Точные решения уравнения Ито — различия между версиями

Материал из synset
Перейти к: навигация, поиск
 
(не показаны 3 промежуточные версии этого же участника)
Строка 5: Строка 5:
 
|}
 
|}
 
----
 
----
<math>\textstyle \bullet</math> Несмотря на простой вид, стохастические уравнения () аналитически интегрировать не так и просто из-за винеровского члена <math>\textstyle \delta W</math>. Это явно видно в случае конечной численной реализации (). Каждое последовательное <math>\textstyle x</math> в итерационной процедуре нелинейным образом зависит от всех предыдущих случайных чисел <math>\textstyle \varepsilon_k</math> (<math>\textstyle \lessdot</math> C). Тем не менее, рассмотрим ситуации, в которых можно получить точные решения.
+
 
 +
<math>\textstyle \bullet</math> Несмотря на простой вид, [[Уравнения Ито|стохастические уравнения (2.4)]] аналитически интегрировать не так и просто из-за винеровского члена <math>\textstyle \delta W</math>. Это явно видно в случае [[Уравнения Ито|конечной численной реализации (2.5)]]. Каждое последовательное <math>\textstyle x</math> в итерационной процедуре нелинейным образом зависит от всех предыдущих случайных чисел <math>\textstyle \varepsilon_k</math> (<math>\textstyle \lessdot</math> C). Тем не менее, рассмотрим ситуации, в которых можно получить точные решения.
  
 
Пусть в процессе Ито функции сноса и волатильности зависят только от времени. Обозначим их через <math>\textstyle f(t)</math> и <math>\textstyle s(t)</math>:
 
Пусть в процессе Ито функции сноса и волатильности зависят только от времени. Обозначим их через <math>\textstyle f(t)</math> и <math>\textstyle s(t)</math>:
  
:<center><math> dx = f(t)\,dt + s(t)\,\delta W. </math></center>
+
{| width="100%"
 +
| width="90%" align="center"|<math> dx = f(t)\,dt + s(t)\,\delta W. </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(2.17)'''</div>
 +
|}
  
Это уравнение легко интегрируется при помощи "дискретной" интерпретации стохастического члена <math>\textstyle \delta W</math>. Рассмотрим итерации, выполняемые по разностной схеме ():
+
Это уравнение легко интегрируется при помощи "дискретной" интерпретации стохастического члена <math>\textstyle \delta W</math>. Рассмотрим итерации, выполняемые по [[Уравнения Ито|разностной схеме (2.5)]]:
  
 
:<center><math>\begin{array}{l} x_1 = x_0 + f_0 \Delta t + s_0 \varepsilon_1 \sqrt{\Delta t},\\ x_2 = x_1 + f_1 \Delta t + s_1 \varepsilon_2 \sqrt{\Delta t} = x_0 + (f_0+f_1) \,\Delta t + (s_0 \varepsilon_1+s_1 \varepsilon_2) \sqrt{\Delta t},\\ ..., \end{array}</math></center>
 
:<center><math>\begin{array}{l} x_1 = x_0 + f_0 \Delta t + s_0 \varepsilon_1 \sqrt{\Delta t},\\ x_2 = x_1 + f_1 \Delta t + s_1 \varepsilon_2 \sqrt{\Delta t} = x_0 + (f_0+f_1) \,\Delta t + (s_0 \varepsilon_1+s_1 \varepsilon_2) \sqrt{\Delta t},\\ ..., \end{array}</math></center>
Строка 21: Строка 25:
 
Скобки в последнем слагаемом содержат сумму независимых гауссовых чисел, каждое из которых имеет волатильность <math>\textstyle s_k</math>. В результате получается гауссово число с волатильностью <math>\textstyle \sqrt{s^2_0+...+s^2_{n-1}}</math>. Поэтому, переходя к непрерывному пределу, получаем (<math>\textstyle \lessdot</math> H):
 
Скобки в последнем слагаемом содержат сумму независимых гауссовых чисел, каждое из которых имеет волатильность <math>\textstyle s_k</math>. В результате получается гауссово число с волатильностью <math>\textstyle \sqrt{s^2_0+...+s^2_{n-1}}</math>. Поэтому, переходя к непрерывному пределу, получаем (<math>\textstyle \lessdot</math> H):
  
:<center><math> x(t) = x(t_0) + \int\limits^t_{t_0} f(\tau) \,d\tau + \left[ \int\limits^t_{t_0} s^2(\tau) \,d\tau \right]^{1/2} \cdot \;\varepsilon. </math></center>
+
{| width="100%"
 +
| width="90%" align="center"|<math> x(t) = x(t_0) + \int\limits^t_{t_0} f(\tau) \,d\tau + \left[ \int\limits^t_{t_0} s^2(\tau) \,d\tau \right]^{1/2} \cdot \;\varepsilon. </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(2.18)'''</div>
 +
|}
  
