Статистическое изучение сезонных колебаний
Под сезонными колебаниями понимаются более или менее устойчивые внутригодовые колебания уровней развития социально-экономических явлений.
Сезонная колеблемость (внутригодовая цикличность) показателей рынка выявляется с помощью индексов сезонности, методов механического сглаживания, гармоникой Фурье. Гармоника Фурье имеет следующие аналитическое выражение:

где k – номер гармоники, определяющей степень точности (адекватности) модели (обычно k берётся в пределах от 1 до 4);
t – время;
а0, аk, bk – параметры уравнения.
При 
При 
Параметры уравнения (
) определяются методом наименьших квадратов:

При анализе внутригодовой динамики n = 12 – по числу месяцев в году. Представляя месячные периоды как части окружности, ряд внутригодовой динамики имеет следующий вид:
| Периоды, t |
|
|
|
|
|
|
|
|
|
|
| |
| Уровни, y | y1 | y2 | y3 | y4 | y5 | y6 | y7 | y8 | y9 | y10 | y11 | y12 |
Для расчёта гармоники используется вспомогательная таблица. Для первой гармоники она имеет следующий вид:
| Месяц | Условное обозначение месяца, t | Уровни ряда, y | cos t | sin t | y · cos t | y · sin t |
|
| . . и т.д. | 1/6 π 1/3 π | y1 y2 y3 | |||||
|
Решение типовых задач
Пример 6.1
Розничный товарооборот в регионе в 2008 – 2012 гг. характеризуется
следующими данными, млн руб:
| 2008 г. | 2009 г. | 2010 г. | 2011 г. | 2012 г. |
| 8 900 | 8 670 | 8 930 | 8 990 |
Для анализа ряда динамики определим:
1) показатели, характеризующие динамику товарооборота по годам;
2) средние показатели динамики за исследуемый период.
Решение
Показатели динамики уровней временного ряда с постоянной (базисные) и переменной (цепные) базой сравнения:
1) абсолютный прирост –
; 
2) коэффициент роста –
; 
3) темп роста –
; 
4) темп прироста –
;
.
5) абсолютное значение 1 % прироста –
; 
6) пункты роста –
.
Результаты расчётов приведены в таблице:
| Год | 2008 г. | 2009 г. | 2010 г. | 2011 г. | 2012 г. | |
| Розничный товарооборот, млн руб. | ||||||
| Абсолютный прирост, млн руб. | Цепной | – | –230 | |||
| Базисный | – | –230 | ||||
| Коэффициент роста | Цепной | – | 0,974 | 1,029 | 1,006 | 1,012 |
| Базисный | – | 0,974 | 1,003 | 1,010 | 1,022 | |
| Темп роста, млн руб. | Цепной | – | 97,4 | 102,9 | 100,6 | 101,2 |
| Базисный | – | 97,4 | 100,3 | 101,0 | 102,2 | |
| Темп прироста, млн руб. | Цепной | – | –2,8 | 2,9 | 0,6 | 1,2 |
| Базисный | – | –2,8 | 0,3 | 1,0 | 2,2 | |
| Абсолютное значение 1% прироста, млн руб. | Цепной | – | 82,14 | 89,66 | 100,00 | 91,67 |
| Базисный | – | 82,14 | 100,00 | 90,00 | 90,91 | |
| Пункты роста | – | – | 2,9 | 0,7 | 1,2 |
Средние характеристики ряда динамики:
1) средний уровень интервального ряда рассчитывается по формуле средней арифметической простой:
млн руб.
2) средний абсолютный прирост (
):
,
млн руб.
3) средний коэффициент роста (
):
;
;
4) средний темп роста:
;
;
5) средний темп прироста:
;
;
6) средняя величина абсолютного значения 1% прироста:
;
млн руб..
Сделаем выводы:
1. Так как темпы роста начиная с 2010 г. имеют значение больше100, то это говорит о тенденции роста объёма розничного товарооборота в регионе.
2. Средний годовой уровень объёма товарооборота в исследуемом периоде составил 8 918 млн руб.
3. В среднем ежегодно товарооборот увеличился на 50 млн руб., или на 0, 3 %
4. Средняя величина абсолютного значения 1% прироста показывает, что при увеличении производства электроэнергии на 1% в абсолютном выражении этот прирост составит 166,7 млн руб.
Пример 6.2
Выполним интервальный прогноз объёма реализации продукции А на 2008 г. с вероятностью 0,99 на основе следующих отчётных данных торгового предприятия:
| Год | ||||||||
| Объём реализа-ции, т |
Для определения формы тренда и расчёта его параметров составим таблицу:
| Год | Объём реализации, т y | Первые разности | t | t2 | y · t | yt | (y – yt)2 |
| - | -7 | -532 | 76,53 | 0,2809 | |||
| -5 | -435 | 87,45 | 0,2025 | ||||
| -3 | -297 | 98,37 | 0,3969 | ||||
| -1 | -109 | 109,25 | 0,0625 | ||||
| +1 | 120,21 | 0,0441 | |||||
| +3 | 131,13 | 0,7569 | |||||
| +5 | 142,05 | 0,0025 | |||||
| +7 | 1 064 | 52,97 | 0,9409 | ||||
| Итого | - | 917,96 | 2,6872 |
Первые разности приблизительно равны между собой, что позволяет в качестве модели тренда выбрать уравнение прямой:
yt = a0 + a1 · t.
Найдем параметры уравнения:

