Ротатабельные центральные композиционные планы

В некоторых случаях ортогональное планирование второго порядка не отвечает потребностям практики – при описании поверхности отклика, особенно в окрестностях точки оптимума, более значимой является оценка дисперсии уравнения в целом, чем оценка дисперсии отдельных коэффициентов полинома. В этом случае обычно стремятся к равномерности распределения информации в уравнении функции отклика по всем направлениям. Такому положению отвечают ротатабельные планы. Подобные планы позволяют минимизировать систематические ошибки, связанные с неадекватностью представления результатов полиномами второго порядка. Но построение такого плана более сложно, чем ортогонального, а сама задача построения не имеет однозначного решения. Один из подходов к построению таких планов состоит в следующем.

Путем специального подбора звездного плеча g ЦКП Бокса можно сделать ротатабельным, иначе говоря, ЦКП Бокса можно сделать или ортогональным или ротатабельным.

Точки ротатабельного ЦКП Бокса второго порядка располагают на концентрических гиперсферах, количество которых не менее двух. Первая гиперсфера может быть вырожденной, т. е. представлять собой центральную точку плана, ее радиус r₁ = 0. Именно такая сфера часто используется на практике.

Вторая гиперсфера соответствует вписанному в нее кубу, выбранному в качестве ядра плана. Для ядра х_i = 1, следовательно, радиус этой гиперсферы r₂ = (х ₁² + х ₂² + … + х_k²)^1/2 = (k)^1/2. Ядро представляет собой ПФЭ вида 2 ^k или ДФЭ вида 2 ^{k – p} , причем должно соблюдаться условие (k – p)/4 > 3/4. Следовательно, с учетом ограничений на ЦКП Бокса, если k ³ 5, то в качестве ядра можно использовать полуреплику, если k ³ 8, ядром может служить четверть реплика.

Третья гиперсфера имеет радиус r₃ = 2 ^{k / 4} для ядра в виде ПФЭ и радиус r₃ = 2 ^{( k - p ) / 4} для ядра в виде ДФЭ.

Таким образом, каждый фактор в ротатабельном ЦКП Бокса варьируется на пяти уровнях. В некоторых случаях радиусы второй и третьей гиперсферы совпадают:

n = 2. r₂ = 2 ^1/2, r₃ = 2 ^2/4 = 2^1/2;

n = 8 и p = 2. r₂ = 8 ^1/2 = 2 ^3/2, r₃ = 2 ^{(8 – 2)/4} = 2^3/2.

Коэффициенты модели и их дисперсии рассчитываются по формулам:

;

,

;

;

;

;

;

.

Представленные формулы справедливы для ротатабельного планирования при любом количестве независимых переменных. Такое планирование не позволяет получить независимые оценки для всех коэффициентов модели, коррелированными оказываются коэффициенты (b₀, b _ii) и (b _ii, b _ij). Взаимную связь этих пар коэффициентов можно охарактеризовать ковариациями:

cov(b₀, b _ii) = – 2s²(ỹ) l₄ A/ N;

cov(b _ii, b _ij) = s² (ỹ) (1–l₄ )A/ N.

Проверка однородности дисперсии воспроизводимости, адекватности модели и значимости коэффициентов модели производится по схеме, рассмотренной в разделе 4.

Если повторные наблюдения имеются только в центре плана, то и величина будет несмещенной оценкой дисперсии ошибок наблюдения. При ненасыщенном планировании остаточная сумма S_R ² = отличается от нуля. Здесь – величина, предсказанная уравнением модели, – найденная экспериментально. Величина s _R ² =S_R / [ N –(k +1)(k +2)/2] характеризует неадекватность модели и также является несмещенной оценкой дисперсии ошибок наблюдения.

На основании рассчитанных величин можно провести все необходимые проверки коэффициентов и модели в целом.

Иногда интерес представляет информация о функции отклика в некоторой окрестности центра плана. В этом случае следует добиться одинаковой погрешности модели внутри гиперсферы единичного радиуса. План, обеспечивающий такое свойство функции отклика, называется униформ-ротатабельным. Для его формирования достаточно обеспечить равенство дисперсии в центре плана (r₀ = 0) и на поверхности гиперсферы радиуса r₂ = 1. Этого добиваются подбором числа наблюдений n ₀ в центре плана, а именно, параметр λ₄ следует взять равным положительному корню квадратного уравнения

2λ₄ (λ₄ – 1)(k + 2) + λ₄ (k + 1) – (k – 1) = 0.

Рассмотренное композиционное планирование представляет собой один из возможных подходов к построению ротатабельных планов второго порядка.

5.4. Композиционные планы типа В_n

Планы типа В_n представляют собой симметричные планы второго порядка с ядром в виде ПФЭ 2 ^k или ДФЭ 2 ^{k – p}, дополненные 2 k звездными точками с плечом g =1 и опытами в центре плана.

Иначе говоря, эти планы состоят из 2 ^k (2 ^{k – p}) вершин k -мерного гиперкуба с координатами ±1, из 2 k центров (n– 1)-мерных граней и некоторого числа опытов в центре гиперкуба.

Количество точек плана с ядром из ПФЭ составляет N = 2 ^k + 2 k +1, для ДФЭ N = 2 ^{k – p} + 2 k +1. В каждой точке проводится равное число опытов.

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

Планы типа В_n близки к D - и G -оптимальным планам.

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

;

;

;

;

где N ₁ – число точек ядра плана; – среднее значение отклика в u -й точке, полученное по r опытам. Если некоторые коэффициенты незначимы, то остальные уточняются по специальным формулам.

Дата добавления: 2014-01-11