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

Оптимизация эксперимента в химии и химической технологии

..pdf
Скачиваний:
38
Добавлен:
15.11.2022
Размер:
14.6 Mб
Скачать

тив из каждой вершины треугольника высоту, разделив к аж дую из них на десять равных по величине отрезков и проведя через полу­ ченные деления прямые, параллельные сторонам треугольника, по­ лучим треугольную сетку. Приближение от каждой данной стороны к противоположной вершине отвечает пропорциональному возра­ станию содержания соответствующего компонента, поэтому пря­ мые, параллельные данной стороне, при последовательном перехо­ де от одной прямой к другой отражают возрастание третьего компонента на 10%. На рис. 44 на .соответствующих высотах тре-

Рис. 44. Концентрационный тре-

Рис.

45.

Концентрационный тре­

угольник Гиббса

 

угольник Розебума

угольника указано содержание

каждого

из

компонентов (в точке

d — 30% А и 70% С). Практически большей частью не прибегают к построению высот, можно откладывать содержание компонентов непосредственно на сторонах треугольника. Такой способ отсчета принят в треугольнике Гиббса. В треугольнике Розебума состав тройной системы отсчитывается по трем отрезкам одной стороны треугольника (рис. 45). В концентрационном треугольнике точки, лежащие на прямой, выходящей из вершины треугольника, соот­ ветствуют смесям с постоянным отношением содержаний компонен­ тов, изображаемых двумя другими вершинами. Свойство (у) обыч­ но представляют проекциями линий равного значения на плоскость концентрационного треугольника.

При <7 = 4 правильный симплекс — тетраэдр, каждая вершина которого соответствует чистым компонентам. Ребро представляет собой двухкомпонентную систему, грань — трехкомпонентную. Точ­ ки внутри тетраэдра соответствуют четырехкомпонентным систе­ мам. Так, компонент Х\ отсутствует на грани х2>х3, Х4 , а по сечени­ ям тетраэдра, приближающимся к вершине Xi, содержание компо­ нента Х\ увеличивается.

Представление кривых свойств этой системы на плоскости не представляется возможным. Поэтому графически такую систему представляют в виде сечений трехмерного симплекса плоскостями,

перпендикулярными одной из его осей. Состав четырехкомпонент­ ных смесей, лежащих в плоскости сечения, определяется уже дву­ мерным симплексом, что позволяет изменение свойств системы представлять в виде контурных кривых. При этом в одном сечении варьируют только тремя компонентами. Переход от одного сечения к другому соответствует изменению четвертого компонента.

При планировании эксперимента для решения задач на диа­ граммах состав — свойство нредполагается, что изучаемое свойство является непрерывной функцией аргументов и может быть с доста­ точной точностью представлено полиномом. Использование методов планирования эксперимента позволяет значительно сократить объем эксперимента при изучении многокомпонентных систем, от­ падает необходимость в пространственном представлении сложных поверхностей, так как свойства можно определять из уравнений. При этом сохраняется возможность графической интерпретации результатов.

Поверхности отклика в многокомпонентных системах имеют, как правило, очень сложный характер. Для адекватного описания таких поверхностей необходимы полиномы высоких степеней и, сле­ довательно, большое количество опытов. Обычный полином степе­

ни п от q переменных имеет С

коэффициентов:

 

 

 

 

у — Ьо + 2 ^1*1

 

2

&ijx ixj +

 

 

 

 

1 <i<q

1 <K]<q

 

 

 

+

2

l>l j kx iX j x k + . . . + ^

b !iU........» * / , * / , * / •

(VI. 2)

 

1 < i < ) < b < q

 

 

n

n

 

Соотношение

Я

 

исключить

q-й

компонент и

2 xi = 1 позволяет

 

 

i=1

 

 

 

 

 

снизить число коэффициентов до C"+n-i. Однако по существу зада­ чи желательно ввести в модель все q компонентов.

Шеффе [36] предложил описывать свойства

смесей приведен­

ными полиномами, получаемыми из

(VI.2) с учетом условия нор­

мированное™ суммы независимых

переменных

(VI. 1). Покажем,

например, как получить такой приведенный полином второй степе­ ни для тройной системы. Общий вид полинома

1/ = *0 + *1*1 + *2 *2 +

