MathLabs

確率と統計

ブラウン運動と確率解析

不規則な運動をモデル化する連続な確率過程で、金融や物理学で使われる独自の解析学を持つ。

直観水の中で踊る花粉の粒は、なぜ独自の種類の微積分を必要とするのだろうか?

1827年、植物学者ロバート・ブラウンは顕微鏡下で水に浮遊する花粉の粒を観察し、それらが絶え間なく震え、瞬間ごとに方向を変え、彼には見分けられる規則性が全くないことに気づいた。この説明は何十年も後になってようやく理解された:各粒は毎秒何十億個もの目に見えない水分子によってあらゆる方向から叩かれており、各瞬間の正味のランダムな衝撃が粒を完全に予測不可能な方向へわずかに押しやる。ブラウン運動とは、この終わりのない、方向性のない震えの数学的な理想化である:連続な——決して跳躍しない——ランダムな経路でありながら、あまりにギザギザであるため、どの瞬間にもきちんと定義された速度を決して持たない。経路が連続でありながらどこも滑らかでないというこの一つの事実こそが、数学者にまったく新しい微積分を構築することを強いる。なぜなら通常の微積分は、曲線を拡大すれば直線のように見え得ることを前提としているからである。

ある時刻におけるブラウン運動の位置のガウス型確率密度曲線。
任意の時刻 t>0t > 0 において、標準ブラウン運動の位置 WtW_t はガウス釣鐘型密度 N(0,t)\mathcal{N}(0, t) に従い、確率 ≈68%\approx 68\% で [−t,+t][-\sqrt{t}, +\sqrt{t}] に入る。

上の分布曲線は、ある時刻に粒子がどこにいるかを表す。2次元ブラウン運動はこの過程を独立に2つ組み合わせたもので、時間を第3の座標にすると、1つの標本経路はギザギザした時空間の曲線になる。

2つの独立なブラウン座標を時間に対して描いたインタラクティブな3D離散標本経路。ステップ数を調節できる。
時間に対する2つの独立なブラウン座標の決定論的な標本経路。x と y は位置、z は時間を表す。ステップ数を変えて、異なる解像度の再現可能な離散標本を比べよう。

大学ウィーナー過程:厳密な定義

定義: 標準ウィーナー過程(ブラウン運動)

標準ウィーナー過程(またはブラウン運動)とは、連続時間確率過程 (Wt)t≥0(W_t)_{t\ge0} であって、以下を満たすものである:(1) ほとんど確実に W0=0W_0=0;(2)独立増分を持つ——任意の 0≤t1<t2<t3<⋯0\le t_1<t_2<t_3<\cdots に対して、増分 Wt2−Wt1, Wt3−Wt2,…W_{t_2}-W_{t_1},\ W_{t_3}-W_{t_2},\dots は互いに独立である;(3) 各増分は正規分布に従い、s<ts<t に対して Wt−Ws∼N(0,t−s)W_t-W_s\sim N(0,t-s);(4) ほとんどすべての標本路 t↦Wtt\mapsto W_t が連続である。性質(3)だけで、この過程はどの時間スケールでも分布として同一に見えることが強いられる——この自己相似性こそが、ブラウン運動が実に多くの異なる離散ランダムウォークの普遍的な極限として現れる理由である。

Wt−Ws∼N(0,t−s)W_t-W_s\sim N(0,t-s)

ブラウン経路の決定的で、ほとんど逆説的な特徴は、その二次変分である:[0,t][0,t] の細かい分割上で WW の増分の二乗を足し合わせると、(確率収束の意味で)分割がどれほど細かくても tt 自身へ収束する ∑i(Wti+1−Wti)2→t\sum_{i}(W_{t_{i+1}}-W_{t_i})^2 \to t——二乗された揺れは決して消えず、ちょうど tt に蓄積される。これを、任意の区間で二次変分が常に 00 になる通常の微分可能な関数 ff と比較してみよ。その増分は Δt\sqrt{\Delta t} ではなく Δt\Delta t のように縮小するからである。この一つのスケーリングの違い——ブラウン運動の増分は Δt\Delta t ではなく Δt\sqrt{\Delta t} のオーダーである——が、確率解析にまつわるあらゆる異常なことの根本原因である。

