
plan_exp
.pdf
4 |
+ |
+ |
- |
- |
- |
- |
+ |
5 |
- |
- |
+ |
+ |
+ |
- |
- |
6 |
+ |
- |
- |
- |
- |
+ |
+ |
7 |
- |
+ |
+ |
+ |
- |
- |
- |
8 |
- |
- |
- |
- |
+ |
+ |
+ |
В этом плане предмет заменен фактором. Нетрудно заметить, что план Хаутелинга обладает свойством симметричности.
Такие планы применяются лишь для отсеивающего эксперимента, так как они не дают возможности оценить линейные взаимодействия между факторами.
3. ПЛАНЫ ВТОРОГО ПОРЯДКА
Если математическая модель, полученная с помощью плана первого порядка неадекватна, то необходимо переходить к планам более высокого, например второго порядка. К таким планам относятся ортогональные [27], (в том числе Хартли [37]) и ротатабельные планы [19].
Наиболее простыми являются ортогональные планы они применяются тогда, когда достаточно точно известна область оптимума. Наиболее универсальны ротатабельные планы. Их применение с наибольшей вероятностью обеспечивает получение адекватных моделей, хотя формулы для определения коэффициентов при использовании этих планов более сложны. Планы Хартли, представляющие собой разновидность ортогональных планов, хотя и с меньшей вероятностью обеспечивают получение

адекватных моделей, используются в тех случаях, когда стремятся получить модель процесса при наименьших затратах времени и средств на проведение эксперимента.
Особенность практически всех планов второго порядка состоит в ,томчто в них полный или дробный факторный эксперимент является ядром планирования, к которому добавляют несколько точек, которые могут находиться как в центре, так и на периферии.
3.1. Ортогональные планы
Начнем с двухфакторного эксперимента, ядром которого является план ПФЭ типа 22. Точки такого плана были показаны на рис.2. Дополним этот план точкой в центре эксперимента и четырьмя точками, расположенными на координатных осях на расстоянииa от начала координат (точки 6, 7, 8, 9), как это показано на рис. 4. Координаты точек плана представлены в таблице 10.
Таблица 10.
|
|
|
|
|
|
Ортогональный план второго |
||||||||
|
|
|
|
|
|
|
порядка для двух факторов |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
№ |
|
|
|
Фактор |
|
|
|
|
|
|
|
|||
опыта |
X1 |
X2 |
X1X2 |
2 |
$ |
|
2 |
$ |
2 |
$ |
2 |
$ |
|
|
X1 |
X |
|
X 2 |
X |
(X1 |
- X)(X |
2 |
- X) |
|
|||||
1 |
+1 |
+1 |
+1 |
1 - X |
|
1 - X |
|
(1 - X) |
2 |
|
|
|||
|
|
|
|
|
$ |
|
|
$ |
|
$ |
|
|
|
|
2 |
-1 |
+1 |
-1 |
1 - X |
|
1 - X |
|
(1 - X) |
2 |
|
|
|||
|
|
|
|
|
$ |
|
|
$ |
|
$ |
|
|
|
|
3 |
-1 |
-1 |
+1 |
1 - X |
|
1 - X |
|
(1 - X) |
2 |
|
|
|||
|
|
|
|
|
$ |
|
|
$ |
|
$ |
|
|
|
|
4 |
+1 |
-1 |
-1 |
1 - X |
|
1 - X |
|
(1 - X) |
2 |
|
|
|||
|
|
|
|
|
$ |
|
|
$ |
|
$ |
|
|
|

