Автор работы: Пользователь скрыл имя, 14 Марта 2012 в 13:33, курсовая работа
Целью данной курсовой работы является подробное раскрытие сущности имитационного моделирования его виды, построение имитационных моделей, проведения расчётов с использованием данного метода, и составление прогнозирования.
Введение ......................................................................................................................4
Глава 1. Общие понятия имитационной модели и имитационного моделирования ……....................................................................................................5
1.1. Моделирование входных данных ............................................................6
1.2. Расчет показателей динамики развития экономических процессов....................................................................................................................8
1.3. Методы корректировки динамического ряда .......................................10
1.4. Определение трендовой компоненты ....................................................11
1.5.Определение сезонной компоненты .......................................................13
1.5.1. Расчёт сезонной волны выработки электроэнергии ...............13
1.5.2. Метод Четвертикова....................................................................14
1.5.3. Нахождение сезонной составляющей с помощью ряда Фурье в качестве аналитической модели сезонности ..........................................................16
1.5.4. Оценка адекватности и точности трендовых моделей.............17
1.5.5. Прогнозирование динамики на основе трендовых моделей.......................................................................................................................19
Глава 2. Построение имитационной модели и составление прогноза................23
2.1. Абсолютный и средний абсолютный прирост.........................24
2.2. Показатели роста и прироста...................................................25
2.3. Выявление аномальных уровней ряда...........................................25
2.4. Определение наличие тренда во временном ряду........................26
2.5. Автокорреляционная функция и временной лаг....................27
2.6. Построение графика сезонной волны........................................28
2.7. Построение аналитической модели ряда с использованием метода Фурье........................................................................................................................30
2.8. Оценка адекватности и точности трендовой модели..................................................................................................................31
2.9. Построение интервального и точного прогноза на 4 уровня вперед...................................................................................................................34
Заключение ...............................................................................................................36
Список используемой литературы .........................................................................37
Проверка равенства математического ожидания случайной компоненты нулю, если она распределена по нормальному закону, осуществляется на основе t-критерия Стьюдента. Расчетное значение этого критерия задается формулой
,
где - среднее арифметическое значение уровней остаточной последовательности , - стандартное (среднеквадратическое отклонение) для этой последовательности.
Если расчетное значение меньше табличного значения статистики Стьюдента с заданным уровнем значимости и числом степеней свободы , то гипотеза о равенстве нулю математического ожидания случайной последовательности принимается; в противном случае эта гипотеза отвергается и модель считается неадекватной.
Проверка независимости значений уровней случайной компоненты, т.е. проверка отсутствия существенной автокорреляции в остаточной последовательности может осуществляться по ряду критериев, наиболее распространенным из которых является d-критерий Дарбина-Уотсона Расчетное значение этого критерия определяется по формуле
. (1.35)
Проверка на точность модели производится по формуле
среднеквадратическое отклонение
средняя относительная ошибка аппроксимации
Прогноз на основании трендовых моделей (кривых роста) содержит два элемента: точечный и интервальный прогнозы. Точечный прогноз – это прогноз, которым называется единственное значение прогнозируемого показателя. Это значение определяется подстановкой в уравнение выбранной кривой роста величины времени , соответствующей периоду упреждения: , и т. д. Такой прогноз называется точечным, так как на графике его можно изобразить в виде точки.
Очевидно, что точное совпадение фактических данных в будущем и прогностических точечных оценок маловероятно. Поэтому точечный прогноз должен сопровождаться двусторонними границами, то есть указанием интервала значений, в котором с достаточной долей уверенности можно ожидать появления прогнозируемой величины. Установление такого интервала называется интервальным прогнозом.
Стандартная (средняя квадратическая) ошибка оценки прогнозируемого показателя определяется по формуле:
,
где - фактическое значение уровня временного ряда для времени t; - расчетная оценка соответствующего показателя по модели (например, по уравнению кривой роста); n – количество уровней в исходном временном ряду; k – число параметров модели.
В случае прямолинейного тренда для расчета доверительного интервала можно использовать аналогичную формулу для парной регрессии, таким образом, доверительный интервал прогноза в этом случае будет иметь вид
, (1.39)
где L – период упреждения; - точечный прогноз по модели на -й момент времени; n – количество наблюдений во временном ряду; стандартная ошибка оценки прогнозируемого показателя, рассчитанная по формуле (1) для числа параметров модели, равного двум; - табличное значение критерия Стьюдента для уровня значимости и числа степеней свободы, равного .
Если выражение
обозначить через K, то формула для доверительного интервала примет вид:
. (1.41)
Значения величины К для оценки доверительных интервалов прогноза относительно линейного тренда табулированы.
Иногда для расчета доверительных интервалов прогноза относительно линейного тренда применяют приведенную выше формулу в несколько преобразованном виде:
(1.42)
Здесь t – порядковый номер уровня ряда (t = 1, 2, . . . , n); - время, для которого делается прогноз; - время, соответствующее середине периода наблюдений для исходного ряда, например =(n+1):2; суммирование ведется по всем наблюдениям.
Эту формулу можно упростить, если, как часто делается на практике, перенести начало отсчета времени на середину периода наблюдений (=0):
. (1.43)
Формула для расчета доверительных интервалов прогноза относительно тренда, имеющего вид полинома второго или третьего порядка, выглядит следующим образом:
, (1.44)
Аналогично вычисляются доверительные интервалы для экспоненциальной кривой роста, а также для кривых роста, имеющих асимптоту (модифицированная экспонента, кривая Гомперца, логистическая кривая), если значение асимптоты известно.
Глава 2. Построение имитационной модели и составление прогноза
По данным таблицы 2 найти:
1. Цепной абсолютный прирост
2. Базисный абсолютный прирост
3. Цепной средний абсолютный прирост
4. Базисный средний абсолютный прирост
5. Прирост за весь период наблюдения
6. Цепной коэффициент роста
7. Базисный коэффициент роста
8. Цепной коэффициент прироста
9. Базисный коэффициент прироста
10. Цепной темп роста
11. Базисный темп роста
12. Цепной темп прироста
13. Базисный темп прироста
14. Выявить аномальные уровни ряда
15. Определить наличие тренда во временном ряду
16. Найти первые 12 коэффициентов автокорреляции, построить автокорреляционную функцию и найти временной лаг
17. Построить график сезонной волны
18. Построить аналитическую модель ряда с использованием метода Фурье
19. Оценить адекватность и точность трендовой модели
20. Построить интервальный и точный прогноз на 4 уровня вперед
Таблица 2.
Первичные данные.
| Yt | t | Yt | t | Yt | t | Yt |
1,00 | 9,00 | 26,00 | 9,00 | 51,00 | 8,00 | 76,00 | 8,00 |
2,00 | 7,00 | 27,00 | 6,00 | 52,00 | 6,00 | 77,00 | 9,00 |
3,00 | 4,00 | 28,00 | 3,00 | 53,00 | 4,00 | 78,00 | 10,00 |
4,00 | 1,00 | 29,00 | 2,00 | 54,00 | 4,00 | 79,00 | 13,00 |
5,00 | -1,00 | 30,00 | 2,00 | 55,00 | 6,00 | 80,00 | 15,00 |
6,00 | -1,00 | 31,00 | 3,00 | 56,00 | 9,00 | 81,00 | 17,00 |
7,00 | 1,00 | 32,00 | 6,00 | 57,00 | 12,00 | 82,00 | 17,00 |
t | Yt | t | Yt | t | Yt | t | Yt |
8,00 | 4,00 | 33,00 | 10,00 | 58,00 | 15,00 | 83,00 | 17,00 |
9,00 | 7,00 | 34,00 | 13,00 | 59,00 | 17,00 | 84,00 | 15,00 |
10,00 | 10,00 | 35,00 | 14,00 | 60,00 | 17,00 | 85,00 | 13,00 |
11,00 | 12,00 | 36,00 | 14,00 | 61,00 | 15,00 | 86,00 | 11,00 |
12,00 | 12,00 | 37,00 | 13,00 | 62,00 | 13,00 | 87,00 | 10,00 |
13,00 | 11,00 | 38,00 | 10,00 | 63,00 | 10,00 | 88,00 | 9,00 |
14,00 | 8,00 | 39,00 | 7,00 | 64,00 | 7,00 | 89,00 | 10,00 |
15,00 | 5,00 | 40,00 | 4,00 | 65,00 | 5,00 | 90,00 | 12,00 |
16,00 | 2,00 | 41,00 | 3,00 | 66,00 | 5,00 | 91,00 | 14,00 |
17,00 | 0,00 | 42,00 | 3,00 | 67,00 | 7,00 | 92,00 | 16,00 |
18,00 | 1,00 | 43,00 | 5,00 | 68,00 | 10,00 | 93,00 | 18,00 |
19,00 | 2,00 | 44,00 | 7,00 | 69,00 | 13,00 | 94,00 | 19,00 |
20,00 | 5,00 | 45,00 | 11,00 | 70,00 | 16,00 | 95,00 | 18,00 |
21,00 | 8,00 | 46,00 | 14,00 | 71,00 | 18,00 | 96,00 | 17,00 |
22,00 | 11,00 | 47,00 | 16,00 | 72,00 | 18,00 | 97,00 | 15,00 |
23,00 | 13,00 | 48,00 | 16,00 | 73,00 | 17,00 | 98,00 | 12,00 |
24,00 | 13,00 | 49,00 | 14,00 | 74,00 | 14,00 | 99,00 | 11,00 |
25,00 | 12,00 | 50,00 | 12,00 | 75,00 | 11,00 | 100,00 | 10,00 |
2.1 Абсолютный и средний абсолютный прирост
Найдём абсолютный прирост и средний абсолютный прирост(1.2-1.4),
Таблица 3.
Абсолютный и средний абсолютный прирост.
Данные | Абсолютный прирост | Средний абсолютный прирост | |||
Цепной | Базисный | Цепной | Базисный | ||
t | Yt | ∆ Yt | ∆ Yt | ∆ Yt | ∆ Yt |
| 9 | - | - | - | - |
2 | 7 | -2 | -2 | -1 | -2 |
3 | 4 | -3 | -5 | -1,5 | -2,5 |
4 | 1 | -3 | -8 | -1,5 | -2,67 |
5 | -1 | -2 | -10 | -1 | -2,5 |
95 | 18 | -1 | 9 | -0,5 | 0,1 |
... | ... | ... | ... | ... | ... |
96 | 17 | -1 | 8 | -0,5 | 0,08 |
97 | 15 | -2 | 6 | -1 | 0,06 |
98 | 12 | -3 | 3 | -1,5 | 0,03 |
99 | 11 | -1 | 2 | -0,5 | 0,02 |
100 | 10 | -1 | 1 | -0,5 | 0,01 |
Прирост за весь период наблюдения для данного временного ряда равен
2.2 Показатели роста и прироста
Для вычисление коэффициент роста используем формулы (1.8), а для коэффициента прироста (1.5-1.7) и темпа прироста (1.9-1.10) . Также найдём темп обоих показателей по формулам
Таблица 4.
Показатели роста и прироста.
Данные | Коэффициент роста | Коэффициент прироста | Темп роста | Темп прироста | |||||
Цепной | Базисный | Цепной | Базисный | Цепной | Базисный | Цепной | Базисный | ||
t | Yt | Ki(р) | Ki(р) | Ki(пр) | Ki(пр) | Тi(р) | Тi(р) | Тi(пр) | Тi(пр) |
1 | 9 | - | - | - | - | - | - | - | - |
2 | 7 | 0,78 | 0,78 | -0,22 | -0,22 | 77,78% | 77,78% | -22,22% | -22,22% |
3 | 4 | 0,57 | 0,44 | -0,43 | -0,56 | 57,14% | 44,44% | -42,86% | -55,56% |
4 | 1 | 0,25 | 0,11 | -0,75 | -0,89 | 25,00% | 11,11% | -75,00% | -88,89% |
5 | -1 | -1 | -0,11 | -2 | -1,11 | -100,00% | -11,11% | -200,00% | -111,11% |
... | ... | ... | ... | ... | ... | ... | ... | ... | ... |
95 | 18 | 0,95 | 2 | -0,05 | 1 | 94,74% | 200,00% | -5,26% | 100,00% |
96 | 17 | 0,94 | 1,89 | -0,06 | 0,89 | 94,44% | 188,89% | -5,56% | 88,89% |
97 | 15 | 0,88 | 1,67 | -0,12 | 0,67 | 88,24% | 166,67% | -11,76% | 66,67% |
98 | 12 | 0,8 | 1,33 | -0,2 | 0,33 | 80,00% | 133,33% | -20,00% | 33,33% |
99 | 11 | 0,92 | 1,22 | -0,08 | 0,22 | 91,67% | 122,22% | -8,33% | 22,22% |
100 | 10 | 0,91 | 1,11 | -0,09 | 0,11 | 90,91% | 111,11% | -9,09% | 11,11% |