Runge–Kutta method (SDE)

From Wikipedia, the free encyclopedia

In mathematics, the Runge - Kutta method is a technique for the approximate numerical solution of a stochastic differential equation. It is a generalization of the Runge-Kutta method for ordinary differential equations to stochastic differential equations.

Consider the Itō diffusion X satisfying the following Itō stochastic differential equation

{\mathrm{d} X_{t}} = a(X_{t}) \, \mathrm{d} t + b(X_{t}) \, \mathrm{d} W_{t},

with initial condition X0 = x0, where Wt stands for the Wiener process, and suppose that we wish to solve this SDE on some interval of time [0, T]. Then the Runge-Kutta approximation to the true solution X is the Markov chain Y defined as follows:

  • partition the interval [0, T] into N equal subintervals of width δ = T ⁄ N > 0:
0 = \tau_{0} < \tau_{1} < \dots < \tau_{N} = T;
  • set Y0 = x0;
  • recursively define Yn for 1 ≤ n ≤ N by
Y_{n + 1} = Y_{n} + a(Y_{n}) \delta + b(Y_{n}) \Delta W_{n} + \frac{1}{2} \left( b(\hat{\Upsilon}_{n}) - b(Y_{n}) \right) \left( (\Delta W_{n})^{2} - \delta \right) \delta^{-1/2},
where
\Delta W_{n} = W_{\tau_{n + 1}} - W_{\tau_{n}}
and
\hat{\Upsilon}_{n} = Y_{n} + a(Y_n) \delta + b(Y_{n}) \delta^{1/2}.

Note that the random variables ΔWn are independent and identically distributed normal random variables with expected value zero and variance δ.

This scheme has strong order 1, meaning that the approximation error of the actual solution at a fixed time scales with the time step δ. It has also weak order 1, meaning that the error on the statistics of the solution scales with the time step δ. See the references for complete and exact statements.

The functions a and b can be time-varying without any complication. The method can be generalized to the case of several coupled equations; the principle is the same but the equations become longer. Higher-order schemes also exist, but become increasingly complex.

[edit] References

  • Kloeden, P.E., & Platen, E. (1999). Numerical Solution of Stochastic Differential Equations. Berlin: Springer. ISBN 3-540-54062-8.