コンテンツにスキップ

多次元伊藤公式と Lévy の特徴付け:共変分・積の公式・Brown 運動の判定法

対象読者: 本テーマの伊藤積分(確率微分方程式)の章とマルチンゲール、ブラウン運動の章を学んだ人。多変数のテイラー展開と特性関数は既知とします。

前提となる記事:確率微分方程式:伊藤積分はなぜ「未来を見ない」のか、マルチンゲールとブラウン運動:公平な賭けから連続時間の確率過程へ

生 Markdown

複数の株価や、3 次元空間を揺れ動く粒子のように、ランダムな揺れが同時にいくつも働く状況を考えます。このとき「二つの量の積は、微小時間にどれだけ変化するか」「ある過程が、ランダムな歩みの基本形である Brown 運動かどうかを、分布を計算せずに見分けられるか」という素朴な疑問が出ます。この記事は、前者に答える多次元の伊藤公式と積の公式、後者に答える Lévy の判定条件を、証明つきで示します。読み終えると、複数成分の確率過程の微小変化を計算でき、与えられた過程が Brown 運動であることを確かめられるようになります。

  • dd 次元 Brown 運動 W=(W1,…,Wd)W=(W^1,\dots,W^d) の成分は互いに独立なので、異なる成分による伊藤積分は直交します。これを条件付きの形で述べたのが共分散等長性 E[ξ∫H dWa∫K dWb]=δabE[ξ∫HK du]\mathbb{E}[\xi\int H\,dW^a\int K\,dW^b]=\delta_{ab}\mathbb{E}[\xi\int HK\,du] で、多次元の計算規則 dWa dWb=δab dtdW^a\,dW^b=\delta_{ab}\,dt の正体です。
  • 伊藤過程 dXi=μi dt+∑aσia dWadX^i=\mu^i\,dt+\sum_a\sigma^{ia}\,dW^a の共変分を ⟨Xi,Xj⟩t=∫0t(σσT)ij du\langle X^i,X^j\rangle_t=\int_0^t(\sigma\sigma^{\mathsf T})^{ij}\,du と定めると、増分の積の和 ∑kΔXkiΔXkj\sum_k\Delta X^i_k\Delta X^j_k はこれに L1L^1 で収束します。証明の要は伊藤積分の 4 次モーメント評価です。
  • 多次元伊藤公式は、2 次の項 12∑i,j∂ijf d⟨Xi,Xj⟩\frac12\sum_{i,j}\partial_{ij}f\,d\langle X^i,X^j\rangle を持つ連鎖律です。特に f(x,y)=xyf(x,y)=xy から積の公式 d(XY)=X dY+Y dX+d⟨X,Y⟩d(XY)=X\,dY+Y\,dX+d\langle X,Y\rangle が出ます。
  • Lévy の特徴付け: 連続なマルチンゲール MM で MiMj−δijtM^iM^j-\delta_{ij}t もマルチンゲールなら、MM は Brown 運動です。証明は eiλxe^{i\lambda x} の 3 次までの Taylor 展開と停止による有界化だけで済みます。
  • 応用として、∫0tsgn⁡(Wu) dWu\int_0^t\operatorname{sgn}(W_u)\,dW_u や直交行列値の過程による回転 ∫O dW\int O\,dW が再び Brown 運動であることがすぐに分かります。

確率微分方程式の章 では、1 本の Brown 運動 BB に関する伊藤積分と伊藤の公式を扱いました。しかし応用で現れる確率過程の多くは多次元です。R3\mathbb{R}^3 の中を拡散する粒子、複数の株価、場の理論で多数のモードを同時に動かす確率量子化など、いずれも複数の独立な雑音で駆動されます。

多次元に移ると、新しい問いが二つ生じます。第一に、二つの成分の積 X1X2X^1X^2 はどう微分されるのか。1 次元の公式 (dB)2=dt(dB)^2=dt の多次元版は dWa dWb=δab dtdW^a\,dW^b=\delta_{ab}\,dt ですが、これは「異なる成分の増分の積の和は 00 に収束する」という主張で、証明が要ります。第二に、ある連続マルチンゲールが Brown 運動であることを、分布を直接計算せずに判定できないか。Lévy の特徴付けはこれに「2 次変分が tt なら Brown 運動」と答えます。その証明で使う「条件付き特性関数による独立性と分布の判定」は、次章の Girsanov の定理 の証明でもそのまま使います。また Lévy の特徴付けは、確率積分による時間変更の議論の土台になります。

この章では、確率微分方程式の章で構成した伊藤積分(定理 4.7[確率微分方程式])だけを道具にして、これらを証明します。一般の連続セミマルチンゲールの理論(Doob–Meyer 分解など)には立ち入らず、係数が有界な伊藤過程に話を絞ります。この範囲で、Girsanov の定理とマルチンゲール表現定理に必要なものはすべてそろいます。

時間区間を [0,T][0,T](T<∞T<\infty)とし、確率空間 (Ω,F,P)(\Omega,\mathcal{F},\mathbb{P}) とフィルトレーション (Ft)(\mathcal{F}_t) を固定します。発展的可測性などの用語は 定義 2.1[確率微分方程式] に従います。フィルトレーションとは時刻ごとに増えていく「その時点までに分かっている情報」の列で、発展的可測とは、各時刻 rr について、rr までの時間と rr までの情報だけで過程の値が決まる、という可測性のことです。

定義 2.1(d 次元 (Ft)(F_t)-Brown 運動)↗

Rd\mathbb{R}^d 値の過程 W=(W1,…,Wd)W=(W^1,\dots,W^d) が (Ft)(\mathcal{F}_t)-Brown 運動であるとは、(1) W0=0W_0=0 ですべての道が連続、(2) WW は (Ft)(\mathcal{F}_t) に適合、(3) s<ts<t のとき増分 Wt−WsW_t-W_s は Fs\mathcal{F}_s と独立で、平均 00、共分散行列 (t−s)Id(t-s)I_d の Gauss 分布 N(0,(t−s)Id)N(0,(t-s)I_d) に従う、の三条件を満たすことをいいます。

共分散行列が対角なので、増分の成分 Wt1−Ws1,…,Wtd−WsdW^1_t-W^1_s,\dots,W^d_t-W^d_s は互いに独立です(命題 5.3[特性関数と Gauss ベクトル])。特に各成分 WaW^a は 1 次元の (Ft)(\mathcal{F}_t)-Brown 運動で、確率微分方程式の章の伊藤積分 ∫H dWa\int H\,dW^a がそのまま使えます。∫stH dWa\int_s^tH\,dW^a は ∫0tH dWa−∫0sH dWa\int_0^tH\,dW^a-\int_0^sH\,dW^a の意味です。

定義 2.2(多次元伊藤過程と有界係数の条件)↗

L2[0,T]\mathcal{L}^2[0,T] を 定義 4.1[確率微分方程式] の空間とします。Rd\mathbb{R}^d 値の過程 H=(H1,…,Hd)H=(H^1,\dots,H^d) で各成分が L2[0,T]\mathcal{L}^2[0,T] に属するもの全体を Ld2\mathcal{L}^2_d と書き、

∫0tHu⋅dWu:=∑a=1d∫0tHua dWua\int_0^t H_u\cdot dW_u:=\sum_{a=1}^d\int_0^t H^a_u\,dW^a_u

とおきます。Rn\mathbb{R}^n 値の過程 XX が

Xti=X0i+∫0tμui du+∑a=1d∫0tσuia dWua(i=1,…,n)X^i_t=X^i_0+\int_0^t\mu^i_u\,du+\sum_{a=1}^d\int_0^t\sigma^{ia}_u\,dW^a_u\qquad(i=1,\dots,n)

と書けるとき(X0X_0 は F0\mathcal{F}_0-可測、μ=(μi)\mu=(\mu^i) と σ=(σia)\sigma=(\sigma^{ia}) は発展的可測で σia∈L2[0,T]\sigma^{ia}\in\mathcal{L}^2[0,T]、∫0T∣μu∣ du<∞\int_0^T|\mu_u|\,du<\infty)、XX を nn 次元の伊藤過程と呼びます。σ\sigma の第 ii 行を σi=(σi1,…,σid)\sigma^i=(\sigma^{i1},\dots,\sigma^{id}) と書きます。さらに定数 KK があって、すべての (u,ω)(u,\omega) と i,ai,a について ∣μui∣≤K|\mu^i_u|\le K、∣σuia∣≤K|\sigma^{ia}_u|\le K となるとき、XX は有界係数の条件 (A) を満たすといいます。

確率積分は連続なバージョンを取ります(定理 4.7[確率微分方程式])。以下で繰り返し使う二つの補題を先に示します。

補題 2.3(有界な単純過程による近似)✓

HH を Rd\mathbb{R}^d 値の発展的可測過程で、定数 LL について ∣Hu(ω)∣≤L|H_u(\omega)|\le L(Euclid ノルム)とします。このとき、[0,T][0,T] の等分割の上の単純過程(各成分が 定義 4.1[確率微分方程式] の意味で単純)H(n)H^{(n)} で、∣H(n)∣≤L|H^{(n)}|\le L かつ E∫0T∣Hu(n)−Hu∣2 du→0\mathbb{E}\int_0^T|H^{(n)}_u-H_u|^2\,du\to0 となるものが存在します。

証明(補題 2.3)

補題 4.4[確率微分方程式] の証明の第 2 段階と第 3 段階を、成分ごとに同じ分割で行えば得られます。確かめるべきは上界 LL が保たれることだけです。第 2 段階の平均化 Ht(m)=m∫t−1/mtHs dsH^{(m)}_t=m\int_{t-1/m}^tH_s\,ds では、ベクトル値積分のノルムが被積分関数のノルムの積分以下であることから ∣Ht(m)∣≤m∫t−1/mt∣Hs∣ ds≤L|H^{(m)}_t|\le m\int_{t-1/m}^t|H_s|\,ds\le L です。第 3 段階の左端離散化は HH の値をそのまま使うので上界を変えません。第 1 段階(切り落とし)は HH がもともと有界なので不要です。各段階の収束は同補題の証明と同じく優収束定理によります。

補題 2.4(確率積分からの既知量の取り出し)✓

H∈L2[0,T]H\in\mathcal{L}^2[0,T]、0≤s<t≤T0\le s<t\le T、ξ\xi を有界な Fs\mathcal{F}_s-可測確率変数とします。このとき ξH1(s,t]∈L2[0,T]\xi H\mathbf{1}_{(s,t]}\in\mathcal{L}^2[0,T] で、各 aa について

