Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ЛЕКЦИЯ 11

.pdf
Скачиваний:
13
Добавлен:
08.05.2015
Размер:
403.43 Кб
Скачать

~(2)

 

 

1

~(2)

 

1

 

 

1

 

 

1 1.

 

 

 

 

 

 

x1

 

1

 

 

0 1, x2

 

 

 

 

 

 

 

 

 

 

 

2

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как точка

~(2)

попадает в область X, то

x

(2)

~(2)

. Так

 

 

x

 

x

как f x(2) f x(1) ,

 

то проверяем критерий останова:

 

 

x(2) x(1)

 

0.5 0.1.

 

Критерий

останова

не

 

выполняется,

 

 

 

 

 

значит, переходим к следующей итерации.

Итерация 3. В соответствии с выражением

 

(6)

определяем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3)

 

 

~(3)

 

 

 

 

 

 

 

 

 

(2)

 

1

 

(2)

 

 

 

 

 

x

 

X

x

X x

 

 

 

 

 

f

x

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

~(3)

 

1

 

 

 

~(3)

 

 

 

 

 

 

 

1

 

1 1.5.

 

 

 

 

x1

1

 

 

0 1, x2

 

1

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(3)

 

выходит за пределы области X.

 

 

Точка x

 

 

 

 

Так как множество X - шар, то можно указать явный вид

проекции точки на множество.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Можно показать, что если

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X x R n :

 

 

 

x x0

 

 

 

 

r -шар, то a Rn :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

x0

 

 

 

 

 

a x0

 

 

 

r .

 

 

 

 

 

 

 

 

 

X

 

 

 

 

x x0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В нашем случае

 

 

 

 

 

 

 

 

 

 

 

~(3)

0

 

 

 

 

~(3)

 

 

 

 

 

 

(3)

 

 

~(3)

 

 

 

 

x

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

X x

0

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

2 .

 

 

~(3)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

~(3)

2

~(3)

2

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

1.5

 

 

 

 

 

 

 

 

 

 

x(3)

 

 

 

 

 

2 0.7845, x(3)

 

 

 

 

2 1.1767.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

3.25

 

 

 

 

 

 

2

 

 

 

 

 

3.25

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как

f x(3) f x(2) ,

то проверяем критерий останова:

x(3) x(2) 0.2786 0.1.

Критерий останова не выполняется, значит, переходим к следующей итерации.

Итерация 4.

~(4)

~(4)

1.6767.

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

1, x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(4)

(1.0,1.6767) выходит за пределы области X.

Точка x

 

 

 

 

 

 

 

 

 

 

 

~(4)

 

 

 

 

 

 

 

 

(4)

~(4)

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

X x

 

 

 

 

 

 

 

 

 

 

2 .

 

 

 

 

 

 

 

~(4)

2

~(4)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

x2

 

 

 

 

 

x(4)

 

1

 

 

0.7245, x(4)

 

1.6767

 

 

 

 

 

2

 

 

2 1.2148.

 

 

 

 

1

 

1.952

2

1.952

 

 

 

 

 

 

 

 

 

 

Так как f x(4) f x(3) , то проверяем критерий останова: x(4) x(3) 0.071 0.1. Критерий останова выполняется, значит, процесс вычислений закончен.

Пример 3. Решить методом проекции градиента:

f x x2

x2

2x

1

8x

2

min,

x X,

(13)

1

2

 

 

 

 

 

X : x1 2x2 12, x1 x2 8, x 0.

(14)

 

 

Решение. Задачу

(13)-

 

(14)

решим методом

проекции градиента, используя соотношение

(6).

 

 

 

 

Выберем начальное приближение x(0)

0,0 .

 

 

 

 

Пусть шаг k , k 0,1,2,... в выражении

(6) k равен

1

.

2

 

 

 

 

 

 

 

 

 

 

Критерий

останова

выберем

следующий:

 

 

x(k 1) x(k)

 

, k 0,1,2,..., где 0.01.

 

 

 

 

 

 

 

 

 

 

 

Итерация 1. В соответствии с выражением (6) определяем

Точка ~(1) x

 

(1)

 

 

 

~(1)

 

(0)

 

1

 

 

 

 

(0)

 

x

 

X

x

X x

 

 

 

 

 

f x

 

.

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(1)

 

 

1

 

 

~(1)

 

 

 

1

 

8 4.

x1

0

 

 

2 1, x2

