- •ISBN
- •Введение в методы оптимизации
- •1.1. Функция спроса и ее эластичность
- •1.2. Функция предложения и рыночное равновесие
- •1.3. Предельные величины в экономике и оптимизация прибыли
- •1.4. Основные виды функций нескольких переменных в экономических задачах
- •1.5. Предельная полезность товара и предельная норма замещения
- •1.6. Критерий оптимального набора товаров и оптимального производственного плана.
- •ТРАНСПОРТНАЯ ЗАДАЧА
- •§ 2.1. Постановка задачи
- •§ 2.2. Построение начального опорного плана транспортной задачи
- •§ 2.3. Решение транспортной задачи методом потенциалов
- •§ 2.4. Открытая модель транспортной задачи
- •§ 2.5. Определение оптимального плана транспортных задач
- •с дополнительными ограничениями
- •Метод искусственного базиса
- •3.1. Симплекс-метод решения задач линейного программирования
- •2.2. Метод искусственного базиса
- •4.1. Постановка задачи. Графический метод решения
- •4.2. Двойственный симплекс-метод
- •4.3. Метод Гомори
- •Задачи многокритериальной
- •оптимизации
- •5.1. Общая постановка задачи многокритериальной оптимизации. Парето-эффективное множество.
- •5.2. Методы решения многокритериальной задачи оптимизации
- •Элементы теории игр
- •6.1. Платежная матрица. Нижняя и верхняя цена игры.
- •6.2. Решение игры в смешанных стратегиях.
- •6.3. Приведение матричной игры к задаче линейного программирования.
- •6.4. Игра с природой.
- •7.1. Выпуклые функции
- •7.2. Теорема Куна-Таккера.
- •8.1. Задача динамического программирования
- •8.2. Задача о распределении средств между предприятиями
- •Рекомендуемая литература
Глава 7 Выпуклые функции и теорема Куна-Таккера
В этой главе изучаются экстремумы выпуклых функций на вы-
пуклых множествах в пространстве Rn . Причем если выпуклое множество задано системой ограничений типа равенств, то решение задачи сводится к поиску критических точек функции Лагранжа. Если же выпуклое множество определено системой неравенств, то оптимальное значение целевой функции находится с помощью решения системы Куна-Таккера, которая и является основной темой данной главы.
7.1. Выпуклые функции
Напомним, что множество M Rn называется выпуклым, если для любых двух точек A M и B M отрезок AB целиком принад-
лежит M , т.е. для любых точек A, B M ,t 0,1 |
точка |
[ ] |
(1−t)A +tB M . Если же найдется пара точек множества M , для ко-
торых это условие не выполняется, то множество выпуклым не является.
Пусть теперь точка X = (x1, x2 , , xn ) принадлежит выпуклому
множеству M Rn , |
а функция |
|
f (X )= f (x1, x2 , , xn ) определена на |
|||||||||||||||||||||||||||||
M . |
|
|
|
|
|
|
|
|
|
|
|
|
|
называется выпуклой на множе- |
||||||||||||||||||
Определение 7.1. Функция |
f |
|
||||||||||||||||||||||||||||||
стве M Rn , |
если для любых точек |
|
|
A, B M ,t 0,1 имеет место |
||||||||||||||||||||||||||||
|
(( |
|
) |
|
|
|
) |
( |
|
|
) |
|
|
( |
|
|
) |
|
|
|
|
( |
|
) |
|
[ |
|
] |
||||
неравенство f |
|
A +tB |
|
−t |
|
f |
A |
+tf |
B |
. |
|
|
|
|||||||||||||||||||
1−t |
|
|
|
≤ |
1 |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Если для любых точек A, B M ,t ( |
0,1) имеет место строгое |
|||||||||||||||||||||||||||||||
неравенство f |
(( |
|
) |
A |
+tB |
) |
|
( |
|
|
) |
f |
|
( |
A |
) |
+tf |
( |
B |
) |
, то функция f на- |
|||||||||||
1−t |
|
|
< 1−t |
|
|
|
|
|
|
|
|
|||||||||||||||||||||
зывается строго выпуклой на множестве M Rn . |
|
|||||||||||||||||||||||||||||||
Если же для любых A, B M ,t |
|
0,1 |
|
имеет место обратное не- |
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
[ |
|
|
|
] |
|
|
|
|
|
|
|
|
|
||
равенство |
|
(( |
|
|
) |
|
|
|
|
) |
|
( |
|
|
|
|
) |
|
|
( |
|
|
) |
|
|
( |
|
) |
|
|||
|
f |
|
|
A +tB |
|
|
|
|
|
f |
|
A |
+tf |
B |
, |
|||||||||||||||||
|
1−t |
|
|
≥ 1−t |
|
|
|
|
|
|
|
|
|
92
то функция f (X ) называется вогнутой (а в случае строгого неравен-
ства – строго вогнутой).
Предположим теперь, что функция y = f (X ) имеет непрерывные частные производные 2-ого порядка на выпуклом множестве
M Rn . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Определение 7.2. Матрицей Гессе функции y = f (X ) |
называ- |
|||||||||||||||||||
ется матрица |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂2 f |
|
∂2 f |
|
|
∂2 f |
|
|
|
|
|
|
|
|
|
||||
|
|
∂x2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
∂x |
∂x |
∂x ∂x |
|
|
|
|
|
|
|
|
||||||||
|
|
|
1 |
1 |
2 |
|
|
1 |
n |
|
|
|
|
|
|
|
|
|||
|
∂2 f |
|
∂2 f |
|
|
∂2 f |
|
|
|
2 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
H ( f )= |
|
∂x |
∂x |
|
∂x2 |
|
∂x |
∂x |
|
= |
∂ |
|
|
f |
|
(X ) |
, |
|||
|
|
|
|
∂x |
|
∂x |
|
|||||||||||||
|
2 |
1 |
|
|
2 |
|
|
2 |
n |
|
|
j |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
∂2 f |
|
∂2 f |
|
|
∂2 f |
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
∂x2 |
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
∂x |
|
∂x |
∂x |
|
|
|
|
|
|
|
|
|
|
||||
|
|
n |
1 |
|
n |
2 |
|
|
|
n |
|
|
|
|
|
|
|
|
|
составленная из вторых частных производных функции y = f (X ).
Определитель матрицы Гессе называется гессианом.
Замечание 7.1. Отметим, что в силу условия непрерывности частных производных 2-ого порядка матрица H является симметричной. Поэтому элементы матрицы Гессе, зависящие от точки X M , можно рассматривать как коэффициенты квадратичной
|
|
n |
|
формы − второго дифференциала |
d 2 |
f = ∑hij (X )dxi dxj |
функции |
i, j=1
y = f (X ).
Теорема 7.1 (достаточное условие выпуклости). Пусть M Rn
– выпуклое открытое множество, y = f (X ) – функция, имеющая в
Mнепрерывные частные производные 2-го порядка. Тогда
1.Функция y = f (X ) является строго выпуклой (строго во-
гнутой) на M , если квадратичная форма d 2 f (X ) положительно (отрицательно) определена.
2.Функция y = f (X ) является выпуклой (вогнутой) на M ,
если квадратичная форма d 2 f (X ) неотрицательно (неположи-
тельно) определена.
Обозначим
93
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂2 f |
|
∂2 f |
|
∂2 |
f |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x2 |
|
∂x ∂x |
∂x ∂x |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
∂2 f |
∂2 f |
|
|
|
|
|
|
1 |
|
|
1 |
|
2 |
|
|
1 |
i |
|
||||
|
|
|
|
|
|
|
|
|
|
|
∂2 f |
|
|
∂2 |
|
2f |
|
∂2 f |
|
|
||||||||
|
∂2 f |
|
|
|
|
∂x12 |
∂x1∂x2 |
|
|
|
|
|
|
|
|
|
||||||||||||
∆ = |
,∆ |
|
= |
|
|
, ,∆ |
i |
= |
|
∂x |
∂x |
|
∂x |
∂x ∂x |
. |
|||||||||||||
∂x12 |
|
|
∂2 f |
∂2 f |
|
|
|
|
|
|
||||||||||||||||||
1 |
|
2 |
|
|
|
|
|
|
2 |
1 |
|
|
|
|
2 |
|
|
2 |
i |
|
||||||||
|
|
|
|
|
|
|
|
|
∂x2 |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
∂x |
∂x |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
2 |
1 |
|
2 |
|
|
|
|
|
∂2 f |
|
∂2 f |
|
∂2 f |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x ∂x |
|
|
∂x ∂x |
|
|
∂x2 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
1 |
|
|
|
i |
|
2 |
|
|
|
i |
|
Тогда в силу критерия Сильвестра достаточное условие выпуклости для строго выпуклых (вогнутых) функций можно переформулировать следующим образом.
Теорема 7.2. При выполнении условий теоремы 7.1 функция y = f (X ) является строго выпуклой в M , если в каждой точке об-
ласти
∆i > 0, i =1,2, n.
Функция y = f (X ) является строго вогнутой в M , если в к а-
ждой точке области
∆1 < 0 , ∆2 > 0 , ∆3 < 0, …, (−1)n ∆n > 0 ,
т.е.
(−1)i ∆i > 0, i =1,2, n.
Для функции двух переменных z = f (x, y) следствие можно
сформулировать в следующей форме.
Следствие 7.1. Пусть M – выпуклое открытое множество на плоскости (x, y) и z = f (x, y) имеет в M непрерывные частные
производные 2-го порядка. Тогда функция z = f (x, y) является строго выпуклой (вогнутой) на множестве M , если в каждой точке
(x, y) M выполняются условия
1.fxx′′ (x, y)> 0 (fxx′′ (x, y)< 0).
∆f (x, y)= |
|
f ′′ (x, y) |
f ′′ |
(x, y) |
|
= fxx′′ (x, y) fyy′′ |
(x, y)−(fxy′′ (x, y)) |
2 |
> 0. |
|
|
||||||||
|
xx |
xy |
(x, y) |
|
|
||||
|
|
fxy′′ (x, y) |
fyy′′ |
|
|
|
|
|
Пример 7.1. Рассмотрим функцию
f (x, y, z)= 2x2 −4xz + 4y2 −8yz +9z2 + 4x +8y −2 z .
Эта функция дважды дифференцируема и ее частные производные fx′ = 4x −4z +4 = 0, fy′ =8y −8z +8, fz′= −4x −8y +18z −20, fxx′′ = 4 ,
94
fxy′′ = 0, |
fxz′′ = −4, |
fyy′′ |
=8, |
fyz′′ = −8 , |
|
|
fzz′′ =18 . |
Матрица |
Гессе равна |
||||||||||||||||||||
|
4 |
0 |
−4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 0 |
|
|
|
|
||||||
H = |
0 8 −8 |
|
|
|
и |
|
∆ = 4 > 0, |
∆ |
2 |
= |
= 32 > 0, |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
0 |
8 |
|
|
|
|
|||||
|
−4 |
−8 |
18 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
4 |
0 |
−4 |
|
=192 > 0. Поэтому функция является строго выпук- |
|||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||
∆3 = |
|
0 |
8 |
−8 |
|
||||||||||||||||||||||||
|
−4 −8 18 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
лой. |
Пример 7.2. Рассмотрим теперь производственную функцию |
||||||||||||||||||||||||||||
|
|
||||||||||||||||||||||||||||
Кобба-Дугласа Q(K, L)= aK β L1−β |
(здесь a > 0 , 0 < β <1). Последова- |
||||||||||||||||||||||||||||
тельно находим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Q′ = aβK β−1L1−β , |
Q′ = a(1 |
−β)K β L−β , Q′′ |
= −aβ (1−β)K β−2 L1−β , |
|
|
||||||||||||||||||||||||
K |
|
|
|
|
|
L |
|
|
|
|
|
|
|
|
KK |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
′′ |
|
|
|
|
|
|
β−1 −β |
|
′′ |
|
|
|
−β)K |
β −β−1 |
. |
|||||
|
|
|
|
|
|
|
|
QKL = aβ (1−β)K |
|
L |
, QLL = −aβ (1 |
L |
|||||||||||||||||
Заметим, что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
Q |
′′ |
Q |
′′ |
|
|
−aβ 1−β |
) |
K |
β−2 1−β |
−aβ 1−β |
) |
K |
β−1 −β |
|
|
|||||||||||
|
|
|
|
|
|
|
= |
L |
|
|
L |
|
|
||||||||||||||||
H (Q)= |
KK |
KL |
|
( |
|
|
|
|
β−1 −β |
( |
|
|
|
|
β −β−1 |
|
|
||||||||||||
|
|
|
|
′′ |
|
′′ |
|
|
−aβ (1−β)K L |
|
−aβ (1−β)K L |
|
|||||||||||||||||
|
|
|
QKL |
QLL |
|
|
|||||||||||||||||||||||
и ∆ = −aβ (1−β)K β−2 L1−β |
< 0, ∆ |
2 |
= 0 . Поэтому достаточное условие |
||||||||||||||||||||||||||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
не дает ответа на вопрос о выпуклости исходной функции. С другой
стороны, второй дифференциал d 2Q = QKK′′ d 2 K +2QKL′′ dKdL +QLL′′ d 2 L легко преобразуется к виду
d 2Q = −aβ (1−β)K β−2 L−β−1 (L2d 2 K −2KLdKdL + K 2d 2 L)=
= −aβ (1−β)K β−2 L−β−1 (LdK − KdL)2 ,
т.е. является неположительно определенной квадратичной формой. Таким образом, функция Кобба-Дугласа является вогнутой функцией.
Пусть функция y = f (X ) определена на некотором множестве
M Rn . Тогда точка X0 называется точкой локального условного максимума (минимума) функции y = f (X ), если f (X )≤ f (X0 ) (или f (X )≥ f (X0 )) для всех точек X M , достаточно близких к X0 .
95