[論文レビュー] On the First Hitting Time Density of an Ornstein-Uhlenbeck Process

Alexander Lipton, Vadim Kaushansky|arXiv (Cornell University)|Oct 4, 2018

Advanced Thermodynamics and Statistical Mechanics参考文献 23被引用数 13

ひとこと要約

本稿では、オーナイズ・ウーレンバック（OU）過程の最初の到達時刻密度を計算するための半解析的手法を提示する。問題は移動境界を伴う熱方程式に変換され、熱ポテンシャルの手法を用いて、第二種の線形ボルテラ積分方程式が導かれる。小時間ではアーベル近似を用いて解析的に解かれ、大時間では数値的に解かれる。これにより、従来の手法に比べて収束性と誤差制御が向上した、高速で正確かつ実装が容易な代替手法が得られる。

ABSTRACT

In this paper, we study the classical problem of the first passage hitting density of an Ornstein--Uhlenbeck process. We give two complementary (forward and backward) formulations of this problem and provide semi-analytical solutions for both. The corresponding problems are comparable in complexity. By using the method of heat potentials, we show how to reduce these problems to linear Volterra integral equations of the second kind. For small values of $t$, we solve these equations analytically by using Abel equation approximation; for larger $t$ we solve them numerically. We also provide a comparison with other known methods for finding the hitting density of interest, and argue that our method has considerable advantages and provides additional valuable insights.

研究の動機と目的

• オーナイズ・ウーレンバック過程の最初の到達時刻密度を、高速で正確かつ実装可能な方法で計算するための手法を開発すること。
• 収束性が悪く、計算コストが高く、実装が複雑であるといった、従来の解析的・数値的手法の限界を克服すること。
• 最初の通過問題を第二種ボルテラ積分方程式に還元する統一的な枠組みを提供する熱ポテンシャル法を用いること。
• 小時間に対しては解析的近似、大時間に対しては数値的解法を組み合わせることで、広範な適用可能性を確保すること。
• ラプラス変換、級数展開、またはモンテカルロシミュレーションなどの先行手法と比較して、より優れた誤差の定量化と制御を実現すること。

提案手法

• 変数変換を用いて、OU過程を移動境界を伴う熱方程式に変換することで問題を変形する。
• 熱ポテンシャルの手法を適用し、境界値問題を第二種の線形ボルテラ積分方程式に還元する。
• 小時間では、アーベル方程式の近似を用いてボルテラ方程式を解析的に解き、到達密度の閉形式漸近式を得る。
• 大時間では、Lipton (2018) と同様の再帰的アルゴリズムを用いてボルテラ方程式を数値的に解く。
• 熱方程式の基本解と境界上の密度関数との畳み込み型積分を用いて解を構成する。
• 数値スキームにおける項の数を制御することで誤差の定量化が可能であり、精度が保証される。

実験結果

リサーチクエスチョン

• RQ1 オーナイズ・ウーレンバック過程の最初の到達時刻密度は、半解析的手法を用いて効率的かつ正確に計算可能か？
• RQ2 ラプラス変換、級数展開、またはモンテカルロ法などの既存手法と比較して、本手法の精度と計算コストはどのように異なるか？
• RQ3 小時間と大時間の時間スケールにおける到達密度の挙動は何か？ そして、その挙動を効果的に捉えるにはどうすればよいか？
• RQ4 熱ポテンシャル法は、拡散過程における移動境界問題に効果的に適用可能か？
• RQ5 先行の解析的・数値的技術と比較して、誤差制御と収束速度はどの程度向上できるか？

主な発見

• 本手法により、OU過程の最初の通過問題が第二種ボルテラ積分方程式に還元され、体系的な解法が可能になる。
• 小時間では、アーベル方程式の近似を用いて解析的に解が得られ、到達密度の明示的漸近式が得られる。
• 大時間では、再帰的公式を用いてボルテラ方程式が数値的に解かれ、安定性と精度が保証される。
• ラプラス逆変換、級数展開、モンテカルロシミュレーションなどの競合手法と比較して、誤差制御と収束性が優れている。
• 数値的実験により、本手法の正確性と効率性が確認され、既存のベンチマークと同等またはそれ以上の結果が得られた。
• 本手法は実装が容易であり、計算の扱いやすさと信頼性の面で、他の解析的手法を上回っている。