A generic construction for high order approximation schemes of semigroups using random grids
Abstract
Our aim is to construct high order approximation schemes for general semigroups of linear operators Pt,t≥ 0. In order to do it, we fix a time horizon T and the discretization steps hl=Tnl,l∈ N and we suppose that we have at hand some short time approximation operators Ql such that Phl=Ql+O(hl1+α) for some α>0. Then, we consider random time grids Π(ω)=\t0(ω)=0<t1(ω)<...<tm(ω)=T\ such that for all 1 k m, tk(ω)-tk-1(ω)=hlk for some lk∈ N, and we associate the approximation discrete semigroup PTΠ(ω)=Qln...Ql1. Our main result is the following: for any approximation order ν, we can construct random grids Πi(ω) and coefficients ci, with i=1,...,r such that \[ Ptf=Σi=1rciE(PtΠi(ω)f(x))+O(n-ν) \]% with the expectation concerning the random grids Πi(ω). Besides, Card(Πi(ω))=O(n) and the complexity of the algorithm is of order n, for any order of approximation ν. The standard example concerns diffusion processes, using the Euler approximation for~Ql. In this particular case and under suitable conditions, we are able to gather the terms in order to produce an estimator of Ptf with finite variance. However, an important feature of our approach is its universality in the sense that it works for every general semigroup Pt and approximations. Besides, approximation schemes sharing the same α lead to the same random grids Πi and coefficients ci. Numerical illustrations are given for ordinary differential equations, piecewise deterministic Markov processes and diffusions.
Turn this paper into a lesson
ArcXiv compiles a structured reading guide from this paper's metadata: plain-English importance, contributions, prerequisite concepts, which sections to read first, flashcards, and a quiz. Grounded in the abstract, never invented.