Журнал вычислительной математики и математической физики, 2023, T. 63, № 3, стр. 390-407

Анализ на основе математической модели механизмов стимулирования производственных инвестиций на несовершенном рынке капитала

Н. К. Обросова 123*, А. А. Шананин 123**

1 Федеральный исследовательский центр “Информатика и управление” РАН
119333 Москва, ул. Вавилова, 44, Россия

2 Московский центр фундаментальной и прикладной математики, МГУ им. М.В. Ломоносова
119992 Москва, Ленинские горы, Россия

3 Московский физико-технический институт (НИУ)
141701 Московская обл., Долгопрудный, Институтский пер., 9, Россия

* E-mail: nobrosova@ya.ru
** E-mail: alexshan@ya.ru

Поступила в редакцию 23.06.2022
После доработки 23.06.2022
Принята к публикации 14.11.2022

Полный текст (PDF)

Аннотация

Проблема возобновления рыночных инвестиций в реальном секторе российской экономики тесно связана с состоянием предпринимательской среды в условиях несовершенного рынка капитала в России и проблемой оценки доходности инвестиционных проектов. Трудности с определением показателя доходности в условиях несовершенной денежно-кредитной системы связаны с существенным расхождением процентных ставок по депозитам и кредитам и могут быть преодолены в рамках подхода Кантора–Липмана, который позволяет вычислить показатель доходности пула инвестиционных проектов, доступных инвестору. С точки зрения собственника производства рыночные инвестиции зависят от состояния предпринимательской среды и конкурируют с инвестициями в потребление. Возникает проблема оценки порогового значения показателя доходности, при котором собственнику выгодно отложить потребление в пользу рыночных инвестиций. Мы предлагаем подход к решению этой проблемы в терминах математической модели инвестиционного поведения собственника производства в условиях несовершенного рынка капитала, формализованной в виде задачи оптимального управления с фазовым ограничением на бесконечном горизонте. Решение задачи основано на построении вязкостного решения уравнения Гамильтона–Якоби–Беллмана. Показано, что инвестиционная стратегия собственника производства может существенно зависеть от состояния предпринимательской среды. Результаты исследования задачи позволили предложить подход к объяснению перехода российской экономики из режима восстановительного роста в режим стагнации в конце 2007 г., сопровождавшийся спадом инвестиционной активности в производственной сфере. Библ. 20. Фиг. 7.

Ключевые слова: модель Кантора–Липмана, доходность инвестиций, оптимальное управление, модель инвестиций, несовершенный рынок, вязкостное решение, уравнение Гамильтона–Якоби–Беллмана.

1. ВВЕДЕНИЕ

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

Ввиду исторически сложившихся особенностей российская экономика характеризуется неоднородностью ее секторов, связанной с высокой конкурентоспособностью на мировом рынке продукции добывающих (капиталоемких) отраслей и технологической отсталостью обрабатывающих (трудоемких) отраслей экономики, проигрывающих в конкуренции импортным аналогам (подробнее см. [1]). В постсоветский период экономика России прошла через несколько эволюционных этапов. Однако неоднородность и технологическая отсталость трудоемкого сектора сохранились. Перечислим кратко особенности основных этапов развития российской экономики за последние три десятилетия.

Этап 1. Трансформационный спад (1990–1998 гг.), сопровождавшийся упрощением структуры производства за счет вытеснения неэффективных технологических цепочек и замены их импортными аналогами, что привело к снижению мультипликатора в производстве с 2 до 1.3.

Этап 2. Восстановительный рост (1998–2008 гг.), сопровождавшийся процессами импортовытеснения в ряде отраслей, ростом мультипликатора до 1.8. В этот период важное значение имело бюджетное финансирование перспективных направлений (национальных проектов), имеющих высокую доходность (в основном, в инфраструктурных отраслях).

Этап 3. Период стагнации экономики, начавшийся к концу 2007 г. и продолжающийся до настоящего времени. В этот период в связи с выросшими транзакционными издержками снизилась доходность крупных национальных проектов в инфраструктурных отраслях, замедлились темпы экономического роста и произошел спад инвестиционной активности в производственной сфере.

На фиг. 1 представлены динамика роста в обрабатывающем секторе (левый график) и доля обрабатывающего сектора в валовой добавленной стоимости (ВДС) в России (правый график) на втором и третьем этапе эволюции экономики (источник: сайт Росстата РФ www.rosstat.gov.ru). Анализ статистических рядов подтверждает замедление темпов роста в обрабатывающем секторе после 2007 г., а также сокращение доли обрабатывающего сектора в ВДС, которая до настоящего времени не вернулась к уровню начала 2000-х годов.

Фиг. 1.

Статистические показатели динамики развития обрабатывающего сектора РФ.

Следует отметить, что тенденции общемировой динамики экономического роста и аналогичных показателей для крупных развивающихся экономик существенно отличаются от показателей роста российской экономики России в период стагнации. По данным Всемирного Банка за период 2008–2019 гг. экономика России выросла на 12.4%, в то время как мировая экономика выросла на 34.9%. Данная ситуация проиллюстрирована на фиг. 2, где приведены статистические данные о динамике валового внутреннего продукта (ВВП) не только для российских данных, но и для общемировых показателей, а также для агрегированного показателя трех крупных развивающихся экономик мира (Бразилии, Индии и Китая).

Фиг. 2.

Динамика ВВП крупных развивающихся экономик и всего мира.

Спад инвестиционной активности в обрабатывающем секторе отечественной экономики в период стагнации требует объяснения, учитывающего особенности российской ситуации. Естественным подходом к решению данного вопроса являются математические модели экономики, основанные на описании интересов и логики поведения основных экономических агентов с учетом особенностей сложившихся экономических отношений и позволяющие анализировать косвенные причины изменения стратегий поведения агентов. Значительный опыт таких исследований накоплен научной школой академиков Н.Н. Моисеева и А.А. Петрова (см. [2], [3]). В настоящей статье мы продолжаем исследования в данном направлении и строим математическую модель инвестиционного поведения собственника в производственной сфере, учитывающую особенности рынка инвестиций в России. Результаты исследования модели позволяют предложить подход к объяснению механизмов восстановления роста в обрабатывающем секторе отечественной экономики.

Одной из важных особенностей российских экономических условий является несовершенство рынка капитала, характеризующееся значительным расхождением процентных ставок по кредитам и депозитам финансовых организаций. Наша гипотеза заключается в том, что причиной сокращения рыночных инвестиций на рубеже 2008 г. и перехода в режим стагнации российской экономики стало изменение состояния предпринимательской среды на несовершенном рынке капитала. Для анализа состоятельности данной гипотезы необходимо ответить на следующие вопросы.

• Что такое предпринимательская среда и как ее моделировать в российских условиях?

• Как описывать изменение состояния предпринимательской среды и как построить числовую характеристику этого состояния?

• Как моделировать поведение собственника капитала в зависимости от характеристики состояния предпринимательской среды?

Ответы на эти вопросы неочевидны и требуют учета отмеченных выше особенностей сложившихся экономических отношений.

Описание предпринимательской среды на несовершенном рынке капитала мы предлагаем проводить в рамках подхода Кантора–Липмана, который впервые был предложен в [4], [5] и развит в более поздних публикациях многих авторов (см. [613]). В разд. 1 мы приводим описание модели Кантора–Липмана и даем краткий обзор результатов исследований и развития модели. Эти результаты обосновывают использование модели Кантора–Липмана для описания предпринимательской среды в условиях несовершенного рынка капитала и позволяют вычислить числовую характеристику предпринимательской среды, являющуюся показателем доходности пула всех доступных инвестору инвестиционных проектов. Такой показатель доходности при некоторых дополнительных предположениях может быть вычислен как минимальный положительный корень инвестиционной функции, которая строится в терминах модели. Тем самым мы предлагаем ответ на первые два поставленных выше вопроса.

В разд. 2 мы переходим к рассмотрению третьего вопроса и показываем, что классические модели потребительского поведения типа модели Рамсея (см. [15]) плохо подходят для целей нашего исследования, а именно, не позволяют описать существенную связь между стратегией поведения инвестора на несовершенном рынке капитала и характеристикой предпринимательской среды, вычисленной в терминах модели Кантора–Липмана. Далее в разд. 3 мы предлагаем подход к моделированию поведения инвестора на несовершенном рынке капитала, который позволяет проанализировать такую связь. Мы строим математическую модель инвестиционного поведения собственника производства в сложившейся предпринимательской среде с показателем доходности $r$, который может быть вычислен в терминах подхода Кантора–Липмана и является входным параметром модели. Для описания потребительского поведения собственника мы используем новый подход, который предполагает заинтересованность собственника в накоплении средств потребления (объектов недвижимости, предметов роскоши и т.п.) наряду с инвестиционной деятельностью. Модель формализована в виде задачи оптимального управления с фазовым ограничением на бесконечном горизонте. Исследование задачи основано на изучении свойств решения аналогичной задачи оптимального управления на конечном горизонте времени и дальнейшем построении обобщенного (вязкостного) решения уравнения Гамильтона–Якоби–Беллмана, соответствующего исходной задаче. Данная задача оптимального управления проанализирована авторами в [16]. В разд. 4 статьи мы приводим схему исследования задачи и основные результаты. Экономическая интерпретация полученных результатов обсуждается в разд. 5. Результаты исследования модели позволяют проявить зависимость инвестиционного поведения собственника предприятия от начального уровня его капитала и характеристики состояния предпринимательской среды. Показано, что активность инвестиционной деятельности малого и среднего бизнеса в производственном секторе зависит от характеристик предпринимательской среды, сложившейся на несовершенном рынке капитала и рост доходности предпринимательской среды приводит к переключению поведения собственника из режима приоритетного потребления в режим приоритетного инвестирования. В то же время представители крупного бизнеса оказываются нейтральны по отношению к изменению состояния предпринимательской среды.

