Ротатабельные планы

Всякий план будет ротатабельным только тогда, когда экспериментальные точки будут лежать не менее чем на двух сферах (в общем случае, для k-мерного пространства,—на двух гиперсферах).

В качестве первой сферы с нулевым радиусом выбирают центральную точку плана — центр эксперимента. Тогда все остальные опыты будут располагаться на второй сфере выбранного радиуса  и будут равно удалены от нулевой точки.

В частном случае, для плана второго порядка при k=2 (плоская задача) вместо сфер можно рассматривать две окружности. Графическое построение такого плана представлено на рис. 1.

Вершины квадрата здесь соответствуют координатам экспе­риментальных точек полного факторного эксперимента типа 2^k, если вокруг данного квадрата описать окружность с центром в нулевой точке, то ее пересечение с центральными осями даст координаты четырех дополнительных звездных точек (±, 0 и 0±). Таким образом, все экспериментальные точки лежат на двух окружностях. Первая - в центре эксперимента с нулевым радиусом, а вторая с радиусом = (величина звездного плеча). Из геометрических соотношений легко определить значение  для данного случая. Как видно, ²=1²+1² и =. В общем случае при числе переменных равном k величина звездного плеча определяется соотношениями ^k/4.

Количество опытов в нулевой точке назначается из условия получения равной информации о процессе как в центре экспе­римента, так к в периферийных точках и зависит от числа вклю­ченных в рассмотрение факторов Такой план называется униформ - планом. Значения величии n_о для ротатабельного планирования табулированы и приведены в литературе (Ю.П. Адлер Введение в планирование эксперимента, М., «Металлургия», 1969; В.В. Налимов и Н.А. Чернова Статистические методы планирования экстремальных экспериментов). Рассмотренный метод построения определяет центральное композиционное ротатабельное униформ планирование второго порядка.

Матрица такого планирования представлена в табл. . В нее включен без изменения ПФЭ типа 2³ , т. е. основные (с 1 по 8) и дополнительные (с 9 по 14) опыты с звезд­ным плечом ^k/43/4=1,682 и опыты (с 15 по 20) в нулевой точке. Для реализации такого плана, используя весь предыду­щий материал, достаточно провести лишь шесть опытов в звезд­ных точках.

Обработка полученных данных для расчета коэффициентов полного квадратного уравнения ведется по сравнительно сложным формулам. Для нашего случая при к=3 эти формулы имеют вид:

b₀=0,166338(0y) - 0,056791 ;

b₁ = 0,073224 (iy)

bii = 0,0625(iiy)+0,006889- 0,056791 (Oy),

bij = 0,125(ijy),

(Oy),=,

(iiy)= ,

(iy)= ,

(ijy),= ,

i, j - номера столбцов матрицы планирования;

u - номер строки,

N - количество опытов в плане.

Статистическй анализ результатов для планов второго прядка

При проведении статистического анализа полученного уравнения регрессии воспользуемся рассмотренной ранее методической схемой для ПФЭ.

Оценка дисперсии, характеризующей ошибку опыта. Определение ошибки опыта проводилось на каждом этапе исследований: для первой и второй полуреплик, для полного факторного эксперимента, для опытов в нулевой точке. Раннее было подтверждено, что все вычисленные дисперсии однородны, т. е. практически равны. Следовательно, для дальнейшего статистического анализа можно использовать любую из рассчитанных дисперсий: построчную S_i², усредненнyю S_ср² дисперсию в нулевой точке S₀² и т. д.

Выбор той или иной оценки в некоторой степени произволен. Однако очевидно, что при ротатабельном планировании эксперимента удобнее оценивать дисперсию по результатам опытов в нулевой точке, так как дублирование этих опытов предусмотрено самим планом эксперимента (опыты с 15 по 20 в матрице планирования). В качестве оценки ошибки опыта примем дисперсию в нулевой точке S_y² f_y=n_o— 1

Определение значимости коэффициентов регрессии. Расчет дисперсий при определении коэффициентов производится до следующим формулам:

S²(b₀ )=

S²(b_i )=

S²(b_ii )=

S²(b_ij )=

где С=; А= ;

.

Значимость коэффициентов определяется с помощью t- критерия Стьюдента. Если коэффициент статистически значим, то влияние соответствующего фактора существенно. Однако обратное утверждение не всегда справедливо. Поэтому при окончательной записи уравнения нельзя исключать значимые коэффициенты и в то же врем по чисто технологическим соображениям допускается введения любого статистически незначимого коэффициента.

При проверке адекватности полученной модели, как и для линейный планов, прежде всего вычисляется остаточная сума квадратов S_R

S_R=

Далее рассчитывается сумма квадратов, связанную с повторением опытов в нулевой точке S_E .

S_Е=

Для расчета дисперсии адекватности необходимо определить число степеней свободы системы.

f_ад=N - m - ( n_o - 1 ),

где N- общее число опытов в плане,

m - число коэффициентов, включенных в модель,

n_о - число опытов в нулевой точке.

Адекватность проверяем по F - критерию Фишера

при f₁ = f_ад , f₂ = f₀