確率微分方程式(SDE )は、1つ以上の項が確率過程である微分方程式であり、[1]結果として得られる解も確率過程です。SDEは純粋数学全体で多くの用途があり、株価、[2]ランダム成長モデル[3]や熱変動の影響を受ける物理システムなどの確率モデルのさまざまな動作をモデル化するために使用されます。
SDE にはランダム微分があり、最も基本的なケースではブラウン運動の分布微分として計算されるランダムホワイトノイズ、またはより一般的にはセミマルチンゲールです。ただし、レヴィ過程[4]のようなジャンプ過程やジャンプを含むセミマルチンゲールなど、他の種類のランダム動作も可能です。
確率微分方程式は一般に微分方程式でもランダム微分方程式でもない。ランダム微分方程式は確率微分方程式と共役である。確率微分方程式は微分多様体にも拡張できる。[5] [6] [7] [8]
背景
確率微分方程式は、 1905年の アルバート・アインシュタインと マリアン・スモルホフスキーの研究における ブラウン運動の理論に由来するが、ルイ・バシュリエは1900年にブラウン運動をモデル化した最初の人物とされ、現在バシュリエ・モデルとして知られる確率微分方程式の非常に初期の例を示した。これらの初期の例のいくつかは、フランスの物理学者ランジュバンにちなんでランジュバン方程式とも呼ばれる線型確率微分方程式であり、ランダムな力を受ける調和振動子の運動を記述する。確率微分方程式の数学的理論は、確率積分の概念を導入し、非線形確率微分方程式の研究を開始した日本の数学者伊藤清の画期的な研究を通じて1940年代に発展した。後にロシアの物理学者ストラトノビッチによって別のアプローチが提案され、通常の微積分に似た微積分につながった。
用語
文献で最も一般的な SDE の形式は、右辺がホワイト ノイズ変数に依存する項によって摂動を受けた常微分方程式です。ほとんどの場合、SDE は対応する確率差分方程式の連続時間制限として理解されています。この SDE の理解は曖昧であり、対応する積分の適切な数学的定義によって補完される必要があります。[1] [3]このような数学的定義は、1940 年代に伊藤清によって最初に提案され、今日伊藤計算として知られるものにつながりました。別の構成は後にロシアの物理学者ストラトノビッチによって提案され、ストラトノビッチ積分として知られるものにつながりました。伊藤積分とストラトノビッチ積分は関連してはいますが異なるオブジェクトであり、どちらを選択するかは検討するアプリケーションによって異なります。伊藤計算は、変数が時間であるアプリケーションでは自然な、非予測性または因果関係の概念に基づいています。一方、ストラトノビッチ計算は、通常の計算に似た規則を持ち、多様体上のランダム運動などの幾何学的問題を扱う際により自然になる固有の幾何学的性質を持っている。ただし、多様体上のランダム運動を伊藤SDE [6]を通じてモデル化することも可能であり、場合によってはその方が望ましい。たとえば、部分多様体上のSDEを最適に近似しようとする場合などである。[9]
SDE に関する別の見方は、微分同相写像の確率フローです。この理解は明確で、確率差分方程式の連続時間制限の Stratonovich バージョンに対応します。SDE に関連付けられているのは、確率分布関数の時間発展を記述する方程式であるスモルホフスキー方程式またはフォッカー・プランク方程式です。フォッカー・プランク発展を微分形式の時間発展に一般化することは、確率発展演算子の概念によって提供されます。
物理科学では、 「ランジュバン SDE」という用語の使用には曖昧さがあります。ランジュバン SDE はより一般的な形式になることもありますが、この用語は通常、勾配フロー ベクトル フィールドを持つ狭いクラスの SDE を指します。このクラスの SDE は、パリシ-スールラス確率量子化手順[10]の出発点であるため、特に人気があり、超対称量子力学に密接に関連する N = 2 超対称モデルにつながります。ただし、物理的な観点からは、このクラスの SDE はあまり興味深いものではありません。なぜなら、位相的超対称性の自発的な破れが決して示されないためです。つまり、(過剰減衰)ランジュバン SDE は決してカオス的ではありません。
確率計算
ブラウン運動あるいはウィーナー過程は、数学的に非常に複雑であることが発見されました。ウィーナー過程は、ほぼ確実にどこでも微分不可能です。[1] [3]そのため、独自の計算規則が必要です。確率微分には、伊藤確率微分法とストラトノビッチ確率微分法という2つの主要なバージョンがあります。それぞれに利点と欠点があり、初心者は特定の状況でどちらが適切であるかについてしばしば混乱します。ガイドラインが存在し (例: Øksendal, 2003) [3]、便利なことに、伊藤 SDE を同等のストラトノビッチ SDE に簡単に変換したり、その逆を行ったりすることができます。[1] [3]それでも、SDE を最初に記述するときには、どちらの計算を使用するか注意する必要があります。
数値解
確率微分方程式を解く数値解析法[11]には、オイラー・丸山法、ミルシュタイン法、ルンゲ・クッタ法(SDE)、ローゼンブロック法[12]、および反復確率積分のさまざまな表現に基づく方法などがある。[13] [14]
物理学での使用
物理学では、SDE は分子動力学から神経力学、天体物理学的物体の力学まで、幅広い応用範囲を持っています。より具体的には、SDE は量子効果が重要でないか、摂動として考慮できるすべての動的システムを記述します。SDE は、動的システム理論をノイズのあるモデルに一般化したものと考えることができます。これは重要な一般化です。なぜなら、実際のシステムは環境から完全に分離することができず、このため常に外部の確率的影響を受けるからです。
新しい未知数を導入することで、高次方程式を複数の結合された 1 次方程式に変換する標準的な手法があります。したがって、SDE の最も一般的なクラスは次のとおりです。
ここで、はシステムの位相(または状態)空間における位置、は微分可能多様体であると仮定され、 は決定論的な進化の法則を表すフロー ベクトル場、 はシステムとガウス白色ノイズ との結合を定義するベクトル場の集合です。が線形空間で が定数である場合、システムは加法性ノイズの影響を受けると言われ、それ以外の場合は乗法性ノイズの影響を受けると言われます。 この用語は、 という限定されたケースを意味しているように見えますが、一般的なケースを意味するようになったため、やや誤解を招きます。
ノイズの固定構成に対して、SDE は初期条件に関して微分可能な一意の解を持ちます。[15]確率的ケースの非自明性は、ノイズ構成にわたってさまざまな対象オブジェクトを平均化しようとすると現れます。この意味で、ノイズが乗法である場合、および SDE が確率差分方程式の連続時間制限として理解されている場合、SDE は一意に定義されたエンティティではありません。この場合、SDE は、伊藤または SDE の Stratonovich 解釈などの「SDE の解釈」と呼ばれるもので補完される必要があります。それでも、SDE を微分同相写像の連続時間確率フローと見なすと、それは確率差分方程式の連続時間制限への Stratonovich アプローチに対応する 一意に定義された数学的オブジェクトです。
物理学では、等価なフォッカー・プランク方程式(FPE)を使用して、確率分布関数を時間の関数として求めるのが主な解決法です。フォッカー・プランク方程式は決定論的偏微分方程式です。シュレーディンガー方程式が量子波動関数の時間発展を与えたり、拡散方程式が化学濃度の時間発展を与えたりするのと同様に、確率分布関数が時間とともにどのように発展するかを示します 。あるいは、モンテカルロシミュレーションによって数値解を得ることもできます。その他の手法には、統計物理学と量子力学の類似性を利用する経路積分(たとえば、フォッカー・プランク方程式は、いくつかの変数を再スケーリングすることでシュレーディンガー方程式に変換できます) や、確率分布関数の統計モーメントの常微分方程式を書き出すことなどがあります。[出典が必要]
確率と数理ファイナンスにおける使用
確率論(およびフィルタリング問題を伴う信号処理や数理ファイナンスなど、確率論の多くの応用)で使用される表記法は若干異なります。これは、確率微分方程式を解く数値法に関する出版物で使用される表記法でもあります。この表記法は、物理学の定式化における時間のランダム関数のエキゾチックな性質をより明確にします。厳密な数学的用語では、 は 通常の関数として選択することはできず、一般化関数としてのみ選択できます。数学的定式化では、この複雑さを物理学の定式化よりも曖昧さが少なく扱います。
典型的な方程式は次の形式である。
ここでウィーナー過程(標準ブラウン運動)を表す。この式は、対応する積分方程式を非公式に表現したものとして解釈される。
上記の式は、連続時間 確率過程 X tの挙動を、通常のルベーグ積分と伊藤積分の和として特徴付けています。確率微分方程式の発見的(しかし非常に役立つ)解釈は、長さδの小さな時間間隔で、確率過程X tの値が、期待値μ ( X t , t ) δと分散σ ( X t , t ) 2 δで正規分布する量だけ変化し、過程の過去の挙動とは無関係であるというものです。これは、ウィーナー過程の増分が独立しており、正規分布しているためです。関数μ はドリフト係数と呼ばれ、σは拡散係数と呼ばれます。確率過程X t は拡散過程と呼ばれ、マルコフ特性を満たします。[1]
SDE の正式な解釈は、SDE の解を構成するものに基づいて行われます。SDE の解には、強い解と弱い解という 2 つの主な定義があります[1]どちらも、 SDE の積分方程式バージョンを解くプロセスX tの存在を必要とします。この 2 つの違いは、基礎となる確率空間( ) にあります。弱い解は確率空間と積分方程式を満たすプロセスで構成されますが、強い解は方程式を満たし、特定の確率空間上で定義されるプロセスです。山田–渡辺の定理は、この 2 つを結び付けています。
重要な例としては、幾何ブラウン運動の方程式が挙げられる。
これは金融数学のブラック・ショールズオプション価格モデル[2]における株価のダイナミクスを表す方程式である。
幾何ブラウン運動を一般化すると、強い解を許容し、分布が異なる幾何ブラウン運動またはブラック・ショールズ・モデルから得られる密度の凸結合であるSDEを定義することも可能であり、解が異なるブラック・ショールズ・モデルの対数正規分布の混合ダイナミクスとして分布する単一のSDEが得られる。[2] [16] [17] [18]これにより、金融数学におけるボラティリティ・スマイルを処理できるモデルが生まれる。
より単純なSDEは算術ブラウン運動と呼ばれる[3]
1900年にルイ・バシュリエによって株価の最初のモデルとして使用され、今日ではバシュリエモデルとして知られています。
より一般的な確率微分方程式もあり、係数μとσ はプロセスX tの現在の値だけでなく、プロセスの以前の値、場合によっては他のプロセスの現在または以前の値にも依存します。その場合、解のプロセスXはマルコフ過程ではなく、拡散過程ではなく伊藤過程と呼ばれます。係数がXの現在および過去の値のみに依存する場合、定義方程式は確率遅延微分方程式と呼ばれます。
フィスク・ストラトノビッチ積分を伴う確率微分方程式をジャンプ付き半マルチンゲールに一般化したもの、それがマーカス型SDEである。マーカス積分はマクシェーンの確率微分法の拡張である。[19]
確率的金融における革新的な応用は、オルンシュタイン・ウーレンベック過程の方程式の使用から生まれる。
これは、株価のリターンが対数正規分布を示すという仮説の下での株価のリターンのダイナミクスを表す方程式である。この仮説の下で、マルチェロ・ミネナによって開発された方法論は、市場の濫用現象を隠す可能性のある異常なリターンを特定できる予測区間を決定する。[20] [21]
多様体上のSDE
より一般的には、確率微分法の理論を微分多様体上に拡張することができ、この目的のためにフィスク・ストラトノビッチ積分が使用される。多様体、有限次元ベクトル空間、通常の条件を満たすフィルタリングされた確率空間を考え、1点コンパクト化を、可測であるとする。書かれた 確率微分方程式は、
はペアであり、
- は連続値セミマルチンゲールであり、
- は上のベクトル束の準同型です。
それぞれについて、マップは線形であり、それぞれについてです。
初期条件を持つのSDEの解は、寿命 まで連続した適応型 値プロセスであり、各テスト関数に対してプロセスは実数値のセミマルチンゲールであり、各停止時間に対して方程式
ほぼ確実に、における微分が成り立ちます。寿命が最大であれば、つまり、
-ほぼ確実です。各テスト関数がセミマルチンゲールであるという事実から、 は上のセミマルチンゲールであることがわかります。最大解が与えられれば、 の時間をまで延長することができ、を まで延長すると、
区別できないプロセスまで。[22] ストラトノビッチSDEは、連鎖律を満たし、ドリフト係数と拡散係数が座標変更時にベクトル場として振舞うことを考えると、多様体上のSDEとしては自然な選択ですが、多様体上の伊藤計算の方が望ましい場合もあります。多様体上の伊藤計算の理論は、シュワルツ写像の概念を通じてローラン・シュワルツによって最初に開発されました。 [6]ジェットバンドルに基づく多様体上の伊藤SDEの関連する2ジェット解釈も参照してください。[8]この解釈は、大きな空間で与えられたSDEの解を、その空間のサブ多様体で与えられたSDEの解で最適に近似しようとするときに役立ちます。 [9]ストラトノビッチベースの射影は最適にならないという点で。これはフィルタリング問題に適用され、最適な射影フィルタにつながっています。[9]
荒れた道として
通常、SDEの解は確率的な設定を必要とする。なぜなら、解に暗黙的に含まれる積分は確率積分だからである。微分方程式をパスごとに扱うことができれば、確率積分を定義する必要はなく、確率理論とは独立に理論を展開できる。これは、SDEを考慮することを示している。
は任意の に対して単一の決定論的微分方程式として定義されます。ここで は与えられた確率空間 ( ) 内のサンプル空間です。しかし、ブラウン運動の経路は無制限に変化し、どこでも確率 1 で微分可能ではないため、SDE の直接的な経路ごとの解釈は不可能であり、 のような項に意味を与える単純な方法はなく、確率積分を任意の単一の に対する積分として単純に経路ごとに定義することも不可能になります。しかし、Wong-Zakai の結果[23]に倣い、正規ノイズを含む SDE の解の極限についてラフパス理論を使用し、ブラウン運動の反復積分の選択定義を追加することで、任意の に対して、たとえば反復ブラウン積分の特定の選択に対して確率 1 で伊藤積分と一致する決定論的な大まかな積分を定義することができます。 [23]反復積分の他の定義は、ストラトノビッチ積分のようなさまざまな確率積分の決定論的経路ごとの同等物につながります。これは例えば金融数学において、確率なしでオプションの価格を決めるために使われてきました。[24]
解決策の存在と一意性
決定論的常微分方程式や偏微分方程式と同様に、与えられた SDE に解があるかどうか、またそれが一意であるかどうかを知ることは重要です。以下は、n次元ユークリッド空間R nに値を持ち、 m次元ブラウン運動Bによって駆動される伊藤 SDE の典型 的な存在定理と一意性定理です。証明は Øksendal (2003、§5.2) にあります。[3]
T > 0とし 、
定数CとDが存在する測定可能な関数であり、
すべてのt ∈ [0, T ] およびすべてのxおよびy ∈ R nに対して、ここで
Z を、 B s、s ≥ 0によって生成されるσ代数から独立し 、有限の2次モーメントを持つランダム変数とします。
すると、確率微分方程式/初期値問題
は、P-ほぼ確実に唯一のt連続解 ( t , ω ) ↦ X t ( ω ) を持ち、X はZとB s、s ≤ tによって生成されるフィルタリングF t Zに適合し、
一般的なケース: 局所リプシッツ条件と最大解
上記の確率微分方程式は、より一般的な形式の特殊なケースに過ぎない。
どこ
- は連続半マルチンゲールであり、は 連続半マルチンゲールである。
- は、ある空でない開集合 からの写像であり、 はからへのすべての線型写像の空間です。
より一般的には、多様体上の確率微分方程式を見ることもできます。
この方程式の解が爆発するかどうかは の選択に依存する。 が何らかの局所リプシッツ条件、すなわちおよび何らかのコンパクト集合と何らかの定数 を満たすと仮定すると、条件
ここで はユークリッドノルムです。この条件は、いわゆる最大解の存在と一意性を保証します。
が連続で、上記の局所リプシッツ条件を満たすと仮定し、 を何らかの初期条件、つまり初期σ代数に関して測定可能な関数とします。 をほぼ確実に予測可能な停止時間とします。値セミマルチンゲールは の 最大解と呼ばれます。
生涯 にわたって
- 停止したプロセスを宣言することは、停止した確率微分方程式の解である。
- 集合上では、ほぼ確実に となる。[25]
いわゆる爆発時間でもあります。
明示的に解ける例
明示的に解けるSDEには以下のものがある: [11]
線形SDE: 一般的なケース
どこ
簡約可能なSDE: ケース1
与えられた微分可能関数に対しては、ストラトノビッチSDEと同等である。
一般解がある
どこ
簡約可能なSDE: ケース2
与えられた微分可能関数に対しては、ストラトノビッチSDEと同等である。
これは次のように要約できる。
ここで、は前と同じように定義されます。その一般解は
SDE と超対称性
SDE の超対称理論では、確率動力学はモデルの位相空間上の微分形式に作用する確率的発展演算子を介して定義されます。確率動力学のこの正確な定式化では、すべての SDE が位相空間の連続性が連続時間フローによって保存されることを表す位相的超対称性を持ちます。この超対称性の自発的な崩壊は、カオス、乱流、自己組織化臨界性などとして分野を超えて知られている普遍的な動的現象の数学的本質であり、ゴールドストーン定理は関連する長距離動的挙動、つまりバタフライ効果、1/fおよびパチパチ音、および地震、神経雪崩、太陽フレアなどのスケールフリー統計を説明します。
参照
参考文献
- ^ abcdef Rogers, LCG ; Williams, David (2000).拡散、マルコフ過程、マルチンゲール、第2巻:伊藤微積分(第2版、ケンブリッジ数学図書館編)。ケンブリッジ大学出版局。doi :10.1017/ CBO9780511805141。ISBN 0-521-77594-9. OCLC 42874839.
- ^ abc Musiela, M.、および Rutkowski, M. (2004)、Martingale Methods in Financial Modelling、第 2 版、Springer Verlag、ベルリン。
- ^ abcdefg Øksendal, Bernt K. (2003).確率微分方程式:応用入門. ベルリン: Springer. ISBN 3-540-04758-1。
- ^ Kunita, H. (2004). レヴィ過程と微分同相写像の確率フローに基づく確率微分方程式。Rao, MM (eds) 実数および確率解析。Trends in Mathematics。Birkhäuser Boston。https://doi.org/10.1007/978-1-4612-2054-1_6
- ^ イムケラー、ピーター; シュマルフス、ビョルン (2001)。「確率微分方程式とランダム微分方程式の共役性とグローバルアトラクターの存在」。ダイナミクスと微分方程式ジャーナル。13 (2): 215–249。doi : 10.1023 /a: 1016673307045。ISSN 1040-7294。S2CID 3120200 。
- ^ abc Michel Emery (1989). 多様体における確率計算. Springer Berlin, ハイデルベルク. Doi https://doi.org/10.1007/978-3-642-75051-9
- ^ Zdzisław Brzeźniak およびKD Elworthy、「バナッハ多様体上の確率微分方程式」、Methods Funct. Anal. Topology 6 (2000)、第1号、43-84ページ。
- ^ ab Armstrong J. およびBrigo D. (2018)。ジェットとしての固有確率微分方程式。Proc. R. Soc. A.、474: 20170559、http://doi.org/10.1098/rspa.2017.0559
- ^ abc Armstrong, J., Brigo, D. and Rossi Ferrucci, E. (2019), 部分多様体上の SDE の最適近似: Itô ベクトル投影と Itô ジェット投影。Proc. London Math. Soc., 119: 176-213. https://doi.org/10.1112/plms.12226.
- ^ Parisi, G.; Sourlas, N. (1979). 「ランダム磁場、超対称性、および負の次元」. Physical Review Letters . 43 (11): 744–745. Bibcode :1979PhRvL..43..744P. doi :10.1103/PhysRevLett.43.744.
- ^ ab Kloeden, PE, Platen E. (1992). 確率微分方程式の数値解法。Springer、ベルリン、ハイデルベルク。DOI: https://doi.org/10.1007/978-3-662-12616-5
- ^ Artemiev, SS, Averina, TA (1997). 常微分方程式と確率微分方程式の数値解析。VSP、ユトレヒト、オランダ。DOI: https://doi.org/10.1515/9783110944662
- ^ Kuznetsov, DF (2023). 反復伊藤およびストラトノビッチ確率積分の強近似: 一般化多重フーリエ級数法。伊藤SDEおよび半線形SPDEの数値積分への応用。Differ. Uravn. Protsesy Upr., no. 1. DOI: https://doi.org/10.21638/11701/spbu35.2023.110
- ^ Rybakov, KA (2023). 反復確率積分のスペクトル表現と非線形確率ダイナミクスのモデリングへの応用。数学、vol. 11、4047。DOI: https://doi.org/10.3390/math11194047
- ^ Slavík, A. (2013). 「一般化微分方程式:初期条件とパラメータに関する解の微分可能性」.数学解析および応用ジャーナル. 402 (1): 261–274. doi : 10.1016/j.jmaa.2013.01.027 .
- ^ Fengler, MR (2005)、インプライドボラティリティのセミパラメトリックモデリング、Springer Verlag、ベルリン。DOI https://doi.org/10.1007/3-540-30591-2
- ^ Brigo, Damiano ; Mercurio, Fabio (2002). 「対数正規混合ダイナミクスと市場ボラティリティスマイルへの調整」International Journal of Theoretical and Applied Finance . 5 (4): 427–446. doi :10.1142/S0219024902001511.
- ^ Brigo, D, Mercurio, F, Sartorelli, G. (2003). 代替資産価格ダイナミクスとボラティリティ スマイル、QUANT FINANC、2003 年、第 3 巻、173 - 183 ページ、ISSN 1469-7688
- ^ スティーブン・マーカス (1981)、「セミマルチゲールによる確率微分方程式のモデリングと近似」、Stochastics、第 4 巻、223 ~ 245 ページ
- ^ 「市場の濫用の検出」。リスク誌。2004年11月2日。
- ^ 「金融市場における市場濫用の検出:定量的アプローチ」。Consob – イタリア証券取引委員会。
- ^ ヴォルフガング、ハッケンブロッホ;タルマイヤー、アントン (1994)。確率分析: Eine Einführung in die Theorie der stetigen Semimartingale (ドイツ語)。 Vieweg+Teubner Verlag ヴィースバーデン。 p. 364-365。ISBN 978-3-519-02229-9。
- ^ ab Friz, P. および Hairer, M. (2020). ラフパスのコースと規則性構造入門、第 2 版、Springer-Verlag、ハイデルベルク、DOI https://doi.org/10.1007/978-3-030-41556-3
- ^ Armstrong, J., Bellani, C., Brigo, D. and Cass, T. (2021). 確率を考慮しないオプション価格設定モデル: ラフパスアプローチ。数学ファイナンス、第31巻、1494~1521ページ。
- ^ ヴォルフガング、ハッケンブロッホ;タルマイヤー、アントン (1994)。確率分析: Eine Einführung in die Theorie der stetigen Semimartingale (ドイツ語)。 Vieweg+Teubner Verlag ヴィースバーデン。 297~299ページ。ISBN 978-3-519-02229-9。
さらに読む
- エヴァンス、ローレンス C (2013)。確率微分方程式入門、アメリカ数学会。
- アドミアン、ジョージ (1983)。確率システム。科学と工学における数学 (169)。フロリダ州オーランド: Academic Press Inc.
- アドミアン、ジョージ (1986)。非線形確率作用素方程式。フロリダ州オーランド:Academic Press Inc. ISBN 978-0-12-044375-8。
- アドミアン、ジョージ (1989)。非線形確率システム理論と物理学への応用。数学とその応用 (46)。ドルドレヒト: Kluwer Academic Publishers Group。
- Calin, Ovidiu (2015)。確率微積分学の非公式入門とその応用。シンガポール:World Scientific Publishing。p. 315。ISBN 978-981-4678-93-3。
- Teugels, J.; Sund, B. 編 (2004).保険数理科学百科事典. チチェスター: Wiley. pp. 523–527.
- CW Gardiner (2004)。確率論的手法ハンドブック:物理学、化学、自然科学向け。Springer。p. 415。
- Thomas Mikosch (1998)。『初等確率微積分学:金融の観点から』シンガポール:World Scientific Publishing。p. 212。ISBN 981-02-3543-7。
- Seifedine Kadry (2007)。「線形確率微分方程式の解」。Wseas Transactions on Mathematics。米国:WSEAS TRANSACTIONS on MATHEMATICS、2007 年 4 月:618。ISSN 1109-2769 。
- Higham., Desmond J. (2001 年 1 月). 「確率微分方程式の数値シミュレーションへのアルゴリズム入門」. SIAM Review . 43 (3): 525–546. Bibcode :2001SIAMR..43..525H. CiteSeerX 10.1.1.137.6375 . doi :10.1137/S0036144500378302.
- Desmond HighamとPeter Kloeden:「確率微分方程式の数値シミュレーション入門」、SIAM、ISBN 978-1-611976-42-7(2021)。
