The Azimuth Project
Milstein scheme

Skip the Navigation Links | Home Page | All Pages | Recently Revised | Authors | Feeds |

Contents

Idea

The Milstein scheme, named after Grigori N. Milstein, is a numerical discrete approximation scheme for stochastic differential equations.

There is an explicit and an implicit version of the Milstein scheme.

The Milstein scheme is the simplest scheme that achieves a higher strong order of convergence than the Euler scheme, namely 1.0.

Details

Let

dX=adt+bdW d X = a \; d t + b \; d W

be an Itô stochastic differential equation with a,ba, b suitable functions. In one dimension, the Milstein scheme is the following approximation scheme:

Y n+1=Y n+αa(τ n+1,Y n+1)(τ n+1−τ n)+(1−α)a(τ n,Y n)(τ n+1−τ n)+b(W τ n+1−W τ n)+12bb′((W τ n+1−W τ n) 2−(τ n+1−τ n)) Y_{n+1} = Y_n + \alpha \; a(\tau_{n+1}, Y_{n+1}) (\tau_{n+1} - \tau_{n}) + (1 - \alpha) a(\tau_{n}, Y_{n}) (\tau_{n+1} - \tau_{n}) + b (W_{\tau_{n+1}} - W_{\tau_n}) + \frac{1}{2} b b' ( (W_{\tau_{n+1}} - W_{\tau_n} )^2 - (\tau_{n+1} - \tau_{n}))

The number α\alpha is chosen to be between 0 and 1.

References

category: computational methods

Revised on May 15, 2012 16:37:20 by Tim van Beek (0.0.0.0)
View Source | Back in time (7 revisions) | See changes | History | Views: Print | TeX | Source | Linked from: Stochastic differential equation, Stochastic delay differential equation