Пример. 1. Требуется найти минимум функции , где А=1; В=2; С=1,2; D=2

1. Требуется найти минимум функции , где А=1; В=2; С=1,2; D=2.

2. Интервал поиска: х1нач = -2, х1кон = 2, х2нач = -2, х2кон = 2.

3. Начальная точка: х10 = 1,4962, х20 = -1.

4. Параметры поиска: шаг h=0,5, погрешность = 0,01.

В окрестности заданной начальной точки (1,4962, -1,0000), R=3,9719 сформируем начальный симплекс, его координатами являются:

вершины х1 х2 R
  1,7462 -0,8557 4,1938
  1,2462 -0,8557 1,7519
  1,4962 -1,2887 6,7368

Конечно, можно и по-другому построить начальный симплекс, например, заданная начальная точка может являться одной из вершин симплекса, а не лежать внутри симплекса, как в данном случае. Это не принципиально и на результат существенно не влияет.

Не трудно заметить, что самой плохой вершиной является третья (самое большое значение критерия). Удалим эту вершину и на противоположной грани построим новый симплекс. Практически это означает, что нужно найти одну вершину нового симплекса. Применим базовый алгоритм для каждой переменной х1 и х2 отдельно.

Получим новую вершину:

вершины х1 х2 R
  1,4962 -0,4226 0,8839

Теперь симплекс образуют вершины 1, 2, 4. Самая плохая вершина – 4. По аналогичной формуле перейдём к следующей новой вершине, «отразив» удаляемую относительно центра грани.

Получим вершину:

вершины х1 х2 R
  0,9962 -0,4226 0,4578

Теперь симплекс образуют вершины 2, 4, 5. Самая плохая вершина – 2. Строим следующую вершину:

вершины х1 х2 R
  1,2462 0,0104 1,5919

Оказалось, что только что полученную вершину нужно будет удалять как самую плохую. Следовательно, необходимо сжимать симплекс, воспользовавшись второй формулой алгоритма.

вершины х1 х2 R
  1,2462 -0,6392 0,8750

Далее приведём без комментариев координаты вершин следующих симплексов, предоставив учащимся самостоятельно определять самую плохую вершину в каждом текущем симплексе.

вершины х1 х2 R  
  0,7462 -0,6392 0,6435 отражение
  0,4962 -0,4226 0,3610 отражение
  0,7462 -0,2061 0,3707 отражение
  0,2462 -0,2061 0,1157 отражение
  0,4962 0,0104 0,2526 отражение
  -0,0663 -0,1520 0,0490 сжатие
  -0,3163 0,0645 0,1238 отражение
  0,1056 -0,1385 0,0458 сжатие
  0,0089 0,0704 0,0100 отражение
  -0,1014 0,0814 0,0256 отражение
  0,0546 0,0199 0,0039 сжатие

Симплекс быстро сжимается вокруг искомой оптимальной точки.

При выходе в район оптимума процесс поиска "зацикливает­ся". Это имеет место тогда, когда приходится выбрасывать только что полученную вершину (эта ситуация возникает в том случае, если значение критерия в новой вершине оказывается самое пло­хое). В этом случае можно "сжимать" симплекс, откладывая но­вую вершину от грани на расстоянии вдвое меньше, чем необхо­димо.

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

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

Основным недостатком метода является невозможность уско­рения поиска вдали от оптимума. Этот недостаток устранен в од­ной из модификаций метода, известной как метод Нелдера — Мида. В этом варианте симплексного метода предусмотрено рас­тяжение симплекса в случае, если новая вершина лучше лучшей в старом симплексе, а также сжатия симплекса, если новая вершина оказывается хуже наихудшей в старом симплексе.

Имеется также операция уменьшения размера симплекса, называемая редукцией, которая вводится в случае зацикливания. При этом уменьшаются все грани симплекса одновременно в два раза. Уменьшение размера симплекса в два раза осуществляется делением пополам расстояния от каждой точки симплекса до точки х 1 определяющей наименьшее значение функции, т.е. будет: xi = xic +0,5(xi - x 1).

КОНТРОЛЬНЫЕ ВОПРОСЫ

Метод Гаусса — Зайделя

1. Достаточно ли провести поиск оптимума поочередно по всем переменным последовательно?

2. Область наиболее предпочтительного использования метода Гаусса — Зайделя.

3. Может ли оказывать влияние на результат поиска (значение оптимума) порядок чередования переменных при поиске?

4. Может ли оказывать влияние на эффективность поиска поря­док чередования переменных?

5. Условие окончания поиска min R (x).

6. Основное достоинство метода.

7. Основной недостаток метода.


Понравилась статья? Добавь ее в закладку (CTRL+D) и не забудь поделиться с друзьями:  



double arrow
Сейчас читают про: