Контрольная работа по "Финансовая математика"

Автор: Пользователь скрыл имя, 09 Января 2012 в 12:42, контрольная работа

Краткое описание

Требуется:
Построить адаптивную мультипликативную модель Хольта-Уинтерса с учетом сезонного фактора, приняв параметры сглаживания α1=0,3; α2=0,6; α3=0,3.
Оценить точность построенной модели с использованием средней относительной ошибки аппроксимации.

Файлы: 1 файл

Финансовая математика_контрольная Вариант 1.doc

— 285.00 Кб (Скачать)

Задание 1.

В каждом варианте приведены поквартальные  данные о кредитах от коммерческого  банка на жилищное строительство (в  условных единицах) за 4 года (всего 16 кварталов, первая строка соответствует первому  кварталу первого года).

Требуется:

  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 шага вперед, т.е. на 1 год.
  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, условие случайности уровней ряда остатков выполнено.

    Проверка  независимости уровней ряда остатков проводим двумя методами:

  1. по d-критерию Дарбина-Уотсона;
  2. по первому коэффициенту автокорреляции 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

Информация о работе Контрольная работа по "Финансовая математика"