в центре плана, а полученную при этом композицию используют для получения математического описания процесса в виде многочлена второй степени. Отсюда и произошло название метода: центральное композиционное планирование.
Значения «звездного» плеча α для ЦКП с различным числом факторов n следующие:
n |
|
2 |
|
3 |
4 |
5 |
|
|
|||||
|
|
|
|
|
|
|
|
|
1,000 |
|
1,215 |
1,414 |
1,547 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Эти значения α выбраны из условия ортогональности матрицы планирования. Уравнение регрессии при ортогональном ЦКП ищут в следующем виде:
y b* b X |
1 |
b |
X |
2 |
... b |
n |
X |
n |
|
|||
0 1 |
|
2 |
|
|
|
|
|
|
||||
b12 X1X2 |
|
... b(n 1)n Xn 1 Xn |
|
(2.23) |
||||||||
b |
X* ... b |
nn |
X *. |
|
|
|
|
|||||
11 |
|
|
1 |
|
|
|
n |
|
|
|
|
|
Переменные величины
|
|
|
1 |
N |
|
|
X*ji |
X 2ji |
|
X 2ji |
|
||
N |
(2.24) |
|||||
|
|
|
j 1 |
(j -номеропыта;i-номерфактора)
введены для того, чтобы матрица планирования была ортогональна и коэффициенты регрессии определялись независимо друг от друга по результатам опытов.
Для того, чтобы получить уравнение регрессии в обычной форме
30
y b0 b1X1 b2X2 ... bnXn
b12X1X2 ... b(n 1)nXn 1Xn |
|
|
(2.25) |
b11X12 b22X22 ... bnn X2n ,
находят величину
|
b |
N |
|
b |
N |
|
|
b0 b0* |
11 |
j 1 |
X2ji ... |
nn |
j 1 |
X 2jn. |
(2.26) |
N |
N |
В табл. 2.6. приведена в качестве примера матрица ортогонального ЦКП для двух факторов, а на рис. 2.3 изображена схема этих опытов.
Таблица 2.6
Ортогональное ЦКП для двух факторов
Система |
Номер |
X1 |
X2 |
X1X2 |
X1* |
X2* |
|
опытов |
опыта |
|
|
|
|
|
|
Полный |
1 |
1 |
1 |
+1 |
+0,33 |
+0,33 |
|
факторный |
2 |
1 |
1 |
1 |
+0,33 |
+0,33 |
|
эксперимент |
3 |
1 |
+1 |
1 |
+0,33 |
+0,33 |
|
Опыты |
4 |
1 |
+1 |
+1 |
+0,33 |
+0,33 |
|
5 |
+1 |
0 |
0 |
+0,33 |
-0,67 |
||
в звездных |
6 |
1 |
0 |
0 |
+0,33 |
-0,67 |
|
точках |
|||||||
7 |
0 |
+1 |
0 |
0,67 |
+0,33 |
||
|
|||||||
Опыт в |
8 |
0 |
1 |
0 |
0,67 |
+0,33 |
|
центре плана |
9 |
0 |
0 |
0 |
0,67 |
0,67 |
|
|
|
|
|
|
|
|
31
|
X2 |
||
|
*+1 |
||
–1* |
|
|
*+1 X 1 |
|
|
||
|
|
||
|
–1* |
||
Рис. 2.3. Схема опытов ортогонального ЦКП для двух факторов:
— опыты полного факторного эксперимента; * – опыты в звездныхточках; 
— опыт в центре плана
Коэффициенты регрессии при ортогональном ЦКП рассчитываются по следующим формулам:
|
|
* |
|
1 |
N |
|
|
|
|
b |
|
|
|
j 1 |
y |
, |
(2.27) |
|
|
N |
||||||
|
0 |
|
j |
|
|
|||
|
N |
|
|
|
|
|
|
|
|
X ji yj |
|
|
|
|
|||
bi |
j 1 |
|
|
|
(гдеi 0), |
(2.28) |
||
N |
|
|
|
|||||
|
(X ji )2 |
|
|
|
|
|||
j 1
32
|
N |
|
|
|
X ji X jk yj |
|
|
bik |
j 1 |
(гдеi k), |
(2.29) |
N |
(X ji X jk )2 j 1
N
X*ji yj
bii |
j 1 |
. |
(2.30) |
N |
(X*ji )2
j 1
Для расчета оценок дисперсий в определении коэффициентов регрессии используют следующие выражения:
|
|
|
|
|
|
|
|
|
|
|
S2 |
|
|
|
|
|
|
|
|
|
|
|
S2* |
|
|
y |
|
, |
|
||||
|
|
|
|
|
|
|
N |
|
||||||||
|
|
|
|
|
|
b0 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
N |
||
|
|
Sb2 Sb2* |
nSb0 |
|
|
X 2ji , |
||||||||||
|
|
|
|
|
|
|
||||||||||
|
|
|
0 |
|
0 |
|
|
|
N |
|
|
j 1 |
||||
|
|
|
|
|
|
S2 |
|
|
|
|
|
|
|
|
|
|
Sb2 |
|
|
|
y |
|
|
|
|
|
(гдеi 0), |
||||||
|
N |
|
|
|
|
|
|
|
||||||||
|
i |
|
|
|
(X ji )2 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
j 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S2 |
|
|
|
|
|
|
|
|
|
|
S2 |
|
|
|
|
y |
|
|
|
|
|
|
|
(гдеi k), |
|||
|
N |
|
|
|
|
|
|
|
|
|
||||||
bjk |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
(X ji X jk )2 |
|
|
|
|
||||||||
|
|
|
|
j 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
S2 |
|||||
|
|
|
|
|
Sb2 |
|
|
|
|
y |
|
. |
||||
|
|
|
|
|
|
N |
|
|
|
|
||||||
|
|
|
|
|
ii |
|
(X |
|
*ji )2 |
|||||||
|
|
|
|
|
|
|
|
|||||||||
j 1
(2.31)
(2.32)
(2.33)
(2.34)
(2.35)
33
Коэффициент bi, считается значимым, если |
bi |
Sb t ,n . |
|
|
i |
Аналогично проверяется значимость остальных коэффициентов регрессии. Проверка адекватности уравнения регрессии осуществляется с помощью критерия Фишера (2.18).
2.8. Метод ротатабельного планирования
Метод ротатабельного планирования эксперимента позволяет получать более точное математическое описание поверхности отклика по сравнению с ортогональным ЦКП, что достигается благодаря увеличению числа опытов в центре плана и специальному выбору величины «звездного» плеча α. В табл. 2.7 приведены основные характеристики матриц ротатабельного планирования.
При ротатабельном ЦКП для вычисления коэффициентов регрессии и соответствующих оценок дисперсий находят следующие константы:
A |
1 |
|
, |
2B (n 2)B n |
|||
B |
nN |
, |
|
|
|
||
(n 2)(N N0 )
C N ,
N N0
где n – число факторов;
N – общее число опытов ротатабельного ЦКП; N0 – число опытов в центре плана.
(2.36)
(2.37)
(2.38)
34