Picture of athlete cycling

Open Access research with a real impact on health...

The Strathprints institutional repository is a digital archive of University of Strathclyde's Open Access research outputs. Strathprints provides access to thousands of Open Access research papers by Strathclyde researchers, including by researchers from the Physical Activity for Health Group based within the School of Psychological Sciences & Health. Research here seeks to better understand how and why physical activity improves health, gain a better understanding of the amount, intensity, and type of physical activity needed for health benefits, and evaluate the effect of interventions to promote physical activity.

Explore open research content by Physical Activity for Health...

Almost sure and moment exponential stability in the numerical simulation of stochastic differential equations

Higham, D.J. and Mao, X. and Yuan, C. (2007) Almost sure and moment exponential stability in the numerical simulation of stochastic differential equations. SIAM Journal on Numerical Analysis, 45 (2). pp. 592-609. ISSN 0036-1429

HighamP70.pdf - Final Published Version

Download (192kB) | Preview


Relatively little is known about the ability of numerical methods for stochastic differential equations (SDEs) to reproduce almost sure and small-moment stability. Here, we focus on these stability properties in the limit as the timestep tends to zero. Our analysis is motivated by an example of an exponentially almost surely stable nonlinear SDE for which the Euler-Maruyama (EM)method fails to reproduce this behavior for any nonzero timestep. We begin by showing that EM correctly reproduces almost sure and small-moment exponential stability for sufficiently small timesteps on scalar linear SDEs. We then generalize our results to multidimensional nonlinear SDEs. We show that when the SDE obeys a linear growth condition, EM recovers almost surely exponential stability very well. Under the less restrictive condition that the drift coefficient of the SDE obeys a one-sided Lipschitz condition, where EM may break down, we show that the backward Euler method maintains almost surely exponential stability.