A mixed-step algorithm for the approximation of the stationary regime of a diffusion

STOCHASTIC PROCESSES AND THEIR APPLICATIONS

卷:124

Article

Pages, Gilles¹  Panloup, Fabien^2,3 

[1] UPMC, UMR 7599, Lab Probabil & Modeles Aleatoires, F-75252 Paris 5, France

[2] Univ Toulouse 3, Inst Math Toulouse, F-31077 Toulouse 4, France

[3] INSA Toulouse, F-31077 Toulouse 4, France

关键词: Stochastic differential equation;    Stationary process;    Steady regime;    Ergodic diffusion;    Central Limit Theorem;    Euler scheme;    Poisson equation;    Richardson Romberg extrapolation;   

DOI  :  10.1016/j.spa.2013.07.011

来源: Elsevier

【 摘 要 】

In some recent papers, some procedures based on some weighted empirical measures related to decreasing-step Euler schemes have been investigated to approximate the stationary regime of a diffusion (possibly with jumps) for a class of functionals of the process. This method is efficient but needs the computation of the function at each step. To reduce the complexity of the procedure (especially for functionals), we propose in this paper to study a new scheme, called the mixed-step scheme, where we only keep some regularly time-spaced values of the Euler scheme. Our main result is that, when the coefficients of the diffusion are smooth enough, this alternative does not change the order of the rate of convergence of the procedure. We also investigate a Richardson-Romberg method to speed up the convergence and show that the variance of the original algorithm can be preserved under a uniqueness assumption for the invariant distribution of the duplicated diffusion, condition which is extensively discussed in the paper. Finally, we conclude by giving sufficient asymptotic confluence conditions for the existence of a smooth solution to a discrete version of the associated Poisson equation, condition which is required to ensure the rate of convergence results. (C) 2013 Published by Elsevier B.V.

【 授权许可】

Free