В разд. 6 мы возвращаемся к описанию предпринимательской среды в терминах подхода Кантора–Липмана применительно к ситуации в России на рубеже 2008 г. На основе результатов исследования модели поведения собственника мы приводим возможное объяснение механизма переключения из режима восстановительного роста в режим стагнации в обрабатывающем (трудоемком) секторе отечественной экономики при ухудшении характеристик предпринимательской среды.

2. ОПИСАНИЕ ПРЕДПРИНИМАТЕЛЬСКОЙ СРЕДЫ В УСЛОВИЯХ НЕСОВЕРШЕННОГО РЫНКА КАПИТАЛА. МОДЕЛЬ КАНТОРА–ЛИПМАНА

В ситуации несовершенного рынка капитала перестают работать традиционно используемые в теории корпоративных финансов подходы к оценке инвестиционных проектов, основанные на расчете доходности с помощью показателя NPV (Net Present Value) – чистой приведенной стоимости денежных потоков, связанных с проектом. В рамках данного подхода инвестиционный проект описывается распределенными во времени финансовыми потоками ${\mathbf{a}} = \left( {{{a}_{0}},{{a}_{1}}, \ldots ,{{a}_{n}}} \right)$, связанными с вложением средств (${{a}_{{\text{i}}}} < 0$) и получением выплат (${{a}_{{\text{i}}}} > 0$) по проекту в $i$-й период времени от начала проекта, $i = 0,~1,...,n$, где $n$ – конечно, т.е. предполагается, что проект полностью завершается через определенный период времени. Если платежи осуществляются через некоторые равные промежутки времени, то NPV проекта с вектором потоков платежей ${\mathbf{a}}$ вычисляется по формуле

(2.1)
$NPV\left( {{\mathbf{a}},r} \right) = {{a}_{0}} + {{a}_{1}}{{e}^{{ - r}}} + \ldots + {{a}_{n}}{{e}^{{ - rn}}},$
где${\text{\;\;}}r{\text{\;}}$ – дефлятор денежных потоков, позволяющий привести все финансовые потоки проекта к начальному моменту времени. В теории корпоративных финансов в качестве такого дефлятора предлагается рассматривать показатель доходности общедоступного для инвесторов альтернативного способа вложения средств. Предполагается, что проект ${\mathbf{a}}$ следует считать рентабельным и поддерживать, если рассчитанный с использованием дефлятора ${\text{\;\;}}r{\text{\;}}$ показатель $NPV\left( {{\mathbf{a}},r} \right) > 0,$ и считать нерентабельным, если $NPV\left( {{\mathbf{a}},r} \right) < 0$. Возникает вопрос о выборе дефлятора денежных потоков $r{\text{\;}}$, так как от его величины существенно зависит оценка рентабельности проектов по критерию NPV. Очевидно, что величина дефлятора $r{\text{\;}}$ заключена между процентными ставками по депозитам ${{r}_{{\text{D}}}}$ и кредитам ${{r}_{{\text{L}}}}$, т.е. ${{r}_{{\text{D}}}} \leqslant r \leqslant $ ${{r}_{{\text{L}}}}$. Поэтому вопрос выбора дефлятора $r{\text{\;}}$ не так важен в условиях совершенного рынка капитала (в странах с развитой рыночной экономикой), так как процентные ставки по кредитам и депозитам в этом случае близки друг к другу и значение $r{\text{\;определяется\;\;достаточно\;точно}}$ $r \approx {{r}_{{\text{D}}}} \approx {{r}_{{\text{L}}}}$ . В ситуации несовершенного рынка капитала, характерного для России, расхождение процентных ставок велико, ${{r}_{{\text{D}}}} < $ ${{r}_{{\text{L}}}}$, и рентабельность проекта существенно зависит от выбора значения дефлятора $r{\text{\;\;}}$ из диапазона ${{r}_{{\text{D}}}} \leqslant r \leqslant $ ${{r}_{{\text{L}}}}$. Поэтому на несовершенном рынке капитала возникает проблема выбора значения дефлятора денежных потоков $r$ инвестиционного проекта. Для определения значения дефлятора следует формализовать понятие общедоступных альтернативных проектов по вложению средств. На несовершенном рынке капитала доходности таких вложений зависят не только от показателей экономической среды, но и от индивидуальных характеристик инвестора, что приводит к их неоднозначности.

Решение этой проблемы было предложено в работах Кантора и Липмана (см. [4], [5]), в которых доходность инвестиционных проектов вычисляется на основе решения задачи управления интенсивностью реализации проекта с целью максимизации дохода инвестора при условии завершения всей инвестиционной деятельности к терминальному моменту времени. Модель Кантора–Липмана позволяет формализовать описание предпринимательской среды на несовершенном рынке капитала как полного пула общедоступных инвесторам проектов по вложению средств и вычислить единую величину дефлятора денежных потоков ${\text{\;}}r,{\text{\;}}$характеризующего этот пул проектов. Значение этого дефлятора отражает доходность альтернативного способа вложения денег в условиях сформировавшейся предпринимательской среды.

В модели Кантора–Липмана предпринимательская среда описывается пулом всех доступных инвестору проектов

${{{\mathbf{a}}}^{m}} = \left( {a_{0}^{m},a_{1}^{m}, \ldots ,a_{n}^{m}} \right),\quad m = 1,{\text{\;}}2,...,M,$
где $n + 1$ – наибольшая продолжительность среди всех проектов (если проект заканчивается раньше, то соответствующие компоненты дополняются нулями). Для решения задачи оценки дефлятора $r$ как доходности альтернативного общедоступного способа вложения денег следует предположить, что все входящие в пул проекты являются стационарными и тиражируемыми. Стационарность проекта означает, что проект является доступным в неизменном виде в любой момент времени. Тиражируемость проекта означает, что в каждый момент времени $t$ проект m может запускаться с произвольной постоянной для данного проекта интенсивностью ${{u}_{m}}\left( t \right) \geqslant 0$. Принимая решение относительно финансовых вложений, инвестор формирует инвестиционную стратегию $\left\{ {{{u}_{m}}\left( t \right)~\,{\text{|}}\,m = 1,2,~.~.~.~,M;~\;t = 0,1,~...,~T - n - 1} \right\}$ при условии, что к моменту времени $T$ вся инвестиционная деятельность должна завершиться. Тогда динамика остатка расчетного счета $s\left( t \right)$ инвестора описывается соотношением
$s\left( {t + 1} \right) = s\left( t \right) + \sum\limits_{m = 1}^M {\sum\limits_{k = 0}^n {a_{k}^{m}{{u}_{m}}} } \left( {t - k} \right),\quad t = 0,1,...,T - 1.$
Заметим, что пул проектов включает, в том числе, проекты депонирования под процент ${{r}_{{\text{D}}}}$ и кредитования под процент ${{r}_{{\text{L}}}}$, если они доступны инвестору. Поэтому дополнительных возможностей занимать и вкладывать деньги у инвестора нет и выполняется условие самофинансирования $s\left( t \right) \geqslant 0,~$ $t = 1,2,...,T$. В модели Кантора–Липмана рассматривается следующая задача об оптимальной стратегии инвестирования:

(2.2)
$s\left( T \right) \to {\text{max,}}$
(2.3)
$s\left( {t + 1} \right) = s\left( t \right) + \sum\limits_{m = 1}^M {\sum\limits_{k = 0}^n {a_{k}^{m}{{u}_{m}}} } \left( {t - k} \right),\quad t = 0,1,...,T - 1,$
(2.4)
$s\left( t \right) \geqslant 0,\quad t = 1,~2,...,T,$
(2.5)
$s\left( 0 \right) = {{x}_{0}} > 0,$
(2.6)
${{u}_{m}}\left( t \right) \geqslant 0,\quad m = 1,~2,...,M;\quad ~t = 0,~1,...,T - n - 1,$
(2.7)
${{u}_{m}}\left( t \right) = 0,\quad m = 1,~2,...,M;\quad t = --n,--n + 1,..., - 1,\quad t = T--n,T--n + 1,...,T - 1.$

