Економіко-математичне моделювання

Економіко-математичне моделювання Економіко-математичне моделювання

library.tneu.edu.ua
from library.tneu.edu.ua More from this publisher
19.07.2013 Views

для цієї групи методів є градієнтні. Непрямі методи полягають у зведенні задачі до такої, знаходження оптимального розв’язку якої вдається спростити. Найпоширенішими методами цього класу є методи квадратичного програмування. Оптимізаційні задачі, на змінні яких накладаються обмеження, розв’язуються методами класичної математики. Оптимізацію з обмеженнями-рівностями виконують методами зведеного градієнта, наприклад, методом множників Лагранжа. У задачах оптимізації з обмеженнями-нерівностями досліджують необхідні та достатні умови існування екстремуму Куна-Таккера. Розглянемо метод множників Лагранжа на прикладі такої задачі нелінійного програмування: , ,..., ) max(min) , x x x f Z = (6.3) ( 1 2 n q ( x , x ,..., x ) i 1 2 n bi ,..., , ( x x x i 1 2 n = , i = 1, m, (6.4) де ( , ,..., ) 1 2 n x x x f та q ) – диференційовані. Ідея методу Лагранжа полягає в заміні окресленої задачі простішою – знаходження екстремуму складнішої функції, але без обмежень. Ця функція називається функцією Лагранжа і записується у вигляді: L( x , x ,..., x ; λ , λ ,..., λ ) = f ( x , x ,..., x ) + λ [ q ( x , x ,..., x ) − b ], (6.5) 1 2 n 1 2 m 1 2 де λ – невизначені поки що величини, так звані множники Лагранжа. i Необхідною умовою екстремуму функції багатьох змінних є рівність нулю частинних похідних стосовно всіх змінних функції. Візьмемо ці частинні похідні і прирівняємо їх до нуля: ⎧ ∂L ⎪ = 0, j = 1, n, ⎪∂x j ⎨ ⎪ ∂L = 0, i = 1, m, ⎪⎩ ∂λi або m ⎧∂f ( x1, x2,..., xn ) ∂qi ( x1, x2,..., xn ) ⎪ + ∑ λi = 0, j = 1, n, ⎨ ∂x j i= 1 ∂x j ⎪ ⎩qi ( x1, x2,..., xn ) − bi = 0 i = 1, m. Отримуємо систему (m+n) рівнянь із (m+n) невідомими, * розв’язавши яку, знайдемо X = ( x , x ,..., x ) та λ = ( λ , λ ,..., λ ) – 1 2 n 0 1 2 m стаціонарні точки. Оскільки їх знайдено з необхідної умови 192 n m ∑ i= 1 i i 1 2 n i

екстремуму, то в них можливий максимум або мінімум. Іноді стаціонарна точка є сідловою (точкою перегину графіка функції). Теорема. Нехай в околі критичної точки (x0; y0) функція f(x, y) має неперервні частинні похідні до другого порядку включно. Розглянемо вираз ( ) ( ) ( ) ( ) 2 2 2 2 ∂ f x, y ∂ f x, y ⎛ ∂ f x, y ⎞ D x, y = ⋅ − ⎜ ⎟ . 2 2 ∂x ∂y ⎝ ∂x∂y ⎠ 1) Якщо D ( x , y ) > 0, то в точці (x0; y0) функція f(x, y) має такі 0 0 екстремуми: 2 ∂ f ( x , y ) 0 0 а) досягає свого максимального значення, якщо < 0; 2 ∂x 2 ∂ f ( x , y ) 0 0 б) досягає свого мінімального значення, якщо > 0; 2 ∂x 2) Якщо D ( x , y ) < 0 , то в точці (x0; y0) функція f(x, y) не має 0 0 екстремумів. 3) Якщо D ( x , y ) = 0, то в точці (x0; y0) для функції f(x, y) 0 0 екстремуми можуть існувати чи не існувати. Для наглядності в деяких випадках цю теорему формулюють інакше. Складаємо матрицю такого виду: 2 2 ⎡ ∂ Z ∂ Z ⎤ ⎢ 2 ∂x ∂x∂y ⎥ H ( x, y) = ⎢ 2 2 ⎥ . ⎢ ∂ Z ∂ Z ⎥ 2 ⎢⎣ ∂y∂x ∂y ⎥⎦ 2 ∂ Z Обчислюємо H ( x, y) = і визначник матриці H ( x, y) : 1 2 ∂x 2 2 ∂ Z ∂ Z 2 2 2 2 2 2 2 2 2 ∂x ∂x∂y ∂ Z ∂ Z ∂ Z ∂ Z ∂ Z ∂ Z ⎛ ∂ Z ⎞ H 2( x, y) = 2 2 = ⋅ − ⋅ = ⋅ − ⎜ ⎟ . 2 2 2 2 ∂ Z ∂ Z ∂x ∂y ∂y∂x ∂x∂y ∂x ∂y ⎝ ∂x∂y ⎠ 2 ∂y∂x ∂y Якщо H ( x , y ) = D( x , y ) > 0, то в точці (x0; y0) досліджувана 2 0 0 0 0 функція має екстремум. Якщо при цьому H ( x , y ) > 0, то в окресленій точці функція 1 0 0 досягає мінімального значення; якщо H ( x , y ) < 0, то – 1 0 0 максимального значення. Розв’яжемо методом множників Лагранжа таку задачу. 193