∑i(Wti+1−Wti)2→t\sum_{i}(W_{t_{i+1}}-W_{t_i})^2 \to t
通常の微積分と確率(伊藤)解析の比較
側面滑らかな関数 f(t)f(t)ブラウン経路 WtW_t
微分可能性微分可能で、接線が存在する至る所微分不可能で、どこにも接線が存在しない
Δt\Delta t にわたる増分の大きさΔt\Delta t のオーダーΔt\sqrt{\Delta t} のオーダー
[0,t][0,t] 上の二次変分00tt(決してゼロにならない)
F(f(t))F(f(t)) の連鎖律dF=F′(f)f′(t) dtdF=F'(f)f'(t)\,dt余分な二階微分の項が必要(伊藤の補題)

標準ウィーナー過程 WW と任意の t>0t>0 について、和が細分の幅が 00 に近づく [0,t][0,t] の分割にわたってとられるとき ∑i(Wti+1−Wti)2→t\sum_{i}(W_{t_{i+1}}-W_{t_i})^2 \to t が成り立ち、この極限は確率収束の意味で成立する(実際には二進分割に沿ってほとんど確実に成り立つ)。結果として、ほとんどすべての標本路 s↦Wss\mapsto W_s は s∈[0,t]s\in[0,t] のどの点においても微分可能ではない。

なぜ正しいのか?

これは確率解析をそもそも必要にする定理である:古典的な微積分で dtdt が扱われるのと同じように dWtdW_t を通常の無限小として扱うことを排除し、以下の伊藤の補題を導く発見的規則 (dWt)2=dt(dW_t)^2=dt を正当化する。

証明

二次変分。[0,t][0,t] を長さ t/nt/n の nn 個の等しい部分に分割し、Qn=∑i=1n(Wit/n−W(i−1)t/n)2Q_n=\sum_{i=1}^n(W_{it/n}-W_{(i-1)t/n})^2 とおく。各項 (Wit/n−W(i−1)t/n)2(W_{it/n}-W_{(i-1)t/n})^2 は平均 t/nt/n(増分が N(0,t/n)N(0,t/n) であるため)を持ち、正規分布の4次モーメントを用いると分散 2(t/n)22(t/n)^2 を持つ。このような独立な nn 個の項を足し合わせると、E[Qn]=tE[Q_n]=t が厳密に成り立ち、n→∞n\to\infty のとき Var(Qn)=2t2/n→0\mathrm{Var}(Q_n)=2t^2/n\to0 となるので、Qn→tQ_n\to t が確率収束(そして平均二乗収束)の意味で成り立つ——これにより二次変分の恒等式が証明される。

至る所微分不可能性。背理法により、WW がある点 ss において有限な微分係数 W′(s)=LW'(s)=L を持って微分可能であると仮定する。すると小さい hh に対して Ws+h−Ws≈LhW_{s+h}-W_s\approx Lh となるので、長さ hh の微小区間にわたる二乗増分は h2h^2 のオーダーになるはずである——しかし上の二次変分の計算は、長さ hh の区間にわたる典型的な二乗増分が hh のオーダー(小さい hh に対して h2h^2 よりずっと大きい)であることを示しており、この hh のオーダーの項を t/ht/h 個の区間にわたって足し合わせることこそが、00 ではなく有限で正の数 tt に収束する合計を生み出す。

この不一致——微分可能性は二次変分が 00 であることを強いるはずだが、実際には t>0t>0 であることが証明されている——は解消不可能であり、微分可能な点は存在し得ない。完全な測度論的議論(1933年のペイリー、ウィーナー、ジグムントによる)は、任意の点の近くで差分商が有界にとどまる確率を直接評価することでこれを厳密にし、この確率が経路上のすべての点について同時に正確にゼロであることを示す。

WtW_t を標準ウィーナー過程とし、F(t,x)F(t,x) を xx について2回、tt について1回連続微分可能な関数とする。このとき過程 F(t,Wt)F(t,W_t) は dF=(∂tF+12∂x2F)dt+∂xF dWtdF=\Big(\partial_tF+\tfrac12\partial_x^2F\Big)dt+\partial_xF\,dW_t を満たす——これは通常の連鎖律と比べて余分な二階微分の項(「伊藤補正」)を持つ確率的連鎖律である。

