[论文解读] A Lie systems approach for the first passage-time of piecewise deterministic processes

Florin Avram, José F. Cariñena+1 more|arXiv (Cornell University)|Aug 16, 2010

Differential Equations and Boundary Problems参考文献 22被引用 3

一句话总结

本文将李系统理论应用于具有相型跳跃的连续时间生灭过程的首达时问题求解，实现了在具有马氏跳跃成分的超几何扩散模型类中 ruin 概率的显式解析解。关键贡献在于通过基于对称性的变换，系统性地将相关积分微分方程约化为可解的里卡蒂方程与常微分方程组，从而在特定漂移与跳跃结构下获得首达时与 ruin 概率的闭式表达式。

ABSTRACT

Our paper illustrates how the theory of Lie systems allows recovering known results and provide new examples of piecewise deterministic processes with phase-type jumps for which the corresponding first-time passage problems may be solved explicitly.

研究动机与目标

• 通过引入相型跳跃分布，将首达时问题的解析可解性从标准仿射模型扩展至更广范围的跳跃扩散过程。
• 探究李系统的代数几何框架是否能为具有非泊松或结构化跳跃动态的过程提供新的可 tractable 首达时问题模型。
• 证明可解李代数与相型分布的结合，可在一类具有马氏跳跃的超几何扩散过程中，为后向柯尔莫哥洛夫方程提供显式解。
• 提出一种系统化方法，通过基于对称性的重参数化与变量替换，将复杂积分微分方程转化为可解的常微分方程。
• 恢复风险理论中的已知结果（如 Cramér-Lundberg 模型），并将其推广至具有广义漂移与跳跃结构的、 ruin 概率具有解析可解性的新模型。

提出的方法

• 将李系统理论应用于控制首达时分布的后向柯尔莫哥洛夫方程，利用无穷小生成元所具有的可解李代数结构。
• 基于生成元导出的里卡蒂方程，通过状态空间的重参数化，实施变量替换，将非自治系统转化为自治系统。
• 引入函数 $ \bar{x}(x) = \frac{1}{2}\textstyle\big(\frac{\nu}{q+\nu}\big)^{1/2} \textstyle\big[\textstyle\text{log}\big(\frac{1-\textstyle\big(1-K\big)^{1/2}}{1+\textstyle\big(1-K\big)^{1/2}}\big)\big] $，实现动力学的线性化，并将问题约化为初等积分。
• 引入新变量 $ \bar{\theta}(x) $，使得里卡蒂方程变为 $ \frac{d\bar{\eta}}{dx} = \sqrt{\frac{-\lambda\mu}{\varphi_K(x)}}(1 - \bar{\eta}^2) $，其解可用双曲正切函数表示。
• 利用变换变量 $ \bar{x}(x) $，求解系统 $ \frac{dM}{dx} = (\mu \eta(x) - \mu)M $ 与 $ \Psi(x) = \eta(x)M(x) $，得到以指数与对数函数表示的显式表达式。
• 施加边界条件 $ \Psi(0) = 1 $，$ M(0) = 1 $，并导出最终解 $ K_1(K) = \frac{\sqrt{\lambda/(q+\lambda)} - \sqrt{1-K}}{\sqrt{\lambda/(q+\lambda)} + \sqrt{1-K}} $，从而实现 ruin 概率的解析计算。

实验结果

研究问题

• RQ1 李系统理论是否可用于求解具有相型跳跃的连续时间生灭过程中首达时问题，超越标准扩散或泊松跳跃模型？
• RQ2 在何种漂移与跳跃强度条件下，相关的后向柯尔莫哥洛夫方程可具备李系统结构，从而实现显式积分？
• RQ3 能否通过变量替换将首达时问题中导出的里卡蒂方程转化为自治系统，从而实现初等函数形式的解析解？
• RQ4 在负漂移与相型跳跃的所提模型中，当初始盈余趋于无穷时，ruin 概率的渐近行为如何？
• RQ5 参数 $ \lambda, \mu, q, K $ 如何影响所导解中 ruin 概率函数的形状与衰减速率？

主要发现

• 对于一类具有相型跳跃的超几何扩散过程，首达时问题可通过李系统实现解析求解，得到 ruin 概率 $ \Psi(x) $ 与相关函数 $ M(x) $ 的显式闭式表达式。
• 解通过基于对称性的变换获得，该变换重参数化状态空间，使里卡蒂方程变为自治形式，从而可通过初等积分实现求解。
• 当 $ x \to \infty $ 时，ruin 概率的渐近行为满足 $ \Psi(x) \sim C \cdot e^{\mu\big(\sqrt{\lambda/(\lambda+q)} - 1\big)x} $，当 $ \sqrt{\lambda/(\lambda+q)} < 1 $ 时指数为负，表明 $ \lim_{x\to\infty} \Psi(x) = 0 $。
• 在 $ \mu = 1.5, \lambda = q = 0.5, K = 0.75 $ 的特定情况下，模型表现出定义良好且衰减的 ruin 概率，与所导解析形式一致。
• 该方法恢复了风险理论中的已知结果（如 Cramér-Lundberg 模型），并将其推广至具有非泊松跳跃动态与结构化相型跳跃分布的新模型。
• 解显式依赖于参数 $ K_1(K) = \frac{\sqrt{\lambda/(q+\lambda)} - \sqrt{1-K}}{\sqrt{\lambda/(q+\lambda)} + \sqrt{1-K}} $，该参数控制 ruin 概率曲线的形状，并确保在 $ x = 0 $ 处满足边界条件。