Лемма Ито — различия между версиями

Материал из synset
Перейти к: навигация, поиск
 
(не показана 1 промежуточная версия этого же участника)
Строка 10: Строка 10:
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> dF = A(x,t)\;dt + B(x,t)\;\delta W </math>
 
  | width="90%" align="center"|<math> dF = A(x,t)\;dt + B(x,t)\;\delta W </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(2.13)'''</div>
 
  |}
 
  |}
  
Строка 19: Строка 19:
 
:<center><math>F(x,t) = F(x_0,t_0) + \frac{\partial F}{\partial x_0}\; \Delta x + \frac{1}{2}\frac{\partial^2 F}{\partial x_0^2}\; (\Delta x)^2 +...+\frac{\partial F}{\partial t_0}\; \Delta t + ...,</math></center>
 
:<center><math>F(x,t) = F(x_0,t_0) + \frac{\partial F}{\partial x_0}\; \Delta x + \frac{1}{2}\frac{\partial^2 F}{\partial x_0^2}\; (\Delta x)^2 +...+\frac{\partial F}{\partial t_0}\; \Delta t + ...,</math></center>
  
где все производные справа вычислены в точке <math>\textstyle x_0, t_0</math>. Для ряда оставлен член второго порядка малости по <math>\textstyle \Delta x</math>. При помощи () мы можем записать <math>\textstyle (\Delta x)^2</math> в следующем виде:
+
где все производные справа вычислены в точке <math>\textstyle x_0, t_0</math>. Для ряда оставлен член второго порядка малости по <math>\textstyle \Delta x</math>. При помощи (2.7) мы можем записать <math>\textstyle (\Delta x)^2</math> в следующем виде:
  
 
:<center><math>(\Delta x)^2 = \bigl( a_0\;\Delta t + b_0\; \varepsilon \sqrt{\Delta t}+ ...\bigr)^2 = b_0^2 \;\varepsilon^2 \;\Delta t +...,</math></center>
 
:<center><math>(\Delta x)^2 = \bigl( a_0\;\Delta t + b_0\; \varepsilon \sqrt{\Delta t}+ ...\bigr)^2 = b_0^2 \;\varepsilon^2 \;\Delta t +...,</math></center>
Строка 27: Строка 27:
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> F=F_0 + \frac{\partial F}{\partial x_0} \cdot (a_0 \Delta t + b_0 \;\varepsilon \sqrt{\Delta t}) + \frac{b^2_0}{2}\;\frac{\partial^2 F}{\partial x_0^2} \;\varepsilon^2 \;\Delta t + \frac{\partial F}{\partial t_0} \;\Delta t+... </math>
 
  | width="90%" align="center"|<math> F=F_0 + \frac{\partial F}{\partial x_0} \cdot (a_0 \Delta t + b_0 \;\varepsilon \sqrt{\Delta t}) + \frac{b^2_0}{2}\;\frac{\partial^2 F}{\partial x_0^2} \;\varepsilon^2 \;\Delta t + \frac{\partial F}{\partial t_0} \;\Delta t+... </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(2.14)'''</div>
 
  |}
 
  |}
  
По ''определению'' () коэффициент сноса в пределе <math>\textstyle \Delta t\to 0</math> равен:
+
По ''определению'' [[Уравнения Ито|(2.6)]] коэффициент сноса в пределе <math>\textstyle \Delta t\to 0</math> равен:
  
 
:<center><math>A(x_0,t_0) \;=\; \frac{\left\langle F-F_0\right\rangle }{\Delta t} \;=\; a_0 \cdot \frac{\partial F}{\partial x_0} + \frac{b^2_0}{2}\cdot \frac{\partial^2 F}{\partial x_0^2} + \frac{\partial F}{\partial t_0},</math></center>
 
:<center><math>A(x_0,t_0) \;=\; \frac{\left\langle F-F_0\right\rangle }{\Delta t} \;=\; a_0 \cdot \frac{\partial F}{\partial x_0} + \frac{b^2_0}{2}\cdot \frac{\partial^2 F}{\partial x_0^2} + \frac{\partial F}{\partial t_0},</math></center>
  
