Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

DYNAMICS OF CHANGE OF NATIONAL INCOME IN THE MODEL OF HARMONIC OSCILLATOR

Gevorkyan E.A. 1 Sinchukov A.V. 1 Tatarnikov O.V. 1
1 Plekhanov Russian University of economics (REU)
In this article the dependence of national income on time Y(t) in the harmonic oscillator model is investigated. It is assumed that national income, as a function of time, satisfies a second-order homogeneous ordinary differential equation with constant real coefficients. It is also shown that such ordinary differential equation describing change of national income with time can be received on the basis of dynamic model of Keynes taking into account that the velocity of change of spendings and final consumption is proportional to national income. The economic meaning of each term in the differential equation of the harmonic oscillator is determined. The solution of this equation by the Euler method leads to an analytical expression of the national income Y(t). It allows to carry out the qualitative and graphic analysis of dynamics of change of national income at various values of the parameters characterizing wave or non-wave process of change of national income. It is shown, that in a certain relation of the parameters, characterizing the process, the change of national income has oscillatory character. At the same time the amplitude of the oscillations decreases with increasing time and the envelope of the maxima changes exponentially. In cases of non-oscillatory nature of changes in national income when time tends to infinity national income tends to zero.
national income
harmonic oscillator
differential equation
oscillatory process

В макроэкономической теории рассматриваются волновые явления в экономике, когда под действием каких-то факторов происходит отклонение макроэкономических показателей от их устойчивого состояния или от их траектории движения [1–2]. Среди таких показателей особое место занимает национальный доход, являющийся одним из обобщающих показателей экономического развития страны. Волновой характер изменения экономики в зависимости от времени является формой прогрессивного развития рыночной экономики. При этом каждая волна характеризует целый цикл экономического развития [3–5]. Отметим, что в основе теории экономического цикла лежит выяснение причин колебаний экономической активности общества во времени. Отметим также, что так как цикличность прямым или косвенным образом действует на рыночную экономику, то она является важной макроэкономической проблемой.

При построении динамической модели волнового изменения национального дохода можно пользоваться различными математическими моделями, описывающимися при помощи обыкновенных дифференциальных уравнений [6].

Цель исследования: найти зависимость национального дохода от времени в рамках модели гармонического осциллятора и провести многосторонний анализ полученного результата.

Постановка задачи, метод решения и результаты исследования

Пусть национальный доход Y(t) как функция от времени t удовлетворяет обыкновенному однородному дифференциальному уравнению второго порядка с постоянными действительными коэффициентами (модель гармонического осциллятора)

gev02.wmf (1)

где первый член gev03.wmf показывает темп изменения национального дохода (в физике это соответствует ускорению движения), второй член gev04.wmf связан с транзакционными издержками (издержки на переговоры, на юридическую защиту и т.д.) (в физике это соответствует силе трения), а третий член gev05.wmf есть рыночная сила (в физике это сила, которая возвращает систему к точке равновесия), γ – коэффициент затухания, ω0 – частота свободных колебаний (γ = 0).

Любопытно отметить, что аналогичное дифференциальное уравнение для функции Y(t) можно получить исходя из динамической модели Кейнса [7]. В основе этой модели лежит основной закон экономического баланса, который имеет вид

gev06.wmf (2)

где Y(t) – национальный доход, E(t) – расходы, S(t) – национальное потребление, I(t) – инвестиции. Представляя функцию потребления в виде суммы

gev07.wmf (3)

где B(t) – конечное потребление, 0 < c < 1 – постоянный коэффициент, и учитывая, что gev08.wmf (m = const), а функции E(t) и B(t) известны, из (2) относительно функции Y(t) получим следующее обыкновенное дифференциальное уравнение первого порядка

gev09.wmf (4)

Дифференцируя обе части уравнения (4) по t и предполагая, что

gev10.wmf

где b – постоянный коэффициент, получим

gev11.wmf (5)

Сравнивая (5) с (1), можно сказать, что национальный доход Y(t) описывается дифференциальным уравнением гармонического осциллятора с частотой свободных колебаний gev12.wmf и коэффициентом затухания gev13.wmf

Решение линейного однородного дифференциального уравнения (1) будем искать методом Эйлера, а именно в виде [8]

