Модель стохастической динамики производственного предприятия, учитывающая оптимальные нормы внутренних инвестиций
- Авторы: Ильина Е.А.1, Сараев Л.А.1
-
Учреждения:
- Самарский национальный исследовательский университет имени академика С.П. Королева
- Выпуск: Том 14, № 3 (2023)
- Страницы: 205-218
- Раздел: МАТЕМАТИЧЕСКИЕ И ИНСТРУМЕНТАЛЬНЫЕ МЕТОДЫ ЭКОНОМИКИ
- URL: https://journals.ssau.ru/eco/article/view/26941
- DOI: https://doi.org/10.18287/2542-0461-2023-14-3-205-218
- ID: 26941
Цитировать
Полный текст
Аннотация
В публикуемой статье предложена новая экономико-математическая модель стохастической динамики развития многофакторного производственного предприятия, восстановление ресурсов которого обеспечивается за счет внутренних инвестиций. Прогнозирование роста объемов выпуска продукции, производственных издержек, прибыли и амортизационных отчислений описывается системой стохастических дифференциальных уравнений. Показано, что динамике оптимального развития предприятия, в рамках которой оно выходит на производственные мощности, отвечающие максимальной прибыли, соответствуют эффективные коэффициенты норм внутренних инвестиций, для вычисления которых установлена специальная система уравнений. Показано, что выбор неэффективных коэффициентов норм внутренних инвестиций не дает возможность предприятию обеспечить свою максимальную прибыль.
Полный текст
Введение
Устойчивый экономический рост производственных предприятий и экономических систем формирует стабильное развитие всей национальной экономики. Прогнозирование экономико-математическими методами показателей динамики развития субъектов экономики является одной из актуальных проблем современной экономической теории, успешное решение которой позволяет проводить полноценную аналитику бизнес-процессов, вычислять оптимальные объемы выпуска продукции, ресурсов, издержек и прибыли, которые способствуют выходу предприятий на эффективные производственные мощности. Классические положения теории экономического роста предприятий и экономических систем подробно представлены в работах [1–7].
На их основе разработан целый спектр моделей роста экономических систем, учитывающий роль технических инноваций и информационных технологий [8–32].
Динамика развития предприятий, которая в силу объективных обстоятельств носит, как правило, стохастический характер, определяется взаимодействием объемов инвестиций в производство и амортизационных отчислений на восстановление объемов ресурсов и затрат на модернизацию средств производства. Поэтому одним из главных математических инструментов для построения моделей экономического развития предприятий является аппарат стохастических дифференциальных уравнений и их систем [33–39].
Целью публикуемой работы является разработка новой экономико-математической модели динамики развития предприятия, которая учитывает случайный характер влияния сопровождающих издержек производства, амортизационных отчислений и внутренних инвестиций. Такая модель позволяет прогнозировать выход мощностей предприятия на эффективное предельное состояние производства, при котором прибыль предприятия становится максимальной.
Пусть предприятие для выпуска своей продукции привлекает объемы факторов производства , представляющие собой капитал, трудовые ресурсы, материалы, технологии, инновации и т. д.
Величины предполагаются непрерывными и непрерывно дифференцируемыми функциями времени.
Начальные значения функций считаются известными, их предельные значения определяется динамикой развития предприятия и подлежат вычислению.
Объем выручки предприятия обеспечивается многофакторной производственной функцией Кобба – Дугласа
. (1)
Здесь показатели степени представляют собой эластичности выпуска по соответствующим ресурсам , коэффициент равен объему выручки, полученной от переработки единичных объемов ресурсов.
Пропорциональные издержки предприятия записываются в виде
. (2)
Здесь – стоимости затрат на единичные объемы ресурсов, – постоянные затраты предприятия.
Формула прибыли предприятия имеет вид
. (3)
Для вычисления максимальной прибыли предприятия необходимо приравнять к нулю все частные производные функции прибыли (3):
,
и составить систему уравнений
. (4)
Решениями системы уравнений (4) являются значения ресурсов
, (5)
соответствующих максимальному значению прибыли :
. (6)
Предельные значения ресурсов , при которых прибыль предприятия обращается в нуль, находятся из уравнения
. (7)
Стохастический характер динамики развития производственного предприятия определяется объемами внутренних инвестиций и объемами амортизаций (износа) ресурсов, которые описываются случайными функциями времени. Поэтому приращение объемов ресурсов за некоторый малый промежуток времени можно выразить суммой трех компонентов:
. (8)
Здесь – частичные амортизационные утраты ресурсов за время , – частичные восстановления ресурсов за время за счет внутренних инвестиций предприятия, – случайные колебания объемов ресурсов предприятия.
Для пропорциональной амортизации величины записываются в виде
. (9)
Здесь – коэффициенты амортизации, выражающие доли утраченных объемов ресурсов за единицу времени.
Величины , выражающие частичные восстановления ресурсов за промежуток времени за счет внутренних инвестиций предприятия, имеют вид
. (10)
Здесь – внутренние инвестиции, восстанавливающие ресурсы; – нормы инвестиций, доли выручки, инвестируемые в ресурсы . Функция времени
описывает пропорциональный рост нормы инвестиций и амортизаций ресурсов на всем рассматриваемом интервале времени .
Имманентная волатильность объемов внутренних инвестиций является причиной случайных колебаний величины , которые могут быть представлены в виде стохастических стандартных винеровских процессов [36]:
. (11)
Здесь – стандартный винеровский процесс, , – показатель волатильности факторов производства , – случайная величина с нормальным законом распределения, нулевым средним значением и единичной дисперсией .
Их формулы (11) следует, что в окрестности начальной точки пространства и в окрестности предельной точки пространства поведение функций становится практически детерминированным.
Подставляя формулы (1), (9)–(11) в уравнения баланса (8), находим
. (12)
Предельный переход в соотношениях (12) при и приводит к системе стохастических дифференциальных уравнений диффузии Ито
. (13)
Здесь
(14)
– коэффициенты сноса системы (13),
(15)
– коэффициенты волатильности системы (13).
Начальные условия для системы уравнений (13) с коэффициентами (14) и (15) имеют вид
. (16)
Уравнения (13)–(15) показывают, что увеличения объемов производственных факторов и соответствующих им объемов выпуска продукции будут продолжаться до тех пор, пока объемы внутренних инвестиций будут превосходить объемы амортизационных отчислений. Рост величин прекратится, и они достигнут своих предельных значений , когда объемы внутренних инвестиций и амортизационных отчислений совпадут. Поскольку в окрестности предельной точки пространства случайный процесс становится практически детерминированным, то величины находятся из уравнения
. (17)
Целью любого производственного предприятия является организация такого режима его работы, при котором прибыль становится максимально возможной. Это достигается только в том случае, если все предельные величины ресурсов будут совпадать со значениями ресурсов , отвечающими максимальной прибыли .
В этом случае оптимальные коэффициенты нормы внутренних инвестиций находятся из соотношения (17)
. (18)
Численное решение системы стохастических дифференциальных уравнений (13) с коэффициентами (14), (15) и начальными условиями (16) строятся на временном отрезке разбитом системой точек методом последовательных приближений Эйлера – Маруямы в соответствии с алгоритмом [25]:
. (19)
Здесь .
Начиная с начальных значений ресурсов на каждом временном интервале с шагом генерируются случайные числа и вычисляются следующие значения . Полученные случайные последовательности и на координатной плоскости образуют системы точек и соответствующие им стохастические траектории. Повторная реализация алгоритма (19) всякий раз генерирует новые случайные значения величины ε и новые стохастические траектории.
Для вычисления математических ожиданий функций необходимо статистически усреднить систему уравнений (13) с коэффициентами (14) и (15):
. (20)
Рассмотрим важный частный случай , согласно которому число ресурсов сводится к одному, а выручка предприятия определяется одним фактором производства . В объем этого ресурса включены капитал, трудовые ресурсы, материалы, технологии, инновации и т. д.
Тогда формулы (1)–(3) принимают вид
, (21)
, (22)
(23)
Максимальное значение прибыли и соответствующее ему значение ресурса находятся из условия
. (24)
Решая уравнение (24), находим значение ресурса :
. (25)
Подставляя выражение (25) в формулу для прибыли (23), получаем максимальное значение прибыли:
(26)
Уравнение (7) для предельных значений ресурсов, при которых прибыль предприятия обращается в нуль, принимает вид
. (27)
Система стохастических уравнений динамики роста производственных факторов (13) с коэффициентами (14), (15) для однофакторного предприятия принимает вид
(28)
Начальное условие для системы (28) имеет вид
. (29)
Формула (18) для оптимального коэффициента нормы внутренних инвестиций записывается в виде
. (30)
Алгоритм (19) численного решения стохастического дифференциального уравнения (28) методом итераций Эйлера – Маруямы принимает вид
. (31)
При реализации алгоритма (31) получаются случайные системы точек и соответствующие им стохастические траектории.
Для вычисления математического ожидания функции необходимо статистически усреднить систему уравнений (28):
. (32)
Точное вычисление статистического момента не представляется возможным, поэтому приближенное значение этого момента вычисляется в предположении, что флуктуации величины относительно ее среднего значения пропорциональны случайной величине :
. (33)
Здесь
(34)
– коэффициент пропорциональности.
Подставляя соотношение (34) в формулу (33), находим
. (35)
В формуле (35) ограничимся тремя слагаемыми в разложении сходящегося биномиального ряда для малых флуктуаций :
. (36)
Усредняя соотношение (36), находим
. (37)
Подстановка выражения (37) в уравнение (32) приводит к дифференциальному уравнению относительно
. (38)
Начальное условие для уравнения (38) имеет вид
. (39)
Оптимальной организацией деятельности производственного предприятия является такой режим его работы, при котором нормы внутренних инвестиций рассчитывается по формуле (30), функция прибыли стремится к своему максимальному значению , а функция ресурса стремиться к значению , соответствующему этому максимальному значению прибыли. При выборе любой другой нормы внутренних инвестиций функция прибыли будет стремиться к другому меньшему предельному значению , соответствующему другому предельному значению ресурса .
На рисунке показано сравнение графиков стохастических траекторий и математических ожиданий функций объемов прибыли , построенных по формуле (23), с результатами численной реализации алгоритма (31) и результатами численного решения задачи Коши (38), (39) для коэффициентов нормы инвестиций и .
Расчетные значения: ; ; ; ; ; ; ; ; ; ; ; ; .
Графики функций объемов прибыли на рисунке показывают, что коэффициент нормы инвестиций выбран неудачно. После достижения максимального значения прибыль предприятия начинает снижаться.
Следует отметить, что кривые математических ожиданий, построенные по результатам решения задачи Коши (38), (39), и средние значения стохастических траекторий, построенных по результатам двухсот реализаций алгоритма (31), практически совпадают.
Рисунок – Сравнение графиков стохастических траекторий и математических ожиданий функций объемов прибыли построенных по формуле (23), с результатами численной реализации алгоритма (31) и результатами численного решения задачи Коши (38), (39) для коэффициентов нормы инвестиций и
Figure – Comparison of graphs of stochastic trajectories and mathematical expectations of functions of profit volumes constructed according to formula (23), the results of the numerical implementation of algorithm (31), and the results of the numerical solution of the Cauchy problem (38), (39) for the investment rate coefficients and
Заключение
- Разработана новая стохастическая модель динамики развития предприятий, восстановление производственных ресурсов которых обеспечивается за счет внутренних инвестиций.
- Для прогнозирования объемов производственных издержек, амортизационных отчислений, выручки и прибыли установлены системы стохастических дифференциальных уравнений.
- Показано, что эффективность динамики развития предприятия зависит от выбора значений коэффициентов норм инвестиций. При неудачном выборе этих коэффициентов производственные мощности предприятий не способны выйти на режим работы с максимальной прибылью.
- Получена система уравнений для вычисления эффективных коэффициентов норм инвестиций, при которых предприятия гарантированно выходят на режим работы с максимальной прибылью.
Об авторах
Елена Алексеевна Ильина
Самарский национальный исследовательский университет имени академика С.П. Королева
Автор, ответственный за переписку.
Email: elenaalex.ilyina@yandex.ru
ORCID iD: 0000-0002-2590-6138
кандидат физико-математических наук, доцент кафедры математики и бизнес-информатики
Россия, 443086, Российская Федерация, г. Самара, Московское шоссе, 34Леонид Александрович Сараев
Самарский национальный исследовательский университет имени академика С.П. Королева
Email: saraev_leo@mail.ru
ORCID iD: 0000-0003-3625-5921
доктор физико-математических наук, профессор, заведующий кафедрой математики и бизнес-информатики
Россия, 443086, Российская Федерация, г. Самара, Московское шоссе, 34Список литературы
- Harrod R.F. The Trade Cycle. Oxford: Clarendon Press, 1936. 234 p. URL: https://archive.org/details/tradecycle0000unse/mode/2up.
- Domar E.D. Capital Expansion, Rate of Growth, and Employment // Econometrica. April 1946, Vol. 14, no. 2. P. 137–147. DOI: https://doi.org/10.2307/1905364.
- Solow R.M. A Contribution to the Theory of Economic Growth // Quarterly Journal of Economics. February 1956. Vol. 70, issue 1. P. 65–94. DOI: https://doi.org/10.2307/1884513.
- Swan T.W. Economic Growth and Capital Accumulation // Economic Record. November 1956. Vol. 32. Issue 2. P. 334–361. URL: https://econpapers.repec.org/scripts/redir.pf?u=http %3A %2F %2Fhdl.handle.net%2F10.1111 %2Fj.1475-4932.1956.tb00434.x;h=repec:bla:ecorec:v:32:y:1956:i:2:p:334-361.
- Kuznets S. Long Swings in the Growth of Population and in Related Economic Variables // Proceedings of the American Philosophical Society. February 1958. Vol. 102, no. 1. P. 25–52. URL: https://www.jstor.org/stable/985303.
- Kuznets S. Quantitative Aspects of the Economic Growth of Nations. Paper VIII: Distribution of Income by Size // Economic Development and Cultural Change. 1963. Vol. 11. No 2. Part 2. P. 1–80. DOI: https://doi.org/10.1086/450006.
- Uzawa H. Optimum Technical Change in an Aggregative Model of Economic Growth // International Economic Review. 1965. Vol. 6, no. 1. P. 18–31. URL: http://links.jstor.org/sici?sici=0020-6598%28196501%296 %3A1 %3C18 %3AOTCIAA %3E2.0.CO %3B2-Y.
- Arrow K.J. The Economic Implications of Learning by Doing // The Review of Economic Studies. 1962. Vol. 29. Issue 3. P. 155–173. URL: https://econpapers.repec.org/scripts/redir.pf?u=http%3A%2F%2Fhdl.handle.net %2F10.2307 %2F2295952;h=repec:oup:restud:v:29:y:1962:i:3:p:155-173.
- Denison E.F. The Contribution of Capital to Economic Growth // Financing Industrial Investment. London: Palgrave Macmillan, pp. 39–87. DOI: https://doi.org/10.1007/978-1-349-04021-6_3.
- Romer P.M. Increasing Returns and Long-run Growth // Journal of Political Economy. October 1986, Vol. 94, number 5. P. 1002–1037. DOI: https://doi.org/10.1086/261420.
- Lucas R.E. On the Mechanics of Economic Development // Journal of Monetary Economics. July 1988, Vol. 22, no. 1. P. 3–42. DOI: https://doi.org/10.1016/0304-3932 %2888 %2990168-7.
- Romer P.M. Endogenous Technological Change // Journal of Political Economy. October 1990, Vol. 98. Number 5. Part 2. P. 71–102. DOI: https://doi.org/10.1086/261725.
- Grossman G.M., Helpman E. Innovation and Growth in the Global Economy. Cambridge, MA: MIT Press. 1991. 376 p. URL: https://books.google.ru/books?id=4ikgmM2vLJ0C&printsec=frontcover#v=onepage&q&f=false.
- Mankiw N., Romer D., Weil D. A Contribution to the Empirics of Economic Growth // Quarterly Journal of Economics. May 1992. Vol. 107, no. 2. P. 407–437. URL: http://links.jstor.org/sici?sici=0033-5533 %28199205 %29107 %3A2 %3C407 %3AACTTEO %3E2.0.CO %3B2-5.
- Grossman G.M., Helpman E. Endogenous Innovation in the Theory of Growth // Journal of Economic Perspectives. 1994. Vol. 8, no. 1. P. 23–44. DOI: http://dx.doi.org/10.1257/jep.8.1.23.
- Barro R.J., Sala-i-Martin X. Economic Growth. 2nd edition. Cambridge MA: MIT Press, 1995. 672 p. URL: http://piketty.pse.ens.fr/files/BarroSalaIMartin2004.pdf.
- Bruno M., Easterly W. Inflation Crises and Long-Run Growth: NBER Working Papers 5209. National Bureau of Economic Research, Inc, 1995. URL: https://www.nber.org/papers/w5209 (дата обращения: 06.03.2012).
- Gong G., Greiner A., Semmler W. The Uzawa – Lucas model without scale effects: theory and empirical evidence // Structural Change and Economic Dynamics. 2004. Vol. 15, issue 4. P. 401–420. DOI: https://doi.org/10.1016/j.strueco.2003.10.002.
- Нижегородцев Р.М. Модели логистической динамики как инструмент экономического анализа и прогнозирования //Моделирование экономической динамики: риск, оптимизация, прогнозирование. Москва, 1997. С. 34–51.
- Бадаш Х.З. Экономико-математическая модель экономического роста предприятия // Вестник Удмуртского университета. Серия Экономика и право. 2009. № 1. С. 5–9. URL: https://elibrary.ru/item.asp?id=11700881. EDN: https://elibrary.ru/jwbhyv.
- Королев А.В., Матвеенко В.Д. О структуре равновесных нестационарных траекторий в модели эндогенного роста Лукаса // Автоматика и телемеханика. 2006. № 4. С. 126–136. URL: https://www.mathnet.ru/rus/at1170.
- Кузнецов Ю.А., Мичасова О.В. Сравнительный анализ применения пакетов имитационного моделирования и систем компьютерной математики для анализа моделей теории экономического роста // Экономический анализ: теория и практика. 2007. № 5 (86). С. 23–30. URL: https://elibrary.ru/item.asp?id=9337066. EDN: https://elibrary.ru/hwikhf.
- Кузнецов Ю.А. Обобщенная модель экономического роста с учетом накопления человеческого капитала // Вестник Санкт-Петербургского университета. Прикладная математика. Информатика. Процессы управления. 2012. № 4. С. 46–57. URL: https://elibrary.ru/item.asp?id=18079557. EDN: https://elibrary.ru/pfqnbt.
- Прасолов А.В. Математические методы экономической динамики. Санкт-Петербург: Лань, 2015. 352 с. URL: https://klex.ru/uzv.
- Ильина Е.А., Сараев Л.А. К теории производственных функций, учитывающей изменение эластичностей выпуска по производственным ресурсам // Экономика и предпринимательство. 2018. № 10 (99). С. 145–150. URL: https://elibrary.ru/item.asp?id=35654399. EDN: https://elibrary.ru/sadnef.
- Сараев А.Л., Сараев Л.А. Показатели нелинейной динамики и предельное состояние производственного предприятия // Экономика и предпринимательство. 2018. № 11 (100). С. 1237–1241. URL: https://elibrary.ru/item.asp?id=36512728. EDN: https://elibrary.ru/ypfjhn.
- Сараев А.Л. Уравнения динамики нестабильных многофакторных экономических систем, учитывающих эффект запаздывания внутренних инвестиций // Казанский экономический вестник. 2015. № 3 (17). С. 68–73. URL: https://elibrary.ru/item.asp?id=24899060. EDN: https://elibrary.ru/uywnhn.
- Ильина Е.А., Сараев А.Л., Сараев Л.А. К теории модернизации производственных предприятий, учитывающей запаздывание внутренних инвестиций // Экономика и предпринимательство, 2017. № 9–4 (86). С. 1130–1134. URL: https://elibrary.ru/item.asp?id=30782945. EDN: https://elibrary.ru/zxqfaf.
- Сараев А.Л., Сараев Л.А. Экономико-математическая модель развития производственных предприятий, учитывающая эффект запаздывания внутренних инвестиций // Экономика и предпринимательство. 2019. № 5 (106). С. 1316–1320. URL: https://elibrary.ru/item.asp?id=39238012. EDN: https://elibrary.ru/aigtur.
- Сараев А.Л., Сараев Л.А. Многофакторная математическая модель развития производственного предприятия за счет внутренних и внешних инвестиций // Вестник Самарского университета. Экономика и управление. 2020. Т. 11, № 2. С. 157–165. DOI: https://doi.org/10.18287/2542-0461-2020-11-2-157-165. EDN: https://elibrary.ru/wdbmkv.
- Ilyina E.A., Saraev L.A. Predicting the dynamics of the maximum and optimal profits of innovative enterprises // Journal of Physics: Conference Series. 2021. Vol. 1784. P. 012002. DOI: https://doi.org/10.1088/1742-6596/1784/1/012002. EDN: https://www.elibrary.ru/xwxltx.
- Saraev A.L., Saraev L.A. Mathematical models of the development of industrial enterprises, with the effect of lagging internal and external investments // Journal of Physics: Conference Series. 2021. Vol. 1784, P. 012010. DOI: https://doi.org/10.1088/1742-6596/1784/1/012010. EDN: https://www.elibrary.ru/qvnrzq.
- Ито К., Маккин Г. Диффузионные процессы и их траектории. Москва: Мир, 1968. 329 с. URL: https://reallib.org/reader?file=507133&ysclid=lnlbqh7jsw514150694.
- Allen E. Modeling with Itô Stochastic Differential Equations. Springer: Mathematical Modelling: Theory and Applications, 2007. Vol. 22. 230 p. DOI: https://doi.org/10.1007/978-1-4020-5953-7.
- Кузнецов Д.С. Стохастические дифференциальные уравнения: теория и практика численного решения. Санкт-Петербург: Издательство Политехнический университет, 2009. 800 с. DOI: https://doi.org/10.18720/SPBPU/2/z17-4.
- Соловьев В.И. Экономико-математическое моделирование рынка программного обеспечения. Москва: ГУУ ВЕГА-ИНФО, 2009. 176 с. URL: https://mpra.ub.uni-muenchen.de/28974/1/MMPO.pdf; https://www.elibrary.ru/item.asp?id=19911710. EDN: https://www.elibrary.ru/qtsxxz.
- Кузнецова И.Ю. Численное решение стохастического дифференциального уравнения методом Эйлера – Маруямы // Международный исследовательский журнал. 2013. № 11–1 (18), С. 8–11. URL: https://www.elibrary.ru/item.asp?id=21041543. EDN: https://www.elibrary.ru/rsltld.
- Ильина Е.А., Парфенова А.Ю., Сараев Л.А. Влияние изменения общего объема рынка на кинетику процесса диффузии инноваций // Вестник Алтайской академии экономики и права. 2019. № 12–1. С. 61–67. DOI: https://doi.org/10.17513/vaael.848. EDN: https://www.elibrary.ru/kfhwfz.
- Сараев А.Л., Сараев Л.А. Математические модели стохастической динамики развития предприятий // Вестник Самарского государственного технического университета. Серия Физико-математические науки. 2020. Т. 24, № 2. С. 343–364. DOI: https://doi.org/10.14498/vsgtu1700. EDN: https://www.elibrary.ru/mltmba.