где подставлено разложение () для <math>\textstyle F</math> и учтено, что <math>\textstyle \left\langle \varepsilon\right\rangle =0</math>, <math>\textstyle \left\langle \varepsilon^2\right\rangle =1</math>. Аналогично, для коэффициента диффузии:
+
где подставлено разложение (2.14) для <math>\textstyle F</math> и учтено, что <math>\textstyle \left\langle \varepsilon\right\rangle =0</math>, <math>\textstyle \left\langle \varepsilon^2\right\rangle =1</math>. Аналогично, для коэффициента диффузии:
  
 
:<center><math>B^2(x_0,t_0) \;=\; \frac{\left\langle (F-F_0)^2\right\rangle }{\Delta t} \;=\; b^2_0 \cdot \left(\frac{\partial F}{\partial x_0}\right)^2.</math></center>
 
:<center><math>B^2(x_0,t_0) \;=\; \frac{\left\langle (F-F_0)^2\right\rangle }{\Delta t} \;=\; b^2_0 \cdot \left(\frac{\partial F}{\partial x_0}\right)^2.</math></center>
Строка 40: Строка 40:
 
Для моментов более высоких порядков в пределе <math>\textstyle \Delta t\to 0</math> получается ноль. Таким образом, это ''действительно'' диффузный процесс.
 
Для моментов более высоких порядков в пределе <math>\textstyle \Delta t\to 0</math> получается ноль. Таким образом, это ''действительно'' диффузный процесс.
  
Выше упоминалось, что для записи стохастического дифференциального уравнения некоторого процесса необходимо вычислить условные средние его изменения первого и второго порядка. При этом следует убедиться, что моменты более высоких порядков при <math>\textstyle \Delta t\to 0</math> стремятся к нулю. Если этого не происходит, то процесс не является диффузным и не может быть записан в форме Ито (). Поэтому необходима полная проверка "диффузности", проведенная выше.
+
Выше упоминалось, что для записи стохастического дифференциального уравнения некоторого процесса необходимо вычислить условные средние его изменения первого и второго порядка. При этом следует убедиться, что моменты более высоких порядков при <math>\textstyle \Delta t\to 0</math> стремятся к нулю. Если этого не происходит, то процесс не является диффузным и не может быть записан в форме Ито (2.13). Поэтому необходима полная проверка "диффузности", проведенная выше.
  
Считая уравнение () первым вычислением в бесконечной итерационной схеме, мы можем также рассуждать аналогично разделу <math>\textstyle \S</math>. Сумма слагаемых вида <math>\textstyle \varepsilon^2 \Delta t</math> приводит к такому же детерминированному результату, как и в отсутствие <math>\textstyle \varepsilon^2</math>. Поэтому можно положить <math>\textstyle \varepsilon^2\to 1</math>.
+
Считая уравнение (2.14) первым вычислением в бесконечной итерационной схеме, мы можем также рассуждать аналогично разделу <math>\textstyle \S</math>. Сумма слагаемых вида <math>\textstyle \varepsilon^2 \Delta t</math> приводит к такому же детерминированному результату, как и в отсутствие <math>\textstyle \varepsilon^2</math>. Поэтому можно положить <math>\textstyle \varepsilon^2\to 1</math>.
  
 
Так как начальный момент был выбран произвольным образом, запишем дифференциал функции <math>\textstyle F(x,t)</math> в форме Ито при помощи бесконечно малой винеровской переменной <math>\textstyle \delta W=\varepsilon\sqrt{dt}</math>:
 
Так как начальный момент был выбран произвольным образом, запишем дифференциал функции <math>\textstyle F(x,t)</math> в форме Ито при помощи бесконечно малой винеровской переменной <math>\textstyle \delta W=\varepsilon\sqrt{dt}</math>:
  
 
{| width="100%"  
 
{| width="100%"  
  | width="90%" align="center"|<math> { \;dF \;= {\Bigr|}}\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)dt + b(x,t)\,\frac{\partial F}{\partial x}\;\delta W\; }. </math>
+
  | width="90%" align="center"|
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
<math> \;dF \;= \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)dt + b(x,t)\,\frac{\partial F}{\partial x}\;\delta W.  
 +
</math>
 +
  | <div width="10%" align="right" style="color:#0000CC">'''(2.15)'''</div>
 
  |}
 
  |}
  
Строка 57: Строка 59:
 
