Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Основы экономической динамики (110

..pdf
Скачиваний:
1
Добавлен:
15.11.2022
Размер:
595.55 Кб
Скачать

где ρ — доля валовых инвестиций в ВВП;

1

 

1

 

 

 

 

 

> 1

— коэффициент усиления (мультипликатор), который показывает,

 

ρ

ρ

 

 

насколько должен быть увеличен ВВП для увеличения валовых инвестиций на единицу.

Таким образом, в широком смысле мультипликатор — усилительное линейное статическое звено, в узком смысле — сам коэффициент усиления.

Акселератор — дифференцирующее звено нулевого порядка, выход которого пропорционален скорости входа.

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

I = rdYdt ,

где r — коэффициент акселерации, т.е. прирост потребности в инвестициях при увеличении ВВП на единицу.

При дискретности времени ∆t или ∆t = 1 (один год) то же уравнение выглядит следующим образом:

It t = r t (Yt Ytt ), It = r (Yt Yt1 ) .

Инерционное звено задается дифференциальным уравнением первого порядка

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 + CeTt .

Константу С находим из начального условия

y (0) = x + C = y0 , C = y0 x ,

12

поэтому окончательно имеем y (t ) = x + (t0 x)eTt .

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

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 (1e1 ) , т.е. постоянную времени Т можно определить как длину промежутка времени, в течение которого переходный процесс проходит основную часть (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 (1c) y (t ) + I

t

,

 

 

 

где (1 - с) — склонность к накоплению.

14

 

При

 

t 0

приходим к уравнению инерционного звена (роль посто-

янной времени выполняет величина

1

, обратная склонности к накопле-

1- c

 

 

 

 

 

 

 

 

1

 

dy

C + I

 

 

нию):

 

 

dt + y

= 1c .

 

 

 

1c

 

 

Последнее уравнение имеет равновесное (стационарное) решение

yE = C + I .

1c

Если в начальный момент спрос на инвестиционные товары изменился с величины I0 до I ( I > I0), то в экономике будет происходить переходный

процесс от значения ВВП

y0

=

C + I0

до значения yE (см. рис. 1.4). При этом

1c

y (t ) = yE + ( y0 yE )e

 

) .

 

 

 

(

 

 

 

 

 

1 1c

 

 

 

 

 

3. Передаточная функция

Понятие передаточной функции динамического элемента связано с операторным методом решения дифференциального уравнения. Суть метода состоит в сведении решения дифференциального уравнения к решению алгебраического уравнения. В основе метода - переход от первоначальных функций времени x(t), y(t) к их образам X(s), Y(s) — преобразованиям Лапласа этих функций. Напомним определение преобразования Лапласа для некоторой функции :

 

 

 

 

F (s) = Lf

(s) = est f (t )dt ,

(3.1)

 

 

 

0

 

а также формулу обратного

перехода от образа к

прообразу:

f (t ) =

1

δ +iest F (s)ds .

(3.2)

 

 

2π i δ i

 

15

Образ производной можно найти по образу функции: поэтому

 

f (t )est dt = f (t )est

 

0+ s f

(t )est dt

 

0

 

 

0

 

Lf (s) = − f (0) + sF (s).

 

 

Lf (n) (s) = sn F (s) sn1 f (0) sn2 f (0) ... s(n1) f (0).

В частности, при f (0) = 0 ,

Lf = sF (s) .

При f (i) (0) = 0,

i = 0,..., n1, 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

 

 

 

 

 

 

 

 

 

 

 

1eα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

n1

где

N (s) = aj y(l) (0) sn1l ,

 

j=0

l =0

m

m1

 

M (s) = bi x(l ) (0) sm1l ,

откуда

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) = est y (t )dt — преобразование Лапласа выхода y(t);

0

x = est (xδ (t ))dt — преобразование Лапласа от правой части (4.5), по-

0

скольку estδ (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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]