*3 *3 +

*12 *1*2 + *13 *1*3 +

 

+ *23*2*3 +

*11*1 +

*22*2 + *33*3-

(VI. 3)

Так к!ак

 

 

 

 

 

*1 + *2 + *з = 1,

(V I.4)

то

 

 

 

 

(V I.5)

*0*1 +

*0*2 + *0*3 =

*0-

Умножив (VI.4) последовательно на Х\,

* 2 и *з, получим

 

*1 = * 1 — *1*2 — * 1 * 3 .

 

* 2 =

*2 — *1*2 — * 2 * 3 .

(VI. 6)

*3 =

*3 — *1*3 — *2*3 -

 

Подставим (VI.5) и (VI.6) в (VI.3) и приведем подобные члены:

У (*о + + Ьп ) х \ 4- ( Ь 0 + 62 + ^22) * 2 + (6о + + *зз) *з-+

+ (^ 12^11 622) * 1*2 + (^13 — ^11 — *33) *1*3 + (^23 — ^22 — ^зз) * 2 * 3 - (V I . 7)

Обозначим

Р/= + */+

Ри = Ь ц Ь ц — b j j .

(VI.8)

Тогда получим приведенный полином второй степени от трех пере­ менных:

У — Pl*l + Рг* 2 + Рз*з + Р12*1*2 + Р13*1* 3 + Ргз*2*з-

( V ( .9 )

Таким образом, число коэффициентов уменьшилось с 10 до 6. Приведенный полином второй степени от q переменных

У =

2

Р /* / +

2

Р/У*/*У

(V I. 10)

 

1 « < ?

 

\ < i < ) < q

 

 

содержит q -\-С1 = С2д+ 1 коэффициентов. Приведенный полином не­ полного третьего порядка для трехкомпонентной смеси

У = Pl*l + Рг* 2 ■+• Рз*3 + Р12*1*2 + Р13*1*2 + Р23*2*3 + Pl23*l*2*3i (VI. 11)

для ^-компонентной смеси

0

=

2

р/ * / +

2

Ри * 1 * } +

2

РU k X i X j X k .

( V I . 12)

 

 

1 < 1 < ?

l < i < } < q

 

1 < i < ) < k < . q

 

 

Приведенный полином

третьего

порядка

для трехкомпонентной

смеси

 

 

 

 

 

 

 

 

 

 

У

= Pl*l + Рг* 2 + Рз* 3 + Pl2*l*2 + Pl3*l* 3 + Р23*2*3 +

 

+ Yl2* l *2 (*1 — * 2) +

Y l3*l*3 (*1 — *3> +

Y23*2*3 (*2 — *3> + P l2 3 * l* 2 * 3 .'

(VI. 13)

для ^-компонентной

 

 

 

 

 

 

 

У =

2

р1*1 +

2

ViJXiXj] +

2

ViJXtXj (хI Xj) +

 

l < K q

1 < l < j < q

 

K i < ] < q

 

 

 

 

 

 

 

+

2

PU k X i X j X k .

 

( V I . 14)

 

 

 

 

l < K ] < k < q

 

 

 

 

 

Приведенный' полином

четвертого

порядка

для

трехкомпонентной

смеси

 

 

 

 

 

 

 

 

 

 

У = 01*1 +

02*2 +

03*3 +

012*1*2 + 013*1*3 +

023*2*3 +

Y l2*l*2 (*1 — * 2 ) +

+Yl3*l* 3 ( * 1 — *3)+ Y23*2 *3 ( * 2 —*3) + ^12*1 *2 ( * 1 ~ *2>2 +

+^13*1*3 (*1 — *3)2 + ^23*2*3(*2 — *3)2 + 01123*1*2*3 +

01223*1*2*3 "Ь 01233*1*2*3»

(V I .1 5 )

для ^-компонентной

 

У =

2 Р ‘х ‘ +

2

$ U x ix J +

2 Y i i x ix J (x i — X J) +

 

 

 

1 <1<я

l < i < ] < q

К

/ < ] < q

 

 

 

 

+

2

l u x ix i ( * / X J ) 2 +

2

$ u ) k x ^x j x k +

2

$ ij } k x ix ) x k +

 

1 </ <]< q

+

2

 