∫stξHu dWua=ξ∫stHu dWuaa.s.\int_s^t\xi H_u\,dW^a_u=\xi\int_s^tH_u\,dW^a_u\quad\text{a.s.}
証明(補題 2.4)

過程 (u,ω)↦ξ(ω)1(s,t](u)(u,\omega)\mapsto\xi(\omega)\mathbf{1}_{(s,t]}(u) は、r≤sr\le s なら [0,r][0,r] 上で 00、r>sr>s なら B([0,r])⊗Fs\mathcal{B}([0,r])\otimes\mathcal{F}_s-可測なので発展的可測です。発展的可測な過程の積は発展的可測で、∣ξH1(s,t]∣≤∥ξ∥∞∣H∣|\xi H\mathbf{1}_{(s,t]}|\le\|\xi\|_\infty|H| なので ξH1(s,t]∈L2[0,T]\xi H\mathbf{1}_{(s,t]}\in\mathcal{L}^2[0,T] です。

HH が単純過程で、その分割点に s,ts,t が含まれる場合(分割を細かくしても単純過程のままなので一般性を失いません)、H1(s,t]=∑ti≥sξi1(ti,ti+1]H\mathbf{1}_{(s,t]}=\sum_{t_i\ge s}\xi_i\mathbf{1}_{(t_i,t_{i+1}]}(和は ti+1≤tt_{i+1}\le t の範囲)で、ξξi\xi\xi_i は ti≥st_i\ge s より Fti\mathcal{F}_{t_i}-可測かつ二乗可積分です。よって 定義 4.2[確率微分方程式] から両辺はともに ∑ξξi(Wti+1a−Wtia)\sum\xi\xi_i(W^a_{t_{i+1}}-W^a_{t_i}) です。一般の HH は 補題 4.4[確率微分方程式] の単純過程 H(n)→HH^{(n)}\to H で近似します。∥ξH(n)1(s,t]−ξH1(s,t]∥L2≤∥ξ∥∞∥H(n)−H∥L2→0\|\xi H^{(n)}\mathbf{1}_{(s,t]}-\xi H\mathbf{1}_{(s,t]}\|_{\mathcal{L}^2}\le\|\xi\|_\infty\|H^{(n)}-H\|_{\mathcal{L}^2}\to0 なので、伊藤等長性(定理 4.7[確率微分方程式])により左辺も右辺も L2(Ω)L^2(\Omega) で収束し、極限でも等式が成り立ちます。

3.1. 異なる成分による積分の直交性

Section titled “3.1. 異なる成分による積分の直交性”

この節のねらいは、異なる成分の Brown 運動による積分どうしの積の期待値が 00 になることを示し、多次元の計算規則 dWa dWb=δab dtdW^a\,dW^b=\delta_{ab}\,dt の中身をはっきりさせることです。証明は、単純過程(区間ごとに値が一定の過程)で確かめてから一般の場合へ近づける、という伊藤積分の定義と同じ流れで進みます。途中の「偏極」は、xy=14((x+y)2−(x−y)2)xy=\frac14\bigl((x+y)^2-(x-y)^2\bigr) という恒等式を使って、2 乗についての等長性の式から積についての式を取り出す手順です。

補題 3.1(条件付き共分散等長性)✓

H,K∈L2[0,T]H,K\in\mathcal{L}^2[0,T]、a,b∈{1,…,d}a,b\in\{1,\dots,d\}、0≤s<t≤T0\le s<t\le T、ξ\xi を有界な Fs\mathcal{F}_s-可測確率変数とします。このとき

E[ξ∫stHu dWua∫stKu dWub]=δab E[ξ∫stHuKu du].\mathbb{E}\left[\xi\int_s^tH_u\,dW^a_u\int_s^tK_u\,dW^b_u\right]=\delta_{ab}\,\mathbb{E}\left[\xi\int_s^tH_uK_u\,du\right].

特に H,G∈Ld2H,G\in\mathcal{L}^2_d について E[∫0TH⋅dW∫0TG⋅dW]=E∫0THu⋅Gu du\mathbb{E}[\int_0^TH\cdot dW\int_0^TG\cdot dW]=\mathbb{E}\int_0^TH_u\cdot G_u\,du です。

証明(補題 3.1)

補題 2.4 により ξ∫stH dWa=∫stξH dWa\xi\int_s^tH\,dW^a=\int_s^t\xi H\,dW^a なので、ξ=1\xi=1 の場合に、HH を ξH\xi H に置き換えて示せば十分です。また H,KH,K を H1(s,t],K1(s,t]H\mathbf{1}_{(s,t]},K\mathbf{1}_{(s,t]} に置き換えます。

a=ba=b のときは、伊藤等長性(定理 4.7[確率微分方程式])を H+KH+K と H−KH-K に適用して差を取る(偏極)と E[I(H)I(K)]=E∫HK du\mathbb{E}[I(H)I(K)]=\mathbb{E}\int HK\,du が出ます。

a≠ba\ne b のときを示します。まず H,KH,K が共通の分割 s=u0<⋯<um=ts=u_0<\dots<u_m=t 上の単純過程 H=∑lξl1(ul,ul+1]H=\sum_l\xi_l\mathbf{1}_{(u_l,u_{l+1}]}、K=∑lηl1(ul,ul+1]K=\sum_l\eta_l\mathbf{1}_{(u_l,u_{l+1}]} の場合です。ΔlWa=Wul+1a−Wula\Delta_lW^a=W^a_{u_{l+1}}-W^a_{u_l} と書くと、示すべきは E[∑l,mξlΔlWa ηmΔmWb]=0\mathbb{E}[\sum_{l,m}\xi_l\Delta_lW^a\,\eta_m\Delta_mW^b]=0 です。l<ml<m の項では ξlΔlWaηm\xi_l\Delta_lW^a\eta_m が Fum\mathcal{F}_{u_m}-可測で可積分、ΔmWb\Delta_mW^b は Fum\mathcal{F}_{u_m} と独立で平均 00 なので、定理 4.4[条件付き期待値] と 命題 4.5[条件付き期待値] から期待値は E[ξlΔlWaηm E[ΔmWb]]=0\mathbb{E}[\xi_l\Delta_lW^a\eta_m\,\mathbb{E}[\Delta_mW^b]]=0 です。l>ml>m も同様です。l=ml=m の項では、(ΔlWa,ΔlWb)(\Delta_lW^a,\Delta_lW^b) が Ful\mathcal{F}_{u_l} と独立で、成分どうしも独立で平均 00 なので、E[ξlηlΔlWaΔlWb]=E[ξlηl] E[ΔlWa] E[ΔlWb]=0\mathbb{E}[\xi_l\eta_l\Delta_lW^a\Delta_lW^b]=\mathbb{E}[\xi_l\eta_l]\,\mathbb{E}[\Delta_lW^a]\,\mathbb{E}[\Delta_lW^b]=0 です(命題 7.5[確率変数と期待値])。

一般の H,KH,K は単純過程 H(n)→HH^{(n)}\to H、K(n)→KK^{(n)}\to K で近似します(補題 4.4[確率微分方程式]。分割は両者の和集合に細分して共通にできます)。Cauchy–Schwarz の不等式と等長性から

∣E[Ia(H)Ib(K)]−E[Ia(H(n))Ib(K(n))]∣≤∥H−H(n)∥L2∥K∥L2+∥H(n)∥L2∥K−K(n)∥L2→0\bigl|\mathbb{E}[I^a(H)I^b(K)]-\mathbb{E}[I^a(H^{(n)})I^b(K^{(n)})]\bigr|\le\|H-H^{(n)}\|_{\mathcal{L}^2}\|K\|_{\mathcal{L}^2}+\|H^{(n)}\|_{\mathcal{L}^2}\|K-K^{(n)}\|_{\mathcal{L}^2}\to0

(Ia(H)=∫H dWaI^a(H)=\int H\,dW^a)なので、00 という値が極限に受け継がれます。最後の主張は、成分ごとの式を a,ba,b について足し合わせたものです。

この補題を形式的に書けば dWua dWub=δab dudW^a_u\,dW^b_u=\delta_{ab}\,du です。条件付きの形(ξ\xi が入っている形)で述べたのは、次の命題と 2 次変分の収束で、マルチンゲール差の直交性として使うためです。

命題 3.2(確率積分の積から作られるマルチンゲール)✓

G,H∈Ld2G,H\in\mathcal{L}^2_d とし、Mt=∫0tG⋅dWM_t=\int_0^tG\cdot dW、Nt=∫0tH⋅dWN_t=\int_0^tH\cdot dW とおきます。このとき

MtNt−∫0tGu⋅Hu du(t∈[0,T])M_tN_t-\int_0^tG_u\cdot H_u\,du\qquad(t\in[0,T])

は (Ft)(\mathcal{F}_t)-マルチンゲールです。

証明(命題 3.2)

Mt,Nt∈L2(Ω)M_t,N_t\in L^2(\Omega) なので積は可積分で、∫0T∣G⋅H∣ du\int_0^T|G\cdot H|\,du も Cauchy–Schwarz の不等式から可積分です。s<ts<t、A∈FsA\in\mathcal{F}_s とし、ΔM=Mt−Ms\Delta M=M_t-M_s、ΔN=Nt−Ns\Delta N=N_t-N_s と書くと

MtNt−MsNs=ΔM ΔN+Ms ΔN+Ns ΔM.M_tN_t-M_sN_s=\Delta M\,\Delta N+M_s\,\Delta N+N_s\,\Delta M.

ΔN\Delta N は 定理 4.7[確率微分方程式] のマルチンゲール性から E[ΔN∣Fs]=0\mathbb{E}[\Delta N\mid\mathcal{F}_s]=0 を満たし、1AMs\mathbf{1}_AM_s は二乗可積分で Fs\mathcal{F}_s-可測なので、定理 4.4[条件付き期待値] により E[1AMsΔN]=0\mathbb{E}[\mathbf{1}_AM_s\Delta N]=0 です。第 3 項も同じです。第 1 項は ΔM ΔN=∑a,b∫stGa dWa∫stHb dWb\Delta M\,\Delta N=\sum_{a,b}\int_s^tG^a\,dW^a\int_s^tH^b\,dW^b と展開し、補題 3.1(ξ=1A\xi=\mathbf{1}_A)を使うと E[1AΔMΔN]=E[1A∫stG⋅H du]\mathbb{E}[\mathbf{1}_A\Delta M\Delta N]=\mathbb{E}[\mathbf{1}_A\int_s^tG\cdot H\,du] です。合わせて E[1A(MtNt−∫0tG⋅H)]=E[1A(MsNs−∫0sG⋅H)]\mathbb{E}[\mathbf{1}_A(M_tN_t-\int_0^tG\cdot H)]=\mathbb{E}[\mathbf{1}_A(M_sN_s-\int_0^sG\cdot H)] を得ます。

