Рис. 3.2. Схема шахтной печи: 1-1 -датчик температуры; 1-2 – регистрирующий прибор
Требуют значительно меньшего числа опытов и быстрее приводят к цели те поисковые методы оптимизации, где шаговое варьирование факторами производится целенаправленно по определенному плану. Поисковые методы оптимизации относятся к классу итерационных процедур, при этом весь процесс разбивается на шаги, на каждом шаге делается ряд опытов и определяется, каким образом нужно изменить факторы, влияющие на процесс, чтобы получить улучшение результата. При этом на каждом очередном шаге получаемая информация используется для выбора последующего шага.
Разработано множество методов пошаговой оптимизации, которые подробно рассматриваются в разделе вычислительной математики – “Численные методы оптимизации”. Мы же рассмотрим только некоторые из них, эффективность использования которых в промышленном и
45
лабораторном эксперименте применительно к металлургическим процессам подтверждена практикой.
3.1. Метод покоординатной оптимизации
Процесс поиска оптимума методом покоординатной оптимизации в графическом виде для двумерного случая представлен на рис. 3.3. По этому методу выбирается произвольная точка М и определяются ее координаты. Поиск оптимума осуществляется поочередным варьированием каждого их факторов. При этом сначала изменяют один фактор
(х ) при фиксированных остальных (х |
|
) до тех пор, |
||||
|
|
|
отклика (точка М ). В |
|||
пока не прекращается прирост функции = |
|
(х ) |
|
|||
дальнейшем |
изменяется |
другой |
фактор |
|
при |
|
фиксированных остальных (х |
= |
t), и далее |
процедура |
|||
повторяется. |
|
|
|
|
|
|
Рис. 3.3. К методу покоординатной оптимизации
46
Данный метод весьма прост, однако при большом числе факторов требуется значительное число опытов, чтобы достичь координат оптимума. Более того, при некоторых зависимостях y=f( ,..., ) этот метод может привести к ложному результату. На рис. 3.3 показан один из таких частных случаев, когда поочередное изменение каждого из факторов в любую сторону вдоль координатных осей и вызывает уменьшение y. В результате решения находится ложный экстремум, находящийся в точке А′ с координатами
, ; в то время как действительное значение максимума
находится в точке А с координатами |
и |
. |
методы. |
В дальнейшем рассмотрим более совершенные |
|||
3.2. Метод крутого восхождения
Известно, что кратчайший, наиболее короткий путь — это движение по градиенту, т.е. перпендикулярно линиям равного уровня, на которых функция отклика принимает постоянные значения ( , ,…, )= . В связи с этим при оптимизации процесса рабочее движение целесообразно совершать в направлении наиболее быстрого возрастания функции отклика, т.е. в направлении градиента функции y.
Существуют различные модификации градиентного метода, одним из них является метод крутого восхождения. Сущность этого метода также рассмотрим на примере двухфакторной задачи (рис. 3.4).
Рис. 3.4. Процедура оптимизации методом крутого восхождения
47
В этом случае шаговое движение осуществляется в направлении наискорейшего возрастания функции отклика, т.е. ( , ). Однако направление корректируют не после каждого следующего шага, а при достижении в некоторой точке на данном направлении частного экстремума функции отклика.
Пусть в окрестности точки М как центра плана поставлен ПФЭ 2 . Координаты отдельных опытов соответствуют точкам 1-4. По результатам ПФЭ можно рассчитать коэффициенты линейного уравнения регрессии.
=+ + .
Градиент функции отклика в этой точке определяется как
= |
|
∙ + |
|
∙ , |
(3.1) |
|
|
где , — единичные векторы в направлении координатных осей.
Следовательно, для движения по градиенту необходимо изменять факторы пропорционально их коэффициентам регрессии и в сторону, соответствующую знаку коэффициента. В процессе поиска двигаются в этом направлении до тех пор, пока не будет обнаружен локальный максимум (точка М на рис. 3.4). В точке последнего находят новое направление градиента (направление М ), осуществляя опять же ПФЭ, и далее процедура повторяется. Стрелками на рис. 3.4 показана траектория движения к оптимуму.
Практически алгоритм сводится к следующей последовательности операций.
1. Планирование и постановка ПФЭ (или ДФЭ) в
окрестности точки начального состояния |
. Расчет |
|
коэффициентов |
линейной математической модели с целью |
|
определения направления градиента. |
|
|
|
48 |
|
2. |
Расчет произведений |
, |
где |
∆ |
— интервалы |
варьирования факторов при ПФЭ (∆ДФЭ). |
|
||||
3. |
Выбор базового фактора |
= |
, у которого | ∆ | = |
||
= |
|
|
|
|
|
4. Выбор шага крутого восхождения для базового фактора ha .
Этот выбор производится на основании имеющейся априорной информации или с учетом опыта исследователя, технологических соображений или других критериев. Относительно выбора шага заметим, что слишком малый шаг потребует значительного числа опытов при движении к оптимуму, а большой шаг создает опасность проскакивания области оптимума.
5. Расчет шагов изменения других факторов по формуле
= ( ∆ ) / . |
(3.2) |
Это соотношение между величинами шагов изменения отдельных факторов обеспечивает движение по градиенту в факторном пространстве.
6. Составление плана движения по градиенту. Для этого в соответствии с определенными значениями шагов изменения факторов и их последовательным алгебраическим суммированием с основным уровнем в точке
= + , = 1,2,…
находят координаты опытов 5, 6, 7, 8, 9, 10 (см. рис. 3.4). Часть этих опытов полагают "мысленными".
"Мысленный" опыт заключается в получении предсказанных (расчетных) значений функции отклика по линейному уравнению регрессии, что позволяет сократить объем реальных опытов, т.е. увеличить скорость продвижения к экстремуму. При "мысленном эксперименте" перевод координат в кодированную форму и подстановка их в уравнение модели объекта должна подтвердить
49