УМК (1013374), страница 26
Текст из файла (страница 26)
Решением системы являются две условно-стационарные точки:A: x1∗ = 1, x 2∗ = 1, λ∗1 = −1;2B: x1∗ = −1, x 2∗ = −1, λ∗1 =1.24. Проверим выполнение достаточных условий экстремума:∂ 2 L (x , λ 1 ) ∂ 2 L (x , λ 1 )а) d 2 L( x ∗ , λ 1∗ ) = 2λ 1∗ dx12 + 2λ 1∗ dx 22 , так как== 2λ∗1 ,22∂x1∂x 2∂ 2 L (x , λ 1 )∂x1∂x 2∂ 2 L (x , λ 1 )== 0;∂x 2 ∂x1∂ g1 ( x)∂ g1 ( x)= 2x2 ;= 2 x1 ,∂ x1∂ x2в) исследуем точку A.
Получаем dg1 ( A ) = 2dx1 + 2dx 2 = 0 , откуда dx1 = − dx 2 .б) dg1 ( x ∗ ) = 2 x1∗ dx1 + 2 x 2∗ dx 2 = 0 , так какС учетом полученного соотношения d 2 L ( A ) = − dx12 − dx 22 = −2dx 22 < 0 при dx 2 ≠ 0 .Поэтому в точке x ∗ = (1,1) – регулярный условный локальный максимум.T154Исследуем точку B. Получаем dg1 (B ) = −2dx1 − 2dx 2 = 0 , откуда dx1 = − dx 2 .С учетом полученного соотношения d 2 L (B ) = dx12 + dx 22 = 2dx 22 > 0 при dx 2 ≠ 0 .Поэтому в точке x ∗ = ( −1, − 1) – регулярный условный локальный минимум.T5.
Подсчитаем значения функции в точках экстремума: f ( A ) = 2, f (B ) = −2 .Графическое решение задачи изображено на рис. 2. x222A1g 1 ( x ) = x1 2 + x 2 2 − 2 = 0−121−22f (x ) = 2−1Bx1f (x ) = 0−2f (x ) = −2Рис. 2Б. ОГРАНИЧЕНИЯ ТИПА НЕРАВЕНСТВПостановка задачиДаны дважды непрерывно дифференцируемые целевая функция f ( x) = f ( x1,…, xn )и функции ограничений g j ( x) = g j ( x1 ,… , x n ) ≤ 0 , j = 1,… , m , определяющие множестводопустимых решений X .Требуется исследовать функцию f ( x ) на экстремум, т.е. определить точки x ∗ ∈ Xее локальных минимумов и максимумов на множестве X :f ( x ∗ ) = min f ( x) ;{}x∈Xf ( x ∗ ) = max f ( x) ,x∈Xгде X = x g j ( x) ≤ 0, j = 1,… , m .Алгоритм решения задачиШаг 1. Составить обобщенную функцию Лагранжа:mL ( x, λ 0 , λ ) = λ 0 f ( x ) + ∑ λ j g j ( x ) .j =1155Шаг 2.
Записать необходимые условия минимума (максимума) первого порядка:∂ L( x ∗ , λ ∗0 , λ ∗ )= 0,i = 1,… , n ;а)∂ xiб) g j ( x ∗ ) ≤ 0 , j = 1, … , m ;в) λ∗j ≥ 0 , j = 1, … , m (для минимума), λ∗j ≤ 0 , j = 1, … , m (для максимума);г) λ ∗j g j ( x ∗ ) = 0 , j = 1, … , m .Шаг 3. Решить систему для двух случаев.Первый случай: λ∗0 = 0 .Второй случай: λ∗0 ≠ 0 (при этом поделить условия, записанные на шаге 2, на λ∗0 изаменитьλ∗jλ∗0на λ∗j ).В результате решения найти условно-стационарные точки x ∗ , выделив из нихполученные при λ∗0 ≠ 0 (они могут быть регулярными точками экстремума).
В каждомиз двух случаев следует начинать с рассмотрения 2 m вариантов удовлетворения условия«г» дополняющей нежесткости.Шаг 4. Для выделенных на шаге 3 точек проверить выполнение достаточныхусловий экстремума первого или второго порядка.Для проверки выполнения достаточных условий первого порядка следует:а) определить число l активных в точке x ∗ ограничений;б) если l = n и λ∗j > 0 для всех j ∈ J a , то в точке x ∗ – условный локальныйминимум. Если l = n и λ∗j < 0 для всех j ∈ J a , то в точке x ∗ – условныйлокальный максимум. Если l < n или соответствующие множители Лагранжане удовлетворяют достаточным условиям первого порядка, проверитьдостаточные условия второго порядка.Для проверки выполнения достаточных условий второго порядка следует:а) записать выражение для второго дифференциала классической функцииЛагранжа в точке ( x ∗ , λ ∗ ) :∂ 2 L( x ∗ , λ ∗ )d L( x , λ ) = ∑ ∑dx i dx j ;∂x i ∂x ji =1 j =1б) записать условия, накладываемые на первые дифференциалы активныхограничений:n ∂ g (x∗ )j∗dg j ( x ) = ∑dx i = 0 , j ∈ J a ; λ∗j > 0 ( λ∗j < 0 );∂ xii =1∗2n∂ g j (x∗ )i =1∂ xidg j ( x ) = ∑∗∗nndx i ≤ 0 , j ∈ J a , λ∗j = 0 ;в) исследовать знак второго дифференциала функции Лагранжа при ненулевыхdx , удовлетворяющих системе, составленной в п.б.
Если d 2 L( x ∗ , λ ∗ ) > 0 , то156в точке x ∗ – условный локальный минимум. Если d 2 L( x ∗ , λ ∗ ) < 0 , то в точкеx ∗ – условный локальный максимум.Если достаточные условия первого и второго порядка не выполняются, следуетпроверить выполнение необходимых условий второго порядка (см. утверждение 3.6 –лекция 2), следуя аналогичной процедуре. Если они выполняются, то требуетсядополнительное исследование, а если нет, то в точке x ∗ нет условного экстремума.Шаг 5. Вычислить значения функции в точках условного экстремума.Пример 2.
Найти условный минимум в задачеf ( x) = x12 + ( x 2 − 2) 2 → min ,g1 ( x) = x12 + x 22 − 1 ≤ 0 ,g 2 ( x) = − x1 ≤ 0 ,g 3 ( x) = − x 2 ≤ 0 . 1. Составим обобщенную функцию Лагранжа:L ( x, λ 0 , λ ) = λ 0 ⎡ x12 + ( x 2 − 2) 2 ⎤ + λ 1 ( x12 + x 22 − 1) + λ 2 (− x1 ) + λ 3 (− x 2 ) .⎣⎦2. Выпишем необходимые условия минимума первого порядка:∂ L (x , λ 0 , λ )а)= 2λ 0 x1 + 2λ1 x1 − λ 2 = 0 ,∂ x1∂ L (x , λ 0 , λ )= 2λ 0 (x 2 − 2) + 2λ1 x 2 − λ 3 = 0 ;∂ x2б) x12 + x 22 − 1 ≤ 0 , − x1 ≤ 0 , − x 2 ≤ 0 ;в) λ1 ≥ 0 , λ 2 ≥ 0 , λ 3 ≥ 0 ;г) λ 1 ( x12 + x 22 − 1) = 0 , λ 2 (− x1 ) = 0 , λ 3 (− x 2 ) = 0 .3. Решим систему для двух случаев.Первый случай: λ 0 = 0 , тогда условия “а” запишутся в виде2λ1x1 − λ 2 = 0 ,2λ1x 2 − λ 3 = 0 .Рассмотрим восемь вариантов выполнения условий «г» дополняющейнежесткости:1) λ1 = 0 , λ 2 = 0 , λ 3 = 0 , при этом не удовлетворяется требование утверждения3.4;2) λ1 ≠ 0 , λ 2 = 0 , λ 3 = 0 , тогда x1 = x 2 = 0 из условия «а», но первое условиедополняющей нежесткости не удовлетворяется;3) λ1 = 0 , λ 2 ≠ 0 , λ 3 = 0 , тогда из первого уравнения в условии «а» имеемλ 2 = 0 , т.е.
имеется противоречие;4) λ1 = 0 , λ 2 = 0 , λ 3 ≠ 0 , тогда из второго уравнения в условии «а» имеем λ 3 = 0 ,т.е. также имеется противоречие;5) λ1 ≠ 0 , λ 2 ≠ 0 , λ 3 = 0 , тогда x1 = 0 и из первого уравнения в условии «а»имеем λ 2 = 0 , т.е. имеется противоречие;1576) λ1 ≠ 0 , λ 2 = 0 , λ 3 ≠ 0 , тогда x 2 = 0 и из второго уравнения в условии «а»имеем λ 3 = 0 , т.е. также имеется противоречие;7) λ1 = 0 , λ 2 ≠ 0 , λ 3 ≠ 0 , тогда не выполняются оба уравнения в условии “а”;8) λ1 ≠ 0 , λ 2 ≠ 0 , λ 3 ≠ 0 , тогда уравнения x1 = x 2 = 0 , x12 + x 22 − 1 = 0 ,следующие из условия “г”, вместе не выполняются.Условно-стационарных точек пока не найдено.Второй случай: λ 0 ≠ 0 . Поделив уравнения приведенной в п.2 системы на λ 0 иλλλзаменив 1 на λ1 , 2 на λ 2 , 3 на λ 3 , получим:λ0λ0λ0∂ L ( x, λ )∂ L (x , λ )= 2 ( x 2 − 2) + 2 λ 1 x 2 − λ 3 = 0 ;а)= 2 x1 + 2λ1 x1 − λ 2 = 0 ,∂ x2∂ x1б) x12 + x 22 − 1 ≤ 0 , − x1 ≤ 0 , − x 2 ≤ 0 ;в) λ1 ≥ 0 , λ 2 ≥ 0 , λ 3 ≥ 0 ;г) λ 1 ( x12 + x 22 − 1) = 0 , λ 2 (− x1 ) = 0 , λ 3 (− x 2 ) = 0 .Рассмотрим восемь вариантов выполнения условий дополняющей нежесткости:1) λ1 = 0 , λ 2 = 0 , λ 3 = 0 , тогда x1 = 0, x 2 = 2 и не выполняется первоеограничение в условии «б»;2) λ1 ≠ 0 , λ 2 = 0 , λ 3 = 0 , тогдаx12 + x 22 − 1 = 0 ,2 x1 (1 + λ1 ) = 0 ,2 (x 2 − 2) + 2λ1 x 2 = 0 .Если λ1 = −1 , то третье уравнение не удовлетворяется.
Если x1 = 0 , то x 2 = ±1 .Ограничениям в условии «б» удовлетворяет x 2 = 1 , при этом λ1 = 1 . Получили условно-стационарную точку А: x1∗ = 0, x 2∗ = 1, λ∗1 = 1, λ∗2 = 0, λ∗3 = 0 ;3) λ1 = 0 , λ 2 ≠ 0 , λ 3 = 0 , тогдаx1 = 0 ,2 x1 − λ 2 = 0 ,2 (x 2 − 2 ) = 0 .Получаем λ 2 = 2 x1 = 0 , что противоречит условию λ 2 ≠ 0 ;4) λ1 = 0 , λ 2 = 0 , λ 3 ≠ 0 , тогдаx2 = 0 ,2 x1 = 0 ,2( x 2 − 2) − λ 3 = 0 .Получаем λ 3 = − 4 < 0 , что противоречит условию «в»;5) λ1 ≠ 0 , λ 2 ≠ 0 , λ 3 = 0 , тогдаx12 + x 22 − 1 = 0 ,x1 = 0 ,2 x1 + 2λ1x1 − λ 2 = 0 ,2 (x 2 − 2) + 2λ1 x 2 = 0 .Из третьего соотношения следует, что λ 2 = 0 , т.е. имеется противоречие;1586) λ1 ≠ 0 , λ 2 = 0 , λ 3 ≠ 0 , тогдаx12 + x 22 − 1 = 0 ,x2 = 0 ,2 x1 + 2λ1x1 = 0 ,2 ( x 2 − 2 ) + 2λ 1 x 2 − λ 3 = 0 .Из последнего соотношения следует, что λ 3 = − 4 < 0 .
Это противоречит условию «в»;7) λ1 = 0 , λ 2 ≠ 0 , λ 3 ≠ 0 , тогдаx1 = x 2 = 0 ,2 x1 − λ 2 = 0 ,2( x 2 − 2) − λ 3 = 0 .Из второго соотношения следует, что λ 2 = 0 . Это противоречит условию λ 2 ≠ 0 ;8) λ1 ≠ 0 , λ 2 ≠ 0 , λ 3 ≠ 0 , тогда x12 + x 22 − 1 = 0 . Из условия «г» следует, чтоx1 = 0, x 2 = 0 . Эта система несовместна.x2f ( x) = 121g1 ( x) = 0Ag 3 ( x) = − x 2 = 01g 2 ( x) = − x1 = 0x1Рис.
34. Проверим выполнение достаточных условий минимума. В точке А имеются дваактивных ограничения, т.е. l = 2 = n = 2 (рис. 3). Так как λ∗1 = 1 > 0, λ∗2 = 0 , тодостаточные условия минимума первого порядка не выполняются ввиду того, чтотребуется строгая положительность соответствующих множителей Лагранжа. Проверимусловия второго порядка: d 2 L ( A ) = (2 + 2λ 1∗ ) dx12 + (2 + 2λ 1∗ ) dx 22 .
Поскольку в точке дваактивных ограничения и для одного из них λ∗1 > 0 , а для другого λ∗2 = 0 , то применимусловия:dg 1 ( A ) = 2 x1∗ dx1 + 2 x 2∗ dx 2 = 2dx 2 = 0, λ∗1 > 0 ;dg 2 ( A ) = − dx1 ≤ 0,159λ∗2 = 0 .В результате d 2 L( A ) = 4dx12 > 0 при dx1 ≥ 0 и dx1 ≠ 0 . Поэтому в точке А – локальныйусловный минимум. С другой стороны, целевая функция и множество допустимыхрешений выпуклые. Поэтому в точке А достигается глобальный минимум.5. Вычислим значение функции в точке глобального минимума: f ( A ) = 1 . В.