Файл: Химмельблау Д. Анализ процессов статистическими методами.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 09.04.2024
Просмотров: 754
Скачиваний: 2
542 |
Глава 8 |
рид, при 228 комбинациях уровней шести независимых перемен ных. Была получена следующая оценка уравнения регрессии:
Т = 349,392 - 2,369 х4 - 3,78 х2 + 6,077 х3 + 7,064 s4 + + l,6795s6 + 3,956 х6 + 3,957z2 + 0,655s2 + 10,019s23 +
+ 7,602 х\ - 0,0428sa + 2,33s2 + 0,564х^ 2 + 0,564st s3 -
- 10,793siS4 + 0,473s!S5 + 0,003^^ + 0,011s2 s3 — 3,852s2s4-r-
+ 0,201s2 s5 — 0,005s2 s6 |
— l,376s3 s4 — 0,010s3 s5 + 0,012s3 s6 |
+ |
+ |
0,942s4 s5 + 2,739s4 x 6 — 1,847 s5 xe . |
(a) |
Для такой регрессии квадрат множественного коэффициента корреляции оказался равным 0,998, однако это уравнение доволь но сложное, и оценить влияние каждой из переменных весьма нелегко. В связи с этим уравнение было обработано на ЭВМ с по мощью программы, специально предназначенной для анализа поверхностей отклика. Уравнение (б) представляет собой полу ченное каноническое уравнение; непосредственно после него запи саны координаты центра системы:
Г - 349,96 = -1,388s2 |
- |
0,064s2 |
+ |
0,604s2 + 1,201s2 + |
|
|
|
|
|
+ |
9,776s2 + 12,049s]; (6) |
координаты центра |
(кодированные переменные): |
||||
sl e |
= |
2,07, |
s4 e |
= |
1,84, |
%2e — 6,57, |
Хье |
= |
4,73, |
||
X 3 e |
— —0,23, S ß e = |
2,87 |
|||
|
|
|
Te |
= 349,96. |
Отметим, что два коэффициента в уравнении (б) отрицательны, а четыре' положительны. Сопоставление с табл. 8.2.1 позволяет сделать заключение, что уравнение (б) представляет седловую гиперповерхность. Рассмотрение координат центра показывает, однако, что он расположен далеко за пределами области экспери мента. Всякий раз, когда возникает подобная ситуация, поверх ность в области эксперимента представляет собой многомерный аналог или наклонного гребня, или наклонной впадины. Более того, поверхность вдоль некоторой оси или всегда поднимается, или всегда опускается; весьма вероятно, что она хорошо описы вается только линейными членами по некоторым переменным.
Из-за удаленности центра системы было решено, что некоторые члены взаимодействия и члены второго порядка в уравнении (а) не дают значимого вклада в отклик. Оказалось, что это дей ствительно так; упрощенное уравнение, полученное в результате
Стратегия |
эффективного |
экспериментирования |
543 |
|
нового регрессионного |
анализа, имело вид |
|
|
|
f = 356,57 - 4,0975а;2 |
- 10,5077*1*4 - |
2,445*! + |
0,6509*2 |
- |
-3,8418*2 z4 — 3,7421*2 -f- 6,1038*3 + 7,5876*2 + |
|
|||
|
+ 7,1016*4 |
+ 1,706*6 |
+ 3,945*в; |
(в) |
множественный коэффициент корреляции оказался равным 0,923.
Уравнение (в) показывает, что зависимость температуры от пере |
|||
менных *з, * 5 |
и * 6 чисто линейная. |
под |
|
На этом |
этапе |
анализа эти три переменные временно |
|
держивались |
равными нулю, чтобы несколько упростить |
ана |
|
лиз. Температура t |
исследовалась как функция оставшихся |
трех |
переменных * ь * 2 и * 4 . Центр новой системы координат |
находился |
||||
в точке |
|
|
|
|
|
|
Хіе |
= |
— 0,643, |
|
|
|
* 2 |
е |
= |
0,707, |
|
|
* 4 |
е |
= |
— 0,734. |
|
Получилось новое |
каноническое уравнение |
|
|||
f — 353,426 = |
- 0J477* 2 + |
11,6025*2 + 1,4812 |
(г) |
где преобразованные переменные связаны с исходными следующим образом:
* і |
= |
0,538 |
(*! + |
0,643) |
+ |
0,681 (*2 - |
0,707) + |
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
0,496 (*4 |
+ |
0,734), |
|
* 2 |
= |
0,568 |
(*! + |
0,643) |
+ |
0,142 (*2 |
- |
0,707) - |
0,811 (*4 |
+ |
||
|
|
|
|
|
|
|
|
|
|
+ 0,734), (д) |
||
*4 = |
0,623 |
(*! + |
0,643) - |
0,718 (*2 |
- |
0,707) + |
0,310 (*4 |
+ |
||||
|
|
|
|
|
|
|
|
|
|
|
+ |
0,734). |
|
В уравнении (г) один член имеет отрицательный коэффициент, |
|||||||||||
а |
два других — положительные; снова получили |
уравнение сед- |
||||||||||
ловой поверхности. Однако заметим, |
что теперь |
центр |
системы |
лежит внутри экспериментальной области и максимальные значе ния канонических переменных * 1 ? * 2 и * 4 представляют собой небольшие величины одного порядка малости. Следовательно, относительная малость коэффициентов при первом и третьем чле
нах уравнения (г) позволяет отбросить эти члены, |
что |
не внесет |
заметной ошибки в рассчитываемый отклик Т. Поэтому |
упрощен |
|
ное уравнение будет иметь вид |
|
|
t - 353,4 = 11,6025 х\, |
|
(е) |
544 |
Глава 8 |
или
х2 |
= ] / |
ГЛ - 3 5 3 , 4 |
(ж) |
|
|
11,60 |
|
Подстановка сюда выражения (д) для х 2 дает уравнение
0,568x^0,142x2—0,811х4 —0,331 = | / Г ~ ^ ' 4 ,
описывающее семейство плоскостей для различных температур Т, имеющих минимальное значение 353,4 °F. Геометрически это
сс^ (скорость подачи)
-Ге=20
х 2 (глубина канала)
(скорости червяка..)
Ф и г . П . 8 . 2 . 3 . Представление системы с помощью упрощенного уравнения
|
|
|
|
(з). |
|
|
|
|
|
|
|
|
Б |
Вершины |
|
|
|
|
|
+1 |
(0,14 |
дюйма) |
—1 (0,10 |
дюйма) |
С |
+1 |
(0,14 |
дюйма) |
|
—1 (50 |
фунт/ч) |
|
—1 (50 |
фунт/ч) |
|
+1 |
(150 |
фунт/ч) |
|
+1 |
(60 |
об/мин) |
D |
—1 (30 |
об/мин) |
|
+1 |
(60 об/мин) |
|
|
|
|
+1 (0,14 |
дюйма) |
|
|
|
|
|
|
|
|
|
—1 (50 |
фунт/ч) |
|
|
|
|
|
|
|
|
—1 (30 |
об/мин) |
|
|
|
|
уравнение представляет стационарный гребень с параллельными плоскостями, распределенными вокруг плоскости, соответствую щей минимальному отклику (фиг. П.8.2.3).
Вспомним, что уравнение (з) и фиг. П.8.2.3 дают лишь упро щенное представление о первоначальной поверхности отклика,
Стратегия |
эффективного |
экспериментирования |
545 |
которая описывается уравнением с шестью независимыми перемен ными. Дл я того чтобы учесть влияние переменных х3, хъ и х6, положенных в уравнении (е) равными нулю, к уравнению (е) были добавлены линейные члены, исключенные из уравнения (в), что дало
Т |
- 353,4 = 11,6025^ + 6,1038^3 + |
1,706х8 |
+ 3,945 хв. |
|
(и) |
|
|
Уравнения (з) и (и) приводят |
к |
|
|
|
|
0,568^ + 0,142^—0,811ж4 —0,331 |
= |
|
|
|
|
|
|
Т — (353,4 + 6,1038ж3 + 1 , 7 0 6 ж 5 + 3,945д:6 ) |
. |
||||
|
/ |
|
ÜT6Ö |
|
• |
( к > |
|
Предварительно было показано, |
что исключенные |
члены х3, |
|||
хъ |
и хв линейны и соответствуют |
наклонной |
впадине или гребню |
на поверхности отклика. Так как все три члена входят в уравнение с положительными коэффициентами, то для минимизации отклика необходимо положить их равными минимальным значениям,
а именно |
— 1 . Это приведет к уменьшению члена в круглых скоб |
|||||||
ках, |
модифицированного |
значения |
Те, от 353,4 до 341,6° F и вы |
|||||
зовет |
дополнительное |
уменьшение |
значения отклика |
на 11,8° F. |
||||
Однако |
уравнение (к) |
все еще описывает семейство |
плоскостей |
|||||
с параметром Т, и поэтому фиг. П.8.2.3 дает точное |
представление |
|||||||
о полной |
поверхности |
отклика с |
минимальным |
значением Те, |
||||
теперь |
равным 341,6° F. |
|
|
|
Следует отметить, что проведенный анализ справедлив лишь для исследованной экспериментальной области, и нет никаких оснований считать, что поверхность отклика будет такой же за пределами этой области.
8.3.М Е Т О Д О П Т И М И З А Ц И И ПРОЦЕССО В
ПО С Р Е Д С Т В О М Э К С П Е Р И М Е Н Т И Р О В А Н И Я
Вэтом разделе рассматриваются вопросы оптимизации устано вившихся процессов с помощью некоторой последовательности планируемых экспериментов, проводимых с целью улучшения характеристик процесса. Процесс оптимизации в трех измерениях напоминает подъем на холм. Хотя целевые функции, содержащие две независимые переменные, весьма редко имеют форму холма, изображенного на фиг. 8.3.1, тем не менее общая процедура заклю чается в том, чтобы, исходя из некоторой точки (набора значений независимых переменных) на холме, подняться на его вершину (или спуститься на дно впадины в процессе минимизации). Необхо димо решить, каким путем идти «наверх» и как, учитывая опреде ленные ограничения, достичь вершины наиболее эффективным спо собом. Важное различие между детерминистической и экспери-