\ < l < ] < k < q

2

 

K i < j < k < q

 

 

 

 

$ i l m x i x l x l +

 

h i k i X i X j X k X i -

(VI. 16)

 

 

 

l < l < ] < k < q

 

K i < J < k < l < q

 

 

 

 

Нелинейная часть этих полиномов называется синергизмом, если она вызывает увеличение отклика по сравнению с откликом, предсказываемым линейной частью уравнения, и антагонизмом — при уменьшении отклика. Например, член fiijXiXj в полиноме вто­ рого порядка называют квадратичным коэффициентом бинарного синергизма компонентов i и /. В полиноме третьего порядка синер­ гизм трехкомпонентной смеси равен

[VtjXiXj + YijXiXj (ЛГ/ - xj)] 4- [tibXiXk + yikXixk(xt xk)] 4-

 

+ [VjbXjXk 4- VkXjXk (Xj — Xk-)] 4- hluXiXjXb,

(VI.17)

где в квадратных скобках — бинарные синергизмы в тройной систе­ ме; Pij/t — кубический коэффициент тройного синергизма компонен­ тов /, /, k.

Возможно другое преобразование исходного полинома степени п от q переменных. Его можно свести к так называемому однород­

ному

полиному

степени

п,

умножая

члены

степени

s<n на

( 2

x l\n~s = l.

Приведем полином

(VI.3) к

однородному:

Ч</«7

/

 

 

 

 

 

 

 

У= fy) (*1 + * 2

4- *з) 2 + bXXi (Xi + х2 + х3) + Ь2Х2 (XI + х2+ х3) +

+ Ь3* 3 ( * 1 + х2 4- *з) + bi2xix2 + ^13^1 ^ 3

+ Ь23Х2 Х3 + Ьпх\ + Ъ22х\ 4-

4- Ъ33х\ = (2*0 “Ь Ь\ 4- Ь2 4- ^12) Х\Х2 4- (260 4- ^1 4- 63 4- 613) Х\Х3 4-

 

4- (2^о 4- Ь2 4 - Ь3 4- Ь23) х2х3 + (6 0 + Ь\ 4- 6 ц) х\ 4-

 

 

 

4- (6 0 + Ь2

+ 622) х\ + (6 о 4 - 6 3 + 633)*3 .

(VI. 18)

Обозначим в (VI.18)

 

 

 

 

 

 

 

?ij= 2^о + */ + Ь] 4- bijt

=

6 0 4 - 6 / 4 - 6 //.

 

Получим однородный полином:

 

 

 

 

 

 

« = 012*1*2+ 013*1*3 +

023*2*3 +

011*1 + 022*2 + 033*3-

(V I. 19)

Число неизвестных коэффициентов в (VI. 19) то же, что и в (VI.9). Приведенные полиномы получили более широкое применение и в дальнейшем будем использовать только их. Таким образом, ми­

нимальное число экспериментальных точек для определения коэф­

фициентов полинома степени п от q переменных составляет Cq+n _i (табл. 63).

Число опытов для полиномов разных степеней

 

 

Степень полинома

 

 

 

Степень полинома

 

Число ком­

 

(непол-

 

 

Число ком­

 

(не^ол-

 

 

понентов

2

3

4

понентов

2

3

4

 

3 ная)

 

•3 ная)

3

6

7

10

15

6

21

41

56

126

4

10

14

20

35

8

36

92

120

330

5

15

25

35

70

10

55

175

220

715

2. Симплекс-решетчатые планы Шеффе. В настоящее время наи­ большее применение получили симплекс-решетчатые планы, пред­ ложенные Шеффе [37, 38]. Эти планы обеспечивают равномерный

разброс экспериментальных точек

по (q—1)-мерному

симплексу.

 

 

 

 

 

 

Экспериментальные

точ­

 

 

 

 

 

 

ки

представляют

{qt п}-

 

 

 

 

 

 

решетку

 

на

 

симплексе,

 

 

 

 

 

 

где q — число

компонен­

 

 

 

 

 

 

тов

смеси;

 

п — степень

 

 

 

 

 

 

полинома.

 

Симплекс-ре­

 

 

 

 

 

 

шетчатые

планы

являют­

 

 

 

 

 

 

ся насыщенными

