《ランダム・ウォークとブラウン運動》

【らんだむ・うぉーくとぶらうんうんどう (random walk and Brownian motion) 】

　ランダム・ウォーク (random walk) とその連続化であるブラウン運動は, でたらめな動きを表現する最も基本的な確率過程で, 幅広い応用がある.

ランダム・ウォーク　$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{X_n\}_{n=1}^\infty\, } $ を互いに独立で同一の分布に従う確率変数の列とするとき,

\begin{equation}\label{B-D-05+RW}

 S_0=s~\mbox{(定数),}\qquad
 S_n = s + \sum_{i=1}^n X_i

\end{equation}

によって定義される確率過程~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{S_n\}_{n=0}^\infty\, } $ をランダム・ウォークと呼ぶ. 特に, ある $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle d>0\, } $ およびすべての $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle n\, } $ に対して, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mathrm{P}(X_n=d)=p$, $\mathrm{P}(X_n=-d)=q=1-p\, } $ であるとき, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{S_n\}_{n=0}^\infty\, } $ は (1次元の) 単純ランダム・ウォークであるといい, さらに $$p=q=1/2\,$$ のとき, 単純ランダム・ウォークは対称であるという. また, 「壁」によって動きが止められたり, 動く範囲が制限されるランダム・ウォークを考えることもできる. $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle X_n\, } $ の独立性より, ランダム・ウォークはマルコフ過程となる.

　初期値~$$s=0\,$$ のランダム・ウォークにおいて, $$n\,$$~ステップ後の位置の期待値と分散は, それぞれ $$\mathrm {E} (S_{n})=n\,\mathrm {E} (X_{1})\,$$, $$\mathrm {V} (S_{n})=n\,\mathrm {V} (X_{1})\,$$ となり, 時間の経過に比例する. 分散が時間の経過に比例することから, ランダム・ウォークは時間が経つにつれて次第に拡散していくことが分かる.

　$$d=1\,$$, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle 0<p<1\, } $ として得られる単純ランダム・ウォーク $$\{S_{n}\}_{n=0}^{\infty }\,$$ は, 整数を状態空間とする周期2の既約なマルコフ連鎖である. このマルコフ連鎖は $$p\neq 1/2\,$$ のとき一時的であり, $$p=q=1/2\,$$ ならば零再帰的となる. たとえば $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle p>1/2\, } $ ならば $$S_{n}\,$$ はだんだん大きくなっていく傾向があり, 正の方へドリフトする. このため出発点に戻ることは保証できなくなり一時的となるのである.

　2次元の対称な単純ランダム・ウォーク~(2次元格子点空間上の4つの隣接点にそれぞれ確率~$$1/4\,$$ で推移する) は零再帰的, 3次元以上の単純ランダム・ウォークはすべて一時的であることも知られている [1].

単純ランダム・ウォークからブラウン運動へ　$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{S_n\}_{n=0}^\infty\, } $ を初期値~$$s=0\$\,$ の対称な単純ランダム・ウォークとする. このランダム・ウォークが1ステップ進むのに $$T\$\,$ だけ時間がかかるとして, $$T\,$$ と $$d\,$$ を同時に0に近づけることを考える. $$t=n\,T\,$$ に対して, 時刻~$$t\,$$ にランダム・ウォークが $$x\,$$ にいる確率を $$v(x,t)\,$$ と表すと, $$v(x,t)\,$$ は差分方程式 $$v(x,t+T)=\{v(x-d,t)+v(x+d,t)\}/2\,$$ を満たすので,

 \frac{v(x,t+T) - v(x,t)}{T}
 = \frac{1}{2}\,\frac{d^2}{T}\,
   \frac{v(x+d,t) - 2\,v(x,t) + v(x-d,t)}{d^2}

が得られる. $$d^{2}/T=\sigma ^{2}\,$$~(定数) を保ったまま $$T\to 0\$(\$d\to 0\$)\,$ とすれば

\begin{equation}\label{B-D-05+Diffusion}

 \frac{\partial v(x,t)}{\partial t}
 = \frac{\sigma^2}{2}\,\frac{\partial^2 v(x,t)}{\partial x^2}

\end{equation}

を得る. 式 (2) は拡散方程式 (diffusion equation) と呼ばれ, その解は初期条件~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle v(0,0)=1\, } $, $$v(x,0)=0\$(\$x\neq 0\$)\,$ のもとで, 正規分布 $$N(0,\sigma ^{2}\,t)\,$$ の密度関数となる. より一般的には, 初期値が0の (必ずしも対称でない) 単純ランダム・ウォークにおいて, $$d^{2}/T=\sigma ^{2}\,$$, $$(p-q)/d=\mu /\sigma ^{2}\,$$ を保ったまま $$T\to 0\,$$ とすると, 時刻~$$t\$\,$ での位置が正規分布~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle N(\mu\,t,\sigma^2\,t)\, } $ に従う確率過程が得られる [1].