Отметим, что условие полного завершения инвестиционной деятельности к терминальному моменту времени приводит к ограничению (2.7). Выполнение условия (2.7) гарантирует, что инвестиционный процесс не является финансовой пирамидой, из которой нельзя выйти, фиксировав доход. Именно это условие гарантирует единственность значения дефлятора денежных потоков проекта в результате решения задачи (подробнее см. [8], [12]).

Обозначим через ${{V}_{T}}$ оптимальное значение функционала в задаче (2.2)–(2.7).

Введем понятие инвестиционной функции $F\left( r \right)$, которая определяется как поточечный максимум показателей NPV (см. (2.1)) всех проектов, входящих в доступный пул проектов инвестора:

(2.8)
$F\left( r \right) = \mathop {{\text{max}}}\limits_{1 \leqslant m \leqslant M} \sum\limits_{k = 0}^n {a_{k}^{m}{{e}^{{ - rk}}}} ,\quad r \in \left[ {0, + \infty } \right)~.$

Справедливо следующее утверждение.

Теорема 1 (Д.Г. Кантор, С.А. Липман, см. [4], [5]).

1. Если $F\left( 0 \right) \leqslant 0$, то для любого $T \geqslant 0~,$ ${{V}_{T}} = {{x}_{0}}$ (пул инвестиционных проектов убыточен).

2. Если $F\left( r \right) > 0$ для всех $r \in \left. {[0,~\infty } \right)$, то существует ${{T}_{0}} > 0$ такое, что для любого $T \geqslant {{T}_{0}}$ справедлива оценка ${{V}_{T}} = + \infty $ (пул инвестиционных проектов допускает арбитраж).

3. Если $F\left( 0 \right) > 0$, существует $\tilde {r} > 0$ такое, что $F\left( {\tilde {r}} \right) < 0$, и $\rho > 0$ – наименьший положительный корень уравнения $F\left( r \right) = 0$, то $\rho $ является доходностью пула инвестиционных проектов, т.е. $\mathop {\lim }\limits_{T \to + \infty } \frac{{\ln \left( {{{V}_{T}}} \right)}}{T} = \rho $.

4. Если $\rho $ – минимальный положительный корень инвестиционной функции $F\left( r \right)$, причем $\rho $ – простой корень, т.е. $\left| {\left\{ {m\left| {\mathop \sum \nolimits_{k = 0}^n a_{k}^{m}{{e}^{{ - \rho k}}} = F\left( \rho \right)} \right.} \right\}} \right| = 1$ и $F{\kern 1pt} '\left( \rho \right) \ne 0$, то существуют положительные константы $0 < c < C$ такие, что $c{{e}^{{\rho T}}} \leqslant {{V}_{T}} \leqslant C{{e}^{{\rho T}}}$.

Из теоремы 1 следует, что в условиях несовершенного рынка капитала темп роста доходов $r$ при реализации пула инвестиционных проектов (дефлятор денежных потоков) может быть найден как минимальный положительный корень инвестиционной функции $F\left( r \right)$. Тем самым теорема Кантора–Липмана позволяет решить вопрос об оценке дефлятора $r$, используемого в (2.1) для вычисления NPV проекта. В рамках этого подхода также удается решить проблему множественности альтернативной характеристики проекта – величины IRR (internal rate of return), ставкой дисконта, при которой NPV проекта обращается в нуль.

Предложенный Д.Г. Кантором и С.А. Липманом подход к построению оценки доходности пула тиражируемых инвестиционных проектов в условиях несовершенного рынка капитала получил развитие в [414]. В работе Л. Адлера и Д. Гейла с использованием теоремы А. Пуанкаре была установлена связь условия оценки Кантора–Липман (см. [4, 5]) доходности пула тиражируемых инвестиционных проектов с теоремой об арбитраже. В [7, 8] И.М. Сонина и Э.Л. Пресмана переформулировали постановку задачи об оптимальной стратегии инвестирования в терминах динамики вектора, задающего финансовую позицию инвестора. Компоненты этого вектора определяют текущие и будущие запасы ликвидных средств инвестора при условии, что он не начинает новых инвестиционных проектов. В терминах динамики вектора финансовых позиций исследованы режимы сбалансированного роста, введено понятие ликвидных финансовых состояний и охарактеризованы режимы, являющиеся финансовыми пирамидами. Подход И.М. Сонина и Э.Л. Пресмана получил развитие в работе В.З. Беленького (см. [9]), в которой модель Кантора–Липмана интерпретировалась как обобщенная модель экономического роста Неймана–Гейла, а результаты ее исследования трактовались как анализ магистральных эффектов в соответствующей краевой задаче.

В [1214] задача об оценке финансового состояния инвестора формализована в виде уравнения Беллмана, решение которого позволяет строить оценки финансовых состояний инвестора и оценивать их ликвидность на основе принадлежности к конусу ликвидных состояний. В [12] показано, что невырожденное решение уравнения Беллмана в конусе ликвидных состояний существует только в том случае, если в качестве дефлятора денежных потоков выбран минимальный положительный корень инвестиционной функции (построенной как поточечный максимум преобразований Лапласа денежных потоков всего пула инвестиционных проектов, доступных инвестору). Классический вариант модели Кантора–Липмана сформулирован для дискретного времени, когда моменты времени поступления денежных потоков синхронизированы. В [11] подход Кантора–Липмана обобщен на случай непрерывного времени, для которого доказан аналог теоремы 1. Это позволяет расширить возможности применения модели Кантора–Липмана в моделях экономической динамики. В [8] исследована более общая постановка задачи инвестора, в которой авторы отказываются от условия (2.7). Доказано, что в этом случае инвестиционная стратегия с доходностью, превышающей минимальный корень инвестиционной функции, требует постоянных реинвестиций для поддержания роста, т.е. из такого инвестиционного процесса нельзя выйти и фиксировать доход (ситуация “финансовой пирамиды”).

Результаты исследования модели и ее модификаций показывают, что подход Кантора–Липмана позволяет

• формализовать описание предпринимательской среды на несовершенном рынке капитала;

• однозначно решить вопрос о выборе дефлятора финансовых потоков $r$ на несовершенном рынке капитала при расчете показателя $NPV\left( {{\mathbf{a}},r} \right)$ инвестиционного проекта a, определенного в (2.1); тем самым дать однозначный ответ о рентабельности проекта a;

• вычислить значение дефлятора финансовых потоков $r$, используемого в выражении (2.1) для расчета показателя $NPV\left( {{\mathbf{a}},r} \right)$, как минимальный положительный корень инвестиционной функции $F\left( r \right)$, определенной в (2.8).

• корректно ввести понятие состояния предпринимательской среды на несовершенном рынке капитала, определяя его как значение дефлятора финансовых потоков $r$, вычисленного в терминах модели.

3. МОДЕЛИРОВАНИЕ ИНВЕСТИЦИОННОГО ПОВЕДЕНИЯ СОБСТВЕННИКА ПРОИЗВОДСТВА НА НЕСОВЕРШЕННОМ РЫНКЕ КАПИТАЛА

Предположим, что собственник капитала принимает инвестиционные решения в производственном секторе в условиях сформировавшейся предпринимательской среды на несовершенном рынке капитала, которую мы будем описывать с помощью модели Кантора–Липмана. Тогда характеристика состояния предпринимательской среды $r$ (доходность пула доступных инвестиционных проектов), определяемая как минимальный положительный корень инвестиционной функции $F\left( r \right)$, является экзогенным входным параметром при описании инвестиционного поведения собственника.

Замечание 1. В условиях существенного расхождения процентных ставок по депозитам и кредитам доходность предпринимательской среды $r$, вычисленная в терминах модели Кантора–Липмана, будет ниже процента по кредитам, т.е. $r < {{r}_{{\text{L}}}}$. Поэтому на несовершенном рынке капитала собственнику невыгодно поддерживать инвестиционную деятельность за счет кредитных средств, а источником рыночных инвестиций в производство является, преимущественно, чистая прибыль собственника от производственной деятельности.

Предположим, что инвестиционное поведение собственника производства определяется выбором одной из двух альтернатив: рыночные инвестиции в производство с доходностью $r$ или расходование полученной от производственной деятельности чистой прибыли на собственные нужды, т.е. на потребление. Целью нашего исследования является анализ гипотезы о том, что состояние предпринимательской среды в условиях несовершенного рынка капитала может влиять на решения собственника о способе вложения свободных средств. Поэтому в математической модели должна быть учтена существенная зависимость между поведением собственника производства и характеристикой состояния предпринимательской среды $r$.

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

Рассмотрим финансовую деятельность собственника производства, который в каждый момент времени $t \geqslant 0~$ получает мгновенную чистую прибыль $\pi $ от производственной деятельности (после уплаты всех налогов и обязательных платежей). Предположим, что собственник может расходовать прибыль либо на конечное потребление в объеме ${\text{С}}\left( t \right) \geqslant 0,$ либо на инвестиции в производство в предпринимательской среде, характеризующейся доходностью $r.$

