5.4.4. Прогнозирование процессов авторегрессии
Рассмотрим процесс порядка
Эвентуальная
прогнозирующая функция в этом случае есть решение уравнения . Она применима ко всем
упреждениям и проходит через последние известных значений ряда. Например,
модель ряда цен акций IBM (ряда ) очень близка к
так
что
Лучший
прогноз на будущее оказывается очень близким к текущему значению цены. Весовая
функция для сосредоточена
в точке , и не происходит
никакого осреднения по прошлому.
Стационарные модели авторегрессии. Процесс порядка , где — стационарный
оператор и с
, будет в общем
случае иметь прогнозирующую функцию, являющуюся совокупностью экспонент и
затухающих синусоид.
В
частности, при =
1 модель порядка (1, 0, 0)
имеет
прогнозирующую функцию, которая при любых является решением уравнения , Отсюда
(5.4.15)
Кроме
того, , так
что
и
Следовательно,
прогноз с минимальной среднеквадратичной ошибкой предсказывает, что текущее
отклонение от среднего должно экспоненциально спадать до нуля. На рис. 5.9, а
показан временной ряд, генерируемый процессом , с его прогнозирующей
функцией для момента . Поведение этой функции целиком
определяется только отклонением . Аналогично минимальная
среднеквадратичная ошибка прогноза для процесса авторегрессии второго порядка
такова, что текущее отклонение от среднего убывает до нуля как затухающая
синусоида
или сумма двух экспонент. На рис. 5.9, б показан временной ряд, генерируемый
процессом ,
и его прогноз для момента . Поведение прогнозирующей функции для
момента целиком
определяется последними двумя отклонениями и .
Рис. 5.9. Прогнозирующие функции процессов
авторегрессии первого и второго порядка:
а — выборка из процесса авторегрессии
первого порядка и
прогнозирующая функция для момента = 14. б — выборка из процесса
авгорегрессии второго порядка и прогнозирующая функция на
момент =14.
1 — ход прогноза, определяемый только по текущему отклонению, 2 — ход
прогноза, определяемый по текущему и предпоследнему отклонениям.
Функция дисперсии для прогноза процесса
(1,0,0).
Чтобы полнее проиллюстрировать использование (5.1.16), выведем функцию
дисперсии для процесса авторегрессии первого порядка. Так как модель в момент может быть
представлена в виде
то из (5.4.15)
следует, что
Отсюда
(5.4.16)
Мы
видим, что для этого стационарного процесса при стремлении к бесконечности дисперсия
растет до постоянного значения , равного среднеквадратичному отклонению
процесса относительно окончательного прогноза . Это поведение существенно иное, чем у
нестационарных моделей, у которых функция дисперсии прогноза при больших
упреждениях неограниченно растет.
Рис. 5.10. Прогнозирующая функция
для двух процессов (1,1,0):
а — прогнозы процесса (1,1,0) (1
— 0,8) , б — прогнозы
процесса (1,1,0) (1 — 0,8) (— 0,2) = .
Нестационарные модели авторегрессии
порядка . Для модели
затухать
до своего среднего значения при прогнозировании на несколько шагов вперед будет
уже -я
разность процесса. Среднее будет, как правило,
предполагаться нулевым, если только не появится каких-либо свидетельств в
пользу противного. Если это необходимо, можно, как указывалось в гл. 4, ввести ненулевое
среднее, заменив в модели на . Например,
рассмотрим модель
(5.4.17)
После
замены на , переходя к условным
математическим ожиданиям при известном на момент прошлом, легко получим (ср. с (5.4.15) и
следующими за ним)
Эта
формула показывает, что прогнозируемая разность экспоненциально убывает
от начального значения до своего среднего значения . Суммируя такие
выражения для ,
изменяющегося от 1 до получаем прогнозирующую функцию
которая
асимптотически приближается к прямой линии
На
рис, 5.10 показаны прогнозы для двух случаев и . Мы увидим в гл. 7, что модель (5.4.17)
с и
хорошо
описывает ряд ;
прогнозы, основанные на этой модели, уже демонстрировались на рис. 5.1 и 5.2.
Рассмотрим теперь прогнозирование некоторых важных смешанных моделей.