Построение, анализ и применение для прогнозирования тренд-сезонной модели временного ряда

Предварительный анализ заданного временного ряда на предмет наличия тренда. Обоснование наличия сезонности по графическому представлению одноименных элементов ряда разных лет. Применение модели для прогноза. Выбор типа остатков и корректировка модели.

Рубрика Программирование, компьютеры и кибернетика
Вид контрольная работа
Язык русский
Дата добавления 12.09.2011
Размер файла 218,8 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ УКРАИНЫ

Факультет экономики и менеджмента

Кафедра Экономической кибернетики

Контрольная работа

По дисциплине

Методы прогнозирования

Тема работы

Построение, анализ и применение для прогнозирования тренд-сезонной модели временного ряда

Симферополь, 2010г.

Задание

1) Предварительный анализ заданного временного ряда на предмет наличия тренда:

a) графическое представление ряда;

b) укрупнение;

c) сглаживание.

2) Обоснование выбора типа модели основного тренда по свойствам кривых роста. Экономический смысл параметров выбранной модели.

3) Обоснование наличия сезонности по графическому представлению одноименных элементов ряда разных лет. Рекомендуется предварительно исключить из исходного ряда основную тенденцию развития (тренд, укрупненные или сглаженные значения) вычитанием или делением. Описать характер сезонных колебаний и по возможности объяснить причины.

4) Построение модели декомпозиции временного ряда. Обоснование формы модели: аддитивная, мультипликативная или смешанная.

5) Анализ точности модели.

6) Проверка адекватности моделей данному временному ряду. Анализ свойств остаточной компоненты.

7) Применение модели для прогноза (точечный и интервальный прогноз на два периода сезонности с заданной доверительной вероятностью Р = 0,95).

8) Верификация прогноза (сопоставление с контрольными данными). Выводы о возможности дальнейшего применения модели.

9) Возможно, изменение модели.

Исходные данные

Житлові будівлі, тис. м2 загальної площі

I кв

 

1234,7

I-II кв

 

2709,4

I-III кв

 

4280,0

за год

 

7566,3

I кв

 

1453,3

I-II кв

 

2893,6

I-III кв

 

4524,0

за год

 

7816,0

I кв

 

1690,9

I-II кв

 

3188,7

I-III кв

 

5046,6

за год

 

8628,4

I кв

 

1761,1

I-II кв

 

3947,6

I-III кв

 

6130,3

за год

 

10243,7

I кв

 

1897,2

I-II кв

 

4044,6

I-III кв

 

6483,4

за год

 

10495,6

Предварительный анализ заданного временного ряда на предмет наличия тренда и сезонности

Тренд - это длительная тенденция изменения экономических показателей, т.е. изменение, определяющее общее направление развития, основную тенденцию временных рядов. Тренд характеризует основные закономерности движения во времени, в некоторой мере свободные от случайных воздействий. При разработке моделей прогнозирования тренд оказывается основной составляющей прогнозируемого временного ряда, на которую накладываются другие составляющие. Результат при этом связывается исключительно с ходом времени. Предполагается, что через время можно выразить влияние всех основных факторов. Если тренд является монотонным (устойчиво возрастает или устойчиво убывает), то анализировать такой ряд обычно нетрудно. Если временные ряды содержат значительную ошибку, то первым шагом выделения тренда является сглаживание.

Представляем исходные значения на точечном графике

Выполним сглаживание

Сглаживание всегда включает некоторый способ локального усреднения данных, при котором несистематические компоненты взаимно погашают друг друга. Самый общий метод сглаживания - скользящее среднее, в котором каждый член ряда заменяется простым или взвешенным средним n соседних членов, где n - ширина окна. Вместо среднего можно использовать медиану значений, попавших в окно.

Для построения модели тренда проведем сглаживание временного ряда с периодом 4. Для этого воспользуемся формулой средней хронологической. В нашем случае данная формула будет выглядеть следующим образом.

Где t = 3,…,n-2; к = 4 - длина цикла сезонности.

Обоснование выбора типа модели основного тренда по свойствам кривых роста

Построив 4 прироста по сглаженному ряду наблюдаем следующую картину

Только четвертые приросты колеблются вокруг постоянного уровня, скорее всего мы можем взять полином третьей степени.

С помощью функции ЛИНЕЙН построим уравнение модели:

ЛИНЕЙН

-0,533959337

19,1

-147,7185881

2244,95

0,114276663

3,6

34,83892337

98,01

0,987638656

34,4

#Н/Д

#Н/Д

319,5894121

12,0

#Н/Д

#Н/Д

1135607,793

