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

книги из ГПНТБ / Кориков А.М. Математические методы планирования эксперимента учеб. пособие

.pdf
Скачиваний:
37
Добавлен:
23.10.2023
Размер:
7.9 Mб
Скачать

-110 -

 

 

Г и с . ИЛ

 

./•игральную

точку

тской поверхности

называют с е д л о м -

или н и н к

и а

к с о м, так как в

направлении одних к а ­

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

Если один из коэффициентов, например Ъ,0 , равен нудв, то под определение центра фигуры здесь подходит любая точка на оси Х г ( - значение параметра оптимизации в любой точке на оси З С а ) . Этот тип поверхности можно рас­ сматривать как предельный тля первых двух случаев, когда

протяженность

по оси

Xстановится бесконечной. Поверх­

ность отклика

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

собой

с т а ц и о н а р н о е

в о э в i а е в и - е ..

 

 

Наконец,

возможен и такой случай, когда один из коэф­

фициентов, например

!Е>г 2 ,

равен

нулю, но центр фигуры на

ходится на бесконечности, перенеся начало координат в ка­ кую-нибудь подходящим образом выбранную точку 0' вблизи центра, получаем уравнение параболы

где & г - крутизна наклона возвышения. Зтот тип поверхностч опять-таки мохно рассматривать как предельный для пер­ вых двух случаев, когда центр фигуры отнесен на бесконеч­ ность. Локально говерхность отклики здесь представляет со­ бой возрастающее розвышение ( " г р е б е н ь " ) .

Прг обработке результатов реального эксперимента вряд ли можно встретиться с последними двумя предельней слу­ чаями. Практически возможен случай, когда це^тр фигуры удален за пределы той Сласти, ь которой варьировались U-J- ременнне, а один из коэффициентов, ЗЦ, или З Ь г 2 , мало отличается от нуля. Тогда поверхность отклика, в зависи­ мости от наклона возвышения, будет аппроксимировать стаци-

 

 

 

 

 

 

 

 

- 1 1 2

-

 

 

 

 

 

 

онарньш или возрастающий

возвышением.

 

'

 

 

 

 

Для

приведения

зависимости

у

( х

) Б каноническому

виду можно воспользоваться достаточно простым

 

м е т о ­

д о м

 

Д а г р а н ж а .

Продемонстрируем

идею этого

мето­

да для

случая

трехмерного

факторного пространства.

 

С

помощью

обозначений

 

 

 

 

 

 

 

 

 

 

1 .

 

х =

 

С - 1 с

 

T " t 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i t

я

±. Ь

*t

t

is

 

 

 

 

 

 

 

 

 

 

 

z

г.

 

 

причем,

c n - c 4 1 ,

c

„ «

C

n ,

 

с , ь

- с и

,

 

 

 

 

 

уравнение регрессии

можно

записать

в следующем

виде :

 

Приведем

к каноническому

виду

квадратичную

форму

 

х С х

.

Если

в матрице

С

элемент

с

0

 

, то

пере­

ходим

к новый

переменным

" t , ,

Хг,

 

в соответствии

с ли­

нейным

преобразованием

 

 

 

 

 

 

 

 

 

 

 

Ч- *i i

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

х 1 = 7* г

1 ~ с

* с

1 '

11

11

t „

 

v>.,

fc„

 

или в матричной форме

где

ис-it

1О

1 t .

Вновых переменных квадратичная форма приобретает слезу* ющия вид:

где

 

| А „

0 ' 0 «

S 5 - J t T C J L - | 0

d

Так как обычно ctu =£0 , то применим ене раз линейное

преобразование

которое в матричной форме имеет следуюищи вид;

 

- м -

 

где

 

 

0

0

 

1

 

 

0

1

 

Запишем квадратичную форму в координатах

, г г , г э

где

Тем самым мн привели квадратичную форму к каноническому виду

хт

С а& =» а., г.

Итак,

имеем

V

 

где l T = ( C 1 t t t , i 4 ) « & T J l r \

Отсюда окончательно можно записать канонический вид для . всего .уравнения регрессии

1-1

—115 -

где

Видповерхности отклика устанавливаем по знакам канони­

ческих коэффициентов. При

& ц >

0

,

t «

1,2,3,

§

( х"

)

имеет

минимум. Если

3biv-c

0 ,

i

«

1,2,3, то

£

( х

)

 

имеет

максимум. Если

vB-b i i

имеют различные

знаки, то

у

(*) -

- поверхность типа мииимахс.

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

Если поверхность отклика - гиперболический параболоид, то определение оптимальных режимов усложняется,

Рассмотрим три метода внбора оптимальных режимов.

 

-

 

ч - 5 . 1 . П е р в ы й

м е т о д - д в и ж е н и е

в д о л ь

к а н о н и ч е с к и х

о с е й

