
- •Общие сведения
- •Постановка задачи безусловной оптимизации функции многих переменных
- •Глобальный оптимум функции
- •Локальный оптимум функции
- •Выпуклость множеств и функций
- •Методы прямого поиска безусловного оптимума функции многих переменных
- •Эвристические методы
- •Метод Хука-Дживса (метод прямого поиска)
- •Метод покоординатного спуска
- •Метод Розенброка (метод вращающихся координат)
- •Метод Нелдера-Мида
- •Теоретические методы прямого поиска
- •Метод сопряженных направлений (метод Пауэлла)
- •Решение задачи безусловной оптимизации функции многих переменных с использованием необходимых и достаточных условий экстремума
- •Методы оптимизации функции одной переменной
- •Метод квадратичной интерполяции
- •Метод Пауэлла (квадратичной интерполяции)
- •Метод первого приемлемого значения (квадратичная интерполяция с использованием производных).
- •Метод кубической интерполяции (Давидона)
- •Методы оптимизации функции многих переменных
- •Общие критерии останова численных методов оптимизации функции многих переменных
- •Модельная схема применения методов оптимизации функции многих переменных
- •Градиентные методы
- •Аффинная (линейная) модель функции многих переменных
- •Метод Коши
- •Расчет направления поиска методами сопряженных градиентов
- •Метод Ньютона-Рафсона
- •Метод Марквардта
- •Методы переменной метрики (квазиньютоновские методы)
- •Соотношение секущих
- •Метод bfgs
- •Метод dfp
Теоретические методы прямого поиска
Метод сопряженных направлений (метод Пауэлла)
На
k-м этапе метода Пауэлла
используютсяnлинейно
независимых направлений поиска; при
этом поиск начинается в некоторой точкеи проводится следующим образом:
Шаг
1. Начиная из
,
с помощью какого-либо одномерного поиска
определяется
так, чтобы
принимала минимальное значение, и
полагается
.
Начиная из
,
определяется
так, чтобы
обращалась в минимум, и полагается
.
Поиск продолжается последовательно в
каждом направлении, всегда начиная из
самой последней точки последовательности,
пока не будут определены все
,i=1, …,n.
Величина
,
полученная при минимизации
в направлении
,
используется на шаге 4 (см. ниже).
Шаг
2. Пауэлл отметил, что поиск, осуществляемый
в соответствии с шагом 1, может привести
к линейно зависимым направлениям поиска,
как, например, в случае, когда одно из
компонент
становится тождественным нулем, поскольку
в этом направлении не может быть движения.
Следовательно, два направления могут
стать коллинеарными. Таким образом, не
следует заменять старое направление
на новое, если после этого направления
нового набора становятся линейно
зависимыми. Пауэлл показал также (на
примере квадратичной функции), что при
нормировании направлений поиска в
соответствии с формулой
=
1,i= 1, …,n,
определитель матрицы, столбцы которой
представляют собой направления поиска,
принимает максимальное значение тогда
и только тогда, когда
взаимно сопряжены относительно Н. Он
пришел к выводу, что направление полного
перемещения наk-м этапеs(k)должно заменять предыдущее направление
только в том случае, когда заменяющий
вектор увеличивает определитель матрицы
направлений поиска, поскольку только
тогда новый набор направлений будет
более эффективным. Следовательно, после
минимизации
в каждом изnнаправлений,
как на шаге 1, проводится один дополнительный
шаг величиной (
),
соответствующий полному перемещению
наk-м этапе и приводящий
в точку (
).
Затем проводится тест (см. шаг 3), чтобы
убедиться, уменьшается ли определитель
матрицы направлений поиска путем
включения нового направления и
отбрасывания старого.
Шаг
3. Обозначим наибольшее уменьшение
в каком-либо направлении поиска наk-м
этапе через
Направление поиска, соответствующее
этому максимальному изменению
,
обозначим через
.
Чтобы сделать обозначения более
компактными, положим
,
и
,
где
и
.
Тогда если
и (или)
,
то следует использовать на (k+1)-м
этапе, т.е
дляi= 1, …,n,
и начать поиск из точки
[или из
в зависимости от того, в какой точке х
функция
принимает наименьшее значение].
Шаг
4. Если тест на шаге 3 не удовлетворен,
то ищется минимум
в направлении вектора
,
проведенного из
в
,
точка этого минимума берется в качестве
начальной для следующего (k+1)-го
этапа. Система направлений, используемых
на (k+1)-м этапе, та же, что
и наk-м этапе, за исключением
направления
,
которое заменяется на
.
Однако
помещают в последний столбец матрицы
направлений, а не на место
.
Следовательно, на (k+1)-м
этапе будут использоваться следующие
направления
Шаг
5. Критерий удовлетворительной сходимости
для метода Пауэлла, используемый для
определения момента окончания поиска
в конце любого этапа, состоит в том, что
изменение по каждой независимой
переменной должно быть меньше, чем
заданная точность
,i= 1, …,n,
или
.