Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Книга_МО.pdf
Скачиваний:
554
Добавлен:
05.06.2015
Размер:
12.36 Mб
Скачать

7.2. Условный экстремум при ограничениях типа неравенств

Пусть требуется исследовать функцию f (x) на экстремум, т.е. определить

точки x U ее локальных минимумов и максимумов

 

f (x ) = min f (x) ,

f (x ) = max f (x) ,

(7.11)

 

U = {x : g j (x) 0,

x U

x U

 

 

 

где

j =1, ..., m}.

Рассмотрим

условный

экстремум

с

ограничениями типа неравенств.

 

 

 

 

 

Стратегия решения задачи

 

 

 

 

 

Находятся точки

x локального

экстремума

с помощью

необходимых

и

достаточных условий минимума и максимума первого и второго порядка при ограничениях типа неравенств. Вычисляются значения f (x ) в найденных точках локального экстремума.

Необходимые условия минимума (максимума) первого порядка. Пусть x

точка локального минимума (максимума) в задаче (7.11). Тогда найдется такое число λ0 0 , и вектор λ = ( λ1 , ..., λm )T , не равные одновременно нулю и такие, что

выполняются следующие условия:

− стационарности обобщенной функции Лагранжа по x

L(x , λ0

, λ )

= 0, i =1, ..., n ;

xi

 

 

 

− допустимости решения

 

 

g j (x ) 0, j =1, ..., m ;

− неотрицательности для условного минимума

(7.12а)

(7.12б)

λj 0, j =1, ..., m

 

(7.12в)

(условие неположительности для условного максимума λj 0,

j =1, ..., m );

− дополняющей нежесткости

 

 

λj g j (x ) = 0, j =1, ..., m .

 

(7.12г)

Если при этом градиенты активных (т.е. g j (x ) = 0 ) в точке

x ограничений

линейно независимы (выполняется условие регулярности), то

λ0

0 . Точки x ,

удовлетворяющие системе (7.12), называются условно-стационарными. Необходимо обратить внимание на следующее.

139

1.В отличие от ограничений типа равенств необходимые условия экстремума первого порядка формулируются отдельно для максимума и минимума. Эти условия впервые были доказаны Куном и Таккером.

2.Если в задаче ограничения записаны в форме g j (x ) 0, то их необходимо

переписать как g j (x ) 0 (7.12).

Обозначим через J a множество индексов ограничений, активных в точке x .

3. При λ0 0 необходимые условия минимума (максимума)

первого порядка

для выпуклых функций f (x), g j (x), j =1, ..., m ( f (x), g j (x),

j =1, ..., m ) будут

одновременно и достаточными условиями глобального минимума (глобального максимума).

4.Из условия дополняющей нежесткости следует, что если ограничение в точке

xпассивное, т.е. g j (x ) < 0 , то λj = 0 , а если активное, т.е. g j (x ) = 0 , то λj 0

(для минимума), и λj 0 (для максимума).

Достаточные условия минимума (максимума) первого порядка. Пусть имеется точка (x , λ ) , удовлетворяющая системе (7.12) при λ0 0 , число активных

ограничений в точке x совпадает с числом

n переменных (при этом условие

регулярности выполняется). Если λj > 0 для всех j J a , то x

− точка условного

локального минимума. Если λj < 0 для всех

j J a , то x

− точка условного

локального максимума в задаче (7.11).

Необходимые условия минимума (максимума) второго порядка. Пусть x

точка локального минимума (максимума) в задаче (7.11) и имеется решение (x , λ )

системы (7.12). Тогда второй дифференциал классической функции Лагранжа,

вычисленный в точке (x , λ ) , неотрицателен (неположителен)

d 2 L(x , λ ) 0

 

 

( d 2 L(x , λ ) 0 )

(7.13)

для всех d x таких, что

 

 

 

 

 

 

d g j (x ) = g j (x

 

) d xi = 0, j J a , λj > 0 (λj < 0)

;

n

 

 

 

 

 

 

 

 

 

 

 

 

 

=

x

i

 

 

 

 

i 1

 

 

 

 

 

d g j (x ) 0, j J a , λj = 0 .

140

Достаточные условия минимума (максимума) второго порядка. Пусть имеется точка (x , λ ) , удовлетворяющая системе (7.12) при λ0 0 . Если в этой точке d 2 L(x , λ ) > 0 ( d 2 L(x , λ ) < 0 ) для всех ненулевых d x таких, что

