|
|
Моделирование сезонных колебанийПростейший подход к моделированию сезонных колебаний – это расчет значений сезонной компоненты методом скользящей средней и построение аддитивной или мультипликативной модели временного ряда. Общий вид аддитивной модели следующий:
Эта модель предполагает, что каждый уровень временного ряда может быть представлен как сумма трендовой ( Общий вид мультипликативной модели выглядит так:
Эта модель предполагает, что каждый уровень временного ряда может быть представлен как произведение трендовой ( Выбор одной из двух моделей осуществляется на основе анализа структуры сезонных колебаний. Если амплитуда колебаний приблизительно постоянна, строят аддитивную модель временного ряда, в которой значения сезонной компоненты предполагаются постоянными для различных циклов. Если амплитуда сезонных колебаний возрастает или уменьшается, строят мультипликативную модель временного ряда, которая ставит уровни ряда в зависимость от значений сезонной компоненты. Построение аддитивной и мультипликативной моделей сводится к расчету значений Процесс построения модели включает в себя следующие шаги. 1) Выравнивание исходного ряда методом скользящей средней. 2) Расчет значений сезонной компоненты 3) Устранение сезонной компоненты из исходных уровней ряда и получение выровненных данных ( 4) Аналитическое выравнивание уровней ( 5) Расчет полученных по модели значений ( 6) Расчет абсолютных и/или относительных ошибок. Если полученные значения ошибок не содержат автокорреляции, ими можно заменить исходные уровни ряда и в дальнейшем использовать временной ряд ошибок Методику построения каждой из моделей рассмотрим на примерах. Пример. Построение аддитивной модели временного ряда. Обратимся к данным об объеме правонарушений на таможне за четыре года, представленным в табл. 4.1. Было показано, что данный временной ряд содержит сезонные колебания периодичностью 4, т.к. количество правонарушений в первый-второй кварталы ниже, чем в третий-четвертый. Рассчитаем компоненты аддитивной модели временного ряда. Шаг 1. Проведем выравнивание исходных уровней ряда методом скользящей средней. Для этого: 1.1. Просуммируем уровни ряда последовательно за каждые четыре квартала со сдвигом на один момент времени и определим условные годовые объемы потребления электроэнергии (гр. 3 табл. 4.5). 1.2. Разделив полученные суммы на 4, найдем скользящие средние (гр. 4 табл. 4.5). Полученные таким образом выровненные значения уже не содержат сезонной компоненты. 1.3. Приведем эти значения в соответствие с фактическими моментами времени, для чего найдем средние значения из двух последовательных скользящих средних – центрированные скользящие средние (гр. 5 табл. 4.5). Таблица 4.5
Шаг 2. Найдем оценки сезонной компоненты как разность между фактическими уровнями ряда и центрированными скользящими средними (гр. 6 табл. 4.5). Используем эти оценки для расчета значений сезонной компоненты Таблица 4.6
Для данной модели имеем:
Корректирующий коэффициент: Рассчитываем скорректированные значения сезонной компоненты ( Проверим равенство нулю суммы значений сезонной компоненты:
Шаг 3. Исключим влияние сезонной компоненты, вычитая ее значение из каждого уровня исходного временного ряда. Получим величины Таблица 4.7
Шаг 4. Определим компоненту
Подставляя в это уравнение значения Шаг 5. Найдем значения уровней ряда, полученные по аддитивной модели. Для этого прибавим к уровням На одном графике отложим фактические значения уровней временного ряда и теоретические, полученные по аддитивной модели.
Рис. 4.6. Для оценки качества построенной модели применим сумму квадратов полученных абсолютных ошибок.
Следовательно, можно сказать, что аддитивная модель объясняет 97% общей вариации уровней временного ряда количества правонарушений по кварталам за 4 года. Шаг 6. Прогнозирование по аддитивной модели. Предположим, что по нашему примеру необходимо дать прогноз об общем объеме правонарушений на I и II кварталы 2003 года. Прогнозное значение
Получим
Значения сезонных компонент за соответствующие кварталы равны:
Т.е. в первые два квартала 2003 г. следовало ожидать порядка 395 и 422 правонарушений соответственно. Построение мультипликативной модели рассмотрим на данных предыдущего примера. Шаг 1. Методика, применяемая на этом шаге, полностью совпадает с методикой построения аддитивной модели. Таблица 4.8
Шаг 2. Найдем оценки сезонной компоненты как частное от деления фактических уровней ряда на центрированные скользящие средние (гр. 6 табл. 4.8). Эти оценки используются для расчета сезонной компоненты Таблица 4.9
Имеем
Определяем корректирующий коэффициент:
Скорректированные значения сезонной компоненты Проверяем условие равенство 4 суммы значений сезонной компоненты:
Шаг 3. Разделим каждый уровень исходного ряда на соответствующие значения сезонной компоненты. В результате получим величины Таблица 4.10
Шаг 4. Определим компоненту
Подставляя в это уравнение значения Шаг 5. Найдем уровни ряда, умножив значения
Рис. 4.7. Расчет ошибки в мультипликативной модели производится по формуле:
Для сравнения мультипликативной модели и других моделей временного ряда можно, по аналогии с аддитивной моделью, использовать сумму квадратов абсолютных ошибок
Сравнивая показатели детерминации аддитивной и мультипликативной моделей, делаем вывод, что они примерно одинаково аппроксимируют исходные данные. Шаг 6. Прогнозирование по мультипликативной модели. Если предположить, что по нашему примеру необходимо дать прогноз об общем объеме правонарушений на I и II кварталы 2003 года, прогнозное значение
Получим
Значения сезонных компонент за соответствующие кварталы равны:
Т.е. в первые два квартала 2003 г. следовало ожидать порядка 409 и 436 правонарушений соответственно. Таким образом, аддитивная и мультипликативная модели дают примерно одинаковый результат по прогнозу. |
|
. (4.3)
), сезонной ( 
) и случайной ( 
) компонент.
. (4.4)
) в аддитивной или ( 
) в мультипликативной модели.
. В моделях с сезонной компонентой обычно предполагается, что сезонные воздействия за период взаимопогашаются. В аддитивной модели это выражается в том, что сумма значений сезонной компоненты по всем кварталам должна быть равна нулю.
.
.
) и заносим полученные данные в таблицу 4.6.
.
(гр. 4 табл. 4.7). Эти значения рассчитываются за каждый момент времени и содержат только тенденцию и случайную компоненту.
.
, найдем уровни 
.
уровня временного ряда в аддитивной модели есть сумма трендовой и сезонной компонент. Для определения трендовой компоненты воспользуемся уравнением тренда
;
.
и 
. Таким образом,
;
.
.
.
.
.
(гр. 4 табл. 4.10), которые содержат только тенденцию и случайную компоненту.
.
.
:
.
;
.
и 
. Таким образом
;
.