0

 

 

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

выходит за пределы области X.

Следовательно, чтобы найти x

(1)

~(1)

, необходимо

 

X x

решить задачу вида (2), то есть:

 

 

 

 

 

 

~(1)

 

 

 

 

2

min, x X,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x x

 

 

 

 

 

 

x1

 

~1

2

x

 

 

 

 

 

~1

 

2

 

x1

 

2 x2

min, x X,

 

x

1

1 2

x

2

4 2

min, x X.

(15)

 

 

 

 

 

 

 

 

 

 

Решением задачи (15) является точка с координатами: x1 0; x2 4.

Следовательно, x(1) 0;4 .

Так как f x(1) f x(0) , то проверяем критерий останова:

x(1) x(0) 4 0.01.

Критерий останова не выполняется, значит, переходим к следующей итерации.

Итерация 2.

 

~(2)

 

 

1

 

 

 

 

~(2)

 

 

 

1

8 4.

 

 

x1

0

 

 

 

 

2 1, x2

 

 

0

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~(2)

выходит за пределы области X.

 

 

 

Точка x

 

 

 

 

Следовательно, чтобы найти x

(2)

 

 

~(2)

,

необходимо

 

X x

 

решить задачу вида (2), то есть:

 

 

 

 

 

 

 

 

 

 

 

x

1 2

x

2

4 2

min, x X.

(16)

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решением задачи

 

 

 

 

 

 

 

(16)

является

точка с

координатами: x

0; x

2

4. Следовательно,

x(2)

0;4 .

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

x(2) x(1) 0 0.01.

Критерий останова выполняется, значит, процесс вычислений закончен.

Задачи для самостоятельного решения.

I. Решить методом проекции градиента следующие задачи:

1. x1x2 12 12 min , 2x1 2x2 x X

X : 0.5 x1 10, 5 x2 10.

x(0) 1,6 .

Ответ: x* 0.5848,5 .

2.

3x

 

2x

 

 

1

x2

x2

x

x

 

min ,

1

2

 

2

 

 

 

 

2

1

2

1

 

x X

 

 

 

 

 

 

 

 

 

 

 

X : 0 x1 3, 0 x2 6.

x(0) 1,2 . Ответ: x* 3,0.5 .

3.

8x

2

4x

2

2x

2

4x

1

x

2

min ,

 

 

 

1

 

2

 

 

x X

 

 

 

 

 

 

 

 

 

 

 

X : 0 x1 4, 0 x2 3.

x(0) 2,1 .

Ответ: x* 1.5,3 .

 

 

4. 4x

1

2x

2

x2

x2

5 min ,

 

 

1

2

x X

 

 

 

 

 

 

X : 1 x1 5, 1 x2 6.

x(0) 4,5 .

Ответ: x* 2,1 .

 

 

 

5.

x2

2x

2

4x

1

2x

2

min ,

 

1

 

2

 

 

x X

 

 

 

 

 

 

 

 

X :1 x1 4,

-2 x2 2.

x(0) 1,0 .

Ответ: x* 2, 0.5 .

6.

x2

2x

2

12x

1

min ,

 

2

 

1

 

x X

 

 

 

 

 

 

X : 0 x1 6,

-1 x2 4.

x(0) 5,3 .

Ответ: x* 3,0 .

7.

4x

1

8x

2

2x

2

2x

2

min ,

 

 

 

 

1

 

2

x X

 

 

 

 

 

 

 

 

 

X : 0 x1 4, 1 x2 3.

x(0) 2,1 . Ответ: x* 1,2 .

8. 6x

1

32x

2

x2

4x

2

min ,

 

 

1

 

2

x X

 

 

 

 

 

 

 

X : 2 x1 7, 3 x2 8.

x(0) 6,6 . Ответ: x* 3,4 .

9. x1x2 50 20 min ,

x1 x2 x X

X : 3 x1 6, 1 x2 4.

x(0) 2,1 . Ответ: x* 5,2 .

10.

1

 

1

 

2x1x2

min ,

 

 

 

 

x1

2x

2

 

x X

X : 1 x1 2,

1 x2 3.

x(0) 0.5,1 .

Ответ: x* 1,0.5 .

II. Решить методом проекции градиента следующие задачи:

1. x2

2x

1

x

2

min ,

 

1

 

 

 

 

 

x X

 

 

 

 

 

 

 

 

 

 

 

 

X : x2

