Полиномиальные модели курса акций
В табл. 3.1 представлен временной ряд xt курса акций компании IBM (в долл.), включающий n =30 наблюдений. Приняв коэффициент адаптации a=0, 5 и период упреждения t=1, требуется аппроксимировать ряд xt с помощью адаптивной полиномиальной модели: а) нулевого порядка (p =0); б) первого порядка (p =1); в) второго порядка (p =2); г) рассчитать прогнозные значения при t =31; д) сравнить адекватность моделей по величине остаточного среднеквадратического отклонения . Таблица 3.1. Курс акций компаний IBM (в долл.) xt и расчеты по адаптивным полиномиальным моделям нулевого (p =0), первого (p =1) и второго (p =2) порядков
Решение. а) Адаптивная полиномиальная модель нулевого порядка (p=0). экспоненциальная средняя имеет вид S t = a xt + bSt-1, b = 1 - a. Начальное условие: где за примем среднее значение, например, первых пяти наблюдений. Для нашего примера Расчетное модельное значение с периодом упреждения t будем определять из соотношения: (3.1) Как уже отмечалось: a=0, 5; t=1. Тогда при t =1 S1 = a× x 1 + (1-a)S0 = 0, 5× 510 + 0, 5× 506 = 508 , при t =2 S1 = 0, 5× 497 + 0, 5× 508 = 502, 5 , при t =3 S1 = 0, 5× 504 + 0, 5× 502, 5 = 503, 65 . Остальные результаты расчетов где t =5, 6,..., 31 получены аналогично и приведены в табл. 3.1 (для модели (p =0)).
б) Адаптивная полиномиальная модель первого порядка (p =1). Первоначально по данным временного ряда xt находим МНК-оценку линейного тренда: и принимаем и Экспоненциальные средние 1-го и 2-го порядков определяются как: (3.2) где b=1‑ a. Отсюда начальные условия: (3.3) Оценка модельного значения ряда с периодом упреждения t равна: (3.4) Для нашего примера при a=0, 5 и t=1 имеем: (3.5) МНК-оценка уравнения тренда имеет вид: откуда и Согласно (3.3) получим: . С учетом (3.5) найдем: При t =1 экспоненциальные средние равны: Откуда: При t =2 имеем: . Остальные результаты расчетов где t =4, 5,..., 31 получены аналогично и приведены в табл. 3.1 (для модели p =1).
в) Адаптивная полиномиальная модель второго порядка (p =2). По данным временного ряда xt находим МНК—оценку параболического тренда и экспоненциальные средние 1-го, 2-го и 3-го порядков: Начальные условия определим как: (3.7) где и . Оценка модельного значения с периодом упреждения t находим из выражения (3.8) Для нашего примера, когда a=b=0, 5 и t=1, выражения (3.7) и (3.8) преобразуются к виду МНК-оценка уравнения тренда равна откуда: . Согласно (3.9) получим начальные условия: Из (3.10) следует: При t =1согласно (3.6) имеем откуда из (3.10): При t =2 с учетом (3.6): откуда согласно (3.10): . При t =3 с учетом (3.6) Остальные результаты расчетов при t =4, 5,..., 30 получены аналогично и приведены в табл. 3.1 (для модели p =2).
г) Прогнозные значения по адаптивным полиномиальным моделям равны (табл. 3.1): как видим они практически совпадают.
д) остаточные среднеквадратические отклонения для моделей соответственно равны: 6, 3; 7, 4 и 8, 1. Из этого следует, что адаптивная полиномиальная модель нулевого порядка наиболее адекватна процессу xt. Рис. 3.1 Исходный (xt) и модельный ряд курса акций IBM (долл.). Адаптивная полиномиальная модель 0-го порядка(p =0).
Рис. 3.2 Исходный (xt) и модельный ряд курса акций IBM (долл.). Адаптивная полиномиальная модель 1-го порядка(p =1). Рис. 3.3 Исходный (xt) и модельный ряд курса акций IBM (долл.). Адаптивная полиномиальная модель 2-го порядка(p =2).
Однако сопоставляя исходные хt и модельные значения временного ряда (рис. 3.1, 3.2 и 3.3) можно предположить, что прогностические свойства адаптивных полиномиальных моделей 1-го и 2-го порядка лучше, чем нулевого. Это следует из близости значений и в последних пяти точках временного ряда (при t =26¸ 30). На этом интервале остаточные среднеквадратические отклонения s для полиномиальных моделей нулевого (p =0), первого (p =1) и второго (p =2) порядков соответственно равны 6, 36; 4, 49 и 4, 45. В этой связи предпочтение следует отдать полиномиальной модели первого порядка, т.к. она проще и практически не отличается по точности от модели второго порядка.
|