Abstract:We address the problem of approximating the moments of the solution, $\boldsymbol{X}(t)$, of an Itô stochastic differential equation (SDE) with drift and a diffusion terms over a time-grid $t_0, t_1, \ldots, t_n$. In particular, we assume an explicit numerical scheme for the generation of sample paths $\hat{\boldsymbol{X}}(t_0), \ldots, \hat{\boldsymbol{X}}(t_n), \ldots$ and then obtain recursive equations that yield any desired non-central moment of $\hat{\boldsymbol{X}}(t_n)$ as a function of the initial condition $\boldsymbol{X}_0$. The core of the methodology is the decomposition of the numerical solution into a "central part" and an "effective noise" term. The central term is computed deterministically from the ordinary differential equation (ODE) that results from eliminating the diffusion term in the SDE, while the effective noise accounts for the stochastic deviation from the numerical solution of the ODE. For simplicity, we describe algorithms based on an Euler-Maruyama integrator, but other explicit numerical schemes can be exploited in the same way. We also apply the moment approximations to construct estimates of the 1-dimensional marginal probability density functions of $\hat{\boldsymbol{X}}(t_n)$ based on a Gram-Charlier expansion. Both for the approximation of moments and 1-dimensional densities, we describe how to handle the cases in which the initial condition is fixed (i.e., $\boldsymbol{X}_0 = \boldsymbol{x}_0$ for some known $\boldsymbol{x_0}$) or random. In the latter case, we resort to polynomial chaos expansion (PCE) schemes to approximate the target moments. The methodology has been inspired by the PCE and differential algebra (DA) methods used for uncertainty propagation in astrodynamics problems. Hence, we illustrate its application for the quantification of uncertainty in a 2-dimensional Keplerian orbit perturbed by a Wiener noise process.

Path Integral Methods for Stochastic Differential Equations

Master equations and the theory of stochastic path integrals

A Primer on Stochastic Partial Differential Equations with Spatially Correlated Noise

Numerical Solution of Stochastic Differential Equations: Diffusion and Jump-Diffusion Processes

Statistical inference for stochastic differential equations

Numerical methods for the deterministic second moment equation of parabolic stochastic PDEs

Construction of stochastic hybrid path integrals using "quantum-mechanical'' operators

A survey of numerical methods for stochastic differential equations

Coupling Methods and Applications on Path Dependent McKean-Vlasov SDEs

A brief review on stability investigations of numerical methods for systems of stochastic differential equations

Polynomial propagation of moments in stochastic differential equations

Path Integral Methods and Applications

Path integral Monte Carlo ground state approach: formalism, implementation, and applications

Wellposedness, exponential ergodicity and numerical approximation of fully super-linear McKean--Vlasov SDEs and associated particle systems

Stochastic solutions and singular partial differential equations

The physical mechanism of stochastic calculus in random walks

Direct Numerical Solutions to Stochastic Differential Equations with Multiplicative Noise

Semi-linear Backward Stochastic Integral Partial Differential Equations Driven by a Brownian Motion and a Poisson Point Process

Field theories and quantum methods for stochastic reaction-diffusion systems

Stochastic methods for solving high-dimensional partial differential equations

Random walk algorithm for the Dirichlet problem for parabolic integro-differential equation