Файл: Химмельблау Д. Анализ процессов статистическими методами.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.04.2024

Просмотров: 754

Скачиваний: 2

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Стратегия эффективного экспериментирования 541

Новые координаты выражаются через первоначальные следую­ щим образом:

 

 

г

1

1

о "

 

 

 

 

Ѵ'2

Ѵ 2

 

 

 

 

 

 

 

и т х

=

1

1

2

 

 

Уе

уе

Уе

 

 

 

 

1

1

i

 

 

 

Іуз

Уз

ysJ

 

Каноническая форма уравнения (а) такова:

 

Y

= 2х\

+

Ъх\ +

Линейные члены.

(г)

Члены первого

порядка

получаются

подстановкой xt,

х2 и х3

в каждый линейный член уравнения (а), причем используется преобразование х = Ux, или

1

У б #2 " Ï7I*3 '

Уз Ж 3 і (д)

1 ~

у!*3 '

Пример 8.2.3. Канонический анализ для экструдера пластмасс

Следующий пример взят из работы [11]. При исследовании про­ цесса прессования пластмассы были получены данные по компози­ ционному плану с шестью независимыми переменными для подгон­ ки уравнения регрессии второго порядка, в котором зависимой

 

 

Таблица

П.8.2.3

 

Планируемые значения

Переменная

кодированные значения

 

- 1

0

+1

Скорость

червяка

х ь

об/мин

30

45

60

Глубина

канала

х2,

дюйм

0,100

0,120

0,140

Температура в камере xs,

°F

330

340

350

Скорость

подачи

xt,

фунт/ч

50

100

150

Температура на

входе хъ,

°F

280

310

340

Длина

измерительной

зоны хв,

10

19

28

дюйм

 

 

 

 

 

 

 

переменной служила температура выходящей массы Т (см. табл. П.8.2.3). Расчеты проводились для экструдера с диаметром 2,5 дюйма, использующего пластифицированный поливинилхло-


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, тем не менее общая процедура заклю­ чается в том, чтобы, исходя из некоторой точки (набора значений независимых переменных) на холме, подняться на его вершину (или спуститься на дно впадины в процессе минимизации). Необхо­ димо решить, каким путем идти «наверх» и как, учитывая опреде­ ленные ограничения, достичь вершины наиболее эффективным спо­ собом. Важное различие между детерминистической и экспери-