d g j (x ) = g j (x

 

) d xi = 0, j J a , λj > 0 (λj < 0) ;

n

 

 

 

 

 

 

 

 

 

 

 

=

x

i

 

 

 

i 1

 

 

 

 

d g j (x ) 0, j J a , λj = 0 ,

то x − точка локального минимума (максимума) в задаче (7.11).

Алгоритм решения задачи.

Шаг 1. Составить обобщенную функцию Лагранжа

λj

λ0

m

L(x, λ0 , λ) = λ0 f (x) + λj g j (x) .

j=1

Шаг 2. Записать необходимые условия минимума (максимума) первого порядка

а)

L(x , λ0 , λ )

= 0, i =1, ..., n ;

б) g j (x

 

) 0, j =1, ..., m ;

xi

 

 

 

 

 

 

в) λj 0, j =1, ..., m (для минимума), λj 0, j =1, ..., m (для максимума);

г) λj g j (x ) = 0, j =1, ..., m .

Шаг 3. Решить систему для двух случаев

1)λ0 = 0 ;

2)λ0 0 (при этом поделить условия, записанные на шаге 2, на λ0 и заменить

на λj ).

Врезультате найти условно-стационарные точки x , выделив из них полученные при λ0 0 (они могут быть регулярными точками экстремума).

Вкаждом из двух случаев следует начинать с рассмотрения 2m вариантов удовлетворения условия "г" дополняющей нежесткости.

Шаг 4. Для выделенных на шаге 3 точек проверить достаточные условия экстремума первого или второго порядка.

Для проверки достаточных условий первого порядка следует: а) определить число l активных в точке x ограничений;

б) если l = n и λj > 0 для всех j J a , то в точке x − локальный минимум;

141

если l = n и λj < 0 для всех j J a , то в точке x − локальный максимум. Если l < n

или соответствующие множители Лагранжа не удовлетворяют достаточным условиям первого порядка, то проверить достаточные условия второго порядка.

Для проверки достаточных условий второго порядка следует:

а) записать выражение для второго дифференциала классической функции Лагранжа в точке (x , λ )

n n

2

L(x

 

, λ

 

) dxi d x j ;

d 2 L(x , λ ) = ∑∑

 

 

 

i=1 j=1

 

xi x j

 

 

б) записать условия, накладываемые на первые дифференциалы активных ограничений

d g j (x ) = g j (x

 

) d xi

= 0, j J a ;

λj > 0 (λj < 0) ;

(7.14)

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i=1

xi

 

 

 

 

 

 

 

 

 

 

 

d g j (x ) = g j (x

 

) d xi 0,

j J a ; λj = 0 ;

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

x

i

 

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

в) исследовать знак второго дифференциала функции Лагранжа для ненулевых

d x , удовлетворяющих системе

(7.14). Если d 2 L(x , λ ) > 0 , то

в

точке x

условный локальный минимум.

Если d 2 L(x , λ ) < 0 , то в точке

x

− условный

локальный максимум.

 

 

 

Если достаточные условия первого и второго порядка не выполняются, следует проверить выполнение необходимых условий второго порядка, следуя аналогичной процедуре. Если они выполняются, то требуется дополнительное исследование, а

если нет, то в точке x

нет условного экстремума.

 

 

 

 

 

Шаг 5. Вычислить значения функции в точках условного экстремума.

 

 

Условия экстремума в исходной задаче (7.11) приведены в табл. 7.2. и 7.3.

 

Необходимые и достаточные условия первого порядка в задаче поиска условного

 

 

экстремума при ограничениях типа неравенств

Таблица 7.2.

 

 

 

 

 

 

 

 

 

 

Необходимые условия первого порядка

Достаточные условия первого порядка (λ0 0)

 

 

 

 

 

 

 

 

x L(x , λ0 , λ );

 

g j (x ),

λ0 0,

Число

λj ,

Тип условно-стационарной

 

 

λj g j (x ),

 

j =1, m

λj , j =1, m

l

j J a

 

точки x

 

 

j =1, ..., m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

 

0

0

n

> 0

условный локальный минимум

 

 

 

 

 

 

 

 

 

2

 

0

 

0

0

n

<0

условный локальный максимум

 

 

 

 

 

 

 

 

 

 

142

Необходимые и достаточные условия второго порядка в задаче поиска условного

 

 

экстремума при ограничениях типа неравенств

Таблица 7.3.

 

 

 

 

 

 

d 2 L(x , λ )

