確率微分方程式

同義語:stochastic differential equationSDE

概要

確率微分方程式(stochastic differential equation)とは、Brown 運動 $B_t$ について $dX_t=b(t,X_t)\,dt+\sigma(t,X_t)\,dB_t$、$X_0=x$ と書かれる方程式で、積分方程式 $X_t=x+\int_0^tb(s,X_s)\,ds+\int_0^t\sigma(s,X_s)\,dB_s$ の略記である(第 2 項は伊藤積分)。係数が連続で $x$ について一様に Lipschitz 条件を満たせば、与えられた Brown 運動に適合した連続な解がただ 1 つ存在し、道ごとに一意である。証明は常微分方程式と同じ逐次近似だが、伊藤積分の等長性と Doob の極大不等式で 2 乗平均を評価する。Ornstein–Uhlenbeck 過程や幾何 Brown 運動が代表例である。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 伊藤の公式, Brown運動, マルチンゲール, Picard–Lindelöfの定理

確率微分方程式(stochastic differential equation, SDE)は、Brown 運動 $B_t$ による「ゆらぎ」を含む微分方程式
$$ dX_t=b(t,X_t)\,dt+\sigma(t,X_t)\,dB_t,\qquad X_0=x $$
である。Brown 運動の道は微分できないので、これは積分方程式
$$ X_t=x+\int_0^tb(s,X_s)\,ds+\int_0^t\sigma(s,X_s)\,dB_s $$
の略記として読む。右辺の第 2 項は 伊藤の公式 の記事で定めた伊藤積分である。$b$ を ドリフト係数、$\sigma$ を 拡散係数 という。$\sigma=0$ なら道ごとの常微分方程式 $x'(t)=b(t,x(t))$ に戻り、$b=0$、$\sigma=1$ なら解は $x+B_t$ である。粒子の拡散、株価のモデル、雑音を受ける力学系など、時間とともにでたらめに揺れながら変化する量を表すのに使われる。
この記事の主定理は、係数が Lipschitz 条件を満たせば解がただ 1 つ存在するという定理である。常微分方程式の Picard–Lindelöfの定理 と同じく逐次近似で解を作るが、収束の評価には伊藤積分の等長性と Doob の極大不等式を使う。この理論は伊藤清(Kiyosi Itô)が確率積分とともに創始した。

設定

以下、$T>0$ を固定し、$(B_t)_{0\le t\le T}$ を確率空間 $(\Omega,\mathcal F,P)$ の上の標準 Brown 運動(Brown運動 の定義「Brown 運動」)で、道はすべての $\omega$ で連続であるとする。$\mathcal N$ を確率 $0$ の事象の全体とし、$\mathcal F_t$ を $B_s$($0\le s\le t$)たちと $\mathcal N$ が生成する $\sigma$-代数とする。
伊藤の公式 の記事では $\mathcal F_t$ を $B_s$ たちだけで生成していたが、確率 $0$ の事象を加えても補題「未来の増分の独立性」はそのまま成り立つ。実際、$\mathcal F_t$ の各事象 $A$ について、$B_s$($s\le t$)たちが生成する $\sigma$-代数の事象 $A'$ で $P(A\triangle A')=0$ となるものがある(そのような $A$ の全体は $\sigma$-代数で、生成元を含む)。よって $B_u-B_t$($u\ge t$)は $\mathcal F_t$ と独立である。同じ記事の伊藤積分の構成と等長性の証明はこの独立性と適合性しか使わないので、同じ記事の定義「適合過程と $L^2$ 連続性」の過程の全体 $\mathcal L_T$ と伊藤積分は、この $\mathcal F_t$ についてもそのまま使える。確率 $0$ の事象を加えておくと、確率 $1$ でしか定まらない極限を確率 $0$ の事象の上で取り替えても適合性が崩れない。
伊藤積分 $\int_0^TH_t\,dB_t$ は $L^2(\Omega)$ の元で、$H$ について線形であり(左端の Riemann 和が線形なので)、伊藤の公式 の定理「伊藤積分の存在と等長性」により期待値 $0$、$E\bigl[(\int_0^TH_t\,dB_t)^2\bigr]=\int_0^TE[H_t^2]\,dt$ を満たす。

伊藤積分の連続な版

