Mean Square Polynomial Stability of Numerical Solutions to a Class of Stochastic Differential Equations

The exponential stability of numerical methods to stochastic differential equations (SDEs) has been widely studied. In contrast, there are relatively few works on polynomial stability of numerical methods. In this letter, we address the question of reproducing the polynomial decay of a class of SDEs using the Euler--Maruyama method and the backward Euler--Maruyama method. The key technical contribution is based on various estimates involving the gamma function

Mean-square convergence and stability of the backward Euler method for stochastic differential delay equations with highly nonlinear growing coefficients

Over the last few decades, the numerical methods for stochastic differential delay equations (SDDEs) have been investigated and developed by many scholars. Nevertheless, there is still little work to be completed. By virtue of the novel technique, this paper focuses on the mean-square convergence and stability of the backward Euler method (BEM) for SDDEs whose drift and diffusion coefficients can both grow polynomially. The upper mean-square error bounds of BEM are obtained. Then the convergence rate, which is one-half, is revealed without using the moment boundedness of numerical solutions. Furthermore, under fairly general conditions, the novel technique is applied to prove that the BEM can inherit the exponential mean-square stability with a simple proof. At last, two numerical experiments are implemented to illustrate the reliability of the theories

Stochastic ordinary differential equations in applied and computational mathematics

Using concrete examples, we discuss the current and potential use of stochastic ordinary differential equations (SDEs) from the perspective of applied and computational mathematics. Assuming only a minimal background knowledge in probability and stochastic processes, we focus on aspects that distinguish SDEs from their deterministic counterparts. To illustrate a multiscale modelling framework, we explain how SDEs arise naturally as diffusion limits in the type of discrete-valued stochastic models used in chemical kinetics, population dynamics, and, most topically, systems biology. We outline some key issues in existence, uniqueness and stability that arise when SDEs are used as physical models, and point out possible pitfalls. We also discuss the use of numerical methods to simulate trajectories of an SDE and explain how both weak and strong convergence properties are relevant for highly-efficient multilevel Monte Carlo simulations. We flag up what we believe to be key topics for future research, focussing especially on nonlinear models, parameter estimation, model comparison and multiscale simulation

Stability of numerical method for semi-linear stochastic pantograph differential equations

Abstract As a particular expression of stochastic delay differential equations, stochastic pantograph differential equations have been widely used in nonlinear dynamics, quantum mechanics, and electrodynamics. In this paper, we mainly study the stability of analytical solutions and numerical solutions of semi-linear stochastic pantograph differential equations. Some suitable conditions for the mean-square stability of an analytical solution are obtained. Then we proved the general mean-square stability of the exponential Euler method for a numerical solution of semi-linear stochastic pantograph differential equations, that is, if an analytical solution is stable, then the exponential Euler method applied to the system is mean-square stable for arbitrary step-size h > 0 $h>0$ . Numerical examples further illustrate the obtained theoretical results

Almost sure exponential stability of the Euler–Maruyama approximations for stochastic functional differential equations

By the continuous and discrete nonnegative semimartingale convergence theorems, this paper investigates conditions under which the Euler–Maruyama (EM) approximations of stochastic functional differential equations (SFDEs) can share the almost sure exponential stability of the exact solution. Moreover, for sufficiently small stepsize, the decay rate as measured by the Lyapunov exponent can be reproduced arbitrarily accurately