dg j (x ),

dg j (x ),

dg j (x ),

Тип условно-стационарной точки x

 

 

j J a ,

j J a ,

j J a ,

 

 

 

 

λj > 0

λj < 0

λj = 0

 

 

 

 

 

 

 

 

1

> 0

0, dx 0

0

условный локальный минимум

 

 

 

 

 

 

2

<0

0, dx 0

0

условный локальный максимум

 

 

 

 

 

 

3

0

0

0

может быть локальный минимум,

 

 

 

 

 

требуется дополнительное исследование

 

 

 

 

 

 

4

0

0

0

может быть локальный максимум,

 

 

 

 

 

требуется дополнительное исследование

 

 

 

 

 

 

5

= 0

0

0

требуется дополнительное исследование

 

 

 

 

 

 

6

= 0

0

0

требуется дополнительное исследование

 

 

 

 

 

 

7

> 0, < 0

0

0

нет экстремума

 

 

 

 

 

 

8

> 0, < 0

0

0

нет экстремума

 

 

 

 

 

 

 

Пример 7.8. Найти условный экстремум

f (x) = x12 + x22 extr , g1 (x) = x1 + x2 2 0 .

□ 1. Составим обобщенную функцию Лагранжа

L(x, λ0 , λ1 ) = λ0 (x12 + x22 ) + λ1 (x1 + x2 2) .

2. Выпишем необходимые условия экстремума первого порядка

а) L(x, λ0 , λ1 ) = 2λ

x + λ = 0,

L(x, λ0 , λ1 ) = 2λ

0

x

2

+ λ = 0

;

x1

0

1

1

x2

 

1

 

 

 

 

 

 

 

 

 

б) x1 + x2 2 0 ;

 

 

 

 

 

 

 

 

 

в) λ1 0 (для минимума),

λ1 0 (для максимума);

 

 

 

 

 

г) λ1 (x1 + x2 2) = 0 .

 

 

 

 

 

 

 

 

3.Решим систему для двух случаев.

Впервом случае λ0 = 0 . Тогда из условия "а" следует, что λ1 = 0 . Это

противоречит требованию существования ненулевого вектора (λ0 , λ1 )T .

143

Во втором

случае

λ

0

0 .

Поделим систему на λ

0

и

заменим

 

λ1

на

λ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Обобщенная функция Лагранжа при этом заменяется классической

 

 

 

 

а) L(x, λ1 ) = 2x + λ = 0,

L(x, λ1 ) = 2x

2

+ λ = 0 ;

 

 

 

 

 

 

 

 

x1

 

1

1

 

 

x2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) x1 + x2 2 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

в) λ1 0 (для минимума),

λ1 0 (для максимума);

 

 

 

 

 

 

г) λ1 (x1 + x2 2) = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

Из условия "г" дополняющей нежесткости следует

 

 

 

 

 

 

 

1) λ1

= 0 (фактически решается задача поиска безусловного экстремума).

 

Тогда x = x = 0,

λ = 0 и условие "б" выполняется. Выполняются необходимые

1

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

условия и для минимума, и для максимума.

 

 

 

 

 

 

 

 

 

2) λ1

0 . Тогда из системы x1 + x2 2 = 0 ,

 

2x1 + λ1 = 0 ,

2x2 + λ1

= 0

 

получаем x = x

 

=1, λ = −2 . Так как λ < 0 , то необходимое условие минимума не

 

1

2

 

1

 

 

1

 

 

 

 

 

 

 

 

 

выполняется (в точке (1, 1)T нет минимума), но выполняется необходимое условие максимума. Таким образом, имеем две условно-стационарные точки.

4.Проверим выполнение достаточных условий экстремума.

Вточке x = (0, 0)T ограничение не является активным, так как g1 (x) = −2 < 0 ,

поэтому достаточные условия первого порядка не удовлетворяются. Проверим

условия второго порядка. Так как d 2 L(x , λ ) = 2dx2

+ 2dx2

> 0 при dx 0 , то в точке

 

1

2

 

 

x = (0, 0)T

регулярный локальный условный минимум,

 

совпадающий в данной

задаче с безусловным.

 

 

 

В точке

x = (1, 1)T ограничение является активным,

но l =1 < n = 2 , поэтому

достаточное условие первого порядка не выполняется. Проверим достаточное условие второго порядка.

 

 

 

d 2 L(x , λ ) = 2dx2 + 2dx2

,

 

 

 

 

 

 

 