この命題が、次の定義を正当化します。

定義 3.3(伊藤過程の共変分)↗

定義 2.2 の伊藤過程 XX について

⟨Xi,Xj⟩t:=∫0tσui⋅σuj du=∫0t(σuσuT)ij du\langle X^i,X^j\rangle_t:=\int_0^t\sigma^i_u\cdot\sigma^j_u\,du=\int_0^t(\sigma_u\sigma_u^{\mathsf T})^{ij}\,du

を XiX^i と XjX^j の共変分、⟨Xi⟩:=⟨Xi,Xi⟩\langle X^i\rangle:=\langle X^i,X^i\rangle を XiX^i の 2 次変分と呼びます。ドリフト μ\mu は共変分に寄与しません。

命題 3.2 により、XX のマルチンゲール部分 Mi=∫σi⋅dWM^i=\int\sigma^i\cdot dW について MiMj−⟨Xi,Xj⟩M^iM^j-\langle X^i,X^j\rangle はマルチンゲールです。たとえば ⟨Wa,Wb⟩t=δabt\langle W^a,W^b\rangle_t=\delta_{ab}t です。

注意 3.4

一般の連続局所マルチンゲール M,NM,N に対しても、MN−AMN-A が局所マルチンゲールとなる連続な有界変動過程 AA(A0=0A_0=0)がただ一つ存在し、これを ⟨M,N⟩\langle M,N\rangle と書きます。存在の証明には Doob–Meyer 分解か、離散 2 次変分の一様収束の議論が要ります(Karatzas–Shreve 第 1.5 節・第 3.2 節、Revuz–Yor 第 IV 章)。この章では伊藤過程に限るので、上の定義で足ります。

2 次変分の収束を示すには、増分の 2 乗の揺らぎを抑える必要があり、そのために 4 次モーメントの評価を使います。

次の補題の証明は二段構えです。まず単純過程について、各区間の寄与 YlY_l が、それまでの情報を固定すると正規分布に従うことを使い、部分和 SlS_l の 4 乗の期待値を区間ごとに足し上げます。そのうえで、一般の場合を単純過程の極限として得ます。

補題 3.5(伊藤積分の 4 次モーメント)✓

HH を Rd\mathbb{R}^d 値の発展的可測過程で ∣Hu∣≤L|H_u|\le L とし、0≤s<t≤T0\le s<t\le T とします。このとき

E[(∫stHu⋅dWu)4]≤6L4(t−s)2.\mathbb{E}\left[\left(\int_s^tH_u\cdot dW_u\right)^4\right]\le6L^4(t-s)^2.
証明(補題 3.5)

単純過程の場合。 H=∑lξl1(ul,ul+1]H=\sum_{l}\xi_l\mathbf{1}_{(u_l,u_{l+1}]}(s=u0<⋯<um=ts=u_0<\dots<u_m=t、ξl\xi_l は Ful\mathcal{F}_{u_l}-可測で ∣ξl∣≤L|\xi_l|\le L)とし、Yl=ξl⋅(Wul+1−Wul)Y_l=\xi_l\cdot(W_{u_{l+1}}-W_{u_l})、Sl=∑k<lYkS_l=\sum_{k<l}Y_k、hl=ul+1−ulh_l=u_{l+1}-u_l とおきます。Wul+1−WulW_{u_{l+1}}-W_{u_l} は Ful\mathcal{F}_{u_l} と独立に N(0,hlId)N(0,h_lI_d) に従うので、補題 2.1[Feller 過程と Feynman–Kac](関数の正部分と負部分に分けて適用)により、Ful\mathcal{F}_{u_l} に関する条件付きで YlY_l は N(0,∣ξl∣2hl)N(0,|\xi_l|^2h_l) に従います。したがって

E[Yl∣Ful]=E[Yl3∣Ful]=0,E[Yl2∣Ful]=∣ξl∣2hl,E[Yl4∣Ful]=3∣ξl∣4hl2.\mathbb{E}[Y_l\mid\mathcal{F}_{u_l}]=\mathbb{E}[Y_l^3\mid\mathcal{F}_{u_l}]=0,\quad\mathbb{E}[Y_l^2\mid\mathcal{F}_{u_l}]=|\xi_l|^2h_l,\quad\mathbb{E}[Y_l^4\mid\mathcal{F}_{u_l}]=3|\xi_l|^4h_l^2 .

∣Sl∣≤L∑k<l∣Wuk+1−Wuk∣|S_l|\le L\sum_{k<l}|W_{u_{k+1}}-W_{u_k}| はすべての次数のモーメントを持つので、(Sl+Yl)4(S_l+Y_l)^4 を展開して 定理 4.4[条件付き期待値] を使うと、奇数次の項が消えて

E[Sl+14]=E[Sl4]+6 E[Sl2∣ξl∣2] hl+3 E[∣ξl∣4] hl2≤E[Sl4]+6L2 E[Sl2] hl+3L4hl2.\mathbb{E}[S_{l+1}^4]=\mathbb{E}[S_l^4]+6\,\mathbb{E}[S_l^2|\xi_l|^2]\,h_l+3\,\mathbb{E}[|\xi_l|^4]\,h_l^2\le\mathbb{E}[S_l^4]+6L^2\,\mathbb{E}[S_l^2]\,h_l+3L^4h_l^2 .

同様に E[Sl2]=∑k<lE[∣ξk∣2]hk≤L2(ul−s)\mathbb{E}[S_l^2]=\sum_{k<l}\mathbb{E}[|\xi_k|^2]h_k\le L^2(u_l-s) です。これを l=0,…,m−1l=0,\dots,m-1 について足すと、∑l(ul−s)hl≤∫st(u−s) du=(t−s)2/2\sum_l(u_l-s)h_l\le\int_s^t(u-s)\,du=(t-s)^2/2 と ∑lhl2≤(t−s)2\sum_lh_l^2\le(t-s)^2 から

E[Sm4]≤6L4⋅(t−s)22+3L4(t−s)2=6L4(t−s)2.\mathbb{E}[S_m^4]\le6L^4\cdot\frac{(t-s)^2}{2}+3L^4(t-s)^2=6L^4(t-s)^2 .

一般の場合。 補題 2.3 により ∣H(n)∣≤L|H^{(n)}|\le L なる単純過程で E∫∣H(n)−H∣2→0\mathbb{E}\int|H^{(n)}-H|^2\to0 となるものを取ります(分割に s,ts,t を加えます)。補題 3.1 から ∫stH(n)⋅dW→∫stH⋅dW\int_s^tH^{(n)}\cdot dW\to\int_s^tH\cdot dW が L2(Ω)L^2(\Omega) で成り立つので、ほとんど確実に収束する部分列が取れます。実際、Yn=∫stH(n)⋅dWY_n=\int_s^tH^{(n)}\cdot dW、Y=∫stH⋅dWY=\int_s^tH\cdot dW として E∣Ynk−Y∣2≤8−k\mathbb{E}|Y_{n_k}-Y|^2\le 8^{-k} となる n1<n2<⋯n_1<n_2<\cdots を選ぶと、マルコフの不等式(命題 3.5[大数の法則と中心極限定理])を ∣Ynk−Y∣2|Y_{n_k}-Y|^2 と a=4−ka=4^{-k} に使って P(∣Ynk−Y∣≥2−k)≤4k⋅8−k=2−k\mathbb{P}(|Y_{n_k}-Y|\ge2^{-k})\le4^k\cdot8^{-k}=2^{-k} で、これは kk について和が有限です。ボレル・カンテリの補題(補題 4.4[大数の法則と中心極限定理])から、ほとんど確実に有限個の kk を除いて ∣Ynk−Y∣<2−k|Y_{n_k}-Y|<2^{-k} となり、Ynk→YY_{n_k}\to Y a.s. です。補題 6.1[ルベーグ積分の定義と収束定理] により、極限でも同じ上界が成り立ちます。

上界が (t−s)2(t-s)^2 の定数倍であることが重要です。積分の 2 乗は t−st-s の大きさなので、4 次モーメントはその 2 乗の大きさに収まり、増分の 2 乗の揺らぎは平均に比べて小さいことになります。

定義 3.3 の共変分が、増分の積の和の極限として実際に現れることを示します。これが 1 次元の 定理 3.1[確率微分方程式] の多次元・伊藤過程版で、伊藤公式の 2 次の項の出どころです。

次の定理の証明では、増分 ΔXkiΔXkj\Delta X_k^i\Delta X_k^j の和を、増分をドリフト部分とマルチンゲール部分に分けて展開し、ドリフトを含む項、マルチンゲール部分の項、Riemann 和の項に整理して、それぞれ 00 または目的の積分に収束することを見ます。ドリフトを含む項は、ドリフトの寄与が時間幅に比例して小さいため、分割を細かくすれば 00 に近づきます。マルチンゲール部分の項は、各区間の誤差 ZkZ_k が互いに直交するので 2 乗の和が小さくなります。残る Riemann 和の項は、連続な被積分関数の積分への通常の収束です。

定理 4.1(共変分の離散近似)✓

XX を有界係数の条件 (A)(定数 KK)を満たす nn 次元伊藤過程、GG を適合した実数値過程で、道が連続かつ ∣G∣≤C|G|\le C なるものとします。t∈(0,T]t\in(0,T]、Π={0=t0<⋯<tm=t}\Pi=\{0=t_0<\dots<t_m=t\} を分割、ΔXki=Xtk+1i−Xtki\Delta X^i_k=X^i_{t_{k+1}}-X^i_{t_k} とおくと

E∣∑kGtk ΔXki ΔXkj−∫0tGu d⟨Xi,Xj⟩u∣⟶0(∥Π∥→0).\mathbb{E}\left|\sum_{k}G_{t_k}\,\Delta X^i_k\,\Delta X^j_k-\int_0^tG_u\,d\langle X^i,X^j\rangle_u\right|\longrightarrow0\qquad(\|\Pi\|\to0).

ここで ∫0tGu d⟨Xi,Xj⟩u=∫0tGu σui⋅σuj du\int_0^tG_u\,d\langle X^i,X^j\rangle_u=\int_0^tG_u\,\sigma^i_u\cdot\sigma^j_u\,du です。

証明(定理 4.1)

hk=tk+1−tkh_k=t_{k+1}-t_k とし、増分をドリフト部分とマルチンゲール部分に分けます。

