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


При оптимизации по методу Гаусса-Зайделя последовательное приближение к экстремуму осуществляется путем поочередного варьирования каждым параметром до достижения частного экстремума выходной величины. Другими словами, изображающая точка clip_image002 перемещается поочередно вдоль каждой из координатных осей факторного пространства Xi (i=1,…,n), причем переход к новой (i+1) координате совершается по достижению частного экстремума целевой функции clip_image004 на предыдущем направлении, то есть в точке clip_image006, где

clip_image008.

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

Направление движения вдоль (i+1)-й координатной оси выбирается по результатам двух пробных экспериментов (измерений отклика clip_image010 и clip_image012 в окрестностях базовой точки clip_image006[1] (то есть точки частного экстремума по предыдущей i-ой переменной).

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

Рис. 8.2. Последовательность экспериментального движения к экстремуму

методом Гаусса-Зайделя

При практическом использовании метода Гаусса-Зайделя (на примере оптимизации двухфакторного процесса, см. рис.8.2) рекомендуется следующая последовательность операций.

1. Определяется начальная точка clip_image029 движения к оптимуму. В реальных условиях отвечает принятому технологическому режиму.

2. Задается шаг варьирования ?Хi по каждой независимой переменной Xi (i=1,2).

3. Для выяснения направления движения в первом рабочем цикле (вдоль оси Х1) совершается пробное движение с центром в начальной точке. С этой целью из начальной базовой точки clip_image029[1]clip_image032 вариацией параметра X1 по DX1 и -DX1 выполняется два пробных шага

clip_image034;

clip_image036,

где производится однократное измерение отклика

clip_image038, g=1,2.

4. Осуществляется сравнение значений отклика в пробных точках и результаты его выражаются с помощью функции y

clip_image040.

5. Осуществляется первый цикл рабочего движения (с тем же шагом DX1) в направлении возрастания уровня clip_image004[1]. Последовательные точки этого движения, очевидно, будут:

clip_image042

6. После каждого рабочего шага производится измерение значения отклика

clip_image044clip_image046clip_image048…..; clip_image050.

7. Первый цикл рабочего движения прекращается по достижении в некоторой точке clip_image052 частного экстремума целевой функции по соответствующей переменной

clip_image054.

Критерием останова служит выполнения неравенств

clip_image056.

8. Точка clip_image058 является исходной для следующего цикла рабочего движения (вдоль X2) и, одновременно, базовой для новых пробных экспериментов в точках clip_image060

clip_image062

clip_image064.

9. Дальнейшая процедура, связанная с выбором направления и организацией второго цикла рабочего движения к оптимуму, целиком аналогична описанной выше. Изображающая точка clip_image002[1] в факторном пространстве занимает при этом последовательные положения

clip_image067

или

clip_image069

10. По окончании второго цикла рабочего движения переходят к третьему (вновь по оси Х1) и т.д.

Поиск прекращается в некоторой точке clip_image071, любое движение из которой приводит к уменьшению значения выходного параметра. С точностью до максимального шага варьирования (DXi)max это и будет точка экстремума целевой функции clip_image073. Считается целесообразным уменьшать шаг DXi по мере приближения к экстремуму.

Следует заметить, что если в пробном движении по i-й переменной оба шага оказались неудачными clip_image075, то переходят к варьированию следующим (i+1) параметром.

Метод Гаусса-Зайделя достаточно прост для практической реализации, достаточно помехоустойчив. Недостатком его можно счесть не самый короткий путь к точке экстремума. Другим недостатком метода является тенденция к ложной остановке процедуры, если в ходе движения поисковая точка окажется на узком “гребне” (см. рис.8.2). Для полной уверенности в достижении именно точки экстремума (если нет теоретических соображений на этот счёт) рекомендуется дважды и трижды повторить всю процедуру поиска из разных исходных точек и прекратить её при достижении одной и той же области clip_image077.

Загрузка...