1

2

 

 

 

dg

(x ) = dx + dx

2

= 0

откуда dx = −dx .

 

1

 

1

 

 

 

1

2

Следовательно, d 2 L(x , λ ) = 4dx22

> 0

при

dx2

0 .

Так как в этой точке

λ = −2 < 0

, то достаточное

условие максимума

не

выполняется. Проверим

1

 

 

 

 

 

 

 

 

 

достаточное условие максимума второго порядка. Так как d 2 L(x , λ ) = 4dx22 0 при

144

любых dx2 , то необходимое условие максимума не выполняется, поэтому в точке x = (1, 1)T максимума нет.

5. Вычислим значение функции в точке условного минимума f (x ) = 0

(рис.7.7). ■

x2

2

g(x) = x1 + x2 - 2 = 0

1

2 x1

x* 1

Рис.7.7. Иллюстрация к решению задачи из примера 7.8 Пример 7.9. Найти условный экстремум в задаче

f (x) = x1 + x2 extr ,

g1 (x) = x12 + x22 1 0 .

□ 1. Составим обобщенную функцию Лагранжа

L(x, λ0 , λ1 ) = λ0 (x1 + x2 ) + λ1 (x12 + x22 1) .

2. Выпишем необходимые условия первого порядка

а) L(x, λ0 , λ1 ) = λ

0

+ 2λ x = 0

,

L(x, λ0 , λ1 ) = λ

0

+ 2λ x

2

= 0

;

 

x1

1

1

 

x2

1

 

 

 

 

 

 

 

 

 

 

 

 

б) x2

+ x2

1 0 ;

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

в) λ1 0 (для минимума),

 

λ1 0 (для максимума);

 

 

 

г) λ (x2 + x2 1) = 0 .

 

 

 

 

 

 

 

 

 

1

1

2

 

 

 

 

 

 

 

 

 

 

3. Решим задачу для двух случаев.

145

В первом случае λ0

= 0 .

Тогда, согласно необходимым условиям экстремума

первого порядка требуется, чтобы λ1 0 . При этом x1

 

= x2 = 0 и не удовлетворяется

условие "г" дополняющей нежесткости.

 

 

 

 

 

 

Во втором случае λ

0

0 . Поделим уравнения на λ

0

и заменим

λ1

на λ

 

 

 

 

 

 

 

 

 

λ0

1

 

 

 

 

 

 

 

 

 

 

 

а) L(x, λ1 ) =1 + 2λ x = 0

L(x, λ1 ) =1 + 2λ x

2

= 0 ;

 

 

 

x1

1

1

 

1

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

б) x2

+ x2

1 0 ;

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

в) λ1 0 (для минимума),

λ1 0 (для максимума);

 

 

г) λ (x2 + x2 1) = 0 .

 

 

 

 

 

 

 

 

 

1

1

2

 

 

 

 

 

 

 

 

 

Из условия "г" дополняющей нежесткости следует два варианта

1)λ1 = 0 . Тогда условие "а" не выполняется.

2)λ1 0 . Тогда x12 + x22 1 = 0 и система имеет два решения (рис. 7.8)

точка A : x

= x = −

2 ,

λ

=

2

(в точке A может быть минимум);

 

 

1

 

2

 

2

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

точка B : x

= x =

 

2 ,

λ = −

2

(в точке B может быть максимум).

 

 

1

 

2

2

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4. Проверим достаточные условия экстремума.

В точках A и B ограничения

являются активными,

но

l =1 < n = 2 .

Поэтому

 

условия

первого

порядка

не

выполняются. Проверим условия второго порядка:

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2 L(x , λ ) = 2λ dx2 + 2λ dx2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

1

1

2

 

 

 

 

 

 

 

 

 

 

 

d g

(x

) = 2x dx + 2x dx

2

= 0 .

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

1

 

2

 

 

 

 

 

 

 

В точках

A и B выполняется dx

= −dx

2

. Так как d 2 L( A) = 4λ dx 2

= 2 2dx 2

> 0

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

2

2

 

при dx2 0 ,

то в точке

A − регулярный локальный условный минимум. Так как

d 2 L(B) = 4λ dx 2 = −2

 

2dx 2 < 0

при

dx

2

0 , то в точке

B − регулярный локальный

1

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

условный максимум. С другой

стороны,

 

функции

f (x)

и

f (x) = −x1 x2

выпуклые и

ограничение

выпуклое,

поэтому в

точках

A

