Ортогональное центральное композиционное планирование


Будем преобразовывать модель следующим образом.

Ортогональное центральное композиционное планирование

clip_image004

т.е. clip_image006

Следует доказать, что при таком преобразовании ортогональность будет соблюдена.

clip_image008

Анализ преобразования в матрице планирования.

G

clip_image010

clip_image012

clip_image014

clip_image016

clip_image012[1]clip_image014[1]

clip_image012[2]clip_image016[1]

clip_image014[2]clip_image016[2]

clip_image019

clip_image021

clip_image023

П

Ф

э

clip_image025

1

+1

-1

-1

-1

+1

+1

+1

clip_image027

clip_image027[1]

clip_image027[2]

2

+1

+1

-1

-1

-1

-1

+1

clip_image027[3]

clip_image027[4]

clip_image027[5]

3

+1

-1

+1

-1

-1

+1

-1

clip_image027[6]

clip_image027[7]

clip_image027[8]

4

+1

+1

+1

-1

+1

-1

-1

clip_image027[9]

clip_image027[10]

clip_image027[11]

5

+1

-1

-1

+1

+1

+1

-1

clip_image027[12]

clip_image027[13]

clip_image027[14]

6

+1

+1

-1

+1

-1

-1

-1

clip_image027[15]

clip_image027[16]

clip_image027[17]

7

+1

-1

+1

+1

-1

+1

+1

clip_image027[18]

clip_image027[19]

clip_image027[20]

8

+1

+1

+1

+1

+1

-1

+1

clip_image027[21]

clip_image027[22]

clip_image027[23]

Звездные точки 2n

9

+1

—clip_image029

0

0

0

0

0

clip_image031

clip_image033

clip_image033[1]

10

+1

+clip_image029[1]

0

0

0

0

0

clip_image031[1]

clip_image033[2]

clip_image033[3]

11

+1

0

—clip_image029[2]

0

0

0

0

clip_image033[4]

clip_image031[2]

clip_image033[5]

12

+1

0

+clip_image029[3]

0

0

0

0

clip_image033[6]

clip_image031[3]

clip_image033[7]

13

+1

0

0

—clip_image029[4]

0

0

0

clip_image033[8]

clip_image033[9]

clip_image031[4]

14

+1

0

0

+clip_image029[5]

0

0

0

clip_image033[10]

clip_image033[11]

clip_image031[5]

(0.0)

15

+1

0

0

0

0

0

0

clip_image033[12]

clip_image033[13]

clip_image033[14]

Нужно найти величину clip_image029[6]. Для этого перемножим clip_image019[1]и clip_image021[1].

clip_image035

Подставим вместо clip_image037ее выражение.

clip_image039

Получаем биквадратное уравнение с неизвестным clip_image029[7]

После решения получаем clip_image041— плечо звёздной точки.

N

clip_image043

clip_image045

clip_image047

clip_image049

clip_image051

clip_image037[1]

1

1.00

4

4

1

9

0.73

2

1.215

8

6

1

15

4

1.414

16

8

1

25

Переходим к ДФЭ

clip_image054clip_image056clip_image054[1] clip_image059

Линейные члены будут иметь одну дисперсию, квадратные другую парные взаимодействия третью.

clip_image061Если получается дробное число, округляем до ближайшего целого.

Получаем модель в следующем виде

clip_image063

Получили модель, полную, квадратичную, но немного не такой свободный член:

clip_image065

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

В результате отсева можем нарваться на дробное число степеней свободы. Раскрывая clip_image056[1]— в модели получим модель, где все наши преобразования содержатся в свободном члене.

Проверка на адекватность — разные дисперсии у разных членов и число степеней свободы N-d.

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

Загрузка...