Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Математические модели в экономике..pdf
Скачиваний:
22
Добавлен:
05.02.2023
Размер:
2.41 Mб
Скачать

35

m

Функция Ф(X , ) F(X ) i fi X называется функцией Лагранжа и

i 1

коэффициенты i множителями Лагранжа.

Можно доказать, что необходимым условием экстремума исходной задачи является обращение в нуль всех частных производных функции Лагранжа.

· · · · · · · · · · · · · · · · · · · · · · ·

Пример 3.2 · · · · · · · · · · · · · · · · · · · · · · ·

Например,

при поиске максимума F X1, X2 X1 X2

при единственном

условии X 2

X 2

1

строится функция Лагранжа:

 

 

1

2

 

 

 

 

 

 

 

 

 

 

 

Ф(X , ) X1 X2 X12 X22 1 .

 

Строим систему уравнений:

 

 

 

 

 

Ф 1 2 X 0,

Ф 1 2 X 0,

Ф

X 2 X 2 1 0,

 

X1

 

 

1

X2

2

 

1

2

 

 

 

 

 

 

решение которой дает 1

2 , X1 X2 и экстремальные значения целе-

вой функции 2 и

2.

 

 

 

 

 

 

Для определения типа найденного экстремума можно построить матрицу вторых производных F X , вычисленных в экстремальной точке, и определить знаки главных ее миноров. Если все они положительны, то найден минимум; если они чередуются, начиная с минуса, то найден максимум (правило Сильве-

стра).

· · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · ·

Сам по себе метод множителей Лагранжа не дает существенного эффекта из-за необходимости решать, как правило, нелинейную систему уравнений; не гарантирует тип отыскиваемого экстремума, кроме глобальных дает и множество локальных экстремумов, но полезен при генерации идей и создании методов нелинейного программирования.

3.2 Теорема Куна – Таккера

Пусть стоит задача минимизации F X при условиях: fi X 0 (i 1...m),

 

X 0,

(3.1)

где X – n-мерный вектор, F X и

fi X – выпуклые функции, что обеспечива-

ет выпуклость множества планов и единственность искомого

минимума

36

(напомним, что мы считаем функцию выпуклой в некоторой точке, если главные миноры матрицы вторых производных положительны) [15].

Введем функцию Лагранжа:

m

Ф(X , ) F(X ) i fi X .

i 1

· · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · ·

Теорема Куна – Таккера утверждает, что вектор X * 0 явля-

ется решением поставленной задачи тогда и только тогда, когда существует вектор λ* 0, такой, что при всех X 0, λ 0

Ф(X *,λ) Ф(X **) Ф(X *).

Так как функция Лагранжа в точке X *, λ* принимает мини-

мум по X и максимум по , то эта точка называется седловой и эта теорема называется теоремой о седловой точке или теоремой о минимаксе.

· · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · ·

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

ности, т. е. существования хотя бы одной допустимой точки

X , где fi X 0

при всех i.

 

 

 

 

 

 

Если F X и

fi X дифференцируемы,

то условия теоремы эквивалент-

ны «локальным» условиям, утверждающим, что в точке X *, λ*

 

 

Ф 0,

X T Ф

0,

X 0,

(3.2)

 

 

X

X

 

 

 

 

 

Ф

Ф

 

 

 

 

 

λ 0,

λT λ

0,

λ 0.

(3.3)

Остановимся на частных случаях задачи.

 

 

Если в поставленной задаче отсутствуют требования неотрицательности

(3.1), то путем замены X разностью двух неотрицательных векторов можно

показать, что условие (3.2) упрощается к виду

Ф 0 (если нет требования не-

 

 

 

 

 

X

 

отрицательности для одной из компонент X , то обратится в нуль производная

Ф(X ) по соответствующей переменной).

 

 

 

Если

fi X 0 при некотором i, то заменой этого равенства системой не-

равенств

fi X 0,

fi X 0 обнаруживаем, что в (3.3)