плана­

 

 

 

 

 

 

ми. По каждому

компо­

 

 

 

 

 

 

ненту'

имеется

(п + 1)

 

 

 

 

 

 

одинаково

 

расположен­

 

 

 

 

 

 

ных

уровней

 

Х\ ~ 0,

1//г,

 

 

 

 

 

 

2//г, ..., 1 и

берутся

все

 

 

 

 

 

 

возможные комбинации с

 

 

 

 

 

 

такими

значениями

кон­

 

 

 

 

 

 

центраций

 

компонентов.

 

 

 

 

 

 

Так, например, для

квад­

 

Рис. 46. {3, п]-решетки:

 

 

ратичной решетки {q, 2},

а — для

полинома

второго

порядка;

б — для

по­

обеспечивающей

прибли­

линома

неполного

третьего

порядка;

в — для

по­

жение

поверхности

от­

линома

третьего порядка;

г — для полинома

чет­

 

вертого порядка

 

 

клика полиномами

вто­

 

 

 

 

 

 

рой

степени

(п = 2),

дол­

жны быть использованы следующие уровни каждого

из фак­

торов: 0, 7г и 1, для

кубической (п = 3) — 0,

!/з, 2/з

и

1 и т. д. Не­

которые {3, п)-решетки представлены на рис. 46, а

{4, п) — на рис.

47. Эти планы частично композиционные. Неполную

 

кубическую

решетку {3, 3*}, например, можно получить

из

{3,

2},

добавив

только одну точку в центре симплекса; решетку

{3, 4}— добавле­

нием точек к

решетке {3,

2}.

 

 

 

 

 

 

 

 

 

Записав координаты точек симплексной решетки, получим мат­ рицу планирования. Построим матрицы планирования для решеток {3,2}, {3,3} и {3,4}.“

Индексы у свойства смеси указывают на относительное содер­ жание каждого компонента в смеси. Например, смесь 1 (табл. 64)

состоит только из компонента Х\, свойство этой смеси обозначается У\, смесь 4 состоит из 72*1 и Ч2 Х2 , свойство обозначается у 12.

Матрица планирования для симплексной' решетки {3, 3} приве­ дена в табл. 65.

*/

ах,

в

Рис. 47. {4, п}-решетки:

а — для

полинома

второго

порядка;

б — для полинома

неполного

третьего

порядка;

в —-для

полинома

третьего

порядка; г — для

полинома

четвер­

 

 

 

того порядка

 

 

В табл. 66 приведена матрица планирования для построения полинома четвертой степени в трехкомпонентной системе.

Коэффициенты приведенных полиномов получают, используя свойство насыщенности плана. Для получения коэффициентов по­ линомаУ

У — P l* l + ?2Х 2

Рз*3 + $12х Iх 2 + P l3 * l* 3

? 2 3 *2 *3

 

 

 

М а т р и ц а п л а н и р о в а н и я

д л я

(3, 2 } -р еш етки

 

 

 

N

 

**

 

хя

^ЭКС

 

N

 

*4

X 2

Ха

^ЭКс

 

 

0

 

0

 

 

4

 

 

1

 

 

1

1

 

У\

 

1 / 2

1 /2

0

 

 

2

0

1

 

0

 

5

 

1 /2

0

1 /2

У?

 

у \

 

 

3

0

0

 

1

 

 

6

 

0

1 /2

1 /2

у

 

*>3

 

 

У23

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

65

 

 

С и м п лекс-реш етчаты й

п л ан

т р етье го

п о р я д к а

 

 

 

 

 

 

д л я тр ех ко м п о н ен тн о й

смеси

 

 

 

 

N

 

х 2

х%

 

у

1

N

 

Хх

Ха

Х$

У

 

 

 

 

 

 

 

|

 

 

 

 

 

 

1

1

0

0

 

У\1

 

6

 

0

2/3

1/3

У223

 

2

0

1

0

 

Уъ

 

7

 

0

1/3

2/3

У233

 

3

0

0

1

 

Уо

 

8

 

2/3

0

1/3

Уц з

 

4

2/3

1/3

0

 

* 4 ,2

 

9

 

1/3

0

2/3

0*133

 

5

1/3

2/3

0

 

У122

 

10

 

1/3

1/3

1/3

