Метод градиента


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

Поскольку координатами вектора clip_image006служат, как известно, коэффициенты при линейных членах разложения функции clip_image002[1] в ряд Тейлора по степеням Хi (i=1,…, n), то соответствующие компоненты вектора градиента могут быть получены как коэффициенты b1, b2,…, bn линейной аппроксимации поверхности отклика вблизи исходной точки clip_image008:

clip_image010. (8.12)

Для получения оценок линейных коэффициентов b0, b1,…, bn можно воспользоваться любым известным способом экспериментального получения математической модели объекта, например, реализовать ПФЭ с центром в точке clip_image008[1].

Более простым, хотя и менее точным, является способ определения каждого из коэффициентов bi по результатам двух пробных движений из точки clip_image012 в точки clip_image014:

clip_image016 (8.13)

Тогда соответствующий коэффициент bi найдется по формуле

clip_image018, clip_image020. (8.14)

Таким образом, процедура оптимизации методом градиента может быть выполнена по следующей схеме.

1. Задается шаг варьирования r и параметр рабочего шага R, единые для всех независимых переменных и всех циклов работы. Напомним, что в стандартизированном масштабе (то есть в пространстве типа clip_image022), как правило, clip_image024.

2. Определяется начальная точка clip_image026 в факторном пространстве clip_image028 (как правило, типовой рабочий режим или вообще любая точка, в которой функция отклика удовлетворяет требованиям исследуемого процесса и не нарушается ни одно его ограничения). В этой точке реализуется пробный эксперимент (8.13) для определения направления первого рабочего шага по вектору gradclip_image030. Объём эксперимента равен 2n, то есть производятся 2n измерений отклика в точках clip_image032.

3. По результатам пробного эксперимента с помощью формулы (8.14) вычисляется вектор

clip_image034.

4. Совершается рабочий шаг в направлении вектора clip_image036

clip_image038.

5. В точке clip_image040 процедура поиска градиента и рабочего движения по нему повторяется. Очевидно, что

clip_image042.

6. Поиск прекращается, когда модуль градиента Y становится малой величиной

clip_image044

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

Загрузка...