解を道ごとの関数として扱うには、各時刻 $t$ の伊藤積分 $\int_0^tH_s\,dB_s$ を、$t$ について連続な 1 つの過程にまとめる必要がある。

伊藤積分の連続な版と極大不等式

$H\in\mathcal L_T$ とする。適合した過程 $(I_t)_{0\le t\le T}$ で、すべての道が連続であり、各 $t$ について確率 $1$ で $I_t=\int_0^tH_s\,dB_s$ となるものがある。このような過程は、確率 $1$ で道全体が一致する意味でただ 1 つである。さらに次が成り立つ。

  1. $\displaystyle E\Bigl[\sup_{0\le t\le T}I_t^2\Bigr]\le4\int_0^TE[H_t^2]\,dt$。
  2. $0\le s\le t\le T$ なら $E[(I_t-I_s)^2]=\int_s^tE[H_r^2]\,dr$。とくに $t\mapsto I_t$ は $L^2(\Omega)$ への連続写像である。

この $I$ も $\int_0^tH_s\,dB_s$ と書く。1 は マルチンゲール の系「$L^p$ での極大不等式」の連続時間版にあたる。

左端の和をマルチンゲールとみて極大不等式を当てる

要点:分割 $\Pi$ ごとに、左端の和を時刻 $t$ で打ち切った連続な過程 $I^\Pi$ を作る。2 つの分割の差は離散時間のマルチンゲールなので、Doob の極大不等式で道全体の差の上限を抑えられる。分割を十分速く細かくすると、$I^\Pi$ は確率 $1$ で一様に収束する。

詳しい証明を開く

打ち切った和。 分割 $\Pi:0=t_0< \cdots< t_k=T$ について $I^\Pi_t:=\sum_{i=1}^kH_{t_{i-1}}(B_{t_i\wedge t}-B_{t_{i-1}\wedge t})$ とおく。道は $t$ について連続で、$I^\Pi_t$ は $\mathcal F_t$-可測である。$t\in(t_{j-1},t_j]$ なら $I^\Pi_t$ は $[0,t]$ の分割 $\{t_0,\dots,t_{j-1},t\}$(幅は $\lvert\Pi\rvert$ 以下)の左端の和なので、$\lvert\Pi\rvert\to0$ で $I^\Pi_t\to\int_0^tH_s\,dB_s$(2 乗平均)である。

