Многие приложения в экономике требуют моделирования систем во времени.
Состояние системы в момент времени t описывается мерным вектором X(t).
X1(t)
X(t) = ….. , X (t) Rn ( R – множество всех вещественных чисел)
Xn(t)
t
Эволюция системы со временем описывается функцией
G (X0, t, ) , где
X0 – начальное состояние системы;
t – время;
- вектор параметров.
Функция g(*) называют также переходной функцией
Функция g(*) – это правило, описывающее текущее состояние как функцию от времени, начальных условий и параметров.
Например: Xt = X0 (1+)t = g (X0, t, )
Функция g(*) как правило не известна. Обычно она задана неявно как решение системы разностных уравнений.
Разностное уравнение или система уравнений – это уравнения в следующей форме : F (t, Xt, Xt+1, …, Xt+m, ) = 0 (1),где
Xt – состояние системы в момент времени t.
Решение уравнения (1) – это последовательность векторов
Xt = X0, X1,…,
Обычно предполагается, что уравнение (1) можно решить аналитически относительно Xt+m и переписать в форме так называемых уравнений – состояний :
Xt+m = f (t, Xt, Xt+1, …,Xt+m-1, ) (2)
Например:
Xt+2 = Xt + Xt+1/2 + t
Любую систему представляют в форме (2) всегда можно ?
Разностное уравнение (2)называется линейным, если F(*) является линейной фуекцией переменных состояний (не обязательно линейно относительно )
В уравнениях (1) и (2) величина m называется порядком системыне является серьезным ограничением, так как системы более высокого порядка путем введения дополнительных переменных и уравнений.
Пример:Xt = f (Xt-1, Yt-1) – система 2-го порядка
Введем Yt = Xt-1
Xt = f(Xt-1, Yt-1)
Yt = Xt-1
Таким образом, мы будем рассматривать только системы 1-го порядка следующего вида:
Xt-1 = f(t, Xt, ) (3)
Уравнение (3) называется автономным, если t не входит в него отдельным аргументом.
Пример:
Рассмотрим динамику основных фондов на предприятии
Kt – стоимость основных фондов предприятия в период t.
- норма амортизации, то есть % основных фондов, которые изъяли на предприятии за год.
It = инвестиции в основные фонды.
Kt+1 = ( 1 - ) Kt + It – уравнение 1-го порядка, линейное, если It = I, тогда
Kt+1 = ( 1 - ) Kt + I – уравнение автономное
Если It = I(t) – неавтономное (зависит от t)
Решение уравнения (3) – это последовательность векторов состояния {Xt}, удовлетворяющих уравнению (3) для всех возможных состояний. Эта последовательность называется траекторией системы. Уравнение (3) показывает, как состояние системы изменяется от периода к периоду, а траектория системы дает ее эволюцию как функцию начальных условий и состояния внешней среды .
Если известно начальное состояние X0, легко получить последовательность решений путем итеративного применения отношения (3), получим переходную функцию следующим образом:
Xt+1 = f (t, Xt, )
{Xt}
X0 = X0
X1 = f ( 0, X0, ) = g (0, X0, )
X2 = f ( 1, X, ) = f (1; f (0, X0, );) = g (1, X0, )
Xt+1 = f (t, Xt, ) = f ( t, g, (t – 1, X0, ),) = g (t, X0, )
Если f (*) однозначная, всюду определенна функция, то существует уникальное решение уравнения (3) для любого X0.
Если функция имеет вид f (t, Xt, ) = / Xt – не всюду опрделенная.
Если f (*) непрерывная дифференциальная функция, то решение также будет гладким относительно и X0
Полученное решение зависит от начального состояния X0.
Задача с граничным условием состоит из уравнения (3) и граничного условия, задаваемого в формуле:
Xs = Xs (4)
Если в уравнении (4) – S = 0 , то оно называется начальным состоянием.
Уравнение (3) имеет много решений, а уравнение (3) + (4) – система – единственное решение, поэтому различают общее и частное решение разностного уравнению (3):
Xtg = X(t, c, ) = {Xt(Xt+1 = f (t, Xt, ))} , где параметр е индексирует частное решение.
Пример:
Xt – размер вклада в момент t
Z - % я ставка
Xt+1 = Xt ( 1+ z ) ; X0 = …
X1 = X0 ( 1 + z)
X2 = X1 ( 1 + z) = X0 ( 1 + z)2 = g (X0, t, z) , где t = 2
Если можно найти общее решение системы (3) . у нас будет полная информация о поведении системы со временем, будет легко определить, как система реагирует на изменение параметров.
К сожалению, общее решение существует только для определенных классов l – го порядка ( в частности для линейных систем )