Автор работы: Пользователь скрыл имя, 12 Ноября 2009 в 15:50, Не определен
Основы прогнозирования
На основе данных итоговой строки определим параметры линеаризированной модели показательного тренда lna0 и lna1:
(34)
(35)
Найдем значение a0, пропотенцировав lna0. Тогда . Найдем a1, пропотенцировав величину lna1. Тогда . В результате расчетов параметров показательный тренд будет иметь следующий конкретный вид:
В
одной системе координат, построим
3 графика: исходный временной ряд, линейную
и показательную модели. Для построения
графиков используем таблицу 15. Рассчитаем
значения линейного и показательного
трендов для каждого момента времени,
подставив соответствующие значения t
в уравнения
и
Таблица 15
t | yt | ||
1 | 2 | 3 | 4 |
1 | 10,2 | 10,326 | 10,407 |
2 | 10,8 | 10,723 | 10,751 |
3 | 10,4 | 11,120 | 11,106 |
4 | 11,9 | 11,517 | 11,472 |
5 | 12,2 | 11,914 | 11,851 |
6 | 12,5 | 12,311 | 12,242 |
7 | 13,1 | 12,708 | 12,646 |
8 | 12,4 | 13,105 | 13,063 |
9 | 13,6 | 13,502 | 13,494 |
10 | 14,3 | 13,899 | 13,939 |
11 | 14,9 | 14,296 | 14,399 |
12 | 13,8 | 14,693 | 14,875 |
Нарисуем
график по исходным данным на рис.6:
На основе визуального анализа можно сделать вывод: что совпадение выбранного нелинейного тренда с реальным трендом временного ряда очевидно.
Выбранные тренды имеют одинаковое число параметров, поэтому для расчета критерия наименьшей суммы квадратов отклонений выберем формулу
Для проведения промежуточных расчетов используем таблицу 16:
Таблица 16
t | yt | ||||
1 | 2 | 3 | 4 | 5 | 6 |
1 | 10,2 | 10,326 | 10,407 | 0,016 | 0,043 |
2 | 10,8 | 10,723 | 10,751 | 0,006 | 0,002 |
3 | 10,4 | 11,120 | 11,106 | 0,518 | 0,498 |
4 | 11,9 | 11,517 | 11,472 | 0,147 | 0,183 |
5 | 12,2 | 11,914 | 11,851 | 0,082 | 0,122 |
6 | 12,5 | 12,311 | 12,242 | 0,036 | 0,067 |
7 | 13,1 | 12,708 | 12,646 | 0,154 | 0,206 |
8 | 12,4 | 13,105 | 13,063 | 0,497 | 0,440 |
9 | 13,6 | 13,502 | 13,494 | 0,010 | 0,011 |
10 | 14,3 | 13,899 | 13,939 | 0,161 | 0,130 |
11 | 14,9 | 14,296 | 14,399 | 0,365 | 0,251 |
12 | 13,8 | 14,693 | 14,875 | 0,797 | 1,156 |
- | - | - | - | 2,789 | 3,109 |
В графу 5 внесем отклонение фактических данных от линейного тренда: для t=1 =(10,2–10,326)2=0,016, для t=2 =(10,8–10,723)2=0,006 и т.д. Аналогичным образом в графу 6 впишем отклонения фактических данных от показательного тренда, и подведем итог.
Сравним значения критерия наименьшей суммы квадратов отклонений для линейного и показательного трендов. Для линейного тренда критерий равен 2,789 (итог графы 5), а для показательного тренда 3,109 (итог графы 6). Коэффициент для линейного тренда меньше, чем для показательного, поэтому линейный тренд лучше аппроксимирует исходные данные. Следовательно, для прогнозирования необходимо взять линейный тренд.
Чтобы оценить адекватность выбранной линейной трендовой модели теоретическому тренду временного ряда, найдем разность еt между исходными данными уt и нашей трендовой моделью по формуле
результаты занесем в таблицу 17.
Таблица 17
t | yt | ||
1 | 2 | 3 | 4 |
1 | 10,2 | 10,326 | -0,126 |
2 | 10,8 | 10,723 | 0,077 |
3 | 10,4 | 11,120 | -0,72 |
4 | 11,9 | 11,517 | 0,383 |
5 | 12,2 | 11,914 | 0,286 |
6 | 12,5 | 12,311 | 0,189 |
7 | 13,1 | 12,708 | 0,392 |
8 | 12,4 | 13,105 | -0,705 |
9 | 13,6 | 13,502 | 0,098 |
10 | 14,3 | 13,899 | 0,401 |
11 | 14,9 | 14,296 | 0,604 |
12 | 13,8 | 14,693 | -0,893 |
- | - | - | - |
Построим график ряда отклонения еt по данным графы 4 на рис.7:
Визуальный анализ показывает, что колебание величины еt не содержит элементов тенденции, т.е. носит случайный характер.
Оценим
адекватность выбранной модели тренда
исходному ряду на основе анализа
данных ряда отклонений еt. Величина еt должна отвечать следующим
четырем условиям (требованиям):
Условие
1. Колебание величины
еt
должно носить случайный
характер.
Проверим условие 1 с помощью критерия поворотных точек.
Величина еt считается поворотной, если она соответствует одному из двух условий:
еt-1< еt >еt+1 или еt-1> еt <еt+1
Для этого на базе данных графы 4 определим поворотные точки и в графе 5 проставим соответствующие им значения. Рассматриваемые точки обозначим Рt. Тем точкам, которые будут поворотными, присвоим значение Рt=1. А тем точкам, которые не будут поворотными, присвоим значение Рt=0. результаты занесем в таблицу 18.
Таблица 18
t | yt | Pt | ||
1 | 2 | 3 | 4 | 5 |
1 | 10,2 | 10,326 | -0,126 | - |
2 | 10,8 | 10,723 | 0,077 | 1 |
3 | 10,4 | 11,120 | -0,72 | 1 |
4 | 11,9 | 11,517 | 0,383 | 1 |
5 | 12,2 | 11,914 | 0,286 | 0 |
6 | 12,5 | 12,311 | 0,189 | 1 |
7 | 13,1 | 12,708 | 0,392 | 1 |
8 | 12,4 | 13,105 | -0,705 | 1 |
9 | 13,6 | 13,502 | 0,098 | 0 |
10 | 14,3 | 13,899 | 0,401 | 0 |
11 | 14,9 | 14,296 | 0,604 | 1 |
12 | 13,8 | 14,693 | -0,893 | - |
- | - | - | - | 7 |
Определим общее число поворотных точек по формуле
Результат занесем в итоговую строку графы 6.
Для проверки выполнения условия 1 выдвинем нулевую гипотезу Н0: колебание величины еt носит случайный характер.
Для проверки этой гипотезы определим математическое ожидание числа поворотных точек по формуле:
и его дисперсию по формуле:
Кроме того, для проверки нулевой гипотезы используем вероятность, равную 95%, которой соответствует коэффициент доверия t=1,96. С помощью формулы
проверим нулевую гипотезу, подставив в нее значения М(Р), D(P), t :
Расчет
показывает, общее
число поворотных
точек – 7 находится
в требуемом интервале.
Это позволяет сделать
следующий вывод: с вероятностью 0,95 (95%)
колебание величины
еt
носит случайный характер
и, следовательно,
отвечает данному условию.
Условие 2.
Распределение величины
еt
соответствует нормальному
распределению.
Проверим распределение еt на соответствие нормальному распределению. Вначале определим среднее квадратическое отклонение по формуле:
для определения расчетного значения критерия RSр из графы 4 таблицы 18 найдем максимальное emax=0,604 и минимальное emin=−0,893 значения. Расчетное значение критерия RSр найдем по формуле:
Следующим шагом проверки условия 2 является нахождение табличного значения RS-критерия – RST по приложению 3 (2, стр. 71). В таблице приводятся нижнее и верхнее значения RS-критерия для n=10 и n=20; а у нас n=12. Для нахождения нижнего и верхнего значений RS-критерия для n=12 используем линейную интерполяцию. Найдем величину RS12н:
Информация о работе Прогнозирование емкости и конъюнктуры рынка