3.8.        ОРТОГОНАЛЬНЫЕ центральные композиционные планы

При использовании ортогональных центральных композиционных планов (ОЦКП) можно получить все коэффициенты квадратичных полиномов (см п.3.2.). Планы второго порядка предполагают исследование не в трех точках (+1, -1, 0), а в пяти (-α, -1, 0, +1, + α); здесь ± α – «звездные точки», т.е. диапазон варьирования факторов разбивается на четыре участка.

Соответственно «+α» обозначает максимальное значение фактора; «-α» – минимальное значение фактора; «+1», «0», «-1» – точки внутри диапазона варьирования.

Кодирование факторов производится по формуле (3.1).

Значения «звездного» плеча (α) для разного количества факторов (К) составляют:

К

2

3

4

5

α

1,0

1,215

1,414

1,547

Для полных планов второго порядка при ином значении К, а также для дробных планов второго порядка значения «звездного» плеча (α) определяется из формулы:

.

Первоначальное преобразование квадратичных членов производится по формуле:

,

где  – сумма i-го столбца матрицы.

Квадратичный член имеет вид:

.

Подсчет коэффициентов (bi) производится по формуле:

.

Первоначально полином ищется в виде:

.

Для перехода к обычному виду (3.5) осуществляется преобразование свободного члена:

.

В табл. 3.14 дан ОЦКП для трех факторов. Число опытов для ОЦКП второго порядка подсчитывается по формуле:

.

Оценка значимости полученных коэффициентов (bi) производится согласно t-критерию (Стьюдента) по формуле:

,                                                   (3.6)

где  – среднеквадратическое отклонение;  – дисперсия коэффициентов.

В свою очередь, определяется по формуле:

,

где  – дисперсия опыта, равная

.

Таблица 3.14

Элементы плана эксперимента

х0

х1

х2

х3

х4= х12 – 0,73

х5 = 22 2- 0,73

х6 = х3 2 -0,73

х7 = х1х2

х 8= х1х3

х9 = х2х3

у

План полного факторного эксперимента типа 23

1

+1

-1

-1

-1

0,27

0,27

0,27

+1

+1

+1

у1

2

+1

+1

-1

-1

0,27

0,27

0,27

-1

-1

+1

у2

3

+1

-1

+1

-1

0,27

0,27

0,27

-1

+1

-1

у3

4

+1

+1

+1

-1

0,27

0,27

0,27

+1

-1

-1

у4

5

+1

-1

-1

+1

0,27

0,27

0,27

+1

-1

-1

у5

6

+1

+1

-1

+1

0,27

0,27

0,27

+1

+1

-1

у6

7

+1

-1

+1

+1

0,27

0,27

0,27

-1

-1

+1

у7

8

+1

+1

+1

+1

0,27

0,27

0,27

+1

+1

+1

у8

«Звездные» точки

9

+1

-1,215

0

0

0,75

-0,73

-0,73

0

0

0

у9

10

+1

+1,215

0

0

0,75

-0,73

-0,73

0

0

0

у10

11

+1

0

-1,215

0

-0,73

0,75

-0,73

0

0

0

у11

12

+1

0

+1,215

0

-0,73

0,75

-0,73

0

0

0

у12

13

+1

0

0

-1,215

-0,73

-0,73

0,75

0

0

0

у13

14

+1

0

0

+1,215

-0,73

-0,73

0,75

0

0

0

у14

Нулевая точка

15

+1

0

0

0

-0,73

-0,73

-0,73

0

0

0

у15

Здесь l – число параллельных опытов, проводимых в нулевой точке плана; - среднее значение .

В случае выполнения неравенства (3.6) коэффициент bi считается значимым, в случае невыполнения коэффициент bi считается незначимым и отбрасывается. Далее проверяется неравенство:

,                                                   (3.7)

где  – табличное значение F-критерия;  – дисперсия адекватности.

Дисперсия адекватности определяется по формуле:

,

где  – наблюденное (определенное из опыта) значение выходной величины;  – предсказанное (рассчитанное по полиному) значение выходной величины; - число степеней свободы (m – число определяемых коэффициентов bi).

В случае выполнения неравенства (3.7) математическая модель процесса считается адекватной при определенном уровне достоверности.