VIII

Интеграл Ито. Стохастические дифференциальные уравнения

[27/0%]
LaTeX
Задача VIII.1

Приведите пример адаптированного, но не прогрессивно измеримого процесса.

Следующие две σ\sigma-алгебры (опииональных и предсказуемых множеств) играют в стохастическом анализе важную роль в самых разнообразных вопросах.

Определение 8. Опциональная σ\sigma-алгебра O\mathscr {O} подмножеств R+×Ω\mathbb {R}_{+} \times \Omega есть σ\sigma-алгебра, порожденная всеми согласованными процессами X={Xt,t⩾0}X = \left\{ X_{t}, t \geqslant 0\right\}, рассматриваемыми как отображения (t,ω)↦Xt(ω)(t, \omega ) \mapsto X_{t}(\omega ) из R+×Ω\mathbb {R}_{+} \times \Omega в (R,B(R))(\mathbb {R}, \mathscr {B}(\mathbb {R})), траектории которых принадлежат пространству Скорохода D[0,∞)D[0, \infty ). Процессы, которые измеримы относительно σ\sigma-алгебры O\mathscr {O}, называются опциональными.

?
Задача VIII.2

Докажите, что равносильным образом опциональная σ\sigma-алгебра O\mathscr {O} может быть определена как σ\sigma-алгебра подмножеств в R+×Ω\mathbb {R}_{+} \times \Omega, порожденная стохастическими интервалами вида