14213,3

#Н/Д

#Н/Д

Уравнение модели:

y = -,53x3+19,1x2-147,72x+2244,95

Выбор типа модели: аддитивная, мультипликативная или смешанная

Общая модель. Основная идея сезонной декомпозиции проста. В общем случае временной ряд типа того, который описан выше, можно представить себе состоящим из трех различных компонент: (1) сезонной компоненты (обозначается St, где t обозначает момент времени), (2) тренда (Tt), (3) циклической компоненты (Ct). Разница между циклической и сезонной компонентой состоит в том, что последняя имеет регулярную (сезонную) периодичность, тогда как циклические факторы обычно имеют более длительный эффект, который к тому же меняется от цикла к циклу. Конкретные функциональные взаимосвязи между этими компонентами могут иметь самый разный вид. Однако, можно выделить два основных способа, с помощью которых они могут взаимодействовать: аддитивно и мультипликативно.

Аддитивная модель:

Xt = TCt + St + It

Мультипликативная модель:

Xt = TCt*St*It

Здесь Xt обозначает значение временного ряда в момент времени t. Если имеются какие-то априорные сведения о циклических факторах, влияющих на ряд (например, циклы деловой конъюнктуры), то можно использовать оценки для различных компонент для составления прогноза будущих значений ряда.

тренд сезонность временный ряд

Так как график не имеет ярко выраженной склонности к возрастанию, а колеблется примерно в одном и том же диапазоне, я взял аддитивную модель.

Построим аддитивную модель.

Для того чтобы определить порядок сезонности, выполним ряд вычислений:

Квартал

Сумма разностей

Число повторений квартала

Средняя разность

Вычетание среднего

I

-3076,17

5

-615,23

-596,00

II

-2462,47

5

-492,49

-473,26

III

-1655,05

5

-331,01

-311,77

IV

6808,96

5

1361,79

1381,03

 

 

Среднее

-19,24

0,00

Сезонная компонента S Квартал

S

I

-596,00

II

-473,26

III

-311,77

IV

1381,03

Построим график Y и T+S чтобы посмотреть на качество аддитивной модели

Выбор типа остатков и корректировка модели

Для каждого значения t определим ошибку модели (остаточную компоненту временного ряда) по следующей формуле:

Для того, чтобы выбрать тип остатков: e = y - y или e = (y/y) - 1, построим оба графика

Так как на первом графике димперсия более постоянна, выберем e = y - y.

С помощью процедуры «Поиск решения» определим новые параметры модели

Новые параметры модели:

 

-0,414011172

13,53396

-76,8736

2000,126

y = -0,41*x3+13,53*x2-76,87*x+2000,13

Проверка гипотезы о нормальном законе распределения остатков по ассиметрии и эксцессу:

0,47

0,22

0,59

-0,76

0,76

H0:

остатки распределены в соответствии с нормальным законом

0,22

<

0,71

0,55

<

1,14

H1:

остатки не распределены в соответствии с нормальным законом

Следовательно, остатки нормально распределены

Проверка постоянства дисперсии остатков:

n1 =

10

D1 =

10892,72291

???

0,05

n2 =

10

D2 =

17834,30861

 

 

F =

0,61077349

F1кр =

4,025994158

Fкр =

3,178893105

F2кр =

0,314574906

 

 

H0:

?????????

H0:

?????????

H1:

??????????

H1:

?????????

Следовательно, остатки гомокедастичны.

Проверка среднего значения остатков:

t

e

e - e1ср

e - e2ср

1

-0,079

-0,086

 

n1 =

10

2

0,036

0,028

 

n2 =

10

3

0,001

-0,006

 

4

0,007

-0,001

 

e1ср =

0,007623145

5

0,112

0,105

 

e2ср =

0,006

6

-0,016

-0,023

 

7

-0,025

-0,032

-0,031

???

0,05

8

-0,038

-0,045

-0,044

9

0,122

0,115

0,116

S =

0,001102894

10

-0,118

-0,125

-0,124

11

-0,037

-0,045

-0,043

H0:

???????

12

-0,030

-0,038

-0,036

H1:

????????

13

-0,012

-0,020

-0,018

14

0,111

0,104

0,105

t =

2,33343192

15

0,000

-0,008

-0,007

16

0,049

0,042

0,043

t1кр =

-2,100922037

17

-0,039

-0,047

-0,046

t2кр =

2,100922037

18

0,016

0,008

0,010

19

0,073

0,065

0,066

20

0,017

0,009

0,011

Суммкв:

0,074666557

0,054814471

Среднее значение остатков непостоянно