ΔXki=Aki+Nki,Aki=∫tktk+1μui du,Nki=∫tktk+1σui⋅dWu.\Delta X^i_k=A^i_k+N^i_k,\qquad A^i_k=\int_{t_k}^{t_{k+1}}\mu^i_u\,du,\qquad N^i_k=\int_{t_k}^{t_{k+1}}\sigma^i_u\cdot dW_u .

∣Aki∣≤Khk|A^i_k|\le Kh_k で、補題 3.1 と ∣σi∣2≤dK2|\sigma^i|^2\le dK^2 から E[(Nki)2]≤dK2hk\mathbb{E}[(N^i_k)^2]\le dK^2h_k です。

ドリフトを含む項。 ∑k∣GtkAkiAkj∣≤CK2∥Π∥t\sum_k|G_{t_k}A^i_kA^j_k|\le CK^2\|\Pi\|t です。また Cauchy–Schwarz の不等式から

E∑k∣GtkAkiNkj∣≤CK∑khk(E[(Nkj)2])1/2≤Cd K2∑khk3/2≤Cd K2∥Π∥1/2t,\mathbb{E}\sum_k|G_{t_k}A^i_kN^j_k|\le CK\sum_kh_k\bigl(\mathbb{E}[(N^j_k)^2]\bigr)^{1/2}\le C\sqrt d\,K^2\sum_kh_k^{3/2}\le C\sqrt d\,K^2\|\Pi\|^{1/2}t,

AjNiA^jN^i の項も同じです。いずれも 00 に収束します。

マルチンゲール部分の項。 Ck=∫tktk+1σi⋅σj duC_k=\int_{t_k}^{t_{k+1}}\sigma^i\cdot\sigma^j\,du、Zk=Gtk(NkiNkj−Ck)Z_k=G_{t_k}(N^i_kN^j_k-C_k) とおきます。∣Ck∣≤dK2hk|C_k|\le dK^2h_k で、補題 3.5(L=dKL=\sqrt dK)と Cauchy–Schwarz の不等式から

E[Zk2]≤2C2(E[(Nki)4]1/2 E[(Nkj)4]1/2+d2K4hk2)≤14 C2d2K4hk2.\mathbb{E}[Z_k^2]\le2C^2\Bigl(\mathbb{E}[(N^i_k)^4]^{1/2}\,\mathbb{E}[(N^j_k)^4]^{1/2}+d^2K^4h_k^2\Bigr)\le14\,C^2d^2K^4h_k^2 .

次に、ξ\xi を有界な Ftk\mathcal{F}_{t_k}-可測確率変数とすると、NkiNkj=∑a,b∫σia dWa∫σjb dWbN^i_kN^j_k=\sum_{a,b}\int\sigma^{ia}\,dW^a\int\sigma^{jb}\,dW^b に 補題 3.1 を ξGtk\xi G_{t_k} について適用して E[ξZk]=0\mathbb{E}[\xi Z_k]=0 を得ます。ξ\xi として Ftk\mathcal{F}_{t_k} の事象の指示関数を取れば E[Zk∣Ftk]=0\mathbb{E}[Z_k\mid\mathcal{F}_{t_k}]=0 です。k<lk<l なら ZkZ_k は Ftl\mathcal{F}_{t_l}-可測で二乗可積分なので、定理 4.4[条件付き期待値] により E[ZkZl]=E[Zk E[Zl∣Ftl]]=0\mathbb{E}[Z_kZ_l]=\mathbb{E}[Z_k\,\mathbb{E}[Z_l\mid\mathcal{F}_{t_l}]]=0 です。よって交差項が消え、

E[(∑kZk)2]=∑kE[Zk2]≤14 C2d2K4∥Π∥t⟶0.\mathbb{E}\Bigl[\Bigl(\sum_kZ_k\Bigr)^2\Bigr]=\sum_k\mathbb{E}[Z_k^2]\le14\,C^2d^2K^4\|\Pi\|t\longrightarrow0 .

Riemann 和の項。 残るのは ∑kGtkCk=∫0tGuΠ σui⋅σuj du\sum_kG_{t_k}C_k=\int_0^tG^\Pi_u\,\sigma^i_u\cdot\sigma^j_u\,du(GuΠ=GtkG^\Pi_u=G_{t_k}, u∈(tk,tk+1]u\in(t_k,t_{k+1}])と ∫0tGu σi⋅σj du\int_0^tG_u\,\sigma^i\cdot\sigma^j\,du の差です。道の連続性から GuΠ→GuG^\Pi_u\to G_u が各 (u,ω)(u,\omega) で成り立ち、被積分関数は 2CdK22CdK^2 で抑えられるので、定理 7.3[ルベーグ積分の定義と収束定理](du⊗dPdu\otimes d\mathbb{P} について)から差の期待値は 00 に収束します。

以上の 4 種類の項を合わせると主張が出ます(L2L^2 収束は L1L^1 収束を導きます)。

G=1G=1、X=WX=W とすれば ∑kΔWkaΔWkb→δabt\sum_k\Delta W^a_k\Delta W^b_k\to\delta_{ab}t で、異なる成分の増分の積の和が 00 に収束することが分かります。符号がランダムな小さい量 ΔWkaΔWkb\Delta W^a_k\Delta W^b_k(大きさ hkh_k)が 1/hk1/h_k 個足されても、互いに打ち消し合って和は ∥Π∥\sqrt{\|\Pi\|} 程度にしかならない、というのが証明の中身です。

通常の微積分の連鎖律では ff の 1 次の項だけが残りますが、確率過程では増分の 2 乗が時間幅と同じ大きさ(h\sqrt{h} の 2 乗)になるため、2 次の項も残ります。次の定理の証明は、時間を細かく刻んで ff の変化を足し合わせ(望遠鏡和)、各区間で 2 次まで Taylor 展開します。1 次の項は確率積分とドリフトの積分に、2 次の項は前節の定理により共変分の積分に収束します。3 次以上の項(剰余)は、増分の 3 乗が h3/2h^{3/2} の大きさなので、足し合わせても 00 に近づきます。

定理 5.1(多次元伊藤公式)✓

XX を有界係数の条件 (A) を満たす nn 次元伊藤過程とし、f∈C1,2([0,T]×Rn)f\in C^{1,2}([0,T]\times\mathbb{R}^n)(tt について 1 回、xx について 2 回連続微分可能)で、xx についての 1 階・2 階の偏導関数 ∂if\partial_if、∂ijf\partial_{ij}f が有界であるとします。このとき、ほとんど確実にすべての t∈[0,T]t\in[0,T] について

f(t,Xt)=f(0,X0)+∫0t∂tf(u,Xu) du+∑i∫0t∂if(u,Xu) μui du+∑i,a∫0t∂if(u,Xu) σuia dWua+12∑i,j∫0t∂ijf(u,Xu) d⟨Xi,Xj⟩u\begin{aligned} f(t,X_t)=f(0,X_0)&+\int_0^t\partial_tf(u,X_u)\,du+\sum_{i}\int_0^t\partial_if(u,X_u)\,\mu^i_u\,du\\ &+\sum_{i,a}\int_0^t\partial_if(u,X_u)\,\sigma^{ia}_u\,dW^a_u+\frac12\sum_{i,j}\int_0^t\partial_{ij}f(u,X_u)\,d\langle X^i,X^j\rangle_u \end{aligned}

が成り立ちます。微分形式では df=∂tf dt+∑i∂if dXi+12∑i,j∂ijf d⟨Xi,Xj⟩df=\partial_tf\,dt+\sum_i\partial_if\,dX^i+\frac12\sum_{i,j}\partial_{ij}f\,d\langle X^i,X^j\rangle です。

証明(定理 5.1)

被積分関数 ∂if(u,Xu)σuia\partial_if(u,X_u)\sigma^{ia}_u は有界で、連続な適合過程と発展的可測過程の積なので発展的可測です。よって右辺の確率積分は定義されます。両辺は tt について連続なので、tt を固定して等式がほとんど確実に成り立つことを示せば、有理数の tt で同時に成り立ち、連続性からすべての tt に及びます。

tt を固定し、分割 Π={0=t0<⋯<tm=t}\Pi=\{0=t_0<\dots<t_m=t\} を取ります。Xk=XtkX_k=X_{t_k}、ΔXk=Xk+1−Xk\Delta X_k=X_{k+1}-X_k、hk=tk+1−tkh_k=t_{k+1}-t_k と書き、望遠鏡和を時間方向と空間方向に分けます。

f(t,Xt)−f(0,X0)=∑k[f(tk+1,Xk+1)−f(tk,Xk+1)]+∑k[f(tk,Xk+1)−f(tk,Xk)].f(t,X_t)-f(0,X_0)=\sum_k\bigl[f(t_{k+1},X_{k+1})-f(t_k,X_{k+1})\bigr]+\sum_k\bigl[f(t_k,X_{k+1})-f(t_k,X_k)\bigr].

以下、道 ω\omega を固定したときの量 R(ω)=sup⁡u≤t∣Xu(ω)∣R(\omega)=\sup_{u\le t}|X_u(\omega)| は有限で、[0,T]×{∣x∣≤R+1}[0,T]\times\{|x|\le R+1\} はコンパクトなので ∂tf\partial_tf と ∂ijf\partial_{ij}f はそこで一様連続です。

時間方向。 第 1 の和の各項は ∫tktk+1∂tf(u,Xk+1) du\int_{t_k}^{t_{k+1}}\partial_tf(u,X_{k+1})\,du です。これと ∫tktk+1∂tf(u,Xu) du\int_{t_k}^{t_{k+1}}\partial_tf(u,X_u)\,du の差の和は、t⋅max⁡ksup⁡u∈[tk,tk+1]∣∂tf(u,Xk+1)−∂tf(u,Xu)∣t\cdot\max_k\sup_{u\in[t_k,t_{k+1}]}|\partial_tf(u,X_{k+1})-\partial_tf(u,X_u)| 以下で、∂tf\partial_tf の一様連続性と道の一様連続性から各 ω\omega で 00 に収束します。

空間方向の Taylor 展開。 積分形の剰余をもつ Taylor の定理から

f(tk,Xk+1)−f(tk,Xk)=∑i∂if(tk,Xk)ΔXki+12∑i,j∂ijf(tk,Xk)ΔXkiΔXkj+ρk,f(t_k,X_{k+1})-f(t_k,X_k)=\sum_i\partial_if(t_k,X_k)\Delta X^i_k+\frac12\sum_{i,j}\partial_{ij}f(t_k,X_k)\Delta X^i_k\Delta X^j_k+\rho_k,∣ρk∣≤12 εΠ ∣ΔXk∣2,εΠ:=max⁡ksup⁡0≤θ≤1∥∇2f(tk,Xk+θΔXk)−∇2f(tk,Xk)∥|\rho_k|\le\frac12\,\varepsilon_\Pi\,|\Delta X_k|^2,\qquad\varepsilon_\Pi:=\max_k\sup_{0\le\theta\le1}\bigl\|\nabla^2f(t_k,X_k+\theta\Delta X_k)-\nabla^2f(t_k,X_k)\bigr\|

