v0.5.7.final / 03
.pdf11 Тема 02: Построение эмпирических статистических моделей ХТП
При этом уравнение регрессии, описывающее эти опытные данные,
записывается с использованием кодированных факторов z |
( j 0,1,2) и, |
||||
|
|
|
j ~ |
~ ~ |
|
соответственно, кодированных коэффициентов регрессии |
a0 |
, a1, a2 |
: |
||
ˆ ~ |
~ |
~ |
|
|
|
y a0 |
a1z1 |
a2 z2 |
|
|
|
В кодированном факторном пространстве в соответствии с указанным планом
проведения эксперимента проведённые опыты представляются точками вершин |
|
квадрата: |
z2 |
|
|
( 1; 1) |
( 1; 1) |
z1
( 1; 1) |
( 1; 1) |
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
12 Тема 02: Построение эмпирических статистических моделей ХТП
Для параметрической идентификации кодированного уравнения регрессии используется метод регрессионного анализа, включающий три этапа:
- определение кодированных коэффициентов уравнения регрессии ~ a
методом наименьших квадратов;
- оценка значимости кодированных коэффициентов регрессии с использованием t – критерия Стьюдента;
- проверка адекватности кодированного уравнения регрессии с использованием F – критерия Фишера.
Реализация двух последних этапов возможна при выполнении свойства однородности дисперсий (одно из требований регрессионного анализа) и проведении параллельных опытов, например, в центре плана
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
13 Тема 02: Построение эмпирических статистических моделей ХТП
При проведении k параллельных опытов в центре плана ( yэксп , s 1, ..., k)среднее |
|
|
0s |
значение ycэксп определяется как среднее арифметическое результатов |
|
измерений во всех параллельных опытах: |
|
k |
|
ycэксп y0экспS |
k |
S 1 |
|
Определение кодированных коэффициентов регрессии
В этом случае используется применяемая при линейном регрессионном анализе матричная формула метода наименьших квадратов (МНК), которая с учётом кодирования факторов имеет вид:
|
|
|
|
|
T |
|
|
|
|
1 |
|
|
|
T |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
~ |
|
~ |
|
~ |
|
|
|
~ |
|
эксп |
|||||
|
|
|
|
|
y |
||||||||||
a |
|
|
|
|
|
|
|||||||||
m 1 1 |
|
m 1 n |
|
|
|
|
|
|
m 1 n |
|
|
||||
|
|
n m 1 |
|
n 1 |
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
14 Тема 02: Построение эмпирических статистических моделей ХТП
где кодированная матрица, зависящая от независимых переменных для двух факторов включает только +1 и -1 и имеет вид:
|
|
|
|
|
z10 |
z11 |
z12 |
|
1 |
1 |
1 |
||
|
|
|
|
|
z20 |
z21 |
z22 |
|
1 1 1 |
||||
|
|
||||||||||||
|
~ |
|
|
|
|||||||||
|
|
|
|
|
|||||||||
|
Z |
|
|
|
|
|
|
|
|||||
4 3 |
|
|
|
z30 |
z31 |
z32 |
|
|
1 |
1 |
|
||
|
|
|
|
|
|
|
|
1 |
|||||
|
|
|
|
|
|
z41 |
|
|
1 |
|
|||
|
|
|
|
|
z40 |
z42 |
1 |
1 |
Матрица z при активном экспериментировании называется матрицей планирования и обладает тремя оптимальными свойствами:
- симметричности: сумма элементов всех столбцов матрицы, кроме первого (точнее, нулевого) равна нулю
n
zij 0 ( j 1, ..., m);
i1
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
15 Тема 02: Построение эмпирических статистических моделей ХТП
- ортогональности: скалярное произведение двух любых столбцов матрицы равно нулю
n |
|
|
|
z Tj zu zij ziu |
0 |
( j, u 0, 1,..., m |
u j); |
i 1
- нормировки: скалярное произведение двух одинаковых столбцов матрицы
равно n ( n 2m в ПФЭ) |
|
|
n |
|
|
z Tj z j zi2j |
n |
( j 0, 1, ..., m). |
i 1 |
|
|
Благодаря перечисленным оптимальным свойствам матрицы планирования информационная матрица в ПФЭ при m=2 равна
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
16 Тема 02: Построение эмпирических статистических моделей ХТП
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
n |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
0 n |
0 , |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
||||||||
|
Z |
|
Z |
|||||||||||||
~ |
~ |
|
|
~ |
|
|
|
|
T |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
3 3 |
3 4 |
|
4 3 |
3 4 |
4 3 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
n |
т.е. она является диагональной с одинаковыми элементами на главной диагонали, равными n = 22 = 4.
Соответственно, корреляционная матрица C также будет диагональной и с одинаковыми элементами главной диагонали:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
1 |
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
T |
|
|
1 |
|
|
|
|
|
1 |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
~ |
|
~ |
~ |
|
|
|
T |
|
|
|
|
|
n 1 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
C |
|
|
|
Z Z |
|
0 |
0 |
|||||||||||||||
3 3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
0 |
n |
1 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
17 Тема 02: Построение эмпирических статистических моделей ХТП
Результатом подстановки последних соотношений в матричную формулу для определения кодированных коэффициентов регрессии будет простая формула:
~ |
n |
|
|
|
|
|
|
|
||
zij |
|
эксп |
n , |
|
j 0, 1, ..., m |
|||||
a j |
yi |
|
||||||||
|
i 1 |
|
|
|
|
|
|
|
||
При учёте взаимодействия двух факторов z1 и |
z2 кодированное уравнение |
|||||||||
регрессии принимает вид: |
|
|
|
|
|
|
|
|||
|
ˆ ~ |
|
~ |
|
|
~ |
|
~ |
||
|
y a0 z0 |
a1z1 |
|
a2 z2 |
a12z1z2 |
|||||
|
|
|
|
|||||||
и в матрицу планирования z |
включается ещё один дополнительный последний |
столбец, каждый элемент которого равен произведению элементов столбцов, соответствующих взаимодействующим факторам:
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
18 Тема 02: Построение эмпирических статистических моделей ХТП
~
4 4
|
|
|
z10 |
|
|
|
|
|
|
|
|
|
z20 |
|
|
|
|
||
Z |
||||
z30 |
||||
|
|
|
|
|
|
|
|
z40 |
z |
|
z |
|
z |
z |
|
|
|
1 |
1 |
1 |
1 |
|
11 |
12 |
|
11 |
12 |
|
|
|
|
|||||
z |
|
z |
|
z |
|
z |
|
|
|
1 |
1 |
1 |
|
|
|
|
|
|
|
1 |
|||||||
|
21 |
|
22 |
|
21 |
|
22 |
|
|
|
|
||
z31 |
z32 |
z31z32 |
|
|
1 |
1 |
1 |
|
|||||
|
|
1 |
|||||||||||
z41 |
z42 |
|
|
|
|
|
1 |
1 |
|
||||
z41z42 |
1 |
1 |
При этом матрица планирования сохраняет все три оптимальных свойства – симметричности, ортогональности и нормировки, а кодированный коэффициент уравнения регрессии при члене, характеризующем взаимодействие факторов, определяется по формуле:
~ |
n |
|
|
|
эксп |
n, |
j, u 1, ..., m |
u j |
|
a ju |
zij ziu yi |
|||
|
i 1 |
|
|
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
19 Тема 02: Построение эмпирических статистических моделей ХТП
В теории ПФЭ доказывается, что при увеличении числа факторов (m>2) матрица планирования z строится с использованием рассмотренной методики, в том числе и с учётом взаимодействия факторов (не только двойного, но и тройного, четверного и т.д.).
В этом случае число столбцов матрицы p зависит от числа учёта взаимодействий факторов (n 2m ) и матрица планирования сохраняет перечисленные оптимальные свойства.
Поэтому для определения кодированных коэффициентов регрессии используются приведённые выше формулы.
Для расчёта натуральных значений коэффициентов в кодированное уравнение |
|
регрессии вместо кодированных факторов z j j 1,...m |
следует подставить |
выражения для последних через натуральные значения факторов x j j 1,...m в соответствии с приведённой выше схемой кодирования.
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |