中文

用递减步长欧拉-丸山格式逼近稳定SDE的不变测度

概率论 2023-10-10 v1 数值分析 动力系统 数值分析

摘要

(Xt)t0(X_t)_{t \ge 0} 为随机微分方程 dXt=b(Xt)dt+AdZt,X0=xdX_t = b(X_t) dt+A dZ_t, \quad X_{0}=x 的解,其中 b:RdRdb: \mathbb{R}^d \rightarrow \mathbb R^d 为Lipschitz函数,ARd×dA \in \mathbb R^{d \times d} 为正定矩阵,(Zt)t0(Z_t)_{t\geq 0}dd 维旋转不变 α\alpha-稳定 L\'evy 过程且 α(1,2)\alpha \in (1,2)xRdx\in\mathbb{R}^{d}。我们使用两种具有递减步长 Γ=(γn)nN\Gamma = (\gamma_n)_{n\in \mathbb{N}} 的欧拉-丸山格式来逼近 (Xt)t0(X_t)_{t \ge 0} 的不变测度:一种以独立同分布的 α\alpha-稳定分布随机变量为其新息,另一种以独立同分布的Pareto分布随机变量为其新息。我们研究了这两种逼近格式在Wasserstein-1距离下的收敛速率。对第一种格式,当函数 bb 为Lipschitz并满足某耗散条件时,我们证明收敛速率为 γn1/α\gamma^{1/\alpha}_n。在对 bb 的二阶方向导数附加假设下,该收敛速率可改进为对任意 κ[1,α)\kappa \in [1,\alpha)γn1+1α1κ\gamma^{1+\frac 1 {\alpha}-\frac{1}{\kappa}}_n。对第二种格式,当函数 bb 二次连续可微时,我们得到收敛速率 γn2αα\gamma^{\frac{2-\alpha}{\alpha}}_n。我们证明速率 γn2αα\gamma^{\frac{2-\alpha}{\alpha}}_n 对一维稳定Ornstein-Uhlenbeck过程是最优的。我们的定理表明,近期关于带加性新息的未调整Langevin算法的显著结果可推广到由 α\alpha-稳定 L\'evy 过程驱动的SDE,且相应收敛速率有类似表现。与先前结果相比,我们对第一种格式将二阶可微条件放宽至Lipschitz条件。

关键词

引用

@article{arxiv.2310.05390,
  title  = {Approximation of the invariant measure for stable SDE by the Euler-Maruyama scheme with decreasing step-sizes},
  author = {Peng Chen and Xinghu Jin and Yimin Xiao and Lihu Xu},
  journal= {arXiv preprint arXiv:2310.05390},
  year   = {2023}
}