Оптимизация эксперимента в химии и химической технологии
..pdfтив из каждой вершины треугольника высоту, разделив к аж дую из них на десять равных по величине отрезков и проведя через полу ченные деления прямые, параллельные сторонам треугольника, по лучим треугольную сетку. Приближение от каждой данной стороны к противоположной вершине отвечает пропорциональному возра станию содержания соответствующего компонента, поэтому пря мые, параллельные данной стороне, при последовательном перехо де от одной прямой к другой отражают возрастание третьего компонента на 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 - |
|
для ^-компонентной
|
У = |
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)
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) |
Аналогично |
выводятся |
соотношения для |
определения коэффи |
циентов любого приведенного полинома при любом числе компо нентов. Так, для полинома неполного третьего порядка для трех компонентной смеси имеем:
Iу = Pi*! 4- Рг*2 4- рз*з 4- Pi2*i-*2 4" Pi3*i*3 4- Ргз-^г-^з 4- Pi23*i*2*3 >
Р/ = |
У\ и т. д. |
|
( V I . 39) |
|
Pi2 = 4 //i2— 2//1 — |
2у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 = 4у 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) > |