Количественные методы в экономических исследованиях. Часть 1. Классические методы оптимизации
.pdf11
2.3. Выбор параметров овражных методов
Начнем с первого овражного метода. Задачей этапа 1 является спуск на дно оврага.
Этот факт можно считать установленным, если на первых шагах этапа имеют место значительные изменения «быстрых» переменных, а на последних шагах их изменения практически отсутствуют. Этап 2 должен обеспечить движение к точке минимума вдоль дна оврага за счет ограниченной расходимости по «быстрым» переменным.
В табл.4 и на рис. 3 приведена траектория первого овражного метода для функции двух переменных, которые четко определены как «быстрая» и «медленная». В табл. 5 и на рис. 4 изображена траектория спуска для функции, обе переменные которой являются
«быстрыми» и обе становятся «медленными» на дне оврага. Хорошо виден спуск на дно оврага на этапе 1 и движение вдоль дна оврага на этапе 2.
Все действия алгоритма ориентированы на наличие оврага. Если овраг отсутствует, то результат этих действий может оказаться отличным от ожидаемого. Это имеет место в непосредственной близости от точки минимума и приводит к тому, что после выполнения очередного этапа 1 фиксируется большее (чем на предыдущем этапе 1) значение функции.
Такое положение означает, что алгоритм уже выполнил свою задачу и получил точку,
близкую к точке минимума. Работу алгоритма можно остановить.
Основная сложность при использовании первого овражного метода заключается в выборе его параметров. Неудачное определение параметров может привести к вырождению его в обычный градиентный или покоординатный спуск. Исходя из изложенного, можно предложить следующую последовательность действий при настройке первого овражного метода.
а) Задать 1 0, 2 максимально возможное значение.
б) Выбором параметров t1 , N1 обеспечить спуск на дно оврага. При этом фиксируется разбиение переменных на «быстрые» и «медленные».
в) Увеличением 1 попытаться остановить движение по «медленным» переменным.
г) Выбором t2 , N 2 обеспечить хорошее изменение «медленных» переменных, не обращая внимание на расходимость по «быстрым» переменным.
д) Уменьшением 2 остановить расходимость по «быстрым» на нужном уровне.
Настройка метода овражного шага заключается в выборе единственного параметра – величины овражного шага t . Из приведенных выше рассуждений следует, что признаком хорошего выбора величины шага может служить значительное уменьшение функции на нескольких первых итерациях с последующим возрастанием функции. В табл. 6 и на рис.5
изображена траектория метода овражного шага. Для наглядности результаты первых трех
12
шагов и последующих шагов изображены на разных рисунках в разном масштабе. Видно,
что увеличение функции (табл. 6) имеет место после того, как алгоритм «побывал» в
непосредственной близости к точке минимума.
13
Приложение 1.
Рис. 1. Градиентный спуск (вариант А): 1а - t=0,001; 1б - t=0,01225.
y |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
5 |
|
2 |
4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
6 |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
10 |
|
|
|
|
|
|
7 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
4,9 |
|
|
|
|
|
|
|
2 |
3 |
4 |
5 |
6 |
7 |
8 |
x |
Рис. 1а
y
5
0
1
2
3
4
4,5 |
5 |
6
7
8
9
10
4
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
x |
Рис. 1б
14
Рис. 2. Наискорейший спуск (вариант В): 2а - итерации 0-2; 2б - итерации 2-10.
y |
|
|
|
|
|
|
|
6,5 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
5,5 |
|
|
|
|
|
|
|
4,5 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
||
3,5 |
|
|
|
|
|
|
|
2,5 |
|
|
|
|
2 |
|
|
1,5 |
|
|
|
|
|
x |
|
1,5 |
2,5 |
3,5 |
4,5 |
5,5 |
6,5 |
||
|
|||||||
|
|
Рис.2а. |
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
5 |
|
|
|
|
2,5 |
|
|
4 |
|
|
|
|
|
|
|
7 |
|
|
||
|
|
|
|
|
|
||
|
|
|
6 |
|
9 |
|
|
|
|
|
|
8 |
|
|
|
|
|
|
|
|
10 |
|
|
1,5 |
|
|
|
|
|
|
|
5,5 |
|
|
6,5 |
|
|
x |
|
|
|
Рис.2б. |
|
|
|
||
15
Рис. 3. Первые два цикла первого овражного метода для варианта А.
y |
|
|
|
|
|
5,6 |
|
|
|
|
|
4,6 |
|
|
|
|
|
|
Первый этап (I) |
|
6 |
|
|
|
|
|
|
|
|
3,6 |
|
|
|
|
(II) |
|
|
|
|
|
|
|
Первый цикл |
|
этап |
||
|
|
|
|
|
|
2,6 |
|
|
|
|
Второй |
|
|
|
|
|
|
1,6 |
|
|
|
|
|
|
Второй цикл |
II |
|
|
|
|
|
|
|
I |
|
0,6 |
|
|
|
|
|
1,6 |
3,6 |
5,6 |
7,6 |
9,6 |
x |
y |
|
I |
|
|
|
|
|
|
|
|
|
1,55 |
|
|
|
|
|
1,35 |
II |
Второй цикл |
|
|
|
|
|
|
|||
1,15
0,95
0,75
7,7 |
8 |
8,3 |
8,6 x |
16
Рис. 4а. Первые два цикла первого овражного метода для варианта B.
y |
|
|
|
|
|
4,9 |
|
|
0 |
|
|
|
|
|
|
|
|
4,4 |
|
|
|
|
|
|
|
|
|
1 цикл |
|
3,9 |
|
|
|
|
|
|
|
|
II |
|
|
3,4 |
|
|
|
|
|
2,9 |
|
|
8 |
|
|
|
|
|
|
|
|
|
1 |
|
|
2 цикл |
|
2,4 |
|
|
|
|
|
1,9 |
|
|
9 |
16 |
|
|
|
|
|
||
1,4 |
|
|
|
|
|
0,9 |
|
|
|
|
|
4,2 |
5,2 |
6,2 |
7,2 |
8,2 |
x |
17
Рис.4б. Второй цикл первого овражного метода для варианта B.
y
2,7
2,5
I
2,3
15
II
2,1
1,9
1,7
1,5
6,1 |
6,6 |
7,1 |
7,6 x |
|
|
|
18 |
|
|
|
|
Рис. 5а. Второй овражный метод (вариант С): итерации 0-3. |
|
|
|||||
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
0,9 |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
0,4 |
|
|
|
|
|
|
|
-0,1 |
|
|
|
|
|
|
|
-0,6 |
|
|
|
|
|
|
|
-1,1 |
|
|
|
|
|
|
|
-1,6 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
-2,1 |
|
|
|
|
|
|
|
2 |
3 |
4 |
5 |
6 |
7 |
8 |
x |
|
|
|
19 |
|
|
|
|
Рис. 5б. Второй овражный метод (вариант С): итерации 3-9. |
|
|
|||||
y |
|
|
|
|
|
|
|
|
|
|
4 |
8 |
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
7 |
|
|
0,99 |
|
|
|
|
9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
|
|
0,98 |
|
|
|
|
|
|
|
0,97 |
|
|
|
|
|
|
|
0,96 |
|
|
|
|
|
|
|
0,95 |
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
0,94 |
|
|
|
|
|
|
|
0,93 |
|
|
|
|
|
|
|
7,5 |
7,6 |
7,7 |
7,8 |
7,9 |
8 |
8,1 |
x |
20
Приложение 2
Табл.1. Траектория спуска метода градиентного спуска для функции
f x, y 70 x - 8 2 y 1 2 1 (вариант А).
Точка минимума: 8,1 .
Уравнение дна оврага:
t=0,01225
t =0,001
x 8 . |
|
|
|
|
|
|
Итераци |
x |
y |
|
f x, y |
|
я |
|
|
|
|
0 |
|
2 |
5 |
2537 |
|
112.29 4.902 1304.51
24.93265 4.8064 674.09
3 |
10.19316 |
4.71314 |
351.48 |
4 |
6.43189 |
4.62217 |
186.25 |
5 |
9.1212 |
4.53343 |
101.48 |
6 |
7.19835 |
4.44686 |
57.87 |
7 |
8.57318 |
4.36241 |
35.30 |
8 |
7.59017 |
4.28003 |
23.52 |
9 |
8.29303 |
4.19967 |
17.25 |
10 |
7.79049 |
4.12128 |
13.81 |
0 |
2 |
5 |
2537 |
1 |
2.84 |
4.992 |
1880.73 |
23.5624 4.984 1395.33
34.1837 4.976 1036.32
44.7179 4.9681 770.77
55.1774 4.9602 574.36
65.5726 4.9522 429.08
75.9124 4.9443 321.61
86.2047 4.9364 242.11
96.456 4.9286 183.30
10 6.6722 4.9207 139.79