Проверка гипотезы об отсутствии автокорелляции:

?

r

n - ?

t

t1кр

t2кр

 

1

-0,27

19

-1,14

-2,11

2,11

+

2

-0,02

18

-0,08

-2,12

2,12

+

3

-0,11

17

-0,43

-2,13

2,13

+

4

-0,09

16

-0,34

-2,14

2,14

+

5

0,17

15

0,63

-2,16

2,16

+

6

-0,19

14

-0,69

-2,18

2,18

+

7

0,17

13

0,59

-2,20

2,20

+

8

-0,54

12

-2,01

-2,23

2,23

+

Автокорелляция присутствует.

Применение модели для прогноза

Прогноз на основании трендовых моделей (кривых роста) содержит два элемента: точечный и интервальный прогнозы. Точечный прогноз -- это прогноз, которым называется единственное значение прогнозируемого показателя. Это значение определяется подстановкой в уравнение выбранной кривой роста величины времени t, соответствующей периоду упреждения: t=n+1; t=n+2 и т. д. Очевидно, что точное совпадение фактических данных в будущем и прогностических точечных оценок маловероятно. Поэтому точечный прогноз должен сопровождаться двусторонними границами, т.е. интервальным прогнозом.

Интервальный прогноз на базе трендовых моделей осуществляется путем расчета доверительного интервала -- такого интервала, в котором с определенной вероятностью можно ожидать появления фактического значения прогнозируемого экономического показателя. Методы, разработанные для статистических совокупностей, позволяют определить доверительный интервал, зависящий от стандартной ошибки оценки прогнозируемого показателя, от времени упреждения прогноза, от количества уровней во временном ряду и от уровня значимости (ошибки) прогноза.

Точечный прогноз на два периода сезонности:

 

t

t2

t3

T

S

y

Y

2004

1

1

1

1936,372498

-596,00

1340,38

1234,7

 

2

4

8

1897,202698

-473,26

1423,94

1474,6

 

3

9

27

1880,132678

-311,77

1568,36

1570,7

 

4

16

64

1882,678372

1381,03

3263,71

3286,3

2005

5

25

125

1902,355711

-596,00

1306,36

1453,3

 

6

36

216

1936,68063

-473,26

1463,42

1440,3

 

7

49

343

1983,16906

-311,77

1671,40

1630,4

 

8

64

512

2039,336936

1381,03

3420,37

3292,0

2006

9

81

729

2102,700189

-596,00

1506,70

1690,9

 

10

100

1000

2170,774753

-473,26

1697,52

1497,8

 

11

121

1331

2241,076561

-311,77

1929,30

1857,9

 

12

144

1728

2311,121546

1381,03

3692,15

3581,8

2007

13

169

2197

2378,425641

-596,00

1782,43

1761,1

 

14

196

2744

2440,504778

-473,26

1967,25

2186,5

 

15

225

3375

2494,874891

-311,77

2183,10

2182,7

 

16

256

4096

2539,051913

1381,03

3920,08

4113,4

2008

17

289

4913

2570,551777

-596,00

1974,56

1897,2

 

18

324

5832

2586,890415

-473,26

2113,63

2147,4

 

19

361

6859

2585,583761

-311,77

2273,81

2438,8

 

20

400

8000

2564,147747

1381,03

3945,18

4012,2

2009

21

441

9261

2520,098307

-596,00

1924,10

 

 

22

484

10648

2450,951373

-473,26

1977,69

 

 

23

529

12167

2354,222879

-311,77

2042,45

 

 

24

576

13824

2227,428757

1381,03

3608,46

 

2010

25

625

15625

2068,08494

-596,00

1472,09

 

 

26

676

17576

1873,707362

-473,26

1400,45

 

 

27

729

19683

1641,811955

-311,77

1330,04

 

 

28

784

21952

1369,914653

1381,03

2750,94

 

Графическое представление:

Интервальный прогноз на два периода сезонности:

L =

8

???

0,05

Sy =

0,064905

n =

20

t? =

2,100922037

 

 

 

t

t2

t3

t4

y

Uy

2004

1

1

1

1

1340,376

1340,376049

1340,376

 

2

4

8

16

1423,945

1423,944764

1423,945

 

3

9

27

81

1568,359

1568,358863

1568,359

 

4

16

64

256

3263,707

3263,706569

3263,707

2005

5

25

125

625

1306,359

1306,359263

1306,359

 

6

36

216

1296

1463,423

1463,422696

1463,423

 

7

49

343

2401

1671,395

1671,395245

1671,395

 

8

64

512

4096

3420,365

3420,365133