gev14.wmf (6)

где λ – пока неизвестные числа (действительные, комплексные или чисто мнимые). Подставляя (6) в (1) для определения характеристических чисел, получим следующее квадратное уравнение

gev15.wmf (7)

решения которого имеют вид

gev16.wmf (8)

Из (8) видно, что характер решений уравнения (7) зависит от знака выражения под квадратным корнем. Рассмотрим случаи:

1. gev17.wmf (коэффициент затухания меньше частоты свободных колебаний).

Очевидно, что в этом случае имеем комплексные корни характеристического уравнения (3) в виде

gev18.wmf (9)

Тогда линейно независимые решения дифференциального уравнения (1) имеют вид

gev19.wmf

gev20.wmf (10)

а общее решение, как суперпозиция независимых решений, запишется в виде

gev21.wmf (11)

где c1 и c2 – произвольные постоянные. Нетрудно заметить, что (11) можно представить в виде

gev22.wmf (12)

где gev23.wmf – амплитуда колебаний, gev24.wmf – начальная фаза. Выражение (8) показывает, что при γ < ω0 и γ > 0 национальный доход в зависимости от времени имеет колебательный характер с резко выраженными максимумами и минимумами. График зависимости национального дохода от времени приведен на рис. 1. Кривые построены согласно формуле 2,5 + Y(t) при изменении t от 0 до 24 (t измеряется в месяцах) и при следующих значениях параметров, характеризующих рассматриваемый процесс: c1 = c2 = 1, ω0 = 1, γ = 0,07; 0,08; 0,1.

gevor1.tif

Рис. 1. Зависимость функции 2,5 + Y(t) от времени t при ω0 = γ:1 – ω0 = 1, γ = 0,07; 2 – ω0 = 1, γ = 0,08; 3 – ω0 = 1, γ = 0,1

Полученные кривые показывают, что огибающие максимумов с возрастанием времени убывают по экспоненциальному закону и при t→∞ стремятся к нулю. Из рис. 1 также следует, что с увеличением коэффициента затухания γ при фиксированном значении частоты собственных колебаний ω0 значения максимумов уменьшаются.

2. gev26.wmf или γкр. > ω0 (затухание критическое, так как при γ > ω0 осциллятор будет совершать колебательное движение).

В этом случае характеристическим числам gev27.wmf соответствуют два независимых решения

gev28.wmf gev29.wmf. (13)

И общее решение дифференциального уравнения (1) представится в виде

gev30.wmf (14)

Из (14) следует, что процесс имеет неколебательный характер, и так как gev31.wmf то осциллятор экспоненциально будет стремиться к положению равновесия. На рис. 2 приведены графики зависимости национального дохода Y как функция от времени t, построенные согласно формуле (14) при следующих значениях параметров: c1 = c2 = 1, ω0 = γ = 0,07; 0,08; 0,1 и 0 ≤ t ≤ 24. Кривые показывают убывающий характер национального дохода как функция от времени, причем при увеличении значений ω0 = γ скорость стремления к нулю при возрастании t увеличивается.

gevor2.tif

Рис. 2. Зависимость национального дохода Y от времени t при ω0 = γ:1 – ω0 = γ = 0,07; 2 – ω0 = γ = 0,08; 3 – ω0 = γ = 1

gevor3.tif

Рис. 3. Зависимость национального дохода Y от времени t при γ > ω0:1 – ω0 = 0,6; γ = 1; 2 – ω0 = 0,5; γ = 1; 3 – ω0 = 0,2; γ = 1

3. γ > ω0 (соответствует большим транзакционным затратам).

В этом случае характеристическое уравнение (3) имеет два различных действительных решения и общее решение дифференциального уравнения (1) будет иметь вид

gev32.wmf (15)

Так как при

gev33.wmf

то имеем неколебательный процесс и осциллятор будет стремиться к положению равновесия. Зависимости национального дохода Y от времени t в этом случае приведены на рис. 3.

Как видно из рисунка, национальный доход как функция от времени при γ > ω0 является убывающей функцией, при этом чем больше ω0, тем быстрее функция стремится к нулю при t→∞.

Заключение

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