です(∥⋅∥\|\cdot\| は行列の作用素ノルム)。

1 次の項。 左端の値で作った過程 YuΠ,i:=∂if(tk,Xk)Y^{\Pi,i}_u:=\partial_if(t_k,X_k)(u∈(tk,tk+1]u\in(t_k,t_{k+1}])を使うと、補題 2.4(ξ=∂if(tk,Xk)\xi=\partial_if(t_k,X_k))により

∑k∂if(tk,Xk)ΔXki=∫0tYuΠ,iμui du+∑a∫0tYuΠ,iσuia dWua.\sum_k\partial_if(t_k,X_k)\Delta X^i_k=\int_0^tY^{\Pi,i}_u\mu^i_u\,du+\sum_a\int_0^tY^{\Pi,i}_u\sigma^{ia}_u\,dW^a_u .

道の連続性から YuΠ,i→∂if(u,Xu)Y^{\Pi,i}_u\to\partial_if(u,X_u) が各 (u,ω)(u,\omega) で成り立ち、∣YΠ,i∣≤sup⁡∣∂if∣|Y^{\Pi,i}|\le\sup|\partial_if| です。ドリフトの積分は優収束定理により各 ω\omega で収束します。確率積分は、E∫0t∣YuΠ,i−∂if(u,Xu)∣2(σuia)2 du→0\mathbb{E}\int_0^t|Y^{\Pi,i}_u-\partial_if(u,X_u)|^2(\sigma^{ia}_u)^2\,du\to0(優収束定理)と伊藤等長性から L2(Ω)L^2(\Omega) で収束します。

2 次の項。 Gu=∂ijf(u,Xu)G_u=\partial_{ij}f(u,X_u) は有界で道が連続な適合過程なので、定理 4.1 により ∑k∂ijf(tk,Xk)ΔXkiΔXkj→∫0t∂ijf(u,Xu) d⟨Xi,Xj⟩u\sum_k\partial_{ij}f(t_k,X_k)\Delta X^i_k\Delta X^j_k\to\int_0^t\partial_{ij}f(u,X_u)\,d\langle X^i,X^j\rangle_u が L1L^1 で成り立ちます。

剰余項。 SΠ:=∑k∣ΔXk∣2=∑i∑k(ΔXki)2S_\Pi:=\sum_k|\Delta X_k|^2=\sum_i\sum_k(\Delta X^i_k)^2 は 定理 4.1(G=1G=1)により ∑i⟨Xi⟩t\sum_i\langle X^i\rangle_t に L1L^1 で収束するので、sup⁡ΠE[SΠ]<∞\sup_\Pi\mathbb{E}[S_\Pi]<\infty です。一方 max⁡k∣ΔXk∣→0\max_k|\Delta X_k|\to0(道の一様連続性)と ∇2f\nabla^2f の一様連続性から εΠ→0\varepsilon_\Pi\to0 が各 ω\omega で成り立ちます。任意の δ,η>0\delta,\eta>0 について、定理 6.1[確率変数と期待値] から

P(∣∑kρk∣>η)≤P(εΠ>δ)+P(SΠ>2η/δ)≤P(εΠ>δ)+δ E[SΠ]2η\mathbb{P}\Bigl(\bigl|\textstyle\sum_k\rho_k\bigr|>\eta\Bigr)\le\mathbb{P}(\varepsilon_\Pi>\delta)+\mathbb{P}\bigl(S_\Pi>2\eta/\delta\bigr)\le\mathbb{P}(\varepsilon_\Pi>\delta)+\frac{\delta\,\mathbb{E}[S_\Pi]}{2\eta}

で、∥Π∥→0\|\Pi\|\to0 の後に δ→0\delta\to0 とすれば右辺は 00 に収束します。

以上より右辺の各項は ∥Π∥→0\|\Pi\|\to0 で確率収束し、その極限の和が公式の右辺です。左辺は Π\Pi によらないので、確率収束の極限の一意性から等式がほとんど確実に成り立ちます。

1 次元の公式(定理 5.2[確率微分方程式])との違いは、2 次の項に ∂ijf\partial_{ij}f と共変分の組が全部現れることだけです。覚え方は dWa dWb=δab dtdW^a\,dW^b=\delta_{ab}\,dt、dt dWa=0dt\,dW^a=0、(dt)2=0(dt)^2=0 のもとで 2 次まで Taylor 展開する、です。

注意 5.2

定理では係数と ff の 1 階・2 階導関数の有界性を仮定しました。これらは停止時刻による局所化で外せます。∫0T(∣μu∣+∣σu∣2) du<∞\int_0^T(|\mu_u|+|\sigma_u|^2)\,du<\infty と f∈C1,2f\in C^{1,2} だけを仮定すると、公式は同じ形で成り立ちます。ただし確率積分は局所マルチンゲールとして定義し直す必要があります(Karatzas–Shreve 第 3.3 節、Revuz–Yor 第 IV.3 節)。下の積の公式は、f(x,y)=xyf(x,y)=xy の導関数が有界でないにもかかわらず、局所化を使わずに直接証明できます。

系 5.3(積の公式)✓

X,YX,Y を、有界係数の条件 (A) を満たす同じ伊藤過程の二つの成分とし(dX=μX dt+σX⋅dWdX=\mu^X\,dt+\sigma^X\cdot dW、dY=μY dt+σY⋅dWdY=\mu^Y\,dt+\sigma^Y\cdot dW)、X0,Y0∈L2(Ω)X_0,Y_0\in L^2(\Omega) とします。このとき XσY,YσX∈Ld2X\sigma^Y,Y\sigma^X\in\mathcal{L}^2_d で、ほとんど確実にすべての tt について

XtYt=X0Y0+∫0tXu dYu+∫0tYu dXu+⟨X,Y⟩t,X_tY_t=X_0Y_0+\int_0^tX_u\,dY_u+\int_0^tY_u\,dX_u+\langle X,Y\rangle_t,

ここで ∫0tXu dYu:=∫0tXuμuY du+∫0tXuσuY⋅dWu\int_0^tX_u\,dY_u:=\int_0^tX_u\mu^Y_u\,du+\int_0^tX_u\sigma^Y_u\cdot dW_u、⟨X,Y⟩t=∫0tσuX⋅σuY du\langle X,Y\rangle_t=\int_0^t\sigma^X_u\cdot\sigma^Y_u\,du です。

証明(系 5.3)

v<wv<w について、補題 3.1 から

E∣Xw−Xv∣2≤2 E(∫vwμX du)2+2 E(∫vwσX⋅dW)2≤2K2(w−v)2+2dK2(w−v).\mathbb{E}|X_w-X_v|^2\le2\,\mathbb{E}\Bigl(\int_v^w\mu^X\,du\Bigr)^2+2\,\mathbb{E}\Bigl(\int_v^w\sigma^X\cdot dW\Bigr)^2\le2K^2(w-v)^2+2dK^2(w-v).

特に E[Xu2]≤2 E[X02]+4K2T2+4dK2T\mathbb{E}[X_u^2]\le2\,\mathbb{E}[X_0^2]+4K^2T^2+4dK^2T は uu について有界で、∣σY∣≤dK|\sigma^Y|\le\sqrt dK なので XσY∈Ld2X\sigma^Y\in\mathcal{L}^2_d です。

分割 Π\Pi について、代数的な恒等式

XtYt−X0Y0=∑kXtkΔYk+∑kYtkΔXk+∑kΔXkΔYkX_tY_t-X_0Y_0=\sum_kX_{t_k}\Delta Y_k+\sum_kY_{t_k}\Delta X_k+\sum_k\Delta X_k\Delta Y_k

が剰余なしに成り立ちます。第 3 の和は 定理 4.1(G=1G=1)により ⟨X,Y⟩t\langle X,Y\rangle_t に L1L^1 で収束します。第 1 の和は、XuΠ:=XtkX^\Pi_u:=X_{t_k}(u∈(tk,tk+1]u\in(t_k,t_{k+1}])とおくと 補題 2.4 の近似版により ∫0tXΠμY du+∫0tXΠσY⋅dW\int_0^tX^\Pi\mu^Y\,du+\int_0^tX^\Pi\sigma^Y\cdot dW に等しいです(XtkX_{t_k} は有界とは限らないので、XtkX_{t_k} を (−N)∨Xtk∧N(-N)\vee X_{t_k}\wedge N で置き換えた等式で N→∞N\to\infty とし、両辺の L2L^2 収束を使います)。ドリフトの積分は道ごとに優収束定理で ∫0tXμY du\int_0^tX\mu^Y\,du に収束します(∣XΠ∣≤sup⁡u≤t∣Xu∣<∞|X^\Pi|\le\sup_{u\le t}|X_u|<\infty)。確率積分については、上の評価から

E∫0t∣XuΠ−Xu∣2∣σuY∣2 du≤dK2 t (2K2∥Π∥2+2dK2∥Π∥)⟶0\mathbb{E}\int_0^t|X^\Pi_u-X_u|^2|\sigma^Y_u|^2\,du\le dK^2\,t\,\bigl(2K^2\|\Pi\|^2+2dK^2\|\Pi\|\bigr)\longrightarrow0

なので、等長性から L2(Ω)L^2(\Omega) で ∫0tXσY⋅dW\int_0^tX\sigma^Y\cdot dW に収束します。第 2 の和も同じです。左辺は Π\Pi によらないので等式が従い、両辺の連続性からすべての tt で同時に成り立ちます。

例 5.4(Brown 運動のノルムの 2 乗)

X=Y=WaX=Y=W^a に 系 5.3 を使うと (Wta)2=2∫0tWa dWa+t(W^a_t)^2=2\int_0^tW^a\,dW^a+t です。aa について足すと

∣Wt∣2=2∫0tWu⋅dWu+d t.|W_t|^2=2\int_0^tW_u\cdot dW_u+d\,t .

右辺の確率積分は平均 00 なので E∣Wt∣2=dt\mathbb{E}|W_t|^2=dt です。また a≠ba\ne b なら ⟨Wa,Wb⟩=0\langle W^a,W^b\rangle=0 なので WtaWtb=∫0tWa dWb+∫0tWb dWaW^a_tW^b_t=\int_0^tW^a\,dW^b+\int_0^tW^b\,dW^a で、WaWbW^aW^b はマルチンゲールです。

