Файл: Кориков А.М. Математические методы планирования эксперимента учеб. пособие.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 ;

 

 

 

 

 

д х,

д х,.