Семестровая 1 (Методы оптимальных решений) / Методы оптимальных решений учебное пособие
.pdfдится сопряженная переменная y(t), строится функция Гамильтона (12.2.7), определяются функции u(t), x(t), y(t), удовлетворяющие условиям (12.2.8)—(12.2.10), и из участков управлений u(t), доставляющих гамильтониану максимум при каждом t, формируется оптимальная управляющая траектория.
Принцип максимума дает необходимые, но, вообще говоря, не достаточные условия для существования максимума в задаче оптимального управления.
Однако в о п р е д е л е н н ы х ч а с т н ы х с л у ч а я х условия принципа максимума являются необходимыми и достаточными услови-
ями оптимальности.
Приведем без доказательства два таких частных случая:
∙когда гамильтониан линеен относительно управляющих параметров
(д о с т а т о ч н ы е у с л о в и я Р о з о н о е р а);
∙когда максимум гамильтониана представляет собой выпуклую
вверх функцию по фазовым координатам (д о с т а т о ч н ы е
ус л о в и я М а н г а с а р я н а).
§12.3. МОДЕЛИРОВАНИЕ ОПТИМАЛЬНОГОЭКОНОМИЧЕСКОГО РОСТА
Опишем модель национальной экономики, предложенную в 1956 г. Р. Солоу, Нобелевским лауреатом 1987 г. в области экономики.
В замкнутой односекторной экономической системе производится один универсальный продукт, который может как потребляться, так и инвестироваться. О с н о в н ы е п р е д п о л о ж е н и я м о д е л и С о л о у состоят в постоянстве темпа прироста числа занятых, износа основных производственных фондов и нормы накопления, отсутствии лага (т. е. запаздывания) капиталовложений.
|
Состояние экономики в момент времени t определяется следую- |
щими показателями: |
|
∙ |
в а л о в ы м в ы п у с к о м X (t) ; |
∙к а п и т а л о м (основными фондами) K (t) ;
∙ |
ч и с л о м з а н я т ы х в производственной сфере L(t) ; |
∙ |
в а л о в ы м и и н в е с т и ц и я м и I (t ) ; |
∙ |
ф о н д о м н е п р о и з в о д с т в е н н о г о п о т р е б л е н и я C (t) . |
Пусть г о д о в о й т е м п п р и р о с т а ч и с л а з а н я т ы х составляет ν, тогда за промежуток времени dt численность занятых изменяется на величину dL = νL(t)dt , значит, для L(t) можно записать дифференциальное уравнение
dL = νL(t) , dt
341
решением которого является функция
L(t) = L0 eνt ,
где L0 — число занятых в начальный момент времени.
Пусть за год выбывает (изнашивается и приходит в негодность) доля μ основных производственных фондов, н о р м а н а к о п л е н и я составляет ρ, а годовой валовый внутренний продукт определяется л и -
н е й н о - о д н о р о д н о й неоклассической |
производственной |
функцией |
||||||||
X = F (K , L) . Тогда и з н о с и и н в е с т и ц и и |
в расчете на год равны |
|||||||||
μK (t) и I (t) = ρX (t) = ρF (K (t), L(t)) |
соответственно, лаг капиталовложе- |
|||||||||
ний отсутствует, значит, п р и р о с т |
ф о н д о в за промежуток времени dt |
|||||||||
составляет dK = −μK (t)dt + I (t)dt или |
|
|
|
|
||||||
|
|
dK |
= −μK (t) + ρF (K (t), L(t)). |
(12.3.1) |
||||||
|
|
|
|
|||||||
|
|
dt |
|
|
|
|
|
|
|
|
Перепишем это уравнение в виде |
|
|
|
|
||||||
|
dK |
|
|
|
K (t) |
|
|
|||
|
|
|
|
= −μK (t) + ρL(t)F |
|
,1 , |
|
|||
|
|
|
|
|
|
|||||
|
dt |
|
|
|
L(t) |
|
|
где |
мы учли, что F (K , L) = LF (K / L, 1) , поскольку производственная |
функция F (K , L) является линейно-однородной. |
|
|
Перейдем теперь к относительным показателям: |
∙ |
ф о н д о в о о р у ж е н н о с т и k (t) = K (t) / L(t) ; |
∙ |
с р е д н е й п р о и з в о д и т е л ь н о с т и т р уд а x(t) = X (t) / L(t) ; |
∙ |
уд е л ь н ы м и н в е с т и ц и я м i(t) = I (t) / L(t) ; |
∙ |
с р е д н е д у ш е в о м у п о т р е б л е н и ю c(t) = C(t) / L(t) . |
Найдем |
|
|
|
||
|
|
|
|
K |
|
|
dk |
|
d |
|
|
|
|
|
|||
|
= |
|
L |
||
|
|
dt |
|||
|
dt |
по формуле производной частного. Имеем:
|
|
|
|
|
|
K |
dK |
L(t) − |
dL |
|
|
|
|
|||||
|
|
|
dk |
|
d |
|
|
|
|
|
K (t) |
|||||||
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
dt |
dt |
||||||||||||
|
|
|
= |
|
L |
= |
|
|
|
|
|
|
, |
|||||
|
|
|
|
|
|
|
|
|
(L(t))2 |
|
|
|
||||||
|
|
|
dt |
dt |
|
|
|
|
|
|
||||||||
при этом |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
dK |
= −μK (t) + ρL(t)F |
K (t) |
,1 , |
|
dL |
= νL(t) , |
|||||||||||
|
|
|
|
|||||||||||||||
|
dt |
|
|
|
|
|
L(t) |
|
|
|
dt |
342
поэтому
|
|
|
|
|
dK |
|
dL |
|
|
|
|
|
|
|
K (t) |
|
|
|
|
|
|
||
|
|
|
|
|
L(t) − |
K (t) |
|
|
−μK (t) |
+ ρL(t)F |
|
,1 L(t) − νL(t)K (t) |
|
||||||||||
|
|
|
|
|
|
|
|||||||||||||||||
dk = dt |
|
= |
|||||||||||||||||||||
|
dt |
= |
|
|
|
L(t) |
|
|
|
|
|
||||||||||||
dt |
|
|
|
|
(L(t))2 |
|
|
|
|
|
|
|
(L(t))2 |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
K (t) |
|
|
|
|
|
|
|
|
|
|
|||
|
|
−(μ + ν)K (t) + ρL(t)F |
|
|
,1 |
L(t) |
|
|
|
K (t) |
K (t) |
|
|
||||||||||
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
L(t) |
|
|
|
|
|
|
||||||||
= |
|
|
|
|
|
|
|
|
|
|
|
|
= −(μ + ν) |
|
+ ρF |
|
,1 |
= |
|||||
|
|
|
|
|
(L(t)) |
2 |
|
|
|
|
|
L(t) |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L(t) |
|
|
= −(μ + ν)k(t) + ρF (k(t),1),
т. е. для фондовооруженности k(t) справедливо дифференциальное уравнение
dk = −(μ + ν)k(t) + ρF (k(t),1). dt
Рассмотрим в качестве производственной функции функцию Кобба — Дугласа F (K , L) = AK αL1−α , при этом F (k ,1) = Ak α .
Введя обозначения λ = μ + ν , f (k ) = F (k,1) = Ak α , получаем окон-
чательно модель Солоу в относительных показателях:
dkt = −λkt + ρf (k (t)),
dt
k0 = K0 / L0 ,
x(t) = f (k (t)), i(t) = ρf (k (t)), c(t) = (1 −ρ) f (k (t)).
Говорят, что экономика находится на стационарной траектории, если относительные показатели не меняются во времени. Поскольку x(t), i(t), и c(t) являются функциями от k(t), то для того, чтобы экономика находилась на стационарной траектории, необходимо и достаточно постоянства во времени фондовооруженности k(t), т. е.
dk = 0 или −λk + ρf (k ) = 0 . dt
Подставим сюда f (k ) = Ak α , получим условие стационарности траектории:
−λk +ρAkα = 0 .
Вынесем k α за скобку:
k α (−λk1−α + ρA)= 0 .
343
Из последнего уравнения видно, что возможны две стационарные
траектории экономики: вырожденная [когда k = 0, при этом x = Ak α = 0 , i = ρAk α = 0 , c = (1 − ρ) Ak α = 0 ] и невырожденная [когда −λk1−α +ρA = 0 ].
Из условия −λk1−α + ρA = 0 k1−α = ρA / λ следует, что на невырожденной стационарной траектории постоянные значения относительных показателей равны
1
= ρA 1−α k λ
|
|
|
|
ρA |
α |
|
i |
|
|
1−α |
|||
|
= ρA(k |
)α = ρA |
λ |
|
||
|
|
|
|
|
,
,
|
|
|
|
|
|
|
|
ρA |
α |
|
|
|
|
||
x |
|
= A(k |
|
|
|
1−α |
|
|
|
|
|
||||
|
|
|
)α = A |
λ |
, |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
α . |
||||
|
|
|
|
|
|
|
|
|
|
|
|
ρA |
|||
c |
|
= (1 − ρ) A(k |
|
|
|
|
|
1−α |
|
||||||
|
|
)α = (1 − ρ) A |
λ |
. |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Исследуем, что произойдет, если экономика отклонится от стационарной траектории. Изобразим на рис. 12.3.1 графики функций λk и ρAk α [здесь α (0,1) ].
Рассмотрим вначале вырожденную стационарную траекторию
(на ней k(t) = 0). Если k(t) станет чуть больше нуля, то, как видно из рис. 12.3.1, ρAk α > λk , поэтому производная
dk = −λk +ρAkα > 0 , dt
откуда следует, что фондовооруженность k(t) будет возрастать. При этом dk / dt остается положительной при всех k (t) (0, k ) , поэтому вырожденная стационарная траектория является неустойчивой: достаточно малейшего возмущения, и k(t) начинает возрастать в сторону k = k ; при k(t) = k производная dk / dt становится равной нулю, т. е. k(t) kt перестает меняться.
λk, ρAkα
λk =A(k )α
λk
k
0 |
k |
Рис. 12.3.1. Исследование устойчивости стационарных траекторий экономики в модели Солоу
344
Если экономика находится на невырожденной стационарной тра-
ектории k(t) = k , и произошло незначительное отклонение фондовооруженности в л е в о от стационарного значения k , то, как мы уже убедились, k(t) начинает возрастать до тех пор, пока вновь не вернется к значению k . Если же k(t) отклонится от k в п р а в о, то, как показывает рис. 12.3.1, ρA(k (t))α < λk (t) , поэтому производная
dk = −λk(t) +ρA(k(t))α < 0 , dt
значит, k(t) будет убывать до тех пор, пока не станет равной k . Таким образом, невырожденная стационарная траектория
|
|
1 |
|
|
|
|
α |
|
|
|
α |
||
|
ρA |
|
|
|
ρA |
|
|
ρA |
|
|
|||
1−α |
, xt |
1−α |
1−α |
||||||||||
k(t) = |
λ |
|
= A |
λ |
|
, i(t) = ρA |
λ |
|
, |
||||
|
|
|
|
|
|
|
|
|
|
ρA |
α |
|
1−α |
|||
c(t) = (1 −ρ) A |
λ |
|
|
|
|
является у с т о й ч и в о й: при любом отклонении от этой траектории экономика стремится к ней вернуться.
Данная невырожденная стационарная траектория носит название
траектории сбалансированного устойчивого экономического роста: чис-
ленность занятых на ней возрастает экспоненциально: L(t) = L0 eνt (конечно, при положительном темпе прироста занятых ν), а все относительные показатели постоянны, значит, все абсолютные показатели возрастают
пропорционально численности занятых L(t).
Рассмотрим теперь простейшую задачу управления экономикой, которая описывается моделью Солоу: попытаемся подобрать такую норму накопления ρ, чтобы удельное потребление на стационарной траектории сбалансированного устойчивого экономического роста было максимальным.
ЗОЛОТОЕ ПРАВИЛО НАКОПЛЕНИЯ. Чтобы удельное потребление на ста-
ционарной траектории сбалансированного экономического роста было максимальным, норма накопления ρ должна быть равна эластичности выпуска по фондам α .
Доказательство. Рассмотрим удельное потребление на стационарной траектории c как функцию нормы накопления:
|
|
|
|
ρA |
α |
|
c |
|
|
1−α |
|||
|
= c |
(ρ) = (1 −ρ) A |
λ |
|
||
|
|
|
|
|
и поставим задачу определения такой нормы накопления ρ, чтобы
c (ρ) → max
или, расписывая c (ρ) подробно,
345
|
ρA |
α |
|
|
1−α |
→ max . |
|||
(1 − ρ) A |
λ |
|
||
|
|
|
В точке максимума первая производная должна быть равна нулю (или не существовать), а вторая производная должна быть отрицательной. В данном случае имеем:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ρA |
|
|
|
|
α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1−α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(ρ) |
|
|
|
|
d |
(1 − ρ) A |
λ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dc |
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dρ |
|
|
|
|
|
|
|
|
) |
|
|
dρ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
) |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
α |
|
|
|
|
|
|
|
|
|
|
|
|
α |
|
|
|
|
|
|
|
|
|
|
|
α |
|
|
|
|
|
|
|
|
|
α |
|
|
|
α |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
A |
|
|
|
|
|
d ((1 − ρ)ρ |
|
|
|
|
|
A |
|
|
|
|
d (ρ |
|
|
|
|
− ρ1+ |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
1−α |
|
|
|
|
1−α |
1−α |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
1−α |
1−α |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
= A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
dρ |
|
|
|
|
|
|
|
λ |
|
|
|
|
|
|
|
|
|
|
dρ |
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
λ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
α |
|
|
|
|
|
α |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
A |
|
d (ρ |
|
|
|
) |
|
|
|
d (ρ |
|
|
) |
|
|
|
|
A |
|
|
|
|
|
|
|
α |
|
|
|
|
α |
−1 |
|
|
|
|
|
1 |
−1 |
||||||||||||||||||||||||||||
|
|
1−α |
|
|
|
1−α |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|||||||||||||||||||||||||||||||||||||||||||||||||
1−α |
|
|
|
|
|
|
1−α |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
= A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
= |
A |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ρ1−α |
|
− |
|
|
|
|
ρ1−α |
= |
|||||||||||||
|
|
|
dρ |
|
|
dρ |
|
|
|
|
|
|
1 − α |
|
|
− α |
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
λ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
λ |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
A |
|
α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
|
α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
1−α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1−α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
|
|
α |
|
|
|
|
|
|
|
α |
A |
|
|
|
|
|
|
α |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
= |
|
|
|
λ |
|
|
(αρ |
|
−1 − ρ |
|
) = |
|
|
|
λ |
|
|
|
ρ |
|
−1 (α − ρ). |
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
1−α |
1−α |
|
|
|
1−α |
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
1 − α |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 − α |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
Видим, что dc (ρ) / dρ = 0 |
при ρ = 0 и при ρ = α . Предлагаем чита- |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
телю самостоятельно убедиться, |
|
что в точке ρ = 0 |
вторая производная |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
d 2c (ρ) / dρ2 |
> 0 , т. е. точка ρ = 0 является точкой м и н и м у м а удельно- |
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
го потребления на стационарной траектории, |
|
|
а в точке ρ = α вторая про- |
изводная d 2c (ρ) / dρ2 < 0 , т. е. точка ρ = α является точкой м а к с и м у - м а удельного потребления, что и требовалось доказать.
Этот результат получен в 1966 г. Э. Фелпсом.
ПРИМЕР 12.3.1. Даны значения параметров A = 103 и α = 0,5 производственной функции Кобба — Дугласа. В модели Солоу с этой производственной функцией требуется рассчитать значения фондовооруженности, производительности труда и удельного потребления на стационарной траектории сбалансированного устойчивого экономического роста, на которой норма накопления равна ρ = 0, 2 , коэффициент выбытия основных производственных фондов за год составляет μ = 0, 2 , а годовой темп прироста численности занятых равен ν = 0, 05 . Сравнить полученное значение удельного потребления с оптимальным.
РЕШЕНИЕ. На стационарной траектории, соответствующей норме накопления ρ = 0, 2 , фондовооруженность
346
|
|
|
|
1 |
|
|
|
1 |
|
|
|
||
|
|
|
rA |
|
|
|
0, 2 ×103 |
|
|
|
|
|
|
|
|
1−α |
|
1−0,5 |
|
4 |
|
||||||
k |
|
= |
|
|
= |
|
|
|
= 64 ×10 |
|
, |
||
|
|
|
|
|
|||||||||
|
|
m + n |
0, 2 + 0,05 |
|
|
|
|
средняя производительность труда
x = A(k )α =10(64 ×104 )0,5 = 8 ×105 ,
удельное потребление
c = (1 - r)x = (1 - 0, 2)8 ×105 = 64 ×104 .
Согласно золотому правилу накопления, для того чтобы на стационарной траектории сбалансированного устойчивого экономического роста удельное потребление было максимальным, нужно выбрать норму накопления ρ равной эластичности выпуска по фондам α , т. е. в рассматриваемом примере максимум удельного потребления на стационарной траектории достигается при r = a = 0,5 . При этом
|
|
1 |
|
|
|
1 |
|
||
|
|
|
|
|
3 |
|
|
|
|
1−α |
0,5 ×10 |
1−0,5 |
|
||||||
k* = |
r A |
|
= |
|
|
= 4 ×106 , |
|||
|
|
|
|||||||
m + n |
0, 2 + 0,05 |
|
|
c = (1 - r )x = (1 - r ) A(k )α = (1 - 0,5) ×103 × (4 ×106 )0,5 =106 64 ×104.
Видим, что оптимальный выбор нормы накопления приводит к существенному увеличению удельного потребления на стационарной траектории — более чем в полтора раза!
Будем теперь считать, что норма накопления не является константой, а изменяется во времени. Для этого удобнее всего представить инвестиции в основные производственные фонды в виде разности валового продукта и валового потребления. Уравнение (12.3.1) в таком случае
превращается в следующее: |
|
dK = -mK (t) + F ( K (t), L(t)) - С(t) . |
(12.3.2) |
dt
Переход в (12.3.2) к относительным показателям дает
dk = -lk(t) + f (k(t)) - c(t). dt
Модель оптимального экономического роста предполагает максимизацию интегрального дисконтированного (по непрерывной ставке δ) удельного потребления
T |
|
J = ∫c(t) e−δt dt ® max |
(12.3.3) |
0 |
|
347
при условиях
|
dk |
= −λk (t) + f (k (t))−c(t), |
|||||
|
|||||||
|
|
|
|
|
|
|
|
dt |
|
|
|
|
|
||
k |
0 |
= K |
0 |
/ L , |
(12.3.4) |
||
|
|
|
0 |
|
|||
x(t) = f (k (t)), |
i(t) = f (k (t))− c(t). |
||||||
|
|
|
|
|
|
|
Горизонт планирования T в данной задаче может быть конечным или бесконечным. В случае планирования на конечный период в максимизируемый целевой функционал (12.3.3) целесообразно добавить слагаемое, накладывающее условие на минимальную фондовооруженность к концу периода [0, T].
Модель (12.3.3)—(12.3.4) представляет собой задачу оптимального управления, фазовой переменной в которой выступает фондовооруженность k(t), а управляющей переменной— удельное потребление c(t). На управление накладывается очевидное ограничение:
|
|
|
|
|
|
|
0 c c(t) f (k (t)) |
|
|||||
(здесь с — |
минимально допустимое удельное потребление). |
|
|||||||||||
Введем одну сопряженную переменную y(t), соответствующую |
|||||||||||||
единственной фазовой переменной k(t), и построим гамильтониан |
|
||||||||||||
|
H (t, k (t), c(t), y(t) ) = c(t) e−δt + y(t) (−λk (t) + f (k (t ))− c(t)). |
|
|||||||||||
Уравнение для сопряженной переменной имеет вид |
|
||||||||||||
|
|
|
|
dy |
= (λ − f ' (k(t)))y(t), |
y(+∞) = 0 . |
(12.3.5) |
||||||
|
|
|
|
|
|||||||||
|
|
|
|
dt |
|
|
|||||||
Сопряженную переменную удобно |
представить в виде |
y(t) = |
|||||||||||
= m(t) e−δt |
и подставить |
|
|
||||||||||
|
|
|
|
|
|
|
dy |
= |
dm |
e−δt −δm(t)e−δt |
|
||
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
dt dt |
|
|
||||
в уравнение (12.3.5): |
|
|
|||||||||||
|
|
dm |
e−δt −δm(t)e−δt = (λ − f ' (k(t)))m(t)e−δt |
|
|||||||||
|
|
|
|
||||||||||
|
|
dt |
|
|
|||||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dm |
= (λ + δ− f ' (k(t)))m(t) . |
(12.3.6) |
|||||
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
dt |
|
|
348
Общее решение уравнения (12.3.6) имеет вид m(t) = e−∫( f ′(k (t ))−(λ+δ))dt , поэтому, очевидно, m(t) > 0.
Гамильтониан
H(t, k (t), c(t), y(t)) = c(t) e−δt + m(t) e−δt (−λk (t) + f (k (t)) − c(t)) =
=((1 − m(t))c(t) + m(t)( f (k (t))− λk (t)))e−δt
при m(t) ¹ 1 линейно зависит от управления c(t), поэтому максимум гамильтониана по управляющей переменной c(t) может достигаться только
на концах отрезка |
c, f (k(t)) |
. Таким образом, оптимальное управление |
||
|
|
|
|
|
при m(t) ¹ 1 определяется так: |
|
|
||
|
|
c (t) = |
c, |
m(t) >1, |
|
|
|
|
|
|
|
|
f (k (t)), m(t) <1. |
При m(t) = 1 гамильтониан
H(t, k (t), c(t), y(t)) = ( f (k (t))− λk (t))e−δt ,
иего максимум по фазовой переменной k(t) достигается при f ′(k (t)) = λ .
Если положить c(t) = f (k (t))− λk (t) , то траектория, соответствующая та-
кому управлению, будет стационарной, так как на ней
dk = −λk(t) + f (k(t))−c(t) = 0 dt
Отсюда следует, что на такой траектории f ′(k (t)) = λ .
Таким образом, при m(t) =1 оптимальное управление определяется как c (t) = f (k (t))−λk (t)
Окончательно получаем:
c, |
|
|
m(t) >1, |
|
|
c (t) = f (k (t))− λk (t), |
m(t) =1, |
(12.3.7) |
|||
f (k (t)), |
m(t) <1. |
|
|||
В случае производственной функции Кобба — |
Дугласа |
||||
|
α |
′ |
αA |
|
|
f (k) = Ak |
|
, f (k) = |
|
, |
|
|
k1−α |
|
Поэтому условие f ′(k ) = λ определяет
349
|
|
|
1 |
|
|
|
|
αA |
|
|
|
k |
1−α |
||||
= |
λ |
. |
|||
|
|
|
Сформулируем теперь правила оптимального управления эко-
номикой в соответствии с (12.3.7). Пока фондовооруженность меньше k , следует ограничить удельное потребление на минимально допустимом уровне c . Как только фондовооруженность достигнет стационарного зна-
чения k , следует скачком увеличить удельное потребление с c до f (k ) − λk . Если же фондовооруженность больше стационарного значе-
ния, то на потребление следует отправлять весь выпуск: c = f (k ) , и когда за счет проедания фондов экономика выйдет на стационарную траекторию k (t) = k , следует уменьшить удельное потребление до f (k ) − λk .
На стационарной траектории при этом обеспечивается поддержание фондовооруженности и удельного потребления на постоянном уровне: k (t) = k , c(t) = f (k ) − λk .
§ 12.4. МОДЕЛИРОВАНИЕДИНАМИКИ ВЗАИМОДЕЙСТВИЯ РАЗРАБОТЧИКОВ КОММЕРЧЕСКОГО
ИНЕКОММЕРЧЕСКОГОПРОГРАММНОГООБЕСПЕЧЕНИЯ
Вэтом параграфе излагаются результаты анализа смешанной дуополии производителей коммерческого программного обеспечения (для определенности, серверной операционной системы Microsoft Windows) и некоммерческого (для определенности, Linux), полученные автором в 2008 г. путем применения аппарата теории оптимального управления.
Будем предполагать рынок программного обеспечения л и н е й н о
ра с т у щ и м с темпом роста a: в единицу времени на рынок приходят a новых пользователей, т. е. суммарное число пользователей на рынке серверных операционных систем к моменту времени t равно N(t) = N0 + at.
Предположим, что каждый новый пользователь выбирает один и только один продукт: или приобретает лицензионную копию Windows, или бесплатно скачивает копию Linux.
Через nW (t) и nL (t) обозначим суммарное число пользователей, ис-
пользующих на момент t операционные системы Windows и Linux соответственно.
Если обозначить q (t ) д о л ю н о в ы х п о л ь з о в а т е л е й, входящих на рынок в момент t и приобретающих Windows, то доля новых пользователей, приобретающих в этот момент Linux, составит 1 − q(t) , поэтому
t |
|
nW (t) = ∫aq(τ)dτ , |
(12.4.1) |
0 |
|
350