Автор: Пользователь скрыл имя, 06 Февраля 2011 в 15:47, контрольная работа
Постройте график временного ряда.
Постройте автокорреляционную функцию временного ряда. Охарактеризуйте структуру ряда.
Постройте аддитивную и мультипликативную модели временного ряда, учитывая, что трендовая компонента имеет линейный вид.
Оцените качество модели через показатели средней абсолютной ошибки и среднего относительного отклонения.
Рассчитайте прогнозные значения на очередной следующий год.
ИСХОДНЫЕ ДАННЫЕ 3
РЕШЕНИЕ 4
1) Построим график временного ряда 4
2) Построим автокорреляционную функцию временного ряда. 4
3 и 4) Построим аддитивную и мультипликативную модели временного ряда 11
5) Рассчитаем прогнозные значения на очередной следующий год 21
5a) Прогнозирование по аддитивной модели 21
5б) Прогнозирование по мультипликативной модели. 22
БИБЛИОГРАФИЧЕСКИЙ СПИСОК 23
Шаг 4. Определим компоненту Т данной модели. Для этого проведём аналитическое выравнивание ряда (Т+Е) с помощью линейного тренда. Линейная функция имеет вид:
= a + b×t.
Параметры a
и b согласно методу наименьших квадратов
находятся решением системы нормальных
уравнений:
где y - фактические
(эмпирические) значения результативного
признака; x - факторный признак.
Разделив обе части уравнений системы
на n, получим систему нормальных уравнений
в виде:
Решение этой системы относительно
искомых параметров дает следующие
выражения:
Для расчета
параметров вычислим необходимые величины
в Таблице 10
a = 8.515-0,199×10.5 = 6,425
Таким образом, имеем следующий линейный тренд:
T = 0.199×t + 6,425
Подставляя в это уравнение значения t = 1,……..,20, найдём уровни T для каждого момента времени (графа 7 Таблицы 10).
Рисунок 1- Потребление электроэнергии (фактические, выравненные и полученные по аддитивной модели значения уровней ряда)
Шаг 5. Найдём значения уровней ряда, полученные по аддитивной модели. Для этого прибавим к уровням T значения сезонной компоненты для соответствующих кварталов. Графически значения (T+S) представлены на Рисунке 1.
Шаг 6. В соответствии с методикой построения аддитивной модели расчёт ошибки производится по формуле
E=Y-(T+S)
Это абсолютная ошибка. Численные значения абсолютных ошибок приведены в графе 9 Таблицы 10.
Теперь найдём
отношение суммы квадратов случайной
компоненты к общей сумме квадратов
отклонений уровней ряда от его среднего
значения:
Вывод:
Построенная аддитивная модель объясняет
99,53% общей вариации уровней временного
ряда потребления электроэнергии за 20
кварталов исследуемых 5-ти лет и её можно
использовать в прогнозах будущего потребления
электроэнергии.
Построим мультипликативную модель.
Шаг 1.
Таблица 11 - Расчёт оценок сезонной компоненты в мультипликативной модели
№ квартала, t |
Потребление электроэнергии, | Итого за четыре квартала | Скользящая средняя за четыре квартала | Центрированная скользящая средняя | Оценка сезонной компоненты |
1 | 2 | 3 | 4 | 5 | 6 |
1 | 7 | - | - | - | - |
2 | 4,3 | - | - | - | - |
3 | 5 | 25,2 | 6,3 | - | - |
4 | 8,9 | 26,3 | 6,58 | 6,44 | 1,38 |
5 | 8,1 | 27,9 | 6,98 | 6,78 | 1,19 |
6 | 5,9 | 29,9 | 7,48 | 7,23 | 0,82 |
7 | 7 | 32,1 | 8,03 | 7,76 | 0,9 |
8 | 11,1 | 33,1 | 8,28 | 8,16 | 1,36 |
9 | 9,1 | 33,9 | 8,48 | 8,38 | 1,09 |
10 | 6,7 | 34,4 | 8,6 | 8,54 | 0,78 |
11 | 7,5 | 35,3 | 8,83 | 8,72 | 0,86 |
12 | 12 | 36,2 | 9,05 | 8,94 | 1,34 |
13 | 10 | 36,9 | 9,23 | 9,14 | 1,09 |
14 | 7,4 | 37,4 | 9,35 | 9,29 | 0,8 |
15 | 8 | 36,9 | 9,23 | 9,29 | 0,86 |
16 | 11,5 | 37,9 | 9,48 | 9,36 | 1,23 |
17 | 11 | 38,7 | 9,68 | 9,58 | 1,15 |
18 | 8,2 | 39,8 | 9,95 | 9,82 | 0,84 |
19 | 9,1 | 40,8 | 10,2 | 10,08 | 0,9 |
20 | 12,5 | - | - | - | - |
Шаг 2. Найдём оценки сезонной компоненты как частное от деления фактических уровней ряда на центрированные скользящие средние (графа 6 Таблицы 11).
Используем эти оценки для расчёта значений сезонной компоненты S (Таблица 12).
Для этого найдём средние за каждый квартал оценки сезонной компоненты .
Взаимопогашаемость сезонных воздействий в мультипликативной модели выражается в том, что сумма значений сезонной компоненты по всем кварталам должна быть равна числу периодов в цикле. В нашем случае число периодов одного цикла (год) равно 4 (четыре квартала).
Таблица 12 - Расчёт значений сезонной компоненты в мультипликативной модели
Показатели | Год | № квартал, i | |||
I | II | III | IV | ||
1 | 2 | 3 | 4 | 5 | 6 |
1
2 3 4 5 |
-
0,82 0,78 0,8 0,84 |
-
0,9 0,86 0,86 0,9 |
1,38
1,36 1,34 1,23 - |
1,19
1,09 1,09 1,15 - | |
Итого за i-й квартал (за 4 года) | 3,24 | 3,52 | 5,31 | 4,52 | |
Средняя оценка сезонной компоненты для i-го квартала, | 0,81 | 0,88 | 1,33 | 1,13 | |
Скорректированная сезонная компонента, | 0,781 | 0,848 | 1,282 | 1,089 |
Для данной модели имеем: 0,81 + 0,88 + 1,33 + 1,13 = 4,15
Определим корректирующий коэффициент : k = = 0,964
Рассчитаем скорректированные значения сезонной компоненты как произведение между её средней оценкой и корректирующим коэффициентом k :
× k,
где i = 1: 4.
Проверим условие равенства нулю суммы значений сезонной компоненты:
0,781 + 0,848 + 1,282 + 1,089 = 4
Таким образом, получены следующие значения сезонной компоненты:
I квартал: = 0,781;
II квартал: = 0,848;
III квартал: = 1,282;
IV квартал:
= 1,089.
Шаг 3. Разделим каждый уровень исходного ряда на соответствующие значения сезонной компоненты. Тем самым получим величины T×E = Y÷S (графа 4 Таблицы 13), которые содержат только тенденцию и случайную компоненту.
Таблица 13 - Расчёт выравненных значений T и ошибок Е в мультипликативной модели
t
(x) |
T×E=
(y) |
() |
t(T×E),
(xy) |
T | T × S | |
||||
1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 |
1 | 7 | 0,781 | 8,96 | 1 | 8,96 | 6,861 | 5,36 | 1,306 | 1,64 | 2,69 |
2 | 4,3 | 0,848 | 5,07 | 4 | 10,14 | 7,066 | 5,99 | 0,718 | -1,69 | 2,86 |
3 | 5 | 1,282 | 3,9 | 9 | 11,7 | 7,271 | 9,32 | 0,536 | -4,32 | 18,67 |
4 | 8,9 | 1,089 | 8,17 | 16 | 32,68 | 7,476 | 8,14 | 1,093 | 0,76 | 0,58 |
5 | 8,1 | 0,781 | 10,37 | 25 | 51,85 | 7,681 | 6,00 | 1,35 | 2,10 | 4,41 |
6 | 5,9 | 0,848 | 6,96 | 36 | 41,76 | 7,886 | 6,69 | 0,882 | -0,79 | 0,62 |
7 | 7 | 1,282 | 5,46 | 49 | 38,22 | 8,091 | 10,37 | 0,675 | -3,37 | 11,37 |
8 | 11,1 | 1,089 | 10,19 | 64 | 81,52 | 8,296 | 9,03 | 1,229 | 2,07 | 4,27 |
9 | 9,1 | 0,781 | 11,65 | 81 | 104,85 | 8,501 | 6,64 | 1,37 | 2,46 | 6,06 |
10 | 6,7 | 0,848 | 7,9 | 100 | 79 | 8,706 | 7,38 | 0,908 | -0,68 | 0,47 |
11 | 7,5 | 1,282 | 5,85 | 121 | 64,35 | 8,911 | 11,42 | 0,657 | -3,92 | 15,40 |
12 | 12 | 1,089 | 11,02 | 144 | 132,24 | 9,116 | 9,93 | 1,208 | 2,07 | 4,30 |
13 | 10 | 0,781 | 12,8 | 169 | 166,4 | 9,321 | 7,28 | 1,374 | 2,72 | 7,40 |
14 | 7,4 | 0,848 | 8,73 | 196 | 122,22 | 9,526 | 8,08 | 0,916 | -0,68 | 0,46 |
15 | 8 | 1,282 | 6,24 | 225 | 93,6 | 9,731 | 12,48 | 0,641 | -4,48 | 20,03 |
16 | 11,5 | 1,089 | 10,56 | 256 | 168,96 | 9,936 | 10,82 | 1,063 | 0,68 | 0,46 |
17 | 11 | 0,781 | 14,08 | 289 | 239,36 | 10,141 | 7,92 | 1,389 | 3,08 | 9,49 |
18 | 8,2 | 0,848 | 9,67 | 324 | 174,06 | 10,346 | 8,77 | 0,935 | -0,57 | 0,33 |
19 | 9,1 | 1,282 | 7,1 | 361 | 134,9 | 10,551 | 13,53 | 0,673 | -4,43 | 19,59 |
20 | 12,5 | 1,089 | 11,48 | 400 | 229,6 | 10,756 | 11,71 | 1,067 | 0,79 | 0,62 |
∑ 210 | 170.3 | - | 176,16 | 2870 | 1986,37 | - | - | - | - | 130,07 |
=10.5 | - | - | - | - | - | - | - | - | - | - |
- | - | 8.808 | - | - | - | - | - | - | - | |
- | - | - | 143.5 | - | - | - | - | - | - | |
- | - | - | - | 99.32 | - | - | - | - | - |
Шаг 4. Определим компоненту Т данной модели. Для этого проведём аналитическое выравнивание ряда (Т×Е) с помощью линейного тренда. Линейная функция имеет вид:
= a + b×t.
Для расчета
параметров вычислим необходимые величины
в Таблице 13
a = 8.808 - 0,205×10.5 = 6,656
Таким образом, имеем следующий линейный тренд:
T = 0.205×t + 6,656
Подставляя в это уравнение значения t = 1,……..,20, найдём уровни T для каждого момента времени (графа 7 Таблицы 13).
Рисунок
2- Потребление электроэнергии
жителями района (фактические,
выравненные и полученные
по мультипликативной
модели значения уровней
ряда)
Шаг 5. Найдём уровни ряда по мультипликативной модели, умножив уровни T на значения сезонной компоненты для соответствующих кварталов. Графически значения (T×S) представлены на Рисунке 2.
Шаг 6. Расчёт
ошибки в мультипликативной модели производится
по формуле:
Численные значения ошибки приведены в графе 9 Таблицы 13. Для того чтобы сравнить мультипликативную модель и другие модели временного ряда используем метод суммы квадратов абсолютных ошибок.
Абсолютные ошибки
в мультипликативной модели определяются
так:
В данной модели сумма квадратов абсолютных ошибок равна 130,07
Теперь найдём
отношение суммы квадратов случайной
компоненты к общей сумме квадратов
отклонений уровней ряда от его среднего
значения:
Вывод: Построенная мультипликативная модель объясняет 99,5% общей вариации уровней временного ряда потребления электроэнергии за 20 кварталов исследуемых 5-ти лет и её можно использовать в прогнозах будущего потребления электроэнергии.
5a) Прогнозирование по аддитивной модели.
Прогнозное значение уровня временного ряда в соответствии с соотношением Y=T+S+E есть сумма трендовой и сезонной компонент.
Объём электроэнергии, потреблённой в течение следующего года, т.е шестого года, рассчитывается как сумма объёмов потребления электроэнергии в 4-х кварталах шестого года, соответственно . Для определения трендовой компоненты воспользуемся уравнением тренда аддитивной модели: T = 0.199×t + 6,425