В соответствии с поставленной задачей в новой системе координат выбирают ось, вдоль которой параметр оптимизации меняется в желаемом направлении и с максимальное скоростью (канонический коэффициент имеет соответствующий знак и мак­ симален по абсолютной величине). Затем, задаваясь различными значениями параметра оптимизации, вычисляют соответствующие им режимы и подвергают их опытной проверке. В связи с сим­ метрией поверхности отклика каждому'значению параметра оп­ тимизации соответствует два режима.

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

Положим, что в задаче требуется определить условия по­

лучения продукта, характеризуемые максимальным значением

 

Л

 

параметра

оптимизации. Пусть if

>t(.f l , а движение осуществ-

х ) §наки

коэффициентов выбраны

произвольно.

 

 

 

-

1П -

ляется

вдольоси

X , . При этом

Х £ = о и формул?'-(4.7) при­

обретает следующий

вид:

 

 

 

«а

- \

x f

 

откуда

 

 

 

 

Ч

Переход от координат Х - к координатам х^, т . е . к коди­ рованным значениям факторов, осуществляют по следующим фор­ мулам перехода:

x 1

= ( X , + x 0 , ) c d s Y - ( x i + x 0 2 , ) S l n , Y ;

( 4 * 8 )

x t

«

( Х , + х в 1 ) s l u y

+ i l j + x

j

соь-у;

JCJ и новой

где л|> - угол между старой

системой

координат

системой координат Х- ,

 

 

 

 

Значения коэффициентов

в формулах

перехода

( s b n f и cosy)

вычисляют

следующим образом:

 

 

 

 

 

л

гг

 

 

 

bin, \ J - 0,5 ( • ) - c o s y ) ;

cos у — 0 , 5 ( 1 + c o s y ) .

В рассмотренном выи1е случае Х г « О, поэтому формулы перехода несколько упрощаются

х , = ( I , + x 0 , ) с о ь у - i s i n , y •, x t « ( X , + x 0 1 ) b i t t y + x o a c o s y .

При решении задачи на минимум параметра оптимизации за­ даются значениями ^ 0 и движутся вдоль оси, соответст­ вующей отрицательному.каноническому коэффициенту. В этой слу­ чае из зависимостей ( 4 . 7 ) и ( 4 . 8 ) получаем

•х, - х и c o s y - ( X l + x o a ) bin \f ,

После расчета кодированных значений факторов для перехо­ да к их натуральным значениям используют формулу кодирования

(см.§ I - I ) .

 

 

 

 

 

 

 

4 - 5 . 2 . В т о р о й

 

м е т о д - "р и д и

а н а л и з "

"Ради анализ" базируется на методе неопределенных множи­

телей ЛагранжА Для выбора оптимальных режимов составляет

 

следующую систему

уравнений:

 

 

 

 

( " Ь я - > • ) * , +0,5

* и. Ч

+

- + 0 , 5 * 1 к х 1 ь

+

0 , 5 6 , - 0 ;

 

• 0 . 5 ^ * 1 + ( U t t r > . ) x t + . . . + 0 , 5 о а х к

+ 0 , 5 D , = 0;

 

0,5 Ь„хл + 0,5-Ьи ж,

4

...

+ 0 , 5 & s K x K

. + 0,5од =0;

( Д . 9)

°>5 4i х,+ 0 , 5 * К 2 х ,

+ . . . .

+ С б к К - г ) х к +

О,5бк -0;

 

где X - неопределенный множитель Лагранжа,

Количество уравнений в системе ( 4 . 9 ) равно числу факто­

ров. Реаение системы ( 4 . . 9) может быть осуществлено только

при заданных значениях X . Выбор значений неопределенных

множителей Лагракка зависит от типа задачи. В случае задачи

на максимум параметра оптимизации рекомендуется выбирать ве­

личину

Л. таким

образом, чтобы она была больше максимального

из

канонических

коэффициентов. В случае

задачи на минимум

выбранное

значение \

должно быть меньше наименьшего из

ка­

нонических

коэффициентов. На величину

Л. накладывается

огра­

ничение, определяемое

параметром Хорля

 

 

 

 

 

 

 

мим

 

 

( * ; 1 0 )

 

 

 

 

 

 

(в зависимости от

где

бЬмоксмаксимальный или ыыкмальный

 

 

мин

 

 

,

 

 

 

задачи)

канонический

коэффициент!

-

коэффициент регрес­

сии

при

Я. -том

квадратичном члене,

 

 

 

использование параметра Хорля дает возможность сузить

интервал изменения значений неопределенного множителя Лагран­

жа до величины,

определяемой

следующим неравенством:

I ? t ' 1

^. "К. >

j j 6

MOKft

 

I

мин I

При этом изменение параметра оптимизации в желаемув сто­ рону соответствует измзнению X. в направлении от параметра Хорля к соответствующему каноническому коэффициенту.

На практике обычно задаются несколькими значениями не­ определенного множителя Лагранжа и затем, после решения

Соседние файлы в папке книги из ГПНТБ