5 |
0 |
0 |
0 |
|
|
$ |
|
|
|
$ |
|
$ |
|
2 |
|
|
- X |
|
|
- X |
|
X |
|
|
|||||||
6 |
-a |
0 |
0 |
a |
2 |
- |
$ |
|
|
$ |
|
$ |
2 |
$ |
|
|
X |
|
- X |
|
- X(a |
|
|
- X) |
|||||||
7 |
+a |
0 |
0 |
a |
2 |
- |
$ |
|
|
$ |
|
$ |
2 |
$ |
|
|
X |
|
- X |
|
- X(a |
|
|
- X) |
|||||||
8 |
0 |
-a |
0 |
|
|
$ |
|
a |
2 |
- |
$ |
$ |
2 |
$ |
|
|
- X |
|
|
X |
- X(a |
|
|
- X) |
|||||||
9 |
0 |
+a |
0 |
|
|
$ |
|
a |
2 |
- |
$ |
$ |
2 |
$ |
|
|
- X |
|
|
X |
- X(a |
|
|
- X) |
|
X2 |
|
|
+a |
|
|
9 |
|
2 |
+1 |
1 |
a |
0 |
|
+ |
a |
X1 |
- |
+1 |
|
|||
6 -1 |
5 |
7 |
|
|
|
3 |
-1 |
4 |
|
|
|
8
-a
Рис. 4
Левая часть плана, ограниченная жирной линией, представляет собой симметричный план. Если такой план дополнить столбцамиХ12 и Х 22 , то он теряет свойство симметричности, т.к. в нем будут
присутствовать только положительные величины. Поэтому симметричность плана создадим искусст-
|
|
2 |
|
и |
2 |
$ |
2 |
$ |
и |
венно, отнимая от координат Х1 |
|
Х 2 какую-то постоянную величину X . Значения величин Х1 |
- X |
||||||
2 |
$ |
|
|
|
|
|
|
|
|
Х2 |
- X представлены в соответствующих столбцах плана. |
|
|
|
|||||
|
Для симметричного плана |
|
сумма величин |
в каждом столбце равна |
нулю. Если это положение |
||||
|
2 |
$ |
|
|
|
|
|
|
|
применить к столбцу Х1 |
- X , то получим |
|
|
|
|
||||
|
4(1 - X) - 3X + 2(a |
2 |
- X) = 0 |
(35) |
|
|
|
||
|
$ |
$ |
|
|
$ |
|
|
|
|
Если план обладает свойством ортогональности, то сумма произведений величин в двух столбцах тоже, равна нулю. Составим столбец из произведений ( Х12 - X$ )( Х22 - X$ ). Тогда
$ |
2 |
$ |
2 |
$ $ |
2 |
= 0 |
(36) |
4(1 - X) |
|
- 4X(a |
|
- X) + X |
|
Решая совместно (35) и (36), получим значения X$ и a, при которых план является и ортогональным, и симметричным.
Для определения коэффициентов регрессии можно воспользоваться следующими формулами
N
å yn
A ¢ = |
n=1 |
, |
(19a) |
|
|||
0 |
N |
|
|
|
|
|