Решение () уравнения () говорит нам, что <math>\textstyle x(t)</math> является нормально распределённым случайным числом со средним и дисперсией, зависящими от времени. Если <math>\textstyle s(t)</math> &mdash; не константа, то будущая неопределённость в значении <math>\textstyle x</math> может увеличиваться уже не как <math>\textstyle \sqrt{t}</math>, а по другому закону.
+
Решение (2.18) уравнения (2.17) говорит нам, что <math>\textstyle x(t)</math> является нормально распределённым случайным числом со средним и дисперсией, зависящими от времени. Если <math>\textstyle s(t)</math> &mdash; не константа, то будущая неопределённость в значении <math>\textstyle x</math> может увеличиваться уже не как <math>\textstyle \sqrt{t}</math>, а по другому закону.
  
Соотношение () позволяет легко вычислить статистические свойства процесса, в частности, его среднее <math>\textstyle \overline{x}(t)</math> и волатильность <math>\textstyle \sigma(t)</math>.
+
Соотношение (2.18) позволяет легко вычислить статистические свойства процесса, в частности, его среднее <math>\textstyle \overline{x}(t)</math> и волатильность <math>\textstyle \sigma(t)</math>.
  
 
<math>\textstyle \bullet</math> Одномерное уравнение Ито для процесса, имеющего произвольный снос <math>\textstyle a(x,t)</math> и волатильность <math>\textstyle b(x,t)</math>
 
<math>\textstyle \bullet</math> Одномерное уравнение Ито для процесса, имеющего произвольный снос <math>\textstyle a(x,t)</math> и волатильность <math>\textstyle b(x,t)</math>
  
:<center><math> dx = a(x,t)\,dt + b(x,t)\,\delta W, </math></center>
+
{| width="100%"
 
+
| width="90%" align="center"|<math> dx = a(x,t)\,dt + b(x,t)\,\delta W, </math>
заменой иногда можно свести к частному случаю (), для которого решение уже известно. Воспользуемся леммой Ито:
+
| <div width="10%" align="right" style="color:#0000CC">'''(2.19)'''</div>
 
+
|}
:<center><math> dF = \underbrace{\left(\frac{\partial F}{\partial t} \,+\, a(x,t) \frac{\partial F}{\partial x} \,+\, \frac{b^2(x,t)}{2}\, \frac{\partial^2 F}{\partial x^2}\right)}_{f(t)}\, dt + \underbrace{b(x,t)\frac{\partial F}{\partial x}}_{s(t)}\, \delta W. </math></center>
 
 
 
Подберем <math>\textstyle F(x,t)</math> таким образом, чтобы множители при <math>\textstyle \delta W</math> и <math>\textstyle dt</math> в () оказались функциями <math>\textstyle s(t)</math> и <math>\textstyle f(t)</math>, зависящими только от времени:
 
 
 
:<center><math> \frac{\partial F}{\partial x} = \frac{s(t)}{b(x,t)},\;\;\;\;\;\;\;\;\;\;\;\;\; \frac{\partial F}{\partial t} + s(t)\cdot \left[\frac{a(x,t)}{b(x,t)}-\frac{1}{2}\frac{\partial b(x,t)}{\partial x} \right] = f(t), </math></center>
 
 
 
где вместо <math>\textstyle \partial F/\partial x</math> в множитель при <math>\textstyle dt</math> подставлено первое уравнение () и его производная по <math>\textstyle x</math> (<math>\textstyle \lessdot</math> H). Возьмём частные производные первого уравнения () по <math>\textstyle t</math> и второго по <math>\textstyle x</math>. Вычитая их, мы придём к ''условию совместности'':
 
 
 
:<center><math> \frac{1}{s(t)}\cdot \frac{\partial }{\partial t}\left\{\frac{s(t)}{b(x,t)} \right\} = \frac{1}{2}\frac{\partial^2 b(x,t)}{\partial x^2} - \frac{\partial }{\partial x}\left\{\frac{a(x,t)}{b(x,t)} \right\}. </math></center>
 
 
 
''Если'' при данных <math>\textstyle a(x,t)</math> и <math>\textstyle b(x,t)</math> можно подобрать такую функцию <math>\textstyle s(t)</math>, при которой уравнение () обратится в тождество, то мы получим решение стохастического уравнения () в следующей неявной форме:
 
 
 