для цієї групи методів є градієнтні. Непрямі методи полягають у<br />

зведенні задачі до такої, знаходження оптимального розв’язку якої<br />

вдається спростити. Найпоширенішими методами цього класу є<br />

методи квадратичного програмування.<br />

Оптимізаційні задачі, на змінні яких накладаються обмеження,<br />

розв’язуються методами класичної математики. Оптимізацію з<br />

обмеженнями-рівностями виконують методами зведеного градієнта,<br />

наприклад, методом множників Лагранжа. У задачах оптимізації з<br />

обмеженнями-нерівностями досліджують необхідні та достатні умови<br />

існування екстремуму Куна-Таккера.<br />

Розглянемо метод множників Лагранжа на прикладі такої<br />

задачі нелінійного програмування:<br />

, ,..., ) max(min) ,<br />

x x x f Z = (6.3)<br />

( 1 2 n<br />

q ( x , x ,..., x )<br />

i 1 2 n bi<br />

,..., , ( x x x<br />

i 1 2 n<br />

= , i = 1,<br />

m,<br />

(6.4)<br />

де ( , ,..., )<br />

1 2 n x x x f та q ) – диференційовані.<br />

Ідея методу Лагранжа полягає в заміні окресленої задачі<br />

простішою – знаходження екстремуму складнішої функції, але без<br />

обмежень. Ця функція називається функцією Лагранжа і записується<br />

у вигляді:<br />

L(<br />

x , x ,..., x ; λ , λ ,..., λ ) = f ( x , x ,..., x ) + λ [ q ( x , x ,..., x ) − b ], (6.5)<br />

1<br />

2<br />

n<br />

1<br />

2<br />

m<br />

1<br />

2<br />

де λ – невизначені поки що величини, так звані множники Лагранжа.<br />

i<br />

Необхідною умовою екстремуму функції багатьох змінних є<br />

рівність нулю частинних похідних стосовно всіх змінних функції.<br />

Візьмемо ці частинні похідні і прирівняємо їх до нуля:<br />

⎧ ∂L<br />

⎪<br />

= 0,<br />

j = 1,<br />

n,<br />

⎪∂x<br />

j<br />

⎨<br />

⎪ ∂L<br />

= 0,<br />

i = 1,<br />

m,<br />

⎪⎩<br />

∂λi<br />

або<br />

m<br />

⎧∂f<br />

( x1,<br />

x2,...,<br />

xn<br />

) ∂qi<br />

( x1,<br />

x2,...,<br />

xn<br />

)<br />

⎪<br />

+ ∑ λi<br />

= 0,<br />

j = 1,<br />

n,<br />

⎨ ∂x<br />

j<br />

i= 1 ∂x<br />

j<br />

⎪<br />

⎩qi<br />

( x1,<br />

x2,...,<br />

xn<br />

) − bi<br />

= 0 i = 1,<br />

m.<br />

Отримуємо систему (m+n) рівнянь із (m+n) невідомими,<br />

*<br />

розв’язавши яку, знайдемо X = ( x , x ,..., x ) та λ = ( λ , λ ,..., λ ) –<br />

1 2 n<br />

0 1 2 m<br />

стаціонарні точки. Оскільки їх знайдено з необхідної умови<br />

192<br />

n<br />

m<br />

∑<br />

i=<br />

1<br />

i<br />

i<br />

1<br />

2<br />

n<br />

i

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!