Контрольная работа по "Финансовая математика"
Контрольная работа, 09 Января 2012, автор: пользователь скрыл имя
Краткое описание
Требуется:
Построить адаптивную мультипликативную модель Хольта-Уинтерса с учетом сезонного фактора, приняв параметры сглаживания α1=0,3; α2=0,6; α3=0,3.
Оценить точность построенной модели с использованием средней относительной ошибки аппроксимации.
Файлы: 1 файл
Финансовая математика_контрольная Вариант 1.doc
— 285.00 Кб (Скачать)В каждом варианте приведены поквартальные данные о кредитах от коммерческого банка на жилищное строительство (в условных единицах) за 4 года (всего 16 кварталов, первая строка соответствует первому кварталу первого года).
Требуется:
- Построить адаптивную мультипликативную модель Хольта-Уинтерса с учетом сезонного фактора, приняв параметры сглаживания α1=0,3; α2=0,6; α3=0,3.
- Оценить точность построенной модели с использованием средней относительной ошибки аппроксимации.
- Оценить адекватность построенной модели на основе исследования:
- случайности остаточной компоненты по критерию пиков;
- независимости уровней ряда остатков по d-критерию (критические значения d1=1,10 и d2=1,37) и по первому коэффициенту автокорреляции при критическом значении r1=0,32;
- нормальности распределения остаточной компоненты по R/S-критерию с критическими значениями от 3 до 4,21.
- Построить точечный прогноз на 4 шага вперед, т.е. на 1 год.
- Отразить на графике фактические, расчетные и прогнозные данные.
Решение:
Мультипликативная модель Хольта-Уинтерса с линейным ростом имеет следующий вид:
Коэффициенты модели рассчитываются по формулам:
Для оценки начальных значений 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, условие случайности уровней ряда остатков выполнено.
Проверка независимости уровней ряда остатков проводим двумя методами:
- по d-критерию Дарбина-Уотсона;
- по первому коэффициенту автокорреляции r1.
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 |