:<center><math> F\bigl(x(t), t\bigr) = F\bigl(x(t_0), t_0\bigr) + \int\limits^t_{t_0} f(\tau) \,d\tau + \left[ \int\limits^t_{t_0} s^2(\tau) \,d\tau \right]^{1/2} \cdot \;\varepsilon, </math></center>
 
 
 
где функция <math>\textstyle f(t)</math> определяется вторым соотношением (), а <math>\textstyle F(x,t)</math> находится из первого уравнения () (<math>\textstyle \lessdot</math> C).
 
 
 
Решение () &mdash; это нестационарный гауссовый процесс для деформации <math>\textstyle x(t)</math> при помощи нелинейной функции <math>\textstyle F(x,t)</math>. Естественно, что разрешимость () позволяет интегрировать уравнение Ито только в ряде частных случаев. Однако, как мы увидим ниже, эти случаи охватывают достаточно широкий класс важных для практических приложений процессов.
 
 
 
==Простые стохастические модели==
 
 
 
<math>\textstyle \bullet</math> ''Логарифмическое блуждание'' определяется уравнением:
 
 
 
:<center><math> { \;\;dx = \mu \,x\,dt + \sigma\,x\,\delta W\; }, </math></center>
 
 
 
где <math>\textstyle \mu</math> и <math>\textstyle \sigma</math> &mdash; константы модели. Часто () называют геометрическим или экспоненциальным броуновским блужданием.
 
 
 
Если стохастического члена нет (<math>\textstyle \sigma=0</math>), то это обычное уравнение экспоненциального роста (<math>\textstyle \mu>0</math>) или снижения (<math>\textstyle \mu<0</math>):
 
 
 
:<center><math>\frac{dx}{dt} = \mu\,x\;\;\;\;\;\;\;\;\;\;\;\;=>\;\;\;\;\;\;\;\;\;\;\;\;\;\;x(t)=x_0\,e^{\mu t}.</math></center>
 
 
 
Подобная зависимость возникает во многих физических, биологических и социальных системах, от радиоактивного распада до роста экономики.
 
 
 
Случайное воздействие вносит в гладкую динамику определённые коррективы. Подставим функции сноса <math>\textstyle a(x,t)=\mu\,x</math> и волатильности <math>\textstyle b(x,t)=\sigma\,x</math> в условие совместности () на стр. \pageref{ito_main_def}. В результате для <math>\textstyle s(t)</math> получается тривиальное уравнение <math>\textstyle \dot{s}(t)=0</math>, где точка сверху обозначает производную по времени. Следовательно, <math>\textstyle s(t)</math> &mdash; это константа, которую удобно выбрать равной <math>\textstyle \sigma</math>. Интегрирование первого уравнения () даёт <math>\textstyle F(x,t)=\ln x</math>, и, соответственно, функция <math>\textstyle f(t)</math> равна <math>\textstyle \mu-\sigma^2/2</math>. Окончательное решение (<math>\textstyle t_0=0</math>) имеет вид:
 
 
 
:<center><math> x(t)= x_0\, e^{\left(\mu -\sigma^2/2 \right)\, t + \sigma\sqrt{t}\, \varepsilon}. </math></center>
 
 
 
Если в процессе Винера <math>\textstyle x</math> может "уползти" при блуждании в область отрицательных значений <math>\textstyle x<0</math>, то для логарифмической модели это невозможно. Подобное свойство можно было ожидать сразу по виду (). По мере приближения к значению <math>\textstyle x=0</math> снос и волатильность уменьшаются. В результате динамика как бы замораживается при <math>\textstyle x\to 0</math>.
 
 
 
Используя интеграл () на стр. \pageref{aver_exp_gauss}, легко вычислить среднее значение и волатильность в произвольный момент времени:
 
 
 
:<center><math>\bar{x}(t)=x_0\,e^{\mu t},\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\sigma_x(t)=\bar{x}(t) \cdot \sqrt{e^{\sigma^2 t}-1} .</math></center>
 
 
 
Заметим, что необходимо решительно бороться с искушением "по обычному" обращаться со стохастическими уравнениями. Например, разделив () на <math>\textstyle x</math>, нельзя внести его под дифференциал: <math>\textstyle dx/x \neq d\ln x</math>. Для подобных действий служит лемма Ито () по которой для процесса логарифмического блуждания <math>\textstyle d(\ln x)=(\mu-\sigma^2/2)\,dt + \sigma \, \delta W</math>. Фактически, при помощи этой замены, найденной по алгоритму стр. \pageref{ito_main_def}, мы и получили решение ().
 
 
 
