Автор: Пользователь скрыл имя, 09 Января 2012 в 12:42, контрольная работа
Требуется:
Построить адаптивную мультипликативную модель Хольта-Уинтерса с учетом сезонного фактора, приняв параметры сглаживания α1=0,3; α2=0,6; α3=0,3.
Оценить точность построенной модели с использованием средней относительной ошибки аппроксимации.
В каждом варианте приведены поквартальные данные о кредитах от коммерческого банка на жилищное строительство (в условных единицах) за 4 года (всего 16 кварталов, первая строка соответствует первому кварталу первого года).
Требуется:
Решение:
Мультипликативная модель Хольта-Уинтерса с линейным ростом имеет следующий вид:
Коэффициенты модели рассчитываются по формулам:
Для оценки начальных значений a(0) и b(0) построим линейную модель тренда временного ряда по первым восьми уровням Y(t).
Линейная модель имеет вид: y(t)=a+bt.
Коэффициенты a и b определяем методом наименьших квадратов по формулам:
t | Y(t) | t*Y(t) | y(t) | |
1 | 28 | 28 | 32,58333 | |
2 | 36 | 72 | 33,45238 | |
3 | 43 | 129 | 34,32143 | |
4 | 28 | 112 | 35,19048 | |
5 | 31 | 155 | 36,05952 | |
6 | 40 | 240 | 36,92857 | |
7 | 49 | 343 | 37,79762 | |
8 | 30 | 240 | 38,66667 | |
Σ | 36 | 285 | 1319 | 63 |
ср.знач. | 4,5 | 35,6 | 164,9 | 35,625 |
b(0)= | 0,9 | F(-3)= | 0,859512 | |
a(0)= | 31,7 | F(-2)= | 1,079664 | |
F(-1)= | 1,27462 | |||
Yp(t)=31,7+0,9*t | F(0)= | 0,785766 |
; ; ;
Далее вычисляем по формулам:
t | Y | a(t) | b(t) | F(t) | Yp(t) |
0 | 31,7 | 0,9 | 0,785766 | ||
1 | 28 | 32,58132 | 0,868443 | 0,859438 | 28,00578 |
2 | 36 | 33,41794 | 0,858897 | 1,078225 | 36,11451 |
3 | 43 | 34,11445 | 0,810181 | 1,266126 | 43,68993 |
4 | 28 | 35,13745 | 0,874026 | 0,792429 | 27,44259 |
5 | 31 | 36,02906 | 0,879301 | 0,860025 | 30,94963 |
6 | 40 | 36,96525 | 0,896369 | 1,080548 | 39,79552 |
7 | 49 | 38,11336 | 0,97189 | 1,277833 | 47,93757 |
8 | 30 | 38,71716 | 0,861464 | 0,781882 | 30,97227 |
9 | 34 | 39,56516 | 0,857423 | 0,859615 | 34,03861 |
10 | 44 | 40,51182 | 0,884197 | 1,083881 | 43,67856 |
11 | 52 | 41,18538 | 0,821004 | 1,268684 | 52,89722 |
12 | 33 | 42,06623 | 0,838959 | 0,783439 | 32,84402 |
13 | 39 | 43,64437 | 1,060714 | 0,879998 | 36,88196 |
14 | 48 | 44,57915 | 1,022933 | 1,079594 | 48,45499 |
15 | 58 | 45,63646 | 1,033246 | 1,270022 | 57,85463 |
16 | 36 | 46,45417 | 0,968585 | 0,77835 | 36,56287 |
Для проверки точности модели вычисляем остатки:
Квартал | Объем кредита | отн. погр-ть % | ||
t | Y | Yp(t) | et | |
0 | ||||
1 | 28 | 28,00578 | -0,00578 | 0,02063387 |
2 | 36 | 36,11451 | -0,11451 | 0,3180925 |
3 | 43 | 43,68993 | -0,68993 | 1,60449915 |
4 | 28 | 27,44259 | 0,557413 | 1,99076144 |
5 | 31 | 30,94963 | 0,050374 | 0,16249833 |
6 | 40 | 39,79552 | 0,204478 | 0,51119446 |
7 | 49 | 47,93757 | 1,062429 | 2,16822255 |
8 | 30 | 30,97227 | -0,97227 | 3,24091111 |
9 | 34 | 34,03861 | -0,03861 | 0,11356266 |
10 | 44 | 43,67856 | 0,321445 | 0,73055605 |
11 | 52 | 52,89722 | -0,89722 | 1,725418 |
12 | 33 | 32,84402 | 0,155981 | 0,47266899 |
13 | 39 | 36,88196 | 2,118044 | 5,43088329 |
14 | 48 | 48,45499 | -0,45499 | 0,94790358 |
15 | 58 | 57,85463 | 0,145375 | 0,25064607 |
16 | 36 | 36,56287 | -0,56287 | 1,56352389 |
Σ | 21,2519759 | |||
ср.знач. | 1,3282485 |
Еотн=1,33<15%, следовательно условие точности выполнено.
Для того, чтобы модель была адекватна исследуемому процессу, ряд остатков et должен обладать свойствами случайности, независимости последовательных уровней, нормальности распределения.
Проверку случайности уровней остаточной компоненты проводим на основе критерия поворотных точек. Для этого каждый уровень ряда et сравниваем с двумя соседними. Если он больше (либо меньше) обоих соседних уровней, то точка считается поворотной и для этой строки ставится 1, в противном случае ставится 0.
t | Y | Yp(t) | E(t) | точки поворота | E2(t) | (E(t)-E(t-1))2 | E(t)*E(t-1) |
0 | |||||||
1 | 28 | 28,00578 | -0,00578 | ххх | 0,0000333793 | ||
2 | 36 | 36,11451 | -0,11451 | 0 | 0,0131132959 | 0,01182348 | 0,000662 |
3 | 43 | 43,68993 | -0,68993 | 1 | 0,4760097997 | 0,33110971 | 0,079007 |
4 | 28 | 27,44259 | 0,557413 | 1 | 0,3107094799 | 1,55587663 | -0,38458 |
5 | 31 | 30,94963 | 0,050374 | 1 | 0,0025375885 | 0,25708826 | 0,028079 |
6 | 40 | 39,79552 | 0,204478 | 0 | 0,0418111646 | 0,02374783 | 0,0103 |
7 | 49 | 47,93757 | 1,062429 | 1 | 1,1287554806 | 0,73608037 | 0,217243 |
8 | 30 | 30,97227 | -0,97227 | 1 | 0,9453154338 | 4,14001378 | -1,03297 |
9 | 34 | 34,03861 | -0,03861 | 0 | 0,0014908328 | 0,87172478 | 0,037541 |
10 | 44 | 43,67856 | 0,321445 | 1 | 0,1033266699 | 0,1296403 | -0,01241 |
11 | 52 | 52,89722 | -0,89722 | 1 | 0,8049989917 | 1,48513712 | -0,28841 |
12 | 33 | 32,84402 | 0,155981 | 0 | 0,0243299999 | 1,1092263 | -0,13995 |
13 | 39 | 36,88196 | 2,118044 | 1 | 4,4861124350 | 3,84969403 | 0,330374 |
14 | 48 | 48,45499 | -0,45499 | 1 | 0,2070192827 | 6,62052558 | -0,9637 |
15 | 58 | 57,85463 | 0,145375 | 1 | 0,0211338100 | 0,36044226 | -0,06614 |
16 | 36 | 36,56287 | -0,56287 | ххх | 0,3168210605 | 0,5016086 | -0,08183 |
Σ | 10 | 8,8835187049 | 21,983739 | -2,26678 |
Общее число поворотных точек p=10. Рассчитаем значение q:
При N=16 – q=6. Так как p>q, условие случайности уровней ряда остатков выполнено.
Проверка независимости уровней ряда остатков проводим двумя методами:
1)
4-d=1,53
Т.к. 1,1<1,37<1,53<2, то уровни ряда остатков являются независимыми.
2) , |r1|=0,255.
Так как |r1|<rтаб=0,32, то уровни ряда остатков независимы.
Нормальность распределения ряда остатков определяем по RS-критерию:
,
Emax= | 2,118044484 |
Emin= | -0,972273333 |
S= | 0,769567788 |
RS= | 4,015653806 |
3<RS<4,21, значит уровни ряда остатков подчиняются нормальному распределению.
Таким образом, все условия адекватности и точности выполнены. Следовательно, можно говорить об удовлетворительном качестве модели и возможности проведения прогноза показателя Yp(t) на четыре квартала вперед.
t | Y | a(t) | b(t) | F(t) | Yp(t) |
0 | 31,7 | 0,9 | 0,785766 | ||
1 | 28 | 32,58132 | 0,868443 | 0,859438 | 28,00578 |
2 | 36 | 33,41794 | 0,858897 | 1,078225 | 36,11451 |
3 | 43 | 34,11445 | 0,810181 | 1,266126 | 43,68993 |
4 | 28 | 35,13745 | 0,874026 | 0,792429 | 27,44259 |
5 | 31 | 36,02906 | 0,879301 | 0,860025 | 30,94963 |
6 | 40 | 36,96525 | 0,896369 | 1,080548 | 39,79552 |
7 | 49 | 38,11336 | 0,97189 | 1,277833 | 47,93757 |
8 | 30 | 38,71716 | 0,861464 | 0,781882 | 30,97227 |
9 | 34 | 39,56516 | 0,857423 | 0,859615 | 34,03861 |
10 | 44 | 40,51182 | 0,884197 | 1,083881 | 43,67856 |
11 | 52 | 41,18538 | 0,821004 | 1,268684 | 52,89722 |
12 | 33 | 42,06623 | 0,838959 | 0,783439 | 32,84402 |
13 | 39 | 43,64437 | 1,060714 | 0,879998 | 36,88196 |
14 | 48 | 44,57915 | 1,022933 | 1,079594 | 48,45499 |
15 | 58 | 45,63646 | 1,033246 | 1,270022 | 57,85463 |
16 | 36 | 46,45417 | 0,968585 | 0,77835 | 36,56287 |
17 | 41,73191 | ||||
18 | 52,24301 | ||||
19 | 62,68817 | ||||
20 | 39,1732 |
Информация о работе Контрольная работа по "Финансовая математика"