Ornstein–Uhlenbeck process

(Redirected from Ornstein-Uhlenbeck process)

In mathematics, the Ornstein–Uhlenbeck process is a stochastic process with applications in financial mathematics and the physical sciences. Its original application in physics was as a model for the velocity of a massive Brownian particle under the influence of friction. It is named after Leonard Ornstein and George Eugene Uhlenbeck.

Five simulations with θ = 1, σ = 1 and μ = 0.
A 3D simulation with θ = 1, σ = 3, μ = (0, 0, 0) and the initial position (10, 10, 10).

The Ornstein–Uhlenbeck process is a stationary Gauss–Markov process, which means that it is a Gaussian process, a Markov process, and is temporally homogeneous. In fact, it is the only nontrivial process that satisfies these three conditions, up to allowing linear transformations of the space and time variables.[1] Over time, the process tends to drift towards its mean function: such a process is called mean-reverting.

The process can be considered to be a modification of the random walk in continuous time, or Wiener process, in which the properties of the process have been changed so that there is a tendency of the walk to move back towards a central location, with a greater attraction when the process is further away from the center. The Ornstein–Uhlenbeck process can also be considered as the continuous-time analogue of the discrete-time AR(1) process.

Definition

edit
 
Simplified formula for the Ornstein–Uhlenbeck process from the mural shown below.
 
Dutch artist collective De Strakke Hand: Leonard Ornstein mural, showing Ornstein as a cofounder of the Dutch Physical Society (Netherlands Physical Society) at his desk in 1921, and illustrating twice the random walk of a drunkard with a simplified formula for the Ornstein–Uhlenbeck process. Oosterkade, Utrecht, The Netherlands, not far from Ornstein's laboratory. Translated text: Prof. Ornstein researches random motion 1930.

The Ornstein–Uhlenbeck process   is defined by the following stochastic differential equation:

 

where   and   are parameters and   denotes the Wiener process.[2][3][4]

An additional drift term is sometimes added:

 

where   is a constant. The Ornstein–Uhlenbeck process is sometimes also written as a Langevin equation of the form

 

where  , also known as white noise, stands in for the supposed derivative   of the Wiener process.[5] However,   does not exist because the Wiener process is nowhere differentiable,[6] and so the Langevin equation only makes sense if interpreted in distributional sense. In physics and engineering disciplines, it is a common representation for the Ornstein–Uhlenbeck process and similar stochastic differential equations by tacitly assuming that the noise term is a derivative of a differentiable (e.g. Fourier) interpolation of the Wiener process.

Fokker–Planck equation representation

edit

The Ornstein–Uhlenbeck process can also be described in terms of a probability density function,  , which specifies the probability of finding the process in the state   at time  .[5] This function satisfies the Fokker–Planck equation

 

where  . This is a linear parabolic partial differential equation which can be solved by a variety of techniques. The transition probability, also known as the Green's function,   is a Gaussian with mean   and variance  :

 

This gives the probability of the state   occurring at time   given initial state   at time  . Equivalently,   is the solution of the Fokker–Planck equation with initial condition  .

Mathematical properties

edit

Conditioned on a particular value of  , the mean is

 

and the covariance is

 

For the stationary (unconditioned) process, the mean of   is  , and the covariance of   and   is  .

The Ornstein–Uhlenbeck process is an example of a Gaussian process that has a bounded variance and admits a stationary probability distribution, in contrast to the Wiener process; the difference between the two is in their "drift" term. For the Wiener process the drift term is constant, whereas for the Ornstein–Uhlenbeck process it is dependent on the current value of the process: if the current value of the process is less than the (long-term) mean, the drift will be positive; if the current value of the process is greater than the (long-term) mean, the drift will be negative. In other words, the mean acts as an equilibrium level for the process. This gives the process its informative name, "mean-reverting."

Properties of sample paths

edit

