关于Euler-Maruyama格式在多维SDEs中处理不连续漂移系数的表现研究

On the performance of the Euler-Maruyama scheme for multidimensional SDEs with discontinuous drift coefficient

摘要 Abstract

本文研究了$d$维随机微分方程(SDEs)在漂移系数不连续情形下的强逼近问题。更具体地讲,我们假设漂移系数在正达到度量的可定向$C^4$超曲面上是分段Lipschitz连续的,且例外集为$\Theta\subset \mathbb{R}^d$,扩散系数为Lipschitz连续,且在$\Theta$的邻域内,两类系数均有界,扩散系数在$\Theta$方向上有非退化的部分。近年来,文献中已经证明了一些关于此类SDEs强逼近的结果,并特别研究了Euler-Maruyama格式的表现。对于$d=1$且$\Theta$有限的情形,已知Euler-Maruyama格式实现了$L_p$误差至少为$1/2$的收敛率(对所有$p\geq 1$),与经典情况下Lipschitz连续系数的情况一致。对于$d>1$,仅知道当$\mu$和$\sigma$系数全局有界时,Euler-Maruyama格式实现$L_2$误差至少为$1/4-$的收敛率。本文证明了在上述设定下,Euler-Maruyama格式实际上对所有$d\in\mathbb{N}$和所有$p\geq 1$实现了至少为$1/2-$的$L_p$误差收敛率。该结果的证明基于经典的变换方法,即将此类SDE转化为具有全局Lipschitz连续系数的SDE,一个适用于非全局$C^2$函数的新It\^{o}公式,以及对时间连续Euler-Maruyama格式的实际位置及其在底层网格的前一时刻位置在超曲面$\Theta$两侧总期望时间的详细分析。

We study strong approximation of $d$-dimensional stochastic differential equations (SDEs) with a discontinuous drift coefficient. More precisely, we essentially assume that the drift coefficient is piecewise Lipschitz continuous with an exceptional set $\Theta\subset \mathbb{R}^d$ that is an orientable $C^4$-hypersurface of positive reach, the diffusion coefficient is assumed to be Lipschitz continuous and, in a neighborhood of $\Theta$, both coefficients are bounded and the diffusion coefficient has a non-degenerate portion orthogonal to $\Theta$. In recent years, a number of results have been proven in the literature for strong approximation of such SDEs and, in particular, the performance of the Euler-Maruyama scheme was studied. For $d=1$ and finite $\Theta$ it was shown that the Euler-Maruyama scheme achieves an $L_p$-error rate of at least $1/2$ for all $p\geq 1$ as in the classical case of Lipschitz continuous coefficients. For $d>1$, it was only known so far, that the Euler-Maruyama scheme achieves an $L_2$-error rate of at least $1/4-$ if, additionally, the coefficients $\mu$ and $\sigma$ are globally bounded. In this article, we prove that in the above setting the Euler-Maruyama scheme in fact achieves an $L_{p}$-error rate of at least $1/2-$ for all $d\in\mathbb{N}$ and all $p\geq 1$. The proof of this result is based on the well-known approach of transforming such an SDE into an SDE with globally Lipschitz continuous coefficients, a new It\^{o} formula for a class of functions which are not globally $C^2$ and a detailed analysis of the expected total time that the actual position of the time-continuous Euler-Maruyama scheme and its position at the preceding time point on the underlying grid are on 'different sides' of the hypersurface $\Theta$.