16
Ротатабелный план - это план, обеспечивающий равноточность исследования по всем направлениям пространства факторов.
А - оптимальный план - это план, обеспечивающий минимальное значение средней дисперсии оценок параметров модели.
Д - оптимальный план – обеспечивает минимальный объем эллипсоида рассеяния оценок параметра.
G - оптимальный план минимизирует максимальные значения дисперсии предсказываемых моделью значений параметра оптимизации. Для G -оптимальных планов в области планирования эксперимента не окажется точек, в которых точность предсказания значений параметра оптимизации будет хуже, чем в точках плана эксперимента.
Всем этим критериям удовлетворяют полные факторные планы (ПФП) и дробные факторные планы (ДФП).
ПФП – план, содержащий все возможные комбинации всех факторов на определенном числе уровней равное число раз.
ДФП – план, содержащий часть комбинаций ПФП.
План эксперимента первого порядка - план с двумя или более уровнями факторов, позволяющий найти раздельные оценки параметров регрессионной модели первого порядка (полинома первого порядка).
Наиболее простыми, но тем не менее достаточно эффективными, являются двухуровневые планы, в которых каждый фактор может находиться только на двух уровнях: верхнем и нижнем.
Двухуровневый ПФП содержит общее количество опытов, выражаемое формулой:
N 2K , |
(8) |
где к – количество факторов.
Этот план пригоден для построения линейчатых моделей вида:
k |
k k |
k k k |
|
y b0 x0 bi xi bij xi x j ... bijk xi x j xn ... |
(9) |
||
i 1 |
i 1 j i 1 |
i 1 j i 1 n j 1 |
|
Количество неизвестных коэффициентов этого уравнения регрессии также определяется по формуле (8). Следовательно, двухуровневый ПФП может быть использован для нахождения всех коэффициентов уравнения (9).
17
Избыточность (насыщенность) плана представляет собой разность между количеством проводимых опытов и количеством определяемых коэффициентов регрессии модели.
Для регрессионных моделей (9), содержащих смешанные произведения факторов хi, xj, xn, …, избыточность ПФП И определяется по формуле:
И 2k (k 1) m , |
(10) |
где m – количество в регрессионной модели членов со смешанными произведениями факторов.
В случае линейной модели неизвестно к+1 коэффициентов. Следовательно, для линейной модели ПФП является избыточным на величину 2к – (к+1).
Избыточные опыты плана могут быть использованы либо для уменьшения случайных погрешностей линейной модели, либо для уменьшения погрешности адекватности модели путем включения в нее произведений факторов: (aij xi x j , aijn xi x j xn , ...), либо и для того и другого вместе.
Планы, у которых количество опытов и количество коэффициентов регрессии равны друг другу называются насыщенными планами.
Методика построения ПФП и ДФП
Для унификации планов факторы в них записываются в безразмерном
нормализованном виде. Для этого вводятся величины: |
|
|||||||||
1) |
Центр эксперимента для данного фактора: |
|
|
|||||||
|
X |
|
|
xmax |
xmin |
. |
|
(11) |
||
|
0 |
|
|
2 |
|
|||||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
2) |
Шаг варьирования данного фактора: |
|
|
|
|
|
||||
|
X |
xmax |
xmin |
, |
(12) |
|||||
|
|
2 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
Тогда безразмерная нормализированная величина фактора х определяется по формуле:
X |
x x0 |
. |
(13) |
|
|||
|
X |
|
|
18
Нетрудно увидеть, что при x=xmin, Xmin= –1; при x=x0, X0=0; при x=xmax, – Xmax=1.Следовательно, в таблицах активного эксперимента используются уровни –1 ,0 , +1, в которых цифру "1" опускают.
Рассмотрим пример, когда количество факторов к=3. В этом случае Non 23 8 , а план эксперимента имеет вид, представленный в табл. 3 [4].
Строки плана представляют собой условия проведения опытов при соответствующих фиксированных значениях уровней факторов Х1, Х2, Х3.
В таблицу вводится столбец Х 0, который всегда равен +1.Это так называемый постоянный фактор, который является сомножителем свободного члена уравнения регрессии (y=b0∙X0+…), где Х0=1.
Правило заполнения столбцов ПФП
Заполнение начинают с левого столбца (X1) сверху вниз чередованием знаков "+" и "–". Следующий столбец заполняют чередованием по 2 знака "+" и "–". Третий столбец – по 4 знака "–" и 4 знака “+” Затем (если бы у нас была n-факторная модель) по 8 знаков "–" и 8 знаков "+", затем по 16 и т.д. Столбец yi – отклик.
Данный ПФП имеет избыточность, равную 23–(3+1)=4. Поэтому для нахождения параметров достаточно использовать половину из опытов ПФП.
Правило заполнения столбцов ДФП
Таблица 4 – Полный 3-факторый план эксперимента
Non |
Х 0 |
Х 1 |
Х 2 |
Х 3 |
yi |
1 |
+ |
– |
– |
– |
y1 |
2 |
+ |
+ |
– |
– |
y2 |
|
|
|
|
|
|
3 |
+ |
– |
+ |
– |
y3 |
4 |
+ |
+ |
+ |
– |
y4 |
|
|
|
|
|
|
5 |
+ |
– |
– |
+ |
y5 |
6 |
+ |
+ |
– |
+ |
y6 |
|
|
|
|
|
|
7 |
+ |
– |
+ |
+ |
y7 |
8 |
+ |
+ |
+ |
+ |
y8 |
|
|
|
|
|
|
Заголовки столбцов здесь те же, а количество опытов – 4 (табл. 4).
Сначала процедура заполнения такая же, как и для случая ПФП. Она проводится до того столбца, где по правилу должны быть поставлены только знаки "–" и чередование знаков отсутствовало бы (в данном случае столбец X3 должен был бы состоять из знаков "–"). В этом и последующих столбцах знаки находятся по правилу перемножения уже имеющихся знаков
данной строки (в данном случае X3= X2∙∙X1).
ДФП называют также дробной репликой ПФП.
Если дробление ПФП производится многократным последовательным делением числа его опытов на 2, то реплика называется регулярной. Число раз "р"
19
такого последовательного деления называется дробностью реплики. При Р =1 ДФП называют полурепликой, при Р =2 –четвертьрепликой, при Р =3 – 1/8 реплики, при Р =4 – 1/16 реплики и т.д.
Таблица 5 Пример дробного факторного плана
Non |
X0 |
X1 |
X2 |
X3 |
yi |
|
|
|
|
(Х3=Х1∙Х2) |
|
1 |
+ |
– |
– |
+ |
y1 |
|
|
|
|
|
|
2 |
+ |
+ |
– |
– |
y2 |
|
|
|
|
|
|
3 |
+ |
– |
+ |
– |
y3 |
4 |
+ |
+ |
+ |
+ |
y4 |
В этом случае расчет количества опытов ДФП осуществляется по формуле
N 2k p , |
(14) |
а избыточность плана определяется по формуле
И 2k p (k 1) m . (15)
Методика построения плана эксперимента второго порядка
Как указывалось ранее, уравнение регрессии 1-го порядка зачастую не является адекватным и для обеспечения адекватности модели к исходному уравнению 1-го порядка добавляют члены со сменными произведениями факторов, типа bi j xi x j , bij xi x j xn и т.д., где bi j , bijn – коэффициенты регрессии смешанных произведений факторов xi , x j , xk и т.д.
Для нахождения этих коэффициентов регрессии используется избыточные опыты того же плана эксперимента 1-го порядка, который использовался в случае уравнения 1-го порядка.
Однако практика показывает, что во многих случаях такое добавление членов
k
смешенных произведений факторов к правой части линейной модели y b0 bi xi
i 1
и не приводит к желаемому результату обеспечения адекватности модели. В этом случае либо уменьшают интервалы варьирования одного из факторов (или, может быть, всех факторов), либо переходят к моделям 2-го, затем 3-го и т.д. порядков, в которые вводят члены, содержащие факторы во 2-й, 3-й и т.д. степенях.
Рассмотрим построение оптимального плана эксперимента 2-го порядка на примере двухфакторной модели:
y f x , x |
2 |
b |
b |
x |
b |
x |
2 |
b |
x2 |
b |
x |
x |
2 |
b |
x2 . |
(16) |
1 |
0 |
1 |
1 |
2 |
|
11 |
1 |
12 |
1 |
|
22 |
2 |
|
20
Функция отклика у этой модели представляет собой поверхность 2-го порядка: цилиндр, конус, эллипсоид, параболоид, гиперболоид, параболический цилиндр, параболический гиперболоид или сферу.
Сечения этой поверхности, параллельные каждой из осей X1 и X2, в общем случае представляют собой кривые также второго порядка (окружности, эллипсы, параболы, гиперболы). Известно, что для однозначного определения такой кривой необходимо располагать координатами не менее 3-х ее точек, т.е. факторы X1 и X2 должны варьироваться не менее чем на 3-х уровнях.
Для удовлетворения этого требования план эксперимента 2-го порядка помимо опытов плана эксперимента 1-го порядка, называемого ядром плана 2-го порядка, должен включать опыты в звездных точках и в центре плана [4]
Центром плана, или центральной точкой плана, называется точка плана соответствующая нулевым нормализованным безразмерным координатам факторов.
Звездной точкой плана эксперимента 2-го порядка называется точка, лежащая на координатной оси факторного пространства.
Звездная точка удалена от центральной точки на расстояние звездного плеча. Для центральных композиционных ротатабельных планов звездное плечо
рассчитывается по формуле: |
|
|
|
|
|
|
|
|
|
||||
Таблица 13 |
|
|
|
|
|
|
|
|
|
|
n p |
|
|
План эксперимента 2-го порядка |
|
|
|
|
|
2 2 . |
|||||||
|
|
|
|
|
|
|
|
|
|
|
(17) |
||
Non |
X0 |
X1 |
X2 |
|
yi |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||||
1 |
+ |
+ |
+ |
|
y1 |
|
|
|
|
|
|
|
|
2 |
+ |
– |
+ |
|
y2 |
Опыты в |
Общее количество звездных |
||||||
3 |
+ |
+ |
– |
|
y3 |
ядре |
точек, |
расположенных |
|
на |
|||
4 |
+ |
– |
– |
|
y4 |
плана |
|
||||||
|
|
|
от центра плана |
||||||||||
5 |
+ |
0 |
|
|
y5 |
|
расстояниях |
||||||
6 |
+ |
0 |
- |
|
y6 |
Опыты в |
равно |
2k. |
Все координаты |
||||
7 |
+ |
|
0 |
|
y7 |
звездны |
|||||||
|
звездной |
точки |
за исключением |
||||||||||
8 |
+ |
- |
0 |
|
y8 |
х |
|||||||
|
|
|
|
|
|
точках |
одной равны |
нулю. Например, |
|||||
7 |
+ |
0 |
0 |
|
y9 |
|
|||||||
|
|
координата |
звездной |
точки, |
|||||||||
8 |
+ |
0 |
0 |
|
y10 |
Опыты в |
|||||||
|
|
|
|
|
|
|
|
||||||
9 |
+ |
0 |
0 |
|
y11 |
центре |
лежащей |
|
на |
положительной |
|||
10 |
+ |
0 |
0 |
|
y12 |
плана |
полуоси Х1, записываются так: (α, |
||||||
11 |
+ |
0 |
0 |
|
y13 |
|
|||||||
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
0, 0, 0, … 0). |
|
|
|
|
|
|
|
Для 2-х факторного плана эксперимента с ядром ПФП |
|
|
|
|
|
|||||||
|
|
|
|
|
|
2 0 |
|
|
|
|
|
|
|
2 4 
2 .