Дробный факторный эксперимент

С увеличением числа факторов количество опытов при ПФЭ растет по показательной функции. В то же время часто заранее известно, что число определяемых коэффициентов в модели объекта намного меньше числа опытов в соответствующем ПФЭ (например, если можно ограничиться линейным приближением модели). В этом случае используя дробный факторный эксперимент (ДФЭ) и получая дробные реплики можно резко сократить число опытов.

ДФЭ позволяет минимизировать число опытов, но при этом оценки коэффициентов получаются смешанными.

Эффективность ДФЭ при определении модели объекта зависит от выбранной системы смешивания. Необходимо выбрать такие дробные реплики, при которых интересующие экспериментатора линейные эффекты смешиваются с теми эффектами взаимодействия, о которых с наибольшей вероятностью можно предположить, что они близки к нулю. Как правило, такое предположение чаще всего справедливо относительно взаимодействий самого высокого порядка.

Для определения системы смешивания оценок используется так называемый определяющий контраст — символическое произведение всех столбцов матрицы планирования равное +1 или -1. Чтобы узнать какие взаимодействия смешаны с данным линейным эффектом (то есть в каких столбцах плана имеются одинаковые элементы) необходимо умножить определяющий контраст на столбец соответствующий данному эффекту.

Положим, что контраст в случае трех факторов равен

. (6.49)

Умножим по очереди определяющий контраст на , получим

(6.50)

поскольку .

Полученным соотношениям (6.50) соответствует система смешанных оценок:

(6.51)

При увеличении числа факторов целесообразность использования дробных реплик для получения линейной модели объекта резко возрастает.

Рассмотрим пример ДФЭ. Положим, что интересующая нас модель объекта является функцией трех переменных (факторов), причем заранее известно, что в рабочей области эта модель близка к линейной, то есть

(6.52)

При ПФЭ и двух уровнях факторов мы должны осуществить опытов. Следовательно, для определения четырех коэффициентов модели проводится вдвое больше число опытов.

Используя ДФЭ для решения поставленной задачи можно уменьшить число опытов. Чтобы дробная реплика представляла собой ортогональный план, в качестве реплики следует брать ПФЭ для меньшего числа факторов. Число опытов при этом должно быть больше (или равно) числа неизвестных коэффициентов модели. Следовательно, для рассматриваемого примера можно ограничиться четырьмя опытами. При этом ортогональность плана достигается если использовать план ПФЭ типа , а в качестве плана для использовать столбец . Матрица такого эксперимента приведена в таблице 3.6.

Таблица 3.6

Полуреплика от ПФЭ типа 23

№ опыта
	+	+	+	+
	+	—	+	—
	+	+	—	—
	+	—	—	+

Такой сокращенный план (половина ПФЭ) называется полурепликой. Пользуясь таким планированием, можно оценить свободный член и три коэффициента модели при линейных членах. Но если модель объекта не строго линейна и коэффициенты при парных взаимодействиях не равны нулю, то найденные коэффициенты будут смешанными в соответствии с системой (3.50), поскольку элементы столбцов , и (эти столбцы в табл. 3.6 не указаны) совпадают с элементами столбцов соответственно. Если модель объекта действительно близка к линейной, то со смешанностью коэффициентов вполне можно примириться.

По указанному в табл. 3.6 плану проводится половина ПФЭ типа и получается полуреплика от , что принято обозначать как . Если знаки одного из элементов столбцов (, или ) изменить на противоположные, то получим вторую половину ПФЭ (еще 4 опыта) и, следовательно, вторую полуреплику от . При осуществлении обоих полуреплик мы провели бы ПФЭ (8 опытов) и получили бы раздельные оценки коэффициентов как для линейных эффектов, так и для эффектов взаимодействий.

Пусть , , тогда план ДФЭ включает матрицу плана ПФЭ для трех факторов , и , а столбцы плана для остальных факторов представляют собой некоторые произведения факторов , и , например, , , , . Такие равенства называют генерируемыми соотношениями, их выбор произволен, однако повторение не допускается. При таком планировании все линейные эффекты смешаны с эффектами взаимодействия.

Число опытов в дробной реплике должно удовлетворять неравенству , где — число факторов, для получения несмешанных оценок линейных эффектов.

Если , то есть число опытов равно числу определяемых коэффициентов линейной модели, дробная реплика представляет собой насыщенный линейный ортогональный план.

Поскольку число опытов в насыщенных планах равно числу определяемых коэффициентов, число степеней свободы дисперсии адекватности . Для поверки адекватности линейной модели, полученной по насыщенному плану, необходим дополнительный эксперимент.

Таким образом оптимальные двухуровневые планы и имеют следующие преимущества:

 планы ортогональны, и поэтому все вычисления просты, все коэффициенты определяются независимо друг от друга;

 каждый коэффициент определяется по результатам всех опытов;

 для данного числа опытов эти планы имеют минимальный определитель ковариационной матрицы ; в связи с этим все коэффициенты определяются с одинаковой и минимальной дисперсией, то есть планы обладают свойством Д-оптимальности;

 линейные планы и обладают свойством ротатабельности.

Покажем справедливость свойства ротатабельности в случае линейной модели

. (6.53)

Вследствие отсутствия корреляции между коэффициентами по закону сложения дисперсий для линейной модели (6.53) имеем:

.

Так как , получим

;

,

где — радиус сферы в К-мерном пространстве.

Рис. 3.3.

Дисперсия выхода растет пропорционально квадрату радиуса сферы и одинаково для всех эквидистантных точек (рис. 3.3).

Планирование обладающее таким свойством называется ротатабельным планированием.