Скачиваний:
12
Добавлен:
13.06.2022
Размер:
1.7 Mб
Скачать

91 Тема 02: Построение эмпирических статистических моделей ХТП

и не учитывают взаимодействия более высоких порядков (тройные, четверные и т.д.), вероятность которых существенно меньше.

Последнее уравнение включает слагаемые с квадратами входных переменных

m

a x2 j j

j 1

и его коэффициенты получаются при обработке результатов активных

 

ˆ

 

ˆ

 

 

 

 

 

 

экспериментов II-го порядка (верхний индекс II при

y

:

 

) – например,

 

y

 

ОЦКП – ортогонального центрального композиционного плана эксперимента.

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

активных экспериментов I-го порядка – верхний индекс I при :

yˆ

 

- например,

 

ПФЭ – полный факторный эксперимент.

 

 

 

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

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

92 Тема 02: Построение эмпирических статистических моделей ХТП

При определении оптимальных условий проведения процесса с использованием эмпирических моделей (например, методом Бокса-Вильсона) выходная

переменная

yˆ

является критерием оптимальности или целевой функцией.

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

Активный эксперимент планируется таким образом, чтобы упростить обработку его результатов методами регрессионного и корреляционного анализа.

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

матрицы

C

при регрессионном анализе и, соответственно, статистическую

независимость коэффициентов регрессии.

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

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

93 Тема 02: Построение эмпирических статистических моделей ХТП

К другим достоинствам активного экспериментирования относятся:

-возможность предсказания количества опытов, которые следуют провести;

-определение точек факторного пространства, где следует проводить опыты;

-отсутствие проблем, связанных с выбором вида уравнения регрессии;

-возможность определения оптимальных параметров процесса экспериментально-статистическим методом;

-сокращение объёма опытных исследований.

ПФЭ и обработка его результатов

Полный факторный эксперимент (ПФЭ) относится к экспериментам I-го порядка, т.к. описывающее его уравнение не включает факторы в квадрате.

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

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

94 Тема 02: Построение эмпирических статистических моделей ХТП

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

ˆ

 

a

 

a x

a

 

x

 

y

 

0

2

2

 

 

 

1

1

 

 

В соответствии с теорией ПФЭ при проведении опытных исследований каждый из факторов варьируется только на двух уровнях – минимальном (кодированное значение -1) и максимальном (кодированное значение +1).

При этом реализуются возможные комбинации минимальных и максимальных значений факторов, в результате чего общее число опытов (n) в ПФЭ равно 2m и

полный факторный эксперимент обычно называется ПФЭ типа

m

2.

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

n 2m

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

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

95 Тема 02: Построение эмпирических статистических моделей ХТП

В последнее уравнение включаются кодированные значения факторов

z

j

 

вместо

x

j , значения которых получаются по следующей схеме кодирования:

 

 

 

 

x

 

x

0

 

z

 

 

j

j

( j 1, ...,m)

 

 

 

 

j

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j

 

 

 

 

 

 

 

 

 

 

xj0 0.5 xminj

xmaxj

 

( j 1, ..., m)

x

 

0.5(x

max

x

min

)

( j 1, ...,m)

j

j

j

 

 

 

 

 

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

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

96 Тема 02: Построение эмпирических статистических моделей ХТП

В результате план проведения эксперимента с учётом вышесказанного и кодирования факторов имеет вид (число факторов равно 2 - m = 2, число опытов n = 2m = 22 = 4):

n

1

2

3

4

p

z

0

 

1

1

1

1

z

1

 

1

1

1

1

z2

1

1

1

1

y

 

 

y

 

 

 

 

 

1

 

y

 

 

2

 

 

 

y

 

 

 

3

 

y

 

 

4

 

 

 

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

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

97 Тема 02: Построение эмпирических статистических моделей ХТП

При этом уравнение регрессии, описывающее эти опытные данные,

 

 

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

z

j

( j 0,1,2)

 

 

~

~

~

 

соответственно, кодированных коэффициентов регрессии

 

 

:

 

 

a

0

, a

, a

2

 

 

 

 

1

 

 

и,

ˆ

~

~

 

~

z

 

y a

a z

a

2

 

0

1

1

2

 

В кодированном факторном пространстве в соответствии с указанным планом проведения эксперимента проведённые опыты представляются точками вершин

квадрата:

z

 

 

2

 

 

( 1; 1)

 

( 1; 1)

z

1

 

( 1; 1)

( 1; 1)

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

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

98 Тема 02: Построение эмпирических статистических моделей ХТП

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

- определение кодированных коэффициентов уравнения регрессии методом наименьших квадратов;

~ a

-оценка значимости кодированных коэффициентов регрессии с использованием t – критерия Стьюдента;

-проверка адекватности кодированного уравнения регрессии с использованием F – критерия Фишера.

Реализация двух последних этапов возможна при выполнении свойства однородности дисперсий (одно из требований регрессионного анализа) и проведении параллельных опытов, например, в центре плана

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

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

99 Тема 02: Построение эмпирических статистических моделей ХТП

При проведении k параллельных опытов в центре плана

( y

эксп

, s 1, ...,

0s

 

 

 

y

эксп

 

 

 

значение

c

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

 

измерений во всех параллельных опытах:

k

)

среднее

k

 

ycэксп y0экспS

k

S 1

 

Определение кодированных коэффициентов регрессии

В этом случае используется применяемая при линейном регрессионном анализе матричная формула метода наименьших квадратов (МНК), которая с учётом кодирования факторов имеет вид:

~

 

 

 

T

 

 

 

1

 

 

T

 

 

 

 

 

 

 

 

 

 

 

~

~

~

y

эксп

a

 

 

 

 

 

 

m 1 1

 

m 1 n

n m 1

m 1 n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

100 Тема 02: Построение эмпирических статистических моделей ХТП

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

 

 

 

 

z

 

z

 

z

 

 

 

1

1

1

 

 

 

10

11

12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

z

 

z

 

1

1

 

 

Z

 

 

 

 

 

1

 

 

 

20

 

21

 

22

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

4 3

 

z

 

z

 

z

 

 

 

1

1

 

 

 

 

 

30

31

32

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

z

 

z

 

1

1

 

 

 

 

40

41

42

 

 

1

 

 

 

 

 

 

 

 

 

Матрица

z

при активном экспериментировании называется матрицей

 

планирования и обладает тремя оптимальными свойствами:

- симметричности: сумма элементов всех столбцов матрицы, кроме первого (точнее, нулевого) равна нулю

n

zij 0 ( j 1, ..., m);

i 1

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

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

Соседние файлы в папке Лекции ХТП