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