3420,365

2006

9

81

729

6561

1506,704

1506,70374

1506,704

 

10

100

1000

10000

1697,517

1697,516819

1697,517

 

11

121

1331

14641

1929,303

1929,302746

1929,303

 

12

144

1728

20736

3692,15

3692,149743

3692,15

2007

13

169

2197

28561

1782,429

1782,429192

1782,429

 

14

196

2744

38416

1967,247

1967,246844

1967,247

 

15

225

3375

50625

2183,101

2183,101076

2183,101

 

16

256

4096

65536

3920,08

3920,080111

3920,08

2008

17

289

4913

83521

1974,555

1974,555328

1974,555

 

18

324

5832

104976

2113,632

2113,632481

2113,632

 

19

361

6859

130321

2273,81

2273,809946

2273,81

 

20

400

8000

160000

3945,176

3945,175944

3945,176

2009

21

441

9261

194481

1924,102

1924,27099

1923,933

 

22

484

10648

234256

1977,693

1977,868154

1977,519

 

23

529

12167

279841

2042,449

2042,630093

2042,268

 

24

576

13824

331776

3608,457

3608,645039

3608,269

2010

25

625

15625

390625

1472,088

1472,284381

1471,893

 

26

676

17576

456976

1400,449

1400,653871

1400,245

 

27

729

19683

531441

1330,038

1330,251881

1329,824

 

28

784

21952

614656

2750,943

2751,166625

2750,719

?

210

2870

44100

722666

 

 

 

 

406

7714

164836

3756718

 

 

 

Графическое представление

Размещено на Allbest.ru


Подобные документы

  • "Наивная" модель прогнозирования. Прогнозирование методом среднего и скользящего среднего. Метод опорных векторов, деревьев решений, ассоциативных правил, системы рассуждений на основе аналогичных случаев, декомпозиции временного ряда и кластеризации.

    курсовая работа [2,6 M], добавлен 02.12.2014

  • Исходный текст программы и ее экранная форма. Программа вычисления и выдачи на печать суммы/произведения элементов бесконечного числового ряда, вычисления числового ряда для известного числа членов ряда. Значение максимального элемента в матрице.

    контрольная работа [29,0 K], добавлен 07.12.2010

  • Последовательность разработки чертежа и модели с типоразмерами из параметрического ряда. Построение таблицы переменных в соответствии с исходными данными. Проектирование параметрической модели в системе Компас-3D, внешние переменные для чертежа детали.

    практическая работа [5,9 M], добавлен 14.04.2016

  • Проверка наличия линейной связи между соответствующими показателями деятельности коммерческих банков Украины в модуле Multiple Regression ППП Statistica. Расчет теоретических значений зависимой переменной и ошибки модели, вид графика линейной функции.

    лабораторная работа [1,5 M], добавлен 19.05.2011

  • Особенности и основные этапы разработка компьютерной модели расчета траектории неуправляемого летательного аппарата, анализ модели динамики. Метод Эйлера как линейное приближение, использующее первые два члена ряда Тейлора, способы решения задач.

    курсовая работа [930,6 K], добавлен 14.01.2013

  • Анализ матрицы коэффициентов парной корреляции. Выбор факторных признаков для построения двухфакторной регрессионной модели. Оценка параметров регрессии по методу наименьших квадратов. Нахождение определителей матриц. Применение инструмента Регрессия.

    контрольная работа [1,0 M], добавлен 13.01.2013

  • Приведена оптимизация расходов и трудозатрат персонала. Реализация модели ARIMA (модели Бокса-Дженкинса), являющейся интегрированной композицией метода авторегрессии и модели скользящего среднего. Применение средств программного продукта Matlab 2013a.

    дипломная работа [876,7 K], добавлен 19.09.2019

  • Построение концептуальной модели и метод имитационного моделирования. Определение переменных уравнений математической модели и построение моделирующего алгоритма. Описание возможных улучшений системы и окончательный вариант модели с результатами.

    курсовая работа [79,2 K], добавлен 25.06.2011

  • Применение комплексного математического моделирования в проектировании. Обзор численных методов в моделировании. Решение дифференциальных уравнений в MathCAD. Анализ исходных и результирующих данных. Описание реализации базовой модели в MathCAD.

    курсовая работа [240,5 K], добавлен 18.12.2011

  • Общая характеристика ателье "Вита", схема модели рабочего процесса. Исследование заданной системы с помощью моделирования динамических рядов, модели типа "система массового облуживания". Построение имитационной модели деятельности данного ателье.

    курсовая работа [1,4 M], добавлен 01.06.2016

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.