CC BY
314
Математическое моделирование
УДК 517.958:330.101.541
МОДЕЛИРОВАНИЕ И АНАЛИЗ ДИНАМИКИ МАКРОЭКОНОМИЧЕСКОЙ СИСТЕМЫ В СРЕДЕ М ATL АВ / SIMULINK
В. С. Лубенцова, Д. А. Андреянов
Самарский государственный технический университет,
443100, Россия, Самара, ул. Молодогвардейская, 244.
E-mails: [email protected], [email protected]
Рассматривается общая теория циклических процессов в экономике. В среде MATLAB/Simulink построена имитационная модель решений дифференциальных уравнений, описывающих циклические колебания деловой активности и трендовую кривую роста ВВП. Исследована устойчивость системы. Приведено графическое изображение решений дифференциальных уравнений, проанализировано влияние мощности акселератора на динамику экономической системы.
Ключевые слова: деловые циклы, метод усреднения, мощность акселератора, уравнение макроэкономической динамики, устойчивость системы, точка бифуркации, долгосрочное прогнозирование.
В настоящее время в макроэкономике утвердилась теория, согласно которой деловые циклы играют значительную роль в формировании тренда экономического роста: эффективность долгосрочного роста является следствием краткосрочных циклических колебаний. Однако характер влияния деловых циклов на формирование траектории долгосрочного экономического роста недостаточно изучен. Последним научным достижением в этой области является создание непрерывных моделей. В работе А. А. Акаева [1] дан вывод общего уравнения макроэкономической динамики, описывающего совместное взаимодействие долгосрочного экономического роста и краткосрочных деловых циклов:
1 d2Y Г к ЛЖ + \1+А-Ь'
dYi
4/ dY^ п AdY' д¥
dY
dt
, ґ-, л / а\ с) У Ь <9У dA
‘-'А-8)жР ("“кУКт -4(1 “Ч1 Ж = Ті +кА' (1)
где Уіїй) —текущий объём (уровень) выпуска валовой продукции (ВВП); У(і) —уровень выпуска, соответствующий траектории долговременного роста; V — мощность акселератора; Л — скорость реакции запаздывания предложения от спроса; к — скорость реакции запаздывания фактических капиталовложений от момента принятия решения об индуцированных инвестициях; в — коэффициент сбережений; 7—параметр Оукена.
Применяя метод усреднения Крылова—Боголюбова—Митропольского [2], автор [1] из уравнения (1) выделил уравнение, представляющее долговременную траекторию экономического роста (2), и уравнение, описывающее циклические колебания
Вера Степановна Лубенцова (к.т.н., доц.), доцент, каф. прикладной математики и информатики. Денис Анатольевич Андреянов, студент, каф. прикладной математики и информатики.
деловой активности вокруг трендовой кривой роста (3):
d2Y _ dY (dA -\
~W + ao^+UJ°Y -X(^+kA)' (2)
d2y
dt2
3 \dtJ
dy_
dt
- Шл
У = v(t), (3)
где <то = —[A + к — kXv — A(1 — s)(3/7]; со2 = Ak; /3 = h/(b — a); у = Y — Y; <p — квази-периодическая функция, выражающая действие внешней вынуждающей силы; /3 — эластичность выпуска по труду; i — норма процента;_ао = А + к — kXv; А — трендовая составляющая независимых инвестиций (A{t) = A{t) + tp(t)). Для дальнейшего анализа были приняты следующие численные значения параметров: А = 4, к = 1, s = 0,25, i = 0,1, /3 = 0,67, 7 = 2,5.
Основным управляющим параметром, оказывающим существенное влияние на динамику исследуемой системы, является мощность v акселератора, которая меняется в определённых пределах. Здесь акселератор выступает как зависимость между объёмом капиталовложений и текущей скоростью изменения выпуска продукции dY /dt.
В уравнении (3) учитываются не только колебания выпуска продукции, но и колебания занятости. Поскольку циклическая безработица имеет негативные экономические и социальные последствия, очевидно, что экономическая политика должна быть направлена на ограничение уровня циклической безработицы и, следовательно, амплитуды колебания выпуска. Таким образом, при сокращении амплитуды колебаний совокупного выпуска продукции одновременно сокращается уровень циклической безработицы.
Количественный анализ динамики предложенной экономической системы проводился путём имитационного моделирования в среде MATLAB/Simulink. Для этого была построена модель заданной экономической системы (см. схему). Данная модель создаётся на основе стандартных блоков среды MATLAB/Simulink [3,4] и позволяет получить искомое значение путём двойного интегрирования блоками Integrator.
Для построения фазовых портретов и проведения качественного анализа в уравнениях (2) и (3) правые части полагаем равными нулю. В процессе моделирования варьировалась мощность акселератора и исследовалось её влияние на состояние равновесия системы. Результаты исследования динамики экономического роста (2) приведены в табл. 1.
Как видно из табл. 1, устойчивый рост возможен в диапазоне мощностей акселератора -у от 0 до 1,25, причём при 0 < v ^ 0,75 имеет место неколебательное (апериодическое) затухающее влияние циклических колебаний, при 0,75 < v < 1,25 — колебательное затухающее влияние, а при v > 1,25 — неустойчивый взрывной рост. В этом случае траектория выпуска стремительно отклоняется от траектории равновесия, но с течением времени взрывной рост заканчивается, начинается спад и устанавливается новое равновесие на новом уровне, которое будет удерживаться до очередного шока. При этом характер влияния циклических колебаний на долговременный рост определяется решением однородной части уравнения (2).
Результаты исследования влияния мощности акселератора на циклические колебания вокруг трендовой кривой приведены в табл. 2.
Как видно из табл. 2, при значениях мощности v > 1,05 любые малые колебания, возникающие в системе, перестают затухать и становятся нарастающими, пока не стабилизируются на уровне ~ у/&о, являющемся следствием нелинейных эффектов. Выходу на предельный цикл предшествует некоторый переходный процесс. При «то ~ 1 и v ~ 1,3 колебания будут существенно нелинейными. Таким образом, экономическая система с нелинейным акселератором является классической автоколебательной системой, в которой роль механизма положительной обратной связи играет нелинейный акселератор, а в качестве коэффициента усиления выступает мощность акселератора v. Если v > 1,05, то в системе будет протекать самоподдерживающийся
263
Схема модели дифференциальных уравнений (2) и (3)
Моделирование и анализ динамики макроэкономической системы .
Интервал изменения мощности акселератора
Характер решения однородного уравнения (2)
Таблица 1 Характерный график решения однородного уравнения (2)
0,75 < V < 1,25
1,25
1,25 < V < 1,75
V > 1,75
неколебательное,
затухающее
колебательное,
затухающее
колебательное,
незатухающее
колебательное,
взрывное
неколебательное,
взрывное
У
0,10
о
-0,10
-0,20
10 15 і
у ■ Ю15
0 5 10 15 і
Таблица 2
Продолжение таблицы 2
колебательный процесс, основные характеристики которого определяются параметрами системы и не зависят от начальных условий.
Результаты качественного анализа решений уравнений (2) и (3), описывающих траекторию долгосрочного экономического роста и циклические колебания деловой активности, позволяют проследить влияние деловых циклов на экономический рост.
Результаты моделирования показывают, что устойчивый долгосрочный экономический рост в динамической модели возникает в диапазоне мощностей акселератора О < V < 1,25. Причём при значениях V < 1,05 любые внешние шоковые воздействия сглаживаются системой и она, будучи выведена из равновесного состояния, асимптотически стремится к устойчивому первоначальному тренду экономического роста.
В интервале 1,05 < V < 1,25 состояние равновесия неустойчиво; в системе происходит бифуркация и рождается устойчивый незатухающий колебательный режим. Для этого режима характерны благоприятствующие структурные преобразования, необходимые для переработки и освоения инноваций. Экономическая система переходит к новому состоянию равновесия, которое фиксируется возникшим в системе автоколебанием. Следовательно, режим автоколебаний в экономической системе чрезвычайно полезен, поскольку способствует структурным преобразованиям, требуемым для восприятия инноваций.
При высоких значениях мощностей акселератора V > 1,25 экономический рост может стать неустойчивым вследствие взрывного характера влияния циклических колебаний системы. Однако благодаря цикличности V и естественному ограничению выпуска, вызванному падением предпринимательской прибыли, последующему спаду производства до наступления нового уровня равновесия, закрепляемого при снижении мощности акселератора, система будет продолжать относительно устойчивый рост, сменяя уровни равновесия и пребывая между ними в неравновесном состоянии, т. е. в состоянии динамического равновесия.
Компьютерное моделирование в среде МАТЬАВ/8л.пш1л.пк эффективно для понимания причинно-следственных связей, анализа и построений процессов экономики. Такое моделирование позволяет сделать исследование системы более наглядным и гибким к изменяющимся параметрам, что при аналитическом решении затруднено. Кроме того, имитационное моделирование позволяет решать задачу прогнозирования развития экономического процесса.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Акаев А. А. Анализ решений общего уравнения макроэкономической динамики // Эко-
номика и математические методы, 2008. Т. 44, №3. С. 62-78. [Akayev A. A. Analysis of solutions of general equations of microeconomic dynamics // Ekonomika i matematicheskie metody, 2008. Vol. 44, no. 3. Pp. 62-78].
2. Акаев А. А. Качественный анализ влияния деловых циклов на экономический рост // Экономика и математические методы, 2009. Т. 45, №3. С. 64-77. [Akayev А. А. Qualitative analysis of the impact of economic cycles on economic growth // Ekonomika i matematicheskie metody, 2009. Vol. 45, no. 3. Pp. 64-77].
3. Цисарь И. Ф. MATLAB Simulink. Компьютерное моделирование экономики. М.: Солон-Пресс, 2008. 256 с. [Cisar I. F. MATLAB Simulink. Computer modeling of economy. Moscow: Solon-Press, 2008. 256 pp.]
4. Hunt B. R., Lipsman R.L., Rosenberg J. M. A Guide to MATLAB: For Beginners And Experienced Users. Cambridge: Cambridge University Press, 2001. 348 pp.
Поступила в редакцию 23/VI/2011; в окончательном варианте — 20/1/2012.
MSC: 35Q91
MODELING AND ANALYSIS OF MACROECONOMIC SYSTEM DYNAMICS IN MATLAB/SIMULINK ENVIRONMENT
V. S. Lubencova, D. A. Andreyanov
Samara State Technical University,
244, Molodogvardeyskaya St., Samara, 443100, Russia.
E-mails: [email protected], [email protected]
The general theory of cyclic processes is examined in an economy. The simulation model of solutions of differential equations describing the cyclic vibrations of business activity and the trend curve of GDP growth is built in MATLAB/Simulink environment. Stability of the system is investigated. The graphic image of solutions of differential equations is brought, the influence of accelerator power on the dynamics of the economic system is analyzed.
Key words: business cycles, method of averaging, power of accelerator, the equation of macroeconomic dynamics, stability of system, bifurcation point, long-term prognostication.
Original article submitted 23/VI/2011; revision submitted 20/1/2012.
Vera S. Lubencova (Ph.D. (Techn.)), Associate Professor, Dept, of Applied Mathematics & Computer Science. Denis A. Andreyanov, Student, Dept, of Applied Mathematics & Computer Science.
