複数の株価や、3 次元空間を揺れ動く粒子のように、ランダムな揺れが同時にいくつも働く状況を考えます。このとき「二つの量の積は、微小時間にどれだけ変化するか」「ある過程が、ランダムな歩みの基本形である Brown 運動かどうかを、分布を計算せずに見分けられるか」という素朴な疑問が出ます。この記事は、前者に答える多次元の伊藤公式と積の公式、後者に答える Lévy の判定条件を、証明つきで示します。読み終えると、複数成分の確率過程の微小変化を計算でき、与えられた過程が Brown 運動であることを確かめられるようになります。
d d d 次元 Brown 運動 W = ( W 1 , … , W d ) W=(W^1,\dots,W^d) W = ( W 1 , … , W d ) の成分は互いに独立なので、異なる成分による伊藤積分は直交します。これを条件付きの形で述べたのが共分散等長性 E [ ξ ∫ H d W a ∫ K d W b ] = δ a b E [ ξ ∫ H K d u ] \mathbb{E}[\xi\int H\,dW^a\int K\,dW^b]=\delta_{ab}\mathbb{E}[\xi\int HK\,du] E [ ξ ∫ H d W a ∫ K d W b ] = δ ab E [ ξ ∫ H K d u ] で、多次元の計算規則 d W a d W b = δ a b d t dW^a\,dW^b=\delta_{ab}\,dt d W a d W b = δ ab d t の正体です。
伊藤過程 d X i = μ i d t + ∑ a σ i a d W a dX^i=\mu^i\,dt+\sum_a\sigma^{ia}\,dW^a d X i = μ i d t + ∑ a σ ia d W a の共変分を ⟨ X i , X j ⟩ t = ∫ 0 t ( σ σ T ) i j d u \langle X^i,X^j\rangle_t=\int_0^t(\sigma\sigma^{\mathsf T})^{ij}\,du ⟨ X i , X j ⟩ t = ∫ 0 t ( σ σ T ) ij d u と定めると、増分の積の和 ∑ k Δ X k i Δ X k j \sum_k\Delta X^i_k\Delta X^j_k ∑ k Δ X k i Δ X k j はこれに L 1 L^1 L 1 で収束します。証明の要は伊藤積分の 4 次モーメント評価です。
多次元伊藤公式は、2 次の項 1 2 ∑ i , j ∂ i j f d ⟨ X i , X j ⟩ \frac12\sum_{i,j}\partial_{ij}f\,d\langle X^i,X^j\rangle 2 1 ∑ i , j ∂ ij f d ⟨ X i , X j ⟩ を持つ連鎖律です。特に f ( x , y ) = x y f(x,y)=xy f ( x , y ) = x y から積の公式 d ( X Y ) = X d Y + Y d X + d ⟨ X , Y ⟩ d(XY)=X\,dY+Y\,dX+d\langle X,Y\rangle d ( X Y ) = X d Y + Y d X + d ⟨ X , Y ⟩ が出ます。
Lévy の特徴付け: 連続なマルチンゲール M M M で M i M j − δ i j t M^iM^j-\delta_{ij}t M i M j − δ ij t もマルチンゲールなら、M M M は Brown 運動です。証明は e i λ x e^{i\lambda x} e iλ x の 3 次までの Taylor 展開と停止による有界化だけで済みます。
応用として、∫ 0 t sgn ( W u ) d W u \int_0^t\operatorname{sgn}(W_u)\,dW_u ∫ 0 t sgn ( W u ) d W u や直交行列値の過程による回転 ∫ O d W \int O\,dW ∫ O d W が再び Brown 運動であることがすぐに分かります。
この記事で使う事前知識
Brown 運動 : 増分が過去と独立で、時間幅 t − s t-s t − s に応じた分散の正規分布に従う、連続な道を持つ過程。→ ブラウン運動
マルチンゲール : 現在までの情報のもとでの将来の期待値が、現在の値に等しい過程。→ マルチンゲール
条件付き期待値 : ある情報が分かっているという前提のもとでの期待値。→ 条件付き期待値
停止時刻 : 「時刻 r r r までに起きたか」が時刻 r r r までの情報で判定できる、ランダムな時刻。→ 停止時刻
特性関数 : 確率変数 Y Y Y に対する E [ e i u ⋅ Y ] \mathbb{E}[e^{iu\cdot Y}] E [ e i u ⋅ Y ] のこと。分布を決めます。→ 特性関数
確率収束 : どんな ε > 0 \varepsilon>0 ε > 0 についても、極限との差が ε \varepsilon ε を超える確率が 0 0 0 に近づくこと。→ 確率収束
確率微分方程式の章 では、1 本の Brown 運動 B B B に関する伊藤積分と伊藤の公式を扱いました。しかし応用で現れる確率過程の多くは多次元です。R 3 \mathbb{R}^3 R 3 の中を拡散する粒子、複数の株価、場の理論で多数のモードを同時に動かす確率量子化など、いずれも複数の独立な雑音で駆動されます。
多次元に移ると、新しい問いが二つ生じます。第一に、二つの成分の積 X 1 X 2 X^1X^2 X 1 X 2 はどう微分されるのか。1 次元の公式 ( d B ) 2 = d t (dB)^2=dt ( d B ) 2 = d t の多次元版は d W a d W b = δ a b d t dW^a\,dW^b=\delta_{ab}\,dt d W a d W b = δ ab d t ですが、これは「異なる成分の増分の積の和は 0 0 0 に収束する」という主張で、証明が要ります。第二に、ある連続マルチンゲールが Brown 運動であることを、分布を直接計算せずに判定できないか。Lévy の特徴付けはこれに「2 次変分が t t t なら Brown 運動」と答えます。その証明で使う「条件付き特性関数による独立性と分布の判定」は、次章の Girsanov の定理 の証明でもそのまま使います。また Lévy の特徴付けは、確率積分による時間変更の議論の土台になります。
この章では、確率微分方程式の章で構成した伊藤積分(定理 4.7[確率微分方程式] )だけを道具にして、これらを証明します。一般の連続セミマルチンゲールの理論(Doob–Meyer 分解など)には立ち入らず、係数が有界な伊藤過程に話を絞ります。この範囲で、Girsanov の定理とマルチンゲール表現定理に必要なものはすべてそろいます。
時間区間を [ 0 , T ] [0,T] [ 0 , T ] (T < ∞ T<\infty T < ∞ )とし、確率空間 ( Ω , F , P ) (\Omega,\mathcal{F},\mathbb{P}) ( Ω , F , P ) とフィルトレーション ( F t ) (\mathcal{F}_t) ( F t ) を固定します。発展的可測性などの用語は 定義 2.1[確率微分方程式] に従います。フィルトレーションとは時刻ごとに増えていく「その時点までに分かっている情報」の列で、発展的可測とは、各時刻 r r r について、r r r までの時間と r r r までの情報だけで過程の値が決まる、という可測性のことです。
定義 2.1 (d 次元 ( F t ) (F_t) ( F t ) -Brown 運動 ) ↗
R d \mathbb{R}^d R d 値の過程 W = ( W 1 , … , W d ) W=(W^1,\dots,W^d) W = ( W 1 , … , W d ) が ( F t ) (\mathcal{F}_t) ( F t ) -Brown 運動 であるとは、(1) W 0 = 0 W_0=0 W 0 = 0 ですべての道が連続、(2) W W W は ( F t ) (\mathcal{F}_t) ( F t ) に適合、(3) s < t s<t s < t のとき増分 W t − W s W_t-W_s W t − W s は F s \mathcal{F}_s F s と独立で、平均 0 0 0 、共分散行列 ( t − s ) I d (t-s)I_d ( t − s ) I d の Gauss 分布 N ( 0 , ( t − s ) I d ) N(0,(t-s)I_d) N ( 0 , ( t − s ) I d ) に従う、の三条件を満たすことをいいます。
共分散行列が対角なので、増分の成分 W t 1 − W s 1 , … , W t d − W s d W^1_t-W^1_s,\dots,W^d_t-W^d_s W t 1 − W s 1 , … , W t d − W s d は互いに独立です(命題 5.3[特性関数と Gauss ベクトル] )。特に各成分 W a W^a W a は 1 次元の ( F t ) (\mathcal{F}_t) ( F t ) -Brown 運動で、確率微分方程式の章の伊藤積分 ∫ H d W a \int H\,dW^a ∫ H d W a がそのまま使えます。∫ s t H d W a \int_s^tH\,dW^a ∫ s t H d W a は ∫ 0 t H d W a − ∫ 0 s H d W a \int_0^tH\,dW^a-\int_0^sH\,dW^a ∫ 0 t H d W a − ∫ 0 s H d W a の意味です。
定義 2.2 (多次元伊藤過程と有界係数の条件 ) ↗
L 2 [ 0 , T ] \mathcal{L}^2[0,T] L 2 [ 0 , T ] を 定義 4.1[確率微分方程式] の空間とします。R d \mathbb{R}^d R d 値の過程 H = ( H 1 , … , H d ) H=(H^1,\dots,H^d) H = ( H 1 , … , H d ) で各成分が L 2 [ 0 , T ] \mathcal{L}^2[0,T] L 2 [ 0 , T ] に属するもの全体を L d 2 \mathcal{L}^2_d L d 2 と書き、
∫ 0 t H u ⋅ d W u : = ∑ a = 1 d ∫ 0 t H u a d W u a \int_0^t H_u\cdot dW_u:=\sum_{a=1}^d\int_0^t H^a_u\,dW^a_u ∫ 0 t H u ⋅ d W u := a = 1 ∑ d ∫ 0 t H u a d W u a とおきます。R n \mathbb{R}^n R n 値の過程 X X X が
X t i = X 0 i + ∫ 0 t μ u i d u + ∑ a = 1 d ∫ 0 t σ u i a d W u a ( 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) X t i = X 0 i + ∫ 0 t μ u i d u + a = 1 ∑ d ∫ 0 t σ u ia d W u a ( i = 1 , … , n ) と書けるとき(X 0 X_0 X 0 は F 0 \mathcal{F}_0 F 0 -可測、μ = ( μ i ) \mu=(\mu^i) μ = ( μ i ) と σ = ( σ i a ) \sigma=(\sigma^{ia}) σ = ( σ ia ) は発展的可測で σ i a ∈ L 2 [ 0 , T ] \sigma^{ia}\in\mathcal{L}^2[0,T] σ ia ∈ L 2 [ 0 , T ] 、∫ 0 T ∣ μ u ∣ d u < ∞ \int_0^T|\mu_u|\,du<\infty ∫ 0 T ∣ μ u ∣ d u < ∞ )、X X X を n n n 次元の伊藤過程 と呼びます。σ \sigma σ の第 i i i 行を σ i = ( σ i 1 , … , σ i d ) \sigma^i=(\sigma^{i1},\dots,\sigma^{id}) σ i = ( σ i 1 , … , σ i d ) と書きます。さらに定数 K K K があって、すべての ( u , ω ) (u,\omega) ( u , ω ) と i , a i,a i , a について ∣ μ u i ∣ ≤ K |\mu^i_u|\le K ∣ μ u i ∣ ≤ K 、∣ σ u i a ∣ ≤ K |\sigma^{ia}_u|\le K ∣ σ u ia ∣ ≤ K となるとき、X X X は有界係数の条件 (A) を満たすといいます。
確率積分は連続なバージョンを取ります(定理 4.7[確率微分方程式] )。以下で繰り返し使う二つの補題を先に示します。
補題 2.3 (有界な単純過程による近似 ) ✓
H H H を R d \mathbb{R}^d R d 値の発展的可測過程で、定数 L L L について ∣ H u ( ω ) ∣ ≤ L |H_u(\omega)|\le L ∣ H u ( ω ) ∣ ≤ L (Euclid ノルム)とします。このとき、[ 0 , T ] [0,T] [ 0 , T ] の等分割の上の単純過程(各成分が 定義 4.1[確率微分方程式] の意味で単純)H ( n ) H^{(n)} H ( n ) で、∣ H ( n ) ∣ ≤ L |H^{(n)}|\le L ∣ H ( n ) ∣ ≤ L かつ E ∫ 0 T ∣ H u ( n ) − H u ∣ 2 d u → 0 \mathbb{E}\int_0^T|H^{(n)}_u-H_u|^2\,du\to0 E ∫ 0 T ∣ H u ( n ) − H u ∣ 2 d u → 0 となるものが存在します。
証明(補題 2.3) 補題 4.4[確率微分方程式] の証明の第 2 段階と第 3 段階を、成分ごとに同じ分割で行えば得られます。確かめるべきは上界 L L L が保たれることだけです。第 2 段階の平均化 H t ( m ) = m ∫ t − 1 / m t H s d s H^{(m)}_t=m\int_{t-1/m}^tH_s\,ds H t ( m ) = m ∫ t − 1/ m t H s d s では、ベクトル値積分のノルムが被積分関数のノルムの積分以下であることから ∣ H t ( m ) ∣ ≤ m ∫ t − 1 / m t ∣ H s ∣ d s ≤ L |H^{(m)}_t|\le m\int_{t-1/m}^t|H_s|\,ds\le L ∣ H t ( m ) ∣ ≤ m ∫ t − 1/ m t ∣ H s ∣ d s ≤ L です。第 3 段階の左端離散化は H H H の値をそのまま使うので上界を変えません。第 1 段階(切り落とし)は H H H がもともと有界なので不要です。各段階の収束は同補題の証明と同じく優収束定理によります。
∎
補題 2.4 (確率積分からの既知量の取り出し ) ✓
H ∈ L 2 [ 0 , T ] H\in\mathcal{L}^2[0,T] H ∈ L 2 [ 0 , T ] 、0 ≤ s < t ≤ T 0\le s<t\le T 0 ≤ s < t ≤ T 、ξ \xi ξ を有界な F s \mathcal{F}_s F s -可測確率変数とします。このとき ξ H 1 ( s , t ] ∈ L 2 [ 0 , T ] \xi H\mathbf{1}_{(s,t]}\in\mathcal{L}^2[0,T] ξ H 1 ( s , t ] ∈ L 2 [ 0 , T ] で、各 a a a について
∫ s t ξ H u d W u a = ξ ∫ s t H u d W u a a.s. \int_s^t\xi H_u\,dW^a_u=\xi\int_s^tH_u\,dW^a_u\quad\text{a.s.} ∫ s t ξ H u d W u a = ξ ∫ s t H u d W u a a.s.
証明(補題 2.4) 過程 ( u , ω ) ↦ ξ ( ω ) 1 ( s , t ] ( u ) (u,\omega)\mapsto\xi(\omega)\mathbf{1}_{(s,t]}(u) ( u , ω ) ↦ ξ ( ω ) 1 ( s , t ] ( u ) は、r ≤ s r\le s r ≤ s なら [ 0 , r ] [0,r] [ 0 , r ] 上で 0 0 0 、r > s r>s r > s なら B ( [ 0 , r ] ) ⊗ F s \mathcal{B}([0,r])\otimes\mathcal{F}_s B ([ 0 , r ]) ⊗ F s -可測なので発展的可測です。発展的可測な過程の積は発展的可測で、∣ ξ H 1 ( s , t ] ∣ ≤ ∥ ξ ∥ ∞ ∣ H ∣ |\xi H\mathbf{1}_{(s,t]}|\le\|\xi\|_\infty|H| ∣ ξ H 1 ( s , t ] ∣ ≤ ∥ ξ ∥ ∞ ∣ H ∣ なので ξ H 1 ( s , t ] ∈ L 2 [ 0 , T ] \xi H\mathbf{1}_{(s,t]}\in\mathcal{L}^2[0,T] ξ H 1 ( s , t ] ∈ L 2 [ 0 , T ] です。
H H H が単純過程で、その分割点に s , t s,t s , t が含まれる場合(分割を細かくしても単純過程のままなので一般性を失いません)、H 1 ( s , t ] = ∑ t i ≥ s ξ i 1 ( t i , t i + 1 ] H\mathbf{1}_{(s,t]}=\sum_{t_i\ge s}\xi_i\mathbf{1}_{(t_i,t_{i+1}]} H 1 ( s , t ] = ∑ t i ≥ s ξ i 1 ( t i , t i + 1 ] (和は t i + 1 ≤ t t_{i+1}\le t t i + 1 ≤ t の範囲)で、ξ ξ i \xi\xi_i ξ ξ i は t i ≥ s t_i\ge s t i ≥ s より F t i \mathcal{F}_{t_i} F t i -可測かつ二乗可積分です。よって 定義 4.2[確率微分方程式] から両辺はともに ∑ ξ ξ i ( W t i + 1 a − W t i a ) \sum\xi\xi_i(W^a_{t_{i+1}}-W^a_{t_i}) ∑ ξ ξ i ( W t i + 1 a − W t i a ) です。一般の H H H は 補題 4.4[確率微分方程式] の単純過程 H ( n ) → H H^{(n)}\to H H ( n ) → H で近似します。∥ ξ H ( n ) 1 ( s , t ] − ξ H 1 ( s , t ] ∥ L 2 ≤ ∥ ξ ∥ ∞ ∥ H ( n ) − H ∥ L 2 → 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 ∥ ξ H ( n ) 1 ( s , t ] − ξ H 1 ( s , t ] ∥ L 2 ≤ ∥ ξ ∥ ∞ ∥ H ( n ) − H ∥ L 2 → 0 なので、伊藤等長性(定理 4.7[確率微分方程式] )により左辺も右辺も L 2 ( Ω ) L^2(\Omega) L 2 ( Ω ) で収束し、極限でも等式が成り立ちます。
∎
この節のねらいは、異なる成分の Brown 運動による積分どうしの積の期待値が 0 0 0 になることを示し、多次元の計算規則 d W a d W b = δ a b d t dW^a\,dW^b=\delta_{ab}\,dt d W a d W b = δ ab d t の中身をはっきりさせることです。証明は、単純過程(区間ごとに値が一定の過程)で確かめてから一般の場合へ近づける、という伊藤積分の定義と同じ流れで進みます。途中の「偏極」は、x y = 1 4 ( ( x + y ) 2 − ( x − y ) 2 ) xy=\frac14\bigl((x+y)^2-(x-y)^2\bigr) x y = 4 1 ( ( x + y ) 2 − ( x − y ) 2 ) という恒等式を使って、2 乗についての等長性の式から積についての式を取り出す手順です。
補題 3.1 (条件付き共分散等長性 ) ✓
H , K ∈ L 2 [ 0 , T ] H,K\in\mathcal{L}^2[0,T] H , K ∈ L 2 [ 0 , T ] 、a , b ∈ { 1 , … , d } a,b\in\{1,\dots,d\} a , b ∈ { 1 , … , d } 、0 ≤ s < t ≤ T 0\le s<t\le T 0 ≤ s < t ≤ T 、ξ \xi ξ を有界な F s \mathcal{F}_s F s -可測確率変数とします。このとき
E [ ξ ∫ s t H u d W u a ∫ s t K u d W u b ] = δ a b E [ ξ ∫ s t H u K u d u ] . \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]. E [ ξ ∫ s t H u d W u a ∫ s t K u d W u b ] = δ ab E [ ξ ∫ s t H u K u d u ] . 特に H , G ∈ L d 2 H,G\in\mathcal{L}^2_d H , G ∈ L d 2 について E [ ∫ 0 T H ⋅ d W ∫ 0 T G ⋅ d W ] = E ∫ 0 T H u ⋅ G u d u \mathbb{E}[\int_0^TH\cdot dW\int_0^TG\cdot dW]=\mathbb{E}\int_0^TH_u\cdot G_u\,du E [ ∫ 0 T H ⋅ d W ∫ 0 T G ⋅ d W ] = E ∫ 0 T H u ⋅ G u d u です。
証明(補題 3.1) 補題 2.4 により ξ ∫ s t H d W a = ∫ s t ξ H d W a \xi\int_s^tH\,dW^a=\int_s^t\xi H\,dW^a ξ ∫ s t H d W a = ∫ s t ξ H d W a なので、ξ = 1 \xi=1 ξ = 1 の場合に、H H H を ξ H \xi H ξ H に置き換えて示せば十分です。また H , K H,K H , K を H 1 ( s , t ] , K 1 ( s , t ] H\mathbf{1}_{(s,t]},K\mathbf{1}_{(s,t]} H 1 ( s , t ] , K 1 ( s , t ] に置き換えます。
a = b a=b a = b のときは、伊藤等長性(定理 4.7[確率微分方程式] )を H + K H+K H + K と H − K H-K H − K に適用して差を取る(偏極)と E [ I ( H ) I ( K ) ] = E ∫ H K d u \mathbb{E}[I(H)I(K)]=\mathbb{E}\int HK\,du E [ I ( H ) I ( K )] = E ∫ H K d u が出ます。
a ≠ b a\ne b a = b のときを示します。まず H , K H,K H , K が共通の分割 s = u 0 < ⋯ < u m = t s=u_0<\dots<u_m=t s = u 0 < ⋯ < u m = t 上の単純過程 H = ∑ l ξ l 1 ( u l , u l + 1 ] H=\sum_l\xi_l\mathbf{1}_{(u_l,u_{l+1}]} H = ∑ l ξ l 1 ( u l , u l + 1 ] 、K = ∑ l η l 1 ( u l , u l + 1 ] K=\sum_l\eta_l\mathbf{1}_{(u_l,u_{l+1}]} K = ∑ l η l 1 ( u l , u l + 1 ] の場合です。Δ l W a = W u l + 1 a − W u l a \Delta_lW^a=W^a_{u_{l+1}}-W^a_{u_l} Δ l W a = W u l + 1 a − W u l a と書くと、示すべきは E [ ∑ l , m ξ l Δ l W a η m Δ m W b ] = 0 \mathbb{E}[\sum_{l,m}\xi_l\Delta_lW^a\,\eta_m\Delta_mW^b]=0 E [ ∑ l , m ξ l Δ l W a η m Δ m W b ] = 0 です。l < m l<m l < m の項では ξ l Δ l W a η m \xi_l\Delta_lW^a\eta_m ξ l Δ l W a η m が F u m \mathcal{F}_{u_m} F u m -可測で可積分、Δ m W b \Delta_mW^b Δ m W b は F u m \mathcal{F}_{u_m} F u m と独立で平均 0 0 0 なので、定理 4.4[条件付き期待値] と 命題 4.5[条件付き期待値] から期待値は E [ ξ l Δ l W a η m E [ Δ m W b ] ] = 0 \mathbb{E}[\xi_l\Delta_lW^a\eta_m\,\mathbb{E}[\Delta_mW^b]]=0 E [ ξ l Δ l W a η m E [ Δ m W b ]] = 0 です。l > m l>m l > m も同様です。l = m l=m l = m の項では、( Δ l W a , Δ l W b ) (\Delta_lW^a,\Delta_lW^b) ( Δ l W a , Δ l W b ) が F u l \mathcal{F}_{u_l} F u l と独立で、成分どうしも独立で平均 0 0 0 なので、E [ ξ l η l Δ l W a Δ l W b ] = E [ ξ l η l ] E [ Δ l W a ] E [ Δ l W b ] = 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 E [ ξ l η l Δ l W a Δ l W b ] = E [ ξ l η l ] E [ Δ l W a ] E [ Δ l W b ] = 0 です(命題 7.5[確率変数と期待値] )。
一般の H , K H,K H , K は単純過程 H ( n ) → H H^{(n)}\to H H ( n ) → H 、K ( n ) → K K^{(n)}\to K K ( n ) → K で近似します(補題 4.4[確率微分方程式] 。分割は両者の和集合に細分して共通にできます)。Cauchy–Schwarz の不等式と等長性から
∣ E [ I a ( H ) I b ( K ) ] − E [ I a ( H ( n ) ) I b ( K ( n ) ) ] ∣ ≤ ∥ H − H ( n ) ∥ L 2 ∥ K ∥ L 2 + ∥ H ( n ) ∥ L 2 ∥ K − K ( n ) ∥ L 2 → 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 E [ I a ( H ) I b ( K )] − E [ I a ( H ( n ) ) I b ( K ( n ) )] ≤ ∥ H − H ( n ) ∥ L 2 ∥ K ∥ L 2 + ∥ H ( n ) ∥ L 2 ∥ K − K ( n ) ∥ L 2 → 0 (I a ( H ) = ∫ H d W a I^a(H)=\int H\,dW^a I a ( H ) = ∫ H d W a )なので、0 0 0 という値が極限に受け継がれます。最後の主張は、成分ごとの式を a , b a,b a , b について足し合わせたものです。
∎
この補題を形式的に書けば d W u a d W u b = δ a b d u dW^a_u\,dW^b_u=\delta_{ab}\,du d W u a d W u b = δ ab d u です。条件付きの形(ξ \xi ξ が入っている形)で述べたのは、次の命題と 2 次変分の収束で、マルチンゲール差の直交性として使うためです。
命題 3.2 (確率積分の積から作られるマルチンゲール ) ✓
G , H ∈ L d 2 G,H\in\mathcal{L}^2_d G , H ∈ L d 2 とし、M t = ∫ 0 t G ⋅ d W M_t=\int_0^tG\cdot dW M t = ∫ 0 t G ⋅ d W 、N t = ∫ 0 t H ⋅ d W N_t=\int_0^tH\cdot dW N t = ∫ 0 t H ⋅ d W とおきます。このとき
M t N t − ∫ 0 t G u ⋅ H u d u ( t ∈ [ 0 , T ] ) M_tN_t-\int_0^tG_u\cdot H_u\,du\qquad(t\in[0,T]) M t N t − ∫ 0 t G u ⋅ H u d u ( t ∈ [ 0 , T ]) は ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲールです。
証明(命題 3.2) M t , N t ∈ L 2 ( Ω ) M_t,N_t\in L^2(\Omega) M t , N t ∈ L 2 ( Ω ) なので積は可積分で、∫ 0 T ∣ G ⋅ H ∣ d u \int_0^T|G\cdot H|\,du ∫ 0 T ∣ G ⋅ H ∣ d u も Cauchy–Schwarz の不等式から可積分です。s < t s<t s < t 、A ∈ F s A\in\mathcal{F}_s A ∈ F s とし、Δ M = M t − M s \Delta M=M_t-M_s Δ M = M t − M s 、Δ N = N t − N s \Delta N=N_t-N_s Δ N = N t − N s と書くと
M t N t − M s N s = Δ M Δ N + M s Δ N + N s Δ M . M_tN_t-M_sN_s=\Delta M\,\Delta N+M_s\,\Delta N+N_s\,\Delta M. M t N t − M s N s = Δ M Δ N + M s Δ N + N s Δ M . Δ N \Delta N Δ N は 定理 4.7[確率微分方程式] のマルチンゲール性から E [ Δ N ∣ F s ] = 0 \mathbb{E}[\Delta N\mid\mathcal{F}_s]=0 E [ Δ N ∣ F s ] = 0 を満たし、1 A M s \mathbf{1}_AM_s 1 A M s は二乗可積分で F s \mathcal{F}_s F s -可測なので、定理 4.4[条件付き期待値] により E [ 1 A M s Δ N ] = 0 \mathbb{E}[\mathbf{1}_AM_s\Delta N]=0 E [ 1 A M s Δ N ] = 0 です。第 3 項も同じです。第 1 項は Δ M Δ N = ∑ a , b ∫ s t G a d W a ∫ s t H b d W b \Delta M\,\Delta N=\sum_{a,b}\int_s^tG^a\,dW^a\int_s^tH^b\,dW^b Δ M Δ N = ∑ a , b ∫ s t G a d W a ∫ s t H b d W b と展開し、補題 3.1 (ξ = 1 A \xi=\mathbf{1}_A ξ = 1 A )を使うと E [ 1 A Δ M Δ N ] = E [ 1 A ∫ s t G ⋅ H d u ] \mathbb{E}[\mathbf{1}_A\Delta M\Delta N]=\mathbb{E}[\mathbf{1}_A\int_s^tG\cdot H\,du] E [ 1 A Δ M Δ N ] = E [ 1 A ∫ s t G ⋅ H d u ] です。合わせて E [ 1 A ( M t N t − ∫ 0 t G ⋅ H ) ] = E [ 1 A ( M s N s − ∫ 0 s G ⋅ 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)] E [ 1 A ( M t N t − ∫ 0 t G ⋅ H )] = E [ 1 A ( M s N s − ∫ 0 s G ⋅ H )] を得ます。
∎
この命題が、次の定義を正当化します。
定義 3.3 (伊藤過程の共変分 ) ↗
定義 2.2 の伊藤過程 X X X について
⟨ X i , X j ⟩ t : = ∫ 0 t σ u i ⋅ σ u j d u = ∫ 0 t ( σ u σ u T ) i j d u \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 ⟨ X i , X j ⟩ t := ∫ 0 t σ u i ⋅ σ u j d u = ∫ 0 t ( σ u σ u T ) ij d u を X i X^i X i と X j X^j X j の共変分 、⟨ X i ⟩ : = ⟨ X i , X i ⟩ \langle X^i\rangle:=\langle X^i,X^i\rangle ⟨ X i ⟩ := ⟨ X i , X i ⟩ を X i X^i X i の 2 次変分 と呼びます。ドリフト μ \mu μ は共変分に寄与しません。
命題 3.2 により、X X X のマルチンゲール部分 M i = ∫ σ i ⋅ d W M^i=\int\sigma^i\cdot dW M i = ∫ σ i ⋅ d W について M i M j − ⟨ X i , X j ⟩ M^iM^j-\langle X^i,X^j\rangle M i M j − ⟨ X i , X j ⟩ はマルチンゲールです。たとえば ⟨ W a , W b ⟩ t = δ a b t \langle W^a,W^b\rangle_t=\delta_{ab}t ⟨ W a , W b ⟩ t = δ ab t です。
2 次変分の収束を示すには、増分の 2 乗の揺らぎを抑える必要があり、そのために 4 次モーメントの評価を使います。
次の補題の証明は二段構えです。まず単純過程について、各区間の寄与 Y l Y_l Y l が、それまでの情報を固定すると正規分布に従うことを使い、部分和 S l S_l S l の 4 乗の期待値を区間ごとに足し上げます。そのうえで、一般の場合を単純過程の極限として得ます。
補題 3.5 (伊藤積分の 4 次モーメント ) ✓
H H H を R d \mathbb{R}^d R d 値の発展的可測過程で ∣ H u ∣ ≤ L |H_u|\le L ∣ H u ∣ ≤ L とし、0 ≤ s < t ≤ T 0\le s<t\le T 0 ≤ s < t ≤ T とします。このとき
E [ ( ∫ s t H u ⋅ d W u ) 4 ] ≤ 6 L 4 ( t − s ) 2 . \mathbb{E}\left[\left(\int_s^tH_u\cdot dW_u\right)^4\right]\le6L^4(t-s)^2. E [ ( ∫ s t H u ⋅ d W u ) 4 ] ≤ 6 L 4 ( t − s ) 2 .
証明(補題 3.5) 単純過程の場合。 H = ∑ l ξ l 1 ( u l , u l + 1 ] H=\sum_{l}\xi_l\mathbf{1}_{(u_l,u_{l+1}]} H = ∑ l ξ l 1 ( u l , u l + 1 ] (s = u 0 < ⋯ < u m = t s=u_0<\dots<u_m=t s = u 0 < ⋯ < u m = t 、ξ l \xi_l ξ l は F u l \mathcal{F}_{u_l} F u l -可測で ∣ ξ l ∣ ≤ L |\xi_l|\le L ∣ ξ l ∣ ≤ L )とし、Y l = ξ l ⋅ ( W u l + 1 − W u l ) Y_l=\xi_l\cdot(W_{u_{l+1}}-W_{u_l}) Y l = ξ l ⋅ ( W u l + 1 − W u l ) 、S l = ∑ k < l Y k S_l=\sum_{k<l}Y_k S l = ∑ k < l Y k 、h l = u l + 1 − u l h_l=u_{l+1}-u_l h l = u l + 1 − u l とおきます。W u l + 1 − W u l W_{u_{l+1}}-W_{u_l} W u l + 1 − W u l は F u l \mathcal{F}_{u_l} F u l と独立に N ( 0 , h l I d ) N(0,h_lI_d) N ( 0 , h l I d ) に従うので、補題 2.1[Feller 過程と Feynman–Kac] (関数の正部分と負部分に分けて適用)により、F u l \mathcal{F}_{u_l} F u l に関する条件付きで Y l Y_l Y l は N ( 0 , ∣ ξ l ∣ 2 h l ) N(0,|\xi_l|^2h_l) N ( 0 , ∣ ξ l ∣ 2 h l ) に従います。したがって
E [ Y l ∣ F u l ] = E [ Y l 3 ∣ F u l ] = 0 , E [ Y l 2 ∣ F u l ] = ∣ ξ l ∣ 2 h l , E [ Y l 4 ∣ F u l ] = 3 ∣ ξ l ∣ 4 h l 2 . \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 . E [ Y l ∣ F u l ] = E [ Y l 3 ∣ F u l ] = 0 , E [ Y l 2 ∣ F u l ] = ∣ ξ l ∣ 2 h l , E [ Y l 4 ∣ F u l ] = 3∣ ξ l ∣ 4 h l 2 . ∣ S l ∣ ≤ L ∑ k < l ∣ W u k + 1 − W u k ∣ |S_l|\le L\sum_{k<l}|W_{u_{k+1}}-W_{u_k}| ∣ S l ∣ ≤ L ∑ k < l ∣ W u k + 1 − W u k ∣ はすべての次数のモーメントを持つので、( S l + Y l ) 4 (S_l+Y_l)^4 ( S l + Y l ) 4 を展開して 定理 4.4[条件付き期待値] を使うと、奇数次の項が消えて
E [ S l + 1 4 ] = E [ S l 4 ] + 6 E [ S l 2 ∣ ξ l ∣ 2 ] h l + 3 E [ ∣ ξ l ∣ 4 ] h l 2 ≤ E [ S l 4 ] + 6 L 2 E [ S l 2 ] h l + 3 L 4 h l 2 . \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 [ S l + 1 4 ] = E [ S l 4 ] + 6 E [ S l 2 ∣ ξ l ∣ 2 ] h l + 3 E [ ∣ ξ l ∣ 4 ] h l 2 ≤ E [ S l 4 ] + 6 L 2 E [ S l 2 ] h l + 3 L 4 h l 2 . 同様に E [ S l 2 ] = ∑ k < l E [ ∣ ξ k ∣ 2 ] h k ≤ L 2 ( u l − s ) \mathbb{E}[S_l^2]=\sum_{k<l}\mathbb{E}[|\xi_k|^2]h_k\le L^2(u_l-s) E [ S l 2 ] = ∑ k < l E [ ∣ ξ k ∣ 2 ] h k ≤ L 2 ( u l − s ) です。これを l = 0 , … , m − 1 l=0,\dots,m-1 l = 0 , … , m − 1 について足すと、∑ l ( u l − s ) h l ≤ ∫ s t ( u − s ) d u = ( t − s ) 2 / 2 \sum_l(u_l-s)h_l\le\int_s^t(u-s)\,du=(t-s)^2/2 ∑ l ( u l − s ) h l ≤ ∫ s t ( u − s ) d u = ( t − s ) 2 /2 と ∑ l h l 2 ≤ ( t − s ) 2 \sum_lh_l^2\le(t-s)^2 ∑ l h l 2 ≤ ( t − s ) 2 から
E [ S m 4 ] ≤ 6 L 4 ⋅ ( t − s ) 2 2 + 3 L 4 ( t − s ) 2 = 6 L 4 ( 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 . E [ S m 4 ] ≤ 6 L 4 ⋅ 2 ( t − s ) 2 + 3 L 4 ( t − s ) 2 = 6 L 4 ( t − s ) 2 . 一般の場合。 補題 2.3 により ∣ H ( n ) ∣ ≤ L |H^{(n)}|\le L ∣ H ( n ) ∣ ≤ L なる単純過程で E ∫ ∣ H ( n ) − H ∣ 2 → 0 \mathbb{E}\int|H^{(n)}-H|^2\to0 E ∫ ∣ H ( n ) − H ∣ 2 → 0 となるものを取ります(分割に s , t s,t s , t を加えます)。補題 3.1 から ∫ s t H ( n ) ⋅ d W → ∫ s t H ⋅ d W \int_s^tH^{(n)}\cdot dW\to\int_s^tH\cdot dW ∫ s t H ( n ) ⋅ d W → ∫ s t H ⋅ d W が L 2 ( Ω ) L^2(\Omega) L 2 ( Ω ) で成り立つので、ほとんど確実に収束する部分列が取れます。実際、Y n = ∫ s t H ( n ) ⋅ d W Y_n=\int_s^tH^{(n)}\cdot dW Y n = ∫ s t H ( n ) ⋅ d W 、Y = ∫ s t H ⋅ d W Y=\int_s^tH\cdot dW Y = ∫ s t H ⋅ d W として E ∣ Y n k − Y ∣ 2 ≤ 8 − k \mathbb{E}|Y_{n_k}-Y|^2\le 8^{-k} E ∣ Y n k − Y ∣ 2 ≤ 8 − k となる n 1 < n 2 < ⋯ n_1<n_2<\cdots n 1 < n 2 < ⋯ を選ぶと、マルコフの不等式(命題 3.5[大数の法則と中心極限定理] )を ∣ Y n k − Y ∣ 2 |Y_{n_k}-Y|^2 ∣ Y n k − Y ∣ 2 と a = 4 − k a=4^{-k} a = 4 − k に使って P ( ∣ Y n k − Y ∣ ≥ 2 − k ) ≤ 4 k ⋅ 8 − k = 2 − k \mathbb{P}(|Y_{n_k}-Y|\ge2^{-k})\le4^k\cdot8^{-k}=2^{-k} P ( ∣ Y n k − Y ∣ ≥ 2 − k ) ≤ 4 k ⋅ 8 − k = 2 − k で、これは k k k について和が有限です。ボレル・カンテリの補題(補題 4.4[大数の法則と中心極限定理] )から、ほとんど確実に有限個の k k k を除いて ∣ Y n k − Y ∣ < 2 − k |Y_{n_k}-Y|<2^{-k} ∣ Y n k − Y ∣ < 2 − k となり、Y n k → Y Y_{n_k}\to Y Y n k → Y a.s. です。補題 6.1[ルベーグ積分の定義と収束定理] により、極限でも同じ上界が成り立ちます。
∎
上界が ( t − s ) 2 (t-s)^2 ( t − s ) 2 の定数倍であることが重要です。積分の 2 乗は t − s t-s t − s の大きさなので、4 次モーメントはその 2 乗の大きさに収まり、増分の 2 乗の揺らぎは平均に比べて小さいことになります。
定義 3.3 の共変分が、増分の積の和の極限として実際に現れることを示します。これが 1 次元の 定理 3.1[確率微分方程式] の多次元・伊藤過程版で、伊藤公式の 2 次の項の出どころです。
次の定理の証明では、増分 Δ X k i Δ X k j \Delta X_k^i\Delta X_k^j Δ X k i Δ X k j の和を、増分をドリフト部分とマルチンゲール部分に分けて展開し、ドリフトを含む項、マルチンゲール部分の項、Riemann 和の項に整理して、それぞれ 0 0 0 または目的の積分に収束することを見ます。ドリフトを含む項は、ドリフトの寄与が時間幅に比例して小さいため、分割を細かくすれば 0 0 0 に近づきます。マルチンゲール部分の項は、各区間の誤差 Z k Z_k Z k が互いに直交するので 2 乗の和が小さくなります。残る Riemann 和の項は、連続な被積分関数の積分への通常の収束です。
定理 4.1 (共変分の離散近似 ) ✓
X X X を有界係数の条件 (A)(定数 K K K )を満たす n n n 次元伊藤過程、G G G を適合した実数値過程で、道が連続かつ ∣ G ∣ ≤ C |G|\le C ∣ G ∣ ≤ C なるものとします。t ∈ ( 0 , T ] t\in(0,T] t ∈ ( 0 , T ] 、Π = { 0 = t 0 < ⋯ < t m = t } \Pi=\{0=t_0<\dots<t_m=t\} Π = { 0 = t 0 < ⋯ < t m = t } を分割、Δ X k i = X t k + 1 i − X t k i \Delta X^i_k=X^i_{t_{k+1}}-X^i_{t_k} Δ X k i = X t k + 1 i − X t k i とおくと
E ∣ ∑ k G t k Δ X k i Δ X k j − ∫ 0 t G u d ⟨ X i , X j ⟩ 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). E k ∑ G t k Δ X k i Δ X k j − ∫ 0 t G u d ⟨ X i , X j ⟩ u ⟶ 0 ( ∥Π∥ → 0 ) . ここで ∫ 0 t G u d ⟨ X i , X j ⟩ u = ∫ 0 t G u σ u i ⋅ σ u j d u \int_0^tG_u\,d\langle X^i,X^j\rangle_u=\int_0^tG_u\,\sigma^i_u\cdot\sigma^j_u\,du ∫ 0 t G u d ⟨ X i , X j ⟩ u = ∫ 0 t G u σ u i ⋅ σ u j d u です。
証明(定理 4.1) h k = t k + 1 − t k h_k=t_{k+1}-t_k h k = t k + 1 − t k とし、増分をドリフト部分とマルチンゲール部分に分けます。
Δ X k i = A k i + N k i , A k i = ∫ t k t k + 1 μ u i d u , N k i = ∫ t k t k + 1 σ u i ⋅ d W u . \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 . Δ X k i = A k i + N k i , A k i = ∫ t k t k + 1 μ u i d u , N k i = ∫ t k t k + 1 σ u i ⋅ d W u . ∣ A k i ∣ ≤ K h k |A^i_k|\le Kh_k ∣ A k i ∣ ≤ K h k で、補題 3.1 と ∣ σ i ∣ 2 ≤ d K 2 |\sigma^i|^2\le dK^2 ∣ σ i ∣ 2 ≤ d K 2 から E [ ( N k i ) 2 ] ≤ d K 2 h k \mathbb{E}[(N^i_k)^2]\le dK^2h_k E [( N k i ) 2 ] ≤ d K 2 h k です。
ドリフトを含む項。 ∑ k ∣ G t k A k i A k j ∣ ≤ C K 2 ∥ Π ∥ t \sum_k|G_{t_k}A^i_kA^j_k|\le CK^2\|\Pi\|t ∑ k ∣ G t k A k i A k j ∣ ≤ C K 2 ∥Π∥ t です。また Cauchy–Schwarz の不等式から
E ∑ k ∣ G t k A k i N k j ∣ ≤ C K ∑ k h k ( E [ ( N k j ) 2 ] ) 1 / 2 ≤ C d K 2 ∑ k h k 3 / 2 ≤ C d K 2 ∥ Π ∥ 1 / 2 t , \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, E k ∑ ∣ G t k A k i N k j ∣ ≤ C K k ∑ h k ( E [( N k j ) 2 ] ) 1/2 ≤ C d K 2 k ∑ h k 3/2 ≤ C d K 2 ∥Π ∥ 1/2 t , A j N i A^jN^i A j N i の項も同じです。いずれも 0 0 0 に収束します。
マルチンゲール部分の項。 C k = ∫ t k t k + 1 σ i ⋅ σ j d u C_k=\int_{t_k}^{t_{k+1}}\sigma^i\cdot\sigma^j\,du C k = ∫ t k t k + 1 σ i ⋅ σ j d u 、Z k = G t k ( N k i N k j − C k ) Z_k=G_{t_k}(N^i_kN^j_k-C_k) Z k = G t k ( N k i N k j − C k ) とおきます。∣ C k ∣ ≤ d K 2 h k |C_k|\le dK^2h_k ∣ C k ∣ ≤ d K 2 h k で、補題 3.5 (L = d K L=\sqrt dK L = d K )と Cauchy–Schwarz の不等式から
E [ Z k 2 ] ≤ 2 C 2 ( E [ ( N k i ) 4 ] 1 / 2 E [ ( N k j ) 4 ] 1 / 2 + d 2 K 4 h k 2 ) ≤ 14 C 2 d 2 K 4 h k 2 . \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 . E [ Z k 2 ] ≤ 2 C 2 ( E [( N k i ) 4 ] 1/2 E [( N k j ) 4 ] 1/2 + d 2 K 4 h k 2 ) ≤ 14 C 2 d 2 K 4 h k 2 . 次に、ξ \xi ξ を有界な F t k \mathcal{F}_{t_k} F t k -可測確率変数とすると、N k i N k j = ∑ a , b ∫ σ i a d W a ∫ σ j b d W b N^i_kN^j_k=\sum_{a,b}\int\sigma^{ia}\,dW^a\int\sigma^{jb}\,dW^b N k i N k j = ∑ a , b ∫ σ ia d W a ∫ σ j b d W b に 補題 3.1 を ξ G t k \xi G_{t_k} ξ G t k について適用して E [ ξ Z k ] = 0 \mathbb{E}[\xi Z_k]=0 E [ ξ Z k ] = 0 を得ます。ξ \xi ξ として F t k \mathcal{F}_{t_k} F t k の事象の指示関数を取れば E [ Z k ∣ F t k ] = 0 \mathbb{E}[Z_k\mid\mathcal{F}_{t_k}]=0 E [ Z k ∣ F t k ] = 0 です。k < l k<l k < l なら Z k Z_k Z k は F t l \mathcal{F}_{t_l} F t l -可測で二乗可積分なので、定理 4.4[条件付き期待値] により E [ Z k Z l ] = E [ Z k E [ Z l ∣ F t l ] ] = 0 \mathbb{E}[Z_kZ_l]=\mathbb{E}[Z_k\,\mathbb{E}[Z_l\mid\mathcal{F}_{t_l}]]=0 E [ Z k Z l ] = E [ Z k E [ Z l ∣ F t l ]] = 0 です。よって交差項が消え、
E [ ( ∑ k Z k ) 2 ] = ∑ k E [ Z k 2 ] ≤ 14 C 2 d 2 K 4 ∥ Π ∥ 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 . E [ ( k ∑ Z k ) 2 ] = k ∑ E [ Z k 2 ] ≤ 14 C 2 d 2 K 4 ∥Π∥ t ⟶ 0. Riemann 和の項。 残るのは ∑ k G t k C k = ∫ 0 t G u Π σ u i ⋅ σ u j d u \sum_kG_{t_k}C_k=\int_0^tG^\Pi_u\,\sigma^i_u\cdot\sigma^j_u\,du ∑ k G t k C k = ∫ 0 t G u Π σ u i ⋅ σ u j d u (G u Π = G t k G^\Pi_u=G_{t_k} G u Π = G t k , u ∈ ( t k , t k + 1 ] u\in(t_k,t_{k+1}] u ∈ ( t k , t k + 1 ] )と ∫ 0 t G u σ i ⋅ σ j d u \int_0^tG_u\,\sigma^i\cdot\sigma^j\,du ∫ 0 t G u σ i ⋅ σ j d u の差です。道の連続性から G u Π → G u G^\Pi_u\to G_u G u Π → G u が各 ( u , ω ) (u,\omega) ( u , ω ) で成り立ち、被積分関数は 2 C d K 2 2CdK^2 2 C d K 2 で抑えられるので、定理 7.3[ルベーグ積分の定義と収束定理] (d u ⊗ d P du\otimes d\mathbb{P} d u ⊗ d P について)から差の期待値は 0 0 0 に収束します。
以上の 4 種類の項を合わせると主張が出ます(L 2 L^2 L 2 収束は L 1 L^1 L 1 収束を導きます)。
∎
G = 1 G=1 G = 1 、X = W X=W X = W とすれば ∑ k Δ W k a Δ W k b → δ a b t \sum_k\Delta W^a_k\Delta W^b_k\to\delta_{ab}t ∑ k Δ W k a Δ W k b → δ ab t で、異なる成分の増分の積の和が 0 0 0 に収束することが分かります。符号がランダムな小さい量 Δ W k a Δ W k b \Delta W^a_k\Delta W^b_k Δ W k a Δ W k b (大きさ h k h_k h k )が 1 / h k 1/h_k 1/ h k 個足されても、互いに打ち消し合って和は ∥ Π ∥ \sqrt{\|\Pi\|} ∥Π∥ 程度にしかならない、というのが証明の中身です。
通常の微積分の連鎖律では f f f の 1 次の項だけが残りますが、確率過程では増分の 2 乗が時間幅と同じ大きさ(h \sqrt{h} h の 2 乗)になるため、2 次の項も残ります。次の定理の証明は、時間を細かく刻んで f f f の変化を足し合わせ(望遠鏡和)、各区間で 2 次まで Taylor 展開します。1 次の項は確率積分とドリフトの積分に、2 次の項は前節の定理により共変分の積分に収束します。3 次以上の項(剰余)は、増分の 3 乗が h 3 / 2 h^{3/2} h 3/2 の大きさなので、足し合わせても 0 0 0 に近づきます。
定理 5.1 (多次元伊藤公式 ) ✓
X X X を有界係数の条件 (A) を満たす n n n 次元伊藤過程とし、f ∈ C 1 , 2 ( [ 0 , T ] × R n ) f\in C^{1,2}([0,T]\times\mathbb{R}^n) f ∈ C 1 , 2 ([ 0 , T ] × R n ) (t t t について 1 回、x x x について 2 回連続微分可能)で、x x x についての 1 階・2 階の偏導関数 ∂ i f \partial_if ∂ i f 、∂ i j f \partial_{ij}f ∂ ij f が有界であるとします。このとき、ほとんど確実にすべての t ∈ [ 0 , T ] t\in[0,T] t ∈ [ 0 , T ] について
f ( t , X t ) = f ( 0 , X 0 ) + ∫ 0 t ∂ t f ( u , X u ) d u + ∑ i ∫ 0 t ∂ i f ( u , X u ) μ u i d u + ∑ i , a ∫ 0 t ∂ i f ( u , X u ) σ u i a d W u a + 1 2 ∑ i , j ∫ 0 t ∂ i j f ( u , X u ) d ⟨ X i , X j ⟩ 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} f ( t , X t ) = f ( 0 , X 0 ) + ∫ 0 t ∂ t f ( u , X u ) d u + i ∑ ∫ 0 t ∂ i f ( u , X u ) μ u i d u + i , a ∑ ∫ 0 t ∂ i f ( u , X u ) σ u ia d W u a + 2 1 i , j ∑ ∫ 0 t ∂ ij f ( u , X u ) d ⟨ X i , X j ⟩ u が成り立ちます。微分形式では d f = ∂ t f d t + ∑ i ∂ i f d X i + 1 2 ∑ i , j ∂ i j f d ⟨ X i , X j ⟩ df=\partial_tf\,dt+\sum_i\partial_if\,dX^i+\frac12\sum_{i,j}\partial_{ij}f\,d\langle X^i,X^j\rangle df = ∂ t f d t + ∑ i ∂ i f d X i + 2 1 ∑ i , j ∂ ij f d ⟨ X i , X j ⟩ です。
証明(定理 5.1) 被積分関数 ∂ i f ( u , X u ) σ u i a \partial_if(u,X_u)\sigma^{ia}_u ∂ i f ( u , X u ) σ u ia は有界で、連続な適合過程と発展的可測過程の積なので発展的可測です。よって右辺の確率積分は定義されます。両辺は t t t について連続なので、t t t を固定して等式がほとんど確実に成り立つことを示せば、有理数の t t t で同時に成り立ち、連続性からすべての t t t に及びます。
t t t を固定し、分割 Π = { 0 = t 0 < ⋯ < t m = t } \Pi=\{0=t_0<\dots<t_m=t\} Π = { 0 = t 0 < ⋯ < t m = t } を取ります。X k = X t k X_k=X_{t_k} X k = X t k 、Δ X k = X k + 1 − X k \Delta X_k=X_{k+1}-X_k Δ X k = X k + 1 − X k 、h k = t k + 1 − t k h_k=t_{k+1}-t_k h k = t k + 1 − t k と書き、望遠鏡和を時間方向と空間方向に分けます。
f ( t , X t ) − f ( 0 , X 0 ) = ∑ k [ f ( t k + 1 , X k + 1 ) − f ( t k , X k + 1 ) ] + ∑ k [ f ( t k , X k + 1 ) − f ( t k , X k ) ] . 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]. f ( t , X t ) − f ( 0 , X 0 ) = k ∑ [ f ( t k + 1 , X k + 1 ) − f ( t k , X k + 1 ) ] + k ∑ [ f ( t k , X k + 1 ) − f ( t k , X k ) ] . 以下、道 ω \omega ω を固定したときの量 R ( ω ) = sup u ≤ t ∣ X u ( ω ) ∣ R(\omega)=\sup_{u\le t}|X_u(\omega)| R ( ω ) = sup u ≤ t ∣ X u ( ω ) ∣ は有限で、[ 0 , T ] × { ∣ x ∣ ≤ R + 1 } [0,T]\times\{|x|\le R+1\} [ 0 , T ] × { ∣ x ∣ ≤ R + 1 } はコンパクトなので ∂ t f \partial_tf ∂ t f と ∂ i j f \partial_{ij}f ∂ ij f はそこで一様連続です。
時間方向。 第 1 の和の各項は ∫ t k t k + 1 ∂ t f ( u , X k + 1 ) d u \int_{t_k}^{t_{k+1}}\partial_tf(u,X_{k+1})\,du ∫ t k t k + 1 ∂ t f ( u , X k + 1 ) d u です。これと ∫ t k t k + 1 ∂ t f ( u , X u ) d u \int_{t_k}^{t_{k+1}}\partial_tf(u,X_u)\,du ∫ t k t k + 1 ∂ t f ( u , X u ) d u の差の和は、t ⋅ max k sup u ∈ [ t k , t k + 1 ] ∣ ∂ t f ( u , X k + 1 ) − ∂ t f ( u , X u ) ∣ t\cdot\max_k\sup_{u\in[t_k,t_{k+1}]}|\partial_tf(u,X_{k+1})-\partial_tf(u,X_u)| t ⋅ max k sup u ∈ [ t k , t k + 1 ] ∣ ∂ t f ( u , X k + 1 ) − ∂ t f ( u , X u ) ∣ 以下で、∂ t f \partial_tf ∂ t f の一様連続性と道の一様連続性から各 ω \omega ω で 0 0 0 に収束します。
空間方向の Taylor 展開。 積分形の剰余をもつ Taylor の定理から
f ( t k , X k + 1 ) − f ( t k , X k ) = ∑ i ∂ i f ( t k , X k ) Δ X k i + 1 2 ∑ i , j ∂ i j f ( t k , X k ) Δ X k i Δ X k j + ρ 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, f ( t k , X k + 1 ) − f ( t k , X k ) = i ∑ ∂ i f ( t k , X k ) Δ X k i + 2 1 i , j ∑ ∂ ij f ( t k , X k ) Δ X k i Δ X k j + ρ k , ∣ ρ k ∣ ≤ 1 2 ε Π ∣ Δ X k ∣ 2 , ε Π : = max k sup 0 ≤ θ ≤ 1 ∥ ∇ 2 f ( t k , X k + θ Δ X k ) − ∇ 2 f ( t k , X k ) ∥ |\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\| ∣ ρ k ∣ ≤ 2 1 ε Π ∣Δ X k ∣ 2 , ε Π := k max 0 ≤ θ ≤ 1 sup ∇ 2 f ( t k , X k + θ Δ X k ) − ∇ 2 f ( t k , X k ) です(∥ ⋅ ∥ \|\cdot\| ∥ ⋅ ∥ は行列の作用素ノルム)。
1 次の項。 左端の値で作った過程 Y u Π , i : = ∂ i f ( t k , X k ) Y^{\Pi,i}_u:=\partial_if(t_k,X_k) Y u Π , i := ∂ i f ( t k , X k ) (u ∈ ( t k , t k + 1 ] u\in(t_k,t_{k+1}] u ∈ ( t k , t k + 1 ] )を使うと、補題 2.4 (ξ = ∂ i f ( t k , X k ) \xi=\partial_if(t_k,X_k) ξ = ∂ i f ( t k , X k ) )により
∑ k ∂ i f ( t k , X k ) Δ X k i = ∫ 0 t Y u Π , i μ u i d u + ∑ a ∫ 0 t Y u Π , i σ u i a d W u a . \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 . k ∑ ∂ i f ( t k , X k ) Δ X k i = ∫ 0 t Y u Π , i μ u i d u + a ∑ ∫ 0 t Y u Π , i σ u ia d W u a . 道の連続性から Y u Π , i → ∂ i f ( u , X u ) Y^{\Pi,i}_u\to\partial_if(u,X_u) Y u Π , i → ∂ i f ( u , X u ) が各 ( u , ω ) (u,\omega) ( u , ω ) で成り立ち、∣ Y Π , i ∣ ≤ sup ∣ ∂ i f ∣ |Y^{\Pi,i}|\le\sup|\partial_if| ∣ Y Π , i ∣ ≤ sup ∣ ∂ i f ∣ です。ドリフトの積分は優収束定理により各 ω \omega ω で収束します。確率積分は、E ∫ 0 t ∣ Y u Π , i − ∂ i f ( u , X u ) ∣ 2 ( σ u i a ) 2 d u → 0 \mathbb{E}\int_0^t|Y^{\Pi,i}_u-\partial_if(u,X_u)|^2(\sigma^{ia}_u)^2\,du\to0 E ∫ 0 t ∣ Y u Π , i − ∂ i f ( u , X u ) ∣ 2 ( σ u ia ) 2 d u → 0 (優収束定理)と伊藤等長性から L 2 ( Ω ) L^2(\Omega) L 2 ( Ω ) で収束します。
2 次の項。 G u = ∂ i j f ( u , X u ) G_u=\partial_{ij}f(u,X_u) G u = ∂ ij f ( u , X u ) は有界で道が連続な適合過程なので、定理 4.1 により ∑ k ∂ i j f ( t k , X k ) Δ X k i Δ X k j → ∫ 0 t ∂ i j f ( u , X u ) d ⟨ X i , X j ⟩ 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 ∑ k ∂ ij f ( t k , X k ) Δ X k i Δ X k j → ∫ 0 t ∂ ij f ( u , X u ) d ⟨ X i , X j ⟩ u が L 1 L^1 L 1 で成り立ちます。
剰余項。 S Π : = ∑ k ∣ Δ X k ∣ 2 = ∑ i ∑ k ( Δ X k i ) 2 S_\Pi:=\sum_k|\Delta X_k|^2=\sum_i\sum_k(\Delta X^i_k)^2 S Π := ∑ k ∣Δ X k ∣ 2 = ∑ i ∑ k ( Δ X k i ) 2 は 定理 4.1 (G = 1 G=1 G = 1 )により ∑ i ⟨ X i ⟩ t \sum_i\langle X^i\rangle_t ∑ i ⟨ X i ⟩ t に L 1 L^1 L 1 で収束するので、sup Π E [ S Π ] < ∞ \sup_\Pi\mathbb{E}[S_\Pi]<\infty sup Π E [ S Π ] < ∞ です。一方 max k ∣ Δ X k ∣ → 0 \max_k|\Delta X_k|\to0 max k ∣Δ X k ∣ → 0 (道の一様連続性)と ∇ 2 f \nabla^2f ∇ 2 f の一様連続性から ε Π → 0 \varepsilon_\Pi\to0 ε Π → 0 が各 ω \omega ω で成り立ちます。任意の δ , η > 0 \delta,\eta>0 δ , η > 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} P ( ∑ k ρ k > η ) ≤ P ( ε Π > δ ) + P ( S Π > 2 η / δ ) ≤ P ( ε Π > δ ) + 2 η δ E [ S Π ] で、∥ Π ∥ → 0 \|\Pi\|\to0 ∥Π∥ → 0 の後に δ → 0 \delta\to0 δ → 0 とすれば右辺は 0 0 0 に収束します。
以上より右辺の各項は ∥ Π ∥ → 0 \|\Pi\|\to0 ∥Π∥ → 0 で確率収束し、その極限の和が公式の右辺です。左辺は Π \Pi Π によらないので、確率収束の極限の一意性から等式がほとんど確実に成り立ちます。
∎
1 次元の公式(定理 5.2[確率微分方程式] )との違いは、2 次の項に ∂ i j f \partial_{ij}f ∂ ij f と共変分の組が全部現れることだけです。覚え方は d W a d W b = δ a b d t dW^a\,dW^b=\delta_{ab}\,dt d W a d W b = δ ab d t 、d t d W a = 0 dt\,dW^a=0 d t d W a = 0 、( d t ) 2 = 0 (dt)^2=0 ( d t ) 2 = 0 のもとで 2 次まで Taylor 展開する、です。
証明(系 5.3) v < w v<w v < w について、補題 3.1 から
E ∣ X w − X v ∣ 2 ≤ 2 E ( ∫ v w μ X d u ) 2 + 2 E ( ∫ v w σ X ⋅ d W ) 2 ≤ 2 K 2 ( w − v ) 2 + 2 d K 2 ( 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 ∣ X w − X v ∣ 2 ≤ 2 E ( ∫ v w μ X d u ) 2 + 2 E ( ∫ v w σ X ⋅ d W ) 2 ≤ 2 K 2 ( w − v ) 2 + 2 d K 2 ( w − v ) . 特に E [ X u 2 ] ≤ 2 E [ X 0 2 ] + 4 K 2 T 2 + 4 d K 2 T \mathbb{E}[X_u^2]\le2\,\mathbb{E}[X_0^2]+4K^2T^2+4dK^2T E [ X u 2 ] ≤ 2 E [ X 0 2 ] + 4 K 2 T 2 + 4 d K 2 T は u u u について有界で、∣ σ Y ∣ ≤ d K |\sigma^Y|\le\sqrt dK ∣ σ Y ∣ ≤ d K なので X σ Y ∈ L d 2 X\sigma^Y\in\mathcal{L}^2_d X σ Y ∈ L d 2 です。
分割 Π \Pi Π について、代数的な恒等式
X t Y t − X 0 Y 0 = ∑ k X t k Δ Y k + ∑ k Y t k Δ X k + ∑ k Δ X k Δ Y k X_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 X t Y t − X 0 Y 0 = k ∑ X t k Δ Y k + k ∑ Y t k Δ X k + k ∑ Δ X k Δ Y k が剰余なしに成り立ちます。第 3 の和は 定理 4.1 (G = 1 G=1 G = 1 )により ⟨ X , Y ⟩ t \langle X,Y\rangle_t ⟨ X , Y ⟩ t に L 1 L^1 L 1 で収束します。第 1 の和は、X u Π : = X t k X^\Pi_u:=X_{t_k} X u Π := X t k (u ∈ ( t k , t k + 1 ] u\in(t_k,t_{k+1}] u ∈ ( t k , t k + 1 ] )とおくと 補題 2.4 の近似版により ∫ 0 t X Π μ Y d u + ∫ 0 t X Π σ Y ⋅ d W \int_0^tX^\Pi\mu^Y\,du+\int_0^tX^\Pi\sigma^Y\cdot dW ∫ 0 t X Π μ Y d u + ∫ 0 t X Π σ Y ⋅ d W に等しいです(X t k X_{t_k} X t k は有界とは限らないので、X t k X_{t_k} X t k を ( − N ) ∨ X t k ∧ N (-N)\vee X_{t_k}\wedge N ( − N ) ∨ X t k ∧ N で置き換えた等式で N → ∞ N\to\infty N → ∞ とし、両辺の L 2 L^2 L 2 収束を使います)。ドリフトの積分は道ごとに優収束定理で ∫ 0 t X μ Y d u \int_0^tX\mu^Y\,du ∫ 0 t X μ Y d u に収束します(∣ X Π ∣ ≤ sup u ≤ t ∣ X u ∣ < ∞ |X^\Pi|\le\sup_{u\le t}|X_u|<\infty ∣ X Π ∣ ≤ sup u ≤ t ∣ X u ∣ < ∞ )。確率積分については、上の評価から
E ∫ 0 t ∣ X u Π − X u ∣ 2 ∣ σ u Y ∣ 2 d u ≤ d K 2 t ( 2 K 2 ∥ Π ∥ 2 + 2 d K 2 ∥ Π ∥ ) ⟶ 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 E ∫ 0 t ∣ X u Π − X u ∣ 2 ∣ σ u Y ∣ 2 d u ≤ d K 2 t ( 2 K 2 ∥Π ∥ 2 + 2 d K 2 ∥Π∥ ) ⟶ 0 なので、等長性から L 2 ( Ω ) L^2(\Omega) L 2 ( Ω ) で ∫ 0 t X σ Y ⋅ d W \int_0^tX\sigma^Y\cdot dW ∫ 0 t X σ Y ⋅ d W に収束します。第 2 の和も同じです。左辺は Π \Pi Π によらないので等式が従い、両辺の連続性からすべての t t t で同時に成り立ちます。
∎
例 5.4 (Brown 運動のノルムの 2 乗 )
X = Y = W a X=Y=W^a X = Y = W a に 系 5.3 を使うと ( W t a ) 2 = 2 ∫ 0 t W a d W a + t (W^a_t)^2=2\int_0^tW^a\,dW^a+t ( W t a ) 2 = 2 ∫ 0 t W a d W a + t です。a a a について足すと
∣ W t ∣ 2 = 2 ∫ 0 t W u ⋅ d W u + d t . |W_t|^2=2\int_0^tW_u\cdot dW_u+d\,t . ∣ W t ∣ 2 = 2 ∫ 0 t W u ⋅ d W u + d t . 右辺の確率積分は平均 0 0 0 なので E ∣ W t ∣ 2 = d t \mathbb{E}|W_t|^2=dt E ∣ W t ∣ 2 = d t です。また a ≠ b a\ne b a = b なら ⟨ W a , W b ⟩ = 0 \langle W^a,W^b\rangle=0 ⟨ W a , W b ⟩ = 0 なので W t a W t b = ∫ 0 t W a d W b + ∫ 0 t W b d W a W^a_tW^b_t=\int_0^tW^a\,dW^b+\int_0^tW^b\,dW^a W t a W t b = ∫ 0 t W a d W b + ∫ 0 t W b d W a で、W a W b W^aW^b W a W b はマルチンゲールです。
例 5.5 (後ろ向き熱方程式の解とマルチンゲール )
u ∈ R d u\in\mathbb{R}^d u ∈ R d を固定し、f ( t , x ) = cos ( u ⋅ x ) e ∣ u ∣ 2 t / 2 f(t,x)=\cos(u\cdot x)\,e^{|u|^2t/2} f ( t , x ) = cos ( u ⋅ x ) e ∣ u ∣ 2 t /2 とおきます。∂ t f = ∣ u ∣ 2 2 f \partial_tf=\frac{|u|^2}{2}f ∂ t f = 2 ∣ u ∣ 2 f 、∑ a ∂ a a f = − ∣ u ∣ 2 f \sum_a\partial_{aa}f=-|u|^2f ∑ a ∂ aa f = − ∣ u ∣ 2 f なので ∂ t f + 1 2 Δ f = 0 \partial_tf+\frac12\Delta f=0 ∂ t f + 2 1 Δ f = 0 です。[ 0 , T ] [0,T] [ 0 , T ] 上で f f f の x x x についての 1 階・2 階導関数は有界なので、定理 5.1 (X = W X=W X = W 、μ = 0 \mu=0 μ = 0 、σ = I d \sigma=I_d σ = I d )から
cos ( u ⋅ W t ) e ∣ u ∣ 2 t / 2 = 1 − ∫ 0 t sin ( u ⋅ W s ) e ∣ u ∣ 2 s / 2 u ⋅ d W s . \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 . cos ( u ⋅ W t ) e ∣ u ∣ 2 t /2 = 1 − ∫ 0 t sin ( u ⋅ W s ) e ∣ u ∣ 2 s /2 u ⋅ d W s . 右辺はマルチンゲールなので E [ cos ( u ⋅ W t ) ] = e − ∣ u ∣ 2 t / 2 \mathbb{E}[\cos(u\cdot W_t)]=e^{-|u|^2t/2} E [ cos ( u ⋅ W t )] = e − ∣ u ∣ 2 t /2 で、Gauss 分布の特性関数が伊藤公式から再現されます。一般に ∂ t f + 1 2 Δ f = 0 \partial_tf+\frac12\Delta f=0 ∂ t f + 2 1 Δ f = 0 を満たす f f f について f ( t , W t ) f(t,W_t) f ( t , W t ) はマルチンゲールになり、これが Feynman–Kac 公式(Feller 過程の章 )の最も簡単な場合です。
マルチンゲールとブラウン運動の章 の 定理 6.1[マルチンゲールとブラウン運動] で、Brown 運動 B B B について B t B_t B t と B t 2 − t B_t^2-t B t 2 − t がマルチンゲールであることを見ました。この節ではその逆を証明します。この節では時間を [ 0 , ∞ ) [0,\infty) [ 0 , ∞ ) とします。証明には停止時刻 (その時刻までの情報で「もう起きたか」を判定できるランダムな時刻)を使います。
補題 6.1 (止めた連続マルチンゲール ) ✓
( M t ) t ≥ 0 (M_t)_{t\ge0} ( M t ) t ≥ 0 を道が連続な ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲール、τ \tau τ を停止時刻とします。このとき ( M t ∧ τ ) t ≥ 0 (M_{t\wedge\tau})_{t\ge0} ( M t ∧ τ ) t ≥ 0 も ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲールです。
証明(補題 6.1) 適合性: t t t を固定し、σ n : = t ∧ ( ⌈ 2 n τ ⌉ / 2 n ) \sigma_n:=t\wedge(\lceil2^n\tau\rceil/2^n) σ n := t ∧ (⌈ 2 n τ ⌉ / 2 n ) とおくと、σ n \sigma_n σ n は t t t 以下の有限個の値(2 進有理数と t t t )しかとらず、各値 q < t q<t q < t について { σ n = q } = { q − 2 − n < τ ≤ q } ∈ F q ⊂ F t \{\sigma_n=q\}=\{q-2^{-n}<\tau\le q\}\in\mathcal{F}_q\subset\mathcal{F}_t { σ n = q } = { q − 2 − n < τ ≤ q } ∈ F q ⊂ F t です。よって M σ n M_{\sigma_n} M σ n は F t \mathcal{F}_t F t -可測で、σ n ↓ t ∧ τ \sigma_n\downarrow t\wedge\tau σ n ↓ t ∧ τ と道の連続性から M t ∧ τ = lim M σ n M_{t\wedge\tau}=\lim M_{\sigma_n} M t ∧ τ = lim M σ n も F t \mathcal{F}_t F t -可測です。可積分性は 定理 6.2[マルチンゲールとブラウン運動] (a) の主張に含まれます。
マルチンゲール性: s < t s<t s < t 、A ∈ F s A\in\mathcal{F}_s A ∈ F s とし、ρ : = τ ∧ s \rho:=\tau\wedge s ρ := τ ∧ s (A A A の上)、ρ : = τ ∧ t \rho:=\tau\wedge t ρ := τ ∧ t (A c A^c A c の上)とおきます。r < s r<s r < s なら { ρ ≤ r } = { τ ≤ r } \{\rho\le r\}=\{\tau\le r\} { ρ ≤ r } = { τ ≤ r } 、s ≤ r < t s\le r<t s ≤ r < t なら { ρ ≤ r } = A ∪ ( A c ∩ { τ ≤ r } ) \{\rho\le r\}=A\cup(A^c\cap\{\tau\le r\}) { ρ ≤ r } = A ∪ ( A c ∩ { τ ≤ r }) 、r ≥ t r\ge t r ≥ t なら { ρ ≤ r } = Ω \{\rho\le r\}=\Omega { ρ ≤ r } = Ω で、いずれも F r \mathcal{F}_r F r に属するので ρ \rho ρ は停止時刻です。定理 6.2[マルチンゲールとブラウン運動] (a) を停止時刻 ρ \rho ρ と τ \tau τ に(定数 t t t で)適用すると E [ M t ∧ ρ ] = E [ M 0 ] = E [ M t ∧ τ ] \mathbb{E}[M_{t\wedge\rho}]=\mathbb{E}[M_0]=\mathbb{E}[M_{t\wedge\tau}] E [ M t ∧ ρ ] = E [ M 0 ] = E [ M t ∧ τ ] です。t ∧ ρ t\wedge\rho t ∧ ρ は A A A の上で τ ∧ s \tau\wedge s τ ∧ s 、A c A^c A c の上で τ ∧ t \tau\wedge t τ ∧ t なので、差を取ると E [ 1 A ( M τ ∧ s − M τ ∧ t ) ] = 0 \mathbb{E}[\mathbf{1}_A(M_{\tau\wedge s}-M_{\tau\wedge t})]=0 E [ 1 A ( M τ ∧ s − M τ ∧ t )] = 0 を得ます。
∎
補題 6.2 (特性関数による独立性と分布の判定 ) ✓
Y Y Y を R d \mathbb{R}^d R d 値確率変数、G \mathcal{G} G を部分 σ-加法族、h ≥ 0 h\ge0 h ≥ 0 とし、すべての A ∈ G A\in\mathcal{G} A ∈ G と u ∈ R d u\in\mathbb{R}^d u ∈ R d について
E [ 1 A e i u ⋅ Y ] = P ( A ) e − ∣ u ∣ 2 h / 2 \mathbb{E}\bigl[\mathbf{1}_Ae^{iu\cdot Y}\bigr]=\mathbb{P}(A)\,e^{-|u|^2h/2} E [ 1 A e i u ⋅ Y ] = P ( A ) e − ∣ u ∣ 2 h /2 が成り立つとします。このとき Y Y Y は G \mathcal{G} G と独立で、N ( 0 , h I d ) N(0,hI_d) N ( 0 , h I d ) に従います。
証明(補題 6.2) A = Ω A=\Omega A = Ω とすると Y Y Y の特性関数は e − ∣ u ∣ 2 h / 2 e^{-|u|^2h/2} e − ∣ u ∣ 2 h /2 で、これは N ( 0 , h I d ) N(0,hI_d) N ( 0 , h I d ) の特性関数です(補題 2.3[特性関数と Gauss ベクトル] を成分ごとに使い、独立な成分の積を取ります)。定理 3.3[特性関数と Gauss ベクトル] から Y ∼ N ( 0 , h I d ) = : γ h Y\sim N(0,hI_d)=:\gamma_h Y ∼ N ( 0 , h I d ) =: γ h です。P ( A ) > 0 \mathbb{P}(A)>0 P ( A ) > 0 なる A ∈ G A\in\mathcal{G} A ∈ G については、確率測度 P ( ⋅ ∣ A ) \mathbb{P}(\cdot\mid A) P ( ⋅ ∣ A ) のもとでの Y Y Y の特性関数も仮定から e − ∣ u ∣ 2 h / 2 e^{-|u|^2h/2} e − ∣ u ∣ 2 h /2 なので、同じ定理により P ( Y ∈ B ∣ A ) = γ h ( B ) \mathbb{P}(Y\in B\mid A)=\gamma_h(B) P ( Y ∈ B ∣ A ) = γ 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 ∩ { Y ∈ B }) = P ( A ) P ( Y ∈ B ) で、P ( A ) = 0 \mathbb{P}(A)=0 P ( A ) = 0 のときも両辺は 0 0 0 です。
∎
次が証明の核心です。伊藤積分を使わず、e i λ x e^{i\lambda x} e iλ x の Taylor 展開とマルチンゲール差の直交性だけを使います。
証明の流れは次のとおりです。過程がまだ大きく動いていない間だけ追いかけ、増分が ± K \pm K ± K の帯から出た時刻 τ \tau τ で止めると、止めた過程は常に [ − K , K ] [-K,K] [ − K , K ] に収まります。有界になれば Taylor 展開の剰余を評価できるので、各小区間で 3 次の項を無視してよいことが示せ、止めた過程の特性関数が正規分布のものになります。最後に K → ∞ K\to\infty K → ∞ として停止を外します。
Lévy の特徴付けの証明で使う停止。増分 N_r − N_s が ±K の帯を出た時刻 τ で止めると、止めた過程は有界になり、分割 r_0 < r_1 < … の各区間で 3 次の Taylor 展開が使える
補題 6.3 (Lévy の特徴付けの 1 次元の核 ) ✓
( N t ) t ≥ 0 (N_t)_{t\ge0} ( N t ) t ≥ 0 を道が連続な ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲールで、各 N t N_t N t が二乗可積分、かつ N t 2 − t N_t^2-t N t 2 − t も ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲールであるとします。このとき 0 ≤ s < t 0\le s<t 0 ≤ s < t 、A ∈ F s A\in\mathcal{F}_s A ∈ F s 、λ ∈ R \lambda\in\mathbb{R} λ ∈ R について
E [ 1 A e i λ ( N t − N s ) ] = 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}. E [ 1 A e iλ ( N t − N s ) ] = P ( A ) e − λ 2 ( t − s ) /2 .
証明(補題 6.3) 停止による有界化。 K > 0 K>0 K > 0 とし、τ : = inf { r ≥ s : ∣ N r − N s ∣ ≥ K } ∧ t \tau:=\inf\{r\ge s:|N_r-N_s|\ge K\}\wedge t τ := inf { r ≥ s : ∣ N r − N s ∣ ≥ K } ∧ t とおきます。道の連続性から、s ≤ r < t s\le r<t s ≤ r < t について { τ ≤ r } = ⋂ m ⋃ q { ∣ N q − N s ∣ > K − 1 / m } \{\tau\le r\}=\bigcap_{m}\bigcup_{q}\{|N_q-N_s|>K-1/m\} { τ ≤ r } = ⋂ m ⋃ q { ∣ N q − N s ∣ > K − 1/ m } (q q q は [ s , r ] [s,r] [ s , r ] の有理数と r r r を動く)は F r \mathcal{F}_r F r に属し、τ \tau τ は停止時刻です。r ∈ [ s , t ] r\in[s,t] r ∈ [ s , t ] について X r : = N r ∧ τ − N s X_r:=N_{r\wedge\tau}-N_s X r := N r ∧ τ − N s とおくと ∣ X r ∣ ≤ K |X_r|\le K ∣ X r ∣ ≤ K です。補題 6.1 を N N N と N 2 − t N^2-t N 2 − t に適用すると、s ≤ r ≤ r ′ ≤ t s\le r\le r'\le t s ≤ r ≤ r ′ ≤ t について Δ : = X r ′ − X r \Delta:=X_{r'}-X_r Δ := X r ′ − X r 、δ : = r ′ ∧ τ − r ∧ τ ∈ [ 0 , r ′ − r ] \delta:=r'\wedge\tau-r\wedge\tau\in[0,r'-r] δ := r ′ ∧ τ − r ∧ τ ∈ [ 0 , r ′ − r ] は
E [ Δ ∣ F r ] = 0 , E [ Δ 2 − δ ∣ F r ] = 0 \mathbb{E}[\Delta\mid\mathcal{F}_r]=0,\qquad\mathbb{E}[\Delta^2-\delta\mid\mathcal{F}_r]=0 E [ Δ ∣ F r ] = 0 , E [ Δ 2 − δ ∣ F r ] = 0 を満たします。後者は Δ 2 = N r ′ ∧ τ 2 − N r ∧ τ 2 − 2 N r ∧ τ Δ \Delta^2=N_{r'\wedge\tau}^2-N_{r\wedge\tau}^2-2N_{r\wedge\tau}\Delta Δ 2 = N r ′ ∧ τ 2 − N r ∧ τ 2 − 2 N r ∧ τ Δ と、N r ∧ τ = N s + X r N_{r\wedge\tau}=N_s+X_r N r ∧ τ = N s + X r が二乗可積分であることから 定理 4.4[条件付き期待値] で従います。
離散的な展開。 Φ ( r ) : = exp ( i λ X r + 1 2 λ 2 ( r ∧ τ − s ) ) \Phi(r):=\exp\bigl(i\lambda X_r+\tfrac12\lambda^2(r\wedge\tau-s)\bigr) Φ ( r ) := exp ( iλ X r + 2 1 λ 2 ( r ∧ τ − s ) ) とおくと、Φ ( r ) \Phi(r) Φ ( r ) は F r \mathcal{F}_r F r -可測で ∣ Φ ( r ) ∣ ≤ e λ 2 ( t − s ) / 2 |\Phi(r)|\le e^{\lambda^2(t-s)/2} ∣Φ ( r ) ∣ ≤ e λ 2 ( t − s ) /2 です。分割 s = r 0 < r 1 < ⋯ < r m = t s=r_0<r_1<\dots<r_m=t s = r 0 < r 1 < ⋯ < r m = t を取り、Δ k = X r k + 1 − X r k \Delta_k=X_{r_{k+1}}-X_{r_k} Δ k = X r k + 1 − X r k 、δ k = r k + 1 ∧ τ − r k ∧ τ \delta_k=r_{k+1}\wedge\tau-r_k\wedge\tau δ k = r k + 1 ∧ τ − r k ∧ τ とおくと Φ ( r k + 1 ) = Φ ( r k ) g ( Δ k , δ k ) \Phi(r_{k+1})=\Phi(r_k)\,g(\Delta_k,\delta_k) Φ ( r k + 1 ) = Φ ( r k ) g ( Δ k , δ k ) 、g ( x , y ) = e i λ x + λ 2 y / 2 g(x,y)=e^{i\lambda x+\lambda^2y/2} g ( x , y ) = e iλ x + λ 2 y /2 です。w = i λ x + 1 2 λ 2 y w=i\lambda x+\frac12\lambda^2y w = iλ x + 2 1 λ 2 y について e w = 1 + w + 1 2 w 2 + R e^w=1+w+\frac12w^2+R e w = 1 + w + 2 1 w 2 + R 、∣ R ∣ ≤ 1 6 ∣ w ∣ 3 e ∣ w ∣ |R|\le\frac16|w|^3e^{|w|} ∣ R ∣ ≤ 6 1 ∣ w ∣ 3 e ∣ w ∣ で、1 2 w 2 = − 1 2 λ 2 x 2 + i 2 λ 3 x y + 1 8 λ 4 y 2 \frac12w^2=-\frac12\lambda^2x^2+\frac{i}{2}\lambda^3xy+\frac18\lambda^4y^2 2 1 w 2 = − 2 1 λ 2 x 2 + 2 i λ 3 x y + 8 1 λ 4 y 2 ですから
g ( x , y ) = 1 + i λ x + 1 2 λ 2 ( y − x 2 ) + ρ ( x , y ) , ∣ ρ ( x , y ) ∣ ≤ c ( ∣ x ∣ 3 + ∣ x ∣ y + y 2 ) 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) g ( x , y ) = 1 + iλ x + 2 1 λ 2 ( y − x 2 ) + ρ ( x , y ) , ∣ ρ ( x , y ) ∣ ≤ c ( ∣ x ∣ 3 + ∣ x ∣ y + y 2 ) が ∣ x ∣ ≤ 2 K |x|\le2K ∣ x ∣ ≤ 2 K 、0 ≤ y ≤ t − s 0\le y\le t-s 0 ≤ y ≤ t − s の範囲で成り立ちます(c c c は λ , K , t − s \lambda,K,t-s λ , K , t − s だけで決まる定数。∣ w ∣ 3 ≤ 4 ( ∣ λ x ∣ 3 + 1 8 λ 6 y 3 ) |w|^3\le4(|\lambda x|^3+\frac18\lambda^6y^3) ∣ w ∣ 3 ≤ 4 ( ∣ λ x ∣ 3 + 8 1 λ 6 y 3 ) と y 3 ≤ ( t − s ) y 2 y^3\le(t-s)y^2 y 3 ≤ ( t − s ) y 2 を使います)。1 A Φ ( r k ) \mathbf{1}_A\Phi(r_k) 1 A Φ ( r k ) は有界で F r k \mathcal{F}_{r_k} F r k -可測なので、上の二つの条件付き期待値から
E [ 1 A Φ ( r k + 1 ) ] − E [ 1 A Φ ( r k ) ] = E [ 1 A Φ ( r k ) ρ ( Δ 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]. E [ 1 A Φ ( r k + 1 )] − E [ 1 A Φ ( r k )] = E [ 1 A Φ ( r k ) ρ ( Δ k , δ k ) ] . k k k について足すと、Φ ( r 0 ) = 1 \Phi(r_0)=1 Φ ( r 0 ) = 1 より
∣ E [ 1 A Φ ( t ) ] − P ( A ) ∣ ≤ c e λ 2 ( t − s ) / 2 E [ ∑ k ( ∣ Δ k ∣ 3 + ∣ Δ k ∣ δ k + δ k 2 ) ] . \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]. E [ 1 A Φ ( t )] − P ( A ) ≤ c e λ 2 ( t − s ) /2 E [ k ∑ ( ∣ Δ k ∣ 3 + ∣ Δ k ∣ δ k + δ k 2 ) ] . 右辺が消えること。 ∥ Π ∥ = max k ( r k + 1 − r k ) \|\Pi\|=\max_k(r_{k+1}-r_k) ∥Π∥ = max k ( r k + 1 − r k ) とします。まず ∑ k δ k ≤ t − s \sum_k\delta_k\le t-s ∑ k δ k ≤ t − s と E ∑ k Δ k 2 = E ∑ k δ k ≤ t − s \mathbb{E}\sum_k\Delta_k^2=\mathbb{E}\sum_k\delta_k\le t-s E ∑ k Δ k 2 = E ∑ k δ k ≤ t − s です。
∑ k δ k 2 ≤ ∥ Π ∥ ( t − s ) \sum_k\delta_k^2\le\|\Pi\|(t-s) ∑ k δ k 2 ≤ ∥Π∥ ( t − s ) 。
Cauchy–Schwarz の不等式から E ∑ k ∣ Δ k ∣ δ k ≤ ( E ∑ k Δ k 2 ) 1 / 2 ( E ∑ k δ k 2 ) 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} E ∑ k ∣ Δ k ∣ δ k ≤ ( E ∑ k Δ k 2 ) 1/2 ( E ∑ k δ k 2 ) 1/2 ≤ ( t − s ) ∥Π ∥ 1/2 。
∑ k ∣ Δ k ∣ 3 ≤ max k ∣ Δ k ∣ ⋅ Q \sum_k|\Delta_k|^3\le\max_k|\Delta_k|\cdot Q ∑ k ∣ Δ k ∣ 3 ≤ max k ∣ Δ k ∣ ⋅ Q 、Q : = ∑ k Δ k 2 Q:=\sum_k\Delta_k^2 Q := ∑ k Δ k 2 なので E ∑ k ∣ Δ k ∣ 3 ≤ ( E max k Δ k 2 ) 1 / 2 ( E Q 2 ) 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} E ∑ k ∣ Δ k ∣ 3 ≤ ( E max k Δ k 2 ) 1/2 ( E Q 2 ) 1/2 。ここで E Q 2 \mathbb{E}Q^2 E Q 2 は分割によらず有界です。実際 Q 2 = ∑ k Δ k 4 + 2 ∑ k Δ k 2 ∑ l > k Δ l 2 Q^2=\sum_k\Delta_k^4+2\sum_k\Delta_k^2\sum_{l>k}\Delta_l^2 Q 2 = ∑ k Δ k 4 + 2 ∑ k Δ k 2 ∑ l > k Δ l 2 で、Δ k 4 ≤ 4 K 2 Δ k 2 \Delta_k^4\le4K^2\Delta_k^2 Δ k 4 ≤ 4 K 2 Δ k 2 、また l ≠ l ′ l\ne l' l = l ′ のとき E [ Δ l Δ l ′ ∣ F r k + 1 ] = 0 \mathbb{E}[\Delta_l\Delta_{l'}\mid\mathcal{F}_{r_{k+1}}]=0 E [ Δ l Δ l ′ ∣ F r k + 1 ] = 0 (l , l ′ > k l,l'>k l , l ′ > k )なので E [ ∑ l > k Δ l 2 ∣ F r k + 1 ] = E [ ( X t − X r k + 1 ) 2 ∣ F r k + 1 ] ≤ 4 K 2 \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 E [ ∑ l > k Δ l 2 ∣ F r k + 1 ] = E [( X t − X r k + 1 ) 2 ∣ F r k + 1 ] ≤ 4 K 2 です。よって E Q 2 ≤ 4 K 2 ( t − s ) + 8 K 2 ( t − s ) = 12 K 2 ( t − s ) \mathbb{E}Q^2\le4K^2(t-s)+8K^2(t-s)=12K^2(t-s) E Q 2 ≤ 4 K 2 ( t − s ) + 8 K 2 ( t − s ) = 12 K 2 ( t − s ) 。一方 max k Δ k 2 ≤ 4 K 2 \max_k\Delta_k^2\le4K^2 max k Δ k 2 ≤ 4 K 2 で、道 r ↦ X r r\mapsto X_r r ↦ X r の一様連続性から各 ω \omega ω で max k ∣ Δ k ∣ → 0 \max_k|\Delta_k|\to0 max k ∣ Δ k ∣ → 0 なので、優収束定理により E max k Δ k 2 → 0 \mathbb{E}\max_k\Delta_k^2\to0 E max k Δ k 2 → 0 。
したがって ∥ Π ∥ → 0 \|\Pi\|\to0 ∥Π∥ → 0 とすると E [ 1 A Φ ( t ) ] = P ( A ) \mathbb{E}[\mathbf{1}_A\Phi(t)]=\mathbb{P}(A) E [ 1 A Φ ( t )] = P ( A ) です。
停止を外す。 K K K を K → ∞ K\to\infty K → ∞ とします。sup r ∈ [ s , t ] ∣ N r − N s ∣ < K \sup_{r\in[s,t]}|N_r-N_s|<K sup r ∈ [ s , t ] ∣ N r − N s ∣ < K となる ω \omega ω では τ = t \tau=t τ = t なので、各 ω \omega ω で K K K が十分大きければ Φ ( t ) = e i λ ( N t − N s ) + λ 2 ( t − s ) / 2 \Phi(t)=e^{i\lambda(N_t-N_s)+\lambda^2(t-s)/2} Φ ( t ) = e iλ ( N t − N s ) + λ 2 ( t − s ) /2 です。∣ Φ ( t ) ∣ ≤ e λ 2 ( t − s ) / 2 |\Phi(t)|\le e^{\lambda^2(t-s)/2} ∣Φ ( t ) ∣ ≤ e λ 2 ( t − s ) /2 なので、優収束定理から E [ 1 A e i λ ( N t − N s ) ] 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) E [ 1 A e iλ ( N t − N s ) ] e λ 2 ( t − s ) /2 = P ( A ) を得ます。
∎
定理 6.4 (Lévy の特徴付け ) ✓
M = ( M 1 , … , M d ) M=(M^1,\dots,M^d) M = ( M 1 , … , M d ) を ( F t ) (\mathcal{F}_t) ( F t ) に適合した R d \mathbb{R}^d R d 値の過程で、M 0 = 0 M_0=0 M 0 = 0 、すべての道が連続、各 M i M^i M i が二乗可積分な ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲールであり、さらにすべての i , j i,j i , j について M t i M t j − δ i j t M^i_tM^j_t-\delta_{ij}t M t i M t j − δ ij t が ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲールであるとします(すなわち ⟨ M i , M j ⟩ t = δ i j t \langle M^i,M^j\rangle_t=\delta_{ij}t ⟨ M i , M j ⟩ t = δ ij t )。このとき M M M は d d d 次元 ( F t ) (\mathcal{F}_t) ( F t ) -Brown 運動です。
証明(定理 6.4) 定義 2.1 の条件 (1)(2) は仮定です。条件 (3) を示します。u ∈ R d u\in\mathbb{R}^d u ∈ R d , u ≠ 0 u\ne0 u = 0 を固定し、N t : = u ⋅ M t / ∣ u ∣ N_t:=u\cdot M_t/|u| N t := u ⋅ M t /∣ u ∣ とおきます。N N N は道が連続な二乗可積分マルチンゲールで、
N t 2 − t = 1 ∣ u ∣ 2 ∑ i , j u i u j ( M t i M t j − δ i j t ) N_t^2-t=\frac{1}{|u|^2}\sum_{i,j}u_iu_j\bigl(M^i_tM^j_t-\delta_{ij}t\bigr) N t 2 − t = ∣ u ∣ 2 1 i , j ∑ u i u j ( M t i M t j − δ ij t ) はマルチンゲールの線形結合なのでマルチンゲールです。補題 6.3 を λ = ∣ u ∣ \lambda=|u| λ = ∣ u ∣ で適用すると、s < t s<t s < t 、A ∈ F s A\in\mathcal{F}_s A ∈ F s について
E [ 1 A e i u ⋅ ( M t − M s ) ] = 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} E [ 1 A e i u ⋅ ( M t − M s ) ] = P ( A ) e − ∣ u ∣ 2 ( t − s ) /2 で、u = 0 u=0 u = 0 でも自明に成り立ちます。補題 6.2 から M t − M s M_t-M_s M t − M s は F s \mathcal{F}_s F s と独立で N ( 0 , ( t − s ) I d ) N(0,(t-s)I_d) N ( 0 , ( t − s ) I d ) に従います。
∎
連続性の仮定は外せません。演習 7.4 で見るように、跳びのあるマルチンゲールでは M t 2 − t M_t^2-t M t 2 − t がマルチンゲールでも Brown 運動にはなりません。証明の中で連続性を使ったのは、停止した過程が有界になる箇所と、max k ∣ Δ k ∣ → 0 \max_k|\Delta_k|\to0 max k ∣ Δ k ∣ → 0 の箇所です。
例 6.5 (符号で折り返した積分 )
d = 1 d=1 d = 1 とし、sgn ( x ) = 1 \operatorname{sgn}(x)=1 sgn ( x ) = 1 (x > 0 x>0 x > 0 )、− 1 -1 − 1 (x ≤ 0 x\le0 x ≤ 0 )として β t : = ∫ 0 t sgn ( W u ) d W u \beta_t:=\int_0^t\operatorname{sgn}(W_u)\,dW_u β t := ∫ 0 t sgn ( W u ) d W u とおきます。sgn ( W u ) \operatorname{sgn}(W_u) sgn ( W u ) は連続な適合過程のボレル関数なので発展的可測で、有界です。β \beta β は連続なマルチンゲールで、命題 3.2 から β t 2 − ∫ 0 t sgn ( W u ) 2 d u = β t 2 − t \beta_t^2-\int_0^t\operatorname{sgn}(W_u)^2\,du=\beta_t^2-t β t 2 − ∫ 0 t sgn ( W u ) 2 d u = β t 2 − t もマルチンゲールです(任意の T T T で成り立つので [ 0 , ∞ ) [0,\infty) [ 0 , ∞ ) 上で成り立ちます)。定理 6.4 により β \beta β は Brown 運動です。β \beta β は W W W の符号を見て向きを変えた過程なので、見かけ上 W W W と強く依存していますが、分布は Brown 運動そのものです。
例 6.6 (ランダムな回転 )
O u O_u O u を発展的可測な過程で、各 ( u , ω ) (u,\omega) ( u , ω ) で d × d d\times d d × d の直交行列であるものとし、W ~ t : = ∫ 0 t O u d W u \tilde W_t:=\int_0^tO_u\,dW_u W ~ t := ∫ 0 t O u d W u (成分は W ~ i = ∑ a ∫ O i a d W a \tilde W^i=\sum_a\int O^{ia}\,dW^a W ~ i = ∑ a ∫ O ia d W a )とおきます。命題 3.2 から W ~ i W ~ j − ∫ 0 t ( O u O u T ) i j d u = W ~ i W ~ j − δ i j t \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 W ~ i W ~ j − ∫ 0 t ( O u O u T ) ij d u = W ~ i W ~ j − δ ij t はマルチンゲールで、定理 6.4 により W ~ \tilde W W ~ は d d d 次元 Brown 運動です。回転を時刻ごとに、過去の道に応じて選び直しても、Brown 運動の等方性は崩れません。
演習 7.1 易
d = 2 d=2 d = 2 とし、W t 1 W t 2 W^1_tW^2_t W t 1 W t 2 を 系 5.3 で確率積分として表し、それを使って E [ ( W t 1 W t 2 ) 2 ] = t 2 \mathbb{E}[(W^1_tW^2_t)^2]=t^2 E [( W t 1 W t 2 ) 2 ] = t 2 を示してください。
解答 ⟨ W 1 , W 2 ⟩ = 0 \langle W^1,W^2\rangle=0 ⟨ W 1 , W 2 ⟩ = 0 なので W t 1 W t 2 = ∫ 0 t W u 2 d W u 1 + ∫ 0 t W u 1 d W u 2 = ∫ 0 t H u ⋅ d W u W^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 W t 1 W t 2 = ∫ 0 t W u 2 d W u 1 + ∫ 0 t W u 1 d W u 2 = ∫ 0 t H u ⋅ d W u 、H = ( W 2 , W 1 ) H=(W^2,W^1) H = ( W 2 , W 1 ) です。補題 3.1 により
E [ ( W t 1 W t 2 ) 2 ] = E ∫ 0 t ( ( W u 2 ) 2 + ( W u 1 ) 2 ) d u = ∫ 0 t 2 u d u = t 2 . \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 [( W t 1 W t 2 ) 2 ] = E ∫ 0 t ( ( W u 2 ) 2 + ( W u 1 ) 2 ) d u = ∫ 0 t 2 u d u = t 2 . 独立性から直接 E [ ( W t 1 ) 2 ] E [ ( W t 2 ) 2 ] = t ⋅ t \mathbb{E}[(W^1_t)^2]\,\mathbb{E}[(W^2_t)^2]=t\cdot t E [( W t 1 ) 2 ] E [( W t 2 ) 2 ] = t ⋅ t と計算した値と一致します。
演習 7.2 標準
a ∈ R a\in\mathbb{R} a ∈ R とし、1 次元 Brown 運動 W W W について m ( t ) : = E [ W t 2 cos ( a W t ) ] m(t):=\mathbb{E}[W_t^2\cos(aW_t)] m ( t ) := E [ W t 2 cos ( a W t )] とおきます。g ( x ) = x 2 cos ( a x ) g(x)=x^2\cos(ax) g ( x ) = x 2 cos ( a x ) は導関数が有界でないので、定理 5.1 をそのままは使えません。χ ∈ C ∞ ( R ) \chi\in C^\infty(\mathbb{R}) χ ∈ C ∞ ( R ) を ∣ x ∣ ≤ 1 |x|\le1 ∣ x ∣ ≤ 1 で 1 1 1 、∣ x ∣ ≥ 2 |x|\ge2 ∣ x ∣ ≥ 2 で 0 0 0 、0 ≤ χ ≤ 1 0\le\chi\le1 0 ≤ χ ≤ 1 なる関数とし、g N ( x ) = g ( x ) χ ( x / N ) g_N(x)=g(x)\chi(x/N) g N ( x ) = g ( x ) χ ( x / N ) に伊藤公式を適用して N → ∞ N\to\infty N → ∞ とすることで
m ′ ( t ) = E [ cos ( a W t ) ] − 2 a E [ W t sin ( a W t ) ] − a 2 2 m ( t ) m'(t)=\mathbb{E}[\cos(aW_t)]-2a\,\mathbb{E}[W_t\sin(aW_t)]-\frac{a^2}{2}m(t) m ′ ( t ) = E [ cos ( a W t )] − 2 a E [ W t sin ( a W t )] − 2 a 2 m ( t ) を導いてください。さらに m ( t ) = ( t − a 2 t 2 ) e − a 2 t / 2 m(t)=(t-a^2t^2)e^{-a^2t/2} m ( t ) = ( t − a 2 t 2 ) e − a 2 t /2 がこれを満たすことを確かめてください。
解答 g N g_N g N は導関数が有界なので、定理 5.1 (X = W X=W X = W 、μ = 0 \mu=0 μ = 0 、σ = 1 \sigma=1 σ = 1 )と確率積分の平均が 0 0 0 であることから E [ g N ( W t ) ] = 1 2 ∫ 0 t E [ g N ′ ′ ( W u ) ] d u \mathbb{E}[g_N(W_t)]=\frac12\int_0^t\mathbb{E}[g_N''(W_u)]\,du E [ g N ( W t )] = 2 1 ∫ 0 t E [ g N ′′ ( W u )] d u です。g N ′ ′ = g ′ ′ χ ( ⋅ / N ) + 2 g ′ χ ′ ( ⋅ / N ) / N + g χ ′ ′ ( ⋅ / N ) / N 2 g_N''=g''\chi(\cdot/N)+2g'\chi'(\cdot/N)/N+g\chi''(\cdot/N)/N^2 g N ′′ = g ′′ χ ( ⋅ / N ) + 2 g ′ χ ′ ( ⋅ / N ) / N + g χ ′′ ( ⋅ / N ) / N 2 で、χ ′ ( ⋅ / N ) \chi'(\cdot/N) χ ′ ( ⋅ / N ) と χ ′ ′ ( ⋅ / N ) \chi''(\cdot/N) χ ′′ ( ⋅ / N ) は N ≤ ∣ x ∣ ≤ 2 N N\le|x|\le2N N ≤ ∣ x ∣ ≤ 2 N の外では 0 0 0 なので、∣ g ∣ ≤ x 2 |g|\le x^2 ∣ g ∣ ≤ x 2 、∣ g ′ ∣ ≤ 2 ∣ x ∣ + ∣ a ∣ x 2 |g'|\le2|x|+|a|x^2 ∣ g ′ ∣ ≤ 2∣ x ∣ + ∣ a ∣ x 2 から ∣ g N ∣ + ∣ g N ′ ′ ∣ ≤ c ( 1 + x 2 ) |g_N|+|g_N''|\le c(1+x^2) ∣ g N ∣ + ∣ g N ′′ ∣ ≤ c ( 1 + x 2 ) (c c c は N N N によらない)です。E [ W u 2 ] = u \mathbb{E}[W_u^2]=u E [ W u 2 ] = u なので、優収束定理により N → ∞ N\to\infty N → ∞ で m ( t ) = 1 2 ∫ 0 t E [ g ′ ′ ( W u ) ] d u m(t)=\frac12\int_0^t\mathbb{E}[g''(W_u)]\,du m ( t ) = 2 1 ∫ 0 t E [ g ′′ ( W u )] d u を得ます。g ′ ′ ( x ) = 2 cos ( a x ) − 4 a x sin ( a x ) − a 2 x 2 cos ( a x ) g''(x)=2\cos(ax)-4ax\sin(ax)-a^2x^2\cos(ax) g ′′ ( x ) = 2 cos ( a x ) − 4 a x sin ( a x ) − a 2 x 2 cos ( a x ) を代入して t t t で微分すると求める式です。
例 5.5 から E [ cos ( a W t ) ] = e − a 2 t / 2 \mathbb{E}[\cos(aW_t)]=e^{-a^2t/2} E [ cos ( a W t )] = e − a 2 t /2 、これを a a a で微分して E [ W t sin ( a W t ) ] = a t e − a 2 t / 2 \mathbb{E}[W_t\sin(aW_t)]=at\,e^{-a^2t/2} E [ W t sin ( a W t )] = a t e − a 2 t /2 です(微分と期待値の交換は優収束定理による)。よって m ′ ( t ) = ( 1 − 2 a 2 t ) e − a 2 t / 2 − a 2 2 m ( t ) m'(t)=(1-2a^2t)e^{-a^2t/2}-\frac{a^2}{2}m(t) m ′ ( t ) = ( 1 − 2 a 2 t ) e − a 2 t /2 − 2 a 2 m ( t ) で、m ( t ) = ( t − a 2 t 2 ) e − a 2 t / 2 m(t)=(t-a^2t^2)e^{-a^2t/2} m ( t ) = ( t − a 2 t 2 ) e − a 2 t /2 を代入すると左辺は ( 1 − 2 a 2 t ) e − a 2 t / 2 − a 2 2 ( t − a 2 t 2 ) e − a 2 t / 2 (1-2a^2t)e^{-a^2t/2}-\frac{a^2}{2}(t-a^2t^2)e^{-a^2t/2} ( 1 − 2 a 2 t ) e − a 2 t /2 − 2 a 2 ( t − a 2 t 2 ) e − a 2 t /2 となり右辺に一致します。m ( 0 ) = 0 m(0)=0 m ( 0 ) = 0 なので、1 階線形常微分方程式の解の一意性からこれが m m m です。検算として、m ( t ) = − d 2 d a 2 e − a 2 t / 2 m(t)=-\frac{d^2}{da^2}e^{-a^2t/2} m ( t ) = − d a 2 d 2 e − a 2 t /2 からも同じ式が出ます。
演習 7.3 標準
d = 2 d=2 d = 2 とし、θ u : = W u 1 \theta_u:=W^1_u θ u := W u 1 として
B t = ∫ 0 t ( cos θ u d W u 1 + sin θ u d W u 2 ) , B t ′ = ∫ 0 t ( − sin θ u d W u 1 + cos θ u d W u 2 ) 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 t = ∫ 0 t ( cos θ u d W u 1 + sin θ u d W u 2 ) , B t ′ = ∫ 0 t ( − sin θ u d W u 1 + cos θ u d W u 2 ) とおきます。( B , B ′ ) (B,B') ( B , B ′ ) が 2 次元 Brown 運動であることを示してください。
解答 O u = ( cos θ u sin θ u − sin θ u cos θ u ) O_u=\begin{pmatrix}\cos\theta_u&\sin\theta_u\\-\sin\theta_u&\cos\theta_u\end{pmatrix} O u = ( cos θ u − sin θ u sin θ u cos θ u ) は各 ( u , ω ) (u,\omega) ( u , ω ) で直交行列で、θ u = W u 1 \theta_u=W^1_u θ u = W u 1 が連続な適合過程なので発展的可測です。( B , B ′ ) = ∫ O d W (B,B')=\int O\,dW ( B , B ′ ) = ∫ O d W なので 例 6.6 から従います。直接には、命題 3.2 により B 2 − t B^2-t B 2 − t 、B ′ 2 − t B'^2-t B ′2 − t (cos 2 + sin 2 = 1 \cos^2+\sin^2=1 cos 2 + sin 2 = 1 )と B B ′ BB' B B ′ (− cos θ sin θ + sin θ cos θ = 0 -\cos\theta\sin\theta+\sin\theta\cos\theta=0 − cos θ sin θ + sin θ cos θ = 0 )がマルチンゲールで、定理 6.4 を使います。
演習 7.4 難
( P t ) (P_t) ( P t ) を強度 1 1 1 の Poisson 過程とし、その自然なフィルトレーションを ( F t ) (\mathcal{F}_t) ( F t ) とします。M t : = P t − t M_t:=P_t-t M t := P t − t について、M t M_t M t と M t 2 − t M_t^2-t M t 2 − t がどちらも ( F t ) (\mathcal{F}_t) ( F t ) -マルチンゲールであることを示し、それでも M M M が Brown 運動でない理由を述べてください。定理 6.4 のどの仮定が破れていますか。
解答 定理 4.5[Markov 連鎖と Poisson 過程] により P t − P s P_t-P_s P t − P s は F s \mathcal{F}_s F s と独立で平均 t − s t-s t − s 、分散 t − s t-s t − s の Poisson 分布に従います。よって E [ M t − M s ∣ F s ] = 0 \mathbb{E}[M_t-M_s\mid\mathcal{F}_s]=0 E [ M t − M s ∣ F s ] = 0 、また
E [ M t 2 − M s 2 ∣ F s ] = E [ ( M t − M s ) 2 ∣ F s ] + 2 M s E [ M t − M s ∣ F s ] = 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 E [ M t 2 − M s 2 ∣ F s ] = E [( M t − M s ) 2 ∣ F s ] + 2 M s E [ M t − M s ∣ F s ] = t − s なので M t 2 − t M_t^2-t M t 2 − t もマルチンゲールです。しかし M t + t = P t M_t+t=P_t M t + t = P t は整数値なので M t M_t M t は正規分布に従いません。破れているのは道の連続性の仮定で、P P P は高さ 1 1 1 の跳びを持ちます。補題 6.3 の証明では max k ∣ Δ k ∣ → 0 \max_k|\Delta_k|\to0 max k ∣ Δ k ∣ → 0 を使いましたが、跳びがあるとこの量は 1 1 1 以上のまま残り、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 — 伊藤の公式と確率積分の章。
執筆: Claude Code(Opus 5.5)、Claude Code(Sonnet 5.5)・監査 3 回(Claude Code)・最終 2026-10-12