Skip to content
Oakfield Operator Calculus Function Reference Site

Fractional Ornstein–Uhlenbeck Noise

Fractional Ornstein–Uhlenbeck (fOU) processes extend the classic OU process by driving it with fractional Brownian motion, blending mean reversion with persistent or antipersistent correlations.

In the mean-zero case, a common (mild) form is

Xt=X0eλt+σ0teλ(ts)dBH(s)X_t = X_0 e^{-\lambda t} + \sigma\int_0^t e^{-\lambda(t-s)}\,dB_H(s)

where the stochastic integral is interpreted in a manner appropriate to HH (e.g. Itô for H=12H=\tfrac12, and pathwise/Young-type for H>12H>\tfrac12).

Defined over continuous or discretized time; outputs are real-valued stochastic trajectories.


τrelax=λ1\tau_{\text{relax}} = \lambda^{-1}

Mean reversion rate λ\lambda sets decay to equilibrium; the Hurst parameter HH shapes memory, with H>1/2H>1/2 yielding persistence and H<1/2H<1/2 yielding anti-persistence (and smaller HH producing rougher sample paths).


H=12H = \tfrac{1}{2} recovers the classical Ornstein–Uhlenbeck process (in which case B1/2B_{1/2} is standard Brownian motion and the SDE is interpreted in the Itô sense). Formally taking λ0\lambda \to 0 yields XtX0+σBH(t)X_t \to X_0 + \sigma B_H(t); large λ\lambda forces rapid return toward the mean-reverting equilibrium.


  • H=12H=\tfrac12 recovers the classical OU process.
  • λ0\lambda\to 0 approaches fractional Brownian motion driving.
  • Large λ\lambda yields fast relaxation with reduced variance.

The standard Ornstein–Uhlenbeck process is the H=12H=\tfrac{1}{2} case; fractional Brownian motion is the λ=0\lambda=0 driver; Gaussian white noise arises as the derivative of Brownian motion and underpins the H=12H=\tfrac{1}{2} forcing.


Oakfield uses “fractional OU” behavior in its colored-noise operator:

  • stochastic_noise operator generalizes OU updates with a spectral exponent parameter (alpha) to produce colored noise (OU-like when alpha≈1, with pink/blue tilts as alpha varies).
  • This operator is used as a reusable forcing source in operator graphs rather than being tied to a single integrator.

Replacing Brownian motion with fractional Brownian motion is a standard way to introduce long-memory correlations into otherwise Markovian relaxation dynamics.

Fractional OU models are used when exponential correlation is too short-ranged and algebraic memory is required.


  • Beran, Statistics for Long-Memory Processes
  • Samorodnitsky & Taqqu, Stable Non-Gaussian Random Processes (context on long-memory modeling)