[论文解读] Randomized First Passage Times

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

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

一句话总结

本文引入了随机首通时（RFPT）问题，其中使用一个随机初始位置 $X$ 来校准布朗运动首次击中边界 $b(t)$ 的 hitting time，使得击中时间分布与目标测度 $u$ 匹配。主要贡献是通过弗雷德霍姆积分方程对解进行表征，并提出一种新颖的矩生成函数的拉普拉斯变换公式，确保线性边界下解的存在性与唯一性，并可推广至伽马分布的混合分布。

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.

研究动机与目标

• 制定并求解随机首通时（RFPT）问题，其中布朗运动在随机初始位置 $X$ 下的首次击中时间被校准为一个目标概率测度 $\mu$。
• 在一般边界函数 $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 是否存在一个随机初始位置 $X$，使得布朗运动对边界 $b(t)$ 的首次击中时间 $\tau_X$ 具有预设分布 $\mu$？
• 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 通过提供矩生成函数 $\tilde{g}(\alpha)$ 的闭式表达式 $\int_0^\infty e^{-\alpha b(t) - \alpha^2 t/2} f(t)\,dt$，建立了解 $X$ 的存在性，该表达式对任意连续可微边界 $b(t)$ 和目标密度 $f(t)$ 均有效。
• 对于线性边界 $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$ 为非随机时，经典首通时问题作为特例被恢复，此时 RFPT 退化为具有确定初始条件的标准 FPT 问题。