Построение математических моделей исследуемых систем на основании экспериментальных данных, страница 14

          Если выбрать соответствующим образом значения  α, то план Бокса может быть сделан ортогональным.

          При соответствующем выборе звездных плеч план Бокса может быть сделан и рототабельным.

          4.2.4.ЦКП Хартли.

          Идея построения планов Хартли такая же как и у Бокса,  но допускается взаимодействие парных и линейных взаимодействий.

          Центральный композиционный план называется планом Хартли, если в качестве его ядра используется дробная реплика с разрешающей способностью III.

         

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

          Недостатки плана Хартли:

1)  смешение воздействий;

2)  они не могут быть сделаны ни ортогональными, ни рототабельными.

Сравнение планов Бокса и Хартли.

Число

факторов

ПБ

ПХ

ядро

число

отсчетов

ядро

число

отсчетов

3

ПФЭ 23

15

ДФЭ 

11

4

ПФЭ 24

25

ДФЭ 

17

5

 ДФЭ      

27

ДФЭ 

27

6

 ДФЭ     

45

ДФЭ 

29

7

 ДФЭ    

79

ДФЭ

47

4.2.5.Ортогональные планы Бокса.

          Для определения параметров квадратичной регрессионной зависимости

можно построить ортогональный план Бокса за счет:

а)  выбора величины звездного плеча;

б)  видоизменения уравнения регресии путем введения дополнительных переменных:     

          При этом уравнение регресии принимает вид:

         

          Коэффициент       при такой замене изменится, а все остальные коэффициенты будут неизменны.

                   

          При  р = 2 матрица плана имеет вид:

план

х0

х1

х2

х1х2

ядро

плана

ПФЭ 22

 

+

+

+

+

+

-

+

-

-

-

+

+

+

-

-

+

1-β

1-β

1-β

1-β

1-β

1-β

1-β

1-β

звезд-

ная

часть

+

+

+

+

α

0

0

0

0

α

0

0

0

0

α2

α2

α2

α2

центр

+

0

0

0

Аналогично строятся планы при  р>2, когда ядро плана определяется ПФЭ или ДФЭ. Коэффициенты  α и β определяются:

 


где N0 – число экспериментов, определяемое ядром плана, N – полное число экспериментов.

          Если, например, N0 = 4, то N = 9 и α = 1, β = 2/3.

          Для рассматриваемого случая матрица соответствующего плана будет иметь такие параметры:

х0

х1

х2

х1х2

+

+

+

+

-

+

-

+

-

-

+

+

+

-

-

+

1/3

1/3

1/3

1/3

1/3

1/3

1/3

1/3

+

+

+

+

-

+

0

0

0

0

-

+

0

0

0

0

1/3

1/3

-2/3

-2/3

-2/3

-2/3

1/3

1/3

+

0

0

0

-2/3

-2/3

          Матрица этого плана получилась ортогональной, поэтому дисперсионная матрица будет диагональной.

где Е – единичная матрица.

          В рассматриваемом нами примере  р = 2

 


Дисперсию оценок коэффициентов а можно определить следующим образом:

          Для определения оценок коэффициентов   ai, aij, aii,      можно использовать  известное соотношение:   

Истинное значение оценки коэффициентов а0:  

          Расчет доверительного интервала для оценок а осуществляется с использованием статистики Стьюдента (число степеней свободы равно 1).

          Для определения коэффициентов α и β используется таблица: