Файл: Мастеров, В. А. Практика статистического планирования эксперимента в технологии биметаллов.pdf

бросим из него худшую точку (пусть, например, это бу­ дет точка 2 ). Построим вместо нее новую точку 5 и но­ вый симплекс 3—5. Удалим из симплекса 3—5 худшую точку ....

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

Экспериментаторы обычно используют правильные симплексы, т. е. фигуры с равными расстояниями вершин от центра симплекса, а начало координат (Х<= 0) сов­ мещают с центром исходного симплекса. Тогда координа­ ты вершин исходного симплекса можно взять из табл. 17. Содержащийся в таблице план обладает свойством орто­ гональности (8), симметричности (12), но

п+1

2 ^ = 0,5,=^=/1+1.

0=1

Поэтому информационная матрица плана имеет вид:

При числе факторов больше двух координаты очеред­ ной точки вычисляют следующим образом:

74

случае возвращаемся к предыдущему симплексу и стро­ им вершину (гг-}-1 )— отражение вершины (и—2).

2. Началось вращение симплексов вокруг одной из вершин (вокруг вершины 11 на рис.21). Это возможно при большой ошибке измерения отклика у или нахожде­ нии симплексов вблизи экстремума у. Ошибки можно устранить дублированием сомнительных опытов. Если подтверждается отсутствие ошибки, в области «зацик­ ливания» симплексов ставится специальный план для изучения области оптимума, например, план второго по­ рядка.

В металлургии метод симплексов нашел широкое при­ менение в изучении свойств сплавов, оптимизации их со­ става [38, 86 и др.].

75

При движении к области оптимума можно использо­ вать и нерегулярные симплексы, а также учитывать зна­ чения отклика во всех точках предыдущего симплекса при выборе новой точки эксперимента [39, 40]. Предпо­ ложим, в вершинах нерегулярного исходного симплекса А, В, С, (рис. 22) значения отклика Уа < У в < У с ■Отбра­ сываем худшую точку А, но при определении координат новой точки Е учтем значения у в точках В и С:

строим такую точку D на противоположной грани, чтобы

BD/DC = ус!ув \

на продолжении отрезка AD строится искомая точка Е так, чтобы A D = D E .

Обобщая рассмотренный метод на случай п-мерного симплекса, получим формулу для расчета координат но­ вой точки:

~( ^ У и + уmin) ^/min •

Полезным приемом ускорения движения к оптимуму является определение по результатам опытов предыду­

76

щего симплекса координат не одной, а нескольких новых точек опытов. Например, в задаче с двумя факторами ис­ ходный симплекс А, В, С дал худшую точку А. Тогда сле­ дующие опыты проводим в точках D и Е, построение ко­ торых ясно из рис. 23: отрезок СЕ равен отрезку ВС, от­ резок CD равен отрезку АС, новые точки лежат на продолжении сторон исходного симплекса. Процедура продолжается, как показано на рисунке.

Таким образом симплекс-планирование: представляет собой хорошо формализованный и про­

стой способ оптимизации; использует при выборе движения к оптимуму самые

свежие наблюдения, что особенно важно при изменении координат оптимума во времени;

позволяет использовать как регулярные (правиль­ ные), так и нерегулярные симплексы; в последнем случае несущественны отклонения координат фактических то­ чек опыта от намеченных.

Однако информация о форме поверхности отклика, получаемая симплекс-методом, ограничена.

ПЛАНИРОВАНИЕ ВТОРОГО ПОРЯДКА

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

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

с увеличенным числом опытов по сравнению с матрица­

77

ми планов первого порядка. Для технолога удобно поль­ зоваться планами второго порядка, получающимися до­ бавлением новых точек к планам первого порядка. Такие планы разработаны и называются композиционными, или последовательными. Если точки плана располагаются

симметрично относительно центра Х = 0 , план называется

центральным.

Центральный композиционный план второго порядка строится следующим образом:

Эти три сферы и образуют план второго порядка. Вы­ бор числа звездных точек и радиуса у оценивается с по­ мощью системы критериев оптимальности планов [41].

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

ОРТОГОНАЛЬНОЕ ЦЕНТРАЛЬНОЕ КОМПОЗИЦИОННОЕ ПЛАНИРОВАНИЕ

На рис. 24 показано расположение точек планов для двух и трех факторов. Величина звездного плеча у равна

1,000 для п = 2, 1,215 для п= 3 и 1,414 для п = 4. Ортого­ нальность плана достигается специальным преобразова­ нием квадратичных переменных и выбором величины

плеча у. В самом деле, если й ^ = ± 1 , то Х1 = + 1 и столб­ цы Х0 и Х\ не ортогональны. Поэтому при расчете коэф­

фициентов регрессии в колонку X£ записывается преоб­ разованная переменная

78

Т а б л и ц а 18

Матрица планирования, результаты опытов и расчетов прочности сварных соединений Х18Н10Т+А99+А М г6

79