У123

 

 

 

 

 

 

 

 

 

 

 

 

Т а б л и ц а

66

 

 

М а т р и ц а

п л а н и р о в а н и я д л я

{3, 4} -р еш етки

 

 

 

N

 

 

 

У

 

N

 

 

**

х »

У

 

i

1

0

0

 

01

 

9

 

3/4

0

I 1/4

У1113

 

2

0

1

0

 

У2

 

10

 

1/4

0

3/4

ioqq

 

 

0

0

1

 

 

11

 

0

3/4

1/4

хооо

 

з

 

Уз

 

 

У0093.

 

4

1/2

1/2

0

 

012

 

12

 

0

1/4

3/4

У ОоОО

 

5

1/2

0

1/2

 

 

13

 

1/2

1/4

1/4

2333

 

 

013

 

 

л 1 1 ПО

 

б

0

1/2

1/2

 

У23

 

14

 

1/4

1/2

1/4

“ 1223

 

7

3/4

1/4

0

 

УИ 12

 

15

 

1/4

1/4

1/2

У 10QQ

 

8

1/4

3/4

0

 

01222

 

 

 

 

 

 

1/0«3

 

 

 

 

 

 

 

 

 

 

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

шести точек матрицы планирования

(табл.

64). Тогда при подста­

новке координат первой точки

(xi = l ,

* 2 = 0 ,

*з = 0 ) получим

0 i = Pi-

 

(VI. 20)

Соответственно

 

 

 

р2 = 02

И Рз =

Д/3.

(VI.21)

При подстановке в уравнение координат четвертой точки получим

У12 = P1V2 + P2V2 + Pl2'/4*

(VI 22)

Но так как р*=уи то

У\2 = У1Ч2 + m 'h + Эj21/4.

Pl2 = 4*/12 —2(/1 2 y 2 .

(VI. 24)

Соответственно

Pi3 = 4(/13 — 2y i — 2yz , Ргз = 4#2з — 2^2 — 2#3.

(VI.25)

Т р и точки, определяющие коэффициент P i j , лежат на одном ребре. Коэффициенты приведенного полинома второго порядка для ^-ком­ понентной смеси

у =

2 *$‘х ‘ + ^ 2

Puxiyt

 

 

\ < i < q

l < l < j < q

 

 

определяются аналогично:

 

 

 

Рi =

Vi>

=

( V I .26)

При последовательной подстановке в полином третьего порядка для трехкомпонентной смеси

V = M l + Р2*2 + Рз*3 +

Pl2*l-*2 +

Pl3*l*3 + р23*2*3 + Yl2*l*2 (*1 “ * 2) +

+ Yl3-*1*3 0*1 — *3) +

V23*2*3 С*2— *3) + Pl23*l*2*3

 

координат точек матрицы планирования

(см. табл. 65)

имеем:

 

 

P l= l//b

р2= //2» Рз=

у,я»

( V I .27)

У т

=

У \ 21з +

У 2 Ч 3 4* Pi22/9 +

Yi22/27»

( V I .28)

У ш

=

У х Ч з +

У221з + Р122/9

Yl22/27

( V I .29)

Просуммировав уравнения

(VI.28) и (VI.29), получим

 

Pl2 = 9А

