plan_exp
.pdf2 |
|
S 2 |
|
S 2 |
2 |
æ |
1 |
|
1 |
ö |
(58) |
SA |
= |
|
+ |
|
= S |
ç |
|
+ |
|
÷ |
|
|
|
|
|
||||||||
|
ij |
N0m¢ 2a 2m |
|
è N0 m¢ |
|
2a 2mø |
|
||||
Для остальных коэффициентов можно пользоваться формулами (53), (54), (55), (56) и (28).
3.4.3. Ротатабельные планы.
Дисперсия в определении коэффициентов регрессии для ротатабельных планов определяется по общей формуле
SA |
= |
bS2 |
(59) |
|
m |
||||
|
|
|
Коэффициенты b приведены в таблице 15. Индекс при b определяется для различных коэффициентов регрессии в соответствии с таблицей 16.
Таблица 16.
Выбор коэффициентов bк формуле (59):
S2 |
S2 |
S2 |
S2 |
S2 |
A |
A0 |
A i |
A ij |
A ii |
b |
b0 |
bi |
bij |
bii |
4.ДРУГИЕ ПЛАНЫ
4.1.Планы типа 3К
Выше отмечалось, что для получения математической модели необходим такой план эксперимента, у которого число опытов N было больше количества коэффициентов математической модели С. Для двух факторов математическая модель третьего порядка имеет вид:
Y = A0 + A1 X1 + A 2 X2 + A12 X1 X2 + A11 X12 + A 22 X22 + A111X13 + A 222 X32 (60)
При этом взаимодействиями третьего порядка пренебрегаем. Условию N>С удовлетворяет план типа 32, т.к. С=8.
В общем для К факторов может быть использован план типа 3к. Соотношения между N и С для различного числа факторов представлены в табл. 17.
Таблица 17.
Характеристики плана 3к.
№ п/п |
К |
N |
С0 |
Сi |
Сii |
Сij |
Сiii |
С |
1 |
1 |
3 |
1 |
1 |
1 |
0 |
1 |
4 |
2 |
2 |
9 |
1 |
2 |
2 |
1 |
2 |
8 |
3 |
3 |
27 |
1 |
3 |
3 |
3 |
3 |
13 |
4 |
4 |
81 |
1 |
4 |
4 |
6 |
4 |
19 |
5 |
5 |
243 |
1 |
5 |
5 |
10 |
5 |
26 |
На рис.5 показаны опытные точки плана 32, а в табл.18 - план, составленный по этим точкам.
|
Х2 |
|
3 |
6 |
9 |
+1 |
|
|
|
|
|
|
|
|
|
|
|
8 |
|
|
|
|
|
2 |
|
|
0 5 |
|
|
Х1 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-1 |
|
|
|
|
|
|
|
|
7 |
+1 |
|
||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
4 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-1 |
|
|
|
|
|||
Рис. 5.
Таблица 18.
План эксперимента типа 32
№ |
|
|
Фактор |
|
|
|
|
опыта |
Х1 |
Х2 |
Х1Х2 |
2 |
€ |
2 |
€ |
|
Х1 - |
Х |
Х2 - |
Х |
|||
1 |
- |
- |
+ |
1- Х |
1- Х |
||
2 |
- |
0 |
0 |
$ |
|
$ |
|
1- Х |
- Х |
|
|||||
3 |
- |
+ |
- |
$ |
|
$ |
|
1- Х |
1- Х |
||||||
|
|
|
|
$ |
|
$ |
|
4 |
0 |
- |
0 |
- Х |
|
1- Х |
|
5 |
0 |
0 |
0 |
$ |
|
$ |
|
- Х |
|
- Х |
|
||||
6 |
0 |
+ |
0 |
$ |
|
$ |
|
- Х |
|
1- Х |
|||||
|
|
|
|
$ |
|
$ |
|
7 |
+ |
- |
- |
1- Х |
1- Х |
||
8 |
+ |
0 |
0 |
$ |
|
$ |
|
1- Х |
- Х |
|
|||||
9 |
+ |
+ |
+ |
$ |
|
$ |
|
1- Х |
1- Х |
||||||
|
|
|
|
$ |
|
$ |
|
å |
0 |
0 |
0 |
0 |
|
0 |
|
Используя прием, изложенный в п. 3.1 дополним план столбцами ( Х12 - Х$ ) и( Х22 - Х$ ), что позволит записать условие симметричности:
6(1 - X) - 3X |
(61) |
|
$ |
$ |
|
Откуда Х$ =2/3. Интересно отметить, что значение Х$ будет одинаково для всех планов типа 3К. Нетрудно доказать, суммируя значения, представленные в последнем столбце, что план 32 (как и все планы 3К), является ортогональным. Поэтому все формулы для определения коэффициентов регрессии, представленные в п.3.1, пригодны и для планов 3К. Необходимо лишь добавить выражение для определения коэффициента при Х 3i .
Такой коэффициент можно определить по формуле:
N
|
|
2 |
$ |
|
|
|
|
å yn ( Xin - X)Xin |
|||||
A iii = |
n=1 |
|
|
|
|
(62) |
N |
$ |
|
2 |
2 |
||
2 |
|
|
||||
|
å( Xin |
- X) |
|
Xin |
||
n=1
4.2. Симплекс-планирование
К-мерным симплексом называют фигуру, которая в К-мерном пространстве имеет К+1 вершину. В частности, симплекс в двухмерном пространстве(т.е. на плоскости) это фигура, имеющая К+1=2+1=3 вершины. Такой фигурой является треугольник. В трехмерном пространстве симплекс имеет четыре вершины, и, следовательно, это пирамида.
Симплекс, все стороны которого равны, называется регулярным. Для плоскости это равносторонний треугольник, для трехмерного пространства - тетраэдр. Если поставить опыты в вершинах симплекса, то число вершин будет равно числу коэффициентов модели первого порядка. Такой план называют насыщенным. Насыщенные планы нашли применение при планировании отсеивающего эксперимента. На рис.6 приведен симплекс-план для двух факторов. Здесь центр эксперимента, расположен-
ный в точке 0, совпадает с центром тяжести треугольника. При этом одна из высот треугольника совпадает с осью Х2. Если принять сторону треугольника а=1, то координаты точек 1 и 2 по оси Х1 будут
|
|
|
|
|
|
|
|
|
|
|
|
|
соответственно -0,5 и +0,5. При этом высота треугольника h= 12 - 0,52 |
= |
|
0,75 = 0,867 . Т.к. центр тяже- |
|||||||||
сти отсекает от высоты третью часть, |
то координата |
точек 1 |
и |
2 по оси Х2 будет равна |
||||||||
h / 3 = - |
0,867 |
= -0,289 , а координата точки 3 по той же оси |
2h |
= |
2 |
×0,867 = 0,578 . |
||||||
|
|
|
||||||||||
3 |
|
3 |
|
3 |
|
|
|
|
||||
I
Х1
Рис.6
Такой план позволяет проводить опыты в условиях ограничений, когда, например, невозможно поставить опыт в области, расположенной выше линииI-I. Ориентируя симплекс можно добиться удобного его расположения в исследуемой области факторного эксперимента.
В табл. 19. представлен универсальный симплекс-план. Таблица построена таким образом, что дает возможность получить план для любого числа факторов для одного(2 опыта), двух (3 опыта) и т.д. В общем случае для К-факторов (К+1 опыт). Каждый такой план выделен двойной линией.
Таблица 19.
Универсальная таблица симплекс-плана.
№ |
|
|
|
Фактор |
|
|
|
|
опыта |
Х1 |
Х2 |
Х3 |
Х4 |
Х5 |
... |
|
ХК |
1 |
-0,5 |
-0,289 |
-0,204 |
-0,158 |
-0,129 |
|
|
|
2 |
+0,5 |
-0,289 |
-0,204 |
-0,158 |
-0,129 |
|
|
|
3 |
0 |
0,578 |
-0,204 |
-0,158 |
-0,129 |
|
|
|
4 |
0 |
0 |
+0,612 |
-0,158 |
-0,129 |
|
- |
1 |
|
|
|
|
|
|
|
2K(K + |
|
|
|
|
|
|
|
|
|
|
5 |
0 |
0 |
0 |
0,632 |
-0,129 |
6 |
0 |
0 |
0 |
0 |
+0,64 |
|
|
|
|
|
5 |
К |
0 |
0 |
0 |
0 |
0 |
Коэффициенты регрессии:
|
|
|
N |
|
N |
|||
|
|
|
å Yn |
|
å Yn |
|||
A 0 |
= |
n=1 |
= |
n=1 |
|
|||
N |
K + 1 |
|||||||
|
|
|
|
|||||
|
N |
|
|
|
||||
|
å Yn Xin |
|
N |
|||||
A i = |
n=1 |
= 2å Yn Xin |
||||||
|
N |
|
||||||
|
|
å Xin2 |
|
n=1 |
||||
n=1
(63)
(64)
4.3. Отсеивающий эесперимент. Метод случайного баланса
По соотношению между числом опытов и числом факторов планы делятся на насыщенные, ненасыщенные и сверхнасыщенные.
Насыщенными называются планы, у которых N=К+1, т.е. планы, которые позволяют определить коэффициенты модели первого порядка.
Планы, у которых N>К+1 называются ненасыщенными. Такие планы, например, рассмотренные выше планы взвешивания, симплекс-планы, планы ДФЭ позволяют получить более сложные модели или более того оценить адекватность моделей первого порядка, могут быть использованы для отсеивающего эксперимента.
Планы, у которых N<К+1 называют сверхнасыщенными. Они также нашли применение при отсеивающем эксперименте. Для этой цели используется метод случайного баланса, полученный Саттеравайтом [38] и обоснованный Мавлотовым [39]. Рекомендуется в тех случаях, когда заведомо известно, что число значимых факторов гораздо меньше, чем число факторов взятых под подозрение.
Этот метод имеет малую чувствительность, т.е. способность отличать малую величину от нуля. Вместе с тем он позволяет при малых затратах труда отделить значимые коэффициенты от незначимых. Метод применяется в следующей последовательности.
1. Все факторы разбиваются на группы, обычно близкие по физической природе. Как правило в каждую группу входит от 2-х до 6-ти факторов. За основу плана берутся планы ПФЭ и ДФЭ, предпочтительнее первые. Рассмотрим применение метода на конкретном примере.
Пусть общее число факторов К=6. Разобьем их на 2 группы, по три фактора в каждой. За основу примем план типа 23, представленный в табл.3. Пусть таблица случайных чисел выдала нам следующую последовательность: 7, 3, 5, 6, 4, 8, 2, 1. Тогда вместо табл.3 мы получим табл.20, которая отличается от исходной только порядком следования опытов.
Воспользуемся таблицей случайных чисел еще раз. Пусть теперь эта последовательность имеет вид 7, 6, 5, 3, 4, 2, 8, 1. Тогда план будет соответствовать табл.21 для факторов Х4, Х5, Х6. Состыкуем эти два плана в один, представленный в табл.22.
Таблица 20.
Вариант плана 23 (последовательность №2)
№ |
|
|
Фактор |
|
|
п/п |
X1 |
|
X2 |
X3 |
|
1 |
- |
|
- |
- |
|
2 |
- |
|
- |
+ |
|
3 |
+ |
|
+ |
- |
|
4 |
- |
|
+ |
- |
|
5 |
+ |
|
- |
+ |
|
6 |
+ |
|
- |
- |
|
7 |
- |
|
+ |
+ |
|
8 |
+ |
|
+ |
+ |
|
|
|
|
|
|
Таблица 21. |
|
|
Вариант плана 23 (последовательность №3) |
|||
|
|
|
|
|
|
№ |
|
|
Фактор |
|
|
п/п |
X4 |
|
X5 |
X6 |
|
