 Ornstein-Uhlenbeck process

An Ornstein-Uhlenbeck process combines mean reversion with continuing random fluctuations. In the scalar model, its stochastic differential equation iswhere $W$ is a standard Brownian motion. The drift points towards $\mu$: it is negative above $\mu$ and positive below it. The coefficient $\theta$ sets the speed of mean reversion, while $\sigma$ sets the strength of the random forcing. We consider a fixed or normal initial value independent of future Brownian increments, giving a Gaussian process and a Markov process.

For a stipulated example, use $\theta =1$, $\mu =2$, $\sigma =\sqrt{2}$ and the fixed initial value $X_0=6$, in arbitrary consistent state and time units. Then the expected value relaxes as $2+4e^{-t}$ and the variance grows as $1-e^{-2t}$. At time one they are about $3.4715$ and $0.8647$; the limiting variance is one. The sample paths keep fluctuating even after the expected value is close to two. The Ornstein-Uhlenbeck transition distribution describes this evolving uncertainty; Ornstein-Uhlenbeck stationary initialization instead starts with that limiting normal distribution already in place.

A physical motivation is the velocity of a particle undergoing Brownian motion: linear drag removes a persistent initial velocity, while random forcing continually renews it. The original calculation gives a normal distribution for the velocity after subtracting its exponentially decaying initial contribution, with decay rate equal to the friction coefficient divided by the particle's mass.^[ref] A position coordinate in an overdamped harmonic optical trap gives another interpretation of the same scalar equation, with $\theta =k/\zeta$ for spring constant $k$ and drag coefficient $\zeta$. Its physical assumptions and connection to measurements are developed in thermal calibration of an optical trap.