Лабораторная работа № 6. РЯДЫ ДИНАМИКИ. АНАЛИТИЧЕСКАЯ МОДЕЛЬ СЕЗОННЫХ И СЛУЧАЙНЫХ ИЗМЕНЕНИЙ. ИНДЕКСЫ СЕЗОННОСТИ.
Цель работы: Выявление аналитической закономерности выровненных рядов динамики, учитывающих сезонную и случайные изменения во времени признаков ряда. Рассчитать индексы сезонности.
Задание. Используя данные о выпуске продукции предприятием за двенадцать месяцев (табл. 1) найти аналитические тригонометрические модели динамического ряда. Соответствующие модели представить графиками. Сравнить полученные результаты. Вычислить индексы сезонности. Сделать выводы.
Условие. Имеются сведения о произведенной продукции предприятием за год (табл. 1).
Таблица 1
Выполнение задания.
Аналитическая модель основного тренда, сезонной и случайной компонент динамического ряда
Анализ колеблемости динамических рядов наряду с выделением случайных колебаний требует изучения периодических колебаний. Изучение периодических («сезонных») колебаний необходимо с целью исключения их влияния на общую динамику для выявления «случайной» колеблемости.
Выделение сезонного периода можно выполнить на основе построения аналитической модели проявления сезонных колебаний. Построение аналитической модели выявляет основной закон колеблемости данного временного ряда в связи с переходом от месяца к месяцу и дает среднюю характеристику внутригодичных колебаний.
При исследовании явлений периодического типа в качестве аналитической формы развития во времени принимается метод Фурье, разложения в ряд по фундаментальной системе тригонометрических функций
:
EMBED Equation.3 = а0 + EMBED Equation.3 (ak cos kt+ bk sin kt). (1)
В этой зависимости величина k определяет гармонику ряда Фурье и может быть взята с разной степенью точности (чаще всего от 1 до 4). Для отыскания параметров нового ряда (1) используется метод наименьших квадратов, т.е.
EMBED Equation.3 = min . (2)
Справедливы следующие формулы для вычисления параметров:
a0 = EMBED Equation.3 ;
ak = EMBED Equation.3 ; (3)
bk = EMBED Equation.3 ;
Для изучения сезонных колебаний на протяжении года необходимо взять n=12 (по числу месяцев в году). Тогда, представляя периоды как части длины окружности, ряд динамики можно записать в следующем виде:
Таблица 2
Для вычисления синусов и косинусов разных гармоник лучше всего пользоваться табл. 3.
Таблица 3
Для k =1 уравнение (1) примет вид:
EMBED Equation.3 = а0 + a1 cos t+ b1 sin t , (4),
в котором параметры а0, a1 и b1 будут найдены из соотношений:
a0 = EMBED Equation.3 ;
a1 = EMBED Equation.3 ; (5)
b1 = EMBED Equation.3 ;
Для k =2 уравнение (1) будет таким:
EMBED Equation.3 = а0 + a1 cos t+ b1 sin t+ a2 cos 2t+ b2 sin 2t , (6),
в котором параметры a2 и b2 определяются из соотношений:
a2 = EMBED Equation.3 ; (7)
b2 = EMBED Equation.3 ;
Применим расчетную схему (4)-(7) к исходной задаче. В результате получим аналитическую модель динамического ряда для четырех гармоник k =1 , k =2 , k =3 и k =4 . Полученные данные представлены в таблице 4, выполненной в среде Ехсеl.
Из этой таблицы найдем коэффициенты ряда Фурье
a0 = EMBED Equation.3 = 63,583 ; a1 = EMBED Equation.3 = ?21,549 ; b1 = EMBED Equation.3 = ?29,910 ;
a2 = EMBED Equation.3 = 0,75 ; b2 = EMBED Equation.3 = 8,227 ;
a3 = EMBED Equation.3 = ?7,167 ; b3 = EMBED Equation.3 = ?4,0 ;
a4 = EMBED Equation.3 = 3,917 ; b4 = EMBED Equation.3 = ?8,516 .
Запишем соответствующие приближения, подставляя найденные коэффициенты в ряд Фурье (1), ограничиваясь числами k =1 , k =2 , k =3 и k =4 , получим следующие четыре модели исходного динамического ряда:
EMBED Equation.3 = 63,583 ? 21,549 cos t ? 29,910 sin t ,
EMBED Equation.3 = 63,583 ? 21,549 cos t ? 29,910 sin t + 0,75 cos 2t+ 8,227 sin 2t , (8)
EMBED Equation.3 = 63,583 ? 21,549 cos t ? 29,910 sin t + 0,75 cos 2t+ 8,227 sin 2t ? 7,167 cos 3t ? 4 sin 3t ,
EMBED Equation.3 = 63,583 ? 21,549 cos t ? 29,910 sin t + 0,75 cos 2t+ 8,227 sin 2t ? 7,167 cos 3t ? 4 sin 3t + 3,917 cos 4t ? 8,516 sin 4t ,
Таблица 4
В последней строке табл. 4 (гр. 2, 5, 8, 11, 14) представлены фактические и расчетные суммарные уровни выпуска продукции, по которым наглядно видно, что модель достаточно точно отражает эмпирические признаки ряда динамики.
Графическое представление аналитических моделей.
Сопоставим выровненные уровни динамического ряда по первой, второй, третьей и четвертой гармониках (см. табл. 4, гр. 5, 8, 11 и 14) графически (рис. 1).
EMBED Excel.Chart.8 \s Рис. 1. График объемов выпуска продукции предприятием за 12 месяцев года. Эмпирический ряд 1, первая гармоника выровненного ряда гр. 5 табл.4 – ряд 2, вторая гармоника выровненного ряда гр. 8 табл.4 – ряд 3, третья гармоника выровненного ряда гр. 11 табл.4 – ряд 4, четвертая гармоника выровненного ряда гр. 14 табл.4 – ряд 5
Анализ результатов приводит к выводу о достаточности использования для выравнивания только первой гармоники. Её значения являются усредненными по сравнению со всеми перечисленными.
3. Индексы сезонности.
Глубину сезонных колебаний месячных данных измеряют индексами сезонности, которые представляют собой отношение средних из фактических уровней одноименных месяцев за рассматриваемый период к средней из выравненных данных по тем же месяцам
yсез. = EMBED Equation.3 , (9)
где EMBED Equation.3 - средняя из фактических уровней i-того месяца за весь рассматриваемый период; EMBED Equation.3 - средний из выравненных уровней i -того месяца, полученные либо применением 12-месячной скользящей средней, либо аналитическим выравниваем.
Вычисление индексов сезонности проиллюстрируем по данным полученной таблицы 4. Будем полагать, что за 5-летний период были рассчитаны средние из фактических уровней выпуска продукции предприятием за одноименные месяцы, затем производилось выравнивание ряда данных по линейному тренду и рассчитывались средние из выравненных данных по каждому месяцу (результаты расчетов приведены в табл. 5, гр. 3 и 4).
Таблица 5
Обобщающим показателем изменений динамического ряда из-за сезонного характера производства является среднее квадратическое отклонение индексов сезонности
? сез.= EMBED Equation.3 .
Расчет ? сез основан па результатах, представленных в гр. 6 табл. 5, и составляет 31,78 %.
Сравнение средних квадратических отклонений, вычисленных за разные периоды, показывает сдвиги в сезонности. Уменьшение ? сез свидетельствует об уменьшении влияния сезонности на динамику анализируемого показателя.
Выводы.
Основные выводы содержатся в каждом пункте выполнения задания.
Варианты заданий. Варианты указаны римскими цифрами в лабораторной № 4. Выбрать двенадцать (и более) признаков динамического ряда из предложенного преподавателем варианта.