{| width="100%"  
 
{| width="100%"  
 
  | width="90%" align="center"|<math> dF \;=\; \frac{\partial F}{\partial t}\,dt + \frac{\partial F}{\partial x}\, dx \;=\; \left(\frac{\partial F}{\partial t} + a(x,t)\, \frac{\partial F}{\partial x}\right)\, dt. </math>
 
  | width="90%" align="center"|<math> dF \;=\; \frac{\partial F}{\partial t}\,dt + \frac{\partial F}{\partial x}\, dx \;=\; \left(\frac{\partial F}{\partial t} + a(x,t)\, \frac{\partial F}{\partial x}\right)\, dt. </math>
  | <div width="10%" align="right" style="color:#0000CC">'''(EQN)'''</div>
+
  | <div width="10%" align="right" style="color:#0000CC">'''(2.16)'''</div>
 
  |}
 
  |}
  
 
В отличие от этого соотношения, в детерминированную часть леммы Ито проникают функция диффузии <math>\textstyle b^2(x,t)</math> и вторая производная по <math>\textstyle x</math>. Происходит это, как мы видели, благодаря корню <math>\textstyle \sqrt{dt}</math>. Это, в свою очередь, связано со свойствами простого аддитивного блуждания, которое является локальным приближением любого процесса Ито.
 
В отличие от этого соотношения, в детерминированную часть леммы Ито проникают функция диффузии <math>\textstyle b^2(x,t)</math> и вторая производная по <math>\textstyle x</math>. Происходит это, как мы видели, благодаря корню <math>\textstyle \sqrt{dt}</math>. Это, в свою очередь, связано со свойствами простого аддитивного блуждания, которое является локальным приближением любого процесса Ито.
  
Для винеровского уравнения <math>\textstyle dx=\mu\, dt+\sigma\delta W</math> с постоянным сносом <math>\textstyle \mu</math> и волатильностью <math>\textstyle \sigma</math> дифференциал ''квадрата'' траектории <math>\textstyle y=x^2</math>, в соответствии с (), удовлетворяет нелинейному уравнению Ито:
+
Для винеровского уравнения <math>\textstyle dx=\mu\, dt+\sigma\delta W</math> с постоянным сносом <math>\textstyle \mu</math> и волатильностью <math>\textstyle \sigma</math> дифференциал ''квадрата'' траектории <math>\textstyle y=x^2</math>, в соответствии с (2.15), удовлетворяет нелинейному уравнению Ито:
  
 
:<center><math>d(x^2) = (2\mu x + \sigma^2)\,dt + 2\sigma\,x\,\delta W\;\;\;\;\;\;=>\;\;\;\;\;\;\;dy = (2\mu \sqrt{y}+\sigma^2)\,dt + 2\sigma\sqrt{y}\,\delta W.</math></center>
 
:<center><math>d(x^2) = (2\mu x + \sigma^2)\,dt + 2\sigma\,x\,\delta W\;\;\;\;\;\;=>\;\;\;\;\;\;\;dy = (2\mu \sqrt{y}+\sigma^2)\,dt + 2\sigma\sqrt{y}\,\delta W.</math></center>

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

Почему Ито << Оглавление >> Точные решения уравнения Ито

Пусть процесс подчиняется уравнению Ито. Рассмотрим обычную гладкую функцию . Если вместо в неё подставить , то станет случайным процессом. Покажем, что он также подчиняется диффузному уравнению Ито:

(2.13)

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

Разложим в ряд Тейлора в окрестности начального фиксированного значения по небольшим и :

где все производные справа вычислены в точке . Для ряда оставлен член второго порядка малости по . При помощи (2.7) мы можем записать в следующем виде:

где оставлено ведущее приближение по . Таким образом, если в начальный момент времени функция равна детерминированному числу , то через малый промежуток времени, в зависимости от значения , это будет случайная величина вида ( C):

(2.14)

По определению (2.6) коэффициент сноса в пределе равен:

где подставлено разложение (2.14) для и учтено, что , . Аналогично, для коэффициента диффузии:

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

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

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

Так как начальный момент был выбран произвольным образом, запишем дифференциал функции в форме Ито при помощи бесконечно малой винеровской переменной :

(2.15)

Это соотношение называется "леммой Ито". Оно играет очень важную роль в теории случайных процессов ( C).

Обратим внимание, что в отсутствие стохастики полный дифференциал функции , в которую подставили решение уравнения , имеет вид:

(2.16)

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

Для винеровского уравнения с постоянным сносом и волатильностью дифференциал квадрата траектории , в соответствии с (2.15), удовлетворяет нелинейному уравнению Ито:

Действуя в обратном направлении, при помощи подходящей замены и леммы Ито можно сводить одни уравнения к другим, решение которых нам известно. Рассмотрим этот подход подробнее.


Почему Ито << Оглавление >> Точные решения уравнения Ито

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