なぜ正しいのか?

これは確率解析において最も頻繁に使われる道具である:ランダムな経路の関数をどのように微分するかを正確に教えてくれ、ブラウン運動に依存する任意の量(オプション価格、物理的な観測量)のダイナミクスを導出する出発点となる。

証明

テイラー展開。滑らかな FF について、小さなステップ Δt\Delta t と ΔW=Wt+Δt−Wt\Delta W=W_{t+\Delta t}-W_t にわたる両変数の通常の2次テイラー展開は ΔF≈∂tF Δt+∂xF ΔW+12∂x2F (ΔW)2+12∂t2F (Δt)2+∂t∂xF Δt ΔW+⋯\Delta F\approx \partial_tF\,\Delta t+\partial_xF\,\Delta W+\tfrac12\partial_x^2F\,(\Delta W)^2+\tfrac12\partial_t^2F\,(\Delta t)^2+\partial_t\partial_xF\,\Delta t\,\Delta W+\cdots と書ける——ここまでは純粋な微積分であり、任意の滑らかな経路に対して成り立つ。

大きさによって項を並べる。ΔW\Delta W は(前の定理で示したとおり Δt\Delta t ではなく)Δt\sqrt{\Delta t} のオーダーであるため、各項は次のようにスケールする:∂tF Δt\partial_tF\,\Delta t は Δt\Delta t のオーダー;∂xF ΔW\partial_xF\,\Delta W は Δt\sqrt{\Delta t} のオーダー(支配的な、主要なランダム項);(Δt)2(\Delta t)^2 と Δt ΔW\Delta t\,\Delta W はそれぞれ (Δt)2(\Delta t)^2 と (Δt)3/2(\Delta t)^{3/2} のオーダーで Δt\Delta t に比べて無視できる;しかし (ΔW)2(\Delta W)^2 は Δt\Delta t のオーダーであり——∂tF Δt\partial_tF\,\Delta t と同じオーダーであって、まったく無視できない。これは (Δx)2(\Delta x)^2 が常に Δx\Delta x に比べて無視できる通常の微積分とは異なる。

(ΔW)2(\Delta W)^2 をその平均で置き換える。二次変分の定理により、多くの小さなステップにわたって足し合わせると (ΔW)2(\Delta W)^2 は Δt\Delta t(その平均であり、ステップが縮小し足し合わされるにつれてその平均の周りの揺らぎは消える)のように振る舞う。これが、無限小ステップの極限における発見的な置き換え規則 (dWt)2=dt(dW_t)^2=dt の非形式的な正当化である。

極限をとる。無視できる高次の項を落とし、Δt→0\Delta t\to0 の極限で (ΔW)2→dt(\Delta W)^2\to dt を代入すると、テイラー展開はまさに主張どおりの微分形 dF=(∂tF+12∂x2F)dt+∂xF dWtdF=\Big(\partial_tF+\tfrac12\partial_x^2F\Big)dt+\partial_xF\,dW_t になる——12∂x2F dt\tfrac12\partial_x^2F\,dt の項こそが、通常の微積分では捨てられるが確率解析では保持しなければならない寄与そのものである。

大学実世界での応用と具体例

物理的な拡散(水中の花粉、固体中を広がる熱、混ざり合う気体分子)を記述するだけでなく、ブラウン運動と伊藤解析は現代の定量金融の数学的基盤となった。株価 StS_t は一般に幾何ブラウン運動、dSt=μSt dt+σSt dWtdS_t=\mu S_t\,dt+\sigma S_t\,dW_t(μ\mu は平均成長率、σ\sigma はボラティリティ)としてモデル化される。金融デリバティブ(オプションなど)の価値 V(t,St)V(t,S_t) に伊藤の補題を適用すると、有名なブラック・ショールズ偏微分方程式、∂tV+12σ2S2∂S2V+rS∂SV−rV=0\partial_tV+\tfrac12\sigma^2S^2\partial_S^2V+rS\partial_SV-rV=0 が導かれ、その解が主要な取引所すべてで取引されるオプションの公正価格を与える。同じ数学——ウィーナー過程と伊藤の補題——は、神経科学におけるニューロン膜電位のモデル、経済学における金利モデル、流体力学における乱流拡散の基盤にもなっている。

