Fractional Brownian motion

From Wikipedia, the free encyclopedia

A normalized fractional Brownian motion (denoted fBm) BH(t) on [0,T], T\in \mathbb{R} is a continuous-time Gaussian process starting at zero, with mean zero, and having the following correlation function:

E[B^H(t) B^H(s)]=\frac{1}{2} (|t|^{2H}+|s|^{2H}-|t-s|^{2H}),

where H, called the Hurst index or Hurst parameter associated to the fractional Brownian motion, is a real number in [0,1].

The value of H determines what kind of process the fBm is:

  • if H = 1 / 2, the process is in fact a regular Brownian motion;
  • if H > 1 / 2, the increments of the process are positively correlated;
  • if H < 1 / 2, the increments of the process are negatively correlated.

Contents

[edit] Properties

[edit] Self-similarity

The process is said to be self-similar, since in terms of distributions:

B^H(at)\sim |a|^{H}B^H(t).

[edit] Stationary increments

It has stationary increments:

BH(t) − BH(s)˜BH(t − s).

[edit] Long-range dependence

For H > 1 / 2, the process exhibits long-range dependence, which means that

\sum_{n=1}^\infty{E[B^H(1)(B^H(n+1)-B^H(n))]}=\infty.


[edit] Regularity

Sample-paths are almost nowhere differentiable. Precisely, almost-all trajectories are Hölder continuous of any order strictly less than H: for each trajectory, there exists a constant c such that

|B^H(t)-B^H(s)| \le c |t-s|^{H-\varepsilon}

for any \varepsilon >0.

[edit] Integration

As for regular Brownian motion, one can define stochastic integrals with respect to fractional Brownian motion, usually called "fractional stochastic integrals". In general though, unlike integrals with respect to regular Brownian motion, fractional stochastic integrals are not Martingales.

[edit] Sample paths

Practical computer realisations of fBm can be generated, although obviously they cannot really be fractal. The sample paths chosen can be thought of as showing discrete sampled points on an fBm process. Three realisations are shown below, each with 1000 points of fBm with Hurst parameter 0.75.
H=0.75 realisation 1 H=0.75 realisation 2 H=0.75 realisation 3

Two realisations are shown below, each showing 1000 points of fBm, the first with Hurst parameter 0.95 and the second with Hurst parameter 0.55.
H=0.95 H=0.55


[edit] Method 1 of simulation :

One can simulate sample-paths of fBm as any Gaussian process of known covariance. Say we aim to have simulated values at t_1, \cdots, t_n.

  • Form the matrix \Gamma=\bigl(R(t_i,\, t_j), i,j=1,\cdots,\, n\bigr) where R(t,s) = (s2H + t2H − | t − s | 2H) / 2.
  • Construct a vector of n numbers drawn according a standard gaussian distribution.
  • Apply σ to this vector yields a sample-path of fBm.

[edit] Method 2 of simulation :

It is also known that B^H(t)=\int_0^t K_H(t,s)dB(s) where B is a standard Brownian motion and K_H(t,s)=\frac{1}{\Gamma(H+1/2)}(t-s)^{H-1/2}\;_2F_1(H-1/2;\, 1/2-H;\; H+1/2;\, 1-t/s). See section on Euler Hypergeometric integral for the definition of 2F1.


Say we want simulate fBm at points 0=t_0< t_1< \cdots < t_n=T.

  • Construct a vector of n numbers drawn according a standard gaussian distribution.
  • Multiply it component-wise by sqrt(T/n) to obtain the increments of a Brownian motion on [0,T]. Denote this vector by (\delta B_1, \cdots, \delta B_n).
  • For each tj, compute

B^H(t_j)=\frac{n}{T}\sum_{i=0}^{j-1} \int_{t_i}^{t_{i+1}} K_H(t_j,\, s)\, ds \ \delta B_i.

The integral may be efficiently computed by Gaussian quadrature. Hypergeometric functions are part of the GNU scientific library .

[edit] Applications

Fractional Brownian motion was initially used in the modelling of hydrological phenomena. It has been shown that fBm can be suitable for the analysis of computer network traffic. By the end of the twentieth century, it was used in as various fields as financial mathematics or random landscape generation. Currently among the popular methods for the analysis of these processes are Wavelets-based methods.