High order weak methods for stochastic differential equations based on modified equations

Inspired by recent advances in the theory of modified differential equations, we propose a new methodology for constructing numerical integrators with high weak order for the time integration of stochastic differential equations. This approach is illustrated with the constructions of new high order weak methods, in particular, implicit integrators well suited for stiff stochastic problems, and integrators that exactly conserve all quadratic first integrals of a stochastic dynamical system. Numerical examples confirm the theoretical results and show the versatility of the methodology.

[1]  G. Maruyama.  Continuous Markov processes and stochastic equations , 1955 .

[2]  G. Mil’shtein.  A Method of Second-Order Accuracy Integration of Stochastic Differential Equations , 1979 .

[3]  B. Øksendal.  Stochastic Differential Equations , 1985 .

[4]  G. Mil’shtein.  Weak Approximation of Solutions of Systems of Stochastic Differential Equations , 1986 .

[5]  D. Talay,et al.  Expansion of the global error for numerical schemes solving stochastic differential equations , 1990 .

[6]  E. Platen,et al.  Option Pricing under Incompleteness and Stochastic Volatility , 1992 .

[7]  E. Platen,et al.  Higher-Order Weak Approximation of Ito Diffusions by Markov Chains , 1992, Probability in the Engineering and Informational Sciences.

[8]  L. Arnold.  Stochastic Differential Equations: Theory and Applications , 1992 .

[9]  W. P. Petersen.  Lagged Fibonacci Series Random Number Generators for the NEC SX-3 , 1994, Int. J. High Speed Comput..

[10]  Vigirdas Mackevičius,et al.  Second-order weak approximations for stratonovich stochastic differential equations , 1994 .

[11]  Peter E. Kloeden,et al.  Extrapolation methods for the weak approximation of Ito diffusions , 1995 .

[12]  Yoshihiro Saito,et al.  Stability Analysis of Numerical Schemes for Stochastic Differential Equations , 1996 .

[13]  Ming-Chen Liao,et al.  Random Motion of a Rigid Body , 1997 .

[14]  Kevin Burrage,et al.  General order conditions for stochastic Runge-Kutta methods for both commuting and non-commuting stochastic ordinary differential equation systems , 1998 .

[15]  Desmond J. Higham.  A-STABILITY AND STOCHASTIC MEAN-SQUARE STABILITY , 2000 .

[16]  Desmond J. Higham,et al.  Mean-Square and Asymptotic Stability of the Stochastic Theta Method , 2000, SIAM J. Numer. Anal..

[17]  G. N. Milstein,et al.  Numerical Methods for Stochastic Systems Preserving Symplectic Structure , 2002, SIAM J. Numer. Anal..

[18]  Desmond J. Higham,et al.  MAPLE and MATLAB for stochastic differential equations in finance , 2002 .

[19]  Jesús Vigo-Aguiar,et al.  Weak Second Order Conditions for Stochastic Runge-Kutta Methods , 2002, SIAM J. Sci. Comput..

[20]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[21]  M. V. Tretyakov,et al.  Stochastic Numerics for Mathematical Physics , 2004, Scientific Computation.

[22]  Tony Shardlow,et al.  Modified Equations for Stochastic Differential Equations , 2006 .

[23]  Ernst Hairer,et al.  Numerical integrators based on modified differential equations , 2007, Math. Comput..

[24]  J. Ortega,et al.  REDUCTION, RECONSTRUCTION, AND SKEW-PRODUCT DECOMPOSITION OF SYMMETRIC STOCHASTIC DIFFERENTIAL EQUATIONS , 2007, 0705.3156.

[25]  Grigorios A. Pavliotis,et al.  Multiscale Methods: Averaging and Homogenization , 2008 .

[26]  Michael B. Giles,et al.  Multilevel Monte Carlo Path Simulation , 2008, Oper. Res..

[27]  Kristian Debrabant,et al.  Classification of stochastic Runge-Kutta methods for the weak approximation of stochastic differential equations , 2008, Math. Comput. Simul..

[28]  Assyr Abdulle,et al.  S-ROCK: Chebyshev Methods for Stiff Stochastic Differential Equations , 2008, SIAM J. Sci. Comput..

[29]  A. Abdulle,et al.  S-ROCK methods for stiff Ito SDEs , 2008 .

[30]  Andreas Rößler.  Second Order Runge-Kutta Methods for Itô Stochastic Differential Equations , 2009, SIAM J. Numer. Anal..

[31]  Cónall Kelly,et al.  Towards a Systematic Linear Stability Analysis of Numerical Methods for Systems of Stochastic Differential Equations , 2010, SIAM J. Numer. Anal..

[32]  Kevin Burrage,et al.  Explicit Stochastic Runge‐Kutta Methods with Large Stability Regions , 2010 .

[33]  B. Rozovskii,et al.  The Oxford Handbook of Nonlinear Filtering , 2011 .

[34]  K. C. Zygalakis,et al.  On the Existence and the Applications of Modified Equations for Stochastic Differential Equations , 2011, SIAM J. Sci. Comput..

[35]  Pol D. Spanos,et al.  Computational Stochastic Mechanics , 2011 .

[36]  Erwan Faou,et al.  Weak Backward Error Analysis for SDEs , 2012, SIAM J. Numer. Anal..

[37]  Assyr Abdulle,et al.  Second weak order explicit stabilized methods for stiff stochastic differential equations , 2012 .