и

B

достигается

глобальный экстремум. Достаточные условия первого и второго порядка проверялись для демонстрации методики.

5. Вычислим значение целевой функции в точках условного экстремума: f ( A) = − 2, f (B) = 2 . ■

146

Рис.7.8. Иллюстрация к решению задачи из примера 7.9

Пример 7.10. Найти условный экстремум

f (x) = (x1 2)2 + (x2 3)2 extr, g1 (x) =x21 +x22 52 0.

□ 1. Составим обобщенную функцию Лагранжа

 

 

L(x, λ

0

, λ ) = λ

0

((x

2)2 + (x

2

3)2 ) + λ (x2

+ x2

52) .

 

 

 

 

 

 

 

1

1

 

1

1

 

 

2

 

 

 

 

 

2. Выпишем необходимые условия первого порядка

 

 

 

 

 

 

 

 

 

а) L(x, λ0 , λ1 ) = 2λ

0

(x 2) + 2λ x = 0 , L(x, λ0 , λ1 ) = 2λ

0

(x

2

3) + 2λ x

2

= 0

;

 

x1

 

1

 

1

1

 

x2

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) x2

+ x2 52 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в) λ1 0 (для минимума),

 

λ1 0 (для максимума);

 

 

 

 

 

 

 

г) λ (x2

+ x2 52) = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3.Решим задачу для двух случаев.

Впервом случае λ0 = 0 . Тогда, согласно необходимым условиям экстремума

первого порядка требуется, чтобы λ1 0 . При этом x1 = x2 = 0 и не выполняется условие "г" дополняющей нежесткости.

147

Во втором случае λ

0

0 . Поделим уравнения на λ

0

и заменим

λ1

на λ

 

 

 

 

 

 

 

 

 

 

 

 

 

λ0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а) L(x, λ1 ) = 2(x 2) + 2λ x = 0

L(x, λ1 ) = 2(x

2

3) + 2λ x

2

= 0

;

 

 

x1

 

1

 

 

1

1

x2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) x2

+ x2

52 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в) λ1 0 (для минимума),

λ1

0 (для максимума);

 

 

 

 

г)

λ (x2 + x2 52) = 0 .

 

 

 

 

 

 

 

 

 

 

 

 

1

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из условия "г" следуют два варианта

 

 

 

 

 

 

 

 

1)

λ1

= 0 .

Тогда

x1 = 2, x2

= 3 и

выполняются необходимые условия и для

минимума,

и

для

 

максимума.

Имеем условно-стационарную

точку A :

x = 2, x

= 3, λ = 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

2) λ1 0 . Тогда x12 + x22 52 = 0 и система имеет решение:

точка B : x

= 4, x = 6,

λ = −

1

; точка C : x = −4, x = −6,

λ = −

3

.

 

 

1

2

1

2

1

2

1

2

 

 

 

 

 

 

 

 

Так как

λ1 < 0 в

обеих

точках, то в них

минимума нет, но может быть

максимум.

 

 

 

 

 

 

 

 

4. Проверим достаточные условия экстремума. В обеих условно-стационарных точках ограничение превращается в равенство, т.е. активно. Так как число активных ограничений l =1 < 2 = n , то условия первого порядка не выполняются.

Так как функция f (x) = −(x1 2)2 (x2 3)2 не является выпуклой, то необходимые условия не являются достаточными. Поэтому проверим условия второго порядка

 

 

 

d 2 L(x , λ ) = (2 + 2λ )dx2

+ (2 + 2λ )dx2 .

 

 

 

1

1

1

1

2

В точке

A

ограничение не

является

активным. Так как λ = 0 , то

 

 

 

 

 

 

 

1

d 2 L(x , λ ) = 2dx2

+ 2dx2

> 0 при dx 0

. Поэтому в точке A − условный локальный

1

1

2

 

 

 

 

 

минимум. Так как целевая функция выпуклая и множество допустимых значений выпукло, то можно заключить, что в данном случае необходимое условие минимума является достаточным. В точке A − глобальный минимум.

В точках B и C ограничение активно. Поэтому d g

(x ) = 2x dx

+ 2x dx

2

= 0 . В

 

 

 

 

 

 

 

 

1

 

 

1

1

 

2

 

точках

B и

C выполняется dx

= −

3

dx

 

. Так как

d 2 L(B) =

(

3

)dx

 

2

+ dx2

> 0 при

 

2

 

2

 

 

 

1

2

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

148

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]