[论文解读] $V$-Integrability, Asymptotic Stability And Comparison Theorem of Explicit Numerical Schemes for SDEs

Łukasz Szpruch, Xīlíng Zhāng|arXiv (Cornell University)|Oct 2, 2013

Stochastic processes and financial applications参考文献 27被引用 14

一句话总结

本文提出了一类针对具有非全局Lipschitz系数的随机微分方程（SDEs）的改进型显式数值格式，使得Khasminski的李雅普诺夫函数技术可应用于离散时间情形，从而建立$V$-可积性、渐近稳定性及比较定理。关键贡献在于：该格式即使在标准Euler格式失效时，仍能继承连续SDE的定性性质，如期望有界性和几乎必然收敛至平衡态。

ABSTRACT

Khasminski's \cite{chas1980stochastic} showed that many of the asymptotic stability and the integrability properties of the solutions to the Stochastic Differential Equations (SDEs) can be obtained using Lyapunov functions techniques. These properties are rarely inherited by standard numerical integrators. In this article we introduce a family of explicit numerical approximations for the SDEs and derive conditions that allow to use Khasminski's techniques in the context of numerical approximations, particularly in the case where SDEs have non globally Lipschitz coefficients. Consequently, we show that it is possible to construct a numerical scheme, that is bounded in expectation with respect to a Lyapunov function, and/or inherit the asymptotic stability property from the SDEs. Finally we show that using suitable schemes it is possible to recover comparison theorem for scalar SDEs.

研究动机与目标

• 将Khasminski的李雅普诺夫函数技术——此前仅适用于连续时间SDE——推广至离散时间数值格式。
• 解决标准显式Euler格式在系数非全局Lipschitz时无法保持可积性与稳定性的缺陷。
• 构造显式数值格式，使其在弱矩条件下方仍能继承底层SDE的$V$-可积性与渐近稳定性。
• 通过适当的格式改进，在数值设置下恢复标量SDE的比较定理。
• 提供一种使用一般李雅普诺夫函数（超越$V(x) = |x|^p$）分析数值解定性性质的框架。

提出的方法

• 引入改进的显式格式$\bar{X}_{k+1} = \bar{X}_k + b^h(t_k, \bar{X}_k)h + \sigma^h(t_k, \bar{X}_k)\Delta W_{k+1}$，通过截断漂移与扩散系数控制增长。
• 对增量项应用投影或截断机制，确保在李雅普诺夫函数$V(x) = \sum_{i=1}^d c_i x_i^{p_i}$下，即使$p_i$不相等，也能保持期望有界。
• 推导无穷小生成元$\mathcal{L}V$的离散时间类比，以模拟连续时间李雅普诺夫分析，从而实现稳定性与可积性证明。
• 利用矩估计及柯西-施瓦茨与切比雪夫不等式，控制格式在$L^{2+\epsilon}$矩条件下的弱误差与强误差。
• 对$\mathds{E}V(\bar{X}_k)$应用Gronwall型估计，依据漂移条件建立指数衰减或有界性。
• 通过误差估计验证收敛性：$\mathds{E}|\bar{X}_{t,x} - \tilde{X}_{t,x}|^2 \leq K(1 + |x|^{2+2/r})h^2$ 与 $\mathds{E}|\bar{X}_{t,x} - \mathds{E}\tilde{X}_{t,x}| \leq K(1 + |x|^{1+3/(2r)})h^{3/2}$。

实验结果

研究问题

• RQ1 能否将李雅普诺夫函数技术适配至具有非全局Lipschitz系数的SDE的显式数值格式中，以确保$V$-可积性？
• RQ2 在Euler格式失效时，显式格式在何种条件下可继承底层SDE的渐近稳定性？
• RQ3 是否可通过改进格式在数值近似设置下恢复标量SDE的比较定理？
• RQ4 如何利用一般李雅普诺夫函数$V(x) = \sum c_i x_i^{p_i}$分析超越$L^p$-范数的数值格式定性性质？
• RQ5 对$b^h$与$\sigma^h$的矩与增长条件需满足何种要求，才能在$V$-李雅普诺夫函数下保证强收敛性与期望有界性？

主要发现

• 所提出的截断显式格式在条件$\mathcal{L}V(x) \leq \rho V(x)$下，满足$\mathds{E}V(\bar{X}_k) \leq e^{\rho T} \mathds{E}V(\bar{X}_0)$，复现了连续时间下的可积性结果。
• 对于渐近稳定性，若$\mathcal{L}V(x) \leq -z(x)$且$z(x) \geq \rho V(x)$，则有$\mathds{E}V(\bar{X}_k) \leq e^{-\rho k h} \mathds{E}V(\bar{X}_0)$，证明了期望下的指数衰减。
• 该格式在$L^{4+4/r}$矩条件下实现$L^2$-收敛，误差界为$\mathds{E}|\bar{X}_{t,x} - \tilde{X}_{t,x}|^2 \leq K(1 + |x|^{2+2/r})h^2$。
• 该方法通过确保格式保持解的单调性与序关系，实现了标量SDE比较定理的恢复。
• 该格式仅需$4\kappa$阶矩即可保证收敛，与平衡格式相当，但因投影机制的存在，实际中对矩的要求更少。
• 对$\mathds{E}X_t^-$应用Gronwall不等式，可得$\mathds{E}X_t^- = 0$对所有$t$成立，表明在适当的漂移结构下，$X_t \geq 0$几乎必然成立。