соответствующая

37

производная обращается в нуль и исчезает условие неотрицательности по i ;

если fi X 0 при всех i, то (3.3) упрощается к виду Ф 0.

 

 

λ

 

Как частный случай условий Куна – Таккера могут быть построены двой-

ственные задачи линейного программирования.

 

Так, при минимизации CT X

при условиях AX B,

X 0 функция Ла-

гранжа имеет вид:

 

 

Ф(X ,λ) CT X λT (AX B)

 

и условия Куна – Таккера:

 

 

Ф C AT λ 0, X T Ф X T С AT λ 0, X 0,

X

X

 

Ф B AX 0,λ

λT

Ф

λT B AX 0, λ 0.

λ

Отсюда получаем с учетом X T C CT X , λT B BT λ, X T AT λ λT AX , что в седловой точке достигается минимум по X и максимум по λ, причем

CT X BT λ,

AX B, AT λ C,

X 0, λ 0.

3.3 Квадратичное программирование. Метод Вулфа – Фрэнка

Рассмотрим задачу минимизации квадратичной функции n переменных:

F(X ) CT X X T DX ,

(3.4)

при линейных ограничениях

 

AX B, X 0,

 

где A – матрица размерности m на n, C, X – n-мерные векторы, B – m-мерный вектор, D – положительно-определенная n-мерная квадратная матрица (матрица называется положительно определенной, если положительны ее главные мино-

ры) [16].

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так, целевая функция

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F(X ) 7X

1

3X

2

22X 2

8X

X

2

X 2

10X

X

3

12X

2

X

3

40X 2

 

 

1

1

 

 

2

 

1

 

 

 

3

представится в виде (3.4), где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

22

4

5

 

 

 

 

 

 

 

 

 

 

 

 

7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

3

, D

 

4

 

1

6

.

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

5

 

6

40

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

38

Положительная определенность матрицы D и линейность ограничений (множество планов выпукло) позволяют использовать теорему Куна – Таккера.

Функция Лагранжа здесь имеет вид:

Ф(X ,λ) CT X X T DX λT (AX B)

и условия Куна – Таккера приводятся к форме:

Ф

 

 

 

 

 

 

 

Ф

 

 

 

 

 

 

 

X C

2DX AT λ

0,

X T

X 0, X 0,

(3.5)

Ф

AX B

0,

λT

Ф

0,

 

λ 0.

(3.6)

λ

λ

 

Если ввести векторы ослабляющих переменных, то (3.5)–(3.6) примут

вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C 2DX AT λ V 0, X TV 0, X 0, V 0;

 

AX B Y 0,

λTY 0, λ 0, Y 0.

 

С учетом неотрицательности переменных можно поставить задачу мини-

мизации (до нуля) функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(X ,λ,Y,V ) X TV λTY

 

 

 

 

 

 

 

при условиях:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2DX AT λ V C,

 

 

 

 

 

 

(3.7)

 

 

AX Y B,

 

 

 

 

 

 

 

(3.8)

 

 

X ,λ,Y,V 0.

 

 

 

 

 

 

 

(3.9)

Если обозначить

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W T X T T ,YT ,V T ,

 

 

 

 

 

 

 

то (3.7)–(3.9) приведутся к виду:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

RW S,

W 0,

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

0

E

 

, S

 

C

 

 

 

.

 

2D AT

 

 

 

 

 

 

 

A 0

E

 

0

 

 

 

 

B

 

 

 

 

 

Метод П. Вулфа и М. Фрэнка сводит решение задачи к форме, допускающей применение симплексной процедуры.

Здесь отыскивается некоторый опорный план W 0. Если g W 0 0, то этот план оптимален. В противном случае отыскивается градиент

grad g W 0 V T ,YT T , X T

39

и его компоненты на один шаг симплексного преобразования (перехода к другому опорному плану) принимаются за коэффициенты «линейной формы». После выбора нового опорного плана выполняются вышеописанные рассуждения.

