|
|
|
|
|
|
|
|
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 |
|
|
|
|
|