
книги из ГПНТБ / Кориков А.М. Математические методы планирования эксперимента учеб. пособие
.pdf-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 |
4ъ 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. в направлении от параметра Хорля к соответствующему каноническому коэффициенту.
На практике обычно задаются несколькими значениями не определенного множителя Лагранжа и затем, после решения