Предположим, что в начальный момент времени собственник обладает положительным начальным капиталом ${{M}_{0}} > 0.$ Тогда с учетом сделанных предположений динамика капитала $M\left( t \right)$ собственника производства определяется следующим уравнением:

(3.1)
$\frac{{dM\left( t \right)}}{{dt}} = \pi - C\left( t \right) + rM\left( t \right),\quad M\left( 0 \right) = {{M}_{0}} > 0,$
где $rM\left( t \right)$ – доход от инвестиций в производство в предпринимательской среде с показателем доходности $r.$

Напомним, что заемные средства не могут обеспечивать инвестиционную деятельность на несовершенном рынке капитала (см. замечание 1), поэтому должно выполняться условие неотрицательности капитала собственника

(3.2)
$M\left( t \right) \geqslant 0.$

Пример. Следуя традиционной схеме Рамсея, предположим, что целью собственника является максимизация дисконтированного с коэффициентом $\delta > 0$ потребления $C\left( t \right)$ за конечный период времени $\left[ {0,T} \right]$. Тогда с учетом (3.1), (3.2) получим, что выбор стратегии $С\left( t \right) \geqslant 0$ потребительского поведения собственника определяется как решение следующей задачи оптимального управления с фазовым ограничением:

(3.3)
$\int\limits_0^T {C\left( t \right){{e}^{{ - \delta t}}}dt} \to \mathop {\max }\limits_{} ,$
(3.4)
$\frac{{dM\left( t \right)}}{{dt}} = \pi - C\left( t \right) + rM\left( t \right),$
(3.5)
$M\left( 0 \right) = {{M}_{0}} > 0,$
(3.6)
$M\left( t \right) \geqslant 0,$
(3.7)
$C\left( t \right) \geqslant 0.$

Предложение 1. Оптимальное поведение собственника в задаче (3.3)–(3.7) в зависимости от соотношения параметров модели определяется одним из двух режимов.

Режим 1. Если $r < \delta $ (низкая доходность инвестиций), то оптимальная стратегия собственника заключается в потреблении всего начального запаса капитала ${{M}_{0}}$ в нулевой момент времени и дальнейшее поддержание потребления на максимальном уровне за счет всей прибыли от производственной деятельности в условиях полного отсутствия инвестиций, т.е.

$C\left( 0 \right) = \delta \left( t \right){{M}_{0}},\quad C\left( t \right) \equiv \pi ,\quad M\left( t \right) \equiv 0~\quad {\text{при\;}}\;{\text{всех}}\quad t \in \left( {0,T} \right],$
где $\delta \left( t \right)$ – дельта-функция.

Режим 2. Если $r > \delta $ (высокая доходность инвестиций), то собственник откладывает потребление до конца временного горизонта планирования своей деятельности, инвестирует всю прибыль в производство, а в конце периода потребляет все накопленные средства, т.е.

$C\left( t \right) \equiv 0~\quad {\text{при}}\;{\text{\;всех}}\quad ~t \in \left[ {0,T} \right)~,\quad C\left( T \right) = \delta \left( {t - T} \right)M\left( T \right),$
где $\delta \left( {t - T} \right)$ – дельта-функция, сосредоточенная в точке $T$.

Доказательство. Рассмотрим значение функционала (3.3) для режимов 1 и 2 и покажем, что для любой другой допустимой траектории $C\left( t \right)$ значение функционала не может увеличиться.

1. Пусть $r < \delta $. Обозначим функционал в задаче (левая часть (3.3)) через $\Phi \left( T \right)$ и найдем его значение ${{\Phi }_{1}}\left( T \right)$ в режиме 1:

${{\Phi }_{1}}\left( T \right) = \int\limits_0^T {C\left( t \right){{e}^{{ - \delta t}}}dt} = {{M}_{0}} + ~\frac{\pi }{\delta }\left( {1 - {{e}^{{ - \delta T}}}} \right).$
Предположим, что существует траектория $C\left( t \right)$, отличная от режима 1, в которой значение функционала $\Phi \left( T \right)$ будет больше, чем ${{\Phi }_{1}}\left( T \right)$. Из (3.4) имеем
(3.8)
$C\left( t \right) = \pi - \frac{{dM\left( t \right)}}{{dt}} + rM\left( t \right).{\text{\;}}$
Подставляя (3.8) в функционал (3.3) и интегрируя по частям, получим
$\begin{gathered} \Phi \left( T \right) = \int\limits_0^T {\left( {\pi - \frac{{dM\left( t \right)}}{{dt}} + rM\left( t \right)} \right){{e}^{{ - \delta t}}}dt} = \\ = {{M}_{0}} + \frac{\pi }{\delta }\left( {1 - {{e}^{{ - \delta T}}}} \right) - ~M\left( T \right){{e}^{{ - \delta T}}} + \int\limits_0^T {\left( {r - \delta } \right)M\left( t \right){{e}^{{ - \delta t}}}dt} \leqslant {{M}_{0}} + ~\frac{\pi }{\delta }\left( {1 - {{e}^{{ - \delta T}}}} \right) = {{\Phi }_{1}}\left( T \right). \\ \end{gathered} $
Последнее неравенство очевидно следует из соотношения $r - \delta < 0$ и условия (3.6). Противоречие.

2. Пусть $r > \delta $. Обозначим значение функционалa $\Phi \left( T \right)$ в режиме 2 через ${{\Phi }_{2}}\left( T \right)$. В режиме 2 уравнение (3.4) принимает вид

$\frac{{dM\left( t \right)}}{{dt}} = \pi + rM\left( t \right),\quad t \in \left[ {0,T} \right).$
Решение последнего уравнения имеет вид
(3.9)
$M\left( t \right) = - ~\frac{\pi }{r} + \left( {{{M}_{0}} + \frac{\pi }{r}} \right){{e}^{{rt}}},$
тогда
${{\Phi }_{2}}\left( T \right) = - ~\frac{\pi }{r}{{e}^{{ - \delta T}}} + \left( {{{M}_{0}} + \frac{\pi }{r}} \right){{e}^{{\left( {r - \delta } \right)T}}}.$
Предположим, что существует траектория $C\left( t \right)$, отличная от режима 2, в которой значение функционала $\Phi \left( T \right)$ будет больше, чем ${{\Phi }_{2}}\left( T \right)$. Аналогично случаю 1, имеем
$\Phi \left( T \right) = \int\limits_0^T {\left( {\pi - \frac{{dM\left( t \right)}}{{dt}} + rM\left( t \right)} \right){{e}^{{ - \delta t}}}dt} = {{M}_{0}} + \frac{\pi }{\delta }\left( {1 - {{e}^{{ - \delta T}}}} \right) - ~M\left( T \right){{e}^{{ - \delta T}}} + \int\limits_0^T {\left( {r - \delta } \right)M\left( t \right){{e}^{{ - \delta t}}}} dt.$
Учитывая (3.9), получим
$\begin{gathered} \int\limits_0^T {\left( {r - \delta } \right)} M\left( t \right){{e}^{{ - \delta t}}}dt = \int\limits_0^T {\left( {r - \delta } \right)} \left( { - ~\frac{\pi }{r} + \left( {{{M}_{0}} + \frac{\pi }{r}} \right){{e}^{{rt}}}} \right){{e}^{{ - \delta t}}}dt = \\ = \frac{{\left( {\delta - r} \right)}}{{\delta r}}\pi \left( {1 - {{e}^{{ - \delta T}}}} \right) + \left( {{{M}_{0}} + \frac{\pi }{r}} \right)\left( {{{e}^{{\left( {{\text{r}} - \delta } \right)T}}} - 1} \right). \\ \end{gathered} $
С учетом (3.6) легко показать, что
$\Phi \left( T \right) = - \frac{\pi }{r}{{e}^{{ - \delta T}}} + \left( {{{M}_{0}} + \frac{\pi }{r}} \right){{e}^{{\left( {r - \delta } \right)T}}} - M\left( T \right){{e}^{{ - \delta T}}} = {{\Phi }_{2}}\left( T \right) - M\left( T \right){{e}^{{ - \delta T}}} \leqslant {{\Phi }_{2}}\left( T \right).$
Получили противоречие. Тем самым режимы 1 и 2 являются оптимальным решением в задаче (3.3)–(3.7). Предположение 1 доказано.

Решение задачи (3.3)–(3.7) показывает, что оптимальные режимы потребления и инвестиций в рассмотренном варианте модели являются тривиальными, нет оптимального режима с одновременным потреблением и инвестициями. При этом поведение системы однотипное при любых значениях начальных условий. При стремлении горизонта планирования к бесконечности интеграл (3.3) расходится, если $r > \delta $, т.е. решение задачи не существует. Рассмотрение более общих вариантов моделей потребительского поведения рамсеевского типа позволяет подробно описывать динамику потребительского поведения домашних хозяйств и их поведения на рынке депозитов и кредитования (см., например, [17], [18]), однако, аналогично рассмотренной в примере модели, не позволяют дать содержательный ответ на вопрос о зависимости инвестиционных решений собственника от доходности предпринимательской среды. Таким образом, в рамках моделей рамсеевского типа не удается объяснить переключение из режима инвестиционного роста в режим стагнации в производственной сфере.

