Задача оптимизации. Метод крутого восхождения (Бокса — Уилсона).


Если модель сделана в виде полинома второго порядка, то найти оптимальные точки очень просто: т.к. сама модель второго порядка символизирует, что мы находимся в области экстремума. Линейная модель свидетельствует:

1) мы находимся далеко от вершины холма

2) в силу малости известной области требуется выработка некоторых правил, согласно которым будем двигаться к экстремуму.

Существует несколько путей нахождения:

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

2) метод градиента

3) метод случайного поиска

4) метод крутого восхождения (наискорейшего спуска)

5) оптимизация в условиях ограничений

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

Важной особенностью метода является также регулярный статистический анализ результатов экспериментов по мере продвижения к экстремуму.

Пусть у нас k факторов, которые независимы друг от друга.

Алгоритм.

1. В центре clip_image004clip_image006 проводим ПФЭ (ДФЭ) с целью получения градиента. При этом результаты эксперимента подвергаются статистическому анализу, который заключается в следующем:

а) проверка воспроизводимости эксперимента (грубые промахи)

б) проверка значимости коэффициентов линейной модели

в) проверка адекватности линейной модели

2. Вычисляются произведения clip_image008 среди них находится max{clip_image009} =clip_image011 (b — величина относительная, clip_image013 — абсолютная).

3. Для базового фактора выбирается шаг крутого восхождения clip_image015 (только для базового). Дальше будем обозначать его clip_image017.

4. Определяются размеры clip_image004[1]clip_image019

5. Проводятся т.н. «мысленные опыты», которые заключаются в вычислении предсказанных значений целевой функции в определенных точках факторного пространства. clip_image021. Для этого независимые переменные линейной модели изменяются с учетом clip_image023 таким образом, чтобы изображающая точка совершала движение в направлении вектора grad y(clip_image025) полученного в пункте 1 занимая последовательно положения clip_image027.

6. «Мысленные опыты» продолжаются до тех пор, пока выполняется неравенство clip_image029

7. Некоторые из «мысленных опытов» (обычно через каждые 2-3 «мысленных шага») реализуются в виде эксперимента один раз. Для проверки соответствия «мысленных опытов» и реальности.

8. Точка clip_image031, где в реальном опыте получено максимальное значение целевой функции, принимается за новую начальную точку нового цикла крутого восхождения

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

10. Поиск прекращается, когда градиент равен нулю, т.е. все bi оказались незначительными.

Загрузка...