A temporally homogeneous Ornstein–Uhlenbeck process can be represented as a scaled, time-transformed Wiener process:

 

where   is the standard Wiener process. This is roughly Theorem 1.2 in Doob 1942. Equivalently, with the change of variable   this becomes

 

Using this mapping, one can translate known properties of   into corresponding statements for  . For instance, the law of the iterated logarithm for   becomes[1]

 

Formal solution

edit

The stochastic differential equation for   can be formally solved by variation of parameters.[7] Writing

 

we get

 

Integrating from   to   we get

 

whereupon we see

 

From this representation, the first moment (i.e. the mean) is shown to be

 

assuming   is constant. Moreover, the Itō isometry can be used to calculate the covariance function by

 

Since the Itô integral of deterministic integrand is normally distributed, it follows that[citation needed]

 

Kolmogorov equations

edit

The infinitesimal generator of the process is[8] If we let  , then the eigenvalue equation simplifies to:  which is the defining equation for Hermite polynomials. Its solutions are  , with  , which implies that the mean first passage time for a particle to hit a point on the boundary is on the order of  .

Numerical simulation

edit

By using discretely sampled data at time intervals of width  , the maximum likelihood estimators for the parameters of the Ornstein–Uhlenbeck process are asymptotically normal to their true values.[9] More precisely,[failed verification] 

 
Four sample paths of different OU-processes with θ = 1, σ =  :
blue: initial value a = 10, μ = 0
orange: initial value a = 0, μ = 0
green: initial value a = −10, μ = 0
red: initial value a = 0, μ = −10

To simulate an OU process numerically with standard deviation   and correlation time  , one method is to apply the finite-difference formula

  where   is a normally distributed random number with zero mean and unit variance, sampled independently at every time-step  .[10]

Scaling limit interpretation

edit

The Ornstein–Uhlenbeck process can be interpreted as a scaling limit of a discrete process, in the same way that Brownian motion is a scaling limit of random walks. Consider an urn containing   blue and yellow balls. At each step a ball is chosen at random and replaced by a ball of the opposite colour. Let   be the number of blue balls in the urn after   steps. Then   converges in law to an Ornstein–Uhlenbeck process as   tends to infinity. This was obtained by Mark Kac.[11]

Heuristically one may obtain this as follows.

Let  , and we will obtain the stochastic differential equation at the   limit. First deduce   With this, we can calculate the mean and variance of  , which turns out to be   and  . Thus at the   limit, we have  , with solution (assuming   distribution is standard normal)  .

Applications

edit

In physics: noisy relaxation

edit

The Ornstein–Uhlenbeck process is a prototype of a noisy relaxation process. A canonical example is a Hookean spring (harmonic oscillator) with spring constant   whose dynamics is overdamped with friction coefficient  . In the presence of thermal fluctuations with temperature  , the length   of the spring fluctuates around the spring rest length  ; its stochastic dynamics is described by an Ornstein–Uhlenbeck process with

 

where   is derived from the Stokes–Einstein equation   for the effective diffusion constant.[12][13] This model has been used to characterize the motion of a Brownian particle in an optical trap.[13][14]

At equilibrium, the spring stores an average energy   in accordance with the equipartition theorem.[15]

In financial mathematics

edit

The Ornstein–Uhlenbeck process is used in the Vasicek model of the interest rate.[16] The Ornstein–Uhlenbeck process is one of several approaches used to model (with modifications) interest rates, currency exchange rates, and commodity prices stochastically. The parameter   represents the equilibrium or mean value supported by fundamentals;   the degree of volatility around it caused by shocks, and   the rate by which these shocks dissipate and the variable reverts towards the mean. One application of the process is a trading strategy known as pairs trade.[17][18][19]

A further implementation of the Ornstein–Uhlenbeck process is derived by Marcello Minenna in order to model the stock return under a lognormal distribution dynamics. This modeling aims at the determination of confidence interval to predict market abuse phenomena.[20][21]

