Polynomial propagation of moments in stochastic differential equations

by   Albert López-Yela, et al.

We address the problem of approximating the moments of the solution, X(t), of an Itô stochastic differential equation (SDE) with drift and a diffusion terms over a time-grid t_0, t_1, …, t_n. In particular, we assume an explicit numerical scheme for the generation of sample paths X̂(t_0), …, X̂(t_n), … and then obtain recursive equations that yield any desired non-central moment of X̂(t_n) as a function of the initial condition 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 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., X_0 = x_0 for some known 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.


page 1

page 2

page 3

page 4


A Strongly Monotonic Polygonal Euler Scheme

Rate of convergence results are presented for a new class of explicit Eu...

Least squares estimators based on the Adams method for stochastic differential equations with small Lévy noise

We consider stochastic differential equations (SDEs) driven by small Lév...

Improving the approximation of the first and second order statistics of the response process to the random Legendre differential equation

In this paper, we deal with uncertainty quantification for the random Le...

On mathematical aspects of evolution of dislocation density in metallic materials

This paper deals with the solution of delay differential equations descr...

Non-linear Gaussian smoothing with Taylor moment expansion

This letter is concerned with solving continuous-discrete Gaussian smoot...

Conditional Karhunen-Loève expansion for uncertainty quantification and active learning in partial differential equation models

We use a conditional Karhunen-Loève (KL) model to quantify and reduce un...

Please sign up or login with your details

Forgot password? Click here to reset