Ниже на левом рисунке приведены логарифмические блуждания с нулевым сносом: <math>\textstyle dx = x\delta W</math>. Видно, как они, прижимаясь к <math>\textstyle x=0</math>, тем не менее, остаются в положительной области. В результате получается несимметричное распределение для <math>\textstyle x</math>, которое в данном случае имеет логнормальный вид. Справа динамика дополнена детерминированным сносом: <math>\textstyle dx = 0.05\cdot x\cdot (dt + \delta W)</math>. Она имеет ярко выраженный экспоненциальный рост со стохастическими колебаниями вокруг экспоненты.
 
 
 
[[File:log_winer.png]]
 
 
 
\\ Эти два примера напоминают нам, что стохастические процессы могут быть как малыми поправками к детерминированной динамике (справа), так и основной сутью исследуемой системы (слева).
 
 
 
Введя винеровский процесс <math>\textstyle W_t=W(t)=\varepsilon\sqrt{t}</math>, решение для логарифмического блуждания можно записать в следующем виде:
 
 
 
:<center><math>x(t) = e^{(\mu-\sigma^2/2)t + \sigma\,W_t}.</math></center>
 
 
 
Действительно, производные для <math>\textstyle x(t)=F(t, W)</math> равны:
 
 
 
:<center><math>\frac{\partial x}{\partial t} = (\mu-\sigma^2/2) \, x,\;\;\;\;\;\;\;\;\frac{\partial x}{\partial W} = \sigma\, x,\;\;\;\;\;\;\;\;\frac{\partial^2 x}{\partial W^2} = \sigma^2\, x.</math></center>
 
 
 
Винеровское блуждание <math>\textstyle W_t</math> имеет нулевой снос <math>\textstyle a=0</math> и единичную волатильность <math>\textstyle b=1</math>. Поэтому по лемме Ито () имеем:
 
 
 
:<center><math>dx = \left( \frac{\partial x}{\partial t} + \frac{1}{2}\,\frac{\partial^2 x}{\partial W^2} \right)\, dt + \frac{\partial x}{\partial W}\,\delta W = \mu \,x\,dt + \sigma\,x\,\delta W.</math></center>
 
 
 
Роль <math>\textstyle x</math> теперь играет процесс <math>\textstyle W</math>, а функция <math>\textstyle F</math> &mdash; это <math>\textstyle x</math>.
 
 
 
Задавая различные функции <math>\textstyle x=F(t, W_t)</math>, удовлетворяющие начальному условию <math>\textstyle x_0=F(0, 0)</math>, можно найти целый класс точно решаемых стохастических уравнений. После подстановки <math>\textstyle F(t, W_t)</math> в лемму Ито необходимо исключить <math>\textstyle W_t</math>, заменив её на <math>\textstyle W_t=G(t, x)</math>, где <math>\textstyle G</math> &mdash; обратная к <math>\textstyle F</math> функция. Кроме этого константа <math>\textstyle x_0</math> должна сократиться, так как это "внешнее" к динамике условие и "порядочное уравнение" не должно зависеть от него. В качестве упражнения стоит проверить решения (<math>\textstyle \mathbf{R}_{}</math>) &mdash; (<math>\textstyle \mathbf{R}_{}</math>) из Справочника (стр. \pageref{r_exact_from_W_1}). К сожалению, чаще таким методом получаются уравнения, в которых снос зависит от волатильности шума, что не очень естественно для практических приложений.
 
 
 
<math>\textstyle \bullet</math> ''Процесс Орнштейна - Уленбека'':
 
 
 
:<center><math> { \;dx = -\beta\cdot (x-\alpha)\,dt + \sigma\,\delta W } </math></center>
 
 
 
описывает блуждание, в котором <math>\textstyle x</math> притягивается к уровню, определяемому константой <math>\textstyle \alpha</math>. При этом волатильность <math>\textstyle \sigma</math> считается постоянной. Если <math>\textstyle x\gg\alpha</math>, то снос становится заметно отрицательным и тянет процесс вниз. При опускании <math>\textstyle x</math> ниже <math>\textstyle \alpha</math> снос оказывается положительным и в среднем поднимает <math>\textstyle x(t)</math> вверх. Параметр <math>\textstyle \beta>0</math> характеризует величину "силы притяжения" к равновесному значению <math>\textstyle \alpha</math>.
 
 
 
