Стохастическая модель динамики относительных приращений цены акции Текст научной статьи по специальности «Математика»

ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА

2011 Математика и механика № 2(14)

УДК 519.23

П.В. Трясучёв СТОХАСТИЧЕСКАЯ МОДЕЛЬ ДИНАМИКИ ОТНОСИТЕЛЬНЫХ ПРИРАЩЕНИЙ ЦЕНЫ АКЦИИ

В работе будет рассматриваться процесс, описывающий относительные приращения цены акции с помощью обобщённого уравнения Ито. Для описания стохастической динамики были взяты цены акции компании Лукойл за период с 18.04.2008 по 17.04.2009, с интервалами Ат = 1 мин, 5 мин, 10 мин, 15 мин, 30 мин и 1 час.

Ключевые слова: стохастический процесс, коэффициент сноса, волатильность, относительные приращения, виннеровский процесс, марковский процесс.

При описании стоимости ценных бумаг на финансовых рынках широко используется модель геометрического броуновского движения [1, 2]. Согласно этой модели, стоимость актива как функция времени ґ подчиняется стохастическому дифференциальному уравнению Ито вида

где постоянные д и с соответственно дрейф и волатильность; - стандартный

винеровский процесс.

В настоящее время становится ясно, что такая модель не является состоятельной [3] и не учитывает многие важные особенности рынка [4, 5]. Чтобы учесть различного рода наблюдаемые закономерности, были введены модели со стохастической волатильностью [1, 2, 6 - 9], согласно которым волатильность рассматривается как случайная переменная и в общем случае как функция с = с (Т(ґ)) некоторого стохастического процесса У(ґ).

В настоящей работе будет рассмотрен стохастический процесс Я(ґ) = (+! - )/ , описывающий относительные приращения цен акций. Суще-

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

Рассмотрим обобщённую модель Ито [1, 2, 8, 9]

= ц(Я, ґ)йґ + ст(ґ)йШ{, (1)

где на функции коэффициента дрейфа и диффузии выдвигаются стандартные условия

-0(ц(Я, ґ)) <ж, ґ ^ ж, ст(ґ) < ж, ґ ^ ж.

В этой связи представляет интерес получение оценок коэффициентов модели стохастической волатильности непосредственно из эмпирических данных, без привлечения метода наименьших квадратов или максимального правдоподобия.

Анализ динамики относительных приращений для акции

Для описания стохастической динамики были взяты цены акции компании Лу-коил за период с 18.04.2008 по 17.04.2009, с интервалами т = 1 мин, 5 мин, 10 мин, 15 мин, 30 мин и 1 ч, и посчитаны относительные приращения для всех наборов данных. Как видно из рис. 1, плотность распределения для близка к нормальной.

8000

«

* 6000

0

1

иЗ Й

® 4000

м н

о

о

5

§ 2000 0

-0,02 -0,01 0 0,01 Я

Рис. 1. Плотность распределения относительных приращений для т = 5 мин. Сплошной линией показана теоретическая плотность распределения, столбцы показывают реальное распределение

Величины Я1 будут зависимыми случайными величинами. Данное обстоятельство существенно усложняет процесс эконометрического анализа, так как делает необходимым рассмотрение совместной многомерной плотности распределения ры (Я1, т1; Я2, т2;...; Яы, %). Тем не менее, при дополнительных предположениях относительно Яь выдвинутых в [10], можно рассматривать данные величины как независимые.

Тем не менее, при дополнительных предположениях относительно Я( удается перейти к последовательности независимых случайных величин (СВ).

Выберем последовательность независимых одинаково распределенных СВ Е |^| < ж, таких, что

Я' =^0 +^1 +... + ь, Я = £„, < > 0.

Известно [11], что она будет мартингал-разностью на семействе / = ст(ю; §!,..., ^ }, те. почти наверное (п.н.) выполнено следующее равенство:

Е &+1 / ) = 0.

Как следствие, Я( будет мартингалом относительно этого семейства с-алгебр /. Поэтому для перехода к исследованию только независимых приращений дос-

таточно найти такую ) 0. Например, в [12] показано, что в случае существования условного среднего а = Е (і^|^) для процесса (Яі ) 0 последний представим в виде процесса

ад

і = а + ХСг^- ,

і=0

где

ад / ч

Сі ; 2|Сг| <ад; а = Е(і|^і) - условное среднее; (^)^0 ~ N(0,ст,)

і=0

- по-

следовательность независимых, нормально распределенных СВ тогда и только тогда, когда а+ р< 1. Поэтому всюду далее будем предполагать, что Я( независимы, переходя, при необходимости, к рассмотрению последовательности ) .

Стохастический процесс Я(т) полностью определяется бесконечным набором совместных плотностей ры(Я1,т1;Я2,т2;...;Яы,), зависящих от N переменных. Существенное упрощение возникает, если Я(т) есть марковский процесс. В этом случае ^точечная плотность распадается на произведение условных плотностей:

р(Я,,Т,. |Яг+1,хг+1) = Р2(Я(,Я;^Т)+1),

Р1( Я+и Т,+1)

где , = 1, 2, ... , N-1 и р(Я,, т, | Я,+1, т,-+1 ) обозначает плотность условной вероятности реализации значения Я, за время т,- при заданном значении Я,+1 за время т,+1. Будем считать, что т,+1 > т,.

Как известно, в случае марковского процесса условные плотности должны удовлетворять уравнению Чэпмена - Колмогорова

р(Я1,т1 | Я2,т2) = |р(Я1,т1 | Я,т)р(Я,т | Я2,т2 )аЯ , (2)

где

р(^ті |R2,Т2) =

(3)

Р(Я2,Т2)

Анализируя исходные данные, легко найти условные плотности распределения; проводя численное интегрирование, можно получить следующие результаты, которые представлены на рис. 2 и 3. Из рис. 3 видно, что процесс Я - марковский случайный процесс.

Рис. 2. Двумерная плотность распределения р(Я1,т1|Я2,т2), где т = 10 мин, т2 = 15 мин

Х^і,Ті|^2,Т2) 0,6