|
N |
|
|
|
å Yn X1n |
, |
|
A1 = |
n=1 |
||
N |
|||
|
å X12n |
|
|
|
n=1 |
|
|
|
N |
|
|
|
å Yn X2n |
|
, |
A 2 = |
n=1 |
|
|
N |
|
||
|
å x2n2 |
|
|
|
n=1 |
|
|
N0
å Yn X1n X2n
A12 = n=1
N
å X1n2 X22n n=1
(37)
(38)
(39)
Здесь N0 - число опытов в ядре планирования (N0 = 4), N - общее число опытов (N = 9).
По аналогии для коэффициентов при квадратичных членах
|
N |
2 |
- X) |
||
|
|
||||
|
å yn (X1n |
||||
|
|
|
|
$ |
|
A11 = |
n=1 |
|
|
|
(40) |
N |
|
$ |
|
||
2 |
|
2 |
|
||
|
å (X1 |
- X) |
|
|
n=1

|
|
|
|
|
|
N |
|
|
|
2 |
|
$ |
|
|
|
|
|
|
|
|
|
|
|
å yn (X |
|
|
|
|
|
||||||
|
|
|
|
|
A 22 = |
2n |
- X) |
|
|
|
|
||||||
|
|
|
|
|
n=1 |
|
|
|
|
|
|
|
|
|
|
(41) |
|
|
|
|
|
|
N |
|
|
2 |
|
$ |
|
2 |
|
|
|
||
|
|
|
|
|
|
å (X |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
2n |
- X) |
|
|
|
|
|
|||||
|
|
|
|
|
|
n=1 |
|
|
|
|
|
|
|
|
|
|
|
Для планов с числом факторов, равным К в общем случае условие симметричности примет вид |
|||||||||||||||||
|
|
K |
|
|
$ |
|
|
$ |
|
|
2 |
|
$ |
|
, |
(35a) |
|
|
|
2 (1 - X) - (2K - |
1)X + 2(a |
|
- X) = 0 |
||||||||||||
а условие ортогональности |
|
|
|
|
|
|
|
|
|
|
|
||||||
2 |
K |
|
$ |
2 |
$ |
|
2 |
|
$ |
|
+ (2K - |
$ 2 |
= 0 , |
(36a) |
|||
|
(1 - X) |
|
- 4X + 2(a |
|
- X) |
3)X |
Выражения (35a) и (36a) позволяют определить X$ и a для любого числа факторов. Эти значения приведем в таблице 11.
Табл. 11.
Значения параметров ортогональных планов в зависимости от числа факторов.
К |
2 |
3 |
4 |
5 |
6 |
7 |

|
N |
|
|
9 |
15 |
25 |
43 |
77 |
143 |
|
a |
|
|
1 |
1,215 |
1,41 |
1,6 |
1,765 |
1,855 |
|
X |
|
|
0,667 |
0,732 |
0,8 |
0,964 |
0,913 |
0,928 |
|
$ |
|
|
|
|
|
|
|
|
N |
|
$ |
|
|
|
|
|
|
|
å (X |
2 |
2 |
2 |
4,3 |
8 |
13,3 |
19,2 |
24,4 |
|
1n |
- X) |
|
|||||||
n=1 |
|
|
|
|
|
|
|
|
|
Помимо a и X$ , в таблице приведены значения суммы в знаменателе выражений для определения квадратичных коэффициентов регрессии по формуле
N
2 |
$ |
|
|
å yn (Xin - X) |
, |
(42) |
|
A ii = n=N1 |
|
å( Xin2 - X$ ) 2
n=1
атакже общее число опытов, составляющих план для различного числа факторов:
N = 2K + 2K + 1 |
(43) |
Формулы для определения линейных коэффициентов регрессии и коэффициентов при взаимодействиях для К факторов имеют вид:
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
å Yn Xin |
, |
|
|
(37а) |
|
|
|
|
||
|
A i |
= |
n=i |
|
|
|
|
|
|
|
||||
|
N |
|
|
|
|
|
|
|||||||
|
|
|
|
å Xin2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n=i |
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
||
|
|
|
å Yn Xin Yjn |
, |
|
|
|
|
|
|
|
|||
|
A ij = |
n=1 |
|
|
(39а) |
|
|
|
|
|||||
|
N |
|
|
|
|
|
|
|||||||
|
|
|
|
å Xin2 X2jn |
|
|
|
|
|
|
|
|
|
|
|
|
|
n=1 |
|
|
|
|
|
|
|
|
|
||
|
Т.к. квадратичные |
коэффициенты регрессии |
А и А |
рассчитывались для столбцов( Х |
2 |
- Х ) и |
||||||||
|
|
|
|
|
|
|
|
11 |
22 |
|
|
|
1 |
$ |
( Х 2 - Х ), то математическая модель для двух факторов имеет вид |
|
|
|
|
||||||||||
2 |
$ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y = A¢0 + A1X1 + A2 X2 + A12 X1X2 |
|
2 |
- X) + A |
2 |
- X) = |
|
|
||||||
|
+ A11 (X1 |
22 (X2 |
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
$ |
|
$ |
|
|
|
|
$ |
|
|
|
|
|
|
2 |
2 |
|
|
||
|
|
|
|
+ A22 ) + A1X1 + A 2X2 + A12 X1X2 + A11X1 |
+ A22 X2 |
|
|
|||||||
|
|
= A¢0 - X(A11 |
|
|
Отсюда следует, что постоянный коэффициент модели для двух факторов:
A |
0 = A |
' |
$ |
+ A22 ) |
(44) |
0 |
- X(A11 |
В общем случае для любого числа факторов К:

A |
|
= A¢ |
- X |
K |
A |
|
|
|
å |
|
|
||||
|
0 |
0 |
$ |
|
ii |
(44а) |
|
|
|
|
i=1
3.2. Планы Хартли
Планы Хартли являются разновидностью ортогональных планов. Однако, в отличие от планов, рассмотренных в предыдущем параграфе, ядром плана здесь служит не полный, а дробный факторный эксперимент, в котором некоторые линейные эффекты коррелируют с линейными взаимодействиями.
Кроме того, в дополнительных точках опыты ставятся вm повторностях, каждая из которых имеет права опыта. В таблице 12 приведен план для трех факторов, составленный на базе ДФЭ типа23-1, в
котором используется генерирующее соотношение Х=Х Х . Опыты в точках с 5-ой по 11-ую прово-
3 1 2
дятся в m повторностях. Поэтому
N = 2K -1 + m(2K + 1) |
(43а) |
Если m=3, то общее число опытов (строчек плана) N = 23-1 + 3(2× 3 + 1) = 25
Таблица 12
План Хартли для трех факторов
№ |
|
Фактор |
|
опыта |
X1 |
X2 |
X3(X1X2) |

1 |
+ |
+ |
+ |
2 |
- |
+ |
- |
3 |
- |
- |
+ |
4 |
+ |
- |
- |
5 |
0 |
0 |
0 |
6 |
-a |
0 |
0 |
7 |
+a |
0 |
0 |
8 |
0 |
-a |
0 |
9 |
0 |
+a |
0 |
10 |
0 |
0 |
-a |
11 |
0 |
0 |
+a |
12 |
0 |
0 |
0 |
13 |
-a |
0 |
0 |
14 |
+a |
0 |
0 |
15 |
0 |
-a |
0 |
16 |
0 |
+a |
0 |
17 |
0 |
0 |
-a |
18 |
0 |
0 |
+a |
19 |
0 |
0 |
0 |
20 |
-a |
0 |
0 |
21 |
+a |
0 |
0 |