В настоящей работе мы предлагаем альтернативный подход, в котором целью собственника является развитие сектора потребления, обладающего своими потребительскими мощностями (недвижимость, предметы роскоши т.п.). Предположим, что собственник заинтересован в максимизации дисконтированного уровня средств потребления $k\left( t \right) \geqslant 0$ в потребительском секторе. Исходя из уровня накопленных к текущему моменту средств потребления $k\left( t \right)$, собственник планирует расходы на приобретение новых средств потребления в объеме $\zeta k\left( {\text{t}} \right)$, где $\zeta \geqslant 0$ – параметр, характеризующий максимальный темп роста средств потребления. Однако в зависимости от состояния конъюнктуры реализация планов может происходить не в полном объеме. Таким образом, в каждый момент времени $t$ собственнику доступны два способа вложения свободных средств:

• реализация доли $v\left( t \right)$ новых потребительских планов $\zeta k\left( {\text{t}} \right)$, где $v\left( t \right) \in \left[ {0,1} \right]$ – параметр управления;

• рыночные инвестиции в производство в предпринимательской среде с доходностью $r$.

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

$\frac{{dk}}{{dt}} = v\zeta k,$
где $k\left( 0 \right) = {{k}_{0}} > 0$, а изменение капитала собственника определяется уравнением
$\frac{{dM}}{{dt}} = \pi - v\zeta k + rM,$
с начальным условием $M\left( 0 \right) = {{M}_{0}} > 0.$

Цель собственника – максимизация дисконтированных с коэффициентом $\delta > 0$ средств потребления на бесконечном горизонте:

(3.10)
${\text{\;}}\int\limits_0^{ + \infty } {k\left( t \right){{e}^{{ - \delta t}}}dt} \to {\text{max}}{\text{.}}$

Замечание 2. 1) Выбор бесконечного горизонта интегрирования в целевом функционале (3.10) имеет естественную интерпретацию. Предположим, что горизонт планов собственника $T > 0$ конечен, но заранее собственнику неизвестен и является случайной величиной, распределенной по экспоненциальному закону с некоторым показателем $\gamma > 0$. В условиях неопределенности окончания деятельности естественно рассматривать задачу о максимизации математического ожидания дисконтированных средств потребления собственника, функционал в которой будет равен

$\int\limits_0^{ + \infty } {\gamma {{e}^{{ - \gamma T}}}} \int\limits_0^T {k\left( t \right){{e}^{{ - \delta t}}}dtdT} = \int\limits_0^{ + \infty } {k\left( t \right){{e}^{{ - \delta t}}}} \int\limits_t^{ + \infty } {\gamma {{e}^{{ - \gamma T}}}dTdt} = \int\limits_0^{ + \infty } {k\left( t \right){{e}^{{ - \left( {\delta + \gamma } \right)t}}}dt.} $
Таким образом, задача по-прежнему формализуется в виде (3.10), где в качестве нового коэффициента дисконтирования следует рассматривать величину $\gamma + \delta $.

2) Будем предполагать, что собственник строит планы новых потребительских проектов $\zeta k\left( {\text{t}} \right)$, исходя из предположения об их “реальности”. Формально это означает, что выполняется условие $\zeta < \delta $, гарантирующее сходимость интеграла (3.10).

В рамках сделанных предположений модель инвестиционного поведения собственника производства формализуется в виде следующей задачи:

(3.11)
$\begin{gathered} \int\limits_0^{ + \infty } {k\left( t \right){{e}^{{ - \delta t}}}dt} \to \mathop {\max }\limits_v {\text{,}}\quad \frac{{dk}}{{dt}} = v\zeta k,\quad k\left( 0 \right) = {{k}_{0}} > 0, \\ \frac{{dM}}{{dt}} = \pi - v\zeta k + rM,\quad M\left( 0 \right) = {{M}_{0}} > 0~,\quad M\left( t \right) \geqslant 0,\quad {\text{\;}}0 \leqslant v \leqslant 1. \\ \end{gathered} $

Задача (4.2), называемая далее ${{P}_{\infty }}$, представляет собой задачу оптимального управления на бесконечном горизонте с ограничением на фазовую переменную $M\left( t \right) \geqslant 0$, которое отражает условие несовершенства рынка капитала в рассматриваемой модели (см. замечание 1). Такая постановка приводит к техническим трудностям в ходе исследования модели. Исследование задачи оптимального управления ${{P}_{\infty }}$ приведено авторами в [16]. В следующем разделе настоящей работы мы приведем краткую схему исследования и сформулируем основные результаты, доказательство которых можно найти в [16].

4. СХЕМА ИССЛЕДОВАНИЯ МОДЕЛИ ИНВЕСТИЦИОННОГО ПОВЕДЕНИЯ СОБСТВЕННИКА. ОСНОВНЫЕ РЕЗУЛЬТАТЫ

Рассмотрим сначала вспомогательную задачу оптимального управления на конечном интервале времени $T$:

(4.1)
$\begin{gathered} \int\limits_0^T {k\left( t \right){{e}^{{ - \delta t}}}dt} \to \mathop {\max }\limits_v {\text{,}}\quad \frac{{dk}}{{dt}} = v\zeta k,\quad \frac{{dM}}{{dt}} = \pi - v\zeta k + rM, \\ k\left( 0 \right) = {{k}_{0}} > 0{\text{,}}\quad M\left( 0 \right) = {{M}_{0}} > 0,\quad M\left( t \right) \geqslant 0,\quad 0 \leqslant v \leqslant 1. \\ \end{gathered} $

Решение задачи (4.1), называемой далее ${{P}_{T}}$, проанализировано в [16]. Схема исследования сводится к применению принципа максимума Понтрягина для задач оптимального управления с фазовым ограничением, приведенного, например, в [19, c. 185], и дальнейшему изучению поведения траекторий прямой и сопряженной систем задачи ${{P}_{T}}$ в зависимости от значений параметров задачи и начальных условий. Приведем формулировки результатов.

Утверждение 1 (см. [16]). Рассмотрим случай $r < \delta $. Справедливы следующие утверждения.

1) В области значений начальных условий $\left( {{{M}_{0}},{{k}_{0}}} \right)$, определяемой соотношением

(4.2)
$0 < {{M}_{0}} < {{k}_{0}}\frac{\zeta }{{\zeta - r}}\left( {{{e}^{{\left( {\zeta - r} \right)T}}} - 1} \right) + \frac{\pi }{r}\left( {{{e}^{{ - rT}}} - 1} \right),$
уравнение
${{M}_{0}}{{e}^{{r\hat {T}}}} + \frac{\pi }{r}\left( {{{e}^{{r\hat {T}}}} - 1} \right) - \frac{{\zeta {{k}_{0}}}}{{\zeta - r}}\left( {{{e}^{{\zeta \hat {T}}}} - {{e}^{{r\hat {T}}}}} \right) = 0$
имеет единственный корень $\hat {T} \in \left. {[0,T} \right)$, и решение задачи ${{P}_{T}}$ определяется следующими соотношениями:

