Optimal simulation schemes for Lévy driven stochastic differential equations

We consider a general class of high order weak approximation schemes for stochastic differential equations driven by L\'evy processes with infinite activity. These schemes combine a compound Poisson approximation for the jump part of the L\'evy process with a high order scheme for the Brownian driven component, applied between the jump times. The overall approximation is analyzed using a stochastic splitting argument. The resulting error bound involves separate contributions of the compound Poisson approximation and of the discretization scheme for the Brownian part, and allows, on one hand, to balance the two contributions in order to minimize the computational time, and on the other hand, to study the optimal design of the approximating compound Poisson process. For driving processes whose L\'evy measure explodes near zero in a regularly varying way, this procedure allows to construct discretization schemes with arbitrary order of convergence.

[1]  Arturo Kohatsu-Higa,et al.  AN OPERATOR APPROACH FOR MARKOV CHAIN WEAK APPROXIMATIONS WITH AN APPLICATION TO INFINITE ACTIVITY L , 2009, 0908.1021.

[2]  D. Applebaum.  Lévy Processes and Stochastic Calculus: Preface , 2009 .

[3]  David Applebaum,et al.  Lévy Processes and Stochastic Calculus by David Applebaum , 2009 .

[4]  Terry Lyons,et al.  Efficient and Practical Implementations of Cubature on Wiener Space , 2011 .

[5]  Ernesto Mordecki,et al.  Adaptive Weak Approximation of Diffusions with Jumps , 2008, SIAM J. Numer. Anal..

[6]  M. Kreĭn,et al.  The Markov Moment Problem and Extremal Problems , 1977 .

[7]  P. Tankov.  High Order Weak Approximation Schemes for Lévy-Driven SDEs , 2010, 1012.5806.

[8]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[9]  Terry Lyons,et al.  Cubature on Wiener space , 2004, Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[10]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[11]  Philip Protter,et al.  The Euler scheme for Lévy driven stochastic differential equations , 1997 .

[12]  Jean Jacod,et al.  The approximate Euler method for Lévy driven stochastic differential equations , 2005 .

[13]  Sylvain Rubenthaler,et al.  Numerical simulation of the solution of a stochastic differential equation driven by a Lévy process , 2003 .

[14]  Josef Teichmann,et al.  A new extrapolation method for weak approximation schemes with applications , 2009, 0911.4380.

[15]  A. Kohatsu-Higa,et al.  Jump-adapted discretization schemes for Lévy-driven SDEs , 2010 .

[16]  S. Ninomiya,et al.  Weak Approximation of Stochastic Differential Equations and Application to Derivative Pricing , 2006, math/0605361.

[17]  Steffen Dereich,et al.  Multilevel Monte Carlo algorithms for L\'{e}vy-driven SDEs with Gaussian correction , 2011, 1101.1369.

[18]  P. Protter.  Stochastic integration and differential equations , 1990 .

[19]  Peter Tankov,et al.  Monte Carlo Option Pricing for Tempered Stable (CGMY) Processes , 2007 .

[20]  M. Meerschaert.  Regular Variation in R k , 1988 .

[21]  Nicola Bruti-Liberati,et al.  Strong approximations of stochastic differential equations with jumps , 2007 .

[22]  S. Asmussen,et al.  Approximations of small jumps of Lévy processes with a view towards simulation , 2001, Journal of Applied Probability.

[23]  J. K. Hunter,et al.  Measure Theory , 2007 .

[24]  P. Protter.  Stochastic integration and differential equations : a new approach , 1990 .