[[0,τ[[:={(t,ω):0⩽t<τ(ω)} [[0, \tau [[: = \left\{ (t, \omega ): 0 \leqslant t < \tau (\omega )\right\}

где τ(ω)\tau (\omega ) — марковские моменты (относительно потока F\mathbb {F} ).

Определение 9. Предсказуемая σ\sigma-алгебра P\mathscr {P} подмножеств R+×Ω\mathbb {R}_{+} \times \Omega есть σ\sigma-алгебра, порожденная всеми согласованными процессами X={Xt,t⩾0}X = \left\{ X_{t}, t \geqslant 0\right\}, рассматриваемыми как отображения (t,ω)↦Xt(ω)(t, \omega ) \mapsto X_{t}(\omega ) из R+×Ω\mathbb {R}_{+} \times \Omega в (R,B(R))(\mathbb {R}, \mathscr {B}(\mathbb {R})), траектории которых принадлежат пространству C=C(R+,R)C = C\left(\mathbb {R}_{+}, \mathbb {R}\right) непрерывных на R+\mathbb {R}_{+}действительных функций. Процессы, которые измеримы относительно σ\sigma-алгебры P\mathscr {P}, называются предсказуемыми.

?
Задача VIII.3

Докажите что равносильным образом предсказуемая σ\sigma-алгебра P\mathscr {P} может быть определена как σ\sigma-алгебра, порожденная всеми согласованными процессами, траектории которых лишь непрерывны слева на (0,∞)(0, \infty ).

Имеются и другие альтернативные (и во многих рассмотрениях более удобные) определения, см., например, [26] и [46].

?
Задача VIII.4

Проверьте, что σ\sigma-алгебра P\mathscr {P} порождается каждым из следующих классов подмножеств:

?
(а)

{0}×A\left\{ 0\right\} \times A с A∈F0A \in \mathscr {F}_{0} и множествами вида (t,∞)×B(t, \infty ) \times B с B∈Ft,t⩾0B \in \mathscr {F}_{t}, t \geqslant 0;

(б)

{0}×A\left\{ 0\right\} \times A с A∈F0A \in \mathscr {F}_{0} и стохастическими интервалами

[[0,τ]]={(t,ω):0⩽t⩽τ(ω)} или ]]τ,∞[[={(t,ω):τ(ω)<t<∞},  [[0, \tau ]] = \left\{ (t, \omega ): 0 \leqslant t \leqslant \tau (\omega )\right\} \text{ или }]] \tau , \infty [[ = \left\{ (t, \omega ): \tau (\omega ) < t < \infty \right\} \text{, }

где τ\tau — марковский момент;

(в)

{0}×A\left\{ 0\right\} \times A с A∈F0A \in \mathscr {F}_{0} и множествами вида (s,t]×B(s, t] \times B с B∈Fs,s⩽tB \in \mathscr {F}_{s}, s \leqslant t.

Задача VIII.5

Убедитесь, что между введенными σ\sigma-алгебрами имеют место следующие включения:

P⊆O⊆Prog⁡⊆B(R+)⊗F∞ \mathscr {P} \subseteq \mathscr {O} \subseteq \operatorname {Prog} \subseteq \mathscr {B}\left(\mathbb {R}_{+}\right) \otimes \mathscr {F} \infty

Напомним, что по теореме 4 всякий согласованный с фильтрацией F=(Ft)t⩾0\mathbb {F} = \left(\mathscr {F}_{t}\right)_{t \geqslant 0} непрерывный слева или непрерывный процесс является прогрессивно измеримым.

?
Задача VIII.6

Пусть процесс XX является прогрессивно измеримым и τ\tau — конечный марковский момент. Покажите, что случайная величина XτX_{\tau } будет Fτ\mathscr {F}_{\tau }-измеримой (это свойство является едва ли не определяющим целесообразность введения понятия прогрессивной измеримости).

?
Задача VIII.7

Проверьте, что простые (кусочно-непрерывные справа) функции f=f(t,ω)f = f(t, \omega ) вида

f(t,ω)=∑k=0m−1ξk(ω)1[tk,tk+1)(t),t∈[0,T] f(t, \omega ) = \sum _{k = 0}^{m-1} \xi _{k}(\omega ) \mathbf{1}_{\left[t_{k}, t_{k+1}\right)}(t), \quad t \in [0, T]

где величины ξk(ω)\xi_{k}(\omega ) — Ftk\mathscr {F}_{t_{k}}-измеримы, 0=t0<⋯<tm=T0 = t_{0} < \cdots < t_{m} = T, являются опциональными.

?
Задача VIII.8

Покажите, что простые (кусочно-непрерывные слева для t∈(0,T]t \in (0, T] ) функции f=f(t,ω)f = f(t, \omega ) вида

f(t,ω)=η(ω)1{0}(t)+∑k=0m−1ηk(ω)1(tk,tk+1](t),t∈[0,T] f(t, \omega ) = \eta (\omega ) \mathbf{1}_{\left\{ 0\right\} }(t)+\sum _{k = 0}^{m-1} \eta _{k}(\omega ) \mathbf{1}_{\left(t_{k}, t_{k+1}\right]}(t), \quad t \in [0, T]

где величина η(ω)\eta (\omega ) — F0\mathscr {F}_{0}-измерима, ηk(ω)\eta_{k}(\omega ) — Ftk\mathscr {F}_{t_{k}}-измеримы, 0=t0<⋯<tm=T0 = t_{0} < \cdots < t_{m} = T, являются предсказуемыми.

?
Задача VIII.9

Проведите построение интеграла Ито, изложенное в § 2 и § 3, отправляясь от класса опциональных простых функций, фигурирующих в упражнении 7 (которым приписывается то же значение интеграла Ито, что и для предсказуемых простых функций из упражнения 8). Можно ли утверждать, что замыкание в LT2L_{T}^{2} данной совокупности простых функций, входящих в LT2L_{T}^{2}, приводит к пространству AT∩LT2\mathscr {A}_{T} \cap L_{T}^{2} ?

?
Задача VIII.10

Пусть τ\tau — марковский момент такой, что τ(ω)⩽T\tau (\omega ) \leqslant T при всех ω∈Ω\omega \in \Omega (T — положительная константа). Докажите, что Iτ(f)=IT(f1(0,τ])I_{\tau }(f) = I_{T}\left(f \mathbf{1}_{(0, \tau ]}\right), где Iτ(f):=:=Iτ(ω)(f)I_{\tau }(f): = : = I_{\tau (\omega )}(f).

?
Задача VIII.11

Докажите, что если f:[0,∞)→Rf:[0, \infty ) \rightarrow \mathbb {R} есть неслучайная функция из пространства Скорохода D[0,∞)D[0, \infty ), то процесс It(f),t⩾0I_{t}(f), t \geqslant 0, является гауссовским. Найдите его среднее и ковариационную функцию.

Обсудим теперь некоторые обобщения интеграла Ито. Интегрирование можно определить для следующего класса функций J1J_{1}.

Определение 10. Будем говорить, что функция ff принадлежит классу J1J_{1}, если f:[0,∞)×Ω→Rf:[0, \infty ) \times \Omega \rightarrow \mathbb {R} прогрессивно измерима и

P(∫0tf2(s,ω)ds<∞)=1,t>0 \mathbb {P}\left(\int _{0}^{t} f^{2}(s, \omega ) d s < \infty \right) = 1, \quad t > 0

Идея построения интеграла It(f)I_{t}(f) для функции f∈J1f \in J_{1} состоит в том, чтобы подобрать последовательность функций fn,n∈Nf_{n}, n \in \mathbb {N}, для которых интеграл It(fn)I_{t}\left(f_{n}\right) уже построен, и таких, что

∫0t(f(s,ω)−fn(s,ω))2ds→P0 при n→∞ \int _{0}^{t}\left(f(s, \omega )-f_{n}(s, \omega )\right)^{2} d s \stackrel{\mathbb {P}}{\rightarrow } 0 \quad \text{ при } n \rightarrow \infty

Отсюда выводится, что последовательность {It(fn)}n∈N\left\{ I_{t}\left(f_{n}\right)\right\}_{n \in \mathbb {N}} является фундаментальной по вероятности. Следовательно, существует случайная величина, обозначаемая It(f)I_{t}(f), такая, что It(fn)→PIt(f),n→∞I_{t}\left(f_{n}\right) \stackrel{\mathbb {P}}{\rightarrow } I_{t}(f), n \rightarrow \infty. Предел It(f)I_{t}(f) записывают, по-прежнему, как ∫(0,t]f(s,ω)dWs\int_{(0, t]} f(s, \omega ) d W_{s} или (f⋅W)t(f \cdot W)_{t}. Часто броуновское движение обозначают {Bt,t⩾0}\left\{ B_{t}, t \geqslant 0\right\}, и тогда пишут (f⋅B)t(f \cdot B)_{t}.

?
Задача VIII.12

Докажите, что существует последовательность {fn,n∈N}\left\{ f_{n}, n \in \mathbb {N}\right\}, описанная выше в (85). Докажите, что если f∈J1f \in J_{1}, то существует п. н. непрерывная модификация процесса It(f),t⩾0I_{t}(f), t \geqslant 0.

Отметим, что для f∈J1f \in J_{1} процесс It(f)I_{t}(f) не обязан быть мартингалом, но является локальным мартингалом (т.е. процессом, для которого существует последовательность моментов остановки τn\tau_{n} таких, что τn↑∞\tau_{n} \uparrow \infty п.н. (n→∞)(n \rightarrow \infty ) и для каждого nn "остановленные" процессы Itτn(f):=It∧τn(f)I_{t}^{\tau_{n}}(f): = I_{t \wedge \tau_{n}}(f), где t⩾0t \geqslant 0, являются мартингалами).

?
Задача VIII.13

Рассмотрим для f∈J1f \in J_{1} процесс

Zt=exp(∫0tf(s,ω)dWs−12∫0tf2(s,ω)ds),t⩾0 Z_{t} = \mathrm{exp}\left(\int _{0}^{t} f(s, \omega ) d W_{s}-\frac{1}{2} \int _{0}^{t} f^{2}(s, \omega ) d s\right), \quad t \geqslant 0

Докажите, что dZt=Ztf(t,ω)dWtd Z_{t} = Z_{t} f(t, \omega ) d W_{t} (по-прежнему, все дифференциальные соотношения означают сокращенную запись соответствующих интегральных).

Подробнее со стохастическим исчислением можно ознакомиться по книгам [26[26, 124,128]124,128]. Кратные стохастические интегралы рассматриваются в [163].

Стохастический интеграл по винеровскому процессу (броуновскому движению) позволяет дать "явное" представление фрактального броуновского движения.

Определение 11. Случайньй процесс X={Xt,t⩾0}X = \left\{ X_{t}, t \geqslant 0\right\} со значениями в Rm\mathbb {R}^{m} называется автомодельным (самоподобным), если для каждого a>0a > 0 можно найти такое b>0b > 0, что

Law⁡(Xat,t⩾0)=Law⁡(bXt,t⩾0) \operatorname {Law}\left(X_{a t}, t \geqslant 0\right) = \operatorname {Law}\left(b X_{t}, t \geqslant 0\right)

Иначе говоря, для таких процессов изменение временно́й шкалы ( t↦at)t \mapsto a t) приводит к тому же самому результату, что и изменение фазовой шкаль ( x↦bx)x \mapsto b x). Если в определении (87) для любого a>0a > 0 параметр b=aHb = a^{H}, то случайный процесс XX называется автомодельным процессом с показателем Харста Н. Величина D=1/HD = 1 / H называется статистической фрактальной размерностью случайного процесса XX.

Напомним, что фрактальное броуновское движение B(H)={B(H)(t),t⩾0}B^{(H)} = \left\{ B^{(H)}(t), t \geqslant 0\right\} для 0<H⩽10 < H \leqslant 1 определяется как центрированный гауссовский процесс с ковариационной функцией

cov⁡(B(H)(s),B(H)(t))=s2H+t2H−∣s−t∣2H,s,t⩾0. \operatorname {cov}\left(B^{(H)}(s), B^{(H)}(t)\right) = s^{2 H}+t^{2 H}-\left|s-t\right|^{2 H}, \quad s, t \geqslant 0.

При H=1/2H = 1 / 2 процесс B(1/2)={B(1/2)(t),t⩾0}B^{(1 / 2)} = \left\{ B^{(1 / 2)}(t), t \geqslant 0\right\} является броуновским движением.

?
Задача VIII.14

Докажите, что B(H)={B(H)(t),t⩾0}B^{(H)} = \left\{ B^{(H)}(t), t \geqslant 0\right\}, где 0<H⩽10 < H \leqslant 1, есть автомодельный процесс (с показателем Харста HH ).

Процессы B(H)B^{(H)} впервые рассматривались А. Н. Колмогоровым в 1940 г. (см. [33]) и были названы им спиралями Винера. Термин фрактальное (или дробное, "fractional") броуновское движение был введен Б. Мандельброттом и Й. Ван Нессом в 1968 г. в работе [164], и там же была получена следующая теорема.

Теорема 10 (Мандельбротт, Ван Несс). Для 0<H<10 < H < 1 и t⩾0t \geqslant 0 фрактальное броуновское движение B(H)B^{(H)} допускает представление

B(H)(t)=cH{∫−∞0[(t−s)H−1/2−(−s)H−1/2]dWs+∫0t(t−s)H−1/2dWs}, B^{(H)}(t) = c_{H}\left\{ \int _{-\infty }^{0}\left[(t-s)^{H-1 / 2}-(-s)^{H-1 / 2}\right] d W_{s}+\int _{0}^{t}(t-s)^{H-1 / 2} d W_{s}\right\} ,

где нормирующая константа cHc_{H} выбрана так, что E[B(H)(1)]2=1\mathbb {E}\left[B^{(H)}(1)\right]^{2} = 1 (см., например, [86; т. 1, с. 281]), a {Ws,s⩾0}u{W~s,s⩾0}\left\{ W_{s}, s \geqslant 0\right\} u\left\{ \widetilde{W}_{s}, s \geqslant 0\right\} — независимые стандартные винеровские процессы, причем Ws=W~−sW_{s} = \widetilde{W}_{-s} для s⩽0s \leqslant 0.

Φ\Phi рактальное броуновское движение (и его дискретные аналоги) используются во многих важных моделях, в частности, для описания динамики финансовых индексов (см. [86]). Сложность изучения процессов B(H)(t)B^{(H)}(t) с 0<H⩽10 < H \leqslant 1 состоит в том, что за исключением случая H=1/2H = 1 / 2 (броуновского движения) и H=1H = 1 эти процессы не являются семимартингалами (т.е. не входят в важный класс процессов, для которых развито стохастическое исчисление, см. [46; гл. 4]). Если в представлении (88) вместо HH взять гёльдеровскую функцию HtH_{t} (т. е. ∣Ht−Hs∣⩽c∣t−s∣α,α>0\left|H_{t}-H_{s}\right| \leqslant c\left|t-s\right|^{\alpha }, \alpha > 0 ) со значениями в (0,1)(0,1), то получается случайный процесс, называемый мультифрактальным броуновским движением.

?
Задача VIII.15

Докажите, что для процесса B(H)B^{(H)} при n→∞n \rightarrow \infty

H^n:=ln⁡(n−1∑k=1n∣B(H)(k/n)−B(H)((k−1)/n)∣)/ln⁡(1/n)→H п.н.  \widehat{H}_{n}: = \ln \left(n^{-1} \sum _{k = 1}^{n}\left|B^{(H)}(k / n)-B^{(H)}((k-1) / n)\right|\right) / \ln (1 / n) \rightarrow H \text{ п.н. }

Пусть J1([0,T])J_{1}([0, T]) — пространство прогрессивно измеримых функций f:[0,T]×Ω→Rf:[0, T] \times \Omega \rightarrow \mathbb {R}, удовлетворяющих (84) при t∈[0,T]t \in [0, T].

Следующая теорема, основанная на понятии стохастического интеграла, описывает структуру броуновских функционалов.

Теорема 11 (Ито, Кларк). Пусть (FtW)t∈[0,T]\left(\mathscr {F}_{t}^{W}\right)_{t \in [0, T]} — естественная фильтрация броуновского движения W={Wt,t∈[0,T]}W = \left\{ W_{t}, t \in [0, T]\right\} и X=X(ω)X = X(\omega ) является FTW∣B(R)\mathscr {F}_{T}^{W} \mid \mathscr {B}(\mathbb {R}) измеримой случайной величиной. Справедливы следующие утверждения.

  1. Если E[X2]<∞\mathbb {E}\left[X^{2}\right] < \infty, то найдется такой стохастический процесс f={f(s),s∈[0,T]}∈AT∩LT2f = \left\{ f(s), s \in [0, T]\right\} \in \mathscr {A}_{T} \cap L_{T}^{2}, что

    X=E[X]+∫0Tf(s,ω)dWs п.н.  X = \mathbb {E}\left[X\right]+\int _{0}^{T} f(s, \omega ) d W_{s} \quad \text{ п.н. }
  2. Если E[∣X∣]<∞\mathbb {E}\left[\left|X\right|\right] < \infty, то представление (89) справедливо с некоторым процессом f∈J1([0,T])f \in J_{1}([0, T]).

  3. Если XX — положительная случайная величина (т.е. P(X>0)=1)\mathbb {P}\left(X > 0\right) = 1) и E[X]<∞\mathbb {E}\left[X\right] < \infty, то найдется процесс f∈J1[0,T]f \in J_{1}[0, T] такой, что XX допускает представление вида X=ZTE[X]X = Z_{T} \mathbb {E}\left[X\right], где ZTZ_{T} определяется согласно (86).

?
Задача VIII.16

Объясните, почему величина XX, фигурирующая в теореме 11 , есть функционал от броуновского движения, т.е. X(ω)=g(W(s,ω),0⩽s⩽T)X(\omega ) = g(W(s, \omega ), 0 \leqslant s \leqslant T), где g:C[0,T]→Rg: C[0, T] \rightarrow \mathbb {R} и g∈B(C[0,T])∣B(R)g \in \mathscr {B}(C[0, T]) \mid \mathscr {B}(\mathbb {R}).

Теорема 11 (Ито, Кларк). Пусть (FtW)t∈[0,T]\left(\mathscr {F}_{t}^{W}\right)_{t \in [0, T]} — естественная фильтрация броуновского движения W={Wt,t∈[0,T]}W = \left\{ W_{t}, t \in [0, T]\right\} и X=X(ω)X = X(\omega ) является FTW∣B(R)\mathscr {F}_{T}^{W} \mid \mathscr {B}(\mathbb {R})-измеримой случайной величиной. Справедливы следующие утверждения.

  1. Если E[X2]<∞\mathbb {E}\left[X^{2}\right] < \infty, то найдется такой стохастический процесс f={f(s),s∈[0,T]}∈AT∩LT2f = \left\{ f(s), s \in [0, T]\right\} \in \mathscr {A}_{T} \cap L_{T}^{2}, что

    X=E[X]+∫0Tf(s,ω)dWsп.н. X = \mathbb {E}\left[X\right]+\int _{0}^{T} f(s, \omega ) d W_{s} \quad \text{п.н.}
  2. Если E[∣X∣]<∞\mathbb {E}\left[\left|X\right|\right] < \infty, то представление выше справедливо с некоторым процессом f∈J1([0,T])f \in J_{1}([0, T]).

  3. Если XX — положительная случайная величина (т.е. P(X>0)=1\mathbb {P}\left(X > 0\right) = 1) и E[X]<∞\mathbb {E}\left[X\right] < \infty, то найдется процесс f∈J1[0,T]f \in J_{1}[0, T] такой, что XX допускает представление вида X=ZTE[X]X = Z_{T} \mathbb {E}\left[X\right], где ZT=exp(∫0Tf(s,ω)dWs−12∫0Tf2(s,ω)ds)Z_{T} = \mathrm{exp}\left(\int_{0}^{T} f(s, \omega ) d W_{s}-\frac{1}{2} \int_{0}^{T} f^{2}(s, \omega ) d s\right).

Доказательства теоремы 11 в различных вариантах даны К. Ито, Дж. Кларком, Дж. Дубом и приведены во многих книгах, см., например, [46,182][46,182].

В качестве примера использования представления функционалов от броуновского движения в виде стохастических интегралов приведем следующую экстремальную задачу:

среди всех моментов остановки τ\tau (относительно естественной фильтрации броуновского движения W={W(t),t∈[0,1]})W = \left\{ W(t), t \in [0,1]\right\} ) найти такой момент τ∗\tau_{*}, для которого достигается

V∗=inf⁡0⩽τ⩽1E[∣Wτ−max⁡0⩽s⩽1Ws∣2] V_{*} = \inf _{0 \leqslant \tau \leqslant 1} \mathbb {E}\left[\left|W_{\tau }-\max _{0 \leqslant s \leqslant 1} W_{s}\right|^{2}\right]

Можно представить себе модель, в которой колебания цен на акции в течение единичного периода времени (например, день) описываются броуновским движением (разумеется, это лишь интерпретация, так как цены не бывают отрицательными) и требуется выбрать момент τ∗\tau_{*} для самой выгодной в среднем квадратическом смысле продажи имеющихся акций. Заметим, что Wτ∗W_{\tau_{*}} дает смещенную оценку для max⁡0⩽s⩽1Ws\max_{0 \leqslant s \leqslant 1} W_{s}, поскольку E[Wτ∗]=0\mathbb {E}\left[W_{\tau_{*}}\right] = 0 и E[max⁡0⩽s⩽1Ws]=2/π\mathbb {E}\left[\max_{0 \leqslant s \leqslant 1} W_{s}\right] = \sqrt{2 / \pi } (поясните в качестве простого упражнения). Поэтому можно обратиться к величине

V~=inf⁡a∈R,τ∈[0,1]E[∣Wτ+a−max⁡0⩽s⩽1Ws∣2] \tilde{V} = \inf _{a \in \mathbb {R}, \tau \in [0,1]} \mathbb {E}\left[\left|W_{\tau }+a-\max _{0 \leqslant s \leqslant 1} W_{s}\right|^{2}\right]

Легко проверить (убедитесь в этом), что

V~=V∗−2/π \widetilde{V} = V_{*}-2 / \pi

Положим

St=max⁡0⩽s⩽tWs,t∈[0,1] S_{t} = \max _{0 \leqslant s \leqslant t} W_{s}, \quad t \in [0,1]

Теорема 12 (Граверсен, Пешкир, Ширяев). Решение приведенной выше задачи дается формулой

τ∗=inf⁡{t∈[0,1]:St−Wt=z∗1−t} \tau _{*} = \inf \left\{ t \in [0,1]: S_{t}-W_{t} = z_{*} \sqrt{1-t}\right\}

где константа z∗z_{*} находится из уравнения

4Φ(z∗)−2z∗ϕ(z∗)−3=0 4 \Phi \left(z_{*}\right)-2 z_{*} \phi \left(z_{*}\right)-3 = 0

здесь Φ\Phi и ф соответственно функция распределения и плотность стандартной нормальной величины. При этом z∗=1.12…,V∗=2Φ(z∗)−1=0.73…z_{*} = 1.12 \ldots , V_{*} = 2 \Phi \left(z_{*}\right)-1 = 0.73 \ldots.

Ключевую роль при доказательстве этого результата играет представление вида

max⁡0⩽s⩽1Ws=b+∫01f(s,ω)dWs \max _{0 \leqslant s \leqslant 1} W_{s} = b+\int _{0}^{1} f(s, \omega ) d W_{s}

где b=b = const (укажите эту константу) и

f(s,ω)=2{1−Φ((St−Wt)/1−t)},s∈[0,1],ω∈Ω f(s, \omega ) = 2\left\{ 1-\Phi \left(\left(S_{t}-W_{t}\right) / \sqrt{1-t}\right)\right\} , \quad s \in [0,1], \quad \omega \in \Omega

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

Определение 12. Согласованный с фильтрацией F=(Ft)t⩾0\mathbb {F} = \left(\mathscr {F}_{t}\right)_{t \geqslant 0} процесс A={At;t∈R+}A = \left\{ A_{t} ; t \in \mathbb {R}_{+}\right\}называется возрастающим, если для п.в. ω\omega имеем A0(ω)=0A_{0}(\omega ) = 0, где At(ω)A_{t}(\omega ) — неубывающая функция по t∈[0,∞)t \in [0, \infty ) и E[At]<∞\mathbb {E}\left[A_{t}\right] < \infty для каждого t∈R+t \in \mathbb {R}_{+}.

Определение 13. Действительный случайный процесс Y={Yt,t⩾0}Y = \left\{ Y_{t}, t \geqslant 0\right\} принадлежит классу Дирихле ( DD ), если семейство величин {Yτ,τ∈S}\left\{ Y_{\tau }, \tau \in S\right\} равномерно интегрируемо, где SS обозначает совокупность всех конечных марковских моментов (относительно фильтрации F\mathbb {F} ).

Справедлив следующий фундаментальный результат (см. [26, т. 1, с. 68]; сравните с теоремой 1 главы IV).

Теорема 13 (разложение Дуба-Мейера). Субмартингал X=(Xt,Ft)t⩾0X = \left(X_{t}, \mathscr {F}_{t}\right)_{t \geqslant 0} класса (D) с траекториями типа càdlàg допускает и притом единственное (с точностью до стохастической неразличимости) разложение вида

Xt=X0+At+Mt X_{t} = X_{0}+A_{t}+M_{t}

где A=(At,Ft)t⩾0A = \left(A_{t}, \mathscr {F}_{t}\right)_{t \geqslant 0} — возрастающий предсказуемый интегрируемый (E[A∞]<∞\mathbb {E}\left[A_{\infty }\right] < \infty, где A∞=lim⁡t→∞AtA_{\infty } = \lim_{t \rightarrow \infty } A_{t}) процесс, а M=(Mt,Ft)t⩾0M = \left(M_{t}, \mathscr {F}_{t}\right)_{t \geqslant 0} есть равномерно интегрируемый мартингал.

Рассмотрим для простоты непрерывный квадратично интегрируемый мартингал M=(Mt)t⩾0M = \left(M_{t}\right)_{t \geqslant 0} на фильтрованном вероятностном пространстве (Ω,F,(Ft)t⩾0,P)\left(\Omega , \mathscr {F},\left(\mathscr {F}_{t}\right)_{t \geqslant 0}, \mathbb {P}\right) такой, что M0=0M_{0} = 0. Класс таких мартингалов обозначается M2c\mathscr {M}_{2}^{c}. Будем предполагать, что фильтрация F=(Ft)t⩾0\mathbb {F} = \left(\mathscr {F}_{t}\right)_{t \geqslant 0} удовлетворяет обычным условиям.

Наметим схему построения интеграла по процессу M=(Mt,Ft)t⩾0∈M2cM = \left(M_{t}, \mathscr {F}_{t}\right)_{t \geqslant 0} \in \mathscr {M}_{2}^{c}. Пусть ⟨M⟩=(⟨M⟩t)t⩾0\left\langle M\right\rangle = \left(\left\langle M\right\rangle_{t}\right)_{t \geqslant 0} — возрастающий предсказуемый процесс, участвующий в разложении Дуба-Мейера субмартингала Mt2,t⩾0M_{t}^{2}, t \geqslant 0, обладающий тем свойством, что E[⟨M⟩∞]<∞\mathbb {E}\left[\left\langle M\right\rangle_{\infty }\right] < \infty.

Как и в случае броуновского движения, интеграл

IT(X)=∫0TXt(ω)dMt(ω) I_{T}(X) = \int _{0}^{T} X_{t}(\omega ) d M_{t}(\omega )

вообще говоря, не может быть определен потраекторно, как интеграл Лебега-Стилтьеса (можно определить в тривиальном случае Mt=0M_{t} = 0 при t⩾0t \geqslant 0 ).

Процесс X={Xt,t⩾0,ω∈Ω}X = \left\{ X_{t}, t \geqslant 0, \omega \in \Omega \right\} назовем простым, если существуют строго возрастающая числовая последовательность {tn}n=0∞\left\{ t_{n}\right\}_{n = 0}^{\infty } с t0=0t_{0} = 0 и lim⁡n→∞tn=∞\lim_{n \rightarrow \infty } t_{n} = \infty, а также последовательность действительных случайных величин {fn(ω)}n=0∞\left\{ f_{n}(\omega )\right\}_{n = 0}^{\infty } и некоторая константа 0<c<∞0 < c < \infty, для которых sup⁡∣fn(ω)∣⩽c\sup \left|f_{n}(\omega )\right| \leqslant c при всех ω∈Ω\omega \in \Omega, причем fn∈Ftn∣B(R)f_{n} \in \mathscr {F}_{t_{n}} \mid \mathscr {B}(\mathbb {R}) и

Xt(ω)=f0(ω)1{0}(t)+∑k=0∞fk(ω)1(tk,tk+1](t),0⩽t<∞ X_{t}(\omega ) = f_{0}(\omega ) \mathbf{1}_{\left\{ 0\right\} }(t)+\sum _{k = 0}^{\infty } f_{k}(\omega ) \mathbf{1}_{\left(t_{k}, t_{k+1}\right]}(t), \quad 0 \leqslant t < \infty

Класс простых процессов обозначим L0\mathscr {L}_{0} и для X∈L0X \in \mathscr {L}_{0} положим

It(X)=∑k=0n−1fk(Mtk+1−Mtk)+fn(Mt−Mtn) при tn⩽t<tn+1 I_{t}(X) = \sum _{k = 0}^{n-1} f_{k}\left(M_{t_{k+1}}-M_{t_{k}}\right)+f_{n}\left(M_{t}-M_{t_{n}}\right) \text{ при } t_{n} \leqslant t < t_{n+1}

Далее доказывается, что определение (92) можно распространить (снова используя должные аппроксимации простыми процессами) с L0\mathscr {L}_{0} на более широкий класс процессов. Интересно отметить, что этот класс будет зависеть от того, является ли функция ⟨M⟩t\left\langle M\right\rangle_{t} для п.в. ω\omega абсолютно непрерывной относительно меры Лебега. Детальное построение и дальнейшие обобщения можно почерпнуть в книгах [46,182][46,182].

?
Задача VIII.17

С помощью формулы Ито докажите, что для детерминированной функции f=f(t)f = f(t) из C1C^{1} справедлива формула "интегрирования по частям":

∫(0,t]f(s)dWs=f(t)Wt−∫(0,t]f′(s)Wsds \int _{(0, t]} f(s) d W_{s} = f(t) W_{t}-\int _{(0, t]} f^{\prime }(s) W_{s} d s

Постройте пример случайной функции f=f(t,ω)f = f(t, \omega ), для которой формула (93) не имеет места.

Пример 2. Рассмотрим стохастически возмущенное уравнение роста популяции

dXtdt=rXt,r= const.  \frac{d X_{t}}{d t} = r X_{t}, \quad r = \text{ const. }

Точнее говоря, рассмотрим стох астическое дифференциальное уравнение вида

dXt=rXtdt+σXtdWt d X_{t} = r X_{t} d t+\sigma X_{t} d W_{t}

с начальным условием X0(X0∈F∣B(R)X_{0}\left(X_{0} \in \mathscr {F} \mid \mathscr {B}(\mathbb {R})\right. и E[X02]<∞)\left.\mathbb {E}\left[X_{0}^{2}\right] < \infty \right), где константы σ>0,r∈R\sigma > 0, r \in \mathbb {R}.

Из формулы Ито непосредственно следует, что процесс

Xt=X0e(r−σ2/2)t+σWt X_{t} = X_{0} e^{\left(r-\sigma ^{2} / 2\right) t+\sigma W_{t}}

является сильным решением уравнения (95). Поскольку X0  ⊥ ⁣ ⁣ ⁣⊥  W={Wt,t⩾0}X_{0} \; \perp \! \! \! \perp \; W = \left\{ W_{t}, t \geqslant 0\right\}, то

E[Xt]=E[X0]exp(rt) \mathbb {E}\left[X_{t}\right] = \mathbb {E}\left[X_{0}\right] \mathrm{exp}\left(r t\right)

Иначе говоря, при E[X0]≠0\mathbb {E}\left[X_{0}\right] \neq 0 средний рост (убывание, если r<0r < 0 ) E[Xt]\mathbb {E}\left[X_{t}\right] будет таким же, как в случае детерминированной модели (94).

Уравнение (95) играет важную роль в финансовой математике. Его решение (96) часто называют "геометрическим броуновским движением". Локальный снос r−σ2/2r-\sigma^{2} / 2 характеризует скорость изменения среднего значения процесса XX, диффузию σ2\sigma^{2} в финансовой литературе часто называют волатильностью. Видимо, П. Самуэльсон был первым, кто осознал важность геометрического броуновского движения для описания эволюции цен, используя для него термин "экономическое броуновское движение".

Определения и результаты, излагавшиеся в главе VIII, ценой некоторых усложнений естественным образом переносятся на многомерный случай. Так, уравнение (59) можно понимать как систему, в которой Xt=(Xt(1),…,Xt(n)),bX_{t} = \left(X_{t}^{(1)}, \ldots , X_{t}^{(n)}\right), b — векторная, σ\sigma — матричная функции, точнее говоря,

b(⋅,⋅):[0,T]×Rn→Rn,σ(⋅,⋅):[0,T]×Rn→Rn×m b(\cdot , \cdot ):[0, T] \times \mathbb {R}^{n} \rightarrow \mathbb {R}^{n}, \quad \sigma (\cdot , \cdot ):[0, T] \times \mathbb {R}^{n} \rightarrow \mathbb {R}^{n \times m}

W={Wt,t⩾0}W = \left\{ W_{t}, t \geqslant 0\right\} — mm-мерное броуновское движение (относительно фильтрации (Ft)t⩾0)\left.\left(\mathscr {F}_{t}\right)_{t \geqslant 0}\right). Считаем также, что

b(t,Xt)dt=(b(1)(t,Xt)dt⋮b(n)(t,Xt)dt),σ(t,Xt)dWt=(σ11(t,Xt)…σ1m(t,Xt)⋮⋮σn1(t,Xt)…σnm(t,Xt))(dWt(1)⋮dWt(m)). \begin{gathered} b\left(t, X_{t}\right) d t =\begin{pmatrix} b^{(1)}\left(t, X_{t}\right) d t \\ \vdots \\ b^{(n)}\left(t, X_{t}\right) d t \end{pmatrix}, \\ \sigma \left(t, X_{t}\right) d W_{t} =\begin{pmatrix} \sigma _{11}\left(t, X_{t}\right) & \ldots & \sigma _{1 m}\left(t, X_{t}\right) \\ \vdots & & \vdots \\ \sigma _{n 1}\left(t, X_{t}\right) & \ldots & \sigma _{n m}\left(t, X_{t}\right) \end{pmatrix}\begin{pmatrix} d W_{t}^{(1)} \\ \vdots \\ d W_{t}^{(m)} \end{pmatrix} . \end{gathered}

Интеграл от векторной функции определяется как вектор из интегралов от компонент, а ∫[0,t]σ(s,Xs)dWs\int_{[0, t]} \sigma \left(s, X_{s}\right) d W_{s} обозначает вектор-функцию с ii-й компонентой, равной

∑k=1m∫[0,t]σik(s,Xs)dWs(k) \sum _{k = 1}^{m} \int _{[0, t]} \sigma _{i k}\left(s, X_{s}\right) d W_{s}^{(k)}

Определение прогрессивной измеримости очевидным образом переносится на многомерный случай. Теорема 13 остается в силе, если в условиях (61) и (62) модуль bb понимать как норму в Rn\mathbb {R}^{n}, а модуль σ\sigma понимать как норму матрицы, например, считать, что ∣σ∣2=∑i=1n∑k=1mσik2\left|\sigma \right|^{2} = \sum_{i = 1}^{n} \sum_{k = 1}^{m} \sigma_{i k}^{2}.

При исследовании систем стохастических дифференциальных уравнений важнейшую роль играют многомерные варианты формулы Ито, т. е. формулы, описывающие стохастический дифференциал процесса H(t,Xt)H\left(t, X_{t}\right), где H=H(t,x1,…,xn)H = H\left(t, x_{1}, \ldots , x_{n}\right), а X=(X(1),…,X(n))X = \left(X^{(1)}, \ldots , X^{(n)}\right) — векторный процесс Ито. Другими словами, XX имеет стохастический дифференциал вида

dXt=f(t,ω)dWt+g(t,ω)dt,t∈[0,∞) d X_{t} = f(t, \omega ) d W_{t}+g(t, \omega ) d t, \quad t \in [0, \infty )

где ff и gg — соответственно матричная и векторная прогрессивно измеримые функции, f=(fik(t,ω),i=1,…,n;k=1,…,m);g=(g1(t,ω),…,gn(t,ω)),Wt=(Wt(1),…,Wt(m))f = \left(f_{i k}(t, \omega ), i = 1, \ldots , n ; k = 1, \ldots , m\right) ; g = \left(g_{1}(t, \omega ), \ldots , g_{n}(t, \omega )\right), W_{t} = \left(W_{t}^{(1)}, \ldots , W_{t}^{(m)}\right). Иначе говоря, пусть

dXt(i)=∑k=1mfik(t,ω)dWt(k)+gi(t,ω)dt,i=1,…,n d X_{t}^{(i)} = \sum _{k = 1}^{m} f_{i k}(t, \omega ) d W_{t}^{(k)}+g_{i}(t, \omega ) d t, \quad i = 1, \ldots , n

где при любом t>0t > 0, всех i=1,…,ni = 1, \ldots , n и k=1,…,mk = 1, \ldots , m

P(∫[0,t]∣gi(s,ω)∣ds<∞)=1 и P(∫[0,t]∣fik(s,ω)∣2ds<∞)=1.  \mathbb {P}\left(\int _{[0, t]}\left|g_{i}(s, \omega )\right| d s < \infty \right) = 1 \text{ и } \mathbb {P}\left(\int _{[0, t]}\left|f_{i k}(s, \omega )\right|^{2} d s < \infty \right) = 1 \text{. }

Теорема 14. Пусть функция H:[0,∞)×Rn→RH:[0, \infty ) \times \mathbb {R}^{n} \rightarrow \mathbb {R} такова, что H∈C1,2H \in C^{1,2}, т. е. существуют непрерывные производные ∂H/∂t\partial H / \partial t и ∂2H/∂xi∂xj,i,j=1,…,n\partial^{2} H / \partial x_{i} \partial x_{j}, i, j = 1, \ldots , n. Тогда процесс Yt=H(t,Xt)Y_{t} = H\left(t, X_{t}\right), дде процесс Ито XtX_{t} определяется согласно (100), (102), имеет стохастический дифференииал, который задается формулой

dYt=∂H∂t(t,Xt)dt+∑i=1n∂H∂xi(t,Xt)dXt(i)+12∑i=1n∑j=1n∂2H∂xi∂xj(t,Xt)dXt(i)dXt(j) d Y_{t} = \frac{\partial H}{\partial t}\left(t, X_{t}\right) d t+\sum _{i = 1}^{n} \frac{\partial H}{\partial x_{i}}\left(t, X_{t}\right) d X_{t}^{(i)}+\frac{1}{2} \sum _{i = 1}^{n} \sum _{j = 1}^{n} \frac{\partial ^{2} H}{\partial x_{i} \partial x_{j}}\left(t, X_{t}\right) d X_{t}^{(i)} d X_{t}^{(j)}

здесь dXt(i)d X_{t}^{(i)} и d Xt(j)X_{t}^{(j)}, фигурирующие в (101), “перемножаются по правилу":

dWt(i)⋅dWt(j)=δijdt,dt⋅dt=dWt(i)⋅dt=dt⋅dWt(i)=0 d W_{t}^{(i)} \cdot d W_{t}^{(j)} = \delta _{i j} d t, \quad d t \cdot d t = d W_{t}^{(i)} \cdot d t = d t \cdot d W_{t}^{(i)} = 0

Доказательство формулы Ито см., например, в [26]. По поводу разнообразных обобщений этой формулы см., например, [11,146][11,146].

Приведем здесь лишь один результат в направлении упомянутых обобщений. Пусть XX — стандартное (одномерное) броуновское движение W={Wt,t⩾0}W = \left\{ W_{t}, t \geqslant 0\right\}. Пусть функция F=F(x)F = F(x) является абсолютно непрерывной,

F(x)=F(0)+∫0xf(y)dy F(x) = F(0)+\int _{0}^{x} f(y) d y

причем функция f=f(y)∈Lloc2(R)f = f(y) \in L_{\mathrm{loc}}^{2}(\mathbb {R}), т.е. для любого c>0c > 0

∫∣y∣⩽cf2(y)dy<∞ \int _{\left|y\right| \leqslant c} f^{2}(y) d y < \infty

Определим [f(W),W][f(W), W] — квадратическую ковариацию процессов f(W)f(W) и WW следующим образом:

[f(W),W]t=P−lim⁡n∑m(f(Wt(n)(m+1)∧t)−f(Wt(n)(m)∧t))××(Wt(n)(m+1)∧t−Wt(n)(m)∧t),t⩾0 \begin{aligned} & {[f(W), W]_{t} = \mathbb {P}-\lim _{n} \sum _{m}\left(f\left(W_{t^{(n)}(m+1) \wedge t}\right)-f\left(W_{t^{(n)}(m) \wedge t}\right)\right) \times } \\ & \times \left(W_{t^{(n)}(m+1) \wedge t}-W_{t^{(n)}(m) \wedge t}\right), \quad t \geqslant 0 \end{aligned}

Здесь P\mathbb {P} - lim⁡\lim обозначает предел по вероятности, а {t(n)(m),m∈N}\left\{ t^{(n)}(m), m \in \mathbb {N}\right\} при каждом n∈Nn \in \mathbb {N} — римановские последовательности детерминированных моментов t(n)(m)t^{(n)}(m), т.е. t(n)(m)⩽t(n)(m+1)t^{(n)}(m) \leqslant t^{(n)}(m+1) для m∈Nm \in \mathbb {N} при каждом n∈Nn \in \mathbb {N} и t(n)(m)→∞t^{(n)}(m) \rightarrow \infty при m→∞m \rightarrow \infty; кроме того, для любого t>0t > 0

sup⁡m(t(n)(m+1)∧t−t(n)(m)∧t)→0,n→∞ \sup _{m}\left(t^{(n)}(m+1) \wedge t-t^{(n)}(m) \wedge t\right) \rightarrow 0, \quad n \rightarrow \infty

Важно подчеркнуть, что поскольку процесс f(W)f(W) не является, вообще говоря, семимартингалом, нетривиален факт существования предела по вероятности, задающего [f(W),W]t[f(W), W]_{t}. Один из результатов работы [129][129] состоит именно в доказательстве существования этого предела.

Теорема 15 (Проттер, Фёльмер, Ширяев; [129]). При сделанных выше предположениях о функции F=F(x)F = F(x) имеет место формула

F(Wt)=F(0)+∫0tf(Ws)dWs+12[f(W),W]t F\left(W_{t}\right) = F(0)+\int _{0}^{t} f\left(W_{s}\right) d W_{s}+\frac{1}{2}[f(W), W]_{t}
?
Задача VIII.18

Пусть функция f∈C2f \in C^{2}. Докажите, что

[f(W),W]t=∫0tf′(Ws)ds,t⩾0 [f(W), W]_{t} = \int _{0}^{t} f^{\prime }\left(W_{s}\right) d s, \quad t \geqslant 0

Тем самым, формула (103) переходит в формулу Ито.

?
Задача VIII.19

Пусть f(x)=∣x∣f(x) = \left|x\right|. Докажите, что [f(W),W]t=2Lt(0)[f(W), W]_{t} = 2 L_{t}(0), где Lt(0)L_{t}(0) — локальное время (П. Леви) броуновского движения в точке 0 , определяемое следующим образом:

Lt=lim⁡ε↓012ε∫0t1{∣Ws∣⩽ε}ds L_{t} = \lim _{\varepsilon \downarrow 0} \frac{1}{2 \varepsilon } \int _{0}^{t} \mathbf{1}\left\{ \left|W_{s}\right| \leqslant \varepsilon \right\} d s

Таким образом, в данном случае формула (103) переходит в формулу Танака для броуновского движения.

Обобщению интеграла Ито на функции многих переменньх посвящена монография [163]. Приведем только один результат этого направления.

Теорема 16 (Ито). При всех t>0t > 0 и n∈Nn \in \mathbb {N} справедлива формула

∫⋯∫0⩽s1⩽⋯⩽sn⩽tdWs1⋯dWsn=tn/2n!Hn(Wtt) \int \cdots \int _{0 \leqslant s_{1} \leqslant \cdots \leqslant s_{n} \leqslant t} d W_{s_{1}} \cdots d W_{s_{n}} = \frac{t^{n / 2}}{n !} H_{n}\left(\frac{W_{t}}{\sqrt{t}}\right)

где HnH_{n} — полином Эрмита степени п, т.е.

Hn(x)=(−1)nex2/2dndxn(e−x2/2),n=0,1,… H_{n}(x) = (-1)^{n} e^{x^{2} / 2} \frac{d^{n}}{d x^{n}}\left(e^{-x^{2} / 2}\right), \quad n = 0,1, \ldots

Обратимся теперь к однородной диффузии, т.е. рассмотрим решение уравнения

dXt=b(Xt)dt+σ(Xt)dWt,t⩾s,Xs=x∈Rn d X_{t} = b\left(X_{t}\right) d t+\sigma \left(X_{t}\right) d W_{t}, \quad t \geqslant s, \quad X_{s} = x \in \mathbb {R}^{n}

где WtW_{t} — mm-мерный винеровский процесс, а функции b:Rn→Rnb: \mathbb {R}^{n} \rightarrow \mathbb {R}^{n} и σ:Rn→Rn×m\sigma : \mathbb {R}^{n} \rightarrow \mathbb {R}^{n \times m} удовлетворяют условиям теоремы 8 , которые в данном случае сводятся к одному требованию: существует L>0L > 0 такое, что

∣b(x)−b(y)∣+∣σ(x)−σ(y)∣⩽L∣x−y∣,x,y∈Rn \left|b(x)-b(y)\right|+\left|\sigma (x)-\sigma (y)\right| \leqslant L\left|x-y\right|, \quad x, y \in \mathbb {R}^{n}

(для вектора ∣⋅∣\left|\cdot \right| обозначает евклидову норму и ∣σ∣2=∑i=1n∑k=1mσik2\left|\sigma \right|^{2} = \sum_{i = 1}^{n} \sum_{k = 1}^{m} \sigma_{i k}^{2} ), из которого вытекает, что для некоторого c>0c > 0

∣b(x)∣2+∣σ(x)∣2⩽c(1+∣x∣2),x∈Rn \left|b(x)\right|^{2}+\left|\sigma (x)\right|^{2} \leqslant c\left(1+\left|x\right|^{2}\right), \quad x \in \mathbb {R}^{n}

Обозначим Xts,xX_{t}^{s, x} единственное сильное решение уравнения (106) при t⩾st \geqslant s.

?
Задача VIII.20

Докажите, что для каждого x∈Rnx \in \mathbb {R}^{n} определенный выше процесс Xts,x,t⩾sX_{t}^{s, x}, t \geqslant s, является однородным марковским процессом.

?
Задача VIII.21

Докажите, что при выполнении условий теоремы 8 единственность сильного решения (на [0,T][0, T] ) будет в более широком классе процессов, у которых существует E[Xt2]\mathbb {E}\left[X_{t}^{2}\right] для t∈[0,T]t \in [0, T].

Теорема 8. Рассмотрим стохастическое дифференциальное уравнение

Xt=X0+∫0tb(s,Xs)ds+∫0tσ(s,Xs)dWs,0⩽t⩽T, X_{t} = X_{0}+\int _{0}^{t} b(s, X_{s}) d s+\int _{0}^{t} \sigma (s, X_{s}) d W_{s}, \quad 0 \leqslant t \leqslant T,

где bb и σ\sigma определены на [0,T]×R[0, T] \times \mathbb {R}. Пусть выполнено условие Липшица: существует константа L=L(T)>0L = L(T) > 0 такая, что

∣b(t,x)−b(t,y)∣+∣σ(t,x)−σ(t,y)∣⩽L∣x−y∣,x,y∈R,t∈[0,T], \left|b(t, x)-b(t, y)\right|+\left|\sigma (t, x)-\sigma (t, y)\right| \leqslant L\left|x-y\right|, \quad x, y \in \mathbb {R}, \quad t \in [0, T],

и условие линейного роста: для некоторого c=c(T)>0c = c(T) > 0

∣b(t,x)∣+∣σ(t,x)∣⩽c(1+∣x∣),x∈R,t∈[0,T]. \left|b(t, x)\right|+\left|\sigma (t, x)\right| \leqslant c(1+\left|x\right|), \quad x \in \mathbb {R}, \quad t \in [0, T].

Пусть величина X0X_{0} является F0\mathscr {F}_{0}-измеримой, причем E[X02]<∞\mathbb {E}\left[X_{0}^{2}\right] < \infty. Тогда существует единственное сильное решение уравнения с F0\mathscr {F}_{0}-измеримым начальным условием X0X_{0} такое, что Xt∈L2(Ω)X_{t} \in L^{2}(\Omega ) для любого t∈[0,T]t \in [0, T] и функция E[Xt2]\mathbb {E}\left[X_{t}^{2}\right] ограничена на [0,T][0, T].

?
Задача VIII.22

Покажите, что если условия теоремы 8 выполнены на полуинтервале [0,T)[0, T), то ее утверждение также справедливо на этом промежутке, только в этом случае рассматривается класс процессов XX, для которых функция E[Xt2]\mathbb {E}\left[X_{t}^{2}\right] будет ограничена на каждом отрезке, вложенном в [0,T)[0, T).

Теорема 8. Рассмотрим стохастическое дифференциальное уравнение

Xt=X0+∫0tb(s,Xs)ds+∫0tσ(s,Xs)dWs,0⩽t⩽T, X_{t} = X_{0}+\int _{0}^{t} b(s, X_{s}) d s+\int _{0}^{t} \sigma (s, X_{s}) d W_{s}, \quad 0 \leqslant t \leqslant T,

где bb и σ\sigma определены на [0,T]×R[0, T] \times \mathbb {R}. Пусть выполнено условие Липшица: существует константа L=L(T)>0L = L(T) > 0 такая, что

∣b(t,x)−b(t,y)∣+∣σ(t,x)−σ(t,y)∣⩽L∣x−y∣,x,y∈R,t∈[0,T], \left|b(t, x)-b(t, y)\right|+\left|\sigma (t, x)-\sigma (t, y)\right| \leqslant L\left|x-y\right|, \quad x, y \in \mathbb {R}, \quad t \in [0, T],

и условие линейного роста: для некоторого c=c(T)>0c = c(T) > 0

∣b(t,x)∣+∣σ(t,x)∣⩽c(1+∣x∣),x∈R,t∈[0,T]. \left|b(t, x)\right|+\left|\sigma (t, x)\right| \leqslant c(1+\left|x\right|), \quad x \in \mathbb {R}, \quad t \in [0, T].

Пусть величина X0X_{0} является F0\mathscr {F}_{0}-измеримой, причем E[X02]<∞\mathbb {E}\left[X_{0}^{2}\right] < \infty. Тогда существует единственное сильное решение уравнения с F0\mathscr {F}_{0}-измеримым начальным условием X0X_{0} такое, что Xt∈L2(Ω)X_{t} \in L^{2}(\Omega ) для любого t∈[0,T]t \in [0, T] и функция E[Xt2]\mathbb {E}\left[X_{t}^{2}\right] ограничена на [0,T][0, T].

?
Задача VIII.23

С помощью формулы Ито докажите, что решением (одномерного) уравнения

dXt=β−XtT−tdt+dWt,t∈[0,T),X0=α d X_{t} = \frac{\beta -X_{t}}{T-t} d t+d W_{t}, \quad t \in [0, T), \quad X_{0} = \alpha

является процесс

Xt=α(1−t/T)+βt/T+(T−t)∫0tdWsT−s X_{t} = \alpha (1-t / T)+\beta t / T+(T-t) \int _{0}^{t} \frac{d W_{s}}{T-s}
?
Задача VIII.24

(продолжение упражнения 23). Докажите, что Xt→βX_{t} \rightarrow \beta п.н. при t→T−t \rightarrow T-. Таким образом, решение уравнения (109) представляет собой броуновский мост над отрезком [0,T][0, T] с закрепленными концами X0=αX_{0} = \alpha и XT=βX_{T} = \beta. Стандартный броуновский мост получается при T=1T = 1 и α=β=0\alpha = \beta = 0.

Из различного рода обобщений теоремы существования и единственности решения стохастического дифференциального уравнения приведем несколько неожиданный результат А.К. Звонкина (см., например, [86; с. 322]), утверждающий, что для существования сильного решения стохастического дифференциального уравнения

dXt=b(t,Xt)dt+dWt d X_{t} = b\left(t, X_{t}\right) d t+d W_{t}

достаточно лишь измеримости по (t,x)(t, x) и равномерной ограниченности b(t,x)b(t, x).

?
Задача VIII.25

Докажите, что стохастическое дифференциальное уравнение

dXt=σ(Xt)dt+dWt d X_{t} = \sigma \left(X_{t}\right) d t+d W_{t}

с "плохим" коэффициентом σ(x)=sgn⁡x\sigma (x) = \operatorname {sgn} x имеет и, притом, единственное сильное решение.

Определение 14. Стохастическое дифференциальное уравнение (59) с начальным условием μ\mu, где μ\mu — заданная мера на B(R)\mathscr {B}(\mathbb {R}), имеет слабое решение на промежутке [0,T][0, T], если найдутся фильтрованное вероятностное пространство (Ω,F,(Ft)t∈[0,T],P)\left(\Omega , \mathscr {F},\left(\mathscr {F}_{t}\right)_{t \in [0, T]}, \mathbb {P}\right), броуновское движение W=(Wt,Ft)t∈[0,T]W = \left(W_{t}, \mathscr {F}_{t}\right)_{t \in [0, T]} и непрерывньй (п.н.) случайньй процесс X=(Xt,Ft)t∈[0,T]X = \left(X_{t}, \mathscr {F}_{t}\right)_{t \in [0, T]} такие, что Law⁡(X0)=μ\operatorname {Law}\left(X_{0}\right) = \mu и Р-п. н. для каждого t>0t > 0 выполнено равенство (60).

Разумеется, вместо отрезка [0,T][0, T] можно рассматривать промежутки вида [u,v][u, v], где 0⩽u<v<∞0 \leqslant u < v < \infty, или [u,v)[u, v) для 0⩽u<v⩽∞0 \leqslant u < v \leqslant \infty.

Важно подчеркнуть, что в отличие от сильного решения, рассматриваемого на заданном фильтрованном вероятностном пространстве с заданным на нем броуновским движением, в определении слабого решения такие объекты (вероятностное пространство и броуновское движение) не фиксируются, а требуется лишь, чтобы они нашлись. Очевидно, сильное решение является и слабым.

Определение 15. Слабая единственность (слабого или сильного) решения уравнения (59) означает, что у любых двух решений (в рассматриваемом смысле) с одинаковыми начальными условиями совпадают распределения, т. е. совпадают конечномерные распределения.

?
Задача VIII.26

Докажите, что уравнение

dXt=σ(Xt)dWt,t∈[0,1],X0=0 d X_{t} = \sigma \left(X_{t}\right) d W_{t}, \quad t \in [0,1], \quad X_{0} = 0

где σ(x)=sgn⁡x\sigma (x) = \operatorname {sgn} x, имеет по крайней мере два решения (но "слабых"). Кроме того, на некоторых вероятностных пространствах у этого уравнения может вовсе и не быть "сильного" решения.

Для выполнения последнего упражнения следует убедиться, что процесс

Bt=∫0tσ(Ws)dWs,t∈[0,1],σ(x)=sgn⁡x B_{t} = \int _{0}^{t} \sigma \left(W_{s}\right) d W_{s}, \quad t \in [0,1], \quad \sigma (x) = \operatorname {sgn} x

является броуновским движением на отрезке [0,1][0,1]. Это утверждение обеспечивает (см., например, [83]) следующая

Теорема 17 (П. Леви). Пусть B=(Bt,Ft)t⩾0B = \left(B_{t}, \mathscr {F}_{t}\right)_{t \geqslant 0} — P\mathbb {P}-п.н. непрерывный квадратично интегрируемый мартингал, заданный на некотором фильтрованном вероятностном пространстве (Ω,F,(Ft)t⩾0,P)\left(\Omega , \mathscr {F},\left(\mathscr {F}_{t}\right)_{t \geqslant 0}, \mathbb {P}\right). Пусть (Bt2−t,Ft)t⩾0\left(B_{t}^{2}-t, \mathscr {F}_{t}\right)_{t \geqslant 0} также является мартингалом, т. е.

E[Bt2−Bs2∣Fs]=t−s,0⩽s⩽t. \mathbb {E}\left[B_{t}^{2}-B_{s}^{2} \mid \mathscr {F}_{s}\right] = t-s, \quad 0 \leqslant s \leqslant t .

Тогда B={Bt,t⩾0}B = \left\{ B_{t}, t \geqslant 0\right\} есть стандартное броуновское движение.

Заметим, что М. Барлоу [92] доказал, что уравнение (113) может не иметь сильного решения даже в случае непрерывных ограниченных функций σ=σ(x)>0\sigma = \sigma (x) > 0.

?
Задача VIII.27

Пусть выполнены условия теоремы 8. Докажите, что решение (слабое или сильное) уравнения (59) слабо единственно.

Упражнение 26 показывает, что могут существовать слабые решения, но не быть сильных. Поэтому естественно ожидать существования слабого решения при менее ограничительных условиях на коэффициенты уравнения (59).

Один из первых результатов в этом направлении (см. об этом [86; т. 1, с. 325]) формулируется следуюшим образом. Рассмотрим стохастическое дифференциальное уравнение

dXt=b(Xt)dt+σ(Xt)dWt d X_{t} = b\left(X_{t}\right) d t+\sigma \left(X_{t}\right) d W_{t}

с начальным распределением μ=Law⁡(X0)\mu = \operatorname {Law}\left(X_{0}\right) таким, что

∫R∣x∣γμ(dx)=E[∣X0∣γ]<∞ для некоторого γ>2.  \int _{\mathbb {R}}\left|x\right|^{\gamma } \mu (d x) = \mathbb {E}\left[\left|X_{0}\right|^{\gamma }\right] < \infty \text{ для некоторого } \gamma > 2 \text{. }

Если коэффициенты b=b(x)b = b(x) и σ=σ(x)\sigma = \sigma (x) являются непрерывными ограниченными функциями, то уравнение (115) имеет слабое решение.

Замечание 3. В предположении ограниченности и невырожденности коэффициента σ(x)\sigma (x) утверждение о существовании и единственности (по распределению) слабого решения уравнения (115) остается в силе, если потребовать лишь ограниченность и измеримость коэффициента b(x)b(x). Приведенные результаты о слабых решениях также допускают обобщения на многомерный случай, на случай коэффициентов, зависящих от прошлого и т. д. (см., например, [86; т. 1, гл. III]).

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

Пусть (Ω,F,(Ft)t⩾0,P)\left(\Omega , \mathscr {F},\left(\mathscr {F}_{t}\right)_{t \geqslant 0}, \mathbb {P}\right) — некоторое фильтрованное вероятностное пространство, W=(Wt,Ft)t⩾0W = \left(W_{t}, \mathscr {F}_{t}\right)_{t \geqslant 0} — mm-мерное броуновское движение, W=(W1,…,Wm)W = \left(W^{1}, \ldots , W^{m}\right). Пусть b=(bt,Ft)t⩾0b = \left(b_{t}, \mathscr {F}_{t}\right)_{t \geqslant 0}, где b=(b1,…,bm)b = \left(b^{1}, \ldots , b^{m}\right) — прогрессивно измеримый mm-мерный случайный процесс такой, что

P(∫0t∥bs∥2ds<∞)=1,t∈[0,T] \mathbb {P}\left(\int _{0}^{t}\left\| b_{s}\right\| ^{2} d s < \infty \right) = 1, \quad t \in [0, T]

здесь ∥bs∥2=(bs1)2+⋯+(bsm)2\left\| b_{s}\right\|^{2} = \left(b_{s}^{1}\right)^{2}+\cdots +\left(b_{s}^{m}\right)^{2} и T<∞T < \infty.

Образуем процесс Z=(Zt,Ft)t∈[0,T]Z = \left(Z_{t}, \mathscr {F}_{t}\right)_{t \in [0, T]}, полагая

Zt=exp(∫0t(bs,dWs)−12∫0t∥bs∥2ds) Z_{t} = \mathrm{exp}\left(\int _{0}^{t}\left(b_{s}, d W_{s}\right)-\frac{1}{2} \int _{0}^{t}\left\| b_{s}\right\| ^{2} d s\right)

где (bs,dWs):=∑k=1mbskdWsk\left(b_{s}, d W_{s}\right): = \sum_{k = 1}^{m} b_{s}^{k} d W_{s}^{k}.

Лемма 10 (см., например, [86; т. 2, с. 326]). Если выполнено условие Новикова:

E[exp(12∫0t∥bs∥2ds)]<∞ \mathbb {E}\left[\mathrm{exp}\left(\frac{1}{2} \int _{0}^{t}\left\| b_{s}\right\| ^{2} d s\right)\right] < \infty

то E[ZT]=1\mathbb {E}\left[Z_{T}\right] = 1 и введенный в (117) процесс Z=(Zt,Ft)t∈[0,T]Z = \left(Z_{t}, \mathscr {F}_{t}\right)_{t \in [0, T]} будет равномерно интегрируемым мартингалом.

В силу положительности ZtZ_{t} (P\mathbb {P}-п.н.) и условия E[ZT]=1\mathbb {E}\left[Z_{T}\right] = 1 на (Ω,FT)\left(\Omega , \mathscr {F}_{T}\right) можнозадать вероятностную меру QT\mathbb {Q}_{T}, полагая

QT(A)=E[1AZT],A∈FT \mathbb {Q}_{T}\left(A\right) = \mathbb {E}\left[\mathbf{1}_{A} Z_{T}\right], \quad A \in \mathscr {F}_{T}

Теорема 18 (Гирсанов). Для введенных выше процессов WW и bb определим

Bt=Wt−∫0tbsds,t∈[0,T] B_{t} = W_{t}-\int _{0}^{t} b_{s} d s, \quad t \in [0, T]

Тогда B=(Bt,Ft)t∈[0,T]B = \left(B_{t}, \mathscr {F}_{t}\right)_{t \in [0, T]} является mm-мерным броуновским движением на фильтрованном пространстве (Ω,FT,(Ft)t∈[0,T],QT)\left(\Omega , \mathscr {F}_{T},\left(\mathscr {F}_{t}\right)_{t \in [0, T]}, \mathbb {Q}_{T}\right).

Эта теорема, как уже было отмечено, дает возможность строить слабые решения стохастических дифференциальных уравнений вида

dXt=b(Xt)dt+dWt,t∈[0,T] d X_{t} = b\left(X_{t}\right) d t+d W_{t}, \quad t \in [0, T]

а также более общих уравнений (см., например, [86; гл. VII, § 36], [148; гл. 5]).

Наряду с интегралом Ито в ряде задач полезно обращение к "симметризованным" интегралам Стратоновича (или фиска-Стратоновича):

∫0tf(s,ω)∘dWs(ω) \int _{0}^{t} f(s, \omega ) \circ d W_{s}(\omega )

где W={Ws,s⩾0}W = \left\{ W_{s}, s \geqslant 0\right\} — броуновское движение, а функция ff входит в определенный класс. Чтобы пояснить суть дела, скажем, что для некоторого класса функций ff интеграл (121) строится с помощью интегральных сумм вида

∑i=0N−1f(ti∗,ω)(Wti+1−Wti) \sum _{i = 0}^{N-1} f\left(t_{i}^{*}, \omega \right)\left(W_{t_{i+1}}-W_{t_{i}}\right)

где 0=t0<⋯<tN=tи ti∗=(ti+ti+1)/2,i=0,…,N−10 = t_{0} < \cdots < t_{N} = t_{\text{и }} t_{i}^{*} = \left(t_{i}+t_{i+1}\right) / 2, i = 0, \ldots , N-1.

Доводы в пользу трактовки стохастического дифференциального уравнения как интегрального уравнения на основе интеграла (121) можно прочитать в [169]. В частности, замена переменных в интеграле Стратоновича не приводит к появлению дополнительных членов второго порядка, возникающих в формуле Ито. Это используется при изучении стохастических дифференциальных уравнений на многообразиях (см. [128]). В то же время при изучении стохастических уравнений вида

Xt=X0+∫0tb(s,Xs)ds+∫0tσ(s,Xs)∘dWs X_{t} = X_{0}+\int _{0}^{t} b\left(s, X_{s}\right) d s+\int _{0}^{t} \sigma \left(s, X_{s}\right) \circ d W_{s}

требуются бо́льшие ограничения на коэффициенты, чем при изучении уравнений вида (60). Если функция σ(t,x)\sigma (t, x) дифференцируема по xx, то уравнение (122) равносильно следующему уравнению Ито:

Xt=X0+∫0tb(s,Xs)ds+12∫0tσx′(s,Xs)σ(s,Xs)ds+∫0tσ(s,Xs)dWs X_{t} = X_{0}+\int _{0}^{t} b\left(s, X_{s}\right) d s+\frac{1}{2} \int _{0}^{t} \sigma _{x}^{\prime }\left(s, X_{s}\right) \sigma \left(s, X_{s}\right) d s+\int _{0}^{t} \sigma \left(s, X_{s}\right) d W_{s}

В частности, если σ(s,x)\sigma (s, x) не зависит от xx, то обе интерпретации уравнения (60) совпадают.

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

Для дальнейшего ознакомления с различными аспектами теории стохастических дифференциальных уравнений можно обратиться, например, к [11,130,154][11,130,154]. Стохастические уравнения в частных производных изучаются в [196]. Квантовая диффузия рассматривается, например, в [88].

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

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

Пусть интересующий нас nn-мерный процесс Xt,t∈[0,T]X_{t}, t \in [0, T], — полезный, но "скрытый" сигнал — описывается стохастическим дифференциальным уравнением (системой) (59) с начальным условием X0=ZX_{0} = Z и пусть "наблюдаемый" qq-мерньй процесс Yt,t∈[0,T]Y_{t}, t \in [0, T], задается стохастическим уравнением

dYt=c(t,Xt)dt+γ(t,Xt)dBt d Y_{t} = c\left(t, X_{t}\right) d t+\gamma \left(t, X_{t}\right) d B_{t}

где c:Rn+1→Rq,γ:Rn+1→Rq×pc: \mathbb {R}^{n+1} \rightarrow \mathbb {R}^{q}, \gamma : \mathbb {R}^{n+1} \rightarrow \mathbb {R}^{q \times p}. Считается, что все рассматриваемые случайные процессы определены на некотором вероятностном пространстве (Ω,F,P)(\Omega , \mathscr {F}, \mathbb {P}), причем pp-мерное броуновское движение B={Bt,t⩾0}B = \left\{ B_{t}, t \geqslant 0\right\} не зависит от mm-мерного броуновского движения W={Wt,t⩾0}W = \left\{ W_{t}, t \geqslant 0\right\}, фигурирующего в (59).

Одна из возможных здесь задач состоит в том, чтобы по наблюдениям {Ys,s∈[0,t]}\left\{ Y_{s}, s \in [0, t]\right\} дать оценку величины XtX_{t} для t∈[0,T]t \in [0, T]. Обозначив Ht=σ{Ys,s∈[0,t]}\mathscr {H}_{t} = \sigma \left\{ Y_{s}, s \in [0, t]\right\}, будем искать оптимальную оценку X~t∈Ht∣B(Rn)\widetilde{X}_{t} \in \mathscr {H}_{t} \mid \mathscr {B}\left(\mathbb {R}^{n}\right) такую, что

E[∣Xt−X~t∣2]=inf⁡{E[∣Xt−U∣2]:U∈L2(Ω,Ht,P)}. \mathbb {E}\left[\left|X_{t}-\tilde{X}_{t}\right|^{2}\right] = \inf \left\{ \mathbb {E}\left[\left|X_{t}-U\right|^{2}\right]: U \in L^{2}\left(\Omega , \mathscr {H}_{t}, \mathbb {P}\right)\right\} .

Легко видеть, что X~t=E[Xt∣Ht]\tilde{X}_{t} = \mathbb {E}\left[X_{t} \mid \mathscr {H}_{t}\right]. Мы будем предполагать, что все рассматриваемые случайные величины обладают конечным вторым моментом (это относится к компонентам векторов), а также, что выполнены условия, обеспечивающие существование решений используемых стохастических уравнений.

Следующий классический результат (см., например, [169]) для линейных систем явился одним из первых результатов, дающих удобные для практической реализации формулы, описывающие динамику оптимального фильтра.

Теорема 19 (фильтр Кальмана-Бьюси). Пусть

dXt=F(t)Xtdt+C(t)dWt,F(t)∈Rn×n,C(t)∈Rn×m,dYt=G(t)Xtdt+D(t)dBt,G(t)∈Rq×n,D(t)∈Rq×p \begin{aligned} & d X_{t} = F(t) X_{t} d t+C(t) d W_{t}, \quad F(t) \in \mathbb {R}^{n \times n}, \quad C(t) \in \mathbb {R}^{n \times m}, \\ & d Y_{t} = G(t) X_{t} d t+D(t) d B_{t}, \quad G(t) \in \mathbb {R}^{q \times n}, \quad D(t) \in \mathbb {R}^{q \times p} \end{aligned}

Пусть матрица D(t)D(t)∗D(t) D(t)^{*} обратима для всех tt и норма матрицы (D(t)D(t)∗)−1\left(D(t) D(t)^{*}\right)^{-1} ограничена на каждом конечном промежутке изменения tt. Тогда решение X~={X~t,t∈[0,T]}\widetilde{X} = \left\{ \widetilde{X}_{t}, t \in [0, T]\right\} задачи (125) удовлетворяет стохастическому дифференциальному уравнению

dX~t=(F−SG∗(DD∗)−1G)X~tdt+SG∗(DD∗)−1dYt d \tilde{X}_{t} = \left(F-S G^{*}\left(D D^{*}\right)^{-1} G\right) \tilde{X}_{t} d t+S G^{*}\left(D D^{*}\right)^{-1} d Y_{t}

с начальным условием X~0=E[X0]\tilde{X}_{0} = \mathbb {E}\left[X_{0}\right], а матричная функция

S(t)=E[(Xt−X~t)(Xt−X~t)∗]∈Rn×n S(t) = \mathbb {E}\left[\left(X_{t}-\tilde{X}_{t}\right)\left(X_{t}-\tilde{X}_{t}\right)^{*}\right] \in \mathbb {R}^{n \times n}

подчиняется матричному дифференциальному уравнению Риккати

dSdt=FS+SF∗−SG∗(DD∗)−1GS+CC∗,S(0)=E[(X0−E[X0])(X0−E[X0])∗]. \begin{gathered} \frac{d S}{d t} = F S+S F^{*}-S G^{*}\left(D D^{*}\right)^{-1} G S+C C^{*}, \\ S(0) = \mathbb {E}\left[\left(X_{0}-\mathbb {E}\left[X_{0}\right]\right)\left(X_{0}-\mathbb {E}\left[X_{0}\right]\right)^{*}\right] . \end{gathered}

Пример 3. Применив теорему 19 к определенному ниже частному случаю скалярных процессов, найдем X~t\tilde{X}_{t}. Пусть

dXt=0, т. е. Xt=X0,E[X0]=0,E[X02]=vdYt=Xtdt+mdBt,Y0=0 \begin{gathered} d X_{t} = 0, \quad \text{ т. е. } X_{t} = X_{0}, \quad \mathbb {E}\left[X_{0}\right] = 0, \quad \mathbb {E}\left[X_{0}^{2}\right] = v \\ d Y_{t} = X_{t} d t+m d B_{t}, \quad Y_{0} = 0 \end{gathered}

Уравнение Риккати (127) для S(t)=E[(Xt−X~t)2]S(t) = \mathbb {E}\left[\left(X_{t}-\tilde{X}_{t}\right)^{2}\right] приобретает вид

откуда

dSdt=−1m2S2,S(0)=v \frac{d S}{d t} = -\frac{1}{m^{2}} S^{2}, \quad S(0) = v S(t)=vm2m2+vt,t⩾0 S(t) = \frac{v m^{2}}{m^{2}+v t}, \quad t \geqslant 0

Для X~={X~t}\tilde{X} = \left\{ \tilde{X}_{t}\right\} из (126) находим, что

dX~t=−vm2+vtX~tdt+vm2+vtdYt,X~0=E[X0]=0 d \tilde{X}_{t} = -\frac{v}{m^{2}+v t} \tilde{X}_{t} d t+\frac{v}{m^{2}+v t} d Y_{t}, \quad \tilde{X}_{0} = \mathbb {E}\left[X_{0}\right] = 0

Отсюда, используя формулу Ито для диффузионных процессов, находим, что

d(X~texp(∫vm2+vtdt))=exp(∫vm2+vtdt)vm2+vtdYt d\left(\tilde{X}_{t} \mathrm{exp}\left(\int \frac{v}{m^{2}+v t} d t\right)\right) = \mathrm{exp}\left(\int \frac{v}{m^{2}+v t} d t\right) \frac{v}{m^{2}+v t} d Y_{t}

следовательно,

X~t=vm2+vtYt,t⩾0 \tilde{X}_{t} = \frac{v}{m^{2}+v t} Y_{t}, \quad t \geqslant 0

Теории фильтрации посвящена, например, книга [160]. В качестве введения в данную область рекомендуем прочитать главу IV в [169].

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

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

Здесь для моделей с непрерывным временем вместо разностных уравнений возникают стохастические дифференциальные уравнения. При этом ценой значительных технических усложнений удается получить аналоги результатов, известных для дискретных моделей. Мы не приводим переформулировок "дискретных" определений, данных в дополнении к главе IV, на случай непрерывного времени (см. [86]). Заметим, что широкую известность получили исследования в рамках модели Блэка-Мертона-Шоулса, используюшие для описания эволюции цен геометрическое броуновское движение. Работы М. Шоулса и Р. Мертона были отмечены в 1997 г. Нобелевской премией по экономике ( Φ\Phi. Блэк скончался, а посмертно Нобелевские премии не присуждаются).

Итак, если теперь (B,S)(B, S)-рынок описывается стохастическими дифференциальными уравнениями

dBt=rBtdt,dSt=St(μdt+σdWt) d B_{t} = r B_{t} d t, \quad d S_{t} = S_{t}\left(\mu d t+\sigma d W_{t}\right)

где B=(Bt)t⩾0B = \left(B_{t}\right)_{t \geqslant 0} — банковский счет, а S=(St)t⩾0S = \left(S_{t}\right)_{t \geqslant 0} — цены акций, то имеет место следующий замечательный результат.

Теорема 20 (формула Блэка-Шоулса). В модели (128) с постоянными r,μ,σr, \mu , \sigma справедливая стоимость стандартного опциона-колл европейского типа с платежной функцией fT=(ST−K)+f_{T} = \left(S_{T}-K\right)^{+}определяется формулой

C(T)=S0Φ(y+)−Ke−rΦ(y−) C(T) = S_{0} \Phi \left(y_{+}\right)-K e^{-r} \Phi \left(y_{-}\right)

где

y±=log⁡(S0/K)+T(r±σ2/2)σT y_{ \pm } = \frac{\log \left(S_{0} / K\right)+T\left(r \pm \sigma ^{2} / 2\right)}{\sigma \sqrt{T}}

а Φ\Phi — функция распределения стандартной нормальной величины. В частности, при S0=KS_{0} = K и r=0r = 0 имеем

C(T)=S0[Φ(σT2)−Φ(−σT2)] C(T) = S_{0}\left[\Phi \left(\frac{\sigma \sqrt{T}}{2}\right)-\Phi \left(-\frac{\sigma \sqrt{T}}{2}\right)\right]

при этом C(T)∼KσT/(2π)C(T) \sim K \sigma \sqrt{T /(2 \pi )}, когда T→0T \rightarrow 0.

Имеются разные подходы к доказательству этой теоремы (см., например, [86]). Отметим, что наряду с мартингальным подходом, опирающимся на теорему Гирсанова, имеется подход (который и был использован Блэком и Шоулсом), состояший в получении при определенных условиях для цены совершенного хеджирования

Y(t,T)=C[t,T]=inf⁡{x>0:∃π с Xtπ=x и XTπ=fT P-п.н. } Y(t, T) = C_{[t, T]} = \inf \left\{ x > 0: \exists \pi \text{ с } X_{t}^{\pi } = x \text{ и } X_{T}^{\pi } = f_{T} \text{ $\mathbb {P}$-п.н. }\right\}

так называемого фундаментального уравнения

∂Y∂t+rS∂Y∂S+12σ2S2∂2Y∂S2=rY \frac{\partial Y}{\partial t}+r S \frac{\partial Y}{\partial S}+\frac{1}{2} \sigma ^{2} S^{2} \frac{\partial ^{2} Y}{\partial S^{2}} = r Y

с краевым условием Y(T,S)=(S−K)+Y(T, S) = (S-K)^{+}. При этом величина C[0,T]=C(T)C_{[0, T]} = C(T) (см. (129)).

Полное представление о проблемах стохастической финансовой математики можно получить, ознакомившись с [86,147,155][86,147,155].

?