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

Материал из synset
Перейти к: навигация, поиск
(Новая страница: «{| width="100%" | width="40%"|Лемма Ито << ! width="20%"|Оглавление | width="40%" align="right"|…»)
 
 
(не показаны 4 промежуточные версии этого же участника)
Строка 5: Строка 5:
 
|}
 
|}
 
----
 
----
 +
 +
<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>:
 +
 +
{| 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>. Рассмотрим итерации, выполняемые по [[Уравнения Ито|разностной схеме (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>
 +
 +
где <math>\textstyle f_k=f(t_k)</math> и <math>\textstyle s_k=s(t_k)</math>. После <math>\textstyle n</math> итераций итоговое значение будет равно:
 +
 +
:<center><math>x = x_0 + (f_0+...+f_{n-1})\,\Delta t + (s_0\varepsilon_1 + ... + s_{n-1}\varepsilon_{n}) \sqrt{\Delta t}.</math></center>
 +
 +
Скобки в последнем слагаемом содержат сумму независимых гауссовых чисел, каждое из которых имеет волатильность <math>\textstyle s_k</math>. В результате получается гауссово число с волатильностью <math>\textstyle \sqrt{s^2_0+...+s^2_{n-1}}</math>. Поэтому, переходя к непрерывному пределу, получаем (<math>\textstyle \lessdot</math> H):
 +
 +
{| 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>
 +
|}
 +
 +
Решение (2.18) уравнения (2.17) говорит нам, что <math>\textstyle x(t)</math> является нормально распределённым случайным числом со средним и дисперсией, зависящими от времени. Если <math>\textstyle s(t)</math> &mdash; не константа, то будущая неопределённость в значении <math>\textstyle x</math> может увеличиваться уже не как <math>\textstyle \sqrt{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>
 +
 +
{| 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>
 +
|}
 +
 +
заменой иногда можно свести к частному случаю (2.17), для которого решение уже известно. Воспользуемся леммой Ито:
 +
 +
{| 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>
 +
|}
 +
 +
Подберем <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>, зависящими только от времени:
 +
 +
{| 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>. Вычитая их, мы придём к ''условию совместности'':
 +
 +
{| 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) в следующей неявной форме:
 +
 +
{| 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 f(t)</math> определяется вторым соотношением (2.21), а <math>\textstyle F(x,t)</math> находится из первого уравнения (2.21) (<math>\textstyle \lessdot</math> C).
 +
 +
Решение (2.23) &mdash; это нестационарный гауссовый процесс для деформации <math>\textstyle x(t)</math> при помощи нелинейной функции <math>\textstyle F(x,t)</math>. Естественно, что разрешимость (2.22) позволяет интегрировать уравнение Ито только в ряде частных случаев. Однако, как мы увидим ниже, эти случаи охватывают достаточно широкий класс важных для практических приложений процессов.
 +
  
 
----
 
----

Текущая версия на 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) позволяет интегрировать уравнение Ито только в ряде частных случаев. Однако, как мы увидим ниже, эти случаи охватывают достаточно широкий класс важных для практических приложений процессов.



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

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