Exact numerical simulation of the Ornstein-Uhlenbeck process and its integral

Daniel T. Gillespie
Phys. Rev. E 54, 2084 – Published 1 August 1996
PDFExport Citation

Abstract

A numerical simulation algorithm that is exact for any time step Δt>0 is derived for the Ornstein-Uhlenbeck process X(t) and its time integral Y(t). The algorithm allows one to make efficient, unapproximated simulations of, for instance, the velocity and position components of a particle undergoing Brownian motion, and the electric current and transported charge in a simple R-L circuit, provided appropriate values are assigned to the Ornstein-Uhlenbeck relaxation time τ and diffusion constant c. A simple Taylor expansion in Δt of the exact simulation formulas shows how the first-order simulation formulas, which are implicit in the Langevin equation for X(t) and the defining equation for Y(t), are modified in second order. The exact simulation algorithm is used here to illustrate the zero-τ limit theorem. © 1996 The American Physical Society.

  • Received 27 November 1995

DOI:https://doi.org/10.1103/PhysRevE.54.2084

©1996 American Physical Society

Authors & Affiliations

Daniel T. Gillespie

  • Research and Technology Division, Naval Air Warfare Center, China Lake, California 93555

References (Subscription Required)

Click to Expand
Issue

Vol. 54, Iss. 2 — August 1996

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×
×

Images

×

Sign up to receive regular email alerts from Physical Review E

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×