(1/112 + //122 “ Ух — У2)

(VI. 30)

Вычитая из (VI.28) уравнение (VI.29), получим Y12:

 

Y12= 9/4 (Зг/112Зе/!22У \ + у %).

(VI. 31)

Подставляя в приведенный полином координаты точек 6 и 7, по­ лучим:

P23 =

91л ( у 223 + У 233

У 2

Уз) >

 

 

 

(VI. 32)

Y23 =

9U (3//223 — &У233 У 2 + У з)

И аналогично после подстановки координат точек 8 и 9

Pis =

9А (1/пз + угзз -

У \ —

У з ) ,

Y13 =

 

 

(V I. 33)

9А (3 i/ii3 — Зугзз у \

+ # 3).

Подставляя координаты последней точки с учетом соотношений

(VI.27) —(VI.33),

получим:

 

 

 

 

 

 

У\23 =

У \ Ч з

+ У 2 Ч 3 +

У з Ч з +

74(1/112 + 1/ 122

У \ “

У 2) +

+ 74 (У П З + У133 — i/1 Уз) +

7 4 (//223 +

У233У2 У з ) + 727р123»

Pl23 =

2 7 у 123 — 27/4 112 +

//122 + УХ13

+ //133 +

//223 +

У 23з) +

 

 

+

9h (//1 +

У2 +

У з )

 

( V I . 34)

Коэффициенты полинома третьего порядка для ^-компонентной смеси

й = 2

$ ‘ х ‘ + 2

h j X i X j + 2 Yt j x i x j ( х , X j ) +

 

1 < / « 7

 

1 < / < У « 7

1 < х < ; « 7

 

 

 

 

+

2

V l J k X l X j X k ,

 

 

 

 

 

l < i < J < q

 

 

 

 

 

 

Ь =

У1,

 

(VI. 35)

 

P/У= 9А (у//у + У Ш — У1 — y j ) ,

(VI.36)

 

У и

= 9/4 (3////У -

3///уу — 0/ +

//у),

(VI. 37)

Р/У/t = 27y u k — 27у4

/zy +

y i j j

+ у i i k 4- y i k k 4- y j j k 4- y j k k ) + 9/ 2 (#/ 4- y j

4- i/ft)-

 

 

 

 

 

 

(V I. 38)

Аналогично

выводятся

соотношения для

определения коэффи­

циентов любого приведенного полинома при любом числе компо­ нентов. Так, для полинома неполного третьего порядка для трех­ компонентной смеси имеем:

= Pi*! 4- Рг*2 4- рз*з 4- Pi2*i-*2 4" Pi3*i*3 4- Ргз-^г-^з 4- Pi23*i*2*3 >

Р/ =

У\ и т. д.

 

( V I . 39)

Pi2 = 4 //i22//1 —

2 и

Т .Д .,

( V I . 40)

Р123 = 27//123 — 12 (//12 4 -

//1з 4-

//2з) 4 -

3 (//14- //2 4- # з)*

( V I . 41)

Для /7 -компонентной смеси

У =

2

РIх i 4"

2

Р i j x l^ J 4"

2

Рi j k * Iх ]X k t

 

l < K q

 

1 < l < j < q

1 < l < j < k < q

 

 

 

 

 

 

P1 = У1.

 

 

( V I . 42)

 

 

 

P/У =

41//У — 2У1 -

2У]>

 

( V I . 43)

Pi j k

= 27//,y/5— 12 (///y 4- ///* 4- y j k ) 4- 3 (/// 4-

y j 4- //ft).

( V I . 44)

Полином четвертого порядка для трехкомпонентной смеси:

У =

P i* i 4 - р2*2 4- Рз*з 4 - Pi2*i* 2 4- Pi3*i*3 4 -

р23*2*3 4 - Yi2*i*2 ( * I * 2) 4 -

 

4- Y i3*i*a C*i — * з )

4- Y23*2*3 ( * 2 — * з )

4 - *12* 1*2 (*1 — * г ) 2 +

4-

*13*1*3 (*1 — *з)2 +

*23*2*3 (*2 *з)2 4- Pll23*?*2*3 4- Pi22:j*i*2*3 +

 

 

4-

Pi233*i*2*f.

( V I . 45)

 

 

Pi =

1/1 И т. д .,

( V I . 46)

Pl2 = 12 — 2х/1 — 2//2,

Vi2 =

8/з ( —*У \

4- 2//Ш 2 —

2 //1222 4- y i ) ,

Y13 =

8/з (

У \

4- 2//Ш З — 2//1ззз +

//3) ,

Y23 =

8/з (

У 2 4* 2//2223—

2//2333 4 -

//з) »

* 1 2 = 8/з ( — У \

4 - 4 - / / 1 1 1 2 6 / / 1 2

4 - 4 / / 1 2 2 2

~ У 2) >

 

 

 

 

®1з =

8/з ( — Ш

+

4 « /ш з —

6i/i3 +

4</1ззз — у з ) ,

 

 

 

 

 

 

®23=

8/з ( ~

У 2

+

41/2223 — 6 1/23 +

41/2333 — //з ) .

 

 

PU 23 = 32 (3//П 23 1/1223 </123з) + 8/з Ф У 1

~

У 2 ~

У з ) ~

I6 (//12 + //13> “

 

 

 

 

16/з

( 5 1 / 1 1 1 2

+

5 1 / 1 1 1 3

31/1222 —

З 1 / 1 3 3 3

1/2223 —

У233з),

 

 

01223 = 32 (31/1223 1/1123 У т з )

+

8/з

( б //2 1/1 //з) — 16 ('1/12 + </1з)

 

 

 

 

16/з (51/1222 +

51/2223 _

3l/1112 —

3//2333 ~ ’ //1113 _

Й133з)*

 

 

01233 = 32 (31/1233 1/1123 У т з )

+ 8/з (6 //3 У 1 У 2) 16 (l/13 + У 2з)

 

 

 

 

16/з (51/1333 +

51/2333 — 3 1 /ш з —

31/2223—

//1112 — Z/l222>-

 

 

Для ^-компонентной смеси:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

£ =

2

М

/ +

1 <

2

 

 

P//-*w +

2

 

Yu

x

i x j i

x i — x j ) +

 

 

 

 

1</<?

 

 

 

K ) < q

 

 

 

l < K

) < q

 

 

 

 

 

 

 

 

 

 

+

 

2

 

Biyjffjey (Jf/ — JC/)2

2

 

$nj kx )x j x k +

 

 

 

 

 

 

l < i < ) < g

 

 

 

 

 

 

 

 

l < l < ) < k < q

 

 

 

 

 

 

+

2

 

$ U l k x i x ) x k +

 

 

2

 

h i k k x i x l x i +

 

 

2

 

$ i i k l x t x i x hx i ,

\ < i < ) < k < q

 

 

 

 

 

K i < j < k < q

 

 

 

 

l< l < j < k < l < q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0/ =*//.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V u ^ i y u — Z y t — Z y j ,

 

 

 

 

 

 

 

 

 

 

 

 

 

Y1) = Ъ ( . — У 1 + 2 У ш ) — 2 У Ш ) + У }).

 

 

 

 

 

 

 

