В дисциплине "Основы эконометрики" тест 6 дается по теме 7. |
Разностные уравнения и их решение
9.1. Уравнения первого и второго порядков
Пусть дана функция . Тогда первая разность
y определяется по формуле
В отличие от дифференциального исчисления мы не предполагаем, что стремится к нулю. Так как большинство экономических данных собираются через определенные равные промежутки времени, полезно считать, что
. Более того, обычно полагают период между наблюдениями
нормализованным, т.е. выбирают
(один месяц, один квартал, один год). Поэтому первые разности можно записать в виде
![\Delta y_{t} = f(t) - f(t - 1) = y_{t} - y_{t - 1},\\
Delta y_{t + 1} = f(t + 1) - f(t) = y_{t + 1} - y_{t},\\
\Delta y_{t + 2} = f(t + 2) - f(t + 1) = y_{t + 2} - y_{t + 1}](/sites/default/files/tex_cache/4576007d164d2a23ff4d55adc6d3c20c.png)
и т.д.
Таким же образом сформируем вторую разность как изменение первой разности:
![\Delta ^{2}y_{t} = \Delta (\Delta y_{t}) = \Delta (y_{t} - y_{t - 1}) = (y_{t} - y_{t - 1}) - (y_{t} - 1 - y_{t - 2}) = y_{t} - 2y_{t - 1} + y_{t - 2},\\
Delta ^{2}y_{t + 1} = y_{t + 1} - 2y_{t} + y_{t - 1}](/sites/default/files/tex_cache/210f51e8c83668168e92741b61278d0c.png)
и т.д.
Аналогично определяется разность n-го порядка:
![\Delta ^{n}y_{t} = \Delta (\Delta ^{n-1}y_{t}).](/sites/default/files/tex_cache/940fe50cff30f64b19dfe264be30d849.png)
Так как в настоящем учебном пособии в основном рассматриваются линейные модели временных рядов, изучим только один специальный случай линейного разностного уравнения -го порядка с постоянными коэффициентами, а именно уравнение вида
Порядок разностного уравнения задается показателем - максимальной величиной шага запаздывания, или, другими словами, максимальным лагом. Уравнение считается линейным потому, что все значения зависимой переменной y входят в него в первой степени.
Коэффициенты называются параметрами уравнения и не зависят от значений
и
. Переменная
представляет возмущающий процесс и может зависеть от времени, текущего и прошлого значений других переменных или иметь стохастический характер изменения. Соответственно, выбирая возмущающий процесс x, можно получить широкий спектр важных макроэкономических моделей.
Рассмотрим, например, стохастическую версию классической кейнсианской модели, предложенную Самуэльсоном:
где
Слагаемые и
имеют нулевые средние и объясняют случайные возмущения в потреблении и инвестировании.
Преобразуем кейнсианскую модель производства и потребления (9.3) к виду
Уравнение (9.4) выражает как функцию от собственных запаздываний (лагов) и возмущающих членов. Выберем возмущающий процесс в виде
![x_{t} = (1 + \beta ) \varepsilon _{ct} + \varepsilon _{it} - \beta \varepsilon _{it}_{-1}.](/sites/default/files/tex_cache/26477997fd7cc2e05c621b96d8e05a25.png)
и получим разностное линейное уравнение второго порядка вида (9.2). Отметим, что в уравнении (9.4) отсутствует свободный член, т.е. .
Важный частный случай для последовательности {x_{t}} получаем при где
- константы,
не зависят от
. При этом можно полагать, что
является последовательностью неопределенных внешних (экзогенных) переменных. Например, если
- последовательность случайных ошибок наблюдений и
, то уравнение (9.2) становится уравнением авторегрессии
![y_{t} = a_{0} + a_{1}y_{t - 1} +\dots + a_{n}y_{t} - n + \varepsilon _{t}.](/sites/default/files/tex_cache/77e3f00e62c5458935ba2a8ac98ddd35.png)
Пусть . Тогда получаем модель случайного блуждания.
Другой полезный пример дает уравнение
Можно проверить, что решением этого разностного уравнения первого порядка является функция
Метод итераций (последовательных приближений). Обозначим известное значение функции в момент времени 0 как
. Обратимся к разностному уравнению первого порядка
Подставляя в (9.7), получаем:
![y_{1} = a_{0} + a_{1}y_{0} + \varepsilon _{1}.](/sites/default/files/tex_cache/af5103cbd395521aa053272e6a00f563.png)
Тем же путем находим :
![y_{2} = a_{0} + a_{1}y_{1} +\varepsilon _{2} = a_{0} + a_{1}(a_{0} + a_{1}y_{0} + \varepsilon _{1}) + \varepsilon _{2} = a_{0} + a_{0}a_{1} + (a_{1})^{2}y_{0} + a_{1}\varepsilon _{2}+ \varepsilon _{2}.](/sites/default/files/tex_cache/2ddad7f9ea979cf125639c2d5a08bd1d.png)
Продолжая процесс, найдем :
![y_{3} = a_{0} + a_{1}y_{2} + \varepsilon _{3} = a_{0}[1 + a_{1} + (a_{1})^{2}] + (a_{1})^{3}y_{0} + (a_{1})^{2} \varepsilon _{1} + a_{1} \varepsilon _{2} + \varepsilon _{3}.](/sites/default/files/tex_cache/06499a025fd6ab7c03358d70ca795418.png)
Для итерации с номером получаем решение уравнения (9.7):
Предположим, что начальное значение неизвестно. Тогда при движении "назад" заменим
на итерацию "назад"
В результате получаем:
Предполагая, что и сдвигаясь назад на m периодов
, получаем:
Можно показать, что полученное решение уравнения (9.7) не единственное. Решением уравнения (9.7) при любом
будет и выражение
Выражение (9.9), как и в случае дифференциальных уравнений первого порядка, представляет сумму общего решения однородного разностного уравнения вида (9.2)
и частного решения неоднородного уравнения (9.7).
Теперь рассмотрим случай . Тогда движением "назад" получаем решение в виде
В отличие от предыдущего случая здесь ошибки
не убывают, а накапливаются.
Рассмотрим теперь уравнение второго порядка
Будем искать решение уравнения второго порядка в том же виде
![y_{t}^{одн} = C\alpha ^{t}, C \ne 0.](/sites/default/files/tex_cache/3691b88d5c78fad2327a76edea4700fc.png)
Подставляя это выражение в (9.10), получаем
Разделив обе части (9.11) на, получим характеристическое уравнение
. При решении характеристического уравнения могут возникнуть три случая.
Случай 1. Дискриминант и существуют два действительных различных корня уравнения. Тогда
![y_{t}^{одн} = C_{1}(a_{1})_{t} + C_{2}(a_{2})^{T} .](/sites/default/files/tex_cache/429ba1281581500f805676ebb15af688.png)
Если абсолютное значение хотя бы одного из корней или
превышает единицу, то однородное решение имеет "взрывной" характер при
.
Случай 2. Если , то общее решение имеет вид
Решение носит "взрывной" характер при . Если
, то из-за слагаемого
поведение общего решения не так ясно. В итоге решение стремится к нулю, но вначале может носить "взрывной" характер.
Случай 3. Если , то
и характеристические корни комплексные. В этом случае однородное решение может быть получено в виде
где
- произвольные постоянные, а \theta удовлетворяет соотношению
Тригонометрическая функция создает волнообразное поведение решения однородного уравнения. Частота колебаний определяется параметром
, а амплитуда колебаний - множителем
. Если
, то амплитуда не меняется. Колебания затухают при
и растут при
.
Характеризация условий устойчивости. В случае существования двух корней характеристического уравнения для устойчивости требуется, чтобы оба корня находились внутри промежутка (-1; 1). В случае совпадающих характеристических корней условие устойчивости принимает вид . Для комплексных корней в случае 3
условие устойчивости следующее: -
.