Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Семестровая 1 (Методы оптимальных решений) / Методы оптимальных решений учебное пособие

.pdf
Скачиваний:
1246
Добавлен:
19.04.2015
Размер:
6.01 Mб
Скачать

дится сопряженная переменная 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