例: ブラウン変位の期待値と分散

ある粒子の位置は W0=0W_0=0 から出発する1次元の標準ウィーナー過程 WtW_t に従う。E[W5]E[W_5] と Var(W5−W2)\mathrm{Var}(W_5-W_2) を求めよ。

解答

最初の量については、s=0s=0 として Wt−Ws∼N(0,t−s)W_t-W_s\sim N(0,t-s) を用いる:W5−W0=W5W_5-W_0=W_5 は N(0,5)N(0,5) であるから、直接 E[W5]=0E[W_5]=0 となる——ブラウン運動にはドリフトがないため、期待される位置は出発点から決して離れない。

2番目の量については、s=2s=2、t=5t=5 として同じ定義性質を適用する:W5−W2∼N(0,5−2)W_5-W_2\sim N(0,5-2)、すなわち N(0,3)N(0,3)。

N(0,σ2)N(0,\sigma^2) に従う確率変数の分散は定義により σ2\sigma^2 であるから、パラメータを直接読み取ると Var(W5−W2)=3\mathrm{Var}(W_5-W_2)=3 となる。

分散が経過時間 t−s=3t-s=3 のみに依存し、出発時刻 s=2s=2 自体には依存しないことに注意せよ——これはウィーナー過程の定義に組み込まれた時間斉次性を反映している。

例: F(t,x)=x2F(t,x)=x^2 への伊藤の補題の適用

伊藤の補題を用いて、ウィーナー過程の二乗の確率微分 d(Wt2)d(W_t^2) を求め、その結果を用いて E[Wt2]=tE[W_t^2]=t を確認せよ。

解答

F(t,x)=x2F(t,x)=x^2 とおくと、∂tF=0\partial_tF=0、∂xF=2x\partial_xF=2x、∂x2F=2\partial_x^2F=2 となる。

これらの微分(x=Wtx=W_t で評価)を伊藤の補題 dF=(∂tF+12∂x2F)dt+∂xF dWtdF=\Big(\partial_tF+\tfrac12\partial_x^2F\Big)dt+\partial_xF\,dW_t に代入すると d(Wt2)=(0+12⋅2)dt+2Wt dWtd(W_t^2)=\big(0+\tfrac12\cdot2\big)dt+2W_t\,dW_t となり、これは d(Wt2)=dt+2Wt dWtd(W_t^2)=dt+2W_t\,dW_t に簡約される。

両辺を 00 から tt まで積分する(W02=0W_0^2=0 を用いて):Wt2=t+2∫0tWs dWsW_t^2=t+2\int_0^tW_s\,dW_s。

期待値をとり、伊藤積分 ∫0tWs dWs\int_0^tW_s\,dW_s は常に平均ゼロであるという重要な事実(過去から独立な増分から構成されているため、蓄積するような系統的なドリフトが存在しない)を用いると:E[Wt2]=t+2E[∫0tWs dWs]=t+0=tE[W_t^2]=t+2E\big[\int_0^tW_s\,dW_s\big]=t+0=t となり、これは Wt∼N(0,t)W_t\sim N(0,t)(分散 tt)からすでに知っていたことを直接確認するものであるが、今回は定義からではなく確率解析の道具立てから純粋に導かれている。

標準ウィーナー過程において、増分 Wt−WsW_t-W_s(s<ts<t に対して)はどの分布に従うか:

[0,t][0,t] 上のブラウン運動の二次変分は:

通常の連鎖律と比べて、伊藤の補題 dF=(∂tF+12∂x2F)dt+∂xF dWtdF=\Big(\partial_tF+\tfrac12\partial_x^2F\Big)dt+\partial_xF\,dW_t に余分な項があるのは:

オプション価格付けのためのブラック・ショールズ偏微分方程式は、何を適用して導かれるか:

参考文献

  1. Ioannis Karatzas, Steven E. Shreve (1991). Brownian Motion and Stochastic Calculus
  2. Kiyosi Itô (1944). Stochastic Integral
  3. Fischer Black, Myron Scholes (1973). The Pricing of Options and Corporate Liabilities