[论文解读] Stability of densities for perturbed Diffusions and Markov Chains

Valentin Konakov, Anna Kozhina+1 more|arXiv (Cornell University)|Jun 29, 2015

Advanced Mathematical Modeling in Engineering参考文献 31被引用 3

一句话总结

本文为非退化扩散过程和马尔可夫链的转移密度在系数扰动下的稳定性建立了定量边界。通过参数化方法和霍尔德连续性假设，推导出以系数差异表示的显式误差估计，收敛速率取决于漂移和扩散系数的正则性，尤其在霍尔德或多项式衰减条件下。

ABSTRACT

We are interested in studying the sensitivity of diffusion processes or their approximations by Markov Chains with respect to a perturbation of the coefficients.

研究动机与目标

• 量化非退化扩散过程与马尔可夫链的转移密度对漂移和扩散系数扰动的敏感性。
• 在系数正则性假设最小的前提下，提供真实密度与扰动密度之间差值的显式误差边界。
• 将稳定性结果扩展至随机微分方程的离散时间马尔可夫链近似，尤其在欧拉格式弱误差分析的背景下。
• 分析系数扰动在模型误设、参数估计与数值格式收敛性等应用中的影响。
• 在霍尔德连续性与多项式衰减假设下，建立密度差值的收敛速率，显式依赖于扰动大小与正则性参数。

提出的方法

• 采用参数化方法，通过包含高斯核与迭代积分的级数展开表示扩散过程与马尔可夫链的转移密度。
• 使用麦克金–辛格方法处理参数化展开，该方法对时间光滑性不足具有鲁棒性，并可适用于离散时间马尔可夫链。
• 通过分析参数化级数项的扰动，推导出密度差值的边界，尤其关注二阶与三阶导数。
• 应用霍尔德不等式与高斯型界来控制漂移与扩散系数的差异，尤其在霍尔德连续性假设下。
• 针对具有多项式衰减的密度（如 $ q_m(t,x) riangleq c^d (t)^{-d/2} Q_{M-(d+5+ u)}(x / t^{1/2}) $）引入卷积分析，以处理卷积积分中对角与非对角区域。
• 利用数学归纳法与伽马函数恒等式控制参数化展开中迭代卷积的增长，从而得到以 $ au^{rac{r heta}{2}} $ 表示的显式收敛速率。

实验结果

研究问题

• RQ1 在漂移与扩散系数发生小扰动时，非退化扩散过程的转移密度如何变化？
• RQ2 当系数被扰动时，马尔可夫链近似密度相对于真实扩散密度的收敛速率如何？
• RQ3 系数的霍尔德连续性与多项式衰减条件如何影响密度在扰动下的稳定性？
• RQ4 参数化方法能否用于推导具有不规则系数的欧拉格式弱误差的显式误差边界？
• RQ5 高阶导数在密度差值中起什么作用？在参数化展开中如何控制浓度损失？

主要发现

• 在系数霍尔德连续性假设下，真实密度与扰动密度之间的差值被有界于 $ rac{ ilde{ ho}_{ ext{sup}}}{(t_j - t_k)^{1 - rac{ u}{2}}} $，其中 $ ilde{ ho}_{ ext{sup}} $ 捕获了系数差异的上确界范数。
• 对于马尔可夫链近似，密度差值满足 $ | ilde{p}^h oxtimes H^{h,(r)} - ilde{p}_ u^h oxtimes H_ u^{h,(r)}| riangleq ext{error} riangleq (r+1) ho_{ u, u, ext{sup}} \frac{igracevert (1 igvee T^{(1- u)/2})c_1 igracevert^{r+1} [ ext{B}( u/2)]^r}{ ext{B}(1 + r u/2)} \frac{c^d}{(t_j - t_i)^{d/2}} Q_{M-(d+5+ u)}\big( \frac{y-x}{(t_j - t_i)^{1/2}/c} \big) (t_j - t_i)^{r u/2} $ ，表明在 $ (t_j - t_i)^{r u/2} $ 下具有显式收敛性。
• 对具有多项式衰减的密度的卷积分析表明，$ I_{t_k}^1(t_i,t_j,x,y) riangleq ig| ig( q_m(t_k - t_i, ullet - x) * q_m(t_j - t_k, ullet - y) ig)(y) ig| riangleq ext{convolution} riangleq \bar{c} q_m(t_j - t_i, y - x) $，在 $ i < k < j $ 下一致成立，条件为 $ m > d $。
• 通过将定义域分割并分别积分 $ |w|^3 $ 与 $ |w|^4 $ 项，该方法控制了参数化展开中的浓度损失，从而要求 $ m > d+5+ u $ 以确保可积性。
• 密度差值的收敛速率为 $ O((t_j - t_i)^{r u/2}) $，在所假设的霍尔德正则性 $ u $ 下为最优，且边界在空间与时间上一致。
• 结果对系数扰动具有鲁棒性，体现在 $ ho_{ u, u, ext{sup}} $ 上，其度量了 $ (b, u) $ 与 $ (b_ u, u_ u) $ 之间的上确界范数差，且边界与该扰动呈线性关系。