Подставим найденные значения параметров в уравнение:
yt = 114,75 + 5,46 t.
Рассчитаем теоретические уровни ряда динамики:
2000 г.: yt = 114,75 + 5,46 · (-7) = 76,53 т.
2001 г.: yt = 114,75 + 5,46 · (-5) = 87,45 т. и т.д.
Выполним точечный прогноз на 2008 г.:
yt+L = 114,75 + 5,46 · 9 = 163,89 т.
Для получения интервального прогноза воспользуемся формулой

Коэффициент доверия
находим по таблице распределения Стьюдента (Приложение Б). При вероятности P = 0,99, (уровень значимости α = 1,0), числе степеней свободы k = n – 2 = 6 коэффициент доверия tα = 3,7.
Рассчитаем стандартную ошибку аппроксимации:

Расчёты квадратов отклонений теоретических уровней ряда от эмпирических приведены в таблице.
Выполним интервальный прогноз

Пример 6.3.Имеются следующие данные о реализации фруктов в магазинах города:
| Месяц | Продано фруктов, т | Расчётная графа | |||||
| Условное обозначение месяца, t | cos t | sin t | y · cos t | y · sin t |
| ||
| 0,0 | 68,00 | 75,03 | |||||
| 1/6 π | 0,866 | 0,5 | 47,63 | 27,50 | 50,11 | ||
| 1/3 π | 0,5 | 0,866 | 25,00 | 43,30 | 34,86 | ||
| 1/2 π | 0,0 | 1,0 | 33,35 | ||||
| 1/3 π | -0,5 | 0,866 | -27,00 | 46,764 | 46,00 | ||
| 5/6 π | -0,866 | 0,5 | -56,29 | 32,50 | 69,41 | ||
| Π | -1,0 | 0,0 | -90,00 | 97,31 | |||
| 7/6 π | -0,866 | -0,5 | -103,92 | -60,00 | 122,23 | ||
| 4/3 π | -0,5 | -0,866 | -72,50 | -125,57 | 137,48 | ||
| 3/2 π | 0,0 | -1,0 | -130,00 | 138,99 | |||
| 5/3 π | 0,5 | -0,866 | 60,00 | -103,92 | 126,34 | ||
| 11/6 π | 0,866 | -0,5 | 82,27 | -47,50 | 102,93 | ||
| Итого | 1 034 | - | - | - | -66,81 | -316,93 | 1 034,04 |
Построим модель сезонных колебаний в объёме реализации фруктов, используя первую гармонику ряда Фурье:
Первая гармоника имеет вид:

Параметры уравнения определяются по формулам:

Для определения параметров уравнения введём в таблицу расчётные графы.
Рассчитаем параметры уравнения:

Модель сезонной волны объёма реализации фруктов примет вид:
= 86,17 – 11,14 cos t – 52,82 sin t.
Определим расчётные (теоретические) уровни для каждого месяца:
январь:
= 86,17 – 11,14 · 1,0 – 52,82 · 0 = 75,03 т.
февраль:
= 86,17 – 11,14 · 0,866 – 52,82 · 0,5 = 50,11 т. и т.д.
Результаты
расчётов занесём в таблицу.