Математическое планирование экспериментов. Планирование экстремальных экспериментов. Объект исследования, критерий оптимизации и факторы, страница 12

Поскольку изменению кодированной переменной на одну единицу соответствует изменение натуральной переменной xi на Ii единиц, изменению xi на bi единиц соответствует изменение xi на biIi единиц. Следовательно, для обеспечения шагового движения вдоль линии крутого восхождения, натуральные переменные  xi  должны меняться пропорционально biIi .

Шаговую процедуру крутого восхождения обычно выполняют в несколько циклов. Последовательность проведения каждого цикла следующая :

1.  строят план факторного эксперимента ПФЭ или ДФЭ с центром    в точке x0(x01; x02; …; x0k);

2.  проверяют гипотезу об однородности дисперсий;

3.  проверяют адекватность модели;

4.  вычисляют коэффициенты линейной модели (4.13) и проверяют их значимость.

Проверив полученное уравнение поверхности отклик на адекватность, можно приступить к движению вдоль линии крутого восхождения. Кординаты  j–ой  точки крутого восхождения определяют по формуле

xij  = xi0 +  (biIi/ bmIm)×jr                                  (4.14)

( при i = 1,2, …, m ; j = 1,2, …, k ), где   r - длина шага вдоль i – ой переменной (Табл. 4.10);

xi0 – основной уровень i – ого параметра (назначается);

Ii – интервал варьирования;

j – порядковый номер точки крутого восхождения.

Опыт крутого восхождения                                                           Табл. 4.10

Показатель

X1

X2

Xk

Основной уровень

X10

X20

Xk0

Интервал варьирования

I1

I2

Ik

Изменение откли-

ка на закодиро-

ванную единицу

b1

b2

bk

Изменение ис-

ходной единицы

b1I1

b2 I2

bkIk

Изменение откли-

ка при изменении

на единиц

r

{(b2 I2)/(b1 I1 )}r

{(bk Ik)/(b1I1 )}r

Округление шаг:

1 - ый

X11 = X10+ r

X21 = X20+{(b2 I2)/(b1 I1 )}r

Xk1 = Xk0+{(bkIk)/(b1 I1 )}r

2 - ой

X12 = X10+ 2r

X22 = X20+{(b2 I2)/(b1 I1 )}2r

Xk2 = Xk0+{(bkIk)/(b1 I1 )}2r

...

i - ый

X1i = X10+ ir

X2i = X20+{(b2 I2)/(b1 I1 )}ir

Xki = Xk0+{(bkIk)/(b1 I1 )}ir

Подставляя полученные координаты (4.14) в уравнение (4.13) гиперплоскости, для каждого j – го шага найдем

Y = b0 + b1{(x1 – x10)/I1} + b2{(x2 – x20)/I2} + … + bk{(xk – xr0)/Ik}      (4.15)

Линейное уравнение (4.15) адекватно описывает результаты опытов лишь в экспериментальной области. При значительном удалении от нее адекватность может нарушиться и движение вдоль направления крутого восхождения после последовательного улучшения может привести к уменьшению значения функции отклика.

Для экспериментальной проверки изменения функции отклика в выбранном направлении градиента производится ее реализация на объекте в некоторых точках. Эти точки могут выбираться через 2-3 шага.

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

Координаты точки с наилучшим значением функции отклика принимают за новые нулевые уровни входных переменных. Вновь проводят ПФЭ или ДФЭ, определяют новое линейное приближение в направлении градиента  и опыт крутого восхождения повторяют.

4.1.11.  Симплекс – метод оптимизации

Рассматриваемые методы планирования экстремальных экспериментов ( покоординатного, градиентного, наискорейшего подъема и т.д.)

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

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