[論文レビュー] Randomized First Passage Times

Sebastian Jaimungal, A. Kreinin+1 more|ArXiv.org|Nov 21, 2009

History and advancements in chemistry参考文献 6被引用数 4

ひとこと要約

本稿は、ランダムな初期位置 $X$ を用いてブラウン運動の境界 $b(t)$ への最初の到達時間をキャリブレーションし、到達時間分布が目的の測度 $u$ に一致するようにするランダム化最初到達時間問題（RFPT）を導入する。主な貢献は、フリードホルム積分方程式による解の特徴付けと、$X$ のモーメント生成関数に対する新規なラプラス変換公式であり、線形境界では存在・一意性が保証され、ガンマ分布の混合分布へも拡張可能である。

ABSTRACT

In this article we study a problem related to the first passage and inverse first passage time problems for Brownian motions originally formulated by Jackson, Kreinin and Zhang (2009). Specifically, define $τ_X = \inf\{t>0:W_t + X \le b(t) \}$ where $W_t$ is a standard Brownian motion, then given a boundary function $b:[0,\infty) o \RR$ and a target measure $μ$ on $[0,\infty)$, we seek the random variable $X$ such that the law of $τ_X$ is given by $μ$. We characterize the solutions, prove uniqueness and existence and provide several key examples associated with the linear boundary.

研究の動機と目的

• ランダム初期位置 $X$ を用いたブラウン運動の最初の到達時間 $\tau_X$ が、指定された確率測度 $\mu$ に一致するように、ランダム化最初到達時間問題（RFPT）を定式化・解明すること。
• 一般の境界関数 $b(t)$ に対して、$\tau_X \sim \mu$ となるような $X$ の分布の存在と一意性を確立すること。
• 目的の到達時間分布が有限個のガンマ分布の混合である場合、初期位置分布 $X$ の解析的解を提供すること。
• 初期条件をランダム化することで、古典的な最初到達時間問題を拡張し、逆問題へのより広範な適用可能性を実現すること。
• 境界関数と目的密度に依存する $X$ のモーメント生成関数のラプラス変換公式を導出し、一般の境界関数に対して数値的逆変換が可能となるようにすること。

提案手法

• RFPT問題を第一種フリードホルム積分方程式として定式化する：$\int_0^\infty f(t|x)g(x)\,dx = f(t)$、ここで $f(t|x)$ は条件付き到達時間密度、$g(x)$ は初期位置 $X$ の密度関数である。
• 定理3.2を用いて $X$ の存在を証明し、$\tilde{g}(\alpha) = \int_0^\infty e^{-\alpha b(t) - \alpha^2 t/2} f(t)\,dt$ として、$X$ の分布を目的密度 $f(t)$ と境界 $b(t)$ に結びつけるモーメント生成関数を提供する。
• 線形境界 $b(t) = \mu t$ の場合、積分変換がラプラス変換に簡略化され、明示的な解の導出が可能になることを示す。
• エルミート多項式の展開と超幾何関数の恒等式を用いて、解カーネルの正値性と可積分性を分析し、混合分布のターゲットに対して存在を保証する。
• 放物型円柱関数とそのエルミート多項式との関係を用いて、特定の境界条件下での到達時間密度の解析的表現を導出する。
• 一般化逆ガウス型密度のラプラス変換を用いて、解の構造を検証し、得られた密度 $g(x)$ が非負かつ可積分であることを保証する。

実験結果

リサーチクエスチョン

• RQ1 境界 $b(t)$ に対して、最初の到達時間 $\tau_X$ が事前に指定された分布 $\mu$ に一致するようなランダム初期位置 $X$ が存在するか？
• RQ2 与えられた境界 $b(t)$ と目的測度 $\mu$ に対して、解 $X$ の分布が一意に定まる条件は何か？
• RQ3 目的分布 $\tau_X$ が有限個のガンマ分布の混合である場合、$X$ の解が解析的に表現可能か？
• RQ4 初期位置 $X$ のラプラス変換は、境界関数 $b(t)$ と目的密度 $f(t)$ とどのように関係するか？
• RQ5 特に線形境界および滑らかな境界に対して、フリードホルム積分方程式 $\int f(t|x)g(x)dx = f(t)$ の解 $g(x)$ の存在を保証する十分条件は何か？

主な発見

• 定理3.2により、任意の連続的に微分可能な境界 $b(t)$ と目的密度 $f(t)$ に対して、$\tilde{g}(\alpha) = \int_0^\infty e^{-\alpha b(t) - \alpha^2 t/2} f(t)\,dt$ として、モーメント生成関数 $\tilde{g}(\alpha)$ の閉形式表現が得られ、解 $X$ の存在が保証される。
• 線形境界 $b(t) = \mu t$ の場合、モーメント生成関数は $\tilde{g}(\alpha) = \mathcal{L}[f](\alpha\mu + \alpha^2/2)$ に簡略化され、ここで $\mathcal{L}$ はラプラス変換を表す。これにより、直接的な数値的逆変換が可能になる。
• 目的分布 $\tau_X$ が有限個のガンマ分布の混合である場合、本稿では $X$ の密度 $g(x)$ に対する明示的な解析的公式を導出しており、解けるRFPT問題のクラスが拡張される。
• 正規性条件の下で、命題3および4により、$g(x)$ がフリードホルム方程式を満たす解として一意に定まることを証明し、分布としての一意性が保証される。
• 正値性が $\mu^2 - 2/a \geq 0$ の条件下で保証され、これにより導出された密度 $g(x)$ が非負かつ可積分であることが保証される。
• 本稿では、$X$ が非確率的（決定的）である場合、古典的な最初到達時間問題が特殊ケースとして回復されることを示している。