In evolutionary biology

edit

The Ornstein–Uhlenbeck process has been proposed as an improvement over a Brownian motion model for modeling the change in organismal phenotypes over time.[22] A Brownian motion model implies that the phenotype can move without limit, whereas for most phenotypes natural selection imposes a cost for moving too far in either direction. A meta-analysis of 250 fossil phenotype time-series showed that an Ornstein–Uhlenbeck model was the best fit for 115 (46%) of the examined time series, supporting stasis as a common evolutionary pattern.[23] This said, there are certain challenges to its use: model selection mechanisms are often biased towards preferring an OU process without sufficient support, and misinterpretation is easy to the unsuspecting data scientist.[24]

Generalizations

edit

It is possible to define a Lévy-driven Ornstein–Uhlenbeck process, in which the background driving process is a Lévy process instead of a Wiener process:[25][26]

 

Here, the differential of the Wiener process   has been replaced with the differential of a Lévy process  .

In addition, in finance, stochastic processes are used where the volatility increases for larger values of  . In particular, the CKLS process (Chan–Karolyi–Longstaff–Sanders)[27] with the volatility term replaced by   can be solved in closed form for  , as well as for  , which corresponds to the conventional OU process. Another special case is  , which corresponds to the Cox–Ingersoll–Ross model (CIR-model).

Higher dimensions

edit

A multi-dimensional version of the Ornstein–Uhlenbeck process, denoted by the N-dimensional vector  , can be defined from

 

where   is an N-dimensional Wiener process, and   and   are constant N×N matrices.[28] The solution is

 

and the mean is

 

These expressions make use of the matrix exponential.

The process can also be described in terms of the probability density function  , which satisfies the Fokker–Planck equation[29]

 

where the matrix   with components   is defined by  . As for the 1d case, the process is a linear transformation of Gaussian random variables, and therefore itself must be Gaussian. Because of this, the transition probability   is a Gaussian which can be written down explicitly. If the real parts of the eigenvalues of   are larger than zero, a stationary solution   moreover exists, given by

 

where the matrix   is determined from the Lyapunov equation  .[5]

See also

edit

Notes

edit
  1. ^ a b Doob 1942.
  2. ^ Karatzas & Shreve 1991, p. 358.
  3. ^ Gard 1988, p. 115.
  4. ^ Gardiner 1985.
  5. ^ a b c Risken 1989.
  6. ^ Lawler 2006.
  7. ^ Gardiner 1985, p. 106.
  8. ^ Holmes-Cerfon, Miranda (2022). "Lecture 12: Detailed balance and Eigenfunction methods" (PDF).
  9. ^ Aït-Sahalia 2002, pp. 223–262.
  10. ^ Kloeden, Platen & Schurz 1994.
  11. ^ Iglehart 1968.
  12. ^ Nørrelykke & Flyvbjerg 2011.
  13. ^ a b Goerlich et al. 2021.
  14. ^ Li et al. 2019.
  15. ^ Nelson 1967.
  16. ^ Björk 2009, pp. 375, 381.
  17. ^ Leung & Li 2016.
  18. ^ Advantages of Pair Trading: Market Neutrality
  19. ^ An Ornstein–Uhlenbeck Framework for Pairs Trading
  20. ^ "Detecting Market Abuse". Risk Magazine. 2 November 2004.
  21. ^ "The detection of Market Abuse on financial markets: a quantitative approach". Consob – The Italian Securities and Exchange Commission.
  22. ^ Martins 1994, pp. 193–209.
  23. ^ Hunt 2007.
  24. ^ Cornuault 2022.
  25. ^ Jespersen, Metzler & Fogedby 1999.
  26. ^ Fink & Klüppelberg 2011.
  27. ^ Chan et al. 1992.
  28. ^ Gardiner 1985, p. 109.
  29. ^ Gardiner 1985, p. 97.

References

edit
edit