logo search
ПОСОБИЕ по численным для издания

4.10 Планы для построения полиномиальных моделей второго порядка

Полином второго порядка от n переменных имеет вид:

Y = b0 + biXi +biiXi2 + bijXiXj (4.65)

Существует большое количество планов для построения данного вида регрессии, наиболее лучшие из них представлены в «Таблицах планов экспериментов» [30].

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

(4.66)

Все остальные называются несимметричными.

Расчет коэффициентов регрессии для симметричных планов осуществляется по ковариационно-корреляционной матрице, которую в общем виде можно представить следующим образом:

(4.67)

где a*; N*; r*;a-1; p-1; p* - элементы матрицы.

Коэффициенты регрессии находят по упрощенной схеме:

(4.68)

Коэффициенты регрессии несимметричных планов находят по матрице для расчета м.н.к.-оценок параметров моделей, которая в общем виде выглядит следующим образом и обозначается L (таблица 7).

Таблица 7 - Матрица для расчета м.н.к. - оценок параметров моделей

Θ 0 Θ ii Θ jj Θ i Θ j Θ ij Θ kl

Θ0

Θii

Θjj

Θi

Θj

Θij

Θi0 Θ0ii Θ0jj Θ0i Θ0j Θ0ij Θ0kl

Θii0 Θiiii Θiijj Θiii Θiij Θiiij Θiikl

… … … … … … …

Коэффициенты регрессии вычисляются по формуле:

(4.69)

Это означает, что оценку каждого параметра получают как скалярное произведение соответствующей строки матрицы L на вектор наблюдений У, который в общем виде может быть представлен следующим образом (таблица 8):

Таблица 8 – План и результаты эксперимента

№ опыта

Факторы

У

Х1

Х2

1

+1

-1

У1

2

-1

0

У2

...

n

+1

0

Уn

Тогда, например, свободный коэффициент регрессии будет вычисляться следующим образом:

В000У1+ Ө0iiУ2+ Ө0jjУ3+ Ө0iУ4+ Ө0jУ5+ ӨijУ6+ Ө0кеУ7+…+ Ө0кеУn (4.70)

Квадратичный коэффициент:

Вii= Өii0У1+ ӨiiiiУ2+ ӨiijjУ3+ ӨiiiУ4+ ӨiijУ5+ ӨiiijУ6+…+ ӨiiкеУn (4.71)