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

    Оглавление 
 

    Задание 1 3

    Задание 2 15

    Задание 3 26

    Список использованной литературы 29 
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     
     

    Задание 1

 

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

    Таблица 1

    Исходные  значения заданного временного ряда

Т 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
Y(t) 30 38 45 30 32 42 51 31 36 46 55 34 41 50 60 37
 

    Требуется:

  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 год.
  2. Отразить на графике фактические, расчетные и прогнозные данные.
 

    Решение:

    Мультипликативная модель Хольта-Уинтерса с линейным ростом имеет следующий общий  вид:

                            

                        (1)

    где k – период упреждения;

     - расчетное значение экономического  показателя для t-го периода;

     , , - коэффициенты модели; они адаптируются, уточняются по мере перехода от членов ряда с номером t-1 к t;

     - значение коэффициента сезонности  того периода, для которого  рассчитывается экономический показатель;

    L – период сезонности (для квартальных данных L = 4).

    Уточнение (адаптация к новому значению параметра  времени t) коэффициентов модели производится с помощью формул:

                            

                 (2) - (4)

    Параметры сглаживания α1, α2 и α3 подбирают путем перебора с таким расчетом, чтобы расчетные данные наилучшим образом соответствовали фактическим (т.е. чтобы обеспечить удовлетворительную адекватность и точность модели).

    Из  формул видно, что для расчета  и необходимо оценить значения этих коэффициентов для предыдущего периода времени (т.е для t=1-1=0).

    Значения  и имеют смысл этих же коэффициентов для четвертого квартала года, предшествующего первому году, для которого имеются данные.

    Для оценки начальных значений и применим линейную модель к первым восьми значениям заданного временного ряда (Таблица 2).

    Линейная  модель имеет вид: . Оценим коэффициенты линейной модели и с помощью метода наименьших квадратов:

    

 
 
 
 

    Таблица 2

    Расчёт  коэффициентов линейной модели

i t Y(t) t2 t·Y(t) Yt
1 1 30 1 30 34,667
2 2 38 4 76 35,440
3 3 45 9 135 36,214
4 4 30 16 120 36,988
5 5 32 25 160 37,762
6 6 42 36 252 38,536
7 7 51 49 357 39,310
8 8 31 64 248 40,083
суммы 36 299 204 1378 299,0
 
 

    

    

    

    

    

    Решение системы:

    

    

    Линейная  модель с учетом полученных коэффициентов  имеет вид:

    

    Коэффициент сезонности для 1-го квартала предыдущего  первому года:

    

    Коэффициент сезонности для 2-го квартала предыдущего  первому года:

    

    Коэффициент сезонности для 3-го квартала предыдущего  первому года:

    

    Коэффициент сезонности для 4-го квартала предыдущего  первому года:

    

    Перейдем к построению адаптивной мультипликативной модели Хольта-Уинтерса. Путём перебора возможных значений параметров сглаживания было установлено, что лучшими являются:

    

    Из  уравнения (1), полагая t=0, k=1, находим Y1:

    

    Из  уравнений (2-4), полагая t=1, находим:

    

    

    

    Из  уравнения (1), полагая t=2, Y2:

    

    Из  уравнений (2-4), полагая t=2, находим:

    

      Из уравнения (1), полагая t=3, Y3:

    

    Из  уравнений (2-4), полагая t=3, находим:

    

    Из  уравнения (1), полагая t=4, Y4:

    

    Из  уравнений (2-4), полагая t=4, находим:

    

    Из  уравнения (1), полагая t=5, Y5:

    

    Из  уравнений (2-4), полагая t=5, находим:

    

    Из  уравнения (1), полагая t=6, Y6:

    

    Из  уравнений (2-4), полагая t=6, находим:

    

    Из  уравнения (1), полагая t=7, Y7:

    

    Из  уравнений (2-4), полагая t=7, находим:

    

    Из  уравнения (1), полагая t=8, Y8:

    

    Из  уравнений (2-4), полагая t=8, находим:

    

    Из  уравнения (1), полагая t=9, Y9:

    

    Из  уравнений (2-4), полагая t=9, находим:

    

    Из  уравнения (1), полагая t=10, Y10:

    

    Из  уравнений (2-4), полагая t=10, находим:

    

    Из  уравнения (1), полагая t=11, Y11:

    

    Из  уравнений (2-4), полагая t=11, находим:

    

    Из  уравнения (1), полагая t=12, Y12:

    

    Из  уравнений (2-4), полагая t=12, находим:

    

    Из  уравнения (1), полагая t=13, Y13:

    

    Из  уравнений (2-4), полагая t=13, находим:

    

    Из  уравнения (1), полагая t=14, Y14:

    

    Из  уравнений (2-4), полагая t=14, находим:

    

    Из  уравнения (1), полагая t=15, Y15:

    

    Из  уравнений (2-4), полагая t=15, находим:

    

    Из  уравнения (1), полагая t=16, Y16:

    

    Из  уравнений (2-4), полагая t=16, находим:

    

    Таблица 3

    