· · · · · · · · · · · · · · · · · · · · · · ·

 

Пример 3.3

· · · · · · · · · · · · · · · · · · · · · · ·

Рассмотрим простой пример.

 

 

 

 

Минимизировать

 

 

 

 

 

 

F(X , X

) 4X 6X

2

X 2

3X 2

1

2

1

1

 

2

при условиях:

X1 2X2 4,

X1, X2 0.

Ставим задачу минимизации до нуля функции

g(X1, X2,λ,Y,V1,V2 ) X1V1 X2V2 λY

при условиях (3.8), где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

0

1

0

–1

0

 

4

R =

0

6

2

0

0

–1

S =

6

 

1

2

0

1

0

0

 

4

Ищем начальный опорный план методом искусственного базиса (не обращая внимания на коэффициенты, меньшие M) (табл. 3.1):

Таблица 3.1 – Начальный опорный план

C

 

Ба-

План

 

 

 

 

 

 

 

M

M

баз

 

зис

W

X1

X 2

λ

Y

V1

V2

Z1

Z2

M

Z1

4

2

 

 

1

0

–1 0

1

0

M

Z2

6

0

6

 

2

0

0 –1

0

1

 

 

Y

4

1

2

 

0

1

0

0

0

0

 

k

10M

2M 6M

 

3M

0

–M –M

0

0

C

 

Ба-

План

 

 

 

 

 

 

 

 

 

M

баз

 

зис

W

X1

X 2

λ

Y

V1

V2

 

 

Z1

M

 

Z1

1

2

–3

0

0

–1

1/2

 

 

1

 

 

λ

3

0

3

 

1

0

0

–1/2

 

 

0

 

 

Y

4

1

2

 

0

1

0

0

 

 

0

 

k

M

2M –3M

0

0

–M

M/2

 

 

0

Находим начальный опорный план исходной задачи:

X1 X2 V1 0, V2 2, Y 4, λ 4,

40

для которого g(W ) 0 0 0 2 4 4 16 0.

Отыскиваем градиент grad g(W ) 0,2,4,4,0,0 и берем его компоненты в качестве коэффициентов «линейной формы» (нормали к поверхности, определяемой функцией g(W ) , в выбранной точке) (табл. 3.2):

Таблица 3.2 – Поиск решения

C

 

Ба-

План

0

2

4

4

0

0

баз

 

зис

W

X1

X 2

λ

Y

V1

V2

0

 

V2

2

4

–6

0

0

–2

1

4

 

λ

4

2

0

1

0

–1

0

4

 

Y

4

1

2

0

1

0

0

 

k

32

12

6

0

0

–4

0

C

 

Ба-

План

0

0

7/2

3

1/2

0

баз

 

зис

W

X1

X 2

λ

Y

V1

V2

0

 

X1

1/2

1

–3/2

0

0

–1/2

1/4

7/2

 

λ

3

0

3

1

0

0

–1/2

3

 

Y

7/2

0

7/2

0

1

1/2

–1/4

 

k

21

0

21

0

0

1

–5/2

C

 

Ба-

План

0

0

0

0

2

1

баз

 

зис

W

X1

X 2

λ

Y

V1

V2

0

 

X1

2

1

0

1/2

0

–1/2

0

0

 

X 2

1

0

1

1/3

0

0

–1/6

0

 

Y

0

0

0

–7/6

1

1/2

1/3

 

k

0

 

 

 

 

 

 

Получен оптимальный план X (2,1).

 

Заметим, что при ограничениях

AX B,

X 0 исчезнет переменная Y ,

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

на λ и

потребуется минимизировать

X TV при условиях:

 

 

2DX AT λ V C,

AX B,

 

X ,V 0.

 

При

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

та можно

элементарными приемами прийти к нему (вместо

максимизации

F(X ) искать минимум F(X ) , условие AX B заменить

на B AX 0

ит. п.).

·· · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · · ·