Strong convergence of split-step backward Euler method for stochastic differential equations with non-smooth drift

In this paper, we are concerned with the numerical approximation of stochastic differential equations with discontinuous/nondifferentiable drifts. We show that under one-sided Lipschitz and general growth conditions on the drift and global Lipschitz condition on the diffusion, a variant of the implicit Euler method known as the split-step backward Euler (SSBE) method converges with strong order of one half to the true solution. Our analysis relies on the framework developed in [D. J. Higham, X. Mao and A. M. Stuart, Strong convergence of Euler-type methods for nonlinear stochastic differential equations, SIAM Journal on Numerical Analysis, 40 (2002) 1041-1063] and exploits the relationship which exists between explicit and implicit Euler methods to establish the convergence rate results.     

Strong convergence rates for backward Euler on a class of nonlinear jump-diffusion problems

We generalise the current theory of optimal strong convergence rates for implicit Euler-based methods by allowing for Poisson-driven jumps in a stochastic differential equation (SDE). More precisely, we show that under one-sided Lipschitz and polynomial growth conditions on the drift coefficient and global Lipschitz conditions on the diffusion and jump coefficients, three variants of backward Euler converge with strong order of one half. The analysis exploits a relation between the backward and explicit Euler methods.     

非自治刚性随机微分方程正则EM分裂方法的收敛性和稳定性

本文研究了数值求解非自治随机微分方程的正则Euler-Maruyama分裂（CEMS）方法,该方程的漂移项系数带有刚性且允许超线性增长,扩散项系数满足全局Lipschitz条件.首先,证明了CEMS方法的强收敛性及收敛速度.其次,证明了在适当条件下CEMS方法是均方稳定的.进一步,利用离散半鞅收敛定理,研究了CEMS方法的几乎必然指数稳定性.结果表明,CEMS方法在漂移系数的刚性部分满足单边Lipschitz条件下可保持几乎必然指数稳定性.最后通过数值实验,检验了CEMS方法的有效性并证实了我们的理论结果.     

Convergence of Numerical Approximation for Jump Models Involving Delay and Mean-Reverting Square Root Process

The mean-reverting square root process with jump has been widely used as a model on the financial market. Since the diffusion coefficient in the model does not satisfy the linear growth condition and local Lipschitz condition, we can not examine its properties by traditional techniques. To overcome the difficulties, we develop several new techniques to examine the numerical method of jump models involving delay and mean-reverting square root. We show that the numerical approximate solutions converge to the true solutions. Finally, we apply the convergence to examine a path-dependent option price and a bond in the financial pricing.     

Strong Convergence of Euler Approximations of Stochastic Differential Equations with Delay Under Local Lipschitz Condition

The strong convergence of Euler approximations of stochastic delay differential equations is proved under general conditions. The assumptions on drift and diffusion coefficients have been relaxed to include polynomial growth and only continuity in the arguments corresponding to delays. Furthermore, the rate of convergence is obtained under one-sided and polynomial Lipschitz conditions. Finally, our findings are demonstrated with the help of numerical simulations.     

Compensated split-step balanced methods for nonlinear stiff SDEs with jump-diffusion and piecewise continuous arguments

This paper deals with numerical solutions of nonlinear stiff stochastic differential equations with jump-diffusion and piecewise continuous arguments. By combining compensated split-step methods and balanced methods, a class of compensated split-step balanced (CSSB) methods are suggested for solving the equations. Based on the one-sided Lipschitz condition and local Lipschitz condition, a strong convergence criterion of CSSB methods is derived. It is proved under some suitable conditions that the numerical solutions produced by CSSB methods can preserve the mean-square exponential stability of the corresponding analytical solutions. Several numerical examples are presented to illustrate the obtained theoretical results and the effectiveness of CSSB methods. Moreover, in order to show the computational advantage of CSSB methods, we also give a numerical comparison with the adapted split-step backward Euler methods with or without compensation and tamed explicit methods.

     

STRONG CONVERGENCE OF A FULLY DISCRETE FINITE ELEMENT METHOD FOR A CLASS OF SEMILINEAR STOCHASTIC PARTIAL DIFFERENTIAL EQUATIONS WITH MULTIPLICATIVE NOISE

This paper develops and analyzes a fully discrete finite element method for a class of semilinear stochastic partial differential equations(SPDEs)with multiplicative noise.The nonlinearity in the diffusion term of the SPDEs is assumed to be globally Lipschitz and the nonlinearity in the drift term is only assumed to satisfy a one-sided Lipschitz condition.These assumptions are the same ones as the cases where numerical methods for general nonlinear stochastic ordinary differential equations(SODEs)under"minimum assumptions"were studied.As a result,the semilinear SPDEs considered in this paper are a direct generalization of these nonlinear SODEs.There are several difficulties which need to be overcome for this generalization.First,obviously the spatial discretization,which does not appear in the SODE case,adds an extra layer of difficulty.It turns out a spatial discretization must be designed to guarantee certain properties for the numerical scheme and its stiffness matrix.In this paper we use a finite element interpolation technique to discretize the nonlinear drift term.Second,in order to prove the strong convergence of the proposed fully discrete finite element method,stability estimates for higher order moments of the H1-seminorm of the numerical solution must be established,which are difficult and delicate.A judicious combination of the properties of the drift and diffusion terms and some nontrivial techniques is used in this paper to achieve the goal.Finally,stability estimates for the second and higher order moments of the L2-norm of the numerical solution are also difficult to obtain due to the fact that the mass matrix may not be diagonally dominant.This is done by utilizing the interpolation theory and the higher moment estimates for the H1-seminorm of the numerical solution.After overcoming these difficulties,it is proved that the proposed fully discrete finite element method is convergent in strong norms with nearly optimal rates of convergence.Numerical experiment results are also presented to validate the theoretical results and to demonstrate the efficiency of the proposed numerical method.     

随机微分方程改进的分裂步单支θ方法的强收敛性

本文提出了一个改进的分裂步单支θ方法,在漂移项系数满足单边Lipschitz条件下,证明了当数值方法的参数θ满足1/2 ≤ θ ≤ 1时,该数值方法对于这类随机微分方程是强收敛的,并在现有文献的基础上将方法的收敛阶从1/2阶提高到1阶;当0 ≤ θ ≤ 1/2时,若漂移项系数进一步满足线性增长条件,该数值方法也是强收敛的,收敛阶为1阶.文末的数值试验验证了理论结果的正确性.     

The Convergence of Truncated Euler-Maruyama Method for Stochastic Differential Equations with Piecewise Continuous Arguments Under Generalized One-Sided Lipschitz Condition

In this paper, we consider the stochastic differential equations with piecewise continuous arguments (SDEPCAs) in which the drift coefficient satisfies the generalized one-sided Lipschitz condition and the diffusion coefficient satisfies the linear growth condition. Since the delay term $t-[t]$ of SDEPCAs is not continuous and differentiable, the variable substitution method is not suitable. To overcome this difficulty, we adopt new techniques to prove the boundedness of the exact solution and the numerical solution. It is proved that the truncated Euler-Maruyama method is strongly convergent to SDEPCAs in the sense of $L^{\bar{q}}(\bar{q}\ge 2)$. We obtain the convergence order with some additional conditions. An example is presented to illustrate the analytical theory.     

Convergence of the Tamed Euler Method for Stochastic Differential Equations with Piecewise Continuous Arguments Under Non-global Lipschitz Continuous Coefficients

In this paper, we deal with the strong convergence of numerical methods for stochastic differential equations with piecewise continuous arguments (SEPCAs) with at most polynomially growing drift coefficients and global Lipschitz continuous diffusion coefficients. An explicit and time-saving tamed Euler method is used to solve this type of SEPCAs. We show that the tamed Euler method is bounded in pth moment. And then the convergence of the tamed Euler method is proved. Moreover, the convergence order is one-half. Several numerical simulations are shown to verify the convergence of this method.     

B-convergence of split-step one-leg theta methods for stochastic differential equations

For stochastic differential equations (SDEs) with a superlinearly growing and globally one-sided Lipschitz continuous drift coefficient, the explicit schemes fail to converge strongly to the exact solution (see, Hutzenthaler, Jentzen and Kloeden in Proc. R. Soc. A, rspa.2010.0348v1?Crspa.2010.0348, 2010). In this article a class of implicit methods, called split-step one-leg theta methods (SSOLTM), are introduced and are shown to be mean-square convergent for such SDEs if the method parameter satisfies $\frac{1}{2}\leq\theta \leq1$ . This result gives an extension of B-convergence from the theta method for deterministic ordinary differential equations (ODEs) to SSOLTM for SDEs. Furthermore, the optimal rate of convergence can be recovered if the drift coefficient behaves like a polynomial. Finally, numerical experiments are included to support our assertions.     

随机微分方程分步单支theta方法的稳定性(英文)

本文研究随机微分方程单支theta方法的均方稳定性.首先,对线性检验方程,当0≤θ<1时,分步单支theta方法在一定的步长限制下能保持原系统的均方稳定性,当θ=1时,方法按任意步长都能保持原系统的稳定性.其次,对满足单边Lipschitz条件的非线性随机微分方程,当1/2<θ0<θ<1时,方法能保持原系统的均方指数稳定性,但对步长有限制,如果θ=1,对步长限制消失.     

Asymptotic error distribution for the Euler scheme with locally Lipschitz coefficients

