前のページ 1/2 次のページ

確率微分方程式(かくりつびぶんほうていしき、英:stochastic differential equation)とは、一つ以上の項が確率過程である微分方程式であって、その結果、解自身も確率過程となるものである。

一般的に、確率微分方程式はブラウン運動(ウィーナー過程)から派生すると考えられる白色雑音を組み込むが、不連続過程の様な他の無作為変動を用いることも可能である。

◆ 背景

確率微分方程式は、ブラウン運動を記述したアインシュタインの有名な論文、および同時期に Smoluchowski により導入された。

しかし、バシュリエ(1900年)の論文「投機の理論」は、ブラウン運動に関連した初期の業績として特筆すべきである。

その後、ランジュバンに引き継がれ、後に伊藤とが確率微分方程式に数学的基礎付けを行った。

◇ 確率解析

ブラウン運動、あるいはウィーナー過程は、数学的には極めて複雑である。

ウィーナー過程の経路は微分不可能であり、従って、微分・積分を行うには、独自の規則が必要となる。

確率解析には、伊藤確率解析、ストラトノビッチ確率解析の 2 つの方法がある。

いつどちらを使うべきかは、曖昧になりがちである。

一般的に、伊藤確率微分方程式を等価なストラトノビッチ確率微分方程式に変換でき、再び戻すことも可能であり、解を得る手助けになる。

しかし、その確率微分方程式を立てた際、どちらの解析によったのか、注意を払わなければならない。

◇ 数値解

確率微分方程式、特に確率偏微分方程式の数値解法は、相対的に未発達な分野である。

常微分方程式の数値解に使用されるアルゴリズムの殆どは、確率微分方程式には殆ど有効に使用できず、数値収束が非常に悪いとされている。

◆ 定義

典型的には、B_t (t≧0) を、 B₀ = 0 を満たす連続時間一次元ブラウン運動(ウィーナー過程)とするとき、積分方程式

:$X\_-\; X\_=\; \backslash int\_t^\backslash mu(X\_u,u)\; du\; +\; \backslash int\_t^\backslash sigma(X\_u,u)\backslash ,\; dB\_u$

を

:$dX\_t\; =\; \backslash mu(X\_t,t)\backslash ,\; dt\; +\; \backslash sigma(X\_t,t)\backslash ,\; dB\_t$

の形に略記したものを、確率微分方程式という。

上記方程式は、連続時間の確率過程 X_t の振る舞いを、一般のルベーグ積分と伊藤積分の和で模している。

確率微分方程式の発見論的だがとても有益な解釈は、

微小時間間隔 δ において、確率過程 X_t の変化が、

期待値 μ(X_t,t)δ、分散 σ²(X_t,t)δ の正規分布に従って変化し、しかも過去の同確率過程の振る舞いと独立である、と見ることである。

ウィーナー過程の変化は互いに独立で正規分布に従うことから、こう考えることができる。

関数 μ(x,t) はドリフト係数(どりふとけいすう、英:drift coefficient)、関数 σ(x,t) は拡散係数(かくさんけいすう、英:diffusion coefficient)という。

確率微分方程式の解として得られる確率過程 X_t は拡散過程(かくさんかてい、英:diffusion process)と呼び、通常はマルコフ過程である。

◇ 強解と弱解

確率微分方程式の理論的解釈は、同方程式の解とは何かによって解釈する。

確率微分方程式の解の主要な定義には、強解(きょうかい、英:strong solution)と弱解(じゃくかい、英:weak solution)の二種類ある。

どちらも、確率微分方程式に対応する積分方程式の解となる確率過程 X_t の存在を要件とする。

両者の違いは、基礎となる確率空間 (Ω, F, P) にある。

弱解とは、確率積分方程式を満たす確率空間と確率過程をいい、

強解は、与えられた確率空間の上で定義され、確率積分方程式を満たす確率過程をいう。

◇ 幾何ブラウン運動

以下の確率微分方程式、

