v0.5.7.final / Тема 2
.pdf91 Тема 02: Построение эмпирических статистических моделей ХТП
и не учитывают взаимодействия более высоких порядков (тройные, четверные и т.д.), вероятность которых существенно меньше.
Последнее уравнение включает слагаемые с квадратами входных переменных
m
a j x2j
j1
иего коэффициенты получаются при обработке результатов активных
экспериментов II-го порядка (верхний индекс II при yˆ : yˆ ) – например, ОЦКП – ортогонального центрального композиционного плана эксперимента.
Предпоследнее уравнение не включает слагаемые с квадратами входных переменных и его коэффициенты получаются при обработке результатов активных экспериментов I-го порядка – верхний индекс I при yˆ : yˆ - например, ПФЭ – полный факторный эксперимент.
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
92 Тема 02: Построение эмпирических статистических моделей ХТП
При определении оптимальных условий проведения процесса с использованием эмпирических моделей (например, методом Бокса-Вильсона) выходная переменная yˆ является критерием оптимальности или целевой функцией.
В теории активного экспериментирования выходную (зависимую) переменную принято называть функцией отклика, а входные (независимые) переменные – факторами. Соответственно - координатное пространство с координатами - факторным пространством, а геометрическое изображение функции отклика в факторном пространстве – поверхностью отклика.
Активный эксперимент планируется таким образом, чтобы упростить обработку его результатов методами регрессионного и корреляционного анализа.
Ортогональные планы экспериментов, используемые при активном экспериментировании, обеспечивают диагональный вид корреляционной матрицы C при регрессионном анализе и, соответственно, статистическую независимость коэффициентов регрессии.
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
93 Тема 02: Построение эмпирических статистических моделей ХТП
К другим достоинствам активного экспериментирования относятся:
-возможность предсказания количества опытов, которые следуют провести;
-определение точек факторного пространства, где следует проводить опыты;
-отсутствие проблем, связанных с выбором вида уравнения регрессии;
-возможность определения оптимальных параметров процесса экспериментально-статистическим методом;
-сокращение объёма опытных исследований.
ПФЭ и обработка его результатов
Полный факторный эксперимент (ПФЭ) относится к экспериментам I-го порядка, т.к. описывающее его уравнение yˆ не включает факторы в квадрате.
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
94 Тема 02: Построение эмпирических статистических моделей ХТП
Для двух факторов и без учёта взаимодействия факторов соответствующая эмпирическая модель может быть записана:
yˆ a0 a1x1 a2 x2
В соответствии с теорией ПФЭ при проведении опытных исследований каждый из факторов варьируется только на двух уровнях – минимальном (кодированное значение -1) и максимальном (кодированное значение +1).
При этом реализуются возможные комбинации минимальных и максимальных значений факторов, в результате чего общее число опытов (n) в ПФЭ равно 2m и
полный факторный эксперимент обычно называется ПФЭ типа |
2.m |
Для определения числа опытов применяется формула:
n 2m
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
95 Тема 02: Построение эмпирических статистических моделей ХТП
В последнее уравнение включаются кодированные значения факторов z j вместо x j , значения которых получаются по следующей схеме кодирования:
|
|
x |
j |
x 0 |
( j 1, ..., m) |
||
z j |
|
|
j |
||||
|
|
x j |
|
||||
|
|
|
|
|
|
|
|
x 0 |
0.5 xmin |
xmax |
( j 1, ..., m) |
||||
j |
|
|
|
j |
|
j |
|
x |
j |
0.5(xmax xmin ) |
( j 1, ..., m) |
|
|
j |
j |
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
96 Тема 02: Построение эмпирических статистических моделей ХТП
В результате план проведения эксперимента с учётом вышесказанного и кодирования факторов имеет вид (число факторов равно 2 - m = 2, число опытов n = 2m = 22 = 4):
n |
p |
z0 |
z1 |
z2 |
y |
|
|
|
|
|
|
1 |
|
1 |
1 |
1 |
y1 |
2 |
|
1 |
1 |
1 |
y2 |
3 |
|
1 |
1 |
1 |
y3 |
4 |
|
1 |
1 |
1 |
y4 |
|
|
|
|
|
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
97 Тема 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) |
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
98 Тема 02: Построение эмпирических статистических моделей ХТП
Для параметрической идентификации кодированного уравнения регрессии используется метод регрессионного анализа, включающий три этапа:
- определение кодированных коэффициентов уравнения регрессии ~ a
методом наименьших квадратов;
- оценка значимости кодированных коэффициентов регрессии с использованием t – критерия Стьюдента;
- проверка адекватности кодированного уравнения регрессии с использованием F – критерия Фишера.
Реализация двух последних этапов возможна при выполнении свойства однородности дисперсий (одно из требований регрессионного анализа) и проведении параллельных опытов, например, в центре плана
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
99 Тема 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 |
|
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |
100 Тема 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
РХТУ им. Д.И. Менделеева |
Кафедра информатики и компьютерного моделирования |