差の評価。 $\omega_H(\eta):=\sup\{\|H_s-H_r\|_{L^2}:\lvert s-r\rvert\le\eta\}$ とおく(伊藤の公式 の定理「伊藤積分の存在と等長性」の証明と同じ記号で、$\eta\to0$ で $0$ に近づく)。2 つの分割 $\Pi,\Pi'$ と有限集合 $D\subset[0,T]$ をとり、これらの点をすべて合わせた分割を $0=s_0< \cdots< s_m=T$ とする。小区間 $(s_{j-1},s_j]$ を含む $\Pi$・$\Pi'$ の小区間の左端を $a_j,a_j'$ とし、$\xi_j:=H_{a_j}-H_{a_j'}$、$\Delta_j:=B_{s_j}-B_{s_{j-1}}$、$N_l:=\sum_{j\le l}\xi_j\Delta_j$ とおくと、$I^\Pi_{s_l}-I^{\Pi'}_{s_l}=N_l$ である。$\xi_j$ は $\mathcal F_{s_{j-1}}$-可測で $L^2$ に属し、$\Delta_j$ は $\mathcal F_{s_{j-1}}$ と独立なので、条件付き期待値 の定理「条件付き期待値の基本性質」の 6・7 から $E[\xi_j\Delta_j\mid\mathcal F_{s_{j-1}}]=\xi_jE[\Delta_j]=0$ であり、$(N_l)$ は $(\mathcal F_{s_l})$ に関するマルチンゲールである。マルチンゲール の系「$L^p$ での極大不等式」($p=2$)と 伊藤の公式 の補題「左端の和の等長性」から

$$E\Bigl[\max_{s\in D}\lvert I^\Pi_s-I^{\Pi'}_s\rvert^2\Bigr]\le E\Bigl[\max_lN_l^2\Bigr]\le4E[N_m^2]=4\sum_jE[\xi_j^2](s_j-s_{j-1})\le4T\,\omega_H(\eta)^2$$

である($\eta:=\max(\lvert\Pi\rvert,\lvert\Pi'\rvert)$)。有限集合 $D$ を増やして $[0,T]$ の稠密な可算集合 $\{T\}\cup\{j2^{-n}T\}$ に近づけると、道の連続性から最大値は $[0,T]$ 全体での上限に増加して近づくので、単調収束定理 により $E[\sup_t\lvert I^\Pi_t-I^{\Pi'}_t\rvert^2]\le4T\omega_H(\eta)^2$ である。同じ議論を $\Pi'$ の代わりに「$H\equiv0$ の和」として $I^\Pi$ 自身に当てると $E[\sup_t(I^\Pi_t)^2]\le4\sum_iE[H_{t_{i-1}}^2](t_i-t_{i-1})$ を得る。

一様収束。 $\lvert\Pi_n\rvert$ が減少して $0$ に近づき $4T\omega_H(\lvert\Pi_n\rvert)^2\le8^{-n}$ となる分割の列をとる。Chebyshev の不等式から $P(\sup_t\lvert I^{\Pi_{n+1}}_t-I^{\Pi_n}_t\rvert>2^{-n})\le4^n\cdot8^{-n}=2^{-n}$ であり、Borel–Cantelliの補題 により、確率 $1$ である番号から先は上限が $2^{-n}$ 以下である。この事象の上で $I^{\Pi_n}$ は一様に Cauchy 列なので、連続関数に一様収束する。その極限を $I$ とし、収束しない確率 $0$ の事象の上では $I\equiv0$ とおく。$\mathcal N\subset\mathcal F_t$ なので $I_t$ は $\mathcal F_t$-可測である。各 $t$ で $I^{\Pi_n}_t$ は $\int_0^tH_s\,dB_s$ に 2 乗平均で、$I_t$ に確率 $1$ で収束するので、両者は確率 $1$ で等しい。

一意性と 1・2。 連続な道をもつ 2 つの過程が各有理数時刻(と $T$)で確率 $1$ で一致すれば、可算個の時刻での一致と連続性から、確率 $1$ で道全体が一致する。1 は Fatouの補題 と上の評価から $E[\sup_tI_t^2]\le\liminf_nE[\sup_t(I^{\Pi_n}_t)^2]\le4\int_0^TE[H_t^2]\,dt$(右辺の和は連続関数 $t\mapsto E[H_t^2]$ の Riemann 和)。2 は、$s$ を分点に含む分割では $[0,t]$ の左端の和が $[0,s]$ の和と $[s,t]$ の和に分かれ、後者は 伊藤の公式 の補題「左端の和の等長性」により 2 乗平均が $\sum E[H^2]\delta$ になることから、極限をとって得る。$\square$

伊藤積分は $H$ について線形なので、$H,H'\in\mathcal L_T$ の連続な版について、確率 $1$ で道全体として $\int_0^t(H_s-H'_s)\,dB_s=\int_0^tH_s\,dB_s-\int_0^tH'_s\,dB_s$ である。

解の定義

以下、係数 $b,\sigma\colon[0,T]\times\mathbb R\to\mathbb R$ は連続関数で、次の Lipschitz 条件(以下 (L) と書く)を満たすとする:定数 $L$ があって
$$ \lvert b(t,x)-b(t,y)\rvert+\lvert\sigma(t,x)-\sigma(t,y)\rvert\le L\lvert x-y\rvert\qquad(t\in[0,T],\ x,y\in\mathbb R). $$
このとき $K:=L+\max_{0\le t\le T}(\lvert b(t,0)\rvert+\lvert\sigma(t,0)\rvert)$ とおくと、$\lvert b(t,x)\rvert+\lvert\sigma(t,x)\rvert\le K(1+\lvert x\rvert)$ である(線形増大条件)。

確率微分方程式の解

$x\in\mathbb R$ とする。確率過程 $(X_t)_{0\le t\le T}$ が次を満たすとき、$X$ を確率微分方程式 $dX_t=b(t,X_t)\,dt+\sigma(t,X_t)\,dB_t$、$X_0=x$ の 解(強い解)という。

  1. $X$ は $(\mathcal F_t)$ に適合し、すべての道 $t\mapsto X_t(\omega)$ は連続である。
  2. 各 $t$ で $E[X_t^2]< \infty$ であり、$t\mapsto X_t$ は $[0,T]$ から $L^2(\Omega)$ への連続写像である。
  3. 確率 $1$ で、すべての $t\in[0,T]$ について
    $$ X_t=x+\int_0^tb(s,X_s)\,ds+\int_0^t\sigma(s,X_s)\,dB_s $$
    が成り立つ。第 1 項は道ごとの Riemann 積分、第 2 項は thm-sde-continuous の連続な版である。

「強い」は、与えられた Brown 運動 $B$ とその情報 $\mathcal F_t$ のもとで解を作ることを指す。3 で第 2 項が意味をもつことは次の補題が保証する。

係数を代入した過程

$X$ が def-sde-solution の 1・2 を満たすとする。このとき $(b(t,X_t))_t$ と $(\sigma(t,X_t))_t$ は $\mathcal L_T$ に属する。さらに
$$ \Phi(X)_t:=x+\int_0^tb(s,X_s)\,ds+\int_0^t\sigma(s,X_s)\,dB_s $$
(第 2 項は連続な版)とおくと、$\Phi(X)$ も 1・2 を満たし、$E[\sup_t\Phi(X)_t^2]< \infty$ である。

要点:連続性・適合性は合成と Riemann 和の極限から、$L^2$ 連続性は Lipschitz 条件と優収束定理から従う。

詳しい証明を開く

$\sigma$ は連続なので $\sigma(t,X_t)$ は $\mathcal F_t$-可測で、道は連続である。$L^2$ 連続性は 2 項に分けて確かめる。$\|\sigma(t,X_t)-\sigma(s,X_s)\|_{L^2}\le\|\sigma(t,X_t)-\sigma(s,X_t)\|_{L^2}+L\|X_t-X_s\|_{L^2}$ の第 2 項は 2 により $s\to t$ で $0$ に近づく。第 1 項の被積分関数は $\sigma$ の連続性により各点で $0$ に近づき、$4K^2(1+\lvert X_t\rvert)^2$ で抑えられるので、優収束定理 により $0$ に近づく。$E[\sigma(t,X_t)^2]\le2K^2(1+E[X_t^2])< \infty$ である。$b$ についても同じである。

$\Phi(X)$ の第 1 項は道ごとに $t$ の連続関数で、$\mathcal F_t$-可測な値の Riemann 和の極限なので $\mathcal F_t$-可測である。$t\mapsto E[X_t^2]$ は連続なので有界であり、$M:=\sup_tE[b(t,X_t)^2]< \infty$ である。Cauchy–Schwarz の不等式から $E[(\int_s^tb(r,X_r)\,dr)^2]\le(t-s)\int_s^tE[b(r,X_r)^2]\,dr\le(t-s)^2M$ なので第 1 項は $L^2$ 連続であり、第 2 項は thm-sde-continuous の 2 により $L^2$ 連続である。最後に $\sup_t\Phi(X)_t^2\le3\bigl(x^2+T\int_0^Tb(s,X_s)^2ds+\sup_tI_t^2\bigr)$($I$ は第 2 項)であり、thm-sde-continuous の 1 から期待値は有限である。$\square$

def-sde-solution の 3 は、$X$ と $\Phi(X)$ の道が確率 $1$ で一致することにほかならない。両方とも道が連続なので、各 $t$ ごとに確率 $1$ で $X_t=\Phi(X)_t$ となることと同値である。

主定理:解の存在と一意性

2 つの過程の差の評価

$X,Y$ が def-sde-solution の 1・2 を満たすとし、$C:=2L^2(T+4)$ とおく。$0\le t\le T$ について
$$ E\Bigl[\sup_{0\le r\le t}\lvert\Phi(X)_r-\Phi(Y)_r\rvert^2\Bigr]\le C\int_0^tE\bigl[\lvert X_s-Y_s\rvert^2\bigr]\,ds $$
が成り立つ。

ドリフトは Cauchy–Schwarz、拡散は極大不等式で抑える

伊藤積分の線形性から、確率 $1$ ですべての $r$ について
$$ \Phi(X)_r-\Phi(Y)_r=\int_0^r\bigl(b(s,X_s)-b(s,Y_s)\bigr)\,ds+\int_0^r\bigl(\sigma(s,X_s)-\sigma(s,Y_s)\bigr)\,dB_s $$
である。$r\le t$ のとき、第 1 項の 2 乗は Cauchy–Schwarz の不等式と (L) により $r\int_0^rL^2\lvert X_s-Y_s\rvert^2ds\le TL^2\int_0^t\lvert X_s-Y_s\rvert^2ds$ 以下である。第 2 項の $r\le t$ での上限の 2 乗の期待値は、thm-sde-continuous の 1 を区間 $[0,t]$ で使い (L) を当てて、$4L^2\int_0^tE[\lvert X_s-Y_s\rvert^2]\,ds$ 以下である。$(u+v)^2\le2u^2+2v^2$ と Fubiniの定理(非負関数の積分の交換)により結論を得る。$\square$

確率微分方程式の解の存在と一意性

$b,\sigma\colon[0,T]\times\mathbb R\to\mathbb R$ を連続で (L) を満たす関数とし、$x\in\mathbb R$ とする。

  1. 存在:$dX_t=b(t,X_t)\,dt+\sigma(t,X_t)\,dB_t$、$X_0=x$ の解 $X$ がある。
  2. 一意性:$X,Y$ がともに解なら、確率 $1$ ですべての $t\in[0,T]$ について $X_t=Y_t$ である。
  3. 解は $E[\sup_{0\le t\le T}X_t^2]< \infty$ を満たす。
Picard の逐次近似を 2 乗平均で評価する

要点:$X^0\equiv x$、$X^{n+1}:=\Phi(X^n)$ とおくと、lem-sde-estimate を繰り返し使って、$n$ 番目の差の 2 乗平均が $(Ct)^n/n!$ の程度に小さくなる。これは Picard–Lindelöfの定理 の命題「逐次近似の誤差」と同じ型の評価であり、$\sum_n(Ct)^n/n!< \infty$ から道ごとの一様収束が従う。一意性も同じ評価で示す。

詳しい証明を開く

逐次近似。 定数過程 $X^0\equiv x$ は def-sde-solution の 1・2 を満たすので、lem-sde-integrand により $X^{n+1}:=\Phi(X^n)$ が帰納的に定まり、どれも 1・2 を満たす。$d_n(t):=E[\sup_{r\le t}\lvert X^{n+1}_r-X^n_r\rvert^2]$ とおく。lem-sde-integrand から $A:=d_0(T)\le2E[\sup_r(X^1_r)^2]+2x^2< \infty$ である。lem-sde-estimate により $d_n(t)\le C\int_0^tE[\lvert X^n_s-X^{n-1}_s\rvert^2]\,ds\le C\int_0^td_{n-1}(s)\,ds$ なので、帰納法で $d_n(t)\le A(Ct)^n/n!$ である。

収束。 $\varepsilon_n:=(A(CT)^n/n!)^{1/2}$ とおくと $\sum_n\varepsilon_n< \infty$ である。Chebyshev の不等式から $P(\sup_r\lvert X^{n+1}_r-X^n_r\rvert>\varepsilon_n^{1/2})\le\varepsilon_n$ であり、Borel–Cantelliの補題 により、確率 $1$ である番号から先は $\sup_r\lvert X^{n+1}_r-X^n_r\rvert\le\varepsilon_n^{1/2}$ である。$\sum_n\varepsilon_n^{1/2}< \infty$(比の判定法)なので、この事象の上で $X^n$ は連続関数に一様収束する。その極限を $X$ とし、残りの確率 $0$ の事象の上では $X\equiv x$ とおくと、$X$ は適合し道は連続である。Minkowski の不等式から $m>n$ について $\|\sup_r\lvert X^m_r-X^n_r\rvert\|_{L^2}\le\sum_{j\ge n}\varepsilon_j$ であり、$m\to\infty$ として Fatouの補題 を使うと $E[\sup_r\lvert X_r-X^n_r\rvert^2]\le(\sum_{j\ge n}\varepsilon_j)^2\to0$ である。よって各 $X_t$ は $L^2$ に属し、$t\mapsto X_t$ は $L^2$ 連続な写像 $t\mapsto X^n_t$ の一様な極限なので $L^2$ 連続である。$X$ は 1・2 を満たす。

極限が解であること。 lem-sde-estimate から $E[\sup_r\lvert\Phi(X)_r-X^{n+1}_r\rvert^2]=E[\sup_r\lvert\Phi(X)_r-\Phi(X^n)_r\rvert^2]\le CT\,E[\sup_r\lvert X_r-X^n_r\rvert^2]\to0$ であり、$X^{n+1}\to X$ も同じ意味で成り立つので、$E[\sup_r\lvert\Phi(X)_r-X_r\rvert^2]=0$ である。よって確率 $1$ で道全体として $X=\Phi(X)$ であり、$X$ は解である。

一意性と 3。 $X,Y$ を解とする。$X=\Phi(X)$、$Y=\Phi(Y)$ なので、lem-sde-integrand から $E[\sup_tX_t^2]< \infty$(3)であり、$D(t):=E[\sup_{r\le t}\lvert X_r-Y_r\rvert^2]$ は有限である。lem-sde-estimate から $D(t)\le C\int_0^tD(s)\,ds$ なので、繰り返し代入して $D(t)\le D(T)(Ct)^n/n!$ がすべての $n$ で成り立ち、$D(T)=0$ である。$\square$

2 の一意性は、同じ Brown 運動に対する 2 つの解の道が一致するという強い意味のもので、道ごとの一意性(pathwise uniqueness)と呼ばれる。証明の評価 $d_n(t)\le A(Ct)^n/n!$ は、常微分方程式の逐次近似の評価と同じ形をしている。違いは、道ごとの差の上限ではなく、その 2 乗の期待値を評価することにある。伊藤積分は道ごとには定まらず 2 乗平均でしか制御できないためである。

例

確定的な被積分関数の部分積分

$0< t\le T$ とし、$g\colon[0,t]\to\mathbb R$ を $C^1$ 級の関数とする。確率 $1$ で
$$ \int_0^tg(s)\,dB_s=g(t)B_t-\int_0^tg'(s)B_s\,ds $$
が成り立つ。左辺は $N\bigl(0,\int_0^tg(s)^2ds\bigr)$ に従う。

和の順序の入れ替え

$[0,t]$ の分割 $0=t_0< \cdots< t_k=t$ について、$B_0=0$ から
$$ \sum_{i=1}^kg(t_{i-1})(B_{t_i}-B_{t_{i-1}})=g(t_{k-1})B_t-\sum_{i=1}^{k-1}\bigl(g(t_i)-g(t_{i-1})\bigr)B_{t_i} $$
である(和の順序を入れ替えた Abel の変形)。平均値の定理から $g(t_i)-g(t_{i-1})=g'(\eta_i)(t_i-t_{i-1})$ となる $\eta_i$ があり、右辺の和は道ごとに連続関数 $g'(s)B_s$ の Riemann 和に近く(差は $g'$ と道の一様連続性で $0$ に近づく)、$\lvert\Pi\rvert\to0$ で $\int_0^tg'(s)B_s\,ds$ に収束する。左辺は 伊藤の公式 の定理「伊藤積分の存在と等長性」により伊藤積分に 2 乗平均で収束するので、部分列をとれば確率 $1$ で収束し、両者の極限は確率 $1$ で等しい。
分布について:左辺の和は独立な正規変数 $g(t_{i-1})(B_{t_i}-B_{t_{i-1}})$ の和なので、正規分布 の定理「独立な正規分布の和」により $N(0,v_\Pi)$、$v_\Pi:=\sum_ig(t_{i-1})^2(t_i-t_{i-1})$ に従う。和を $Y_\Pi$、極限を $Y$ とすると $\lvert E[e^{i\lambda Y_\Pi}]-E[e^{i\lambda Y}]\rvert\le\lvert\lambda\rvert\,E\lvert Y_\Pi-Y\rvert\to0$ であり、$v_\Pi\to v:=\int_0^tg^2ds$ なので $E[e^{i\lambda Y}]=e^{-\lambda^2v/2}$ である。正規分布 の定理「積率母関数と特性関数」と 特性関数(確率論) の系「特性関数は分布を決める」により $Y\sim N(0,v)$ である。$\square$

Ornstein–Uhlenbeck 過程

$a>0$、$\sigma_0>0$ を定数とし、$dX_t=-aX_t\,dt+\sigma_0\,dB_t$、$X_0=x$ を考える。係数 $b(t,x)=-ax$、$\sigma(t,x)=\sigma_0$ は (L) を満たすので、thm-sde-main により解はただ 1 つである。その解は
$$ X_t=e^{-at}x+\sigma_0\int_0^te^{-a(t-s)}\,dB_s $$
であり、$X_t\sim N\Bigl(e^{-at}x,\ \dfrac{\sigma_0^2}{2a}(1-e^{-2at})\Bigr)$ である。
確かめる。拡散係数が定数なので $\int_0^t\sigma_0\,dB_s=\sigma_0B_t$ であり、$Y_t:=X_t-\sigma_0B_t$ とおくと、方程式は道ごとの線形常微分方程式 $Y_t'=-a(Y_t+\sigma_0B_t)$、$Y_0=x$ になる。定数変化法で解くと
$$ X_t=e^{-at}x+\sigma_0B_t-a\sigma_0\int_0^te^{-a(t-s)}B_s\,ds $$
で、これは道が連続で適合し、方程式を道ごとに満たす。$\lvert X_t\rvert\le\lvert x\rvert+(1+aT)\sigma_0\sup_{s\le T}\lvert B_s\rvert$ で、右辺は thm-sde-continuous の 1($H\equiv1$)により 2 乗可積分なので、優収束定理 から $L^2$ 連続性も従う。lem-sde-parts を $g(s)=e^{as}$ に当てて $e^{-at}$ を掛けると、右辺の後ろ 2 項は $\sigma_0\int_0^te^{-a(t-s)}dB_s$ に等しく、分散は $\sigma_0^2\int_0^te^{-2a(t-s)}ds=\frac{\sigma_0^2}{2a}(1-e^{-2at})$ である。$t\to\infty$ で分布は $N(0,\sigma_0^2/2a)$ に近づく。ドリフト $-aX_t$ が値を $0$ に引き戻し、雑音と釣り合った分布に落ち着く。

幾何 Brown 運動

$\mu\in\mathbb R$、$\sigma_0>0$、$x>0$ とし、$dX_t=\mu X_t\,dt+\sigma_0X_t\,dB_t$、$X_0=x$ を考える。係数 $\mu x$、$\sigma_0x$ は (L) を満たすので解はただ 1 つで、それは
$$ X_t=x\exp\Bigl(\sigma_0B_t+\Bigl(\mu-\frac{\sigma_0^2}2\Bigr)t\Bigr) $$
である。実際、$f(y,s):=x\exp(\sigma_0y+(\mu-\sigma_0^2/2)s)$ は $\partial_yf=\sigma_0f$、$\partial_y^2f=\sigma_0^2f$、$\partial_sf=(\mu-\sigma_0^2/2)f$ を満たすので、時刻にもよる伊藤の公式(MP10 Theorem 7.14(p. 197、証明 pp. 198–199)で増加過程を $s$ とした場合。この記事では証明しない。同定理の仮定 $E\int_0^t\partial_yf(B_s,s)^2\,ds=\sigma_0^2\int_0^tE[X_s^2]\,ds< \infty$ は、$E[X_s^2]=x^2e^{(2\mu+\sigma_0^2)s}$ から従う)から $X_t-x=\int_0^t\sigma_0X_s\,dB_s+\int_0^t\mu X_s\,ds$ である(同書の確率積分は、この記事の被積分過程では左端の和の極限として 伊藤の公式 の伊藤積分と一致する)。$X_t$ は対数正規分布に従う。$s\le t$ なら $B_s+B_t=2B_s+(B_t-B_s)\sim N(0,3s+t)$ なので、正規分布 の定理「積率母関数と特性関数」から $E[X_sX_t]=x^2\exp\bigl(\frac{\sigma_0^2}2(3s+t)+(\mu-\frac{\sigma_0^2}2)(s+t)\bigr)$ であり、これは $s,t$ について連続で $s=t$ のとき $E[X_t^2]$ に等しい。よって $E[(X_t-X_s)^2]=E[X_t^2]-2E[X_sX_t]+E[X_s^2]\to0$($s\to t$)で、$X$ は def-sde-solution の 1・2 も満たす。
$\sigma_0$ を含む補正 $-\sigma_0^2t/2$ は、伊藤の公式 の 2 次の項から来る。この過程は株価のモデルとして使われ、$E[X_t]=xe^{\mu t}$ なので $\mu$ は平均の成長率を表す。

弱い解

thm-sde-main は、与えられた Brown 運動 $B$ から解を作る。これに対し、Brown 運動も含めて組を作ればよいとする解の概念がある。

弱い解

確率空間 $(\Omega',\mathcal F',Q)$、フィルトレーション $(\mathcal G_t)_{0\le t\le T}$、$(\mathcal G_t)$ に関する標準 Brown 運動 $W$($W$ は $(\mathcal G_t)$ に適合し、各 $t\le u$ で $W_u-W_t$ が $\mathcal G_t$ と独立)、適合した過程 $X$ の組で、$X$ が $W$ と $(\mathcal G_t)$ について def-sde-solution の 1–3 を満たすもの(3 の伊藤積分とその連続な版は、伊藤の公式 の伊藤積分と thm-sde-continuous を、$B$ と $\mathcal F_t$ の代わりに $W$ と「$\mathcal G_t$ に $Q$ で確率 $0$ の事象を加えた $\sigma$-代数」で考えたもの。設定の段落と同じ理由でそのまま使える)を、$dX_t=b(t,X_t)\,dt+\sigma(t,X_t)\,dW_t$、$X_0=x$ の 弱い解 という。

強い解は、$W=B$、$\mathcal G_t=\mathcal F_t$ とした弱い解である。弱い解では、$X$ が $W$ の過去だけから決まることを要求しない。係数が (L) を満たさなくても弱い解が作れる場合がある。たとえば、$b\colon\mathbb R\to\mathbb R$ が有界で連続なら、Girsanovの定理 の系「ドリフトを付けた方程式の弱い解」により、確率測度を取り替えて $dX_t=b(X_t)\,dt+dW_t$、$X_0=0$ の弱い解が得られる。この構成では、もとの Brown 運動 $B$ をそのまま解 $X$ とし、$W_t:=B_t-\int_0^tb(B_s)\,ds$ が新しい確率測度のもとで Brown 運動になる。逆に、弱い解が存在して道ごとの一意性が成り立てば、強い解が存在し、解の分布もただ 1 つに決まる(山田–渡辺の定理。解と一意性の正確な意味は同記事の定義による)。道ごとの一意性の仮定を外せないことは、同記事の例「田中の方程式」が示している。

条件を外したときの反例

拡散係数を $\sigma=0$ とすると、方程式は道ごとの常微分方程式になる。thm-sde-main の仮定の役割は、すでにこの場合に現れる。

外す条件反例成り立たなくなること
(L)($0$ の近くで Lipschitz でない)$b(t,x)=3\lvert x\rvert^{2/3}$、$\sigma=0$、$x=0$一意性:$X\equiv0$ と $X_t=t^3$ がともに解
(L) の大域性(各有界集合の上では Lipschitz)$b(t,x)=x^2$、$\sigma=0$、$x=1$、$T\ge1$存在:$[0,T]$ 上の解がない

1 行目:$X_t=t^3$ について $X_t'=3t^2=3\lvert X_t\rvert^{2/3}$ である。$b$ は連続だが、$x\to0$ で $\lvert b(t,x)-b(t,0)\rvert/\lvert x\rvert=3\lvert x\rvert^{-1/3}\to\infty$ なので (L) を満たさない。2 行目:$b$ は有界集合の上では Lipschitz で、Picard–Lindelöfの定理 の系「局所 Lipschitz 条件の下での存在と一意性」により、$[0,T]$ 上の解があれば $[0,1)$ の上で $1/(1-t)$ と一致する。これは $t\to1$ で発散するので、$[0,1]$ で連続な解はない。大域的な (L) から出る線形増大条件が、解が有限時間で無限大に飛ぶこと(爆発)を防いでいる。

補足

  • 多次元の場合:応用では、$X_t\in\mathbb R^d$、$b$ を $\mathbb R^d$ 値、$\sigma$ を $d\times m$ 行列値とし、$m$ 次元 Brown 運動(成分が独立な Brown 運動)で $dX_t=b\,dt+\sigma\,dB_t$ を考えることが多い。この記事では 1 次元の場合だけを扱う。
  • 偏微分方程式との関係:$b=0$、$\sigma=1$ のとき解は $x+B_t$ で、その期待値 $E[f(x+B_t)]$ が熱方程式を満たすことが 伊藤の公式 の系「期待値と熱方程式」である。
  • 名前:伊藤清は、ランダムな被積分過程をもつ確率積分を初めて定義し(1944 年。MP10 第 7 章の Notes and comments、p. 222)、それにもとづいて確率微分方程式の理論を築いた。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する