$v\left( t \right) = \left\{ \begin{gathered} 1,\quad 0 \leqslant t < \widehat {T,} \hfill \\ \frac{\pi }{{\zeta k\left( t \right)}},\quad \hat {T} < t \leqslant T; \hfill \\ \end{gathered} \right.{\text{\;}}$
$k\left( t \right) = \left\{ \begin{gathered} {{k}_{0}}{{e}^{{\zeta t}}},\quad 0 \leqslant t < \widehat {T,} \hfill \\ {{k}_{0}}{{e}^{{\zeta \hat {T}}}} + \pi \left( {t - \hat {T}} \right),\quad \hat {T} < t \leqslant T; \hfill \\ \end{gathered} \right.$
$M\left( t \right) = \left\{ \begin{gathered} {{M}_{0}}{{e}^{{rt}}} + \frac{\pi }{r}\left( {{{e}^{{rt}}} - 1} \right) - \frac{{\zeta {{k}_{0}}}}{{\zeta - r}}\left( {{{e}^{{\zeta t}}} - {{e}^{{rt}}}} \right),\quad 0 \leqslant t < \widehat {T,} \hfill \\ 0,\quad \hat {T} < t \leqslant T. \hfill \\ \end{gathered} \right.$

2) В области значений начальных условий $\left( {{{M}_{0}},{{k}_{0}}} \right)$ , удовлетворяющих неравенству

(4.3)
${{M}_{0}} \geqslant {{k}_{0}}\frac{\zeta }{{\zeta - r}}\left( {{{e}^{{\left( {\zeta - r} \right)T}}} - 1} \right) + \frac{\pi }{r}\left( {{{e}^{{ - rT}}} - 1} \right),$
решение задачи ${{P}_{T}}$ при любых $t \in \left[ {0,T} \right]$ определяется соотношениями

$v\left( t \right) = 1,$
$k\left( t \right) = {{k}_{0}}{{e}^{{\zeta t}}},$
$M\left( t \right) = {{M}_{0}}{{e}^{{rt}}} + \frac{\pi }{r}\left( {{{e}^{{rt}}} - 1} \right) - \frac{{\zeta {{k}_{0}}}}{{\zeta - r}}\left( {{{e}^{{\zeta t}}} - {{e}^{{rt}}}} \right).$

Утверждение 2. Рассмотрим случай $r > \delta $. Справедливы следующие утверждения.

1) В области значений начальных условий $\left( {{{M}_{0}},{{k}_{0}}} \right)$ , определяемой неравенствами (4.2), т.е.

$0 < {{M}_{0}} < {{k}_{0}}\frac{\zeta }{{\zeta - r}}\left( {{{e}^{{\left( {\zeta - r} \right)T}}} - 1} \right) + \frac{\pi }{r}\left( {{{e}^{{ - rT}}} - 1} \right),$
существует единственное решение $\left( {\tilde {T},\hat {T}} \right){\text{\;}}$ системы уравнений
$\left\{ \begin{gathered} {{M}_{0}}{{e}^{{r\hat {T}}}} + \frac{\pi }{r}\left( {{{e}^{{r\hat {T}}}} - 1} \right) - \frac{{\zeta {{k}_{0}}}}{{\zeta - r}}\left( {{{e}^{{\zeta \left( {\hat {T} - \tilde {T}} \right)}}} - {{e}^{{r\left( {\hat {T} - \tilde {T}} \right)}}}} \right) = 0, \hfill \\ \frac{r}{{\delta \left( {r - \zeta } \right)}}\left( {1 - {{e}^{{r\left( {\hat {T} - T} \right)}}}} \right)\left( {{{e}^{{\left( {\delta - \zeta } \right)\left( {\tilde {T} - \hat {T}} \right)}}} - {{e}^{{\left( {\delta - r} \right)\left( {\tilde {T} - \hat {T}} \right)}}}} \right) + \frac{1}{{\delta - \zeta }}\left( {1 - {{e}^{{\left( {\delta - \zeta } \right)\left( {\tilde {T} - \hat {T}} \right)}}}} \right) = 0, \hfill \\ \end{gathered} \right.$
причем $0 \leqslant \tilde {T} < \hat {T} \leqslant T$.

Решение задачи ${{P}_{T}}$ определяется соотношениями

${\text{\;}}v\left( t \right) = \left\{ \begin{gathered} \mathop 0\limits_{}^{} ,\quad 0 \leqslant t \leqslant \tilde {T}, \hfill \\ \mathop 1\limits_{}^{} ,\quad \tilde {T} \leqslant t \leqslant \hat {T}, \hfill \\ \frac{\pi }{{\zeta \left( {{{k}_{0}}{{e}^{{\zeta \left( {\hat {T} - \tilde {T}} \right)}}} + \pi \left( {t - \hat {T}} \right)} \right)}},\quad \hat {T} \leqslant t \leqslant T; \hfill \\ \end{gathered} \right.$
$k\left( t \right) = \left\{ \begin{gathered} {{k}_{0}},\quad 0 \leqslant t \leqslant \widetilde {T,} \hfill \\ {{k}_{0}}{{e}^{{\zeta \left( {t - \tilde {T}} \right)}}},\quad \tilde {T} \leqslant t \leqslant \widehat {T,} \hfill \\ {{k}_{0}}{{e}^{{\zeta \left( {\hat {T} - \tilde {T}} \right)}}} + \pi \left( {t - \hat {T}} \right),\quad \hat {T} \leqslant t \leqslant T; \hfill \\ \end{gathered} \right.$
$M\left( t \right) = \left\{ \begin{gathered} {{M}_{0}}{{e}^{{rt}}} + \frac{\pi }{r}\left( {{{e}^{{rt}}} - 1} \right),\quad 0 \leqslant t \leqslant \widetilde {T,} \hfill \\ {{M}_{0}}{{e}^{{rt}}} + \frac{\pi }{r}\left( {{{e}^{{rt}}} - 1} \right) - \frac{{\zeta {{k}_{0}}}}{{\zeta - r}}\left( {{{e}^{{\zeta \left( {t - \tilde {T}} \right)}}} - {{e}^{{r\left( {t - \tilde {T}} \right)}}}} \right),\quad \tilde {T} < t \leqslant \widehat {T,} \hfill \\ 0,\quad \hat {T} < t \leqslant T. \hfill \\ \end{gathered} \right.$

2) В области значений начальных условий $\left( {{{M}_{0}},{{k}_{0}}} \right)$, удовлетворяющих неравенству (4.3), т.е.

${{M}_{0}} \geqslant {{k}_{0}}\frac{\zeta }{{\zeta - r}}\left( {{{e}^{{\left( {\zeta - r} \right)T}}} - 1} \right) + \frac{\pi }{r}\left( {{{e}^{{ - rT}}} - 1} \right),$
решение задачи ${{P}_{T}}{\text{\;при\;любых\;}}t \in \left[ {0,T} \right]{\text{\;}}$ определяется соотношениями

$v\left( t \right) = 1,$
$k\left( t \right) = {{k}_{0}}{{e}^{{\zeta t}}},$
$M\left( t \right) = {{M}_{0}}{{e}^{{rt}}} + \frac{\pi }{r}\left( {{{e}^{{rt}}} - 1} \right) - \frac{{\zeta {{k}_{0}}}}{{\zeta - r}}\left( {{{e}^{{\zeta t}}} - {{e}^{{rt}}}} \right).$

Замечание 3. 1) Утверждения 1, 2 позволяют разделить фазовую плоскость ${\text{(}}{{M}_{0}}{\text{,\;}}{{k}_{0}}{\text{)}}$ на две области с разными режимами оптимального управления для задачи ${{P}_{T}}.$ Граница между областями определяется равенством

(4.4)
${{M}_{0}} = {{k}_{0}}\frac{\zeta }{{\zeta - r}}\left( {{{e}^{{\left( {\zeta - r} \right)T}}} - 1} \right) + \frac{\pi }{r}\left( {{{e}^{{ - rT}}} - 1} \right).$

2) При $T \to \infty $ линейная граница (4.4) на плоскости $({{M}_{0}},{{k}_{0}})$ стремится к прямой:

(4.5)
${{k}_{0}} = \left( {r{{M}_{0}} + \pi } \right)\frac{{{{{\left( {r - \zeta } \right)}}_{ + }}}}{{r\zeta }}.$

Вернемся теперь к исходной задаче ${{P}_{\infty }}$ с бесконечным горизонтом. Необходимое и достаточное условие оптимальности в задаче ${{P}_{\infty }}$ – существование функции цены, удовлетворяющей уравнению Гамильтона–Якоби–Беллмана для задачи ${{P}_{\infty }}$, в котором закон управления, соответствующий оптимальной траектории, является максимизатором.

Задаче оптимального управления ${{P}_{\infty }}$ соответствует уравнение Гамильтона–Якоби–Беллмана вида (см. [20, гл. 1, ч. 3])

(4.6)
$\delta V = k + \zeta k{{\left( {\frac{{\partial V}}{{\partial k}} - \frac{{\partial V}}{{\partial M}}} \right)}_{ + }} + \frac{{\partial V}}{{\partial M}}\left( {\pi + rM} \right).$

Анализ решения задачи ${{P}_{T}}$ на конечном горизонте позволяет найти вязкостное (обобщенное) (см. [20]) решение уравнения Гамильтона–Якоби–Беллмана для задачи ${{P}_{\infty }}$ и построить позиционное оптимальное управление в исходной модели. В [16] доказаны следующие основные результаты.

Теорема 2 (см. [16]). Пусть $r < $ $\delta .$ Тогда для задачи ${{P}_{\infty }}$ справедливы следующие утверждения.

1) Оптимальное управление определяется выражением

$v\left( {k,M} \right) = \left\{ \begin{gathered} 1,\quad M > 0, \hfill \\ {\text{min}}\left( {\frac{\pi }{{\zeta k}},1} \right),\quad M = 0. \hfill \\ \end{gathered} \right.$

2) Функция цены

${{V}_{1}}\left( {k,M{\text{\;}}} \right) = \left\{ \begin{gathered} \frac{k}{{\delta - \zeta }} - \frac{{\zeta k}}{{\delta \left( {\delta - \zeta } \right)}}{{e}^{{ - \left( {\delta - \zeta } \right)\hat {T}\left( {M,k} \right)}}} + \frac{\pi }{{{{\delta }^{2}}}}{{e}^{{ - \delta \hat {T}\left( {M,k} \right)}}},\quad {\text{если}}\quad rM + \pi < \frac{{r\zeta k}}{{r - \zeta }},{\text{\;\;}} \hfill \\ \frac{k}{{\delta - \zeta }},\quad {\text{если}}\quad rM + \pi \geqslant \frac{{r\zeta k}}{{r - \zeta }}, \hfill \\ \end{gathered} \right.$
является дифференцируемой функцией и удовлетворяет уравнению Гамильтона–Якоби–Беллмана (4.6). Величина $\hat {T}\left( {M,k} \right) > 0$ является единственным корнем уравнения
$\left( {rM + \pi - \frac{{r\zeta k}}{{r - \zeta }}} \right){{e}^{{r\hat {T}}}} = \pi - \frac{{r\zeta k}}{{r - \zeta }}{{e}^{{\zeta \hat {T}}}}$
при $rM + \pi < \frac{{r\zeta k}}{{r - \zeta }}$, $r - \zeta > 0$. В противном случае положим $\hat {T}\left( {M,k} \right) = + \infty $.

3) Фазовые траектории системы, соответствующие оптимальному управлению, имеют вид

$k\left( t \right) = \left\{ \begin{gathered} {{k}_{0}}{{e}^{{\zeta t}}},\quad 0 \leqslant t < \hat {T}\left( {{{M}_{0}},{{k}_{0}}} \right), \hfill \\ {{k}_{0}}{{e}^{{\zeta \hat {T}\left( {{{M}_{0}},{{k}_{0}}} \right)}}} + \pi \left( {t - \hat {T}\left( {{{M}_{0}},{{k}_{0}}} \right)} \right),\quad t \geqslant \hat {T}\left( {{{M}_{0}},{{k}_{0}}} \right); \hfill \\ \end{gathered} \right.$
$M\left( t \right) = \left\{ \begin{gathered} {{M}_{0}}{{e}^{{rt}}} + \frac{\pi }{r}\left( {{{e}^{{rt}}} - 1} \right) - \frac{{\zeta {{k}_{0}}}}{{\zeta - r}}\left( {{{e}^{{\zeta t}}} - {{e}^{{rt}}}} \right),\quad 0 \leqslant t < \hat {T}\left( {{{M}_{0}},{{k}_{0}}} \right), \hfill \\ 0,\quad t \geqslant \hat {T}\left( {{{M}_{0}},{{k}_{0}}} \right). \hfill \\ \end{gathered} \right.$

Теорема 3 (см. [16]). Пусть $r > $ $\delta .$ Тогда для задачи ${{P}_{\infty }}$ справедливы следующие утверждения.

1) Оптимальное управление определяется выражением

$v\left( {k,M~} \right) = \left\{ \begin{gathered} 0,\quad {\text{если}}\quad rM + \pi < \frac{{r\zeta k}}{{r - \zeta }}, \hfill \\ 1,\quad {\text{если}}\quad rM + \pi \geqslant \frac{{r\zeta k}}{{r - \zeta }}. \hfill \\ \end{gathered} \right.$

2) Функция цены

$V\left( {k,M} \right) = \left\{ \begin{gathered} \frac{k}{\delta }\left( {1 + \frac{\zeta }{{\left( {\delta - \zeta } \right)}}{{{\left( {\frac{{\left( {r - \zeta } \right)\left( {rM + \pi } \right)}}{{r\zeta k}}} \right)}}^{{\frac{\delta }{r}}}}} \right),\quad {\text{если\;}}\quad rM + \pi \leqslant \frac{{r\zeta k}}{{r - \zeta }}, \hfill \\ \frac{k}{{\delta - \zeta }},\quad {\text{если}}\quad rM + \pi > \frac{{r\zeta k}}{{r - \zeta }} \hfill \\ \end{gathered} \right.$
является вязкостным решением уравнения Гамильтона–Якоби–Беллмана (4.6).

3) Фазовые траектории системы, соответствующие оптимальному управлению, имеют вид

$k\left( t \right) = \left\{ \begin{gathered} {{k}_{0}},\quad 0 \leqslant t \leqslant \widetilde {T,} \hfill \\ {{k}_{0}}{{e}^{{\zeta \left( {t - \tilde {T}} \right)}}},\quad t \geqslant \tilde {T}; \hfill \\ \end{gathered} \right.$
$M\left( t \right) = \left\{ \begin{gathered} {{M}_{0}}{{e}^{{rt}}} + \frac{\pi }{r}\left( {{{e}^{{rt}}} - 1} \right),\quad 0 \leqslant t \leqslant \widetilde {T,} \hfill \\ {{M}_{0}}{{e}^{{rt}}} + \frac{\pi }{r}\left( {{{e}^{{rt}}} - 1} \right) - \frac{{\zeta {{k}_{0}}}}{{\zeta - r}}\left( {{{e}^{{\zeta \left( {t - \tilde {T}} \right)}}} - {{e}^{{r\left( {t - \tilde {T}} \right)}}}} \right),\quad t \geqslant \tilde {T}, \hfill \\ \end{gathered} \right.$
где

(4.7)
$~\tilde {T} = {{\left( {\frac{1}{r}\ln \frac{{\zeta {{k}_{0}}}}{{\left( {r - \zeta } \right)\left( {r{{M}_{0}} + \pi } \right)}}} \right)}_{ + }}.$

Замечание 4. Пусть $r > $ $\delta $. Если $\tilde {T} > 0$ (см. (4.7)), то фазовые траектории задачи ${{P}_{\infty }}$, соответствующие оптимальному управлению, при $t \geqslant \tilde {T}$ имеют вид

$k\left( t \right) = {{k}_{0}}{{e}^{{\zeta \left( {t - \tilde {T}} \right)}}},\quad M\left( t \right) = \frac{{\zeta {{k}_{0}}}}{{r - \zeta }}{{e}^{{\zeta \left( {t - \tilde {T}} \right)}}} - \frac{\pi }{r}$
и связаны линейным соотношением

$rM\left( t \right) + \pi = \frac{{r\zeta k\left( t \right)}}{{r - \zeta }}.$

5. ЭКОНОМИЧЕСКАЯ ИНТЕРПРЕТАЦИЯ РЕЗУЛЬТАТОВ

Результаты исследования модели инвестиционного поведения собственника производства в предпринимательской среде с доходностью $r$ на несовершенном рынке капитала позволяют проанализировать зависимость оптимальных стратегий собственника от параметров модели и начальных условий.

1. В условиях низкой доходности инвестиций ($r < \delta $) оптимальной стратегией поведения собственника является приоритетное расходование прибыли от производственной деятельности на личное потребление (теорема 2). На фиг. 3 представлен типичный фазовый портрет траекторий системы в задаче $({{P}_{\infty }})$ при условиях $r > \delta ,$ $~r > \zeta $. Заметим, что на фазовой плоскости существуют две области значений фазовых переменных исходной задачи, в которых наблюдается разный тип поведения собственников. Граница между областями – прямая, определенная в (4.5) (при $r < \zeta $ прямая (4.5) вырождается в прямую $k\left( t \right) = 0$ и нижняя область исчезает). Траектории системы, стартующие в верхней области на фиг. 3, соответствуют собственникам, которые обладают меньшим начальным капиталом ${{M}_{0}}$, чем агенты, чьи траектории расположены в нижней области на фиг. 3. Поэтому естественно интерпретировать агентов из верхней области как собственников малого и среднего бизнеса, а агентов из нижней области – как крупных собственников. Анализ полученных в теореме 2 результатов позволяет утверждать, что для собственников малого и среднего бизнеса при низкой доходности предпринимательской среды характерно интенсивное накопление средств потребления до полного исчерпания капитала $M$, а далее расходование на приобретение средств потребления всего текущего дохода от производственной деятельности (верхняя область на фиг. 3).

Фиг. 3.

Фазовый портрет системы в случае $r < $ $\delta $.

Представители крупного бизнеса имеют другой тип поведения. Их финансовое состояние позволяет наращивать средства потребления с максимально возможным темпом $\zeta $ , а оставшуюся часть прибыли направлять на инвестиции в производство независимо от прочих соотношений параметров модели (нижняя область на фиг. 3).

2. В условиях высокой доходности инвестиций ($r > \delta $) представителям малого и среднего бизнеса оказывается выгодно отсрочить личное потребление и наращивать капитал за счет доходов от рыночных инвестиций (теорема 3). На фиг. 4 представлен типичный фазовый портрет траекторий системы в задаче ${{P}_{\infty }}$ при условии $r > \delta $ . Аналогично п. 1, на фазовой плоскости $\left( {M,k{\text{\;}}} \right)$ можно выделить две области: верхняя соответствует представителям малого и среднего бизнеса, а нижняя – представителям крупного бизнеса. Заметим, что поведение крупных собственников (нижняя область на фиг. 4) не отличается от рассмотренного в п. 1 в ситуации низкой доходности предпринимательской среды. Таким образом, на крупных инвесторов не оказывает влияние состояние предпринимательской среды и целесообразно использовать другие методы стимулирования инвестиционной активности. В отличие от этого представители малого и среднего бизнеса инвестируют в развитие производства до определенного уровня своего благосостояния и только потом начинают вкладывать прибыль от основной деятельности в накопление средств потребления и одновременно инвестировать в производство с одинаковым максимально возможным темпом $\zeta $ (верхняя область на фиг. 4).

Фиг. 4.

Фазовый портрет системы в случае $r > $ $\delta $.

6. ОБЪЯСНЕНИЕ ПЕРЕКЛЮЧЕНИЯ ИЗ РЕЖИМА ВОССТАНОВИТЕЛЬНОГО РОСТА В РЕЖИМ СТАГНАЦИИ НА ОСНОВЕ ПОДХОДА КАНТОРА–ЛИПМАНА

В результате исследования модели инвестиционного поведения собственника на несовершенном рынке капитала показано, что улучшение состояния предпринимательской среды (рост доходности $r$) может привести к переключению поведения собственников малого и среднего бизнеса из режима потребления в режим инвестирования. Наоборот, ухудшение характеристик инвестиционной среды (снижение доходности $r$) приводит к сворачиванию инвестиционной деятельности представителей малого и среднего бизнеса.

В заключение работы на основе подхода Кантора–Липмана опишем возможный механизм переключения из режима восстановительного роста в режим стагнации российской экономики на рубеже 2008 г. Напомним, что в терминах модели Кантора–Липмана величина доходности $r$ на несовершенном рынке капитала может быть найдена как минимальный подложительный корень инвестиционной функции $F\left( r \right)$, определенной в (2.8).

На фиг. 5 приведена типичная для рассматриваемого периода схема взаимодействия экономических агентов, поведение которых могло повлиять на прекращение роста в производственной сфере российской экономики.

Фиг. 5.

Схема взаимодействия агентов в российской экономике.

Период восстановительного роста сопровождался крупными высокодоходными (процент по доходу ${{r}_{0}} > 0$) национальными проектами по реформированию инфраструктурных отраслей (реформа электроэнергетики, связь, транспортные проекты). Государство выделяло средства на эти проекты при условии рыночного софинансирования. Источником такого софинансирования являлась система коммерческих банков. В ситуации высокой доходности ${{r}_{0}}$ национальных проектов банки могли назначать достаточно высокие проценты по кредитам $0 < {{r}_{{\text{L}}}} < {{r}_{0}}$. Партнером банков являлись собственники предприятий обрабатывающей промышленности, которые пользовались услугами кредитования под процент ${{r}_{{\text{L}}}}$ и депонирования средств под процент $0 < {{r}_{{\text{D}}}} < {{r}_{{\text{L}}}}$, который также назначался банками. При этом расхождение этих процентных ставок оказывалось значительным в условиях несовершенного рынка капитала, однако величина процента по депозитам ${{r}_{{\text{D}}}}$ могла поддерживаться банками на достаточно высоком уровне при высокой процентной ставке ${{r}_{0}}$. В рамках подхода Кантора–Липмана это означает, что инвестиционный полином (2.8), описывающий предпринимательскую среду собственника в производственной сфере, мог иметь несколько корней (далее мы рассматриваем ситуацию трех корней ${{r}_{1}},\;~{{r}_{2}},\;~{{r}_{3}}$, см. фиг. 6, 7), но высокий процент ${{r}_{D}}$ экранировал два из трех его положительных корней (фиг. 6), что приводило к удержанию относительно высокой ставки доходности ${{r}_{3}}$ предпринимательской среды в производственной сфере. Эта ситуация проиллюстрирована на фиг. 6.

Фиг. 6.

Ситуация высокой доходности предпринимательской среды: ${{r}_{2}} < {{r}_{{\text{D}}}} < {{r}_{{\text{L}}}} < {{r}_{0}}$.

Фиг. 7.

Ситуация низкой доходности предпринимательской среды: ${{r}_{{\text{D}}}} < {{r}_{2}} < {{r}_{{\text{L}}}} < {{r}_{0}}$.

Со временем рост транзакционных издержек в инфраструктурных отраслях привел к снижению доходности ${{r}_{0}}$ национальных проектов. Соответственно, банки снизили процент по кредитам ${{r}_{{\text{L}}}} < {{r}_{0}}$. Следствием этого явилось снижение процента по депозитам ${{r}_{{\text{D}}}}$, который не может превышать ${{r}_{{\text{L}}}}$. Снижение ${{r}_{{\text{D}}}}$ было более значительным, поскольку система коммерческих банков к этому моменту накопила большой объем ликвидности и не была более заинтересована в привлечении депозитов. В результате падения процента ${{r}_{D}}$ минимальный положительный корень инвестиционной функции в производственной сфере перескочил из значения ${{r}_{3}}$ в значение ${{r}_{1}} < {{r}_{3}}$. Данная ситуация проиллюстрирована на фиг. 7. В результате доходность предпринимательской среды в обрабатывающем секторе производства резко снизилась. Собственники предприятий могли либо продолжать инвестиции под невыгодные проценты ${{r}_{1}},$ либо положить средства на депозит под более высокий процент ${{r}_{{\text{D}}}} > {{r}_{1}}$ (режим рантье). Можно предположить, что коэффициент дисконтирования доходов собственника производства фактически приблизился к значению процента по депозиту $\delta \approx {{r}_{{\text{D}}}}$, поэтому проект депонирования стал альтернативой производственной деятельности. Такая ситуация могла привести к спаду инвестиционной активности в производственной сфере и способствовать переходу в режим стагнации российской экономики.

Список литературы

  1. Шананин А.А. Анализ финансового состояния инвестора на основе модели Кантора–Липмана // Тр. ИММ УрО РАН. 2020. Т. 26. № 1. С. 293–306.

  2. Петров А.А., Поспелов И.Г., Шананин А.А. Опыт математического моделирования экономики. М.: Энергоатомиздат, 1996.

  3. Петров А.А., Поспелов И.Г., Шананин А.А. От госплана к неэффективному рынку: математический анализ эволюции российских экономических структур. UK: The Edwin Mellen Press, 1999.

  4. Cantor D.G., Lipman S.A. Investment selection with imperfect capital markets // Econometrica. 1983. V. 51. № 4. P. 1121–1144.

  5. Cantor D.G., Lipman S.A. Optimal investment selection with a multitude of projects // Econometrica. 1995. V. 63. № 5. P. 1231–1240.

  6. Adler L., Gale D. Arbitrate and growth rate for riskless investments in a stationary economy // Mathematical F-inance. 1997. V. 7. № 1. P. 73–81.

  7. Sonin I.M. Growth rate, internal rates of return and turn pikes in an investment model // Economic Theory. 1996. V. 5. P. 383–400.

  8. Presman E.L., Sonin I.M. Growth rate, internal rates of return and financial bubbles. Moscow: CEMI Rus. Acad. Sci., 2000.

  9. Беленький В.З. Экономическая динамика: анализ инвестиционных проектов в рамках линейной модели Неймана–Гейла. М.: Рос. экон. шк., 2002.

  10. Ващенко М.П. Оценка доходности инвестиционных проектов в условиях неопределенности // Матем. моделирование. 2009. Т. 21. № 3. С. 18–30.

  11. Ващенко М.П., Шананин А.А. Оценка доходности пула инвестиционных проектов в модели оптимального инвестирования в непрерывном времени // Матем. моделирование. 2012. Т. 24. № 3. С. 70–86.

  12. Shananin A.A., Vashchenko M.P., Zhang Sh. Financial bubbles existence in the Cantor–Lippman model for continuous time // Lobachevskii J. Math. 2018. V. 39. № 7. P. 929–935.

  13. Шананин А.А. Математическое моделирование инвестиций на несовершенном рынке капитала // Тр. ИММ УрО РАН. 2019. Т. 25. № 4. С. 265–274.

  14. Шананин А.А. Анализ финансового состояния инвестора на основе модели Кантора–Липмана // Тр. ИММ УрО РАН. 2020. Т. 26. № 1. С. 293–306.

  15. Ramsey F.P. A mathematical theory of savings // Econ. J. 1928. № 38. P. 543–559.

  16. Obrosova N.K., Shananin A.A., Spiridonov A.A. A Model of investment behavior of enterprise owner in an imperfect capital market // Lobachevskii J. Math. 2022. V. 43. № 4. P. 1023–1036.

  17. Рудева А.В., Шананин А.А. Синтез управления в модифицированной модели Рамсея с учетом ограничения ликвидности. // Дифференц. уравнения. 2009. Т. 45. № 12. С. 1799–1803.

  18. Тарасенко М.В., Трусов Н.В., Шананин А.А. Математическое моделирование экономического положения домашних хозяйств в России // Ж. вычисл. матем. и матем. физ. 2021. Т. 61. № 6. С. 1034–1056.

  19. Кларк Ф. Оптимизация и негладкий анализ. М.: Наука, 1988.

  20. Bardi M. Capuzzo-Dolcetta I. Optimal control and viscosity solutions of Hamilton–Jacobi–Bellman equations. Boston: Birkhauser, 1997.

Дополнительные материалы отсутствуют.