Материал: 862

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

 

 

 

 

 

 

 

 

R(xi )

 

1

1

 

 

 

 

 

 

 

 

 

 

 

При этом величина

 

 

 

 

определяется при х1

= 2 и х2

= 1,9.

x2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

X

 

 

Имеем

 

R(xi )

 

= 4

х2 + х1 = 4 1,9 2 = 9,6.

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

Отсюда

х22 = 1,9 – 0,02 9,6 = 1,708.

7.Вычисляем величину целевой функции на 2-м шаге при х12=2

их22=1,708.

R(xi) 22 2 (1,708)2 2 1,708 13,25.

8. Сопоставляя величины целевой функции в точке, откуда шагали (15,02), с точкой, куда шагали (13,25), делаем вывод о том, что шаг оказался удачным.

Теперь по рассмотренному алгоритму продолжаем движение до тех пор, пока функция будет уменьшаться (т.е. изменяться желательным образом), а результаты этих расчетов для краткости изложения материала заносим в табл. 5.

Таким образом, удваивая фактор шага, дошли до 7-го шага, когда целевая функция возросла (3,579). Шаг оказался неудачным. Условимся неудачные шаги обозначать со штрихом (этот шаг будет 7'), чтобы отличить их от шагов удачных. Последним удачным был 6- й шаг (х16 = 2; х26 = –0,627). Уменьшаем фактор шага в два раза (вместо 0,64 берем 0,32) и осуществляем очередной шаг из 6-й точки.

Имеем х17 = 2; х27 = –0,464; R(xi) = 3,502. На 6-м шаге R(x) = 3,532. Это значит, что шаг оказался удачным. Обозначаем его как шаг 7 и продолжаем движение. Уменьшаем фактор шага еще в два раза (0,16) и осуществляем очередную попытку движения из 7-й точки (х18=2; х28=–0,484). При этом целевая функция равна 3,5, а в 7-й точке R(xi)=3,502. Значит, шаг удачный. Уменьшаем еще в два раза фактор шага (0,08) и осуществляем движение из последней (8-й) удачной точки. Девятая точка будет иметь координаты х19 = 2; х29 = –0,489

R(2; –0,489) = 3,5, R(2; –0,489) = 3,5. Целевая функция не ухудшилась. Как видно из результатов расчетов (см. табл. 5), последующее уменьшение фактора шага в два раза (0,04) не дало желаемого результата – функция по сравнению с точкой 9 (3,5) не уменьшилась (3,5). Поэтому определяем частные производные по осям, чтобы выбрать новое направление движения в последней удачной точке – точке 10 (2; –0,49).

40

Отсюда имеем

R(xi )

2 x

x

2

2 2 ( 0,49) 3,51;

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

R(xi )

 

 

 

 

x1

 

 

 

 

 

 

4 x

2

x 4 ( 0,49) 2 0,04.

 

 

 

 

 

 

 

x2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 5

 

 

 

Результаты расчетов(метод релаксаций)

 

 

 

 

Номер шага

 

 

 

х1

 

 

х2

R(х1; х2)

Фактор шага

0

 

 

 

 

2

 

 

 

 

2

16

0,01

1

 

 

 

 

2

 

 

 

 

1,9

15,02

0,02

2

 

 

 

 

2

 

 

 

1,708

13,25

0,04

3

 

 

 

 

2

 

 

 

1,354

10,375

0,08

4

 

 

 

 

2

 

 

 

0,76

6,675

0,16

5

 

 

 

 

2

 

 

-0,046

3,912

0,32

6

 

 

 

 

2

 

 

-0,627

3,532

0,64

7'

 

 

 

 

2

 

 

-0,301

3,579

0,32

7

 

 

 

 

2

 

 

-0,464

3,502

0,16

8

 

 

 

 

2

 

 

-0,484

3,5

0,08

9

 

 

 

 

2

 

 

-0,489

3,5

0,04

10

 

 

 

 

2

 

 

 

-0,49

3,5

0,01

11

 

 

 

 

1,965

 

 

 

-0,49

3,379

0,02

12

 

 

 

 

1,896

 

 

 

-0,49

3,145

0,04

13

 

 

 

 

1,763

 

 

 

-0,49

2,725

0,08

14

 

 

 

 

1,52

 

 

 

-0,49

2,04

0,16

15

 

 

 

 

1,112

 

 

 

-0,49

1,171

0,32

16

 

 

 

 

0,557

 

 

 

-0,49

0,517

0,64

17

 

 

 

 

0,157

 

 

 

-0,49

0,427

1,28

18'

 

 

 

 

0,382

 

 

 

-0,49

0,439

0,64

18

 

 

 

 

0,205

 

 

 

-0,49

0,422

0,01

19

 

 

 

 

0,206

 

 

 

-0,49

0,421

0,02

20

 

 

 

 

0,206

 

 

 

-0,49

0,421

0,01

21

 

 

 

 

0,206

 

 

-0,472

0,39

0,02

22

 

 

 

 

0,206

 

 

-0,438

0,385

0,04

23

 

 

 

 

0,206

 

 

-0,376

0,246

0,08

24

 

 

 

 

0,206

 

 

-0,272

0,114

0,16

25

 

 

 

 

0,206

 

 

 

-0,13

0,049

0,32

26

 

 

 

 

0,206

 

 

 

0,029

0,048

0,32

27

 

 

 

 

0,206

 

 

-0,074

0,037

0,16

28

 

 

 

 

0,206

 

 

-0,045

0,037

0,08

29

 

 

 

 

0,206

 

 

 

0,043

0,037

0,01

30'

 

 

 

 

0,206

 

 

-0,043

0,037

0,01

30

 

 

 

 

0,202

 

 

-0,043

0,036

0,02

31

 

 

 

 

0,195

 

 

-0,043

0,033

004

32

 

 

 

 

0,181

 

 

-0,043

0,028

0,32

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

41

 

 

 

Окончание таблицы 5

Номер шага

х1

х2

R(х1; х2)

Фактор шага

33

0,079

-0,043

0,006

0,32

34

0,042

-0,043

0,003

0,32

35

0,028

-0,043

0,003

0,16

36'

0,028

-0,066

0,007

0,08

36

0,028

-0,024

0,0012

0,04

37'

0,028

-0,038

0,0026

0,01

37''

0,028

-0,041

0,0026

0,01

38'

0,0295

-0,041

0,003

 

Наибольшая производная — по оси х1 (3,51), поэтому выбираем новое направление движения по координате х1, а х2 фиксируем на уровне –0,49 и продолжаем вычисления.

11

10

 

R(x

i

)

 

х1

= х1

- h

 

 

 

2 0,01 3,51 1,965,

x

 

 

 

 

 

 

 

 

10

 

 

 

1

 

 

X

1

 

 

 

 

 

 

 

R(1,965; -0,49) = 3,379.

Поскольку на 10-м шаге R(x1)=3,5, то шаг оказался удачным. Удваивая фактор шага, мы доходим до шага 18, при котором исследуемая функция увеличилась по сравнению с предыдущим (17- м) шагом (0,439>0,427). Обозначаем неудачный шаг номером 18', уменьшаем фактор шага вдвое (вместо 1,28 берем 0,64) и продолжаем движение из последней удачной точки (т. 17). Шаг оказался удачным (0,427>0,422). Из точки 18 осуществляем 19 и 20-й шаги с фактором шага 0,02 и 0,01 соответственно. Поскольку шаги 19 и 20 дали одинаковые отклики (0,421), а фактор шага 20 был равен заданной точности поиска (0,01), то в точке 20 вновь определяем частные производные по осям координат:

R(xi ) 2 x1 x2 2 0,206 ( 0,49) 0,078;

x1

R(xi) 4 x2 x1 4 ( 0,49) 0,206 1,754.

x2

На этом основании фиксируем координату х1=0,206 и продолжаем движение из точки 20 с фактором шага h= =0,01, который удваиваем на каждом шаге. Доходим до 26-го шага, затем уменьшаем фактор шага, начиная с 27-го шага до 30-го. Точки 27, 28, 29 и 30 имеют одинаковые отклики (0,037) при факторе шага в точке

30h= =0,01.

Вточке 30 вновь определяем частные производные по осям:

42

R(xi ) 2 x1 x2 2 0,206 ( 0,043) 0,369;

x1

R(xi ) 4 x2 x1 4 ( 0,043) 0,206 0,034.

x2

Считаем точку реверса (т. 30) как неудачную, присваиваем ей номер 30' и движемся далее, вначале увеличивая, а с 35-го шага уменьшая фактор шага. Последний удачный шаг 36, R(xi)=0,0012. Изменение фактора шага и направления движения (шаги 37', 37'', 38') не дает желаемого результата, поэтому в качестве решения принимается последняя удачная точка (точка 36).

Ответ: R(0,028 0,01; -0,024 0,01)=0,0012.

Метод крутого восхождения.

Сущность этого метода сводится к следующему.

1. Вычисляем величины первой производной по осям в исходной точке, т.е. при х10 = х20 = 2 (см. табл. 6).

 

 

R(x)

 

 

 

 

 

 

 

 

2x1 x2

2 2 2 6;

 

 

 

 

 

x

 

 

0

 

 

1

 

X

 

 

 

 

R(x)

 

 

 

 

 

 

 

4 2 2 10.

 

 

 

 

x2

 

 

 

 

 

 

0

 

 

X

 

2. Вычисляем соотношение первых производных по осям

 

 

R(x)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

0

 

 

 

 

 

 

 

X

 

 

6

 

 

 

1

1

 

0,6.

 

 

 

R(x)

 

 

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

0

 

 

 

 

 

 

 

1

X

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Вычисляем значение целевой функции в исходной точке

R(xi) 22 2 22 16.

 

Принимаем шаг по оси х2

 

равным

х2 = 0,5 (по мере

приближения к оптимуму величины производных будут уменьшаться, а поэтому будем уменьшать шаги поиска), тогда х1 = 0,5 0,6 0,3.

4. Осуществляем первый шаг, т.е. вычисляем координаты

x11 и x12 ;

43

x11 x10 Δx1 2 0,3 1,7;

х12 x20 Δx2 2 0,5 1,5.

5. Вычисляем величину целевой функции в новой точке

R(x11, x12 ) (1,7)2 2 (1,5)2 1,5 1,7 9,94.

Дальнейшие вычисления заносим в табл. 6. Движение осуществляем до локального минимума, т.е. до минимальной величины функции на данном направлении. Исходя из этого, на 7-м шаге наблюдается увеличение целевой функции (см. табл. 6), поэтому возвращаемся к данным шага 6 и вычисляем вновь первые производные для х15 = 0,5 и х25 = –0,5. Имеем

 

R(x)

 

 

 

 

0,5;

 

 

x

 

 

 

5

 

 

 

1

X

1

 

 

 

 

 

 

 

R(x)

 

 

 

 

1,5;

 

 

x2

 

 

 

 

 

5

X

 

 

 

 

 

 

2

 

R(x5)

1

 

 

отношение

1

 

 

 

, т.е. изменение величин по осям как 1: 3.

R(x25)

 

 

3

 

R(x)

Здесь следует учесть, что

x2

1,5, т.е. имеем отрицательный

X 52

знак, значит, х2 будем суммировать с х2 в отличие от предыдущего направления, где эта величина вычиталась.

Принимаем х2 = –0,3, тогда х1 = 0,1. Присваиваем неудачному шагу 7 штриховой индекс (так мы условились обозначать неудачные шаги) и движемся дальше по вышеизложенной методике, занося результаты расчетов в табл. 6. Как видно из табл. 6, 9-й шаг оказался неудачным, поскольку 0,44> 0,4.

Рассчитываем первые производные в последней удачной точке, т.е. в точке 8:

 

R(x)

 

 

0,7;

 

R(x)

 

0,7.

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

x2

 

 

 

 

 

7

 

 

 

 

7

X

 

X

 

 

 

 

 

44