x2

4x 6x

2

7 0

,

1

2

 

 

 

1

 

 

 

x1 0, x2

0.

 

 

 

 

 

 

Ответ: x* 1.178,5.307 .

 

2. 2x

1

3x

2

x2

3x2

min ,

 

 

 

 

 

1

 

2

x X

 

 

 

 

 

 

 

 

 

 

X : x2

x2

2x 2x

2

2 0

,

1

2

 

 

 

1

 

 

 

x1 0, x2

0.

 

 

 

 

 

 

Ответ: x* 1,3 .

 

 

 

 

3. 4x

1

8x

2

x2

2x

2

min ,

 

 

1

 

2

x X

 

 

 

 

 

 

 

X : x

2

x

2

4x 2x

2

3 0,

1

 

 

 

 

2

 

 

1

 

 

 

 

 

x1 0, x2

 

0.

 

 

 

 

 

 

 

 

Ответ: x* 1.4718,2.3119 .

4. x2

x

1

 

2x

2

min ,

1

 

 

 

 

 

 

 

 

 

x X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X : x2

x2

4x 4x

2

4 0,

1

 

 

 

2

 

1

 

 

 

 

x1 0, x2

 

0.

 

 

 

 

 

 

 

 

Ответ: x* 1.044,3.7568 .

5. x2

2x

 

4x

2

min ,

1

 

 

 

 

 

1

 

 

 

 

x X

 

 

 

 

 

 

 

 

 

 

 

 

 

X : x

2

x

2

8x 8x

2

28 0,

1

 

 

 

 

2

 

1

 

 

 

 

x1 0, x2

 

0.

 

 

 

 

 

 

 

 

Ответ: x* 2.703,5.5227 .

6. x2

2x

2

3x2

min ,

2

 

 

 

 

 

 

 

 

1

 

 

x X

 

 

 

 

 

 

 

 

 

 

 

 

 

X : x

2

x

2

2x 4x

2

0,

1

 

 

 

 

2

 

 

1

 

 

 

 

 

x1 0, x2

 

0.

 

 

 

 

 

 

 

 

Ответ: x* 3.1656,1.443 .

7. 6x

1

2x2 x

x

2

 

min ,

 

 

 

 

 

 

 

1

 

 

1

 

 

 

x X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X : x2

x2

4x 4x

2

4 0,

1

 

 

 

2

 

1

 

 

 

 

x1 0, x2

 

0.

 

 

 

 

 

 

 

 

Ответ: x* 2.377,3.964 .

8. x2

x

2

 

x x

2

3x

 

1 min ,

1

 

 

 

 

2

 

 

1

 

 

 

 

1

x X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X : x

2

x

2

2x 8x

 

2

15 0,

1

 

 

 

 

2

 

 

1

 

 

 

 

 

x1 0, x2

 

0.

 

 

 

 

 

 

 

 

Ответ: x* 1.7317,2.7898 .

9. x

1

6x

2

x2

3x

2

3x

1

x

2

min ,

 

 

1

 

2

 

 

x X

 

 

 

 

 

 

 

 

 

 

 

X : x2

x

2

4x 8x

2

18 0,

1

 

2

 

 

1

 

 

x1 0, x2

0.

 

 

 

 

Ответ: x* (2.828,2.8535) .

10. x 2x

2

x

2

min ,

1

 

 

 

2

x X

 

 

 

 

 

 

X : x2

x2

2x 6x

2

9 0,

1

2

 

1

 

 

x1 0, x2

0.

 

 

 

 

Ответ: x* 5.8335,2.0909 .

III. Решить методом проекции градиента следующие

задачи:

 

 

 

 

 

 

 

 

1. x2 x2

x x

2

3x

1

1 min ,

1

2

 

 

1

 

x X

 

 

 

 

 

 

 

 

X : x1 x2

4,

 

 

 

 

4x1 4x2 7,

x1 4x2 2,

x0.

Ответ: x* 1.375,0.375 .

2.

 

1

x2

2x

2

x

x

 

3x

 

2 min ,

 

2

2

2

 

8

1

 

1

 

 

x X

 

 

 

 

 

 

 

 

 

X : x1 x2 3, x1 2x2 4,

x 0.

Ответ: x* 2.571,0.7143 .

3. x12 x22 8x1 10x2 min , x X

X : 3x1 2x2 6 , x 0 .

Ответ: x* 0.3077,2.538 .