Файл: Кориков А.М. Математические методы планирования эксперимента учеб. пособие.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 27.06.2024
Просмотров: 168
Скачиваний: 0
|
Если |
\j ( х г |
, |
х а |
) |
имеех |
минимум, ю |
максимум |
tj. ( |
х г |
, |
х ^ |
|||||||||||
следует искать в вершинах квадрата ограничений, го есть в |
|
||||||||||||||||||||||
точках: |
( |
I . l ) |
, |
|
( i . - D , |
|
|
+ 1 ) > |
( . Х | |
_ I |
} > |
И з |
з н |
а ч е н ] |
й |
|
|
||||||
Ч |
t 5 г » |
в |
э |
т |
и х |
|
" ч к а х |
выбирается наибольшее. |
Зю |
и |
бу |
||||||||||||
дет |
максимальное |
значение |
£ |
( х , , ^ , х^ |
на грани |
х , |
» |
+ 1 . |
|||||||||||||||
|
Если |
( Х у ) |
|
внутри области |
|
|
|
|
|
|
,• |
|
|
||||||||||
ь= |
2,3 |
имеет |
максимум, |
то это максимальное |
значение |
и |
|
бу |
|||||||||||||||
дет |
решением на |
грани |
х , |
|
* |
+ 1 . |
|
|
|
|
|
|
|
|
|
|
|||||||
|
В остальных случаях необходимо максимум ^ |
( ссг , х 3 ) |
|
ис |
|||||||||||||||||||
кать на поверхности квадрата ограничений |
|
- I s |
Х ^ £ - И |
|
|
, |
|||||||||||||||||
Ъ= 2 , 3 . Фиксируем |
|
осг |
= |
+1 и ищем экстремум |
i | ( |
х ч ) |
|
|
|||||||||||||||
при наличии ограничения |
- |
1 |
* х & 4 + 1 |
|
.Решение |
этой |
|
од |
|||||||||||||||
номерной |
задачи ухе |
не |
представляет |
труда. Запоминаем получен- |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
л |
|
|
|
|
|
|
|
|
|
|
|
ное |
максимальное |
|
значение |
tj |
и соответствующие |
ему |
значения |
||||||||||||||||
координат. Затем |
последовательно фиксируем |
х |
|
= - |
I |
, х |
* + 1 , |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
л |
|
X j = |
- I , отыскиваем |
и запоминаем максимальное |
значение |
^ . |
|||||||||||||||||||
Из полученных четырех максимальных |
значений |
£ |
выбираем наи |
||||||||||||||||||||
большее |
- |
решение |
задачи |
на |
грани |
х , » + I . |
|
|
|
|
|
|
|
||||||||||
|
2 . фиксируем |
|
х , = - I и вновь повторяем все действия, |
|
|||||||||||||||||||
описанные в предыдущем пункте I . |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
Затем |
последовательно |
исследуем,грани |
х г |
» |
+ 1 , х ^ » - |
I |
, |
|||||||||||||||
X j = |
+ 1 , |
х 3 = |
|
- |
I . Окончательно из решений на |
каждой |
грани |
||||||||||||||||
выбираем |
то, которое |
соответствует |
наибольшему |
значению |
|
|
|||||||||||||||||
у |
( |
I |
) . |
Заметим, |
что этот |
способ |
отыскания экстремальных |
||||||||||||||||
значений |
функции |
|
^ |
( х ) |
можно применять |
и без |
предваритель |
||||||||||||||||
ного |
приведения |
|
ее |
к каноническому |
виду. Для |
предварительной |
- 1 2 2 -
классификации поверхностей оптика в отом случае следует ис пользовать определители Сильвестра, полученные иг матрицы
соли все определителиСильвестра положительны
|
|
" |
1 1 |
^ |
1 2 . |
> 0 , - д . - |
|
|
|
||||
то • ^ |
( х ) имеет миаиму_1} |
|
|
|
||
если нечетные определители Сильвестра отрицатзльны, а |
||||||
чехные |
положительны |
|
|
|
|
|
то |
^ х ) имеет |
максимум |
|
|
|
|
в |
остальных |
случаях |
tj |
( |
х ) это поверхность т-ипа мини- |
|
макс, |
|
|
|
|
|
|
4 - 5 . * , П р и м е р ы В качестве первого примера исследуем поверхность откли
ка, описываемую уравнением ( 4 . 4 ) . Для задачи с двумя факто рами уравнение ( 4 . 6 ) примет -следующий вид
Яосле подстановки в это уравнение значений коэффициен
тов регрессии получим
Ъг~ |
1,41-5» - 5,544 = 0 . |
Откуда
Так как канонические коэффициенты имеют разные знаки, поверхность отклика является гиперболическим параболоидом.
Координаты центра поверхности найдены из системы уравне
ний
После подстановки значений коэффициентов получим
5,2 |
x t , - t 3 , 0 i e |
+ ' 3 , 4 5 - 0 ; |
|
ь,о |
z w |
+ г,ьбх 0 & - 1, 32- 0 . |
|
Откуда |
|
|
|
х и |
- - |
0,195, |
х 6 У - - 0 , 8 . |
Зная координаты центра, рассчитаем по формуле ( 4 . 4 ) со |
|||
ответствующее им значение параметра оптимизации, которое |
|||
оказывается |
равным 85,34 . Отсюда уравнение ( 4 . 4 ) в канони |
ческой форме приобретает вид
$ - 65, 54 - 5 , 1 2 I * |
- 1,72 X J . |
( 4 , 1 2 ) |
|
|
|
|
|
|
угол поворота новой системы координат относительно ста |
||
рой |
равен 19,2° . |
|
|
|
Для определения оптимальных режимов воспользуемся мето |
||
дом |
"ридх анализ". |
|
|
|
Вычисляем значение |
параметра Юрля по формуле (4 . 10) |
X ' « 2 ( £ „ - 6 J - 2 ( 3 , 1 * + 1 , 1 9 ) « б , 62.
- 1 2 4 —
Отсвда допустимые значения X лежат в пределах
8,62 > А > 3 , 1 2 .
Затем задаенся несколькими значениями Л. , постепенно приближаясь х величине, равной 3 , 1 2 , и вычисляем по форму
лах ( 4 . I I ) оптимальные режимы.
Пусть X ш 4 , 0 , тогда
_ - |
0,25 • 3,0 • 1,52. + 0,5- 0,43(1,19+4,0) |
1 |
^ . |
|
|
U ,6 - 4,0)С - 1,19 - 4,0) - 0,25(5,0) г |
|
|
|
_0,&5 6 1 S & , - |
0,5 о » ^ - 7 1 ) ^ |
|
|
|
^ |
(b„-\)(ibi-\)- |
0,25 о* |
|
|
|
0,25 - 5,0 • 5 . 4S + 0,3 - %Ы ( 2 , 6 0 - 4 , 0 ? |
Q |
^ |
|
|
( 2 , 6 - 4 , 0 ) ( - 1 , 1 9 - 4 , 0 } - 0 , 2 5 ( М ) г |
|
|
После |
подстановки |
значений |
х , |
= 1,577 и |
хъ » 0,324 в |
||
уравнение |
( 4 . 4 ) находим, что при этом режиме |
параметр оп |
|||||
тимизации равен 9 6 , 8 5 |
$ . |
|
|
|
А |
||
Аналогично вычисляем значения |
i , , i , |
i |
|||||
ij , задавшись |
|||||||
X - 3,95 и 3,90 |
|
|
|
|
|
||
X |
ш 3,95 ; |
. |
% - 3,90 |
; |
|
||
. . х, = 1,669 ; |
|
х,= 1,827 |
; |
|
|||
х 3 = 0,359 ; |
|
X j - 0,405 |
; |
|
|||
/ } - 9 9 , 2 8 ; |
|
9 = 1 0 1 , 6 0 . |
|
— 125 —
Дальнейшее уменьшение величины неопределенного множите
ля Лагранка будет вызывать дальнейший рост величины парамет ра оптимизации. Наряду с этим значения параметра оптимизации
больше 100$ не имеют смысла. Это значит, что области фактор
ного пространства, соответствующие значениям неопределенного
множителя Лагранжа в пределах |
от 3,12 до 3 , 9 0 , по-видимому, |
|||||||
плохо |
описываются |
уравнением ( 4 . 4 ) . |
|
|||||
|
Во втором примере мы продемонстрируем применение метода |
|||||||
перебора |
граней гиперкуба ограничений. |
|
||||||
|
Пусть |
уравнение |
регрессии |
имеет следующий |
канонический |
|||
ШЛу |
- |
5 0 , 9 |
4 1 V |
|
4,270 X ? |
+ 4, 5в7 X * |
|
|
Видно, |
что £ |
( х, |
, |
х а ) это поверхность типа минимакс. |
||||
|
1) |
Полагаем |
х 1 |
= +1 и в сечении получаем |
параболу |
J |
U * b 4 |
i ' 9 5 1 |
+4,592, х г + 0,216 х\ |
, |
|
||||||||
которая |
максимальное |
значение |
принимает |
при |
x t - I : |
||||||||
£ ( 1 , 1 ) - 4 8 , 5 6 5 . |
|
|
|
|
|
|
|
|
|||||
2 ) |
При |
х^ * |
- I |
функция |
|
|
|
|
|
||||
% |
ij t i j ) - |
50,761 + t , 1 H x t |
+ 0,216 |
х\ |
|
||||||||
и имеет |
максимум |
при |
|
х г . |
+1 |
: |
^ ( - 1,1) = |
52,1165. |
|||||
3 ) |
При х г |
= |
+1 |
функция |
|
|
|
|
|
||||
|
$ ( . х , ) « 51,534 |
- |
1,777 х , - |
1, 195 |
х * |
|
|||||||
имеет максимум |
в |
точке |
х , |
* |
- 0 , 7 4 1 , |
то есть |
• |
||||||
^ |
(-0,741? |
I ) |
- 52,217 . |
|
|
|
|
|
|
- 1 2 6 -
4 ) И, наконец, в сечении хг = - I
^ ( х , ) * 45,024 - 5 , 0 5 2 х , - -1,195
максимальное значение достигается при х , = - I и равно
£ ( - 1 , - 1 ) = 4 9 , 8 8 1 .
Из полученных результатов видно, что наибольшее значение
параметра оптимизации |
= 52,217 достигается в точке |
фак |
||
торного пространства с |
координатами хл = - 0 , 7 4 1 ; |
Х г = |
I . |
|
§ 4 - 6 . О т ы с к а н и е |
у с л о в н о г о |
э к с т |
||
р е м у м а |
п р и |
н а л и ч и и |
н е с |
|
к о л ь к и х |
п о в е р х н о с т е й |
|
о т к |
ли к а
Впрактической работе часто приходится отыскивать услов
ный |
экстремум функции |
отклика |
«• у ( |
( х , |
, хг |
, . . . , х ^ ) |
||||||
при |
ограничениях,накладываемых |
другой |
функцией |
|
|
|
||||||
угтЧг |
( х , , х г , . . . , |
Хк) |
|
.Примером |
такой |
задачи |
||||||
моке^т быть определение режима получения продукта |
заданного |
|||||||||||
состава |
при максимальном |
выходе. Обозначим |
черзз |
^ |
выход |
|||||||
продукта, а через |
у 4 |
его |
качество. При k |
= 2. эта |
задача |
|||||||
решается просто графически.' При большем числе |
факторов за |
|||||||||||
дачу решают с помощью метода неопределенных множителей |
||||||||||||
Лагранка. |
|
|
|
|
|
|
|
|
|
|
||
|
Метод неопределенных |
множителей Лагранжа сводится ж |
||||||||||
решению |
системы |
уравнений |
|
|
|
|
|
|
|
|||
|
|
— — |
+ |
\ |
|
=* 0 ; |
|
|
|
|
|
д х, |
д х,. |