книги / Методы оптимизации эксперимента в химической технологии
..pdfфактора на /-м уровне; I,L y2i...jq... |
/ —сумма квадратов |
результатов |
всех N=$?k опытов; |
|
|
55кор = --------- |
~ Ч |
(V. 139) |
— средний квадрат суммы результатов всех N=22k опытов.
Пример 7. Исследовался одностадийный процесс получения водорастворимых поли электролитов путем радикальной полимеризации винилпиридиновых солей без их про межуточного выделения. Процесс зависит от большого числа количественных и каче ственных факторов. Необходимо определить оптимальные условия процесса.
Р е ш е н и е . Для оптимизации процесса синтеза водорастворимых полиэлектролитов на основе 2-метил-5-винилпиридина (2,5 МВП) был использован сложный план —дробная реплика 2 6-2, совмещенная с двумя латинскими квадратами (табл. 51).
|
|
|
Т а б л и ц а |
51. |
Матрица планирования |
|
|
||
Номер |
|
* 2 |
*3 |
х4 |
*5 |
|
*7 |
|
У |
опыта |
|
* 6 |
* 8 |
||||||
1 |
_ |
_ |
_ |
_ |
_ |
+ |
А |
0 |
30,155 |
2 |
+ |
— |
— |
— |
+ |
— |
В |
1 |
31,761 |
3 |
— |
+ |
— |
— |
+ |
- — |
С |
2 |
12,899 |
4 |
+ |
+ |
— |
— |
— |
+ |
D |
3 |
80,715 |
5 |
— |
— |
+ |
— |
+ |
+ |
В |
3 |
88,876 |
6 |
+ |
— |
+ |
— |
— |
— |
А |
2 |
16,806 |
7 |
— |
+ |
+ |
— |
— |
— |
D |
1 |
12,463 |
8 |
+ |
+ |
+ |
— |
+ |
+ |
С |
0 |
81,443 |
9 |
— |
— |
— |
+ |
— |
+ |
С |
1 |
24,959 |
1 0 |
+ |
— |
— |
+ |
+ |
— |
D |
0 |
19,602 |
1 1 |
— |
+ |
— |
+ |
+ ’ |
— |
А |
3 |
21,761 |
1 2 |
+ • |
+ |
— |
+ |
— |
+ |
В |
2 |
52,953 |
13 |
— |
— |
+ |
+ |
+ |
+ |
D |
2 |
78,220 |
14 |
+ |
— |
+ |
+ |
— |
— |
С |
3 |
20,953 |
15 |
— |
+ |
+ |
+ |
— |
— |
В |
0 |
1,684 |
16 |
+ |
+ |
+ |
+ |
+ |
+ |
А |
1 |
48,910 |
Рассматривалось влияние на выход полимера (у) восьми факторов, из которых два
качественных - галогенакилы RX |
и растворители —менялись на |
четырех |
уровнях, |
а один качественный —инициатор —на двух. Значения выбранных |
уровней |
для всех |
|
исследуемых факторов приведены в табл. 52. |
|
|
|
Вследствие ценных комбинаторных |
свойств, плана и рационального |
выбора |
факторов |
уже в процессе реализации матрицы планирования определены условия (опыты 4 , 5 , 8 , 13), в которых выходы полимеров удовлетворяют технологическим требованиям.
Для определения интенсивности влияния различных параметров на выход полимера был проведен факторный и дисперсионный анализы полученных результатов. Ошибка вос производимости 5Воспр “ 1 ,9 4 , используемая в факторном анализе, определена из пред варительных опытов. Число степеней свободы/ ВОспр " 6 .
В табл. 52 приведены эффекты факторов, введенных в планирование на двух уровнях, полученные по формуле (V.133). Значимость этих эффектов проверялась по критерию Стьюдента. Табличное значение критерия Стьюдента г0,0б(6 ) —2,45. Эффект фактора * 2 (соотношение реагирующих компонентов) оказался незначимым. Таким об разом, избыток галоидного алкила не влияет на выход полимера. Незначимый эффект в табл. 52 заменен нулем. Значимость главных эффектов факторов, введенных в план как на двух, так и на четырех уровнях, проверялась при помощи многофакторного дисперсионного анализа. Для оценки значимости эффектов в дисперсионном анализе было использовано отношение средних квадратов, обусловленных действием соответ ствующих факторов, к среднему квадрату,' связанному с ошибкой опыта, имеющее рас пределение Фишера. При этом к сумме квадратов, связанной с ошибкой опыта, отне сена с соответствующим числом степеней свободы сумма квадратов, обусловленная
Н о |
|
|
|
|
|
|
Выход Но |
|
|
|
|
|
|
Выход |
|
мер |
х, |
*2 |
*3 |
*4 |
*6 |
*в |
продук |
мер |
х \ |
*2 |
*3 |
*4 |
*6 |
*в |
продук |
опы |
|
|
|
|
|
|
та у, % |
опы |
|
|
|
|
|
|
та у, % |
та |
|
|
|
|
|
|
|
та |
|
|
|
|
|
|
|
1 |
-1 |
-1 |
-1 |
-1 |
- 1 |
0 |
59,03 |
9 |
- 1 |
-1 |
-1 |
+ 1 |
- 1 |
2 |
60,08 |
2 |
+ 1 |
-1 |
- 1 |
-1 |
+ 1 |
1 |
91,69 |
10 |
+ 1 |
-1 |
- 1 |
+ 1 |
+ 1 |
3 |
87,9 |
3 |
- 1 |
+ 1 |
- 1 |
-1 |
+ 1 |
2 |
70,32 |
11 |
-1 |
+ 1 |
- 1 |
+ 1 |
+ 1 |
0 |
96,29 |
4 |
+ 1 |
+ 1 |
- 1 |
-1 |
- 1 |
3 |
80,73 |
12 |
+ 1 |
+ 1 |
- 1 |
+ 1 |
-1 |
1 |
87,5 |
5 |
- 1 |
-1 |
+ 1 |
-1 |
+ 1 |
1 |
68,14 |
13 |
- 1 |
- 1 |
+ 1 |
+ 1 |
+ 1 |
3 |
77,58 |
6 |
-И |
-1 |
+ 1 |
- 1 |
-1 |
0 |
91,93 |
14 |
+ 1 |
- 1 |
+ 1 |
+ 1 |
- 1 |
2 |
65,32 |
7 |
- 1 |
+ 1 |
+ 1 |
-1 |
-1 |
3 |
50,16 |
15 |
-1 |
+ 1 |
+ 1 |
+ 1 |
- 1 |
1 |
79,59 |
8 |
+ 1 |
+ 1 |
+ 1 |
- 1 |
+ 1 |
2 |
97,09 |
16 |
+ 1 |
+ 1 |
+ 1 |
+ 1 |
+ 1 |
0 |
97,58 |
Дисперсия воспроизводимости, определенная по параллельным опытам, ^ оспр —19,637; /воспр * 16- Ошибка эффекта
Значимость эффектов факторов, введенных в план на двух уровнях, проверялась по критерию Стьюдента для уровня значимости д —0,05 и числа степеней свободы / —16;
/о,о5(16) —2,12.
Значение /-критерия для каждого фактора:
8,659
=7" Q= 7,815 > 2,12 — эффект значим 1»108
3,537 |
= 3,192 > 2,12 — эффект значим |
t2= 7— |
|
1,108 |
|
0,384 |
< 2,12 — эффект незначим |
*8 = Г 108 = |
2 671
1,108 = 2,411 >2,12 — эффект значим
h 7,015
= 7 ™ тя= 6*331 > 2,12 — эффект значим
1 у I UO
Эффект фактора хз оказался незначим, таким образом, избыток растворителя не влияет на выход феллавина. Увеличение времени экстракций, увеличение количества экстоакций и повышение температуры, т.е. оставшиеся четыре фактора приводят к увеличению выхода феллавина.
Значимость различия между эффектами фактора хв, веденном в плане на четырех уровнях, проверялась с помощью множественного рангового критерия Дункана. В таб лице приведены значения среднего выхода для различных уровней фактора хв:
В и д р а с т в о р и т е л я У р о в н и ф а к т о р а С р е д н и й в ы х о д , %
Метанол |
0 |
86,207 |
50%-ный этанол |
1 |
81,73 |
50%-ный изопропиловый спирт |
2 |
73,202 |
50%-ный метанол |
3 |
74,092 |
Проранжируем средние значения выхода в порядке возрастания:
73,202 |
74,092 |
81,73 |
|
86,207 |
у (2) |
J?(3) |
J/(0 |
|
yi°) |
Ошибка воспроизводимости |
результатов |
4,432%. Ошибка среднего равна |
||
*7 = У $/п = |
У 19,637/4 = 2,216%. |
|||
Выпишем значимые ранги из таблицы Дункана |
для |
уровня значимости р —0,05 и |
||
числа степеней свободы nD^ f B0СП|Г |
|
|
|
|
Р . . . . . . |
2 |
3 |
4 |
|
Ранги . . . . |
3,01 |
3,15 |
3,23 |
Наименьшие значимые ранги (НЗР), умноженные н& ошибку среднего sy, равны:
|
Р ................................. 2 |
3 |
4 |
||
|
НЗР X Зу . . . . |
6,67 7,003 |
7,202 |
||
Проверим значимость различия между средними: |
|
||||
JT(°) |
= 86,207 |
— 73,202 = |
13,005 > |
7,202 — различие значимое |
|
^ (3) = 76,207 — 74,092= 12,115 > 7,003 —различие значимое |
|||||
—i^1) == 86,207 |
— 81,73=4,477< 6,67 — различие незначимое |
||||
— у |
= 81,73 — 73,202 = 8,528 > 7,003 — различие значимое |
||||
—.jjT<3>== 81,73 — 74,092 = 7,638 > 6,67 — различие значимое |
|||||
у (3) |
<2> = 74,092 |
— 73,202= 0,89 < 6,67 — различие незначимое |
|||
Из анализа |
результатов |
следует, что |
лучшими |
растворителями являются метанол |
и 50%-ный этанол. Выбираем в качестве растворителя для процесса твердофазной экстракции 50%-ный этанол, который является наиболее эффективным, удешевляет процесс получения препарата и менее ядовит.
Таким образом, оптимальные условия процесса экстракции феллавина из листьев
бархата |
амурского |
получились следующие: число экстракций — 3; время первой экстрак |
ции—6 |
ч; время |
второй и третьей экстракции —6 ч; соотношение сырье : раствори |
тель —1:6; растворитель —50%-ный этанол; температура —80°С (кипение). |
11. Метод последовательного симплекс-планирования. В рассмотрен ных планах типа 2к и 2к~р экспериментальные точки располагались в вершинах многомерного куба. В качестве экспериментального плана можно также использовать регулярный симплекс. Симплексом в /r-мерном пространстве называют выпуклый многогранник, имеющий ровно к + 1 вершину, каждая из которых определяется пересечением к гиперплоскостей данного пространства. Примером симплекса в двумерном пространстве, т. е. на плоскости, служит треугольник. В трехмерном пространстве симплексом будет любая четырехгранная пирамида, имеющая четыре вершины, каждая из которых образована пересечением трех плоскостей —граней пирамиды.
Симплекс называется регулярным, если расстояния между всеми его вершинами равны. Так, регулярными симплексами являются правильный треугольник (двумерный симплекс), тетраэдр (трехмерный симплекс). При планировании экспериментов обычно используют регулярные симплексы. Однако регулярность симплекса, как и на правление градиента в методе крутого восхождения, и свойство
ротатабельности планов не будут инвариантными к масштабу коорди нат факторного пространства. При изменении масштаба регулярный симплекс может стать нерегулярным. С другой стороны, всегда мож но подобрать соответствующее преобразование системы координат, делающее нерегулярный симплекс регулярным.
В экспериментальной практике симплексные планы наиболее широко используются для решения задач оптимизации на стадии движения к почти стационарной области. При этом, чтобы сделать симплекс регулярным, используется линейное преобразование
где . z] - у'-я координата центра плана; Azj—интервал варьирования по у-му фактору.
Для оптимизации используется следующее важное свойство симплек са: против любой из его вершин Aj расположена только одна грань, на которой можно построить новый симплекс, отличающийся от прежнего расположением новой вершины Aj, тогда как остальные вершины обоих симплексов совпадают. Последовательным отбрасыва нием вершин осуществляется перемещение исходного симплекса в факторном пространстве.
Метод последовательного симплекс-планирования состоит в следую щем: планируют исходную серию опытов так, чтобы точки, соответ ствующие условиям проведения этих опытов, образовывали регуляр ный симплекс в факторном пространстве. После проведения опытов выявляется вершина, отвечающая условиям, при которых получаются наихудшие результаты. Далее строится новый симплекс, для чего наихудшая точка исходного симплекса заменяется новой, расположен ной симметрично относительно центра грани симплекса, находящейся против наихудшей точки. Новая точка вместе с оставшимися снова образует регулярный симплекс, центр тяжести которого смещен по сравнению с исходным в направлении: худшая точка —центр тяжести остальных точек. Это направление в общем случае не является наи более крутым, однако оно обращено в сторону повышения качества процесса.
После реализации опыта в дополнительной точке опять произво дится сопоставление результатов, снова выявляется наихудшая точка, которая также заменяется ее зеркальным отражением, и т. д. Шаговое восхождение с последовательным отбрасыванием наихудших точек повторяется до области, близкой к экстремуму.
На рис. 53 показаны схемы достижения экстремума одной и той же поверхности отклика методами крутого восхождения и симплекспланирования. Рассмотрим движение к экстремуму на примере зада чи отыскания наибольшего значения целевой функции двух факторов. Для достижения экстремума методом крутого восхождения (рис. 53, а) в окрестности точки М с известным значением целевой функции был поставлен полный факторный эксперимент 22 (точки 1—4), движение по градиенту осуществлялось в опытах 5—9 до тех пор, пока значение целевой функции не начало ухудшаться. С центром в лучшей точ-
Рис. 53. Симплексный метод и крутое восхождение по поверхности отклика
кё 7 пришлось вновь реализовать план 22 (точки 10—13). Новое движе ние по градиенту (точки 14, 15) приводит к экстремальному значению целевой функции. При использовании симплекс-планирования (рис. 53, б)
в исходном симплексе (точки 1—3) худшей оказалась |
точка 2. Точ |
ка 4 является зеркальным отражением худшей точки |
относительно |
С1—центра грани 1—3. В новом симплексе 1, 3, 4 худшей оказалась точка 1. В результате применения симплексного метода достигли области оптимума (симплекс 9, 10, 11). Таким образом, оба метода потребовали примерно одинакового числа опытов. Из рис. 53 видно, что вблизи оптимума при применении симплексного метода может возникнуть зацикливание. Достигнув области оптимума, симплекс начинает враще ние вокруг вершины с максимальным значением отклика. Если симплекс располагается относительно поверхности отклика таким образом, что значение отклика в новой точке опять получается самым плохим, необходимо вернуться к предыдущему симплексу и попробовать следующее благоприятное направление. Наличие ошибок в определе нии отклика снижает скорость движения к экстремуму.
Исходный симплекс может быть по-разному ориентирован в фактор ном пространству Если центр симплекса совпадает с началом коорди нат, одна из вершин лежит на координатной оси, а остальные располага ются симметрично относительно координатных осей, плоскостей и гипер
плоскостей (в многомерном случае), то координаты вершин симплекса задаются матрицей X:
Х1 |
* 2 |
. •• |
XJ |
. |
Xk |
* h |
|
— *1 |
* 2 |
. . . |
XJ . . . |
|
x k |
||
0 |
— 2 * 2 . . . X) . . . |
Xk. |
Xk |
||||
х = |
|
|
|
|
|
(V. 140) |
|
0 |
0 . . . - j x j |
. . • |
* f c - 1 |
||||
Xk |
|||||||
о |
о |
|
0 |
|
(k I) |
Xfr |
|
о |
о |
|
0 |
|
0 |
—kxh |