Одним из методов моделирования периодических колебаний является ряд Фурье. Построение ряда Фурье для временного ряда зависит от присутствия или отсутствия тенденции во ряду динамики. Если временной ряд стационарен, т.е. тенденция отсутствует, методика построения ряда Фурье применяется непосредственно к уровням временного ряда; если ВР не стационарен, т.е. тенденция присутствует, методика построения ряда Фурье к отклонениям от тенденции.
Теоретически стационарный временной ряд конечной длины n с периодическими колебаниями можно представить, как
, (1)
где ω – угловая частота, измеряемая в радианах в единицу времени (ω=2πf, 0 ≤ ω ≤ 2π);
N= n /2;
.
Параметры уравнения (1) можно оценить МНК.
При прогнозировании по временному ряду без тенденции прогноз выполняется по ряду Фурье. При прогнозировании по временному ряду с тенденцией выполняется суммарный прогноз, т.е. строится прогноз по уравнению тренда и к нему прибавляется прогноз по ряду Фурье отклонений от тренда.
Построим модель динамики производства яиц в Волгоградской области на основе данных за период 2016-2019 гг. (млн. шт.).
Анализ графика ВР (рис. 1) показывает, что данный ряд содержит периодические колебания, но не имеет трендовой компоненты (т.е. ряд стационарен). Уровни ряда колеблются вокруг среднего значения
Рисунок 1 – График временного ряда
Анализ коррелограммы подтверждает вывод о наличии сезонных колебаний периодичностью в четыре квартала.
Построим ряд Фурье с четырьмя гармониками.
где t - принимает значения от 0 с увеличением на 2π/n.
Расчеты для определения параметров ряда Фурье выполним в Excel, используя инструмент анализа данных Регрессия.
Таблица 1 - Расчет параметров по ряду Фурье
№ |
yt |
t |
cos t |
sin t |
cos 2t |
sin 2t |
cos 3t |
sin 3t |
cos 4t |
sin 4t |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
I.2016 |
1 |
145,9 |
0,0000 |
1,0000 |
0,0000 |
1,0000 |
0,0000 |
1,0000 |
0,0000 |
1,0000 |
II.2016 |
2 |
231,8 |
0,3927 |
0,9239 |
0,3827 |
0,7071 |
0,7071 |
0,3827 |
0,9239 |
0,0000 |
III.2016 |
3 |
254,3 |
0,7854 |
0,7071 |
0,7071 |
0,0000 |
1,0000 |
-0,7071 |
0,7071 |
-1,0000 |
IV.2016 |
4 |
164,7 |
1,1781 |
0,3827 |
0,9239 |
-0,7071 |
0,7071 |
-0,9239 |
-0,3827 |
0,0000 |
I.2017 |
5 |
155,4 |
1,5708 |
0,0000 |
1,0000 |
-1,0000 |
0,0000 |
0,0000 |
-1,0000 |
1,0000 |
II.2017 |
6 |
214 |
1,9635 |
-0,3827 |
0,9239 |
-0,7071 |
-0,7071 |
0,9239 |
-0,3827 |
0,0000 |
III.2017 |
7 |
218,7 |
2,3562 |
-0,7071 |
0,7071 |
0,0000 |
-1,0000 |
0,7071 |
0,7071 |
-1,0000 |
IV.2017 |
8 |
171,8 |
2,7489 |
-0,9239 |
0,3827 |
0,7071 |
-0,7071 |
-0,3827 |
0,9239 |
0,0000 |
I.2018 |
9 |
155,4 |
3,1416 |
-1,0000 |
0,0000 |
1,0000 |
0,0000 |
-1,0000 |
0,0000 |
1,0000 |
II.2018 |
10 |
235,7 |
3,5343 |
-0,9239 |
-0,3827 |
0,7071 |
0,7071 |
-0,3827 |
-0,9239 |
0,0000 |
III.2018 |
11 |
247,2 |
3,9270 |
-0,7071 |
-0,7071 |
0,0000 |
1,0000 |
0,7071 |
-0,7071 |
-1,0000 |
IV.2018 |
12 |
162,9 |
4,3197 |
-0,3827 |
-0,9239 |
-0,7071 |
0,7071 |
0,9239 |
0,3827 |
0,0000 |
I.2019 |
13 |
156,6 |
4,7124 |
0,0000 |
-1,0000 |
-1,0000 |
0,0000 |
0,0000 |
1,0000 |
1,0000 |
II.2019 |
14 |
230,6 |
5,1051 |
0,3827 |
-0,9239 |
-0,7071 |
-0,7071 |
-0,9239 |
0,3827 |
0,0000 |
III.2019 |
15 |
241,1 |
5,4978 |
0,7071 |
-0,7071 |
0,0000 |
-1,0000 |
-0,7071 |
-0,7071 |
-1,0000 |
IV.2019 |
16 |
182,7 |
5,8905 |
0,9239 |
-0,3827 |
0,7071 |
-0,7071 |
0,3827 |
-0,9239 |
0,0000 |
Сумма |
136 |
3169 |
6,2832 |
0,00 |
0,00 |
0,00 |
0,00 |
0,00 |
0,00 |
0,00 |
Результат применения инструмента Регрессия представлен на рис. 2.
Рисунок 2 - Результат применения инструмента Регрессия
Модель с четырьмя гармониками хорошо описывает исходный ВР, однако параметры модели для 1, 2 и 3 гармоник по t-критерию Стьюдента статистически незначимы.
Модель рассматриваемого временного ряда представим в виде ряда Фурье с одной (четвертой) гармоникой:
. (2)
На рис. 3 представлен поквартальный график фактического и расчетного объема производства яиц в Волгоградской области за 2016-2019 гг.
Рисунок 3 – Динамика объема производства продукции
Список используемых источников:
Заяц, О. А. Эконометрика: Учебное пособие / О. А. Заяц; О. А. Заяц; ФГБОУ ВО Волгоградский ГАУ. – Волгоград: Волгоградский государственный аграрный университет, 2021. – 140 с.
Афанасьев, В.Н. Анализ временных рядов и прогнозирование: учебник / В.Н. Афанасьев, М.М. Юзбашев. - 2-е изд., перераб. и доп. – М.: Финансы и статистика, 2010. – 320 с. - Текст: непосредственный.
Федеральная служба государственной статистики. Краткосрочные экономические индикаторы по Волгоградской области. – Режим доступа: https://gks.ru/region/ind1118/Main.htm