18 世紀のフランスで流行した賭け方に、こういうものがあります。まず 1 単位を賭ける。負けたら次は 2 単位、また負けたら 4 単位、と倍々に増やしていく。勝率が 1 / 2 1/2 1/2 の公平なゲームなら、いつかは必ず勝つ。そして最初に勝った時点で、それまでの損失 1 + 2 + ⋯ + 2 k − 1 = 2 k − 1 1 + 2 + \cdots + 2^{k-1} = 2^k - 1 1 + 2 + ⋯ + 2 k − 1 = 2 k − 1 に対して 2 k 2^k 2 k の払い戻しを受けるので、差し引き必ず 1 単位の利益が出る。この賭け方は当時のフランス語で martingale と呼ばれました。今日の数学用語「マルチンゲール」の語源です。
もちろんこれは必勝法ではありません。破綻する理由は「所持金が有限」だからですが、では所持金が無限にあれば本当に必勝なのでしょうか。あるいは、もっと巧妙な賭け方なら公平なゲームを有利にできるのでしょうか。この素朴な疑問に完全な答えを与えるのが、この章で扱う理論です。答えを先に言うと、過去の情報だけを見て賭け金を決めるかぎり、有界な時間内では絶対に有利にできません (Theorem 3.7 と Theorem 4.3 )。そして倍賭け戦略が「必勝」に見えたのは、停止時刻 τ \tau τ が有界でも可積分でもなかったからです。
もう一つの出発点は物理です。1827 年、植物学者ロバート・ブラウンは水に浮かべた花粉から出た微粒子が不規則に運動し続けるのを観察しました。1905 年にアインシュタインがこれを分子の熱運動による衝突として説明し、変位の分散が時間に比例することを導きます。数学的には、これに先立つ 1900 年にバシュリエが株価の変動モデルとして本質的に同じ過程を扱っていました。
この二つの流れ、賭けの理論と物理の理論は、20 世紀半ばに一つになりました。ブラウン運動はマルチンゲールの典型例であり、逆にレヴィの特徴づけによれば、連続なマルチンゲールでその二次変分が t t t であるものはブラウン運動しかありません(Remark 6.5 )。この章はその合流点を見るための章です。
この章では、時間とともに増えていく情報を σ \sigma σ -加法族の増大列で表します。以下 ( Ω , F , P ) (\Omega, \mathcal{F}, P) ( Ω , F , P ) は確率空間とします(確率空間とコルモゴロフの公理 の Definition 4.1[Probability Spaces and Kolmogorov's Axioms] を参照)。
Definition 2.1 (フィルトレーションと適合過程 )
添字集合を T = { 0 , 1 , 2 , … } T = \{0, 1, 2, \ldots\} T = { 0 , 1 , 2 , … } (離散時間)または T = [ 0 , ∞ ) T = [0, \infty) T = [ 0 , ∞ ) (連続時間)とする。
F \mathcal{F} F の部分 σ \sigma σ -加法族の族 ( F t ) t ∈ T (\mathcal{F}_t)_{t \in T} ( F t ) t ∈ T が、s ≤ t s \le t s ≤ t ならば F s ⊂ F t \mathcal{F}_s \subset \mathcal{F}_t F s ⊂ F t を満たすとき、( F t ) t ∈ T (\mathcal{F}_t)_{t\in T} ( F t ) t ∈ T をフィルトレーション という。
確率過程 ( X t ) t ∈ T (X_t)_{t \in T} ( X t ) t ∈ T が、すべての t ∈ T t \in T t ∈ T について X t X_t X t が F t \mathcal{F}_t F t -可測であるとき、( X t ) (X_t) ( X t ) は ( F t ) (\mathcal{F}_t) ( F t ) に適合している という。
特に F t X : = σ ( X s : s ≤ t , s ∈ T ) \mathcal{F}_t^X := \sigma(X_s : s \le t, \ s\in T) F t X := σ ( X s : s ≤ t , s ∈ T ) を ( X t ) (X_t) ( X t ) の自然なフィルトレーション という。
F t \mathcal{F}_t F t は「時刻 t t t までに観測できた出来事の全体」を表します。適合性は「X t X_t X t の値は時刻 t t t には確定している」ことを意味し、未来を先読みしないという要請です。
証明で繰り返し使う条件付き期待値の性質を挙げておきます。Z , W Z, W Z , W を可積分な確率変数、G ⊂ F \mathcal{G} \subset \mathcal{F} G ⊂ F を部分 σ \sigma σ -加法族とします(詳細は 条件付き期待値 の Definition 3.1[条件付き期待値] 以降を参照)。
名称 主張 線形性 E [ a Z + b W ∣ G ] = a E [ Z ∣ G ] + b E [ W ∣ G ] E[aZ + bW \mid \mathcal{G}] = aE[Z\mid\mathcal{G}] + bE[W\mid\mathcal{G}] E [ a Z + bW ∣ G ] = a E [ Z ∣ G ] + b E [ W ∣ G ] a.s.タワー性 H ⊂ G \mathcal{H} \subset \mathcal{G} H ⊂ G なら E [ E [ Z ∣ G ] ∣ H ] = E [ Z ∣ H ] E\bigl[E[Z\mid\mathcal{G}]\bigm\vert \mathcal{H}\bigr] = E[Z\mid\mathcal{H}] E [ E [ Z ∣ G ] H ] = E [ Z ∣ H ] a.s.既知量の取り出し Y Y Y が G \mathcal{G} G -可測で Y Z YZ Y Z が可積分なら E [ Y Z ∣ G ] = Y E [ Z ∣ G ] E[YZ\mid\mathcal{G}] = Y\,E[Z\mid\mathcal{G}] E [ Y Z ∣ G ] = Y E [ Z ∣ G ] a.s.独立性 Z Z Z が G \mathcal{G} G と独立なら E [ Z ∣ G ] = E [ Z ] E[Z\mid\mathcal{G}] = E[Z] E [ Z ∣ G ] = E [ Z ] a.s.イェンセンの不等式 φ \varphi φ が凸で φ ( Z ) \varphi(Z) φ ( Z ) が可積分なら φ ( E [ Z ∣ G ] ) ≤ E [ φ ( Z ) ∣ G ] \varphi(E[Z\mid\mathcal{G}]) \le E[\varphi(Z)\mid\mathcal{G}] φ ( E [ Z ∣ G ]) ≤ E [ φ ( Z ) ∣ G ] a.s.
以下、条件付き期待値に関する等式・不等式はすべて「ほとんど確実に」の意味で述べます。毎回 a.s. と書くのは煩雑なので、断りのないかぎり省略します。
Definition 3.1 (マルチンゲール )
( F n ) n ≥ 0 (\mathcal{F}_n)_{n \ge 0} ( F n ) n ≥ 0 をフィルトレーション、( X n ) n ≥ 0 (X_n)_{n\ge0} ( X n ) n ≥ 0 を実数値確率過程とする。次の三条件が成り立つとき、( X n ) (X_n) ( X n ) は ( F n ) (\mathcal{F}_n) ( F n ) に関するマルチンゲール であるという。
( X n ) (X_n) ( X n ) は ( F n ) (\mathcal{F}_n) ( F n ) に適合している。
すべての n ≥ 0 n \ge 0 n ≥ 0 について E [ ∣ X n ∣ ] < ∞ E[|X_n|] < \infty E [ ∣ X n ∣ ] < ∞ 。
すべての n ≥ 0 n \ge 0 n ≥ 0 について E [ X n + 1 ∣ F n ] = X n E[X_{n+1} \mid \mathcal{F}_n] = X_n E [ X n + 1 ∣ F n ] = X n 。
条件 3 を E [ X n + 1 ∣ F n ] ≥ X n E[X_{n+1}\mid\mathcal{F}_n] \ge X_n E [ X n + 1 ∣ F n ] ≥ X n で置き換えたものを劣マルチンゲール 、E [ X n + 1 ∣ F n ] ≤ X n E[X_{n+1}\mid\mathcal{F}_n] \le X_n E [ X n + 1 ∣ F n ] ≤ X n で置き換えたものを優マルチンゲール という。連続時間 T = [ 0 , ∞ ) T = [0,\infty) T = [ 0 , ∞ ) の場合は、条件 3 を「0 ≤ s ≤ t 0 \le s \le t 0 ≤ s ≤ t に対し E [ X t ∣ F s ] = X s E[X_t \mid \mathcal{F}_s] = X_s E [ X t ∣ F s ] = X s 」に読み替える。
条件 2 を落とすと条件付き期待値そのものが定義できないので、可積分性は飾りではありません。名前の直感は次のとおりです。X n X_n X n を「時刻 n n n での持ち金」と読むと、マルチンゲールは公平な賭け、劣マルチンゲールは自分に有利な賭け、優マルチンゲールは不利な賭けに対応します。日本語の「劣」「優」と有利・不利が逆に見えるのは、劣マルチンゲールが優調和関数ではなく劣調和関数に対応するという解析学側の命名に合わせているためです。
Proposition 3.2 (マルチンゲールの基本性質 )
( X n ) n ≥ 0 (X_n)_{n\ge0} ( X n ) n ≥ 0 を ( F n ) (\mathcal{F}_n) ( F n ) に関するマルチンゲールとする。このとき、0 ≤ n ≤ m 0 \le n \le m 0 ≤ n ≤ m を満たす任意の整数 n , m n, m n , m について
E [ X m ∣ F n ] = X n , E [ X m ] = E [ X 0 ] E[X_m \mid \mathcal{F}_n] = X_n, \qquad E[X_m] = E[X_0] E [ X m ∣ F n ] = X n , E [ X m ] = E [ X 0 ] が成り立つ。
Proof(Proposition 3.2) n n n を固定し、m m m に関する帰納法で第一式を示します。m = n m = n m = n のときは X n X_n X n が F n \mathcal{F}_n F n -可測(Definition 3.1 の条件 1)なので E [ X n ∣ F n ] = X n E[X_n\mid\mathcal{F}_n] = X_n E [ X n ∣ F n ] = X n です。
m ≥ n m \ge n m ≥ n で E [ X m ∣ F n ] = X n E[X_m\mid\mathcal{F}_n] = X_n E [ X m ∣ F n ] = X n が成り立つとします。F n ⊂ F m \mathcal{F}_n \subset \mathcal{F}_m F n ⊂ F m なのでタワー性が使えて、
E [ X m + 1 ∣ F n ] = E [ E [ X m + 1 ∣ F m ] ∣ F n ] = E [ X m ∣ F n ] = X n . E[X_{m+1}\mid\mathcal{F}_n] = E\bigl[E[X_{m+1}\mid\mathcal{F}_m]\bigm|\mathcal{F}_n\bigr] = E[X_m \mid \mathcal{F}_n] = X_n . E [ X m + 1 ∣ F n ] = E [ E [ X m + 1 ∣ F m ] F n ] = E [ X m ∣ F n ] = X n . 二つ目の等号で Definition 3.1 の条件 3 を、三つ目で帰納法の仮定を使いました。
第二式は、第一式で n = 0 n = 0 n = 0 とし、両辺の期待値を取ればよいです。E [ E [ X m ∣ F 0 ] ] = E [ X m ] E\bigl[E[X_m\mid\mathcal{F}_0]\bigr] = E[X_m] E [ E [ X m ∣ F 0 ] ] = E [ X m ] なので E [ X m ] = E [ X 0 ] E[X_m] = E[X_0] E [ X m ] = E [ X 0 ] が出ます。
∎ Proposition 3.2 は「公平なゲームを何回続けても、平均持ち金は初期値のまま」という主張です。ここで重要なのは、m m m が定数 であることです。m m m を確率変数(たとえば「勝ち越したところでやめる」)に置き換えると結論は一般には壊れます。どこまでなら壊れないかが Theorem 4.3 の内容です。
Example 3.3 (ランダムウォークと二次のマルチンゲール )
ξ 1 , ξ 2 , … \xi_1, \xi_2, \ldots ξ 1 , ξ 2 , … を独立同分布な確率変数列で、E [ ξ 1 ] = 0 E[\xi_1] = 0 E [ ξ 1 ] = 0 、σ 2 : = E [ ξ 1 2 ] < ∞ \sigma^2 := E[\xi_1^2] < \infty σ 2 := E [ ξ 1 2 ] < ∞ を満たすものとします。S 0 = x S_0 = x S 0 = x (定数)、S n = x + ξ 1 + ⋯ + ξ n S_n = x + \xi_1 + \cdots + \xi_n S n = x + ξ 1 + ⋯ + ξ n とおき、F n = σ ( ξ 1 , … , ξ n ) \mathcal{F}_n = \sigma(\xi_1,\ldots,\xi_n) F n = σ ( ξ 1 , … , ξ n ) (F 0 \mathcal{F}_0 F 0 は自明な σ \sigma σ -加法族)とします。
まず ( S n ) (S_n) ( S n ) はマルチンゲールです。適合性と可積分性(E ∣ S n ∣ ≤ ∣ x ∣ + n E ∣ ξ 1 ∣ < ∞ E|S_n| \le |x| + nE|\xi_1| < \infty E ∣ S n ∣ ≤ ∣ x ∣ + n E ∣ ξ 1 ∣ < ∞ )は明らかなので、条件 3 を確かめます。ξ n + 1 \xi_{n+1} ξ n + 1 は F n \mathcal{F}_n F n と独立なので
E [ S n + 1 ∣ F n ] = E [ S n ∣ F n ] + E [ ξ n + 1 ∣ F n ] = S n + E [ ξ n + 1 ] = S n + 0 = S n . E[S_{n+1}\mid\mathcal{F}_n] = E[S_n \mid \mathcal{F}_n] + E[\xi_{n+1}\mid\mathcal{F}_n] = S_n + E[\xi_{n+1}] = S_n + 0 = S_n . E [ S n + 1 ∣ F n ] = E [ S n ∣ F n ] + E [ ξ n + 1 ∣ F n ] = S n + E [ ξ n + 1 ] = S n + 0 = S n . 既知量の取り出し(S n S_n S n は F n \mathcal{F}_n F n -可測)と独立性の性質を使いました。
次に M n : = S n 2 − n σ 2 M_n := S_n^2 - n\sigma^2 M n := S n 2 − n σ 2 もマルチンゲールです。S n + 1 = S n + ξ n + 1 S_{n+1} = S_n + \xi_{n+1} S n + 1 = S n + ξ n + 1 を展開すると
E [ S n + 1 2 ∣ F n ] = E [ S n 2 + 2 S n ξ n + 1 + ξ n + 1 2 ∣ F n ] = S n 2 + 2 S n E [ ξ n + 1 ] + E [ ξ n + 1 2 ] = S n 2 + σ 2 E[S_{n+1}^2 \mid \mathcal{F}_n] = E[S_n^2 + 2S_n\xi_{n+1} + \xi_{n+1}^2 \mid \mathcal{F}_n] = S_n^2 + 2S_n E[\xi_{n+1}] + E[\xi_{n+1}^2] = S_n^2 + \sigma^2 E [ S n + 1 2 ∣ F n ] = E [ S n 2 + 2 S n ξ n + 1 + ξ n + 1 2 ∣ F n ] = S n 2 + 2 S n E [ ξ n + 1 ] + E [ ξ n + 1 2 ] = S n 2 + σ 2 となります。ここでも S n S_n S n と S n 2 S_n^2 S n 2 が F n \mathcal{F}_n F n -可測であることと、ξ n + 1 \xi_{n+1} ξ n + 1 が F n \mathcal{F}_n F n と独立であることを使いました。両辺から ( n + 1 ) σ 2 (n+1)\sigma^2 ( n + 1 ) σ 2 を引けば E [ M n + 1 ∣ F n ] = M n E[M_{n+1}\mid\mathcal{F}_n] = M_n E [ M n + 1 ∣ F n ] = M n です。可積分性は E [ S n 2 ] = x 2 + n σ 2 < ∞ E[S_n^2] = x^2 + n\sigma^2 < \infty E [ S n 2 ] = x 2 + n σ 2 < ∞ から従います。
Example 3.4 (ドゥーブのマルチンゲール )
Z Z Z を E ∣ Z ∣ < ∞ E|Z| < \infty E ∣ Z ∣ < ∞ なる確率変数、( F n ) (\mathcal{F}_n) ( F n ) を任意のフィルトレーションとし、X n : = E [ Z ∣ F n ] X_n := E[Z \mid \mathcal{F}_n] X n := E [ Z ∣ F n ] とおきます。条件付き期待値の定義から X n X_n X n は F n \mathcal{F}_n F n -可測、また E ∣ X n ∣ ≤ E [ E [ ∣ Z ∣ ∣ F n ] ] = E ∣ Z ∣ < ∞ E|X_n| \le E\bigl[E[|Z| \mid \mathcal{F}_n]\bigr] = E|Z| < \infty E ∣ X n ∣ ≤ E [ E [ ∣ Z ∣ ∣ F n ] ] = E ∣ Z ∣ < ∞ (イェンセンの不等式を φ ( u ) = ∣ u ∣ \varphi(u)=|u| φ ( u ) = ∣ u ∣ に適用)なので可積分です。タワー性から
E [ X n + 1 ∣ F n ] = E [ E [ Z ∣ F n + 1 ] ∣ F n ] = E [ Z ∣ F n ] = X n E[X_{n+1}\mid\mathcal{F}_n] = E\bigl[E[Z\mid\mathcal{F}_{n+1}]\bigm|\mathcal{F}_n\bigr] = E[Z\mid\mathcal{F}_n] = X_n E [ X n + 1 ∣ F n ] = E [ E [ Z ∣ F n + 1 ] F n ] = E [ Z ∣ F n ] = X n となり、マルチンゲールです。これは「固定された量 Z Z Z についての知識が、情報が増えるにつれて更新されていく」過程であり、ベイズ推定や機械学習の逐次学習の枠組みと同じ形をしています。
Example 3.5 (離散時間の指数マルチンゲール )
Example 3.3 の設定で、さらに φ ( λ ) : = E [ e λ ξ 1 ] < ∞ \varphi(\lambda) := E[e^{\lambda\xi_1}] < \infty φ ( λ ) := E [ e λ ξ 1 ] < ∞ となる λ \lambda λ を取ります(たとえば ξ 1 = ± 1 \xi_1 = \pm1 ξ 1 = ± 1 が等確率なら φ ( λ ) = cosh λ \varphi(\lambda) = \cosh\lambda φ ( λ ) = cosh λ ですべての λ \lambda λ で有限です)。このとき
Y n : = e λ S n φ ( λ ) n Y_n := \frac{e^{\lambda S_n}}{\varphi(\lambda)^n} Y n := φ ( λ ) n e λ S n はマルチンゲールです。実際、e λ S n + 1 = e λ S n e λ ξ n + 1 e^{\lambda S_{n+1}} = e^{\lambda S_n}e^{\lambda\xi_{n+1}} e λ S n + 1 = e λ S n e λ ξ n + 1 と分解し、e λ S n e^{\lambda S_n} e λ S n が F n \mathcal{F}_n F n -可測、e λ ξ n + 1 e^{\lambda\xi_{n+1}} e λ ξ n + 1 が F n \mathcal{F}_n F n と独立であることを使うと
E [ e λ S n + 1 ∣ F n ] = e λ S n E [ e λ ξ n + 1 ] = e λ S n φ ( λ ) E[e^{\lambda S_{n+1}}\mid\mathcal{F}_n] = e^{\lambda S_n}E[e^{\lambda\xi_{n+1}}] = e^{\lambda S_n}\varphi(\lambda) E [ e λ S n + 1 ∣ F n ] = e λ S n E [ e λ ξ n + 1 ] = e λ S n φ ( λ ) となり、両辺を φ ( λ ) n + 1 \varphi(\lambda)^{n+1} φ ( λ ) n + 1 で割れば主張を得ます。λ \lambda λ を自由に動かせるので、この族は到達時刻のラプラス変換を計算する道具になります(Exercise 7.4 )。
「どんな賭け方をしても公平なゲームは公平なまま」を定式化します。鍵は、時刻 n n n の賭け金は時刻 n − 1 n-1 n − 1 までの情報だけで決まる、という制約です。
Definition 3.6 (可予測過程とマルチンゲール変換 )
確率過程 ( H n ) n ≥ 1 (H_n)_{n\ge1} ( H n ) n ≥ 1 が、すべての n ≥ 1 n \ge 1 n ≥ 1 について H n H_n H n が F n − 1 \mathcal{F}_{n-1} F n − 1 -可測であるとき、( H n ) (H_n) ( H n ) は可予測 であるという。( X n ) n ≥ 0 (X_n)_{n\ge0} ( X n ) n ≥ 0 に対し
( H ⋅ X ) 0 : = 0 , ( H ⋅ X ) n : = ∑ k = 1 n H k ( X k − X k − 1 ) ( n ≥ 1 ) (H\cdot X)_0 := 0, \qquad (H\cdot X)_n := \sum_{k=1}^{n} H_k\,(X_k - X_{k-1}) \quad (n \ge 1) ( H ⋅ X ) 0 := 0 , ( H ⋅ X ) n := k = 1 ∑ n H k ( X k − X k − 1 ) ( n ≥ 1 ) で定まる過程 ( H ⋅ X ) n (H\cdot X)_n ( H ⋅ X ) n を、H H H による X X X のマルチンゲール変換 という。
H k H_k H k は「k k k 回目の勝負に賭ける単位数」、X k − X k − 1 X_k - X_{k-1} X k − X k − 1 は「1 単位あたりの損益」、( H ⋅ X ) n (H\cdot X)_n ( H ⋅ X ) n は「n n n 回目までの累積損益」です。倍賭け戦略も、直前までの勝敗を見て賭け金を決めるので可予測です。
Proof(Theorem 3.7) 適合性を確かめます。H k H_k H k は F k − 1 ⊂ F n \mathcal{F}_{k-1} \subset \mathcal{F}_n F k − 1 ⊂ F n -可測(k ≤ n k \le n k ≤ n )、X k , X k − 1 X_k, X_{k-1} X k , X k − 1 は F n \mathcal{F}_n F n -可測なので、有限和である ( H ⋅ X ) n (H\cdot X)_n ( H ⋅ X ) n は F n \mathcal{F}_n F n -可測です。
可積分性を確かめます。三角不等式から
E [ ∣ ( H ⋅ X ) n ∣ ] ≤ ∑ k = 1 n C k ( E ∣ X k ∣ + E ∣ X k − 1 ∣ ) < ∞ E\bigl[|(H\cdot X)_n|\bigr] \le \sum_{k=1}^{n} C_k\bigl(E|X_k| + E|X_{k-1}|\bigr) < \infty E [ ∣ ( H ⋅ X ) n ∣ ] ≤ k = 1 ∑ n C k ( E ∣ X k ∣ + E ∣ X k − 1 ∣ ) < ∞ であり、右辺の有限性は Definition 3.1 の条件 2 によります。
最後に条件 3 です。( H ⋅ X ) n + 1 − ( H ⋅ X ) n = H n + 1 ( X n + 1 − X n ) (H\cdot X)_{n+1} - (H\cdot X)_n = H_{n+1}(X_{n+1}-X_n) ( H ⋅ X ) n + 1 − ( H ⋅ X ) n = H n + 1 ( X n + 1 − X n ) であり、H n + 1 H_{n+1} H n + 1 は F n \mathcal{F}_n F n -可測で有界ですから、既知量の取り出しが使えて
E [ H n + 1 ( X n + 1 − X n ) ∣ F n ] = H n + 1 E [ X n + 1 − X n ∣ F n ] = H n + 1 ⋅ 0 = 0. E\bigl[H_{n+1}(X_{n+1}-X_n)\bigm|\mathcal{F}_n\bigr] = H_{n+1}\,E[X_{n+1}-X_n\mid\mathcal{F}_n] = H_{n+1}\cdot 0 = 0 . E [ H n + 1 ( X n + 1 − X n ) F n ] = H n + 1 E [ X n + 1 − X n ∣ F n ] = H n + 1 ⋅ 0 = 0. 最後の等号は Definition 3.1 の条件 3 です。よって E [ ( H ⋅ X ) n + 1 ∣ F n ] = ( H ⋅ X ) n E[(H\cdot X)_{n+1}\mid\mathcal{F}_n] = (H\cdot X)_n E [( H ⋅ X ) n + 1 ∣ F n ] = ( H ⋅ X ) n となります。( H ⋅ X ) 0 = 0 (H\cdot X)_0 = 0 ( H ⋅ X ) 0 = 0 と Proposition 3.2 から E [ ( H ⋅ X ) n ] = 0 E[(H\cdot X)_n] = 0 E [( H ⋅ X ) n ] = 0 です。
∎ Definition 4.1 (停止時刻 )
写像 τ : Ω → { 0 , 1 , 2 , … } ∪ { + ∞ } \tau : \Omega \to \{0,1,2,\ldots\} \cup \{+\infty\} τ : Ω → { 0 , 1 , 2 , … } ∪ { + ∞ } が、すべての n ≥ 0 n \ge 0 n ≥ 0 について
{ τ ≤ n } ∈ F n \{\tau \le n\} \in \mathcal{F}_n { τ ≤ n } ∈ F n を満たすとき、τ \tau τ を ( F n ) (\mathcal{F}_n) ( F n ) に関する停止時刻 という。連続時間の場合は、すべての t ≥ 0 t \ge 0 t ≥ 0 について { τ ≤ t } ∈ F t \{\tau \le t\} \in \mathcal{F}_t { τ ≤ t } ∈ F t を要請する。
離散時間では { τ ≤ n } ∈ F n \{\tau \le n\} \in \mathcal{F}_n { τ ≤ n } ∈ F n と { τ = n } ∈ F n \{\tau = n\}\in\mathcal{F}_n { τ = n } ∈ F n は同値です。実際、{ τ = n } = { τ ≤ n } ∖ { τ ≤ n − 1 } \{\tau = n\} = \{\tau\le n\}\setminus\{\tau\le n-1\} { τ = n } = { τ ≤ n } ∖ { τ ≤ n − 1 } で、{ τ ≤ n − 1 } ∈ F n − 1 ⊂ F n \{\tau\le n-1\}\in\mathcal{F}_{n-1}\subset\mathcal{F}_n { τ ≤ n − 1 } ∈ F n − 1 ⊂ F n ですから前者から後者が出ます。逆は { τ ≤ n } = ⋃ k ≤ n { τ = k } \{\tau\le n\} = \bigcup_{k\le n}\{\tau = k\} { τ ≤ n } = ⋃ k ≤ n { τ = k } です。
定義の意味は「時刻 n n n の時点で、もう止まったかどうかが判定できる」ことです。「株価が最高値を付けた日に売る」は最高値がいつだったか将来にならないと分からないので停止時刻ではありません。「株価が初めて 1000 円を超えた日に売る」は停止時刻です。
Lemma 4.2 (止めた過程もマルチンゲール )
( X n ) n ≥ 0 (X_n)_{n\ge0} ( X n ) n ≥ 0 を ( F n ) (\mathcal{F}_n) ( F n ) に関するマルチンゲール、τ \tau τ を ( F n ) (\mathcal{F}_n) ( F n ) に関する停止時刻とする。X n τ : = X n ∧ τ X^\tau_n := X_{n\wedge\tau} X n τ := X n ∧ τ (n ∧ τ = min ( n , τ ) n\wedge\tau = \min(n,\tau) n ∧ τ = min ( n , τ ) )とおくと、( X n τ ) n ≥ 0 (X^\tau_n)_{n\ge0} ( X n τ ) n ≥ 0 もマルチンゲールであり、特に E [ X n ∧ τ ] = E [ X 0 ] E[X_{n\wedge\tau}] = E[X_0] E [ X n ∧ τ ] = E [ X 0 ] がすべての n n n で成り立つ。
Proof(Lemma 4.2) H k : = 1 { τ ≥ k } H_k := \mathbf{1}_{\{\tau \ge k\}} H k := 1 { τ ≥ k } とおきます。{ τ ≥ k } = { τ ≤ k − 1 } c \{\tau \ge k\} = \{\tau \le k-1\}^{c} { τ ≥ k } = { τ ≤ k − 1 } c であり、Definition 4.1 より { τ ≤ k − 1 } ∈ F k − 1 \{\tau\le k-1\}\in\mathcal{F}_{k-1} { τ ≤ k − 1 } ∈ F k − 1 ですから H k H_k H k は F k − 1 \mathcal{F}_{k-1} F k − 1 -可測、すなわち可予測です。また 0 ≤ H k ≤ 1 0 \le H_k \le 1 0 ≤ H k ≤ 1 で有界です。
このとき
( H ⋅ X ) n = ∑ k = 1 n 1 { τ ≥ k } ( X k − X k − 1 ) = ∑ k = 1 n ∧ τ ( X k − X k − 1 ) = X n ∧ τ − X 0 (H\cdot X)_n = \sum_{k=1}^{n}\mathbf{1}_{\{\tau\ge k\}}(X_k - X_{k-1}) = \sum_{k=1}^{n\wedge\tau}(X_k - X_{k-1}) = X_{n\wedge\tau} - X_0 ( H ⋅ X ) n = k = 1 ∑ n 1 { τ ≥ k } ( X k − X k − 1 ) = k = 1 ∑ n ∧ τ ( X k − X k − 1 ) = X n ∧ τ − X 0 が成り立ちます。二つ目の等号は、k ≤ τ k \le \tau k ≤ τ のときだけ項が残ることと、k k k は n n n 以下を動くことによります。三つ目は望遠鏡和です。
Theorem 3.7 より ( H ⋅ X ) n (H\cdot X)_n ( H ⋅ X ) n はマルチンゲールで、定数 X 0 X_0 X 0 を足しても(X 0 X_0 X 0 は F 0 \mathcal{F}_0 F 0 -可測で可積分なので)マルチンゲール性は保たれます。実際 E [ X 0 ∣ F n ] = X 0 E[X_0\mid\mathcal{F}_n] = X_0 E [ X 0 ∣ F n ] = X 0 です。よって X n τ = X 0 + ( H ⋅ X ) n X^\tau_n = X_0 + (H\cdot X)_n X n τ = X 0 + ( H ⋅ X ) n はマルチンゲールであり、E [ X n τ ] = E [ X 0 τ ] = E [ X 0 ] E[X^\tau_n] = E[X^\tau_0] = E[X_0] E [ X n τ ] = E [ X 0 τ ] = E [ X 0 ] を得ます。
∎ Lemma 4.2 は E [ X n ∧ τ ] = E [ X 0 ] E[X_{n\wedge\tau}] = E[X_0] E [ X n ∧ τ ] = E [ X 0 ] を与えます。欲しいのは E [ X τ ] = E [ X 0 ] E[X_\tau] = E[X_0] E [ X τ ] = E [ X 0 ] ですから、n → ∞ n \to \infty n → ∞ の極限を取る必要があります。極限と期待値の交換に条件が要る、というのがこの定理の全体像です。
Theorem 4.3 (ドゥーブの任意停止定理 )
( X n ) n ≥ 0 (X_n)_{n\ge0} ( X n ) n ≥ 0 を ( F n ) (\mathcal{F}_n) ( F n ) に関するマルチンゲール、τ \tau τ を停止時刻とする。次のいずれかが成り立つならば、X τ X_\tau X τ は可積分で E [ X τ ] = E [ X 0 ] E[X_\tau] = E[X_0] E [ X τ ] = E [ X 0 ] である。
(a) ある定数 K K K が存在して τ ≤ K \tau \le K τ ≤ K が確率 1 で成り立つ。
(b) P ( τ < ∞ ) = 1 P(\tau < \infty) = 1 P ( τ < ∞ ) = 1 であり、かつ可積分な確率変数 Y Y Y (E [ Y ] < ∞ E[Y] < \infty E [ Y ] < ∞ 、Y ≥ 0 Y \ge 0 Y ≥ 0 )が存在して、すべての n n n について ∣ X n ∧ τ ∣ ≤ Y |X_{n\wedge\tau}| \le Y ∣ X n ∧ τ ∣ ≤ Y が確率 1 で成り立つ。
Proof(Theorem 4.3) (a) の場合。n = K n = K n = K とすると K ∧ τ = τ K \wedge \tau = \tau K ∧ τ = τ なので、Lemma 4.2 の等式 E [ X K ∧ τ ] = E [ X 0 ] E[X_{K\wedge\tau}] = E[X_0] E [ X K ∧ τ ] = E [ X 0 ] がそのまま E [ X τ ] = E [ X 0 ] E[X_\tau] = E[X_0] E [ X τ ] = E [ X 0 ] です。可積分性も X τ = X K ∧ τ X_\tau = X_{K\wedge\tau} X τ = X K ∧ τ が可積分であることから従います。
(b) の場合。P ( τ < ∞ ) = 1 P(\tau<\infty) = 1 P ( τ < ∞ ) = 1 なので、確率 1 の事象の上で、n n n が十分大きければ n ∧ τ = τ n \wedge \tau = \tau n ∧ τ = τ となります。したがって
X n ∧ τ → n → ∞ X τ a.s. X_{n\wedge\tau} \xrightarrow{n\to\infty} X_\tau \quad \text{a.s.} X n ∧ τ n → ∞ X τ a.s. です。仮定より ∣ X n ∧ τ ∣ ≤ Y |X_{n\wedge\tau}| \le Y ∣ X n ∧ τ ∣ ≤ Y で Y Y Y は可積分ですから、優収束定理(ルベーグ積分の定義と収束定理 の Theorem 7.3[ルベーグ積分の定義と収束定理] を参照)が適用できて、X τ X_\tau X τ は可積分であり
E [ X τ ] = lim n → ∞ E [ X n ∧ τ ] = E [ X 0 ] E[X_\tau] = \lim_{n\to\infty} E[X_{n\wedge\tau}] = E[X_0] E [ X τ ] = n → ∞ lim E [ X n ∧ τ ] = E [ X 0 ] となります。最後の等号は Lemma 4.2 です。
∎ 実際に使いやすいのは次の形です。
Corollary 4.4 (増分が有界な場合 )
( X n ) (X_n) ( X n ) をマルチンゲール、τ \tau τ を停止時刻とする。E [ τ ] < ∞ E[\tau] < \infty E [ τ ] < ∞ であり、さらにある定数 C C C が存在してすべての n ≥ 1 n \ge 1 n ≥ 1 で ∣ X n − X n − 1 ∣ ≤ C |X_n - X_{n-1}| \le C ∣ X n − X n − 1 ∣ ≤ C が確率 1 で成り立つならば、E [ X τ ] = E [ X 0 ] E[X_\tau] = E[X_0] E [ X τ ] = E [ X 0 ] である。
Proof(Corollary 4.4) E [ τ ] < ∞ E[\tau]<\infty E [ τ ] < ∞ より P ( τ < ∞ ) = 1 P(\tau<\infty) = 1 P ( τ < ∞ ) = 1 です。また
∣ X n ∧ τ − X 0 ∣ = ∣ ∑ k = 1 n ∧ τ ( X k − X k − 1 ) ∣ ≤ ∑ k = 1 n ∧ τ ∣ X k − X k − 1 ∣ ≤ C ( n ∧ τ ) ≤ C τ |X_{n\wedge\tau} - X_0| = \Bigl|\sum_{k=1}^{n\wedge\tau}(X_k - X_{k-1})\Bigr| \le \sum_{k=1}^{n\wedge\tau} |X_k - X_{k-1}| \le C\,(n\wedge\tau) \le C\tau ∣ X n ∧ τ − X 0 ∣ = k = 1 ∑ n ∧ τ ( X k − X k − 1 ) ≤ k = 1 ∑ n ∧ τ ∣ X k − X k − 1 ∣ ≤ C ( n ∧ τ ) ≤ C τ なので、Y : = ∣ X 0 ∣ + C τ Y := |X_0| + C\tau Y := ∣ X 0 ∣ + C τ とおけば ∣ X n ∧ τ ∣ ≤ Y |X_{n\wedge\tau}| \le Y ∣ X n ∧ τ ∣ ≤ Y がすべての n n n で成り立ちます。E [ Y ] ≤ E ∣ X 0 ∣ + C E [ τ ] < ∞ E[Y] \le E|X_0| + CE[\tau] < \infty E [ Y ] ≤ E ∣ X 0 ∣ + C E [ τ ] < ∞ ですから、Theorem 4.3 の (b) が適用できます。
∎ ギャンブラーの破産問題。所持金 x から出発し、0 または N に到達した時刻が停止時刻 τ Example 4.5 (破産確率と平均ゲーム時間 )
ξ 1 , ξ 2 , … \xi_1,\xi_2,\ldots ξ 1 , ξ 2 , … を P ( ξ i = 1 ) = P ( ξ i = − 1 ) = 1 / 2 P(\xi_i = 1) = P(\xi_i = -1) = 1/2 P ( ξ i = 1 ) = P ( ξ i = − 1 ) = 1/2 の独立同分布列とし、整数 N ≥ 2 N \ge 2 N ≥ 2 と 0 < x < N 0 < x < N 0 < x < N を固定して S n = x + ξ 1 + ⋯ + ξ n S_n = x + \xi_1+\cdots+\xi_n S n = x + ξ 1 + ⋯ + ξ n とします。F n = σ ( ξ 1 , … , ξ n ) \mathcal{F}_n = \sigma(\xi_1,\ldots,\xi_n) F n = σ ( ξ 1 , … , ξ n ) とし、
τ : = inf { n ≥ 0 : S n = 0 または S n = N } \tau := \inf\{n \ge 0 : S_n = 0 \ \text{または}\ S_n = N\} τ := inf { n ≥ 0 : S n = 0 または S n = N } とおきます(そのような n n n がなければ τ = + ∞ \tau = +\infty τ = + ∞ )。{ τ ≤ n } \{\tau \le n\} { τ ≤ n } は S 0 , … , S n S_0,\ldots,S_n S 0 , … , S n だけで決まるので τ \tau τ は停止時刻です。
第一段:E [ τ ] < ∞ E[\tau] < \infty E [ τ ] < ∞ を示す。 N N N 回連続で + 1 +1 + 1 が出れば、出発点が ( 0 , N ) (0,N) ( 0 , N ) のどこであっても値は N N N 以上に到達するはずですが、歩幅が 1 なので途中で必ず N N N に一致します。つまりその区間内で吸収されます。各長さ N N N のブロックで + 1 +1 + 1 が N N N 回連続する確率は 2 − N 2^{-N} 2 − N で、ブロックごとに独立ですから
P ( τ > ( m + 1 ) N ) ≤ ( 1 − 2 − N ) P ( τ > m N ) , P(\tau > (m+1)N) \le \bigl(1 - 2^{-N}\bigr)P(\tau > mN), P ( τ > ( m + 1 ) N ) ≤ ( 1 − 2 − N ) P ( τ > m N ) , したがって P ( τ > m N ) ≤ ( 1 − 2 − N ) m P(\tau > mN) \le (1-2^{-N})^m P ( τ > m N ) ≤ ( 1 − 2 − N ) m です。q : = 1 − 2 − N < 1 q := 1-2^{-N} < 1 q := 1 − 2 − N < 1 とおくと
E [ τ ] = ∑ k ≥ 0 P ( τ > k ) ≤ ∑ m ≥ 0 N P ( τ > m N ) ≤ N ∑ m ≥ 0 q m = N 1 − q = N 2 N < ∞ E[\tau] = \sum_{k \ge 0} P(\tau > k) \le \sum_{m\ge0} N\,P(\tau > mN) \le N\sum_{m\ge0}q^m = \frac{N}{1-q} = N\,2^{N} < \infty E [ τ ] = k ≥ 0 ∑ P ( τ > k ) ≤ m ≥ 0 ∑ N P ( τ > m N ) ≤ N m ≥ 0 ∑ q m = 1 − q N = N 2 N < ∞ となり、特に P ( τ < ∞ ) = 1 P(\tau<\infty)=1 P ( τ < ∞ ) = 1 です。
第二段:破産確率。 ( S n ) (S_n) ( S n ) は Example 3.3 よりマルチンゲールです。n ≤ τ n \le \tau n ≤ τ の間 0 ≤ S n ≤ N 0 \le S_n \le N 0 ≤ S n ≤ N なので ∣ S n ∧ τ ∣ ≤ N |S_{n\wedge\tau}| \le N ∣ S n ∧ τ ∣ ≤ N であり、定数 Y = N Y = N Y = N が Theorem 4.3 の (b) の条件を満たします。よって
x = E [ S 0 ] = E [ S τ ] = 0 ⋅ P ( S τ = 0 ) + N ⋅ P ( S τ = N ) x = E[S_0] = E[S_\tau] = 0\cdot P(S_\tau = 0) + N\cdot P(S_\tau = N) x = E [ S 0 ] = E [ S τ ] = 0 ⋅ P ( S τ = 0 ) + N ⋅ P ( S τ = N ) となり、目標額 N N N に到達する確率は
P ( S τ = N ) = x N , P ( S τ = 0 ) = 1 − x N P(S_\tau = N) = \frac{x}{N}, \qquad P(S_\tau = 0) = 1 - \frac{x}{N} P ( S τ = N ) = N x , P ( S τ = 0 ) = 1 − N x と決まります。所持金 1 万円で 100 万円を目指すなら、成功確率はちょうど 1 / 100 1/100 1/100 です。
第三段:平均ゲーム時間。 σ 2 = E [ ξ 1 2 ] = 1 \sigma^2 = E[\xi_1^2] = 1 σ 2 = E [ ξ 1 2 ] = 1 なので、Example 3.3 より M n : = S n 2 − n M_n := S_n^2 - n M n := S n 2 − n はマルチンゲールです。∣ M n ∧ τ ∣ ≤ N 2 + ( n ∧ τ ) ≤ N 2 + τ |M_{n\wedge\tau}| \le N^2 + (n\wedge\tau) \le N^2 + \tau ∣ M n ∧ τ ∣ ≤ N 2 + ( n ∧ τ ) ≤ N 2 + τ であり、第一段より E [ N 2 + τ ] < ∞ E[N^2+\tau]<\infty E [ N 2 + τ ] < ∞ ですから、再び Theorem 4.3 の (b) が使えて E [ M τ ] = E [ M 0 ] = x 2 E[M_\tau] = E[M_0] = x^2 E [ M τ ] = E [ M 0 ] = x 2 です。すなわち
x 2 = E [ S τ 2 ] − E [ τ ] = N 2 ⋅ x N − E [ τ ] = N x − E [ τ ] x^2 = E[S_\tau^2] - E[\tau] = N^2\cdot\frac{x}{N} - E[\tau] = Nx - E[\tau] x 2 = E [ S τ 2 ] − E [ τ ] = N 2 ⋅ N x − E [ τ ] = N x − E [ τ ] となり、
E [ τ ] = N x − x 2 = x ( N − x ) E[\tau] = Nx - x^2 = x(N-x) E [ τ ] = N x − x 2 = x ( N − x ) を得ます。x = N / 2 x = N/2 x = N /2 のとき最大で N 2 / 4 N^2/4 N 2 /4 回。N = 100 N = 100 N = 100 、x = 50 x = 50 x = 50 なら平均 2500 回の勝負が必要という、直感よりかなり大きい数字です。
ランダムウォーク S n S_n S n の分散は n n n に比例します。時間を 1 / n 1/n 1/ n 倍に縮めると同時に空間を 1 / n 1/\sqrt{n} 1/ n 倍に縮めると分散が保存されるので、S ⌊ n t ⌋ / n S_{\lfloor nt\rfloor}/\sqrt{n} S ⌊ n t ⌋ / n という縮尺のもとで極限を取るのが自然です。中心極限定理(Theorem 5.6)[The Law of Large Numbers and the Central Limit Theorem] から、固定した t t t でこの量は N ( 0 , t ) N(0,t) N ( 0 , t ) に分布収束します。その極限を「過程として」定義したものが次です。
Definition 5.1 (ブラウン運動(ウィーナー過程) )
確率空間 ( Ω , F , P ) (\Omega,\mathcal{F},P) ( Ω , F , P ) 上の実数値確率過程 ( B t ) t ≥ 0 (B_t)_{t\ge0} ( B t ) t ≥ 0 が次の四条件を満たすとき、( B t ) (B_t) ( B t ) を標準ブラウン運動 (ウィーナー過程)という。
P ( B 0 = 0 ) = 1 P(B_0 = 0) = 1 P ( B 0 = 0 ) = 1 。
(独立増分)任意の 0 ≤ t 0 < t 1 < ⋯ < t k 0 \le t_0 < t_1 < \cdots < t_k 0 ≤ t 0 < t 1 < ⋯ < t k に対し、増分 B t 1 − B t 0 , B t 2 − B t 1 , … , B t k − B t k − 1 B_{t_1}-B_{t_0},\ B_{t_2}-B_{t_1},\ \ldots,\ B_{t_k}-B_{t_{k-1}} B t 1 − B t 0 , B t 2 − B t 1 , … , B t k − B t k − 1 は独立である。
(正規増分)任意の 0 ≤ s < t 0 \le s < t 0 ≤ s < t に対し、B t − B s B_t - B_s B t − B s は平均 0 0 0 、分散 t − s t-s t − s の正規分布 N ( 0 , t − s ) N(0,t-s) N ( 0 , t − s ) に従う。
(連続性)確率 1 で、写像 t ↦ B t ( ω ) t \mapsto B_t(\omega) t ↦ B t ( ω ) は [ 0 , ∞ ) [0,\infty) [ 0 , ∞ ) 上連続である。
条件 3 で分散が t − s t-s t − s であること、すなわち s s s と t t t の差だけ で決まることが、増分の定常性を含んでいます。条件 4 は他の三条件から導かれるものではありません。同じ有限次元分布を持ちながら道が至るところ不連続になるような過程も作れてしまうので、連続性は独立に要請する必要があります。
Definition 5.4 (中心化ガウス過程 )
確率過程 ( X t ) t ≥ 0 (X_t)_{t\ge0} ( X t ) t ≥ 0 が、任意の有限個の時刻 t 1 , … , t k ≥ 0 t_1,\ldots,t_k \ge 0 t 1 , … , t k ≥ 0 と任意の実数 a 1 , … , a k a_1,\ldots,a_k a 1 , … , a k に対し ∑ i = 1 k a i X t i \sum_{i=1}^{k} a_i X_{t_i} ∑ i = 1 k a i X t i が正規分布(分散 0 0 0 の退化した場合も含む)に従うとき、( X t ) (X_t) ( X t ) をガウス過程 という。さらにすべての t t t で E [ X t ] = 0 E[X_t]=0 E [ X t ] = 0 のとき中心化ガウス過程 という。関数 ( s , t ) ↦ E [ X s X t ] (s,t)\mapsto E[X_sX_t] ( s , t ) ↦ E [ X s X t ] をその共分散関数 という。
Proposition 5.5 (ブラウン運動のガウス過程による特徴づけ )
連続な見本道を持ち P ( X 0 = 0 ) = 1 P(X_0=0)=1 P ( X 0 = 0 ) = 1 を満たす確率過程 ( X t ) t ≥ 0 (X_t)_{t\ge0} ( X t ) t ≥ 0 について、次は同値である。
(i) ( X t ) (X_t) ( X t ) は標準ブラウン運動である。
(ii) ( X t ) (X_t) ( X t ) は中心化ガウス過程であり、共分散関数が E [ X s X t ] = min ( s , t ) E[X_sX_t] = \min(s,t) E [ X s X t ] = min ( s , t ) である。
Proof(Proposition 5.5) (i) ならば (ii)。 時刻 0 = t 0 < t 1 < ⋯ < t k 0 = t_0 < t_1 < \cdots < t_k 0 = t 0 < t 1 < ⋯ < t k を取り、増分ベクトルを Δ = ( X t 1 − X t 0 , … , X t k − X t k − 1 ) T \boldsymbol{\Delta} = (X_{t_1}-X_{t_0}, \ldots, X_{t_k}-X_{t_{k-1}})^{\mathsf{T}} Δ = ( X t 1 − X t 0 , … , X t k − X t k − 1 ) T とします。Definition 5.1 の条件 2, 3 より Δ \boldsymbol{\Delta} Δ の成分は独立な正規確率変数なので、Δ \boldsymbol{\Delta} Δ は多変量正規分布に従います。X t i = ∑ j ≤ i ( X t j − X t j − 1 ) X_{t_i} = \sum_{j\le i}(X_{t_j}-X_{t_{j-1}}) X t i = ∑ j ≤ i ( X t j − X t j − 1 ) ですから ( X t 1 , … , X t k ) T = A Δ (X_{t_1},\ldots,X_{t_k})^{\mathsf{T}} = A\boldsymbol{\Delta} ( X t 1 , … , X t k ) T = A Δ (A A A は下三角の 0 0 0 -1 1 1 行列)と書けます。多変量正規分布の線形像は多変量正規分布ですから、( X t 1 , … , X t k ) (X_{t_1},\ldots,X_{t_k}) ( X t 1 , … , X t k ) は多変量正規、すなわち ( X t ) (X_t) ( X t ) はガウス過程です。平均は E [ X t ] = E [ X t − X 0 ] = 0 E[X_t] = E[X_t - X_0] = 0 E [ X t ] = E [ X t − X 0 ] = 0 です。
共分散は s < t s < t s < t として
E [ X s X t ] = E [ X s ( X t − X s ) ] + E [ X s 2 ] = E [ X s ] E [ X t − X s ] + s = 0 + s = s E[X_sX_t] = E\bigl[X_s(X_t - X_s)\bigr] + E[X_s^2] = E[X_s]\,E[X_t-X_s] + s = 0 + s = s E [ X s X t ] = E [ X s ( X t − X s ) ] + E [ X s 2 ] = E [ X s ] E [ X t − X s ] + s = 0 + s = s と計算されます。第二の等号で、X s = X s − X 0 X_s = X_s - X_0 X s = X s − X 0 と X t − X s X_t - X_s X t − X s が独立(条件 2)であることと、Var ( X s ) = s \operatorname{Var}(X_s) = s Var ( X s ) = s (条件 3)を使いました。s = t s = t s = t の場合も E [ X t 2 ] = t E[X_t^2] = t E [ X t 2 ] = t なので、まとめて E [ X s X t ] = min ( s , t ) E[X_sX_t] = \min(s,t) E [ X s X t ] = min ( s , t ) です。
(ii) ならば (i)。 条件 1 は仮定、条件 4 も仮定です。0 ≤ t 0 < t 1 < ⋯ < t k 0 \le t_0 < t_1 < \cdots < t_k 0 ≤ t 0 < t 1 < ⋯ < t k を取ります。各増分 X t i − X t i − 1 X_{t_i}-X_{t_{i-1}} X t i − X t i − 1 は ( X t 0 , … , X t k ) (X_{t_0},\ldots,X_{t_k}) ( X t 0 , … , X t k ) の線形結合ですから、増分ベクトル全体は多変量正規分布に従います。s < t ≤ u < v s < t \le u < v s < t ≤ u < v のとき
E [ ( X t − X s ) ( X v − X u ) ] = min ( t , v ) − min ( t , u ) − min ( s , v ) + min ( s , u ) = t − t − s + s = 0 E\bigl[(X_t-X_s)(X_v-X_u)\bigr] = \min(t,v)-\min(t,u)-\min(s,v)+\min(s,u) = t - t - s + s = 0 E [ ( X t − X s ) ( X v − X u ) ] = min ( t , v ) − min ( t , u ) − min ( s , v ) + min ( s , u ) = t − t − s + s = 0 です(s < t ≤ u < v s < t \le u < v s < t ≤ u < v より min ( t , v ) = t \min(t,v)=t min ( t , v ) = t 、min ( t , u ) = t \min(t,u)=t min ( t , u ) = t 、min ( s , v ) = s \min(s,v)=s min ( s , v ) = s 、min ( s , u ) = s \min(s,u)=s min ( s , u ) = s )。よって相異なる増分は無相関で、多変量正規分布において無相関は共分散行列が対角であることを意味し、そのとき同時密度が積に分解するので独立です。これが条件 2 です。
最後に条件 3 です。E [ X t − X s ] = 0 E[X_t-X_s]=0 E [ X t − X s ] = 0 であり、
E [ ( X t − X s ) 2 ] = E [ X t 2 ] − 2 E [ X s X t ] + E [ X s 2 ] = t − 2 s + s = t − s E\bigl[(X_t-X_s)^2\bigr] = E[X_t^2] - 2E[X_sX_t] + E[X_s^2] = t - 2s + s = t-s E [ ( X t − X s ) 2 ] = E [ X t 2 ] − 2 E [ X s X t ] + E [ X s 2 ] = t − 2 s + s = t − s なので、X t − X s X_t - X_s X t − X s は N ( 0 , t − s ) N(0,t-s) N ( 0 , t − s ) に従います。
∎ Proposition 5.5 の利点は、共分散関数さえ計算すればブラウン運動かどうかが判定できることです。実際に使ってみます。
Proposition 5.6 (ブラウン運動の不変性 )
( B t ) t ≥ 0 (B_t)_{t\ge0} ( B t ) t ≥ 0 を標準ブラウン運動とする。次の各過程もまた標準ブラウン運動である。
(a) (対称性)X t : = − B t X_t := -B_t X t := − B t 。
(b) (スケーリング)c > 0 c > 0 c > 0 を定数として X t : = c − 1 / 2 B c t X_t := c^{-1/2}B_{ct} X t := c − 1/2 B c t 。
(c) (マルコフ性・平行移動)s ≥ 0 s \ge 0 s ≥ 0 を定数として X t : = B t + s − B s X_t := B_{t+s}-B_s X t := B t + s − B s 。
(d) (時間反転)X 0 : = 0 X_0 := 0 X 0 := 0 、t > 0 t > 0 t > 0 に対して X t : = t B 1 / t X_t := t\,B_{1/t} X t := t B 1/ t 。
Proof(Proposition 5.6) いずれも Proposition 5.5 を使います。(a)(b)(c) では道の連続性は明らかです(連続関数の定数倍・合成・平行移動だから)。
(a) − B -B − B の有限次元分布は B B B のそれの線形像なのでガウス過程であり、E [ ( − B s ) ( − B t ) ] = E [ B s B t ] = min ( s , t ) E[(-B_s)(-B_t)] = E[B_sB_t] = \min(s,t) E [( − B s ) ( − B t )] = E [ B s B t ] = min ( s , t ) です。
(b) X t 1 , … , X t k X_{t_1},\ldots,X_{t_k} X t 1 , … , X t k は B c t 1 , … , B c t k B_{ct_1},\ldots,B_{ct_k} B c t 1 , … , B c t k の定数倍なのでガウス過程です。平均は 0 0 0 、共分散は
E [ X s X t ] = c − 1 E [ B c s B c t ] = c − 1 min ( c s , c t ) = min ( s , t ) E[X_sX_t] = c^{-1}E[B_{cs}B_{ct}] = c^{-1}\min(cs,ct) = \min(s,t) E [ X s X t ] = c − 1 E [ B cs B c t ] = c − 1 min ( cs , c t ) = min ( s , t ) です。
(c) X t 1 , … , X t k X_{t_1},\ldots,X_{t_k} X t 1 , … , X t k は B B B の値の線形結合なのでガウス過程、平均は 0 0 0 です。u ≤ v u \le v u ≤ v とすると
E [ X u X v ] = E [ ( B u + s − B s ) ( B v + s − B s ) ] = min ( u + s , v + s ) − s − s + s = ( u + s ) − s = u E[X_uX_v] = E\bigl[(B_{u+s}-B_s)(B_{v+s}-B_s)\bigr] = \min(u+s,v+s) - s - s + s = (u+s) - s = u E [ X u X v ] = E [ ( B u + s − B s ) ( B v + s − B s ) ] = min ( u + s , v + s ) − s − s + s = ( u + s ) − s = u となり min ( u , v ) \min(u,v) min ( u , v ) に一致します。
(d) t > 0 t > 0 t > 0 では t ↦ t B 1 / t t\mapsto tB_{1/t} t ↦ t B 1/ t は連続関数の積なので連続です。ガウス過程であることは B B B の値の定数倍だからで、平均は 0 0 0 、0 < s ≤ t 0 < s \le t 0 < s ≤ t に対し
E [ X s X t ] = s t E [ B 1 / s B 1 / t ] = s t ⋅ min ( 1 / s , 1 / t ) = s t ⋅ 1 t = s = min ( s , t ) E[X_sX_t] = st\,E[B_{1/s}B_{1/t}] = st\cdot\min(1/s,1/t) = st\cdot\frac{1}{t} = s = \min(s,t) E [ X s X t ] = s t E [ B 1/ s B 1/ t ] = s t ⋅ min ( 1/ s , 1/ t ) = s t ⋅ t 1 = s = min ( s , t ) です。残るのは t ↓ 0 t \downarrow 0 t ↓ 0 での連続性、すなわち lim t ↓ 0 X t = 0 \lim_{t\downarrow0}X_t = 0 lim t ↓ 0 X t = 0 を確率 1 で示すことです。X X X と B B B は ( 0 , ∞ ) (0,\infty) ( 0 , ∞ ) 上で同じ有限次元分布を持ちますから、可算個の時刻だけで決まる事象については両者の確率が一致します。事象
A : = { lim t ↓ 0 , t ∈ Q > 0 X t = 0 } A := \Bigl\{\lim_{t\downarrow0,\ t\in\mathbb{Q}_{>0}} X_t = 0\Bigr\} A := { t ↓ 0 , t ∈ Q > 0 lim X t = 0 } は正の有理数時刻という可算集合の値だけで決まるので、P ( A ) = P ( lim t ↓ 0 , t ∈ Q > 0 B t = 0 ) = 1 P(A) = P\bigl(\lim_{t\downarrow0,\,t\in\mathbb{Q}_{>0}}B_t = 0\bigr) = 1 P ( A ) = P ( lim t ↓ 0 , t ∈ Q > 0 B t = 0 ) = 1 です。最後の等号は B B B の道の連続性と B 0 = 0 B_0=0 B 0 = 0 によります。A A A の上では、X X X が ( 0 , ∞ ) (0,\infty) ( 0 , ∞ ) で連続であることと合わせて lim t ↓ 0 X t = 0 \lim_{t\downarrow0}X_t = 0 lim t ↓ 0 X t = 0 が従います。実際、任意の t > 0 t>0 t > 0 に対し t t t に十分近い有理数 r r r を取れば ∣ X t − X r ∣ |X_t - X_r| ∣ X t − X r ∣ を小さくできるからです。
∎ Corollary 5.7 (ブラウン運動の大数の法則 )
標準ブラウン運動 ( B t ) (B_t) ( B t ) について lim t → ∞ B t t = 0 \displaystyle\lim_{t\to\infty}\frac{B_t}{t} = 0 t → ∞ lim t B t = 0 が確率 1 で成り立つ。
Proof(Corollary 5.7) Proposition 5.6 (d) の過程 X u = u B 1 / u X_u = uB_{1/u} X u = u B 1/ u はブラウン運動であり、特に u ↓ 0 u \downarrow 0 u ↓ 0 で X u → X 0 = 0 X_u \to X_0 = 0 X u → X 0 = 0 が確率 1 で成り立ちます。u = 1 / t u = 1/t u = 1/ t とおけば X 1 / t = B t / t X_{1/t} = B_t/t X 1/ t = B t / t であり、t → ∞ t\to\infty t → ∞ は u ↓ 0 u\downarrow0 u ↓ 0 に対応するので、B t / t → 0 B_t/t \to 0 B t / t → 0 が確率 1 で従います。
∎ ブラウン運動の道がどれほど粗いかを、二次変分が定量的に教えてくれます。
Theorem 5.8 (ブラウン運動の二次変分 )
( B t ) (B_t) ( B t ) を標準ブラウン運動、t > 0 t > 0 t > 0 を固定する。二進分割 t k n : = k t / 2 n t^n_k := kt/2^n t k n := k t / 2 n (k = 0 , 1 , … , 2 n k = 0,1,\ldots,2^n k = 0 , 1 , … , 2 n )に沿う二次変分を
Q n : = ∑ k = 1 2 n ( B t k n − B t k − 1 n ) 2 Q_n := \sum_{k=1}^{2^n}\bigl(B_{t^n_k} - B_{t^n_{k-1}}\bigr)^2 Q n := k = 1 ∑ 2 n ( B t k n − B t k − 1 n ) 2 と定める。このとき E [ Q n ] = t E[Q_n] = t E [ Q n ] = t かつ Var ( Q n ) = 2 t 2 / 2 n \operatorname{Var}(Q_n) = 2t^2/2^n Var ( Q n ) = 2 t 2 / 2 n であり、Q n → t Q_n \to t Q n → t が L 2 L^2 L 2 の意味でも確率 1 でも成り立つ。
Proof(Theorem 5.8) h : = t / 2 n h := t/2^n h := t / 2 n 、Δ k : = B t k n − B t k − 1 n \Delta_k := B_{t^n_k}-B_{t^n_{k-1}} Δ k := B t k n − B t k − 1 n とおきます。Definition 5.1 の条件 2, 3 より Δ 1 , … , Δ 2 n \Delta_1,\ldots,\Delta_{2^n} Δ 1 , … , Δ 2 n は独立で、各々 N ( 0 , h ) N(0,h) N ( 0 , h ) に従います。
期待値は E [ Δ k 2 ] = h E[\Delta_k^2] = h E [ Δ k 2 ] = h なので E [ Q n ] = 2 n h = t E[Q_n] = 2^n h = t E [ Q n ] = 2 n h = t です。
分散を計算します。Z ∼ N ( 0 , h ) Z \sim N(0,h) Z ∼ N ( 0 , h ) に対し E [ Z 4 ] = 3 h 2 E[Z^4] = 3h^2 E [ Z 4 ] = 3 h 2 (正規分布の 4 次モーメント)ですから
Var ( Δ k 2 ) = E [ Δ k 4 ] − ( E [ Δ k 2 ] ) 2 = 3 h 2 − h 2 = 2 h 2 . \operatorname{Var}(\Delta_k^2) = E[\Delta_k^4] - \bigl(E[\Delta_k^2]\bigr)^2 = 3h^2 - h^2 = 2h^2 . Var ( Δ k 2 ) = E [ Δ k 4 ] − ( E [ Δ k 2 ] ) 2 = 3 h 2 − h 2 = 2 h 2 . Δ k \Delta_k Δ k たちは独立なので Δ k 2 \Delta_k^2 Δ k 2 たちも独立であり、分散は加法的です。よって
Var ( Q n ) = 2 n ⋅ 2 h 2 = 2 n ⋅ 2 ⋅ t 2 2 2 n = 2 t 2 2 n . \operatorname{Var}(Q_n) = 2^n\cdot 2h^2 = 2^n\cdot 2\cdot\frac{t^2}{2^{2n}} = \frac{2t^2}{2^n} . Var ( Q n ) = 2 n ⋅ 2 h 2 = 2 n ⋅ 2 ⋅ 2 2 n t 2 = 2 n 2 t 2 . したがって E [ ( Q n − t ) 2 ] = Var ( Q n ) → 0 E[(Q_n - t)^2] = \operatorname{Var}(Q_n) \to 0 E [( Q n − t ) 2 ] = Var ( Q n ) → 0 、すなわち Q n → t Q_n \to t Q n → t が L 2 L^2 L 2 ノルムの意味で成り立ちます(L^p 空間と関数解析への導入 を参照)。
確率 1 での収束を示します。任意の ε > 0 \varepsilon > 0 ε > 0 に対しチェビシェフの不等式より
P ( ∣ Q n − t ∣ > ε ) ≤ Var ( Q n ) ε 2 = 2 t 2 ε 2 2 n P(|Q_n - t| > \varepsilon) \le \frac{\operatorname{Var}(Q_n)}{\varepsilon^2} = \frac{2t^2}{\varepsilon^2 2^n} P ( ∣ Q n − t ∣ > ε ) ≤ ε 2 Var ( Q n ) = ε 2 2 n 2 t 2 であり、右辺は n n n について総和可能です。ボレル–カンテリの補題(Lemma 6.5)[Probability Spaces and Kolmogorov's Axioms] より P ( ∣ Q n − t ∣ > ε が無限個の n で起こる ) = 0 P(|Q_n-t|>\varepsilon \text{ が無限個の } n \text{ で起こる}) = 0 P ( ∣ Q n − t ∣ > ε が無限個の n で起こる ) = 0 です。ε = 1 / m \varepsilon = 1/m ε = 1/ m (m ∈ N m\in\mathbb{N} m ∈ N )について可算個の共通部分を取れば、確率 1 で Q n → t Q_n \to t Q n → t を得ます。
∎ 滑らかな関数なら二次変分は 0 0 0 になります。実際 f f f が [ 0 , t ] [0,t] [ 0 , t ] でリプシッツ連続(定数 L L L )なら ∑ ( f ( t k n ) − f ( t k − 1 n ) ) 2 ≤ 2 n L 2 h 2 = L 2 t 2 / 2 n → 0 \sum(f(t^n_k)-f(t^n_{k-1}))^2 \le 2^n L^2h^2 = L^2t^2/2^n \to 0 ∑ ( f ( t k n ) − f ( t k − 1 n ) ) 2 ≤ 2 n L 2 h 2 = L 2 t 2 / 2 n → 0 です。ブラウン運動の二次変分が 0 0 0 でない正の値 t t t に収束するというのは、道が微分可能な関数とは質的に違うことを示しています。その最も重要な帰結が次です。
Corollary 5.9 (道は有界変動でない )
t > 0 t>0 t > 0 を固定する。確率 1 で、見本道 s ↦ B s ( ω ) s\mapsto B_s(\omega) s ↦ B s ( ω ) は [ 0 , t ] [0,t] [ 0 , t ] 上で有界変動ではない。すなわち
sup { ∑ k = 1 m ∣ B s k − B s k − 1 ∣ : m ∈ N , 0 = s 0 < s 1 < ⋯ < s m = t } = + ∞ \sup\Bigl\{\sum_{k=1}^{m}|B_{s_k}-B_{s_{k-1}}| \ :\ m\in\mathbb{N},\ 0 = s_0 < s_1 < \cdots < s_m = t\Bigr\} = +\infty sup { k = 1 ∑ m ∣ B s k − B s k − 1 ∣ : m ∈ N , 0 = s 0 < s 1 < ⋯ < s m = t } = + ∞ が確率 1 で成り立つ。
Proof(Corollary 5.9) Ω 0 \Omega_0 Ω 0 を「道が連続であり、かつ Theorem 5.8 の意味で Q n → t Q_n\to t Q n → t となる」ω \omega ω の全体とします。P ( Ω 0 ) = 1 P(\Omega_0) = 1 P ( Ω 0 ) = 1 です。
ω ∈ Ω 0 \omega\in\Omega_0 ω ∈ Ω 0 を取り、この道が [ 0 , t ] [0,t] [ 0 , t ] 上で有界変動であると仮定して矛盾を導きます。全変動を V < ∞ V < \infty V < ∞ とします。道は [ 0 , t ] [0,t] [ 0 , t ] 上連続なので、コンパクト集合上の連続関数として一様連続です(連続関数と一様連続性 の Theorem 7.1[Continuous Functions and Uniform Continuity] を参照)。したがって分割の幅 t / 2 n → 0 t/2^n \to 0 t / 2 n → 0 に対して
m n : = max 1 ≤ k ≤ 2 n ∣ B t k n ( ω ) − B t k − 1 n ( ω ) ∣ ⟶ 0 ( n → ∞ ) m_n := \max_{1\le k\le 2^n}\bigl|B_{t^n_k}(\omega)-B_{t^n_{k-1}}(\omega)\bigr| \longrightarrow 0 \quad (n\to\infty) m n := 1 ≤ k ≤ 2 n max B t k n ( ω ) − B t k − 1 n ( ω ) ⟶ 0 ( n → ∞ ) です。すると
Q n ( ω ) = ∑ k ∣ B t k n − B t k − 1 n ∣ 2 ≤ m n ∑ k ∣ B t k n − B t k − 1 n ∣ ≤ m n V ⟶ 0 Q_n(\omega) = \sum_{k}\bigl|B_{t^n_k}-B_{t^n_{k-1}}\bigr|^2 \le m_n\sum_{k}\bigl|B_{t^n_k}-B_{t^n_{k-1}}\bigr| \le m_n V \longrightarrow 0 Q n ( ω ) = k ∑ B t k n − B t k − 1 n 2 ≤ m n k ∑ B t k n − B t k − 1 n ≤ m n V ⟶ 0 となります。ここで中央の不等式は各項の絶対値の一方を最大値で抑えたもの、右の不等式は全変動の定義(二進分割も分割の一つ)です。一方 ω ∈ Ω 0 \omega\in\Omega_0 ω ∈ Ω 0 より Q n ( ω ) → t > 0 Q_n(\omega)\to t > 0 Q n ( ω ) → t > 0 でしたから矛盾します。よって Ω 0 \Omega_0 Ω 0 上のすべての ω \omega ω で道は有界変動ではありません。
∎ 以下 ( F t ) t ≥ 0 (\mathcal{F}_t)_{t\ge0} ( F t ) t ≥ 0 は F t = σ ( B s : 0 ≤ s ≤ t ) \mathcal{F}_t = \sigma(B_s : 0\le s\le t) F t = σ ( B s : 0 ≤ s ≤ t ) 、すなわちブラウン運動の自然なフィルトレーションとします。Definition 5.1 の条件 2 から、t > s t > s t > s のとき増分 B t − B s B_t - B_s B t − B s は F s \mathcal{F}_s F s と独立です。これは、F s \mathcal{F}_s F s が B u B_{u} B u (u ≤ s u\le s u ≤ s )で生成され、それらが増分 B u i − B u i − 1 B_{u_i}-B_{u_{i-1}} B u i − B u i − 1 (u i ≤ s u_i \le s u i ≤ s )の線形結合として書けることによります。
Theorem 6.1 (ブラウン運動から作られるマルチンゲール )
( B t ) t ≥ 0 (B_t)_{t\ge0} ( B t ) t ≥ 0 を標準ブラウン運動、( F t ) (\mathcal{F}_t) ( F t ) をその自然なフィルトレーションとする。次の三つはいずれも ( F t ) (\mathcal{F}_t) ( F t ) に関するマルチンゲールである。
(a) B t B_t B t 自身。
(b) B t 2 − t B_t^2 - t B t 2 − t 。
(c) 任意に固定した実数 λ \lambda λ に対し、Z t λ : = exp ( λ B t − 1 2 λ 2 t ) Z^\lambda_t := \exp\bigl(\lambda B_t - \tfrac{1}{2}\lambda^2 t\bigr) Z t λ := exp ( λ B t − 2 1 λ 2 t ) (指数マルチンゲール )。
Proof(Theorem 6.1) 0 ≤ s < t 0 \le s < t 0 ≤ s < t とし、Δ : = B t − B s \Delta := B_t - B_s Δ := B t − B s とおきます。Definition 5.1 の条件 3 より Δ ∼ N ( 0 , t − s ) \Delta \sim N(0,t-s) Δ ∼ N ( 0 , t − s ) で、上で述べたとおり Δ \Delta Δ は F s \mathcal{F}_s F s と独立です。適合性はいずれの過程も B t B_t B t の可測関数なので明らかです。
(a) 可積分性は E ∣ B t ∣ = 2 t / π < ∞ E|B_t| = \sqrt{2t/\pi} < \infty E ∣ B t ∣ = 2 t / π < ∞ です。B t = B s + Δ B_t = B_s + \Delta B t = B s + Δ と分解し、B s B_s B s が F s \mathcal{F}_s F s -可測、Δ \Delta Δ が F s \mathcal{F}_s F s と独立であることを使うと
E [ B t ∣ F s ] = B s + E [ Δ ∣ F s ] = B s + E [ Δ ] = B s + 0 = B s . E[B_t\mid\mathcal{F}_s] = B_s + E[\Delta\mid\mathcal{F}_s] = B_s + E[\Delta] = B_s + 0 = B_s . E [ B t ∣ F s ] = B s + E [ Δ ∣ F s ] = B s + E [ Δ ] = B s + 0 = B s . (b) 可積分性は E [ B t 2 ] = t < ∞ E[B_t^2] = t < \infty E [ B t 2 ] = t < ∞ です。B t 2 = ( B s + Δ ) 2 = B s 2 + 2 B s Δ + Δ 2 B_t^2 = (B_s+\Delta)^2 = B_s^2 + 2B_s\Delta + \Delta^2 B t 2 = ( B s + Δ ) 2 = B s 2 + 2 B s Δ + Δ 2 を条件付けると
E [ B t 2 ∣ F s ] = B s 2 + 2 B s E [ Δ ] + E [ Δ 2 ] = B s 2 + 0 + ( t − s ) . E[B_t^2\mid\mathcal{F}_s] = B_s^2 + 2B_s\,E[\Delta] + E[\Delta^2] = B_s^2 + 0 + (t-s) . E [ B t 2 ∣ F s ] = B s 2 + 2 B s E [ Δ ] + E [ Δ 2 ] = B s 2 + 0 + ( t − s ) . 既知量の取り出し(B s B_s B s と B s 2 B_s^2 B s 2 は F s \mathcal{F}_s F s -可測)と独立性を使いました。両辺から t t t を引くと E [ B t 2 − t ∣ F s ] = B s 2 − s E[B_t^2 - t\mid\mathcal{F}_s] = B_s^2 - s E [ B t 2 − t ∣ F s ] = B s 2 − s です。
(c) 正規分布 N ( 0 , σ 2 ) N(0,\sigma^2) N ( 0 , σ 2 ) のモーメント母関数は E [ e λ Z ] = e λ 2 σ 2 / 2 E[e^{\lambda Z}] = e^{\lambda^2\sigma^2/2} E [ e λ Z ] = e λ 2 σ 2 /2 です。まず可積分性は
E [ Z t λ ] = e − λ 2 t / 2 E [ e λ B t ] = e − λ 2 t / 2 e λ 2 t / 2 = 1 < ∞ E\bigl[Z^\lambda_t\bigr] = e^{-\lambda^2t/2}E[e^{\lambda B_t}] = e^{-\lambda^2t/2}e^{\lambda^2t/2} = 1 < \infty E [ Z t λ ] = e − λ 2 t /2 E [ e λ B t ] = e − λ 2 t /2 e λ 2 t /2 = 1 < ∞ から従います(被積分関数が正なので絶対値を付けても同じです)。次に e λ B t = e λ B s e λ Δ e^{\lambda B_t} = e^{\lambda B_s}e^{\lambda\Delta} e λ B t = e λ B s e λ Δ と分解して
E [ e λ B t ∣ F s ] = e λ B s E [ e λ Δ ] = e λ B s e λ 2 ( t − s ) / 2 E\bigl[e^{\lambda B_t}\bigm|\mathcal{F}_s\bigr] = e^{\lambda B_s}E\bigl[e^{\lambda\Delta}\bigr] = e^{\lambda B_s}e^{\lambda^2(t-s)/2} E [ e λ B t F s ] = e λ B s E [ e λ Δ ] = e λ B s e λ 2 ( t − s ) /2 となります。両辺に e − λ 2 t / 2 e^{-\lambda^2t/2} e − λ 2 t /2 を掛けると
E [ Z t λ ∣ F s ] = e λ B s e λ 2 ( t − s ) / 2 − λ 2 t / 2 = e λ B s − λ 2 s / 2 = Z s λ E\bigl[Z^\lambda_t\bigm|\mathcal{F}_s\bigr] = e^{\lambda B_s}e^{\lambda^2(t-s)/2 - \lambda^2t/2} = e^{\lambda B_s - \lambda^2s/2} = Z^\lambda_s E [ Z t λ F s ] = e λ B s e λ 2 ( t − s ) /2 − λ 2 t /2 = e λ B s − λ 2 s /2 = Z s λ を得ます。
∎ 任意停止定理は連続時間でもほぼ同じ形で成り立ちます。証明は離散時間の結果を二進時刻に制限して適用し、道の右連続性で極限を取る、という筋です。
Theorem 6.2 (任意停止定理(連続時間版) )
( M t ) t ≥ 0 (M_t)_{t\ge0} ( M t ) t ≥ 0 を ( F t ) (\mathcal{F}_t) ( F t ) に関するマルチンゲールで、確率 1 で道が連続であるものとする。τ \tau τ を ( F t ) (\mathcal{F}_t) ( F t ) に関する停止時刻とする。このとき、
(a) 任意の定数 T ≥ 0 T \ge 0 T ≥ 0 について E [ M T ∧ τ ] = E [ M 0 ] E[M_{T\wedge\tau}] = E[M_0] E [ M T ∧ τ ] = E [ M 0 ] が成り立つ。
(b) さらに P ( τ < ∞ ) = 1 P(\tau<\infty)=1 P ( τ < ∞ ) = 1 であり、可積分な確率変数 Y ≥ 0 Y \ge 0 Y ≥ 0 が存在してすべての T ≥ 0 T \ge 0 T ≥ 0 で ∣ M T ∧ τ ∣ ≤ Y |M_{T\wedge\tau}| \le Y ∣ M T ∧ τ ∣ ≤ Y となるならば、E [ M τ ] = E [ M 0 ] E[M_\tau] = E[M_0] E [ M τ ] = E [ M 0 ] である。
Example 6.4 (ブラウン運動の出口問題 )
a , b > 0 a, b > 0 a , b > 0 を定数とし、
τ : = inf { t ≥ 0 : B t = − a または B t = b } \tau := \inf\{t \ge 0 : B_t = -a \ \text{または}\ B_t = b\} τ := inf { t ≥ 0 : B t = − a または B t = b } とおきます。道が連続なので τ \tau τ は停止時刻です。
第一段:E [ τ ] ≤ max ( a , b ) 2 E[\tau] \le \max(a,b)^2 E [ τ ] ≤ max ( a , b ) 2 。 c : = max ( a , b ) c := \max(a,b) c := max ( a , b ) とおくと、t ≤ τ t \le \tau t ≤ τ の間 ∣ B t ∣ ≤ c |B_t| \le c ∣ B t ∣ ≤ c です。Theorem 6.1 (b) と Theorem 6.2 (a) から、任意の T T T について
0 = E [ B T ∧ τ 2 − ( T ∧ τ ) ] , すなわち E [ T ∧ τ ] = E [ B T ∧ τ 2 ] ≤ c 2 . 0 = E\bigl[B_{T\wedge\tau}^2 - (T\wedge\tau)\bigr], \quad\text{すなわち}\quad E[T\wedge\tau] = E\bigl[B_{T\wedge\tau}^2\bigr] \le c^2 . 0 = E [ B T ∧ τ 2 − ( T ∧ τ ) ] , すなわち E [ T ∧ τ ] = E [ B T ∧ τ 2 ] ≤ c 2 . T ∧ τ T\wedge\tau T ∧ τ は T T T について単調増加で τ \tau τ に収束するので、単調収束定理から E [ τ ] ≤ c 2 < ∞ E[\tau] \le c^2 < \infty E [ τ ] ≤ c 2 < ∞ です。特に P ( τ < ∞ ) = 1 P(\tau<\infty)=1 P ( τ < ∞ ) = 1 です。
第二段:到達確率。 ∣ B T ∧ τ ∣ ≤ c |B_{T\wedge\tau}| \le c ∣ B T ∧ τ ∣ ≤ c なので Theorem 6.2 (b) が Y = c Y = c Y = c で使えて、E [ B τ ] = E [ B 0 ] = 0 E[B_\tau] = E[B_0] = 0 E [ B τ ] = E [ B 0 ] = 0 です。p : = P ( B τ = b ) p := P(B_\tau = b) p := P ( B τ = b ) とおくと B τ B_\tau B τ は − a -a − a か b b b しか取らないので
0 = − a ( 1 − p ) + b p = − a + ( a + b ) p ⟹ p = a a + b . 0 = -a(1-p) + bp = -a + (a+b)p \ \Longrightarrow\ p = \frac{a}{a+b} . 0 = − a ( 1 − p ) + b p = − a + ( a + b ) p ⟹ p = a + b a . 上に行くか下に行くかの確率が、逆側の壁までの距離に比例するわけです。
第三段:平均到達時間。 M t : = B t 2 − t M_t := B_t^2 - t M t := B t 2 − t に対し ∣ M T ∧ τ ∣ ≤ c 2 + τ |M_{T\wedge\tau}| \le c^2 + \tau ∣ M T ∧ τ ∣ ≤ c 2 + τ であり、第一段より c 2 + τ c^2+\tau c 2 + τ は可積分です。よって Theorem 6.2 (b) から E [ B τ 2 ] − E [ τ ] = 0 E[B_\tau^2] - E[\tau] = 0 E [ B τ 2 ] − E [ τ ] = 0 、すなわち
E [ τ ] = E [ B τ 2 ] = a 2 ( 1 − p ) + b 2 p = a 2 ⋅ b a + b + b 2 ⋅ a a + b = a b ( a + b ) a + b = a b . E[\tau] = E[B_\tau^2] = a^2(1-p) + b^2p = a^2\cdot\frac{b}{a+b} + b^2\cdot\frac{a}{a+b} = \frac{ab(a+b)}{a+b} = ab . E [ τ ] = E [ B τ 2 ] = a 2 ( 1 − p ) + b 2 p = a 2 ⋅ a + b b + b 2 ⋅ a + b a = a + b ab ( a + b ) = ab . a = b = 1 a = b = 1 a = b = 1 なら平均 1 単位時間で ± 1 \pm1 ± 1 のどちらかに到達します。Example 4.5 の答え x ( N − x ) x(N-x) x ( N − x ) で x = a x = a x = a 、N = a + b N = a+b N = a + b と置いたものに一致しており、ランダムウォークの極限としてブラウン運動を見る見方と整合しています。
Exercise 7.1 易
( X n ) n ≥ 0 (X_n)_{n\ge0} ( X n ) n ≥ 0 を ( F n ) (\mathcal{F}_n) ( F n ) に関するマルチンゲール、φ : R → R \varphi : \mathbb{R}\to\mathbb{R} φ : R → R を凸関数とし、すべての n n n について E [ ∣ φ ( X n ) ∣ ] < ∞ E[|\varphi(X_n)|] < \infty E [ ∣ φ ( X n ) ∣ ] < ∞ が成り立つとする。このとき ( φ ( X n ) ) n ≥ 0 (\varphi(X_n))_{n\ge0} ( φ ( X n ) ) n ≥ 0 は劣マルチンゲールであることを示せ。さらにこれを使って、∣ X n ∣ |X_n| ∣ X n ∣ と X n 2 X_n^2 X n 2 (後者は E [ X n 2 ] < ∞ E[X_n^2]<\infty E [ X n 2 ] < ∞ のとき)が劣マルチンゲールであることを示せ。
Solution 適合性は、φ \varphi φ が凸ゆえ連続、したがってボレル可測であることから従います。X n X_n X n が F n \mathcal{F}_n F n -可測なら φ ( X n ) \varphi(X_n) φ ( X n ) も F n \mathcal{F}_n F n -可測です。可積分性は仮定です。
条件付き期待値に対するイェンセンの不等式を G = F n \mathcal{G}=\mathcal{F}_n G = F n に適用すると
φ ( E [ X n + 1 ∣ F n ] ) ≤ E [ φ ( X n + 1 ) ∣ F n ] \varphi\bigl(E[X_{n+1}\mid\mathcal{F}_n]\bigr) \le E\bigl[\varphi(X_{n+1})\bigm|\mathcal{F}_n\bigr] φ ( E [ X n + 1 ∣ F n ] ) ≤ E [ φ ( X n + 1 ) F n ] です。左辺は Definition 3.1 の条件 3 より φ ( X n ) \varphi(X_n) φ ( X n ) に等しいので、E [ φ ( X n + 1 ) ∣ F n ] ≥ φ ( X n ) E[\varphi(X_{n+1})\mid\mathcal{F}_n] \ge \varphi(X_n) E [ φ ( X n + 1 ) ∣ F n ] ≥ φ ( X n ) 、すなわち劣マルチンゲールです。
φ ( u ) = ∣ u ∣ \varphi(u)=|u| φ ( u ) = ∣ u ∣ は凸で、E ∣ X n ∣ < ∞ E|X_n|<\infty E ∣ X n ∣ < ∞ は Definition 3.1 の条件 2 そのものなので、( ∣ X n ∣ ) (|X_n|) ( ∣ X n ∣ ) は劣マルチンゲールです。φ ( u ) = u 2 \varphi(u)=u^2 φ ( u ) = u 2 も凸で、E [ X n 2 ] < ∞ E[X_n^2]<\infty E [ X n 2 ] < ∞ を仮定すれば ( X n 2 ) (X_n^2) ( X n 2 ) も劣マルチンゲールです。なお ( X n 2 ) (X_n^2) ( X n 2 ) がマルチンゲールになるとは限りません。Example 3.3 では E [ S n + 1 2 ∣ F n ] = S n 2 + σ 2 > S n 2 E[S_{n+1}^2\mid\mathcal{F}_n] = S_n^2+\sigma^2 > S_n^2 E [ S n + 1 2 ∣ F n ] = S n 2 + σ 2 > S n 2 (σ > 0 \sigma>0 σ > 0 のとき)で、真に劣マルチンゲールです。
Exercise 7.2 標準
ξ 1 , ξ 2 , … \xi_1,\xi_2,\ldots ξ 1 , ξ 2 , … を P ( ξ i = 1 ) = P ( ξ i = − 1 ) = 1 / 2 P(\xi_i=1)=P(\xi_i=-1)=1/2 P ( ξ i = 1 ) = P ( ξ i = − 1 ) = 1/2 の独立同分布列、S 0 = 0 S_0=0 S 0 = 0 、S n = ξ 1 + ⋯ + ξ n S_n=\xi_1+\cdots+\xi_n S n = ξ 1 + ⋯ + ξ n 、F n = σ ( ξ 1 , … , ξ n ) \mathcal{F}_n=\sigma(\xi_1,\ldots,\xi_n) F n = σ ( ξ 1 , … , ξ n ) とする。M n : = S n 3 − 3 n S n M_n := S_n^3 - 3nS_n M n := S n 3 − 3 n S n が ( F n ) (\mathcal{F}_n) ( F n ) に関するマルチンゲールであることを示せ。
Solution ∣ S n ∣ ≤ n |S_n|\le n ∣ S n ∣ ≤ n なので ∣ M n ∣ ≤ n 3 + 3 n 2 < ∞ |M_n| \le n^3 + 3n^2 < \infty ∣ M n ∣ ≤ n 3 + 3 n 2 < ∞ であり、有界な確率変数として可積分です。適合性も M n M_n M n が S n S_n S n と定数の多項式であることから明らかです。
ξ : = ξ n + 1 \xi := \xi_{n+1} ξ := ξ n + 1 とおくと ξ 2 = 1 \xi^2 = 1 ξ 2 = 1 (確率 1)、E [ ξ ] = 0 E[\xi]=0 E [ ξ ] = 0 、E [ ξ 3 ] = E [ ξ ] = 0 E[\xi^3] = E[\xi] = 0 E [ ξ 3 ] = E [ ξ ] = 0 です。S n + 1 = S n + ξ S_{n+1}=S_n+\xi S n + 1 = S n + ξ を三乗すると
S n + 1 3 = S n 3 + 3 S n 2 ξ + 3 S n ξ 2 + ξ 3 = S n 3 + 3 S n 2 ξ + 3 S n + ξ . S_{n+1}^3 = S_n^3 + 3S_n^2\xi + 3S_n\xi^2 + \xi^3 = S_n^3 + 3S_n^2\xi + 3S_n + \xi . S n + 1 3 = S n 3 + 3 S n 2 ξ + 3 S n ξ 2 + ξ 3 = S n 3 + 3 S n 2 ξ + 3 S n + ξ . ξ \xi ξ は F n \mathcal{F}_n F n と独立、S n S_n S n と S n 2 S_n^2 S n 2 は F n \mathcal{F}_n F n -可測なので
E [ S n + 1 3 ∣ F n ] = S n 3 + 3 S n 2 E [ ξ ] + 3 S n + E [ ξ ] = S n 3 + 3 S n . E[S_{n+1}^3\mid\mathcal{F}_n] = S_n^3 + 3S_n^2E[\xi] + 3S_n + E[\xi] = S_n^3 + 3S_n . E [ S n + 1 3 ∣ F n ] = S n 3 + 3 S n 2 E [ ξ ] + 3 S n + E [ ξ ] = S n 3 + 3 S n . 一方、Example 3.3 より E [ S n + 1 ∣ F n ] = S n E[S_{n+1}\mid\mathcal{F}_n]=S_n E [ S n + 1 ∣ F n ] = S n なので E [ 3 ( n + 1 ) S n + 1 ∣ F n ] = 3 ( n + 1 ) S n E[3(n+1)S_{n+1}\mid\mathcal{F}_n] = 3(n+1)S_n E [ 3 ( n + 1 ) S n + 1 ∣ F n ] = 3 ( n + 1 ) S n です。差を取ると
E [ M n + 1 ∣ F n ] = S n 3 + 3 S n − 3 ( n + 1 ) S n = S n 3 − 3 n S n = M n E[M_{n+1}\mid\mathcal{F}_n] = S_n^3 + 3S_n - 3(n+1)S_n = S_n^3 - 3nS_n = M_n E [ M n + 1 ∣ F n ] = S n 3 + 3 S n − 3 ( n + 1 ) S n = S n 3 − 3 n S n = M n となり、マルチンゲールです。
Exercise 7.3 標準
P ( ξ i = 1 ) = p P(\xi_i=1)=p P ( ξ i = 1 ) = p 、P ( ξ i = − 1 ) = q = 1 − p P(\xi_i=-1)=q=1-p P ( ξ i = − 1 ) = q = 1 − p (0 < p < 1 0 < p < 1 0 < p < 1 、p ≠ 1 / 2 p \ne 1/2 p = 1/2 )の独立同分布列に対し S n = x + ξ 1 + ⋯ + ξ n S_n = x + \xi_1+\cdots+\xi_n S n = x + ξ 1 + ⋯ + ξ n (x x x は 0 < x < N 0 < x < N 0 < x < N なる整数)とする。r : = q / p r := q/p r := q / p とおくとき、Y n : = r S n Y_n := r^{S_n} Y n := r S n がマルチンゲールであることを示し、Example 4.5 の τ \tau τ に対して P ( S τ = N ) P(S_\tau = N) P ( S τ = N ) を求めよ。
Solution まず E [ r ξ 1 ] = p r + q r − 1 = p ⋅ q p + q ⋅ p q = q + p = 1 E[r^{\xi_1}] = pr + qr^{-1} = p\cdot\frac{q}{p} + q\cdot\frac{p}{q} = q + p = 1 E [ r ξ 1 ] = p r + q r − 1 = p ⋅ p q + q ⋅ q p = q + p = 1 です。可積分性は、∣ S n − x ∣ ≤ n |S_n - x| \le n ∣ S n − x ∣ ≤ n より 0 < Y n = r S n ≤ r x max ( r , 1 / r ) n 0 < Y_n = r^{S_n} \le r^{x}\max(r, 1/r)^{n} 0 < Y n = r S n ≤ r x max ( r , 1/ r ) n と、各 n n n について Y n Y_n Y n が定数で抑えられることから従います。適合性は Y n Y_n Y n が S n S_n S n の関数であることから明らかです。r S n + 1 = r S n r ξ n + 1 r^{S_{n+1}} = r^{S_n}r^{\xi_{n+1}} r S n + 1 = r S n r ξ n + 1 と分解し、r S n r^{S_n} r S n が F n \mathcal{F}_n F n -可測、ξ n + 1 \xi_{n+1} ξ n + 1 が F n \mathcal{F}_n F n と独立であることから
E [ Y n + 1 ∣ F n ] = r S n E [ r ξ n + 1 ] = r S n ⋅ 1 = Y n E[Y_{n+1}\mid\mathcal{F}_n] = r^{S_n}E[r^{\xi_{n+1}}] = r^{S_n}\cdot 1 = Y_n E [ Y n + 1 ∣ F n ] = r S n E [ r ξ n + 1 ] = r S n ⋅ 1 = Y n となり、マルチンゲールです。
Example 4.5 の第一段の議論(+ 1 +1 + 1 が N N N 回連続する確率は p N > 0 p^N > 0 p N > 0 に置き換わる)から P ( τ < ∞ ) = 1 P(\tau<\infty)=1 P ( τ < ∞ ) = 1 が同様に従います。n ≤ τ n\le\tau n ≤ τ の間 0 ≤ S n ≤ N 0\le S_n\le N 0 ≤ S n ≤ N なので ∣ Y n ∧ τ ∣ ≤ max ( 1 , r N ) |Y_{n\wedge\tau}| \le \max(1, r^N) ∣ Y n ∧ τ ∣ ≤ max ( 1 , r N ) と有界で、Theorem 4.3 の (b) が定数 Y = max ( 1 , r N ) Y=\max(1,r^N) Y = max ( 1 , r N ) で使えます。よって E [ Y τ ] = E [ Y 0 ] = r x E[Y_\tau] = E[Y_0] = r^x E [ Y τ ] = E [ Y 0 ] = r x 、すなわち u : = P ( S τ = N ) u := P(S_\tau=N) u := P ( S τ = N ) とおくと
r x = r 0 ( 1 − u ) + r N u = 1 + ( r N − 1 ) u . r^{x} = r^{0}(1-u) + r^{N}u = 1 + (r^N-1)u . r x = r 0 ( 1 − u ) + r N u = 1 + ( r N − 1 ) u . r ≠ 1 r \ne 1 r = 1 より r N ≠ 1 r^N \ne 1 r N = 1 なので
P ( S τ = N ) = u = r x − 1 r N − 1 = ( q / p ) x − 1 ( q / p ) N − 1 . P(S_\tau = N) = u = \frac{r^{x}-1}{r^{N}-1} = \frac{(q/p)^{x}-1}{(q/p)^{N}-1} . P ( S τ = N ) = u = r N − 1 r x − 1 = ( q / p ) N − 1 ( q / p ) x − 1 . p → 1 / 2 p \to 1/2 p → 1/2 の極限ではロピタルの定理により x / N x/N x / N に収束し、Example 4.5 の結果と整合します。p < 1 / 2 p < 1/2 p < 1/2 (不利なゲーム)で N N N が大きいと r > 1 r > 1 r > 1 より u ≈ r x − N → 0 u \approx r^{x-N} \to 0 u ≈ r x − N → 0 となり、成功確率は指数関数的に小さくなります。
Exercise 7.4 難
( B t ) (B_t) ( B t ) を標準ブラウン運動、a > 0 a > 0 a > 0 とし、τ a : = inf { t ≥ 0 : B t = a } \tau_a := \inf\{t\ge0 : B_t = a\} τ a := inf { t ≥ 0 : B t = a } とおく。θ > 0 \theta > 0 θ > 0 に対して
E [ e − θ τ a 1 { τ a < ∞ } ] = e − a 2 θ E\bigl[e^{-\theta\tau_a}\mathbf{1}_{\{\tau_a<\infty\}}\bigr] = e^{-a\sqrt{2\theta}} E [ e − θ τ a 1 { τ a < ∞ } ] = e − a 2 θ を示し、さらに P ( τ a < ∞ ) = 1 P(\tau_a<\infty)=1 P ( τ a < ∞ ) = 1 を導け。
Solution λ : = 2 θ > 0 \lambda := \sqrt{2\theta} > 0 λ := 2 θ > 0 とし、Theorem 6.1 (c) の指数マルチンゲール Z t = exp ( λ B t − λ 2 t / 2 ) Z_t = \exp(\lambda B_t - \lambda^2t/2) Z t = exp ( λ B t − λ 2 t /2 ) を考えます。λ 2 / 2 = θ \lambda^2/2 = \theta λ 2 /2 = θ なので Z t = exp ( λ B t − θ t ) Z_t = \exp(\lambda B_t - \theta t) Z t = exp ( λ B t − θ t ) です。
t ≤ τ a t \le \tau_a t ≤ τ a の間 B t ≤ a B_t \le a B t ≤ a なので、λ > 0 \lambda>0 λ > 0 と θ > 0 \theta>0 θ > 0 より
0 < Z t ∧ τ a = exp ( λ B t ∧ τ a − θ ( t ∧ τ a ) ) ≤ e λ a 0 < Z_{t\wedge\tau_a} = \exp\bigl(\lambda B_{t\wedge\tau_a} - \theta(t\wedge\tau_a)\bigr) \le e^{\lambda a} 0 < Z t ∧ τ a = exp ( λ B t ∧ τ a − θ ( t ∧ τ a ) ) ≤ e λa と一様に有界です。定数 Y = e λ a Y = e^{\lambda a} Y = e λa は可積分なので、τ a < ∞ \tau_a<\infty τ a < ∞ が確率 1 とはまだ分かっていませんが、有界収束定理を直接使えます。極限を調べます。
τ a < ∞ \tau_a<\infty τ a < ∞ となる ω \omega ω では t → ∞ t\to\infty t → ∞ のとき Z t ∧ τ a → exp ( λ a − θ τ a ) Z_{t\wedge\tau_a} \to \exp(\lambda a - \theta\tau_a) Z t ∧ τ a → exp ( λa − θ τ a ) です(B τ a = a B_{\tau_a}=a B τ a = a は道の連続性による)。τ a = ∞ \tau_a=\infty τ a = ∞ となる ω \omega ω では Z t ∧ τ a = Z t = exp ( t ( λ B t / t − θ ) ) Z_{t\wedge\tau_a}=Z_t = \exp(t(\lambda B_t/t - \theta)) Z t ∧ τ a = Z t = exp ( t ( λ B t / t − θ )) であり、Corollary 5.7 より B t / t → 0 B_t/t\to0 B t / t → 0 が確率 1 で成り立つので、指数の中身は − ∞ -\infty − ∞ に発散し Z t → 0 Z_t \to 0 Z t → 0 です。まとめると確率 1 で
Z t ∧ τ a ⟶ e λ a e − θ τ a 1 { τ a < ∞ } . Z_{t\wedge\tau_a} \longrightarrow e^{\lambda a}e^{-\theta\tau_a}\mathbf{1}_{\{\tau_a<\infty\}} . Z t ∧ τ a ⟶ e λa e − θ τ a 1 { τ a < ∞ } . Theorem 6.2 (a) より各 t t t で E [ Z t ∧ τ a ] = E [ Z 0 ] = 1 E[Z_{t\wedge\tau_a}] = E[Z_0] = 1 E [ Z t ∧ τ a ] = E [ Z 0 ] = 1 なので、有界収束定理から
1 = e λ a E [ e − θ τ a 1 { τ a < ∞ } ] , すなわち E [ e − θ τ a 1 { τ a < ∞ } ] = e − λ a = e − a 2 θ . 1 = e^{\lambda a}\,E\bigl[e^{-\theta\tau_a}\mathbf{1}_{\{\tau_a<\infty\}}\bigr],
\qquad\text{すなわち}\qquad
E\bigl[e^{-\theta\tau_a}\mathbf{1}_{\{\tau_a<\infty\}}\bigr] = e^{-\lambda a} = e^{-a\sqrt{2\theta}} . 1 = e λa E [ e − θ τ a 1 { τ a < ∞ } ] , すなわち E [ e − θ τ a 1 { τ a < ∞ } ] = e − λa = e − a 2 θ . 次に θ ↓ 0 \theta \downarrow 0 θ ↓ 0 とします。被積分関数 e − θ τ a 1 { τ a < ∞ } e^{-\theta\tau_a}\mathbf{1}_{\{\tau_a<\infty\}} e − θ τ a 1 { τ a < ∞ } は θ \theta θ について単調増加で 1 1 1 で抑えられ、各点で 1 { τ a < ∞ } \mathbf{1}_{\{\tau_a<\infty\}} 1 { τ a < ∞ } に収束します。有界収束定理より左辺は P ( τ a < ∞ ) P(\tau_a<\infty) P ( τ a < ∞ ) に収束し、右辺 e − a 2 θ e^{-a\sqrt{2\theta}} e − a 2 θ は 1 1 1 に収束します。よって P ( τ a < ∞ ) = 1 P(\tau_a<\infty)=1 P ( τ a < ∞ ) = 1 です。
なお、この結果から E [ τ a ] = + ∞ E[\tau_a] = +\infty E [ τ a ] = + ∞ も従います。P ( τ a < ∞ ) = 1 P(\tau_a<\infty)=1 P ( τ a < ∞ ) = 1 が分かったので示した等式は E [ e − θ τ a ] = e − a 2 θ E[e^{-\theta\tau_a}] = e^{-a\sqrt{2\theta}} E [ e − θ τ a ] = e − a 2 θ と書けます。u ↦ e − θ u u \mapsto e^{-\theta u} u ↦ e − θ u は凸なのでイェンセンの不等式より E [ e − θ τ a ] ≥ e − θ E [ τ a ] E[e^{-\theta\tau_a}] \ge e^{-\theta E[\tau_a]} E [ e − θ τ a ] ≥ e − θ E [ τ a ] であり(E [ τ a ] = ∞ E[\tau_a]=\infty E [ τ a ] = ∞ のときは右辺を 0 0 0 と読みます)、両辺の対数を取ると − a 2 θ ≥ − θ E [ τ a ] -a\sqrt{2\theta} \ge -\theta E[\tau_a] − a 2 θ ≥ − θ E [ τ a ] 、すなわち
E [ τ a ] ≥ a 2 θ θ = a 2 θ E[\tau_a] \ge \frac{a\sqrt{2\theta}}{\theta} = \frac{a\sqrt{2}}{\sqrt{\theta}} E [ τ a ] ≥ θ a 2 θ = θ a 2 です。θ ↓ 0 \theta\downarrow0 θ ↓ 0 とすれば右辺は + ∞ +\infty + ∞ に発散するので E [ τ a ] = + ∞ E[\tau_a]=+\infty E [ τ a ] = + ∞ です。確率 1 で到達するのに平均到達時間は無限大、というのがブラウン運動の再帰性の特徴です。
D. Williams, Probability with Martingales , Cambridge University Press, 1991 — 第 10〜14 章。離散時間マルチンゲールの標準的な入門書で、任意停止定理と収束定理の扱いが丁寧です。
R. Durrett, Probability: Theory and Examples , 5th ed., Cambridge University Press, 2019 — 第 4 章(マルチンゲール)、第 7 章(ブラウン運動)。
I. Karatzas and S. E. Shreve, Brownian Motion and Stochastic Calculus , 2nd ed., Springer, 1991 — 第 1 章(連続時間マルチンゲール)、第 2 章(ブラウン運動の構成)、第 3 章(レヴィの特徴づけ)。
P. Mörters and Y. Peres, Brownian Motion , Cambridge University Press, 2010 — 第 1 章(構成と道の性質)、第 2 章(強マルコフ性)。
舟木直久『確率論』朝倉書店、2004 — 測度論的確率論からマルチンゲール、ブラウン運動までを日本語で通した標準的な教科書です。
N. Wiener, “Differential-Space”, Journal of Mathematics and Physics 2 (1923), 131–174 — ブラウン運動の測度論的構成を与えた原論文。
方針。 Theorem 5.2 の証明として標準的なレヴィ–チースィエルスキの構成を、[ 0 , 1 ] [0,1] [ 0 , 1 ] 上に限って概観します。考え方は「独立な正規確率変数の列を係数とする関数級数を作り、それが一様収束することを示す」というものです。一様収束極限は連続関数なので、Definition 5.1 の条件 4 が構成から自動的に出る点が要点です。
ハール関数とシャウダー関数。 H 0 , 0 ( t ) : = 1 H_{0,0}(t) := 1 H 0 , 0 ( t ) := 1 とし、n ≥ 0 n \ge 0 n ≥ 0 と 1 ≤ k ≤ 2 n + 1 1 \le k \le 2^{n+1} 1 ≤ k ≤ 2 n + 1 が奇数のときに、区間 [ ( k − 1 ) 2 − n − 1 , k 2 − n − 1 ) [(k-1)2^{-n-1}, k2^{-n-1}) [( k − 1 ) 2 − n − 1 , k 2 − n − 1 ) で 2 n / 2 2^{n/2} 2 n /2 、次の同じ長さの区間で − 2 n / 2 -2^{n/2} − 2 n /2 、他で 0 0 0 を取る関数を並べると、これらは L 2 ( [ 0 , 1 ] ) L^2([0,1]) L 2 ([ 0 , 1 ]) の完全正規直交系をなします(ハール系)。その原始関数 S n , k ( t ) : = ∫ 0 t H n , k ( u ) d u S_{n,k}(t) := \int_0^t H_{n,k}(u)\,du S n , k ( t ) := ∫ 0 t H n , k ( u ) d u をシャウダー関数といい、S n , k S_{n,k} S n , k は台が長さ 2 − n 2^{-n} 2 − n の区間に含まれる三角形状の関数で、その最大値は 2 − n / 2 − 1 2^{-n/2-1} 2 − n /2 − 1 です。
級数の定義。 ( Z n , k ) (Z_{n,k}) ( Z n , k ) を独立な標準正規確率変数の族とし、
B t : = ∑ n , k Z n , k S n , k ( t ) B_t := \sum_{n,k} Z_{n,k}\,S_{n,k}(t) B t := n , k ∑ Z n , k S n , k ( t ) と定めます。この級数が確率 1 で [ 0 , 1 ] [0,1] [ 0 , 1 ] 上一様収束することが証明の山場です。同じ n n n に属する S n , k S_{n,k} S n , k たちは台が交わらないので、第 n n n 段の寄与の一様ノルムは ( max k ∣ Z n , k ∣ ) ⋅ 2 − n / 2 − 1 \bigl(\max_{k}|Z_{n,k}|\bigr)\cdot 2^{-n/2-1} ( max k ∣ Z n , k ∣ ) ⋅ 2 − n /2 − 1 で抑えられます。標準正規分布の裾の評価 P ( ∣ Z ∣ > u ) ≤ e − u 2 / 2 P(|Z| > u) \le e^{-u^2/2} P ( ∣ Z ∣ > u ) ≤ e − u 2 /2 (u ≥ 1 u \ge 1 u ≥ 1 )と、第 n n n 段の項数が 2 n 2^{n} 2 n 程度であることから、u n : = c n u_n := c\sqrt{n} u n := c n を十分大きな c c c で取れば
P ( max k ∣ Z n , k ∣ > c n ) ≤ 2 n e − c 2 n / 2 P\Bigl(\max_{k}|Z_{n,k}| > c\sqrt{n}\Bigr) \le 2^{n}e^{-c^2n/2} P ( k max ∣ Z n , k ∣ > c n ) ≤ 2 n e − c 2 n /2 が総和可能になります。ボレル–カンテリの補題より、確率 1 で有限個の n n n を除いて第 n n n 段の寄与は c n 2 − n / 2 − 1 c\sqrt{n}\,2^{-n/2-1} c n 2 − n /2 − 1 以下であり、この列は総和可能です。したがって ワイエルシュトラスの M 判定法(Corollary 6.4)[関数列と一様収束] により級数は確率 1 で一様収束し(関数列と一様収束 を参照)、極限 B t B_t B t は連続関数になります。
性質の確認。 残りは Proposition 5.5 を使うのが早道です。B t B_t B t は独立な正規確率変数の(L 2 L^2 L 2 収束する)線形結合なのでガウス過程であり、平均は 0 0 0 です。共分散はハール系の完全性(パーセバルの等式)から
E [ B s B t ] = ∑ n , k S n , k ( s ) S n , k ( t ) = ∫ 0 1 1 [ 0 , s ] ( u ) 1 [ 0 , t ] ( u ) d u = min ( s , t ) E[B_sB_t] = \sum_{n,k}S_{n,k}(s)S_{n,k}(t) = \int_0^1 \mathbf{1}_{[0,s]}(u)\,\mathbf{1}_{[0,t]}(u)\,du = \min(s,t) E [ B s B t ] = n , k ∑ S n , k ( s ) S n , k ( t ) = ∫ 0 1 1 [ 0 , s ] ( u ) 1 [ 0 , t ] ( u ) d u = min ( s , t ) と計算されます。ここで S n , k ( s ) = ⟨ 1 [ 0 , s ] , H n , k ⟩ S_{n,k}(s) = \langle \mathbf{1}_{[0,s]}, H_{n,k}\rangle S n , k ( s ) = ⟨ 1 [ 0 , s ] , H n , k ⟩ であることを使いました。よって Proposition 5.5 より ( B t ) t ∈ [ 0 , 1 ] (B_t)_{t\in[0,1]} ( B t ) t ∈ [ 0 , 1 ] は標準ブラウン運動です。[ 0 , ∞ ) [0,\infty) [ 0 , ∞ ) 上のブラウン運動は、独立なコピーを [ 0 , 1 ] [0,1] [ 0 , 1 ] 、[ 1 , 2 ] [1,2] [ 1 , 2 ] 、… と端点をつないで並べれば得られます。