Условие совместности () даёт уравнение <math>\textstyle \dot{s}(t)=\beta s(t)</math>. Решая его и первое уравнение () для <math>\textstyle F(x,t)</math>, мы каждый раз выбираем константы интегрирования наиболее "удобным способом", так как начальные условия уже учтены в (), а нам необходимо найти ''простейшую'' замену, исключающую <math>\textstyle x</math> из сноса и волатильности:
 
 
 
:<center><math>s(t)=\sigma e^{\beta t},\;\;\;\;\;\;\;\;\;F(x,t)=x e^{\beta t},\;\;\;\;\;\;\;\;\;\;f(t)=\alpha\beta e^{\beta t}.</math></center>
 
 
 
В результате решение записывается в следующем виде (<math>\textstyle t_0=0</math>):
 
 
 
:<center><math> x(t) = \alpha + \bigl(x_0-\alpha\bigr) e^{-\beta t} + \frac{\sigma}{\sqrt{2\beta}}\, \sqrt{1-e^{-2\beta t}}\cdot \; \varepsilon. </math></center>
 
 
 
Несложно увидеть, что <math>\textstyle x(t)</math> оказывается гауссово распределённой величиной со средним и дисперсией, зависящими от времени.
 
 
 
Если <math>\textstyle \beta>0</math>, то среднее при больших временах стремится к равновесному уровню <math>\textstyle \alpha</math>. Волатильность становится равной <math>\textstyle \sigma/\sqrt{2\beta}</math>. При винеровском или логарифмическом блуждании <math>\textstyle x(t)</math> может уйти как угодно далеко от своего начального значения <math>\textstyle x_0</math>. Для процесса () <math>\textstyle x(t)</math> "заперта" в статистическом коридоре с шириной, равной двойной волатильности <math>\textstyle \sigma/\sqrt{2\beta}</math>.
 
 
 
При малых <math>\textstyle \beta</math> процесс Орнштейна-Уленбека по своему поведению становится очень близким к обычному винеровскому блужданию. Траектория <math>\textstyle x(t)</math> достаточно долго блуждает выше или ниже <math>\textstyle \alpha</math>, не уходя, тем не менее, на бесконечность. Волатильность стремится к <math>\textstyle \sigma/\sqrt{2\beta}</math>, и тем больше, чем меньше <math>\textstyle \beta</math>. Следовательно, характерный коридор, в котором происходит блуждание, при малых <math>\textstyle \beta</math> расширяется. Если и <math>\textstyle \sigma</math>, и <math>\textstyle \beta</math> достаточно большие, <math>\textstyle x(t)</math> часто пересекает равновесный уровень, начиная напоминать обычный белый шум.
 
 
 
Наличие равновесного уровня в модели Орнштейна-Уленбека полезно для различных финансовых приложений. Например, в случае курсов валют <math>\textstyle \alpha</math> может быть паритетом покупательной способности (<math>\textstyle \lessdot</math> C), а для процентной ставки - её долгосрочным значением.
 
 
 
Примеры реализаций блуждания Орнштейна-Уленбека при различных параметрах приведены ниже. На левом рисунке <math>\textstyle \beta=0.1</math>, <math>\textstyle \sigma=0.1</math>. На правом &mdash; <math>\textstyle \beta=1</math>, <math>\textstyle \sigma=0.5</math>. Величина <math>\textstyle \alpha</math> в обоих случаях равна единице.
 
 
 
[[File:ol.png]]
 
 
 
\\
 
 
 
Необходимо помнить, что, если решение выражено через винеровскую переменную <math>\textstyle W_t</math>, её всегда можно переписать через гауссову случайную величину, заменив <math>\textstyle W_t=\varepsilon\sqrt{t}</math>. Обратное, вообще говоря, не верно. Если в решении есть <math>\textstyle \varepsilon</math>, нельзя его выразить через <math>\textstyle W_t</math>, подставив <math>\textstyle \varepsilon\to W_t/\sqrt{t}</math>. В качестве упражнения имеет смысл проверить, что подобная замена в (<math>\textstyle </math>) приводит к случайной функции, не удовлетворяющей ().
 
 
 
Можно объединить положительность <math>\textstyle x</math> и его притяжение к равновесному уровню в следующей логарифмической модели с притяжением:
 
 
 
:<center><math> dx = -\beta \cdot x\cdot\left(\ln \frac{x}{\alpha} - 1\right)\, dt + \sigma \,x \,\delta W. </math></center>
 
 
 
Если <math>\textstyle x>\alpha</math>, то снос отрицательный, а при <math>\textstyle x<\alpha</math> &mdash; положительный. Множитель <math>\textstyle x</math> "замораживает" динамику при приближении к <math>\textstyle x=0</math>. Для этой модели несложно найти точное решение (<math>\textstyle \lessdot</math> H).
 
 
 