In traditional works on numerical schemes for solving stochastic differential equations (SDEs), the globally Lipschitz assumption is often assumed to ensure different types of convergence. In practice, this is often too strong a condition. Brownian motion driven SDEs used in applications sometimes have coefficients which are only Lipschitz on compact sets, but the paths of the SDE solutions can be arbitrarily large. In this paper, we prove convergence in probability and a weak convergence result under a less restrictive assumption, that is, locally Lipschitz and with no finite time explosion. We prove if a numerical scheme converges in probability uniformly on any compact time set (UCP) with a certain rate under a global Lipschitz condition, then the UCP with the same rate holds when a globally Lipschitz condition is replaced with a locally Lipschitz plus no finite explosion condition. For the Euler scheme, weak convergence of the error process is also established. The main contribution of this paper is the proof of $\sqrt{n}$ weak convergence of the normalized error process and the limit process is also provided. We further study the boundedness of the second moments of the weak limit process and its running supremum under both global Lipschitz and locally Lipschitz conditions.     

Stability analysis of block boundary value methods for neutral pantograph equation

This paper deals with the convergence and stability properties of block boundary value methods (BBVMs) for the neutral pantograph equation. Due to its unbounded time lags and limited computer memory, a change in the independent variable is used to transform a pantograph equation into a non-autonomous differential equation with a constant delay but variable coefficients. It is shown under the classical Lipschitz condition that a BBVM is convergent of order p if the underlying boundary value method is consistent with order p. Furthermore, it is proved under a certain condition that BBVMs can preserve the asymptotic stability of exact solutions for the neutral pantograph equation. Meanwhile, some numerical experiments are given to confirm the main conclusions.     

Almost sure exponential stability of numerical solutions for stochastic delay differential equations with jumps

This paper deals with the almost sure exponential stability of the Euler-type methods for nonlinear stochastic delay differential equations with jumps by using the discrete semimartingale convergence theorem. It is shown that the explicit Euler method reproduces the almost sure exponential stability under an additional linear growth condition. By replacing the linear growth condition with the one-sided Lipschitz condition, the backward Euler method is able to reproduce the stability property.     

Discrete approximations, relaxation, and optimization of one-sided Lipschitzian differential inclusions in Hilbert spaces

We study discrete approximations of nonconvex differential inclusions in Hilbert spaces and dynamic optimization/optimal control problems involving such differential inclusions and their discrete approximations. The underlying feature of the problems under consideration is a modified one-sided Lipschitz condition imposed on the right-hand side (i.e., on the velocity sets) of the differential inclusion, which is a significant improvement of the conventional Lipschitz continuity. Our main attention is paid to establishing efficient conditions that ensure the strong approximation (in the W^1,p-norm as p1) of feasible trajectories for the one-sided Lipschitzian differential inclusions under consideration by those for their discrete approximations and also the strong convergence of optimal solutions to the corresponding dynamic optimization problems under discrete approximations. To proceed with the latter issue, we derive a new extension of the Bogolyubov-type relaxation/density theorem to the case of differential inclusions satisfying the modified one-sided Lipschitzian condition. All the results obtained are new not only in the infinite-dimensional Hilbert space framework but also in finite-dimensional spaces.     

带Poisson跳和Markovian调制的中立型随机延迟微分方程的数值解的收敛性

本文研究带Poisson跳和Markovian调制的中立型随机微分方程的数值解的收敛性质.用数值逼近方法求此微分方程的解,并证明了Euler近似解在此线性增长条件和全局Lipschitz条件更弱的条件下仍均方收敛于此方程的解析解.     

The tamed Milstein method for commutative stochastic differential equations with non-globally Lipschitz continuous coefficients

For stochastic differential equations (SDEs) with a superlinearly growing and globally one-sided Lipschitz continuous drift coefficient, the classical explicit Euler scheme fails to converge strongly to the exact solution. Recently, an explicit strongly convergent numerical scheme, called the tamed Euler method, has been proposed in [8] for such SDEs. Motivated by their work, we here introduce a tamed version of the Milstein scheme for SDEs with commutative noise. The proposed method is also explicit and easily implementable, but achieves higher strong convergence order than the tamed Euler method does. In recovering the strong convergence order one of the new method, new difficulties arise and kind of a bootstrap argument is developed to overcome them. Finally, an illustrative example confirms the computational efficiency of the tamed Milstein method compared to the tamed Euler method.     

Convergence and Stability of the Truncated Euler-Maruyama Method for Stochastic Differential Equations with Piecewise Continuous Arguments

In this paper, we develop the truncated Euler-Maruyama (EM) method for stochastic differential equations with piecewise continuous arguments (SDEPCAs), and consider the strong convergence theory under the local Lipschitz condition plus the Khasminskii-type condition. The order of convergence is obtained. Moreover, we show that the truncated EM method can preserve the exponential mean square stability of SDEPCAs. Numerical examples are provided to support our conclusions.     

On delay-dependent stability for a class of nonlinear stochastic systems with multiple state delays

Global asymptotic stability conditions for nonlinear stochastic systems with multiple state delays are obtained based on the convergence theorem for semimartingale inequalities, without assuming the Lipschitz conditions for nonlinear drift functions. The Lyapunov–Krasovskii and degenerate functionals techniques are used. The derived stability conditions are directly expressed in terms of the system coefficients. Nontrivial examples of nonlinear systems satisfying the obtained stability conditions are given.