Абс.погр.,

Отн.погр. в %

1 2 3 4 5 6 7 8
0 - 33,893 0,774 - - - -
1 30 34,776 0,807 0,8602 29,689 0,311 1,04
2 38 35,453 0,768 1,0755 38,467 -0,467 1,23
3 45 35,984 0,697 1,2583 46,000 -1,000 2,22
4 30 37,037 0,804 0,8029 29,060 0,940 3,13
5 32 37,649 0,746 0,8540 32,549 -0,549 1,72
6 42 38,592 0,805 1,0832 41,295 0,705 1,68
7 51 39,737 0,907 1,2734 49,574 1,426 2,80
8 31 40,034 0,724 0,7858 32,633 -1,633 5,27
9 36 41,176 0,850 0,8662 34,809 1,191 3,31
10 46 42,158 0,889 1,0880 45,523 0,477 1,04
11 55 43,091 0,902 1,2752 54,817 0,183 0,33
12 34 43,776 0,837 0,7803 34,568 -0,568 1,67
13 41 45,430 1,082 0,8880 38,643 2,357 5,75
14 50 46,345 1,032 1,0825 50,602 -0,602 1,20
15 60 47,280 1,003 1,2715 60,415 -0,415 0,69
16 37 48,023 0,925 0,7744 37,676 -0,676 1,83
Сумма - - - - - - 34,90

    Модель  Хольта-Уинтерса 
 

    Модель  Хольта-Уинтерса:

    

    Проверка  качества модели:

    Абсолютная погрешность (должна удовлетворять определённым условиям точности и адекватности).

    Таблица 4

    Промежуточные расчёты для оценки адекватности модели