0,4

-0,2

-0,1

і

і

■Ю,2

1

і

1

і

\

\

Ч

0

0,1

R

Рис. 3. Условная плотность распределения р(Я1,т1|Я2,т2), Я! = -0,022.

Сплошной линией показана условная плотность, вычисленная по формуле (3), точками - по формуле (2)

Таким образом, с учётом вышесказанного, будем считать, что процесс Я подчиняется уравнению (1), где ц(Я,/) - коэффициент сноса, с(0 - волатильность, а ёШ - приращения для стандартного винеровского процесса.

Коэффициент ц(Я,0 найдём как коэффициент Б1(Я,() уравнения Фоккера-Планка, которое соответствует дифференциальному уравнению (1):

dp(R, t) 1 d

dt

, 2D2(R,t)P(R,t)Di(R,t)p(R,t),

2 dR2

d_

dR

где

D—(R,t) = lim— f (R— - R2 )p(R—,t + Ax | R2,t)dR1.

X^Ü T R

Численные значения коэффициента D—(R,t) представлены на рис. 4, из которого видно, что коэффициент сноса линеен.

Di(R,t)

Ü 1Ü 2Ü 3Ü 4Ü R Рис. 4. Численные значения коэффициента сноса

Как не раз оказывалось на практике, волатильность является стохастическим процессом, который описывается следующим уравнением [13]:

йст = а(ст)йї + Р(ст)йТ , (4)

где а(о) - коэффициент сноса для волатильности; Р(о) - волатильность волатильности, а йТ - приращения винеровского случайного процесса. Будем полагать процесс о(ї) стационарным, поэтому параметры уравнения (4) не зависят от времени.

Будем искать волатитльность волатильности в виде vоr [13], для этого возведём левую и правую часть (4) в квадрат и получим

(йст)2 = р2(ст)йї ,

так как (йТ)2 = йї, (йї)2 = 0, йї-йТ = 0. Из имеющихся значений волатильности для всех шести временных серий найдём (До)2, а также среднее значение М[(До)2]. Таким образом, можно получить следующую зависимость:

М[(Дст)2 ] = р2(ст)Дї .

Положим Р(о) = vоY, а также напомним, что Д = т. Затем, логарифмируя, получим

1п(М[(Дст)2 ]) = 1п^2т) + 2у 1п(ст).

На основе линейной регрессионной модели находим

V = 6,178; у = 0,509, результаты представлены на рис. 5.

Таким образом, Р(с) = voy.

Чтобы определить а(с), рассмотрим плотность вероятности р(с, /). Плотность этой вероятности определяет уравнение Фоккера-Планка

$Р = ^(руи-^ар). (5)

д/ 2 дст2' ' да ’

Теперь предположим, что нам известна плотность вероятности р(с, /) в стационарном состоянии, обозначим её Роо(о), тогда уравнение (5) перепишется следующим образом:

0 Ч^21-)-^40. <б)