ь i )

=

8/з

(

-

/// +

t y u i j — 6 у и

+ Ь у ш )

— У ] ) ,

 

 

 

? l l j k

= 32 (3lU l j k

У \ Ш к

y i j b k )

+ 8/з

(6y . l У )

у к )

—/16 ( У 1) +

 

 

 

+

У i k

— 16/з (5 у ш )

+

5 it u t k

3 y i j j j

S y i k k k

 

y j j j k

y j k k k ) >

 

h j j k

32 (3y i j j k

1l i i j k

 

V i i k k ) + 8/з Ф У ) ~

У 1 У к)

16 (.УЦ + У,]к) —

 

 

 

16/з

Ф У л ) ) )

+ S i j j j j k ?>Уa i ] S y j k k k У ш к

У i k k k ) 1

 

Vljkk

= 3 2 (3iH j k k —

y i l i k

 

y i j j k )

+ 8/з

Ф У к — У l

У ) ) —

16( у ik

+ V j k )

 

 

 

1б/з фу, ikkk + 5 y j kk k 3 y u i k — Z y j j j k

-

У i l l)

y i j j j ) ,

 

h l k i

= 256in j k i

— 32 (1l i i j k

+ U u j i + iи

м

+ y i j j k

+ У Ц ] 1

+ y j j k i

+ y i l k k

+

+

yikkl

+ 1и ш

+

У ti l t

+

y j k i i

+ y i k i i )

+ 32/з {.Уi n ) +

У т ь

+ У m i

+ y i j j j

+

 

 

+ У 1 ) ) к + y j j j l + У I k k k + V j k k k + y k h k l + У l l l l + У J i l l + У k i l l ) •

 

После определения коэффициентов уравнения регрессии необ­ ходимо провести статистический анализ полученных результатов: проверить адекватность уравнения и построить доверительные ин­ тервалы значений отклика, предсказываемые по уравнению регрес­ сии. При постановке эксперимента по симплекс-решетчатым планам нет степеней свободы для проверки адекватности уравнения, так как эти планы насыщенные. Для проверки адекватности ставят опы­ ты в дополнительных, так называемых контрольных точках. Число контрольных точек и их координаты связаны с постановкой задачи и особенностями эксперимента. При этом стараются предусмотреть

9*