На самом деле логарифмическая модель с притяжением является простой деформацией процесса Орнштейна-Уленбека. Действительно, если <math>\textstyle x</math> удовлетворяет уравнению (), то несложно проверить, что <math>\textstyle y=e^x</math> будет удовлетворять (). Уравнение () так же соотносится с (), как логарифмическое блуждание с процессом Винера.
 
 
 
Ещё одну модель уместно назвать ''броуновской ловушкой'':
 
 
 
:<center><math> dx = -\beta \cdot (x-\alpha)\,dt \;+\; \sigma \cdot(x-\alpha) \, \delta W. </math></center>
 
 
 
Член со сносом обеспечивает притяжение к уровню <math>\textstyle x=\alpha</math>, в окрестности которого волатильность становится очень маленькой, а динамика &mdash; детерминированной. В результате процесс рано или поздно гарантированно притягивается к значению <math>\textstyle x=\alpha</math> (<math>\textstyle \lessdot</math> H).
 
 
 
<math>\textstyle \bullet</math> Можно рассмотреть общее ''стационарное уравнение'', снос и волатильность которого не зависят от времени:
 
 
 
:<center><math>dx = a(x)\,dt + b(x)\,\delta W.</math></center>
 
 
 
Условие совместности записывается следующим образом:
 
 
 
:<center><math> \frac{\dot{s}(t)}{s(t)} = \frac{1}{2} \,b\cdot b'' - b \cdot \left(\frac{a}{b}\right)' = \gamma, </math></center>
 
 
 
где штрих производная по <math>\textstyle x</math>, точка &mdash; по времени, и опущены аргументы у функций. Левая часть зависит только от времени, правая &mdash; только от <math>\textstyle x</math>, поэтому это выражение равно некоторой константе, которую мы обозначили через <math>\textstyle \gamma</math>. Проинтегрировав это уравнение, найдём связь между сносом и волатильностью:
 
 
 
