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

v0.5.7.final / 09

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

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

ОЦКП и обработка его результатов

Ортогональный центральный композиционный эксперимент (ОЦКП) относится к экспериментам II – го порядка, так как описывающее его уравнение yˆ включает факторы в квадрате и поэтому может описывать поверхности функций отклика в окрестности их экстремальных значений.

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

yˆ II a0 a1x1 a2 x2 a12x1x2 a11x12 a22x22

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

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

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

В соответствии с методикой ортогонального центрального композиционного плана эксперимента (ОЦКП) здесь, также как и для ПФЭ, осуществляется кодирование факторов по приведённой выше схеме (3 уравнения написанные выше!!), и для обеспечения ортогонального свойства матрицы планирования эксперимента в уравнение регрессии yˆ включается некоторая постоянная S.

В результате уравнение регрессии при m = 2 принимает вид:

ˆ ~

~

~

~

~

2

 

S

~

2

 

S

 

y a0 z0

a1z1

a2 z2

a12z1z2

a11

z1

 

a22

z2

 

 

Для определения большего числа кодированных коэффициентов, чем при обработке ПФЭ, и описания поверхности функции отклика вблизи её экстремума («почти стационарной области»), количество опытов в этом случае увеличивается.

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

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

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

При этом опыты, проводимые при ПФЭ (n 2m ) , дополняются опытами в

«звёздных» точках факторного пространства (n 2m) и опытами в центре

плана (nc ) с координатами z1 0 и

z2 0

«Звёздные» точки в факторном пространстве располагаются на осях координат на расстоянии +α и –α от центра плана эксперимента; причём величина α называется «звёздным» плечом и её значения, так же как величина S, определяются из условия ортогональности матрицы планирования z для ОЦКП.

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

N n n nc ,

или с учётом приведённых выше равенств:

N 2m 2m nc .

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

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

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

Для случая двух факторов (m = 2):

N 8 nc

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

( 1; 1)

( 1; 1)

 

 

 

 

 

(0;0)

 

 

 

z1

 

 

 

( 1; 1)

( 1; 1)

 

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

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

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

План проведения экспериментов в этом случае может быть представлен:

Матрица планирования z представляет собой часть плана проведения эксперимента без горизонтальных и вертикальных заголовков таблицы и вектора наблюдения y (правого столбца).

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

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

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

2m

2m

nc

n

p

z0

z1

z2

z1z2

z12 S z22 S

y

 

1

+1

-1

-1

+1

1-S

1-S

y1

2

+1

+1

-1

-1

1-S

1-S

y2

3

+1

-1

+1

-1

1-S

1-S

y

3

 

 

 

 

 

 

 

4

+1

+1

+1

+1

1-S

1-S

y4

5

+1

-α

0

0

α2-S

-S

y5

6

+1

+α

0

0

α2-S

-S

y6

7

+1

0

-α

0

-S

α2-S

y7

8

+1

0

+α

0

-S

α2-S

y8

9

+1

0

0

0

-S

-S

y9

 

 

 

 

 

 

 

 

N

+1

0

0

0

-S

-S

yN

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

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

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

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

В соответствии с методом наименьших квадратов эти коэффициенты определяются по матричной формуле:

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

1

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

~

 

~

 

 

 

 

~

 

эксп

 

 

 

 

 

 

 

 

 

 

 

y

 

a

 

 

 

 

 

 

 

 

 

m 1 1

 

 

m 1 n

 

 

 

 

 

 

 

m 1 n

 

 

 

 

 

 

n m 1

 

 

n 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

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

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

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

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

~

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C z

y

, где

a

 

 

 

 

 

T

 

1

 

 

 

 

 

~

 

 

 

 

 

 

 

 

C

z

z

 

 

 

 

 

 

 

 

Из-за свойства ортогональности матрицы планирования z необходимо определить только диагональные элементы информационной матрицы:

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

n

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 n

0 ,

~

 

~

 

 

~

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Z

 

Z

 

 

 

3 3

 

3 4

 

4 3

 

3 4

4 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

n

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

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

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

 

 

 

 

T

 

 

~

 

 

 

 

 

 

 

 

 

 

z,

z

 

а затем диагональные элементы корреляционной матрицы:

 

 

 

 

 

 

T

 

 

1

 

 

 

 

 

 

 

 

 

1

n 1

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

~

~

 

 

 

 

T

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

Z Z

 

 

0 n

 

0

3 3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

 

 

 

 

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

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

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

Определение экстремума функции отклика

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

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

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