Основы экономической динамики (110
..pdfгде ρ — доля валовых инвестиций в ВВП;
1 |
|
1 |
|
|
|
|
|
> 1 |
— коэффициент усиления (мультипликатор), который показывает, |
|
ρ |
|||
ρ |
|
|
насколько должен быть увеличен ВВП для увеличения валовых инвестиций на единицу.
Таким образом, в широком смысле мультипликатор — усилительное линейное статическое звено, в узком смысле — сам коэффициент усиления.
Акселератор — дифференцирующее звено нулевого порядка, выход которого пропорционален скорости входа.
Например, инвестиции могут быть выражены через скорость изменения ВВП следующим образом:
I = rdYdt ,
где r — коэффициент акселерации, т.е. прирост потребности в инвестициях при увеличении ВВП на единицу.
При дискретности времени ∆t или ∆t = 1 (один год) то же уравнение выглядит следующим образом:
It t = r t (Yt − Yt− t ), It = r (Yt − Yt−1 ) .
Инерционное звено задается дифференциальным уравнением первого порядка
a1 |
dy |
+ a0 y = x (t ) , a0 ≠ 0. |
(2.2) |
|
dt |
|
|
Уравнение (2.2) можно привести к стандартному виду путем деления его на a0.
T dy |
+ y = x (t ) |
(2.3) |
dt |
, |
где T = |
a1 |
, |
x (t ) = |
x (t ) |
. |
a0 |
|
||||
a |
|||||
|
|
|
|
0 |
|
11
Инерционное звено описывает процесс «отработки» заданного входного воздействия x(t) (значок «~» опустим), таким образом, что скорость «отработки» пропорциональна разности между входом и выходом:
dy = |
1 |
x (t ) − y (t ) |
|
|
|||
dt |
T |
|
. |
|
|
Пример 2.1. Рассмотрим модель освоения введенных производственных мощностей. Обозначим через х (х = const) введенную производственную мощность, а через y(t) — фактическое производство на базе этой мощности в момент t (фактическое использование мощности, y(t) < x). Сделав предположение, что прирост производства пропорционален недоиспользованной мощности:
y = y(x - y) t ,
приходим к уравнению инерционного звена:
T dy + y = x |
, |
T = |
1 |
, |
y (0) = y |
, |
y < x |
. |
(2.4) |
|
|||||||||
at |
|
y |
0 |
0 |
|||||
|
|
|
|
|
|
|
|
Общее решение однородного уравнения имеет вид:
T dy |
+ y = 0 , y = ceλt . |
(2.5) |
dt |
|
|
Подставив его в (2.5), получим:
(Tλ + 1) y = 0 ,
но у ≠0, поэтому приходим к характеристическому уравнению (относительно X):
Tλ + 1 = 0 , или |
λ = − |
1 |
. |
|
|||
|
|
T |
Поскольку частным решением неоднородного уравнения (2.4) является y = x , то общее решение этого уравнения примет вид:
y = x + Ce− Tt .
Константу С находим из начального условия
y (0) = x + C = y0 , C = y0 − x ,
12
поэтому окончательно имеем y (t ) = x + (t0 − x)e− Tt .
Переходный процесс освоения производственных мощностей, описываемый этим решением, завершается выходом на заданный размер мощности
lim y (t ) = x .
t→∞
Общая картина переходного процесса показана на рис. 1. 4.
y
x
|
|
|
|
|
|
y (t ) = x + ( y0 − x )e− |
t |
|
|
|
|
|
|
|
T |
||
y0 |
|
|
|
|
|
|
||
|
1 |
− e |
− |
t |
|
|
||
|
y (t ) = x |
T |
|
|
||||
|
|
|
|
|
|
|
|
0 |
|
|
|
|
t |
Рис. 1.4. Переходный процесс освоения производственных мощностей |
|||||
При у0 = 0 решение примет вид: |
|
|
|
|
|
y (t ) |
|
− |
t |
|
, |
|
|||||
= x 1 |
− e T |
||||
|
|
|
|
|
|
поэтому y (t ) = x (1− e−1 ) , т.е. постоянную времени Т можно определить как длину промежутка времени, в течение которого переходный процесс проходит основную часть (2/3) своего пути от 0 до х.
Пример 2.2. Модель установления равновесной цены. В модели рассматривается рынок одного товара, время считается непрерывным, спрос d и предложение s линейно зависят от цены:
d = a − bp , s = α + β p , a > 0, b > p, α > 0, β > 0, a > α .
Основное предположение модели состоит в том, что изменение цены пропорционально превышению спроса над предложением:
13
p = γ (d − s) t , γ > 0 ,
т.е. в случае действительного превышения спроса над предложением цена возрастает, в противном случае — падает.
Из основного предположения модели вытекает следующее дифференциальное уравнение для цены:
1 |
|
dp |
+ (b + β ) p = a − α , p (0) = p0 , |
γ |
|
dt |
|
т.е. процесс описывается уравнением инерционного звена с
T = |
1 |
и |
pE = a − α |
, |
|
γ (b + β ) |
|||||
|
|
b + β |
|
где pE — равновесная цена (точка пересечения прямых спроса и предложения). Таким образом, цена как выход инерционного звена ведет себя так, как это показано на рис. 1.4.
В модели Кейнса предполагается, что ВВП y(t + 1) в следующем году равен совокупному спросу предыдущего (текущего) года, а совокупный спрос, состоящий из спроса на потребительские (С) и инвестиционные (I) товары, зависит только от ВВП текущего года:
y (t + 1) = C y (t ) + I (t ) .
При линейной зависимости спроса на потребительские товары от ВВП и примерном постоянстве спроса на инвестиционные товары приходим к соотношению
y (t + 1) = C + cy (t ) + I , |
(2.6) |
где C — минимальный объем фонда потребления; с(0< с < 1) — склонность к потреблению.
Соотношение, действующее при дискретности времени в один год, примет форму:
y (t + t ) − y (t ) = C − (1− c) y (t ) + I |
t |
, |
|
|
|
|
где (1 - с) — склонность к накоплению.
14
|
При |
|
t → 0 |
приходим к уравнению инерционного звена (роль посто- |
||||
янной времени выполняет величина |
1 |
, обратная склонности к накопле- |
||||||
1- c |
||||||||
|
|
|
|
|
|
|
||
|
1 |
|
dy |
C + I |
|
|
||
нию): |
|
|
dt + y |
= 1− c . |
|
|
||
|
1− c |
|
|
Последнее уравнение имеет равновесное (стационарное) решение
yE = C + I .
1− c
Если в начальный момент спрос на инвестиционные товары изменился с величины I0 до I ( I > I0), то в экономике будет происходить переходный
процесс от значения ВВП |
y0 |
= |
C + I0 |
до значения yE (см. рис. 1.4). При этом |
|||
1− c |
|||||||
y (t ) = yE + ( y0 − yE )e |
|
) . |
|
|
|
||
( |
|
|
|
|
|||
|
−1 1−c |
|
|
|
|
|
3. Передаточная функция
Понятие передаточной функции динамического элемента связано с операторным методом решения дифференциального уравнения. Суть метода состоит в сведении решения дифференциального уравнения к решению алгебраического уравнения. В основе метода - переход от первоначальных функций времени x(t), y(t) к их образам X(s), Y(s) — преобразованиям Лапласа этих функций. Напомним определение преобразования Лапласа для некоторой функции :
|
|
|
∞ |
|
F (s) = Lf |
(s) = e− st f (t )dt , |
(3.1) |
||
|
|
|
0 |
|
а также формулу обратного |
перехода от образа к |
прообразу: |
||
f (t ) = |
1 |
δ +i∞ est F (s)ds . |
(3.2) |
|
|
||||
|
2π i δ −i∞ |
|
15
Образ производной можно найти по образу функции: поэтому
∞ |
∞ |
|
||
f ′ (t )e− st dt = f (t )e− st |
|
0∞ + s f |
(t )e− st dt |
|
|
||||
0 |
|
|
0 |
|
Lf ′ (s) = − f (0) + sF (s). |
|
|
||
Lf (n) (s) = sn F (s) − sn−1 f ′ (0) − sn−2 f ′ (0) − ... − s(n−1) f ′ (0). |
||||
В частности, при f (0) = 0 , |
Lf ′ = sF (s) . |
|||
При f (i) (0) = 0, |
i = 0,..., n−1, Lf (n) (s) = sn F (s) . |
В табл. 1.1 приведены преобразования Лапласа некоторых функций. Таблица 1.1. Преобразования Лапласа типовых функций
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f (t ),t > 0 |
|
|
|
F(s) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
δ (t) |
1 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
χ (t) |
1 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
s |
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
1 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
s2 |
|
|
|
|
|
|
|
|
tn |
1 |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
sn+1 |
|
|
|
|
|
|
||
|
n! |
|
|||||||||||||
e−αt |
1 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
s + a |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|||||
te−αt |
1 |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
(s + a)2 |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|||||||
1− e−αt |
|
|
|
|
|
a |
|||||||||
|
|
|
|
|
s(s + a) |
|
|
||||||||
|
|
|
|
|
|
|
|
||||||||
sinωt |
|
|
|
|
|
ω |
|||||||||
|
|
|
|
|
s2 + ω 2 |
|
|
||||||||
|
|
|
|
|
|
|
|
||||||||
cosωt |
|
|
|
|
|
s |
|||||||||
|
|
|
|
|
s2 + ω 2 |
|
|
||||||||
|
|
|
|
|
|
|
|||||||||
e−α t sinωt |
|
|
|
|
|
ω |
|||||||||
|
|
|
|
(s + a)2 + ω 2 |
|
||||||||||
|
|
|
|
|
|||||||||||
e−α t cosωt |
|
|
|
s + a |
|||||||||||
|
|
|
|
(s + a)2 + ω 2 |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16
Применим преобразование Лапласа к обеим частям уравнения динамического элемента (2.1) (пользуясь формулой для образа производных):
|
n |
|
|
m |
|
X (s) − M (s), |
|
aj s j Y (s) − N (s) = |
ai si |
||||
|
j=0 |
|
i=0 |
|
|
|
n |
n−1 |
где |
N (s) = aj y(l) (0) sn−1−l , |
|
|
j=0 |
l =0 |
m |
m−1 |
|
M (s) = bi x(l ) (0) sm−1−l , |
откуда |
|
i=0 |
l =0 |
|
|
|
m |
|
|
|
|
|
|
Y (s) = G (s) X (s) + R (s), где |
G (s) = |
bi si |
|
|
R (s) = |
N (s) − M (s) |
|
R(s) = 0 при |
i=0 |
|
, |
n |
, |
||||
n |
j |
|
a j s j |
|||||
|
|
aj s |
|
|
|
j=0 |
|
|
|
|
j=0 |
|
|
|
|
|
|
y( j) (0) = x(i) (0) = 0, j = 0,1,..., n −1, i = 0,1,..., m −1.
Передаточной функцией G(s) динамической системы (подсистемы, элемента) называется отношение образа выхода к образу входа при нулевых условиях.
В передаточной функции динамической системы (подсистемы, звена) содержатся все сведения о ее поведении при нулевых начальных условиях. В самом деле, по входу x(t) находим его образ X(s), затем умножаем этот образ на передаточную функцию, тем самым получаем образ выхода Y (s) = G(s)X (s) и, наконец, пользуясь либо табл. 1.1, либо непосредственно по определению находим выход y(t). Если начальные условия ненулевые, то к этому решению еще добавится «шлейф», образ которого — R(s).
4. Колебательное звено
Колебательное звено задается дифференциальным уравнением второго порядка
|
d |
2 |
2y + a1 |
dy |
m |
|
a2 |
|
+ a0 y = bi x(i) (t ) |
(4.1) |
|||
|
dt |
|
dt |
i=0 |
|
17
с отрицательным дискриминантом, составленным из коэффициентов в левой части уравнения (4.1).
Колебательное звено описывает циклические процессы в экономике. Пример 4.1. Однономенклатурная система управления запасами может быть рассмотрена как колебательное звено. Пусть x(t) и x(t) — фактические интенсивности расхода и поступления товара в систему управления запасами в момент t. Поскольку интенсивность расхода заранее неизвестна, то всегда будет образовываться запас y(t) (если y(t) > 0, то это действительно запас, если y(t)<0, то это дефицит). Изменение запаса следующим обра-
зом связано с интенсивностями расхода и поставок:
t = (x − x) t или |
dy |
= x − x . |
(4.2) |
|
dt |
|
|
Управлять интенсивностью поставок можно только по известному значению запаса y(t) (ведь интенсивность расхода неизвестна).
Имеется два варианта управления:
1) изменение поставок пропорционально (с обратным знаком) величине запаса (при положительном запасе интенсивность поставок уменьшается, при отрицательном — увеличивается):
x = −a0 y t , a0 > 0 .
2) изменение интенсивности поставок пропорционально (с обратным знаком) как запасу, так и скорости его изменения:
|
dy |
t , |
a0 > 0, |
a1 > 0, |
a0 |
> |
a12 |
x = − a0 y + a1 |
|
4 |
|||||
|
dt |
|
|
|
|
|
(при положительном запасе интенсивность поставок уменьшается, при отрицательном — увеличивается, при положительной скорости роста запаса интенсивность поставок уменьшается, при отрицательной — увеличивается).
18
Первый случай. Взяв производную от обеих частей (4.2) и подста-
вив в это выражение |
dy |
= −a0 y , получаем дифференциальное уравнение вто- |
||
|
dt |
|
|
|
рого порядка для запаса: |
|
|
||
|
|
d 2 2y |
+ a0 y = − dx . |
(4.3) |
|
|
dt |
dt |
|
Это уравнение колебательного звена с a2 = 1, |
a1 = 0 и дискриминан- |
том d = −4a0 < 0 . Характеристическое уравнение имеет вид (подставляем в однородное уравнение y = Ceλt ): λ 2 + a0 = 0 . Его корни взаимно сопряженные мнимые: λ1 = i a0 , λ2 = −i a0 .
Пусть на вход системы, находившейся в начальный момент в состоянии равновесия х = 0, у = 0, у'(0) = 0, начали поступать заявки на товар с интенсивностью x(t) =x = const.
Или алгебраически:
x (t ) = x χ (t ) , |
χ (t ) = 1, t > 0 |
, |
|
0, t < 0 |
|
где χ (t ) — функция Хэвисайда. |
|
|
Производная от функции Хэвисайда равна обобщенной функции Дирака, которая принимает бесконечно большое значение в точке t = 0 , равна
∞ |
|
|
|
|
|
нулю при t ≠ 0 и δ (t )dt = 1. |
|
|
|
|
|
−∞ |
|
|
|
|
|
Поскольку χ ′ (t ) = δ (t ) , то |
в этой ситуации |
dy |
= xδ (t ) , |
и уравнение (4.3) |
|
|
|
dt |
|
|
|
принимает вид: |
|
|
|
|
|
d 2 2y |
+ a0 y = − xδ (t ), y (0) = 0, |
y′ (0) = 0 |
. |
(4.4) |
|
dt |
|
|
|
|
|
Решим это уравнение операторным методом, применив преобразова- |
|||||
ние Лапласа к обеим частям уравнения: |
|
|
|
|
|
|
(s2 + a0 )Y (s) = − x , |
|
(4.5) |
||
|
19 |
|
|
|
|
∞
где Y (s) = e− st y (t )dt — преобразование Лапласа выхода y(t);
0
∞
− x = e− st (− xδ (t ))dt — преобразование Лапласа от правой части (4.5), по-
0
∞
скольку e− stδ (t )dt = 1.
0
Из (4.5) находим преобразование Лапласа выхода y(t):
Y (s) = − |
|
|
x |
|
|
, ω = a0 . |
s |
2 |
+ |
ω |
2 |
||
|
|
|
|
И, наконец, по табл. 1.1 восстанавливаем выход: v (t ) = ωx sin ωt .
Таким образом, в первом случае при постоянной интенсивности расхода х запас y(t) будет испытывать незатухающие гармонические колебания
с амплитудой ωx (рис. 1.5).
При таких незатухающих колебаниях промежутки, когда имеется действительный запас y(t) > 0, будут чередоваться с промежутками дефицита y(t) < 0, что крайне отрицательно скажется на финансовом положении организации, отвечающей за систему управления запасами. Для того чтобы система управления запасами снова вошла в состояние равновесия, необходимо учитывать не только величину запаса y(t), но и скорость его изменения, как это и предусмотрено во втором случае.
20