Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

v0.5.7.final / 03

.pdf
Скачиваний:
7
Добавлен:
13.02.2015
Размер:
404.69 Кб
Скачать

11 Тема 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 в соответствии с приведённой выше схемой кодирования.

РХТУ им. Д.И. Менделеева

Кафедра информатики и компьютерного моделирования

Соседние файлы в папке v0.5.7.final