Основные понятия и планирование эксперимента (стр. 6 из 7)

Например, для ОЦКП при числе факторов n=3 имеем следующие параметры плана

,

,

.

Сам план принимает вид

Очевидно, что план является ортогональным. В отличие от планов ПФЭ для ОЦКП сумма квадратов факторов разных столбцов не является одинаковой.

По результатам опытов плана формируется полином

.

Коэффициенты полинома

определяется как

.

Можно преобразовать полином к виду

,

где

.

Значения параметров ОЦКП при числе факторов n

При n =2 ОЦКП совпадает с планом ПФЭ 23. Звездные точки ОЦКП в этом случае лежат на границах варьирования факторов. Если точки плана ПФЭ 2n всегда лежат на окружности (поверхности шара, гипершара), то точки плана ОЦКП не лежат на какой-либо одной окружности (поверхности шара, гипершара). План ОЦКП не является насыщенным. Так, например, для n = 3 полином имеет одиннадцать членов со своими коэффициентами, но для их определения используются пятнадцать опытов.

Пример плана ОЦКП для n = 2.

Параметры плана N0=4, N=9,

= 1, а = 2/3, 1-а=1/3, -а=-2/3, .Использован рассмотренный ранее план ПФЭ 22 с добавленными опытами 5-9.

Коэффициенты полинома составляют

;

;

;

;

;

.

Полином принимает вид

(Ранее по плану ПФЭ 22 был сформирован полином

Рассчитанные значения

по полиному приведены в плане. Также приведены величины , подтверждающие достаточно высокую точность полинома. Так в центральной точке плана, в отличие от случая применения плана ПФЭ 22, расхождений нет.

Лекция 8. Рототабельные планы

Рототабельные планы – это планы, у которых точки плана располагаются на окружностях (сферах, гиперсферах). У рототабельного плана первого порядка точки плана располагаются на одной окружности (сфере, гиперсфере) с радиусом R

,

где V=1,…, N - номер точки плана, i =1,…, n – номер фактора.

В таком случае точность оценивания функции отклика по любому направлению факторного пространства (для всех точек плана) одинаковая.

Рототабельный план может быть симметричным, когда точки плана располагаются симметрично друг друга. Рассмотренный ранее план ПФЭ 2n – рототабельный симметричный план первого порядка.

У рототабельных планов второго порядка точки плана располагаются на двух концентрических гиперсферах с радиусами R1 и R2 . В таких планах

,

для V =1,…, N0 и

,

для W=1,…, n0,

где V и W – текущие номера точек плана в двух подмножествах опытов N0 и n0 из их общего количества N, относящихся к двум разным концентрическим сферам. Одна из сфер может быть вырожденной, когда R2=0. Рассмотренный ранее ортогональный центрально-композиционный план второго порядка (ОЦКП) не является рототабельным планом, так как его точки лежат на трех концентрических окружностях (сферах, гиперсферах). При n=2 это очевидно из рис. 14. “Звездные” точки плана и точки плана ПФЭ 2n лежат на разных окружностях.

Рис. 14. Расположение точек ОЦКП на трех окружностях

Рототабельный план может быть ортогональным, если выполняется условие

,

где

, , , - номера столбцов плана.

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

Рототабельный ортогональный центрально-композиционный план (РОЦКП) строится аналогично рассмотренному ранее ОЦКП. К использованному в качестве ядра плану ПФЭ 2n добавляются “звездные” точки - по две на каждый фактор и несколько точек в центре плана. “Звездные” точки должны располагаться на поверхности гиперсферы с радиусом R, на которой лежат и точки плана ПФЭ 2n, то есть величина плеча “звездных” точек

должна равняться радиусу R. Это может быть обеспечено, при выполнении условия ортогональности, только при соответствующем выборе числа наблюдений в центральной (нулевой) точке плана n0. Для РОЦКП n0 зависит от числа факторов n. Напомним, что в ОЦКП n0 = 1 для любого числа n.

Радиус сферы, на которой лежат точки плана ПФЭ 2n при двух уровнях варьирования факторов с диапазоном

1 составляет (рис. 15)

Рис. 15. Радиус окружности (сферы), на которой лежат точки плана ПФЭ 2n при диапазоне варьирования факторов от –1 до +1:

а) - n=1,

;

б) - n=2,

;

в) - n=3,

.

Таким образом, при построении РОЦКП с ядром из плана ПФЭ 2n плечо “звездных” точек определяется числом факторов

.

Раннее при определении параметров ортогонального композиционного плана второго порядка с ядром из плана ПФЭ 2n было получено

,

где

- число точек плана ПФЭ,

- полное число точек композиционного плана второго порядка,

- константа преобразования элементов столбцов, соответствующих квадратам факторов.

В этом случае для РОЦКП число наблюдений в центре плана

.

Если n0 не целое, то при практическом построении плана его округляют до целого, но свойство ортогональности плана нарушается.

Параметры РОЦКП в зависимости от числа факторов

В [1] без вывода для РОЦКП рекомендуется принимать

.

Тогда

.

Параметры РОЦКП по [1]

Пример рототабельного ортогонального центрально-композиционного плана для n = 2.

Параметры плана:

Нет необходимости проводить восемь раз (точки с 9 по 16) опыты в центре плана. Достаточно провести этот опыт один раз и записать результат во все восемь строк. Строки сокращать нельзя, так как нарушается свойство ортогональности, и коэффициенты полинома будут определены неверно.

Коэффициенты квадратичного полинома рассчитаются, как и ранее.

Использован рассмотренный ранее план ПФЭ 22 с добавленными опытами 5-16.

,

, ,

,