Landau distribution

From testwiki
Jump to navigation Jump to search

Template:Probability distribution

In probability theory, the Landau distribution[1] is a probability distribution named after Lev Landau. Because of the distribution's "fat" tail, the moments of the distribution, such as mean or variance, are undefined. The distribution is a particular case of stable distribution.

Definition

The probability density function, as written originally by Landau, is defined by the complex integral:

p(x)=12πiaia+ieslog(s)+xsds,

where a is an arbitrary positive real number, meaning that the integration path can be any parallel to the imaginary axis, intersecting the real positive semi-axis, and log refers to the natural logarithm. In other words it is the Laplace transform of the function ss.

The following real integral is equivalent to the above:

p(x)=1π0etlog(t)xtsin(πt)dt.

The full family of Landau distributions is obtained by extending the original distribution to a location-scale family of stable distributions with parameters α=1 and β=1,[2] with characteristic function:[3]

φ(t;μ,c)=exp(itμ2ictπlog|t|c|t|)

where c(0,) and μ(,), which yields a density function:

p(x;μ,c)=1πc0etcos(t(xμc)+2tπlog(tc))dt,

Taking μ=0 and c=π2 we get the original form of p(x) above.

Properties

The approximation function for μ=0,c=1
  • Translation: If XLandau(μ,c) then X+mLandau(μ+m,c).
  • Scaling: If XLandau(μ,c) then aXLandau(aμ2aclog(a)π,ac).
  • Sum: If XLandau(μ1,c1) and YLandau(μ2,c2) then X+YLandau(μ1+μ2,c1+c2).

These properties can all be derived from the characteristic function. Together they imply that the Landau distributions are closed under affine transformations.

Approximations

In the "standard" case μ=0 and c=π/2, the pdf can be approximated[4] using Lindhard theory which says:

p(x+log(x)1+γ)exp(1/x)x(1+x),

where γ is Euler's constant.

A similar approximation [5] of p(x;μ,c) for μ=0 and c=1 is:

p(x)12πexp(x+ex2).
  • The Landau distribution is a stable distribution with stability parameter α and skewness parameter β both equal to 1.

References

Template:Reflist

Template:ProbDistributions