:<center><math>a= \frac{\bigl(b^2\bigr)'}{4} + \eta \cdot b - \gamma b\cdot \int \frac{dx}{b},</math></center>
 
 
 
где <math>\textstyle \eta</math> &mdash; ещё один параметр.
 
 
 
Если <math>\textstyle b(x)=\sigma=const</math> &mdash; мы приходим к уравнению Орнштейна-Уленбека (), стр. \pageref{equat_OU}. Для <math>\textstyle b(x)=\sigma x</math> точно решаемой задачей является логарифмическая модель с притяжением (), частным случаем которой является логарифмическое блуждание. При <math>\textstyle b(x)=\sigma \sqrt{x}</math> снос должен явным образом зависеть от <math>\textstyle \sigma</math>:
 
 
 
:<center><math>a(x) = \frac{\sigma^2}{4} + \alpha \sqrt{x} + 2\beta x.</math></center>
 
 
 
Решение такого уравнения имеет вид (<math>\textstyle x_0=x(0)</math>, <math>\textstyle \beta>0</math>):
 
 
 
:<center><math>x(t) =\left[ \sqrt{x_0}\,e^{\beta t} + \frac{\alpha}{2\beta}\, \left(e^{\beta t} - 1\right) + \frac{\sigma}{\sqrt{8\beta}} \,\sqrt{e^{2\beta t}-1}\cdot \varepsilon \right]^2.</math></center>
 
 
 
Если <math>\textstyle a(x)/b(x)=const</math>, или сноса при блуждании нет <math>\textstyle a(x)=0</math>, то условие совместности () упрощается:
 
 
 
:<center><math>\frac{b''}{2} = \frac{\gamma}{b}.</math></center>
 
 
 
Умножая его на интегрирующий множитель <math>\textstyle b'</math>, получаем решение в неявной форме:
 
 
 
:<center><math>x-\alpha = \int \frac{db}{\sqrt{\beta + 4 \gamma \ln b}},</math></center>
 
 
 
где <math>\textstyle \alpha</math> и <math>\textstyle \beta</math> &mdash; константы интегрирования.
 
 
 
<math>\textstyle \bullet</math> ''Броуновский мост''. Рассмотрим теперь уравнение Ито со сносом, зависящим не только от <math>\textstyle x</math>, но и от времени <math>\textstyle t</math>:
 
 
 
:<center><math> { \;dx = -\frac{x-\alpha}{T-t} \, dt + \sigma\,\delta W\; }. </math></center>
 
 
 
Константа <math>\textstyle T</math> &mdash; это выделенное время в будущем (<math>\textstyle t<T</math>), когда снос становится бесконечным. Условие совместности даёт:
 
 
 
:<center><math> s(t)=\frac{\sigma}{T-t},\;\;\;\;\;\;F(x,t)= \frac{x}{T-t},\;\;\;\;\;\;\;f(t)=\frac{\alpha}{(T-t)^2}. </math></center>
 
 
 
В результате получаем решение в следующем виде (<math>\textstyle x_0=x(t_0)</math>):
 
  
:<center><math>x(t) = \alpha + (x_0-\alpha)\, \frac {T-t}{T-t_0} + \sigma\, \sqrt{\frac{(t-t_0)(T-t)}{T-t_0}}\cdot \varepsilon.</math></center>
+
заменой иногда можно свести к частному случаю (2.17), для которого решение уже известно. Воспользуемся леммой Ито:
  
Среднее процесса при <math>\textstyle t\to T</math> стремится к <math>\textstyle \alpha</math>. При этом волатильность оказывается равной нулю. Это означает, что <math>\textstyle x(t)</math> ''гарантированно'' в процессе блуждания достигает равновесного значения <math>\textstyle x(T)=\alpha</math>:
+
{| width="100%"
 +
| width="90%" align="center"|<math> dF = \underbrace{\left(\frac{\partial F}{\partial t} \,+\, a(x,t) \frac{\partial F}{\partial x} \,+\, \frac{b^2(x,t)}{2}\, \frac{\partial^2 F}{\partial x^2}\right)}_{f(t)}\, dt + \underbrace{b(x,t)\frac{\partial F}{\partial x}}_{s(t)}\, \delta W. </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(2.20)'''</div>
 +
|}
  
[[File:bridge.png]]
+
Подберем <math>\textstyle F(x,t)</math> таким образом, чтобы множители при <math>\textstyle \delta W</math> и <math>\textstyle dt</math> в (2.20) оказались функциями <math>\textstyle s(t)</math> и <math>\textstyle f(t)</math>, зависящими только от времени:
  
\\ На рисунках в обоих случаях <math>\textstyle \alpha=1</math>. Слева <math>\textstyle \sigma=0.1</math>, справа <math>\textstyle \sigma=0.05</math>. Соединение начального условия <math>\textstyle x_0=x(0)</math> и "конечного" <math>\textstyle x(T)=\alpha</math> стохастическими траекториями и дало живописное название процессу.
+
{| width="100%"
 +
| width="90%" align="center"|<math> \frac{\partial F}{\partial x} = \frac{s(t)}{b(x,t)},\;\;\;\;\;\;\;\;\;\;\;\;\; \frac{\partial F}{\partial t} + s(t)\cdot \left[\frac{a(x,t)}{b(x,t)}-\frac{1}{2}\frac{\partial b(x,t)}{\partial x} \right] = f(t), </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(2.21)'''</div>
 +
|}
  
Можно рассмотреть броуновский мост в более общем случае, с произвольными коэффициентами, зависящими от времени:
+
где вместо <math>\textstyle \partial F/\partial x</math> в множитель при <math>\textstyle dt</math> подставлено первое уравнение (2.21) и его производная по <math>\textstyle x</math> (<math>\textstyle \lessdot</math> H). Возьмём частные производные первого уравнения (2.21) по <math>\textstyle t</math> и второго по <math>\textstyle x</math>. Вычитая их, мы придём к ''условию совместности'':
  
:<center><math>dx = -\beta(t)\cdot \Bigl(x-\alpha(t)\Bigr) \, dt + \sigma(t)\,\delta W.</math></center>
+
{| width="100%"
 +
| width="90%" align="center"|<math> \frac{1}{s(t)}\cdot \frac{\partial }{\partial t}\left\{\frac{s(t)}{b(x,t)} \right\} = \frac{1}{2}\frac{\partial^2 b(x,t)}{\partial x^2} - \frac{\partial }{\partial x}\left\{\frac{a(x,t)}{b(x,t)} \right\}. </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(2.22)'''</div>
 +
|}
  
Условия совместности дают:
+
''Если'' при данных <math>\textstyle a(x,t)</math> и <math>\textstyle b(x,t)</math> можно подобрать такую функцию <math>\textstyle s(t)</math>, при которой уравнение (2.22) обратится в тождество, то мы получим решение стохастического уравнения (2.19) в следующей неявной форме:
  
:<center><math>s(t)=\sigma(t)e^{\,\int\limits^t_{t_0} \beta(t)dt},\;\;\;F(x,t)=x\, \frac{s(t)}{\sigma(t)},\;\;\;f(t)=\alpha(t) \beta(t)\,\frac{s(t)}{\sigma(t)}.</math></center>
+
{| width="100%"
 +
| width="90%" align="center"|<math> F\bigl(x(t), t\bigr) = F\bigl(x(t_0), t_0\bigr) + \int\limits^t_{t_0} f(\tau) \,d\tau + \left[ \int\limits^t_{t_0} s^2(\tau) \,d\tau \right]^{1/2} \cdot \;\varepsilon, </math>
 +
| <div width="10%" align="right" style="color:#0000CC">'''(2.23)'''</div>
 +
|}
  
Для частного выбора <math>\textstyle \beta(t)=\beta/(T-t)</math>, <math>\textstyle \alpha(t)=\alpha</math>, <math>\textstyle \sigma(t)=\sigma</math>, где <math>\textstyle \alpha, \beta</math>, <math>\textstyle T</math> и <math>\textstyle \sigma</math> &mdash; константы модели, получаем решение в следующем виде (<math>\textstyle t_0=0</math>):
+
где функция <math>\textstyle f(t)</math> определяется вторым соотношением (2.21), а <math>\textstyle F(x,t)</math> находится из первого уравнения (2.21) (<math>\textstyle \lessdot</math> C).
  
:<center><math>x(t) = \alpha +\frac{x_0-\alpha}{T^\beta} (T-t)^\beta + \sigma\cdot \left[\frac{(T-t)}{2\beta -1} \left(1- \frac{(T-t)^{2\beta-1}}{T^{2\beta-1}} \right) \right]^{1/2} \cdot\;\varepsilon.</math></center>
+
Решение (2.23) &mdash; это нестационарный гауссовый процесс для деформации <math>\textstyle x(t)</math> при помощи нелинейной функции <math>\textstyle F(x,t)</math>. Естественно, что разрешимость (2.22) позволяет интегрировать уравнение Ито только в ряде частных случаев. Однако, как мы увидим ниже, эти случаи охватывают достаточно широкий класс важных для практических приложений процессов.
  
Заданием функции <math>\textstyle \alpha(t)</math> можно добиться произвольного выгиба "моста" вверх или вниз.
 
  
 
----
 
----

Текущая версия на 18:32, 9 марта 2010

Лемма Ито << Оглавление >> Простые стохастические модели

Несмотря на простой вид, стохастические уравнения (2.4) аналитически интегрировать не так и просто из-за винеровского члена . Это явно видно в случае конечной численной реализации (2.5). Каждое последовательное в итерационной процедуре нелинейным образом зависит от всех предыдущих случайных чисел ( C). Тем не менее, рассмотрим ситуации, в которых можно получить точные решения.

Пусть в процессе Ито функции сноса и волатильности зависят только от времени. Обозначим их через и :

(2.17)

Это уравнение легко интегрируется при помощи "дискретной" интерпретации стохастического члена . Рассмотрим итерации, выполняемые по разностной схеме (2.5):

где и . После итераций итоговое значение будет равно:

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

(2.18)

Решение (2.18) уравнения (2.17) говорит нам, что является нормально распределённым случайным числом со средним и дисперсией, зависящими от времени. Если — не константа, то будущая неопределённость в значении может увеличиваться уже не как , а по другому закону.

Соотношение (2.18) позволяет легко вычислить статистические свойства процесса, в частности, его среднее и волатильность .

Одномерное уравнение Ито для процесса, имеющего произвольный снос и волатильность

(2.19)

заменой иногда можно свести к частному случаю (2.17), для которого решение уже известно. Воспользуемся леммой Ито:

(2.20)

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

(2.21)

где вместо в множитель при подставлено первое уравнение (2.21) и его производная по ( H). Возьмём частные производные первого уравнения (2.21) по и второго по . Вычитая их, мы придём к условию совместности:

(2.22)

Если при данных и можно подобрать такую функцию , при которой уравнение (2.22) обратится в тождество, то мы получим решение стохастического уравнения (2.19) в следующей неявной форме:

(2.23)

где функция определяется вторым соотношением (2.21), а находится из первого уравнения (2.21) ( C).

Решение (2.23) — это нестационарный гауссовый процесс для деформации при помощи нелинейной функции . Естественно, что разрешимость (2.22) позволяет интегрировать уравнение Ито только в ряде частных случаев. Однако, как мы увидим ниже, эти случаи охватывают достаточно широкий класс важных для практических приложений процессов.



Лемма Ито << Оглавление >> Простые стохастические модели

Стохастический мир - простое введение в стохастические дифференциальные уравнения