例 5.5(後ろ向き熱方程式の解とマルチンゲール)

u∈Rdu\in\mathbb{R}^d を固定し、f(t,x)=cos⁡(u⋅x) e∣u∣2t/2f(t,x)=\cos(u\cdot x)\,e^{|u|^2t/2} とおきます。∂tf=∣u∣22f\partial_tf=\frac{|u|^2}{2}f、∑a∂aaf=−∣u∣2f\sum_a\partial_{aa}f=-|u|^2f なので ∂tf+12Δf=0\partial_tf+\frac12\Delta f=0 です。[0,T][0,T] 上で ff の xx についての 1 階・2 階導関数は有界なので、定理 5.1(X=WX=W、μ=0\mu=0、σ=Id\sigma=I_d)から

cos⁡(u⋅Wt) e∣u∣2t/2=1−∫0tsin⁡(u⋅Ws) e∣u∣2s/2 u⋅dWs.\cos(u\cdot W_t)\,e^{|u|^2t/2}=1-\int_0^t\sin(u\cdot W_s)\,e^{|u|^2s/2}\,u\cdot dW_s .

右辺はマルチンゲールなので E[cos⁡(u⋅Wt)]=e−∣u∣2t/2\mathbb{E}[\cos(u\cdot W_t)]=e^{-|u|^2t/2} で、Gauss 分布の特性関数が伊藤公式から再現されます。一般に ∂tf+12Δf=0\partial_tf+\frac12\Delta f=0 を満たす ff について f(t,Wt)f(t,W_t) はマルチンゲールになり、これが Feynman–Kac 公式(Feller 過程の章)の最も簡単な場合です。

マルチンゲールとブラウン運動の章 の 定理 6.1[マルチンゲールとブラウン運動] で、Brown 運動 BB について BtB_t と Bt2−tB_t^2-t がマルチンゲールであることを見ました。この節ではその逆を証明します。この節では時間を [0,∞)[0,\infty) とします。証明には停止時刻(その時刻までの情報で「もう起きたか」を判定できるランダムな時刻)を使います。

補題 6.1(止めた連続マルチンゲール)✓

(Mt)t≥0(M_t)_{t\ge0} を道が連続な (Ft)(\mathcal{F}_t)-マルチンゲール、τ\tau を停止時刻とします。このとき (Mt∧τ)t≥0(M_{t\wedge\tau})_{t\ge0} も (Ft)(\mathcal{F}_t)-マルチンゲールです。

証明(補題 6.1)

適合性: tt を固定し、σn:=t∧(⌈2nτ⌉/2n)\sigma_n:=t\wedge(\lceil2^n\tau\rceil/2^n) とおくと、σn\sigma_n は tt 以下の有限個の値(2 進有理数と tt)しかとらず、各値 q<tq<t について {σn=q}={q−2−n<τ≤q}∈Fq⊂Ft\{\sigma_n=q\}=\{q-2^{-n}<\tau\le q\}\in\mathcal{F}_q\subset\mathcal{F}_t です。よって MσnM_{\sigma_n} は Ft\mathcal{F}_t-可測で、σn↓t∧τ\sigma_n\downarrow t\wedge\tau と道の連続性から Mt∧τ=lim⁡MσnM_{t\wedge\tau}=\lim M_{\sigma_n} も Ft\mathcal{F}_t-可測です。可積分性は 定理 6.2[マルチンゲールとブラウン運動] (a) の主張に含まれます。

マルチンゲール性: s<ts<t、A∈FsA\in\mathcal{F}_s とし、ρ:=τ∧s\rho:=\tau\wedge s(AA の上)、ρ:=τ∧t\rho:=\tau\wedge t(AcA^c の上)とおきます。r<sr<s なら {ρ≤r}={τ≤r}\{\rho\le r\}=\{\tau\le r\}、s≤r<ts\le r<t なら {ρ≤r}=A∪(Ac∩{τ≤r})\{\rho\le r\}=A\cup(A^c\cap\{\tau\le r\})、r≥tr\ge t なら {ρ≤r}=Ω\{\rho\le r\}=\Omega で、いずれも Fr\mathcal{F}_r に属するので ρ\rho は停止時刻です。定理 6.2[マルチンゲールとブラウン運動] (a) を停止時刻 ρ\rho と τ\tau に(定数 tt で)適用すると E[Mt∧ρ]=E[M0]=E[Mt∧τ]\mathbb{E}[M_{t\wedge\rho}]=\mathbb{E}[M_0]=\mathbb{E}[M_{t\wedge\tau}] です。t∧ρt\wedge\rho は AA の上で τ∧s\tau\wedge s、AcA^c の上で τ∧t\tau\wedge t なので、差を取ると E[1A(Mτ∧s−Mτ∧t)]=0\mathbb{E}[\mathbf{1}_A(M_{\tau\wedge s}-M_{\tau\wedge t})]=0 を得ます。

補題 6.2(特性関数による独立性と分布の判定)✓

YY を Rd\mathbb{R}^d 値確率変数、G\mathcal{G} を部分 σ-加法族、h≥0h\ge0 とし、すべての A∈GA\in\mathcal{G} と u∈Rdu\in\mathbb{R}^d について

E[1Aeiu⋅Y]=P(A) e−∣u∣2h/2\mathbb{E}\bigl[\mathbf{1}_Ae^{iu\cdot Y}\bigr]=\mathbb{P}(A)\,e^{-|u|^2h/2}

が成り立つとします。このとき YY は G\mathcal{G} と独立で、N(0,hId)N(0,hI_d) に従います。

証明(補題 6.2)

A=ΩA=\Omega とすると YY の特性関数は e−∣u∣2h/2e^{-|u|^2h/2} で、これは N(0,hId)N(0,hI_d) の特性関数です(補題 2.3[特性関数と Gauss ベクトル] を成分ごとに使い、独立な成分の積を取ります)。定理 3.3[特性関数と Gauss ベクトル] から Y∼N(0,hId)=:γhY\sim N(0,hI_d)=:\gamma_h です。P(A)>0\mathbb{P}(A)>0 なる A∈GA\in\mathcal{G} については、確率測度 P(⋅∣A)\mathbb{P}(\cdot\mid A) のもとでの YY の特性関数も仮定から e−∣u∣2h/2e^{-|u|^2h/2} なので、同じ定理により P(Y∈B∣A)=γh(B)\mathbb{P}(Y\in B\mid A)=\gamma_h(B) です。よって P(A∩{Y∈B})=P(A) P(Y∈B)\mathbb{P}(A\cap\{Y\in B\})=\mathbb{P}(A)\,\mathbb{P}(Y\in B) で、P(A)=0\mathbb{P}(A)=0 のときも両辺は 00 です。

次が証明の核心です。伊藤積分を使わず、eiλxe^{i\lambda x} の Taylor 展開とマルチンゲール差の直交性だけを使います。

証明の流れは次のとおりです。過程がまだ大きく動いていない間だけ追いかけ、増分が ±K\pm K の帯から出た時刻 τ\tau で止めると、止めた過程は常に [−K,K][-K,K] に収まります。有界になれば Taylor 展開の剰余を評価できるので、各小区間で 3 次の項を無視してよいことが示せ、止めた過程の特性関数が正規分布のものになります。最後に K→∞K\to\infty として停止を外します。

K−Kτst
Lévy の特徴付けの証明で使う停止。増分 N_r − N_s が ±K の帯を出た時刻 τ で止めると、止めた過程は有界になり、分割 r_0 < r_1 < … の各区間で 3 次の Taylor 展開が使える

補題 6.3(Lévy の特徴付けの 1 次元の核)✓

(Nt)t≥0(N_t)_{t\ge0} を道が連続な (Ft)(\mathcal{F}_t)-マルチンゲールで、各 NtN_t が二乗可積分、かつ Nt2−tN_t^2-t も (Ft)(\mathcal{F}_t)-マルチンゲールであるとします。このとき 0≤s<t0\le s<t、A∈FsA\in\mathcal{F}_s、λ∈R\lambda\in\mathbb{R} について

E[1A eiλ(Nt−Ns)]=P(A) e−λ2(t−s)/2.\mathbb{E}\bigl[\mathbf{1}_A\,e^{i\lambda(N_t-N_s)}\bigr]=\mathbb{P}(A)\,e^{-\lambda^2(t-s)/2}.
証明(補題 6.3)

停止による有界化。 K>0K>0 とし、τ:=inf⁡{r≥s:∣Nr−Ns∣≥K}∧t\tau:=\inf\{r\ge s:|N_r-N_s|\ge K\}\wedge t とおきます。道の連続性から、s≤r<ts\le r<t について {τ≤r}=⋂m⋃q{∣Nq−Ns∣>K−1/m}\{\tau\le r\}=\bigcap_{m}\bigcup_{q}\{|N_q-N_s|>K-1/m\}(qq は [s,r][s,r] の有理数と rr を動く)は Fr\mathcal{F}_r に属し、τ\tau は停止時刻です。r∈[s,t]r\in[s,t] について Xr:=Nr∧τ−NsX_r:=N_{r\wedge\tau}-N_s とおくと ∣Xr∣≤K|X_r|\le K です。補題 6.1 を NN と N2−tN^2-t に適用すると、s≤r≤r′≤ts\le r\le r'\le t について Δ:=Xr′−Xr\Delta:=X_{r'}-X_r、δ:=r′∧τ−r∧τ∈[0,r′−r]\delta:=r'\wedge\tau-r\wedge\tau\in[0,r'-r] は

E[Δ∣Fr]=0,E[Δ2−δ∣Fr]=0\mathbb{E}[\Delta\mid\mathcal{F}_r]=0,\qquad\mathbb{E}[\Delta^2-\delta\mid\mathcal{F}_r]=0