ブラウン運動　イギリスの植物学者ブラウン (R. Brown) は, 水面に浮く花粉中の微粒子が極めて不規則な動きをすることを見いだした. アインシュタイン (A. Einstein) は, この運動が拡散方程式 (2) によって特徴づけられることを示し, その後ウィナー (N. Wiener) らによって確率過程としての基盤が築かれた. この確率過程をブラウン運動 (Brownian motion) またはウィーナー過程 (Wiener process) と呼ぶ.

　(1次元の) ブラウン運動~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B(t)\}_{t\ge0}\, } $ は次の性質を満たす実数値確率過程である:

1. 独立増分過程である.

2. 任意の $$s\$\,$, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t>0\, } $ に対して $$B(s+t)-B(s)\,$$ は正規分布~$$N(0,\sigma ^{2}\,t)\,$$ に従う.

3. $$B(0)=0\,$$ かつ $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle B(t)\, } $ は $$t=0\,$$ で連続.

1. より, 時刻 $$s\,$$ 以降の $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B(t)\}_{t\ge s}\, } $ の振る舞いは $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle s\, } $ までの履歴には依存しないため, ブラウン運動はマルコフ過程である. さらに, ブラウン運動が強マルコフ性を持つこと, 標本路が連続となることも知られている [2].

　$$\sigma ^{2}\,$$ を拡散係数と呼び, 特に $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \sigma^2=1\, } $ のブラウン運動を標準ブラウン運動と呼ぶ. また, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle B_d(t) = \mu\,t + B(t)\, } $ によって定まる $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B_d(t)\}_{t\ge0}\, } $ をドリフトを持つブラウン運動と呼び, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mu\, } $ をドリフト係数と呼ぶ.

鏡像原理 ドリフトのないブラウン運動 $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B(t)\}_{t\ge0}\, } $ に対して $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \tau_a\, } $ を $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B(t)\}_{t\ge0}\, } $ が初めて $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle a\, } $ を横切る時刻とすると, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \tau_a\, } $ は停止時 (stopping time) となる. $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t\ge\tau_a\, } $ において $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B(t)\}_{t\ge\tau_a}\, } $ と $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle a\, } $ に関して対称な標本路を持つ確率過程~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{\bar{B}(t)\}_{t\ge0}\, } $ を

 \bar{B}(t) = \left\{\begin{array}{ll}
                B(t),       &\quad t<\tau_a, \\
                2\,a - B(t), &\quad t\ge\tau_a,
              \end{array}\right.

で定める. $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B(t)\}_{t\ge0}\, } $ が強マルコフ性を持つことと, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B(t)\}\, } $ と $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{\bar{B}(t)\}\, } $ の対称性から, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{B(t)\}\, } $ と $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{\bar{B}(t)\}\, } $ は同じ確率法則に従うことがわかる. 一般にこのような性質を鏡像原理 (reflection principle) と呼び, 初到達時間の分布などを求める際に利用される.

拡散過程 ドリフト係数や拡散係数が位置~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle x\, } $ や時刻~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle t\, } $ に依存した値~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mu(x,t)\, } $, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \sigma^2(x,t)\, } $ をとるように一般化して得られる確率過程~$構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \{D(t)\}_{t\ge0}\, } $ を拡散過程 (diffusion process) と呼び, $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \mu(x,t)\, } $ と $構文解析に失敗 (MathML、ただし動作しない場合はSVGかPNGで代替（最新ブラウザーや補助ツールに推奨）: サーバー「https://en.wikipedia.org/api/rest_v1/」から無効な応答 ("Math extension cannot connect to Restbase."):): {\displaystyle \sigma^2(x,t)\, } $ を, それぞれドリフト関数, 拡散関数と呼ぶ. 拡散過程は強マルコフ性を持ち, その標本路は連続である. 逆に, 連続な標本路を持つマルコフ過程は拡散過程となることが知られている.

　ブラウン運動や拡散過程の標本路は, 連続であるがいたるところで微分不可能という性質を持っている. このため拡散過程の解析においては, 確率積分や確率微分方程式といった通常の微分や積分とは異なる概念が必要となる [3, 4].

---

参考文献

[1] W. Feller,　An Introduction to Probability Theory and Its Applications, Volume 1, 2nd Ed., John Wiley & Sons, 1957. 河田龍夫監訳, 『確率論とその応用 I』, 紀伊国屋書店, 1960 (上巻), 1961 (下巻).

[2] K. Itô and H. P. McKean, Diffusion Processes and Their Sample Paths, Second Printing, Springer-Verlag, 1996.

[3] 木島正明, 『ファイナンス工学入門 第I部 ランダムウォークとブラウン運動』, 日科技連, 1994.

[4] 渡辺信三, 『確率微分方程式』, 産業図書, 1975.