:$dX\_t\; =\; \backslash mu\; X\_t\backslash ,\; dt\; +\; \backslash sigma\; X\_t\backslash ,\; dB\_t$

は重要な例であり、この解を幾何ブラウン運動(きかぶらうんうんどう、英:geometric Brownian motion)という。

これは、数理ファイナンスにおいて、ブラック・ショールズ・オプション価格モデルで、株式価格の動きを模す方程式である。

◇ 伊藤過程

係数関数 μ と σ が、解確率過程 X_t の現在の値のみならず、同過程の過去の値、または他の確率過程の現在と過去の値にも依存する、さらに一般的な確率微分方程式が考えられる。

この場合、解確率過程 X_t はマルコフ過程ではなく、その解は拡散過程ではなく伊藤過程(いとうかてい、英:It? process)と呼ばれる。

係数関数が現在と過去の X_t の値のみに依存する場合、定義する確率微分方程式は、確率遅延微分方程式(かくりつちえんびぶんほうていしき、英:stochastic delay differential equation)という。

◆ 解の存在と一意性

決定論的な常微分方程式や偏微分方程式と同様、与えられた確率微分方程式の解が存在するか、存在するとして一意か否かを知ることは、重要である。

下記は、n 次元ユークリッド空間 Rⁿ に値を取り、m 次元ブラウン運動 B を無作為項とする伊藤確率微分方程式の解の存在および一意性に関する一般的定理である。

参考文献に記したエクセンダールの本の §5.2 には、証明が記載されている。

T > 0 とする。

:$\backslash mu\; :\; \backslash mathbb^\backslash times[外部リンク]\; T\backslash to\; \backslash mathbb^$

:$\backslash sigma\; :\; \backslash mathbb^\backslash times[外部リンク]\; T\backslash to\; \backslash mathbb^$

は可測関数で、適当な定数 C、D が存在し、

任意の t ∈ [外部リンク] T、任意の x, y ∈ Rⁿ に対し、次の 2 条件を満たすとする。

:$\backslash big|\; \backslash mu\; (x,\; t)\; \backslash big|\; +\; \backslash big|\; \backslash sigma\; (x,\; t)\; \backslash big|\; \backslash leq\; C\; \backslash big(\; 1\; +\; |\; x\; |\; \backslash big)$

:$\backslash big|\; \backslash mu\; (x,\; t)\; -\; \backslash mu\; (y,\; t)\; \backslash big|\; +\; \backslash big|\; \backslash sigma\; (x,\; t)\; -\; \backslash sigma\; (y,\; t)\; \backslash big|\; \backslash leq\; D\; |\; x\; -\; y\; |$

ここで、

:$|\; \backslash sigma\; |^=\; \backslash sum\_^|\; \backslash sigma\_|^$

である。

確率変数 Z は、_s≧0 により生成される σ 加法族と独立であり、かつ、

:

を満たすとする。

このとき、確率微分方程式、

:$$

dX_ = \mu (X_, t) dt + \sigma (X_, t) dB_ \ , 0 \le t \le T

:$X\_=\; Z\; \backslash ,$

は、以下の 2 つの性質を有する t に関して連続な解

を、P に関して殆ど確実に一意に有する。

・ X は、Z と B_s (s≦t) により生成される増大情報系

可測空間 (Ω, F) において、t ∈[0, ∞) を助変数とする部分集合族 が増大情報系(ぞうだいじょうほうけい、英:filtration)であるとは、F_t が F の部分 σ 加法族であって、かつ 0≦s<t ⇒ F_s⊆F_t を満たすことをいう。

に適合する。

増大情報系 が与えられた確率空間 (Ω, F, P) 上の確率過程 が に適合する(英:adapted)とは、任意の t に対して X_t が F_t 可測になることをいう。

・

◆ 脚注

◆ 参考文献

・ 舟木直久(2005)、確率微分方程式、岩波書店、ISBN 4-00-005196-2

前のページ 1/2 次のページ

・確率微分方程式 page1

・確率微分方程式 page2