An optimal polynomial approximation of Brownian motion

Abstract

In this paper, we will present a strong (or pathwise) approximation of standard Brownian motion by a class of orthogonal polynomials. The coefficients that are obtained from the expansion of Brownian motion in this polynomial basis are independent Gaussian random variables. Therefore it is practical (requires $N$ independent Gaussian coefficients) to generate an approximate sample path of Brownian motion that respects integration of polynomials with degree less than $N$. Moreover, since these orthogonal polynomials appear naturally as eigenfunctions of an integral operator defined by the Brownian bridge covariance function, the proposed approximation is optimal in a certain weighted $L^2(\mathbb{P})$ sense. In addition, discretizing Brownian paths as piecewise parabolas gives a locally higher order numerical method for stochastic differential equations (SDEs) when compared to the standard piecewise linear approach. We shall demonstrate these ideas by simulating Inhomogeneous Geometric Brownian Motion (IGBM). This numerical example will also illustrate the deficiencies of the piecewise parabola approximation when compared to a new version of the asymptotically efficient log-ODE (or Castell-Gaines) method.

arXiv comment27 pages, 8 figures

References in corpus (1)

• A Multiscale Guide to Brownian Motion

Cited by in corpus (8)

• First order convergence of Milstein schemes for McKean-Vlasov equations and interacting particle systems
• Brownian bridge expansions for Lévy area approximations and particular values of the Riemann zeta function
• An asymptotic radius of convergence for the Loewner equation and simulation of $SL{E}_k$ traces via splitting
• High order splitting methods for SDEs satisfying a commutativity condition
• A semicircle law and decorrelation phenomena for iterated Kolmogorov loops
• Modelling pathwise uncertainty of Stochastic Differential Equations samplers via Probabilistic Numerics
• Generative Modelling of Lévy Area for High Order SDE Simulation
• Simple third order operator-splitting schemes for stochastic mechanics and field theory