Точки поворота
1 0,311 - 0,097 - - -
2 -0,467 0 0,218 -0,778 0,605 -0,145
3 -1,000 1 1,000 -0,533 0,285 0,467
4 0,940 1 0,884 1,940 3,764 -0,940
5 -0,549 1 0,302 -1,489 2,218 -0,516
6 0,705 0 0,496 1,254 1,572 -0,387
7 1,426 1 2,032 0,721 0,520 1,004
8 -1,633 1 2,666 -3,058 9,354 -2,328
9 1,191 1 1,419 2,824 7,975 -1,945
10 0,477 0 0,228 -0,714 0,510 0,569
11 0,183 0 0,034 -0,294 0,086 0,088
12 -0,568 1 0,323 -0,751 0,565 -0,104
13 2,357 1 5,554 2,925 8,553 -1,339
14 -0,602 1 0,363 -2,959 8,756 -1,420
15 -0,415 1 0,172 0,187 0,035 0,250
16 -0,676 - 0,457 -0,261 0,068 0,280
Сумма 1,681 10 16,243 -0,987 44,865 -6,466
 

    Проверка  точности модели:

    Условие точности выполнено, если средняя относительная  погрешность  (в табл.  Модель Хольта-Уинтерса в столб.8) в среднем не превышает 5%. У нас условие точности выполнено, т.к.

     %

     %

    Проверка  условия адекватности:

    Для адекватности модели ряд остатков E(t) должен обладать свойствами случайности, независимости последовательных уровней, нормальности распределения.

  • Проверка случайности уровней ряда остатков.

    

    При N = 16:

    

    Р = 10 – число поворотных точек. Так как p > q, то условие случайности выполнено.

  • Проверка независимости уровней ряда остатков (отсутствие автокорреляции:

    а) по d-критерию Дарбина-Уотсона,

    б) по первому коэффициенту автокорреляции r(1).

    а). 

    Так как d > 2, то присутствует отрицательная автокорреляция, тогда d надо уточнить, вычитая полученное из 4:

    

    Так как  , то критерий ответа не дает.

    б). 

    Так как > , то уровни ряда остатков зависимы, в ряду остатков присутствует автокорреляция.

  • Проверка соответствия ряда остатков нормальному распределению по RS – критерию:

    

      = 2,356 - максимальное значение уровней ряда остатков ,

      = -1,633 - минимальное значение уровней ряда остатков ,

     - среднеквадратическое  отклонение. 

    

    Так как RS попадает в интервал от 3 до 4,21, то уровни ряда  подчиняются нормальному распределению.

    Таким образом, почти все условия адекватности и точности выполнены. Т.е. модель удовлетворительна и можно проводить прогноз показателя Y(t) на 4 квартала вперед.

    Расчёт  прогнозных значений экономического показателя (на год вперёд, с t=17 по t=20):

    с помощью а(16) и b(16) (см. табл. Модель Хольта-Уинтерса) по формуле (1) найдём прогнозные Y(t):

    

    

    Рис.1. График фактических, расчетных и  прогнозных данных по модели Хольта-Уинтерса 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

    Задание 2

 

    Даны  цены (открытия, максимальная, минимальная  и закрытия) за 10 дней. Интервал сглаживания  принять равным пяти дням. Рассчитать:

    - экспоненциальную  скользящую среднюю;

    - момент;

    - скорость  изменения цен;

    - индекс  относительной силы;

    - %R, %K, %D.

    Расчеты проводить для всех дней, для которых  эти расчеты можно выполнить на основании имеющихся данных.

    Таблица 5

    Исходные  данные

Дни Цены
Макс. Мин. Закр.
1 765 685 750
2 792 703 733
3 740 706 733
4 718 641 666
5 680 600 640
6 693 638 676
7 655 500 654
8 695 630 655
9 700 640 693
10 755 686 750
 

    Решение:

    Рассчитаем  экспоненциальную скользящую среднюю:

    

    где - значение экспоненциальной скользящей средней текущего дня t;

     - цена закрытия t- го дня,

     - коэффициент. 

    Интервал  сглаживания n=5.

    Тогда коэффициент К будет равен:

    Вычислим  простую среднюю  для первых  5 дней:

    

    

    

    

    

    

    

    Таблица 6

Дни Цена закрытия EMAt MOMt ROCt, % Повышение цены Понижение цены Сумма повышений Сумма понижений RSI
1 2 3 4 5 6 7 8 9 10
1 750       -      
2 733         17      
3 733       0        
4 666         67      
5 640 704,40 -110 85,33   26 0 110 0,00
6 676 694,93 -57 92,22 36   36 110 24,66
7 654 681,29 -79 89,22   22 36 115 23,84
8 655 672,53 -11 98,35 1   37 115 24,34
9 693 679,35 53 108,28 38   75 48 60,98
10 750 702,90 74 110,95 57   75 22 85,71
Контрольная работа по "Финансовой математике". 19