Проинтегрировав один раз уравнение (6), получаем выражение для коэффициента сноса волатильности

а(ст) = ———(2рш) .

^ 2Р- дст' '

Распределение рю(а), как видно из рис. 6, логнормальное. Тогда аналитическое выражение для данной плотности примет вид

Ры =■

1

~е 2а2

урЫаст

где 1пСТ характеризует математическое ожидание 1пст, а а - дисперсию. Графики, полученные для а(о), можно увидеть на рис. 7.

80

Ж

1 60

2

иЗ

л

X

§ 40

20

1,7199-10-6 0,0002 0,0004 0,0005 0,0007 0,0009

Рис. 6. Плотность распределения волатильности (т=15мин)

1п2[СТ

ст

Рис. 7. Коэффициент сноса для волатильности. На графике снизу-вверх представлены коэффициенты сноса для т = 1, 5, 10, 15, 30, 60 мин

В настоящее время знания о случайном и детерминированном характере процессов, лежащих в основе эволюции финансовых рынков, в значительной мере ограничены. При условии, что относительные ценовые приращения R представляют собой марковский процесс на временной шкале t, было показано, что непосредственно из эмпирических данных могут быть получены коэффициенты обобщённого уравнения Ито (1). В частности, были посчитаны коэффициенты а и в для процесса стохастической волатильности, что позволяет строить краткосрочные прогнозы и доверительные интервалы для волатильности, которые широко используются для оценивания опционов и прогнозирования цены активов.

ЛИТЕРАТУРА

1. Hull J., White A. The Pricing of Options on Assets with Stochastic Volatilities // J. Finance. 1987. V. XLII. P. 281-300.

2. Dragulescu A., Yakovenko V.M. Probability distribution of returns in the Heston model with stochastic volatility // Quant. Finance. 2002. V. 2. P. 443-453.

3. ContR. Empirical properties of asset returns: stylized facts and statistical issues // Quant. Finance. 2001. V. 1. P.223-236.

4. Friedrich R., Peinke J. .Renner Ch. How to quantify deterministic and random influences on the statistics of the foreign exchange market // Phys. Rev. Lett. 2000. V. 84. P. 5224-5227.

5. Ivanova K.,Ausloos M., and Takayasu H. Deterministic and stochastic influences on Japan and US stock and foreign exchange markets. A Fokker-Planck approach //arXiv:cond-mat/0301268.

6. Stein E.M. and Stein J.C. Stock Price Distributions with Stochastic Volatility: An Analytic Approach // Rev. Financial Studies. 1991. V. 4. P. 727-752.

7. Heston S.L. A closed-form solution for options with stochastic volatility with applications to bond and currency options // Rev. Financial Studies. 1993. V. 6. P. 327-343.

8. Бухбиндер Г.Л., Чистилин К.М. Описание российского фондового рынка в рамках модели Гестона // Математическое моделирование. 2005. Т. 17. № 10. С. 31-38.

9. Remer R., Mahnke R. Application of Heston model and its solution to German DAX data // Physica A. 2004. V. 344. P. 236-239.

10. ШиряевА.Н. Вероятность. М.: Наука, 1980. 574 с.

11. McNeil A.J., Frey R., Embrechts P. Quantitative Risk Management. Concepts, Techniques and Tools. Princeton, NJ: Princeton University Press, 2005.

12. Бухбиндер Г.Л., Чистилин К.М. Стохастическая динамика котировок акций РАО ЕЭС // Математическое моделирование. 2005. Т. 17. № 2. С. 119-125.

13. Wilmott P., Oztukel A. Uncertain parameters, an empirical stochastic colatility model and confidance limits // Int. J. Theor. Appl. Fin. 1998. V. 1. P. 175-198.

Статья поступила 01.12.2010 г.

Tryasuchev P.V. STOCHASTIC MODEL OF DYNAMIC RELATIVE INCREMENTS STOCK PRICE. In this paper, the process of relative increment of stock price is considered. The process is described using the generalized Ito equation. Stochastic dynamics was described with Lukoil stock prices during the period of 18.04.2008 up to 17.04.2009, with intervals Дт = 1 min, 5 min, 10 min, 15 min, 30 min, and 60 min.

Keywords: Stochastic process, drift, volatility, relative increments, Wiener process, Markov process.

TRYASUCHEV Petr Vladimirovich (Tomsk Polytechnic University)

E-mail: pet3001@yandex.ru