2.6. Метод дробных реплик
С увеличением количества факторов резко возрастает количество опытов полного факторного эксперимента. Это видно из уравнения N=2n. Однако для нахождения коэффициентов регрессии не всегда требуется много опытов. В таких случаях можно уменьшить объем экспериментальных работ, воспользовавшись методом дробных реплик.
Рассматриваемый метод заключается в том, что для нахождения математического описания процесса используется определенная часть полного факторного эксперимента: 1/2, 1/4 и т. д. Эти системы опытов называются дробными репликами
(табл. 2.4).
Таблица 2.4
Полный трехфакторный эксперимент и его дробные реплики
Номер |
|
Факторы |
|
Функция |
Дробные |
|
||||
опыта |
X1 |
|
x2 |
|
x3 |
отклика |
Реплики |
|
|
|
1 |
1 |
|
1 |
|
1 |
y1 |
|
|
|
|
2 |
+1 |
|
1 |
|
1 |
y2 |
¼ 1 |
|
||
3 |
1 |
|
|
|
1 |
y3 |
|
|
|
|
|
+1 |
|
|
2 |
|
|||||
|
|
|
|
|
||||||
4 |
+1 |
|
+1 |
|
1 |
y4 |
¼ |
|
|
|
5 |
1 |
|
1 |
|
+1 |
y5 |
|
|
|
|
6 |
+1 |
|
1 |
|
+1 |
y6 |
¼ |
1 |
|
|
7 |
1 |
|
+1 |
|
y7 |
|
|
|||
|
|
+1 |
|
|
|
|
||||
|
2 |
|
||||||||
8 |
+1 |
|
+1 |
|
+1 |
y8 |
¼ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Расчет коэффициентов регрессии, проверка значимости коэффициентов и адекватности математического описания в данном случае производятся так же, как и при полном факторном эксперименте.
25
Пусть, например, требуется найти коэффициенты уравнения регрессии
y b0 b1X1 b2 X2 b3 X3
Здесь четыре неизвестных коэффициента.
Если для их определения воспользоваться полным трехфакторным экспериментом, то необходимо провести 8 опытов. Однако эту задачу можно решить и с помощью меньшего количества опытов.
Например, возьмем матрицу полного двухфакторного эксперимента (табл. 2.5) и приравняем произведение X1X2 к фактору X3.
Пользуясь формулой (2.14), вычислим:
b0 |
|
|
1 |
|
(y1 y2 y3 y4 ) ; |
||||||||||||
4 |
|
||||||||||||||||
b |
|
|
1 |
|
( y |
|
y |
|
y |
|
y |
|
) ; |
||||
1 |
4 |
1 |
2 |
3 |
4 |
||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||
b2 |
|
1 |
|
( y1 |
y2 |
y3 y4 ) . |
|||||||||||
|
4 |
|
|||||||||||||||
|
|
|
|
|
Таблица 2.5 |
|
|
|
Планирование типа 23 1 |
|
|||
|
|
|
|
|
|
|
Номер |
X1 |
X2 |
X1X2 |
X3 |
Функция |
|
опыта |
отклика |
|||||
|
|
|
|
|||
1 |
1 |
1 |
+1 |
+1 |
y1 |
|
2 |
+1 |
1 |
1 |
1 |
y2 |
|
3 |
1 |
+1 |
1 |
1 |
y3 |
|
4 |
+1 |
+1 |
+1 |
+1 |
y4 |
|
|
|
|
|
|
|
|
|
|
|
26 |
|
|
|
Отметим, что в табл. 2.5 столбцы для произведения X1X2 и фактора X3 полностью совпадают. Поэтому коэффициенты b12 и b3 не могут быть определены раздельно. С помощью формулы (2.14) может быть найдена только их сумма:
1
b12 b3 4 (y1 y2 y3 y4 ).
Этот недостаток рассматриваемого плана является своеобразной «платой» за уменьшение общего количества опытов с восьми до четырех.
Такое планирование эксперимента, когда некоторые из факторов приравнивают к произведениям нескольких факторов, называется планированием со смешиванием. Его обозначают символом 2n p , где n – общее число факторов, а р – число факторов, приравненных к произведениям. С этой точки зрения в табл. 2.5 приведена матрица планирования типа
23 1.
Существует правило, позволяющее определить, какие коэффициенты регрессии определяются совместно при планировании со смешиванием. Рассмотрим это правило на конкретном примере.
Пример. Методом дробных реплик будем искать математическое описание процесса в виде уравнения регрессии:
y b0 b1X1 b2 X2 B3 X3 b4 X4 b5 X5 .
Воспользуемся планированием типа 25 2 и примем:
X4 X1X2;
X5 X1X2 X3.
27
Такие равенства в методе дробных реплик называют
генерирующими соотношениями.
Следует отметить, что выбор генерирующих соотношений в общем случае произволен. Однако он существенно влияет на характер совместных оценок коэффициентов регрессии.
Правило определения совместных оценок коэффициентов состоит в следующем.
1.Примем во внимание, что
Xi2 1;
Xi 1 Xi.
2. Умножив обе части генерирующих соотношений соответственно на X4 и X5, получим:
1 X1X2 X4; 1 X1X2 X3 X5.
Эти равенства называют определяющими контрастами. Перемножив их почленно, получим новые определяющие контрасты. В данном случае это
1 X3 X4 X5.
3. Составим алгебраическую сумму из единицы и правых частей всех полученных определяющих контрастов:
S1 X1X2 X4 X1X2 X3 X5 X3 X4 X5
4.Умножив каждый из факторов на S и заменив факторы соответствующими коэффициентами разложения в ряд Тейлора (2.5), получим
28
b1 1 24 235 1345 ; b2 2 14 135 2345 ; b3 3 1234 125 45 ;
b4 4 12 12345 35; b5 5 1245 123 34.
2.7. Метод ортогонального центрального композиционного планирования
В том случае, когда поверхность отклика не может быть описана многочленом вида (2.9), для адекватного математического описания используется многочлен более высокой степени, например, отрезок ряда Тейлора (2.5), содержащий члены с квадратами переменных. С этой целью используют центральное композиционное планирование эксперимента (ЦКП). Различают два вида ЦКП: ортогональное и ротатабельное.
Количество опытов при ортогональном ЦКП определяется по формуле
N 2n 2n 1, |
(2.22) |
где 2n – количество опытов, образующих полный факторный эксперимент;
2n – число так называемых «звездных» точек в факторном пространстве, имеющих координаты (±α, 0, 0, ..., 0); (0, ±α, 0,
..., 0), ..., (0, 0, ..., ±α). Здесь величина α называется «звездным» плечом;
1 – опыт в центре планирования, т. е. в точке факторного пространства с координатами (0, 0, ..., 0).
Если с помощью полного факторного эксперимента не удается получить адекватного математического описания в форме (2.9), то к нему добавляют опыты в «звездных» точках и
29