を満たします。後者は Δ2=Nr′∧τ2−Nr∧τ2−2Nr∧τΔ\Delta^2=N_{r'\wedge\tau}^2-N_{r\wedge\tau}^2-2N_{r\wedge\tau}\Delta と、Nr∧τ=Ns+XrN_{r\wedge\tau}=N_s+X_r が二乗可積分であることから 定理 4.4[条件付き期待値] で従います。

離散的な展開。 Φ(r):=exp⁡(iλXr+12λ2(r∧τ−s))\Phi(r):=\exp\bigl(i\lambda X_r+\tfrac12\lambda^2(r\wedge\tau-s)\bigr) とおくと、Φ(r)\Phi(r) は Fr\mathcal{F}_r-可測で ∣Φ(r)∣≤eλ2(t−s)/2|\Phi(r)|\le e^{\lambda^2(t-s)/2} です。分割 s=r0<r1<⋯<rm=ts=r_0<r_1<\dots<r_m=t を取り、Δk=Xrk+1−Xrk\Delta_k=X_{r_{k+1}}-X_{r_k}、δk=rk+1∧τ−rk∧τ\delta_k=r_{k+1}\wedge\tau-r_k\wedge\tau とおくと Φ(rk+1)=Φ(rk) g(Δk,δk)\Phi(r_{k+1})=\Phi(r_k)\,g(\Delta_k,\delta_k)、g(x,y)=eiλx+λ2y/2g(x,y)=e^{i\lambda x+\lambda^2y/2} です。w=iλx+12λ2yw=i\lambda x+\frac12\lambda^2y について ew=1+w+12w2+Re^w=1+w+\frac12w^2+R、∣R∣≤16∣w∣3e∣w∣|R|\le\frac16|w|^3e^{|w|} で、12w2=−12λ2x2+i2λ3xy+18λ4y2\frac12w^2=-\frac12\lambda^2x^2+\frac{i}{2}\lambda^3xy+\frac18\lambda^4y^2 ですから

g(x,y)=1+iλx+12λ2(y−x2)+ρ(x,y),∣ρ(x,y)∣≤c (∣x∣3+∣x∣y+y2)g(x,y)=1+i\lambda x+\tfrac12\lambda^2(y-x^2)+\rho(x,y),\qquad|\rho(x,y)|\le c\,\bigl(|x|^3+|x|y+y^2\bigr)

が ∣x∣≤2K|x|\le2K、0≤y≤t−s0\le y\le t-s の範囲で成り立ちます(cc は λ,K,t−s\lambda,K,t-s だけで決まる定数。∣w∣3≤4(∣λx∣3+18λ6y3)|w|^3\le4(|\lambda x|^3+\frac18\lambda^6y^3) と y3≤(t−s)y2y^3\le(t-s)y^2 を使います)。1AΦ(rk)\mathbf{1}_A\Phi(r_k) は有界で Frk\mathcal{F}_{r_k}-可測なので、上の二つの条件付き期待値から

E[1AΦ(rk+1)]−E[1AΦ(rk)]=E[1AΦ(rk) ρ(Δk,δk)].\mathbb{E}[\mathbf{1}_A\Phi(r_{k+1})]-\mathbb{E}[\mathbf{1}_A\Phi(r_k)]=\mathbb{E}\bigl[\mathbf{1}_A\Phi(r_k)\,\rho(\Delta_k,\delta_k)\bigr].

kk について足すと、Φ(r0)=1\Phi(r_0)=1 より

∣E[1AΦ(t)]−P(A)∣≤c eλ2(t−s)/2 E[∑k(∣Δk∣3+∣Δk∣δk+δk2)].\bigl|\mathbb{E}[\mathbf{1}_A\Phi(t)]-\mathbb{P}(A)\bigr|\le c\,e^{\lambda^2(t-s)/2}\,\mathbb{E}\Bigl[\sum_k\bigl(|\Delta_k|^3+|\Delta_k|\delta_k+\delta_k^2\bigr)\Bigr].

右辺が消えること。 ∥Π∥=max⁡k(rk+1−rk)\|\Pi\|=\max_k(r_{k+1}-r_k) とします。まず ∑kδk≤t−s\sum_k\delta_k\le t-s と E∑kΔk2=E∑kδk≤t−s\mathbb{E}\sum_k\Delta_k^2=\mathbb{E}\sum_k\delta_k\le t-s です。

  • ∑kδk2≤∥Π∥(t−s)\sum_k\delta_k^2\le\|\Pi\|(t-s)。
  • Cauchy–Schwarz の不等式から E∑k∣Δk∣δk≤(E∑kΔk2)1/2(E∑kδk2)1/2≤(t−s)∥Π∥1/2\mathbb{E}\sum_k|\Delta_k|\delta_k\le(\mathbb{E}\sum_k\Delta_k^2)^{1/2}(\mathbb{E}\sum_k\delta_k^2)^{1/2}\le(t-s)\|\Pi\|^{1/2}。
  • ∑k∣Δk∣3≤max⁡k∣Δk∣⋅Q\sum_k|\Delta_k|^3\le\max_k|\Delta_k|\cdot Q、Q:=∑kΔk2Q:=\sum_k\Delta_k^2 なので E∑k∣Δk∣3≤(Emax⁡kΔk2)1/2(EQ2)1/2\mathbb{E}\sum_k|\Delta_k|^3\le(\mathbb{E}\max_k\Delta_k^2)^{1/2}(\mathbb{E}Q^2)^{1/2}。ここで EQ2\mathbb{E}Q^2 は分割によらず有界です。実際 Q2=∑kΔk4+2∑kΔk2∑l>kΔl2Q^2=\sum_k\Delta_k^4+2\sum_k\Delta_k^2\sum_{l>k}\Delta_l^2 で、Δk4≤4K2Δk2\Delta_k^4\le4K^2\Delta_k^2、また l≠l′l\ne l' のとき E[ΔlΔl′∣Frk+1]=0\mathbb{E}[\Delta_l\Delta_{l'}\mid\mathcal{F}_{r_{k+1}}]=0(l,l′>kl,l'>k)なので E[∑l>kΔl2∣Frk+1]=E[(Xt−Xrk+1)2∣Frk+1]≤4K2\mathbb{E}[\sum_{l>k}\Delta_l^2\mid\mathcal{F}_{r_{k+1}}]=\mathbb{E}[(X_t-X_{r_{k+1}})^2\mid\mathcal{F}_{r_{k+1}}]\le4K^2 です。よって EQ2≤4K2(t−s)+8K2(t−s)=12K2(t−s)\mathbb{E}Q^2\le4K^2(t-s)+8K^2(t-s)=12K^2(t-s)。一方 max⁡kΔk2≤4K2\max_k\Delta_k^2\le4K^2 で、道 r↦Xrr\mapsto X_r の一様連続性から各 ω\omega で max⁡k∣Δk∣→0\max_k|\Delta_k|\to0 なので、優収束定理により Emax⁡kΔk2→0\mathbb{E}\max_k\Delta_k^2\to0。

したがって ∥Π∥→0\|\Pi\|\to0 とすると E[1AΦ(t)]=P(A)\mathbb{E}[\mathbf{1}_A\Phi(t)]=\mathbb{P}(A) です。

停止を外す。 KK を K→∞K\to\infty とします。sup⁡r∈[s,t]∣Nr−Ns∣<K\sup_{r\in[s,t]}|N_r-N_s|<K となる ω\omega では τ=t\tau=t なので、各 ω\omega で KK が十分大きければ Φ(t)=eiλ(Nt−Ns)+λ2(t−s)/2\Phi(t)=e^{i\lambda(N_t-N_s)+\lambda^2(t-s)/2} です。∣Φ(t)∣≤eλ2(t−s)/2|\Phi(t)|\le e^{\lambda^2(t-s)/2} なので、優収束定理から E[1Aeiλ(Nt−Ns)] eλ2(t−s)/2=P(A)\mathbb{E}[\mathbf{1}_Ae^{i\lambda(N_t-N_s)}]\,e^{\lambda^2(t-s)/2}=\mathbb{P}(A) を得ます。

定理 6.4(Lévy の特徴付け)✓

M=(M1,…,Md)M=(M^1,\dots,M^d) を (Ft)(\mathcal{F}_t) に適合した Rd\mathbb{R}^d 値の過程で、M0=0M_0=0、すべての道が連続、各 MiM^i が二乗可積分な (Ft)(\mathcal{F}_t)-マルチンゲールであり、さらにすべての i,ji,j について MtiMtj−δijtM^i_tM^j_t-\delta_{ij}t が (Ft)(\mathcal{F}_t)-マルチンゲールであるとします(すなわち ⟨Mi,Mj⟩t=δijt\langle M^i,M^j\rangle_t=\delta_{ij}t)。このとき MM は dd 次元 (Ft)(\mathcal{F}_t)-Brown 運動です。

証明(定理 6.4)

定義 2.1 の条件 (1)(2) は仮定です。条件 (3) を示します。u∈Rdu\in\mathbb{R}^d, u≠0u\ne0 を固定し、Nt:=u⋅Mt/∣u∣N_t:=u\cdot M_t/|u| とおきます。NN は道が連続な二乗可積分マルチンゲールで、

Nt2−t=1∣u∣2∑i,juiuj(MtiMtj−δijt)N_t^2-t=\frac{1}{|u|^2}\sum_{i,j}u_iu_j\bigl(M^i_tM^j_t-\delta_{ij}t\bigr)

はマルチンゲールの線形結合なのでマルチンゲールです。補題 6.3 を λ=∣u∣\lambda=|u| で適用すると、s<ts<t、A∈FsA\in\mathcal{F}_s について

E[1Aeiu⋅(Mt−Ms)]=P(A) e−∣u∣2(t−s)/2\mathbb{E}\bigl[\mathbf{1}_Ae^{iu\cdot(M_t-M_s)}\bigr]=\mathbb{P}(A)\,e^{-|u|^2(t-s)/2}

で、u=0u=0 でも自明に成り立ちます。補題 6.2 から Mt−MsM_t-M_s は Fs\mathcal{F}_s と独立で N(0,(t−s)Id)N(0,(t-s)I_d) に従います。

連続性の仮定は外せません。演習 7.4 で見るように、跳びのあるマルチンゲールでは Mt2−tM_t^2-t がマルチンゲールでも Brown 運動にはなりません。証明の中で連続性を使ったのは、停止した過程が有界になる箇所と、max⁡k∣Δk∣→0\max_k|\Delta_k|\to0 の箇所です。

例 6.5(符号で折り返した積分)

d=1d=1 とし、sgn⁡(x)=1\operatorname{sgn}(x)=1(x>0x>0)、−1-1(x≤0x\le0)として βt:=∫0tsgn⁡(Wu) dWu\beta_t:=\int_0^t\operatorname{sgn}(W_u)\,dW_u とおきます。sgn⁡(Wu)\operatorname{sgn}(W_u) は連続な適合過程のボレル関数なので発展的可測で、有界です。β\beta は連続なマルチンゲールで、命題 3.2 から βt2−∫0tsgn⁡(Wu)2 du=βt2−t\beta_t^2-\int_0^t\operatorname{sgn}(W_u)^2\,du=\beta_t^2-t もマルチンゲールです(任意の TT で成り立つので [0,∞)[0,\infty) 上で成り立ちます)。定理 6.4 により β\beta は Brown 運動です。β\beta は WW の符号を見て向きを変えた過程なので、見かけ上 WW と強く依存していますが、分布は Brown 運動そのものです。

例 6.6(ランダムな回転)

OuO_u を発展的可測な過程で、各 (u,ω)(u,\omega) で d×dd\times d の直交行列であるものとし、W~t:=∫0tOu dWu\tilde W_t:=\int_0^tO_u\,dW_u(成分は W~i=∑a∫Oia dWa\tilde W^i=\sum_a\int O^{ia}\,dW^a)とおきます。命題 3.2 から W~iW~j−∫0t(OuOuT)ij du=W~iW~j−δijt\tilde W^i\tilde W^j-\int_0^t(O_uO_u^{\mathsf T})^{ij}\,du=\tilde W^i\tilde W^j-\delta_{ij}t はマルチンゲールで、定理 6.4 により W~\tilde W は dd 次元 Brown 運動です。回転を時刻ごとに、過去の道に応じて選び直しても、Brown 運動の等方性は崩れません。

演習 7.1易

d=2d=2 とし、Wt1Wt2W^1_tW^2_t を 系 5.3 で確率積分として表し、それを使って E[(Wt1Wt2)2]=t2\mathbb{E}[(W^1_tW^2_t)^2]=t^2 を示してください。

解答

⟨W1,W2⟩=0\langle W^1,W^2\rangle=0 なので Wt1Wt2=∫0tWu2 dWu1+∫0tWu1 dWu2=∫0tHu⋅dWuW^1_tW^2_t=\int_0^tW^2_u\,dW^1_u+\int_0^tW^1_u\,dW^2_u=\int_0^tH_u\cdot dW_u、H=(W2,W1)H=(W^2,W^1) です。補題 3.1 により

E[(Wt1Wt2)2]=E∫0t((Wu2)2+(Wu1)2) du=∫0t2u du=t2.\mathbb{E}[(W^1_tW^2_t)^2]=\mathbb{E}\int_0^t\bigl((W^2_u)^2+(W^1_u)^2\bigr)\,du=\int_0^t2u\,du=t^2 .

独立性から直接 E[(Wt1)2] E[(Wt2)2]=t⋅t\mathbb{E}[(W^1_t)^2]\,\mathbb{E}[(W^2_t)^2]=t\cdot t と計算した値と一致します。

演習 7.2標準

a∈Ra\in\mathbb{R} とし、1 次元 Brown 運動 WW について m(t):=E[Wt2cos⁡(aWt)]m(t):=\mathbb{E}[W_t^2\cos(aW_t)] とおきます。g(x)=x2cos⁡(ax)g(x)=x^2\cos(ax) は導関数が有界でないので、定理 5.1 をそのままは使えません。χ∈C∞(R)\chi\in C^\infty(\mathbb{R}) を ∣x∣≤1|x|\le1 で 11、∣x∣≥2|x|\ge2 で 00、0≤χ≤10\le\chi\le1 なる関数とし、gN(x)=g(x)χ(x/N)g_N(x)=g(x)\chi(x/N) に伊藤公式を適用して N→∞N\to\infty とすることで

m′(t)=E[cos⁡(aWt)]−2a E[Wtsin⁡(aWt)]−a22m(t)m'(t)=\mathbb{E}[\cos(aW_t)]-2a\,\mathbb{E}[W_t\sin(aW_t)]-\frac{a^2}{2}m(t)

を導いてください。さらに m(t)=(t−a2t2)e−a2t/2m(t)=(t-a^2t^2)e^{-a^2t/2} がこれを満たすことを確かめてください。

解答

gNg_N は導関数が有界なので、定理 5.1(X=WX=W、μ=0\mu=0、σ=1\sigma=1)と確率積分の平均が 00 であることから E[gN(Wt)]=12∫0tE[gN′′(Wu)] du\mathbb{E}[g_N(W_t)]=\frac12\int_0^t\mathbb{E}[g_N''(W_u)]\,du です。gN′′=g′′χ(⋅/N)+2g′χ′(⋅/N)/N+gχ′′(⋅/N)/N2g_N''=g''\chi(\cdot/N)+2g'\chi'(\cdot/N)/N+g\chi''(\cdot/N)/N^2 で、χ′(⋅/N)\chi'(\cdot/N) と χ′′(⋅/N)\chi''(\cdot/N) は N≤∣x∣≤2NN\le|x|\le2N の外では 00 なので、∣g∣≤x2|g|\le x^2、∣g′∣≤2∣x∣+∣a∣x2|g'|\le2|x|+|a|x^2 から ∣gN∣+∣gN′′∣≤c(1+x2)|g_N|+|g_N''|\le c(1+x^2)(cc は NN によらない)です。E[Wu2]=u\mathbb{E}[W_u^2]=u なので、優収束定理により N→∞N\to\infty で m(t)=12∫0tE[g′′(Wu)] dum(t)=\frac12\int_0^t\mathbb{E}[g''(W_u)]\,du を得ます。g′′(x)=2cos⁡(ax)−4axsin⁡(ax)−a2x2cos⁡(ax)g''(x)=2\cos(ax)-4ax\sin(ax)-a^2x^2\cos(ax) を代入して tt で微分すると求める式です。

例 5.5 から E[cos⁡(aWt)]=e−a2t/2\mathbb{E}[\cos(aW_t)]=e^{-a^2t/2}、これを aa で微分して E[Wtsin⁡(aWt)]=at e−a2t/2\mathbb{E}[W_t\sin(aW_t)]=at\,e^{-a^2t/2} です(微分と期待値の交換は優収束定理による)。よって m′(t)=(1−2a2t)e−a2t/2−a22m(t)m'(t)=(1-2a^2t)e^{-a^2t/2}-\frac{a^2}{2}m(t) で、m(t)=(t−a2t2)e−a2t/2m(t)=(t-a^2t^2)e^{-a^2t/2} を代入すると左辺は (1−2a2t)e−a2t/2−a22(t−a2t2)e−a2t/2(1-2a^2t)e^{-a^2t/2}-\frac{a^2}{2}(t-a^2t^2)e^{-a^2t/2} となり右辺に一致します。m(0)=0m(0)=0 なので、1 階線形常微分方程式の解の一意性からこれが mm です。検算として、m(t)=−d2da2e−a2t/2m(t)=-\frac{d^2}{da^2}e^{-a^2t/2} からも同じ式が出ます。

演習 7.3標準

d=2d=2 とし、θu:=Wu1\theta_u:=W^1_u として

Bt=∫0t(cos⁡θu dWu1+sin⁡θu dWu2),Bt′=∫0t(−sin⁡θu dWu1+cos⁡θu dWu2)B_t=\int_0^t\bigl(\cos\theta_u\,dW^1_u+\sin\theta_u\,dW^2_u\bigr),\qquad B'_t=\int_0^t\bigl(-\sin\theta_u\,dW^1_u+\cos\theta_u\,dW^2_u\bigr)

とおきます。(B,B′)(B,B') が 2 次元 Brown 運動であることを示してください。

解答

Ou=(cos⁡θusin⁡θu−sin⁡θucos⁡θu)O_u=\begin{pmatrix}\cos\theta_u&\sin\theta_u\\-\sin\theta_u&\cos\theta_u\end{pmatrix} は各 (u,ω)(u,\omega) で直交行列で、θu=Wu1\theta_u=W^1_u が連続な適合過程なので発展的可測です。(B,B′)=∫O dW(B,B')=\int O\,dW なので 例 6.6 から従います。直接には、命題 3.2 により B2−tB^2-t、B′2−tB'^2-t(cos⁡2+sin⁡2=1\cos^2+\sin^2=1)と BB′BB'(−cos⁡θsin⁡θ+sin⁡θcos⁡θ=0-\cos\theta\sin\theta+\sin\theta\cos\theta=0)がマルチンゲールで、定理 6.4 を使います。

演習 7.4難

(Pt)(P_t) を強度 11 の Poisson 過程とし、その自然なフィルトレーションを (Ft)(\mathcal{F}_t) とします。Mt:=Pt−tM_t:=P_t-t について、MtM_t と Mt2−tM_t^2-t がどちらも (Ft)(\mathcal{F}_t)-マルチンゲールであることを示し、それでも MM が Brown 運動でない理由を述べてください。定理 6.4 のどの仮定が破れていますか。

解答

定理 4.5[Markov 連鎖と Poisson 過程] により Pt−PsP_t-P_s は Fs\mathcal{F}_s と独立で平均 t−st-s、分散 t−st-s の Poisson 分布に従います。よって E[Mt−Ms∣Fs]=0\mathbb{E}[M_t-M_s\mid\mathcal{F}_s]=0、また

E[Mt2−Ms2∣Fs]=E[(Mt−Ms)2∣Fs]+2Ms E[Mt−Ms∣Fs]=t−s\mathbb{E}[M_t^2-M_s^2\mid\mathcal{F}_s]=\mathbb{E}[(M_t-M_s)^2\mid\mathcal{F}_s]+2M_s\,\mathbb{E}[M_t-M_s\mid\mathcal{F}_s]=t-s

なので Mt2−tM_t^2-t もマルチンゲールです。しかし Mt+t=PtM_t+t=P_t は整数値なので MtM_t は正規分布に従いません。破れているのは道の連続性の仮定で、PP は高さ 11 の跳びを持ちます。補題 6.3 の証明では max⁡k∣Δk∣→0\max_k|\Delta_k|\to0 を使いましたが、跳びがあるとこの量は 11 以上のまま残り、Taylor 展開の 3 次の剰余が消えません。

  • I. Karatzas, S. E. Shreve, Brownian Motion and Stochastic Calculus, 2nd ed., Springer, 1991 — 第 3.2 節(確率積分と 2 次変分)、第 3.3 節(伊藤の公式と Lévy の特徴付け)。
  • D. Revuz, M. Yor, Continuous Martingales and Brownian Motion, 3rd ed., Springer, 1999 — 第 IV 章(確率積分・伊藤の公式)、第 IV.3 節(Lévy の特徴付け)。
  • J.-F. Le Gall, Brownian Motion, Martingales, and Stochastic Calculus, Springer, 2016 — 第 5 章(伊藤の公式とその応用)。
  • B. Øksendal, Stochastic Differential Equations: An Introduction with Applications, 6th ed., Springer, 2003 — 第 4.2 節(多次元伊藤公式)、第 8.6 節(Lévy の特徴付けの応用)。
  • 舟木直久『確率微分方程式』岩波書店、2005 — 伊藤の公式と確率積分の章。

このテーマの 100 本ノック(100 問)

執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12

この記事の誤りを報告する ・運営: 夢現技研合同会社 ・料金 ・利用条件 ・特定商取引法に基づく表記

© 2026 夢現技研合同会社 ・コード例は MIT ライセンスです。