Автор: Пользователь скрыл имя, 04 Апреля 2012 в 14:24, контрольная работа
1) Построить адаптивную мультипликативную модель Хольта-Уинтерса с учетом сезонного фактора, приняв параметры сглаживания α1 = 0,3, α2 = 0,6, α3 = 0,3.
2) Оценить точность построенной модели с использованием средней относительной ошибки аппроксимации.
3) Оценить адекватность построенной модели на основе исследования:
Случайности остаточной компоненты по критерию пиков;
Независимости уровней ряда остатков по d- критерию (критические значения d1 = 1,10 и d2 = 1,37) и по первому коэффициенту автокорреляции при критическом значении r1 = 0,32;
Нормальности распределения остаточной компоненты по R / S критерию с критическими значениями от 3 до 4,21.
ЗАДАНИЕ 1
Приведены данные о кредитах от коммерческого банка на жилищное строительство за 4 года:
t | Y(t) |
1 | 36 |
2 | 46 |
3 | 55 |
4 | 35 |
5 | 39 |
6 | 50 |
7 | 61 |
8 | 37 |
9 | 42 |
10 | 54 |
11 | 64 |
12 | 40 |
13 | 47 |
14 | 58 |
15 | 70 |
16 | 43 |
Требуется:
1) Построить адаптивную мультипликативную модель Хольта-Уинтерса с учетом сезонного фактора, приняв параметры сглаживания α1 = 0,3, α2 = 0,6, α3 = 0,3.
2) Оценить точность построенной модели с использованием средней относительной ошибки аппроксимации.
3) Оценить адекватность построенной модели на основе исследования:
Случайности остаточной компоненты по критерию пиков;
Независимости уровней ряда остатков по d- критерию (критические значения d1 = 1,10 и d2 = 1,37) и по первому коэффициенту автокорреляции при критическом значении r1 = 0,32;
Нормальности распределения остаточной компоненты по R / S критерию с критическими значениями от 3 до 4,21.
4) Построить точечный прогноз на 4 шага вперед, то есть на 1 год.
5) Отразить на графике фактические, расчетные и прогнозные данные.
Решение:
1) Строим адаптивную мультипликативную модель Хольта-Уинтерса с учетом сезонного фактора
Общий вид модели:
где – расчетное значение уровня для момента времени t с периодом упреждения k;
k – период упреждения;
L – период сезонности;
(t – L) - индекс сезонного коэффициента за аналогичный период прошлого года;
Ft – мультипликативный индекс сезонности;
a0(t); a1(t) – параметры модели;
Находим начальные оценки параметров a0(0) и а1(0), построив линейную трендовую модель при n = 8. Для нахождения параметров строим систему нормальных уравнений:
Получили модель:
t | Y(t) | t *Y(t) | t2 | Yрасч(t) |
1 | 36 | 36 | 1 | 41,917 |
2 | 46 | 92 | 4 | 42,762 |
3 | 55 | 165 | 9 | 43,607 |
4 | 35 | 140 | 16 | 44,452 |
5 | 39 | 195 | 25 | 45,297 |
6 | 50 | 300 | 36 | 46,142 |
7 | 61 | 427 | 49 | 46,987 |
8 | 37 | 296 | 64 | 47,832 |
36 | 359 | 1651 | 204 | - |
Находим начальные мультипликативные индексы сезонности:
Производим корректировку параметров по формулам:
, , - параметры адаптации экспоненциального сглаживания.
Расчет скорректированных параметров, мультипликативных индексов сезонности, расчетных значений представлен в таблице:
Год | Квартал | t | Y(t) | a0 | a1 | Yрасч | F |
0 | 1 | - | - | - | - | - | 0,860 |
2 | - | - | - | - | - | 1,080 | |
3 | - | - | - | - | - | 1,280 | |
4 | - | - | 41,071 | 0,845 | - | 0,780 | |
1 | 1 | 1 | 36 | 41,901 | 0,841 | 36,045 | 0,859 |
2 | 2 | 46 | 42,701 | 0,828 | 46,146 | 1,078 | |
3 | 3 | 55 | 43,364 | 0,779 | 55,706 | 1,273 | |
4 | 4 | 35 | 44,354 | 0,842 | 34,451 | 0,786 | |
2 | 1 | 5 | 39 | 45,250 | 0,858 | 38,844 | 0,861 |
2 | 6 | 50 | 46,188 | 0,882 | 49,715 | 1,081 | |
3 | 7 | 61 | 47,326 | 0,959 | 59,915 | 1,283 | |
4 | 8 | 37 | 47,928 | 0,852 | 37,934 | 0,777 | |
3 | 1 | 9 | 42 | 48,781 | 0,852 | 41,995 | 0,861 |
2 | 10 | 54 | 49,732 | 0,882 | 53,645 | 1,084 | |
3 | 11 | 64 | 50,401 | 0,818 | 64,914 | 1,275 | |
4 | 12 | 40 | 51,288 | 0,839 | 39,820 | 0,779 | |
4 | 1 | 13 | 47 | 52,866 | 1,060 | 44,879 | 0,878 |
2 | 14 | 58 | 53,803 | 1,023 | 58,446 | 1,080 | |
3 | 15 | 70 | 54,850 | 1,031 | 69,898 | 1,276 | |
4 | 16 | 43 | 55,678 | 0,970 | 43,527 | 0,775 |
2) Оценим точность модели
Рассчитываем среднюю относительную ошибку аппроксимации:
S = 1 / n · ∑ | E(t) / y(t)| · 100% = 1/16*0,182*100%=1,14%
Ошибка меньше 7%, модель считается точной.
3) Оценим адекватность построенной модели
- Свойство случайности. используем критерий поворотных точек.
Строим график остатков. Всего 10 поворотных точек.
10 > 6, следовательно, критерий поворотных точек выполняется и остатки имеют случайный характер распределения.
t | Y(t) | Yрасч(t) | E(t) | IE(t)/y(t)I | E(t)2 | (Et-Et-1)2 |
1 | 36 | 36,049 | -0,049 | 0,001 | 0,002 |
|
2 | 46 | 46,159 | -0,159 | 0,003 | 0,025 | 0,012 |
3 | 55 | 55,709 | -0,709 | 0,013 | 0,503 | 0,303 |
4 | 35 | 34,424 | 0,576 | 0,016 | 0,332 | 1,651 |
5 | 39 | 38,825 | 0,175 | 0,004 | 0,031 | 0,161 |
6 | 50 | 49,717 | 0,283 | 0,006 | 0,080 | 0,012 |
7 | 61 | 59,939 | 1,061 | 0,017 | 1,126 | 0,605 |
8 | 37 | 37,912 | -0,912 | 0,025 | 0,832 | 3,893 |
9 | 42 | 42,034 | -0,034 | 0,001 | 0,001 | 0,771 |
10 | 54 | 53,681 | 0,319 | 0,006 | 0,102 | 0,125 |
11 | 64 | 64,9 | -0,9 | 0,014 | 0,810 | 1,486 |
12 | 40 | 39,803 | 0,197 | 0,005 | 0,039 | 1,203 |
13 | 47 | 44,962 | 2,038 | 0,043 | 4,153 | 3,389 |
14 | 58 | 58,603 | -0,603 | 0,010 | 0,364 | 6,975 |
15 | 70 | 70,039 | -0,039 | 0,001 | 0,002 | 0,318 |
16 | 43 | 43,659 | -0,659 | 0,015 | 0,434 | 0,384 |
ИТОГО |
|
| 0,585 | 0,182 | 8,835 | 21,288 |
Свойство независимости. Применим критерий Дарбина-Уотсона
d1 = 1,10 d2 = 1,37 2 < dрасч < 4 Следовательно находим d’
d’ = 4 – dрасч = 4 – 2,409 = 1,591
Т.к. 1,37 < d’ = 1,591 < 2 , свойство выполняется, остатки независимы, автокорреляция в ряду остатков отсутствует.
Свойство подчинения нормальному закону распределения. Используем R/S критерий.
Расчетное значение R/S критерия попадает внутрь критических значений (3 < 3,851 < 4,21), следовательно, свойство выполняется.
Вывод: поскольку рядом остатков выполняются все свойства, то модель считается адекватной. Модель качественная.
4) Построим точечный прогноз на 4 шага вперед, то есть на один год.
Yрасч (t+k) = [55,678 + 0,970 · k] · F(t-L+k) – адаптивная мультипликативная модель Хольта-Уинтерса
Ft : F(4;1) = 0,878
F(4;2) = 1,080
F(4;3) = 1,276
F(4;4) = 0,775
t = 17; k = 1 Yрасч(17) = [a0(16)+ a1(16)*1]*F(4;1)=(55,678 + 0,970*1)*0.878 = 49,726
t = 18; k = 2 Yрасч(18) = [a0(16)+ a1(16)*1]*F(4;2)=(55,678 + 0,970*2)*1,080 = 62,246
t = 19; k = 3 Yрасч(19) = [a0(16)+ a1(16)*1]*F(4;3)=(55,678 + 0,970*3)*1,276 = 49,726
t = 20; k = 4 Yрасч(20) = [a0(16)+ a1(16)*1]*F(4;4)=(55,678 + 0,970*4)*0.775 = 49,726
ЗАДАНИЕ 2
Даны цены (открытия, максимальная, минимальная и закрытия) за 10 дней. Интервал сглаживания принять равным 5 дням.
Дни | Цены | ||
МАКС | МИН | ЗАКР | |
1 | 600 | 550 | 555 |
2 | 560 | 530 | 530 |
3 | 536 | 501 | 524 |
4 | 545 | 521 | 539 |
5 | 583 | 540 | 569 |
6 | 587 | 562 | 581 |
7 | 582 | 561 | 562 |
8 | 573 | 556 | 573 |
9 | 610 | 579 | 592 |
10 | 645 | 585 | 645 |
Информация о работе Контрольная работа по "Финансовой математике"