Разделы

Авто
Бизнес
Болезни
Дом
Защита
Здоровье
Интернет
Компьютеры
Медицина
Науки
Обучение
Общество
Питание
Политика
Производство
Промышленность
Спорт
Техника
Экономика

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

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

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

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

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

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

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

n = 2. r2 = 2 1/2, r3 = 2 2/4 = 21/2;

n = 8 и p = 2. r2 = 8 1/2 = 2 3/2, r3 = 2 (8 – 2)/4 = 23/2.

Пример 5.3.1.Построить матрицу ротатабельного ЦКП Бокса второго порядка для трех факторов.

Решение.

Ядром плана является ПФЭ вида 23 (радиус соответствующей гиперсферы r2 = 31/2 = 1,732). Звездные точки располагаются на гиперсфере с радиусом r3 = 23/4 = 1,682 и имеют координаты ( 1,682; 0; 0), (0; 1,682; 0), (0; 0; 1,682). Матрица планирования включает три гиперсферы и соответствует табл. 5.3, в которой g = 1,682. План содержит 15 точек и является ненасыщенным – количество оцениваемых коэффициентов 10.

В табл. 5.5 приведены минимально необходимые сведения для составления рассмотренного вида ротатабельных ЦКП.

Таблица 5.5

Количество факторов Число точек ПФЭ Число звездных точек Значение g
1,414
1,682
2,000
2,378
5, полуреплика 2,000
2,828
6, полуреплика 2,378
3,364
7, полуреплика 2,828

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

;

,

 

 

;

;

;

;

;

.

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

cov(b0, bii) = – 2s2(ỹ) l4 A/N ;

cov(bii, bij) = s2 (ỹ) (1–l4 )A/N.

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

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

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

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

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

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

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

Планы типа Вn представляют собой симметричные планы второго порядка с ядром в виде ПФЭ 2k или ДФЭ 2kp, дополненные 2k звездными точками с плечом g =1 и опытами в центре плана. Иначе говоря, эти планысостоят из 2k (2kp) вершин k-мерного гиперкуба с координатами ±1, из 2k центров (n–1)-мерных граней и некоторого числа опытов в центре гиперкуба. Количество точек плана с ядром из ПФЭ составляет N = 2k + 2k +1, для ДФЭ N = 2kp + 2k +1. В каждой точке проводится равное число опытов. Планы этого типа имеют минимально количество уровней варьирования факторов, равное трем, что позволяет более точно выдерживать режимы работы изделий при натурных испытаниях по сравнению с планами, в которых требуется большее число уровней изменения управляемых переменных. Планы типа Вn близки к D- и G-оптимальным планам.

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

;

;

;

;

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

Дата публикации:2014-01-23

Просмотров:628

Вернуться в оглавление:

Комментария пока нет...


Имя* (по-русски):
Почта* (e-mail):Не публикуется
Ответить (до 1000 символов):







 

2012-2019 lekcion.ru. За поставленную ссылку спасибо.