Circular uniform distribution

From Wikipedia, the free encyclopedia
Jump to: navigation, search

In probability theory and directional statistics, a circular uniform distribution is a probability distribution on the unit circle whose density is uniform for all angles.

Description[edit]

The probability density function (pdf) of the circular uniform distribution is:


f_{UC}(\theta)=\frac{1}{2\pi}.

In terms of the circular variable z=e^{i\theta} the circular moments of the circular uniform distribution all zero, except for m_0:

\langle z^n\rangle=\delta_n

where \delta_n is the Kronecker delta symbol.

The mean angle is undefined, and the length of the mean resultant is zero.


R=|\langle z^n\rangle|=0\,

Distribution of the mean[edit]

The sample mean of a set of N measurements z_n=e^{i\theta_n} drawn from a circular uniform distribution is defined as:


\overline{z} = \frac{1}{N}\sum_{n=1}^N z_n = \overline{C}+i\overline{S} = \overline{R}e^{i\overline{\theta}}

where the average sine and cosine are:


\overline{C}=\frac{1}{N}\sum_{n=1}^N \cos(\theta_n)\qquad\qquad\overline{S}=\frac{1}{N}\sum_{n=1}^N \sin(\theta_n)

and the average resultant length is:


\overline{R}^2=|\overline{z}|^2=\overline{C}^2+\overline{S}^2

and the mean angle is:


\overline{\theta}=\mathrm{Arg}(\overline{z}). \,

The sample mean for the circular uniform distribution will be concentrated about zero, becoming more concentrated as N increases. The distribution of the sample mean for the uniform distribution is given by:[1]


\frac{1}{(2\pi)^N}\int_\Gamma \prod_{n=1}^N d\theta_n = P(\overline{R})P(\overline{\theta})\,d\overline{R}\,d\overline{\theta}

where \Gamma\, consists of intervals of 2\pi in the variables, subject to the constraint that \overline{R} and \overline{\theta} are constant, or, alternatively, that \overline{C} and \overline{S} are constant. The distribution of the angleP(\overline{\theta}) is uniform


P(\overline{\theta})=\frac{1}{2\pi}

and the distribution of \overline{R} is given by:[1]


P_N(\overline{R})=N^2\overline{R}\int_0^\infty J_0(N\overline{R}\,t)J_0(t)^Nt\,dt
A 10,000 point Monte Carlo simulation of the distribution of the sample mean of a circular uniform distribution for N = 3

where J_0 is the Bessel function of order zero. There is no known general analytic solution for the above integral, and it is difficult to evaluate due to the large number of oscillations in the integrand. A 10,000 point Monte Carlo simulation of the distribution of the mean for N=3 is shown in the figure.

For certain special cases, the above integral can be evaluated:


P_2(\overline{R})=\frac{2}{\pi \sqrt{1-\overline{R}^2}}.

For large N, the distribution of the mean can be determined from the central limit theorem for directional statistics. Since the angles are uniformly distributed, the individual sines and cosines of the angles will be distributed as:


P(u)du=\frac{1}{\pi}\,\frac{du}{\sqrt{1-u^2}}

where u=\cos\theta_n\, or \sin\theta_n\,. It follows that they will have zero mean and a variance of 1/2. By the central limit theorem, in the limit of large N, \overline{C}\, and \overline{S}\,, being the sum of a large number of i.i.d's, will be normally distributed with mean zero and variance 1/2N. The mean resultant length \overline{R}\,, being the square root of the sum of two normally distributed variables, will be Chi-distributed with two degrees of freedom (i.e.Rayleigh-distributed) and variance 1/2N:


\lim_{N\rightarrow\infty}P_N(\overline{R})=2N\overline{R}\,e^{-N\overline{R}^2}.

Entropy[edit]

The differential information entropy of the uniform distribution is simply

H_U=-\int_\Gamma \frac{1}{2\pi}\ln\left(\frac{1}{2\pi}\right)\,d\theta = \ln(2\pi)

where \Gamma is any interval of length 2\pi. This is the maximum entropy any circular distribution may have.

See also[edit]

References[edit]

  1. ^ a b Jammalamadaka, S. Rao; Sengupta, A. (2001). Topics in Circular Statistics. World Scientific Publishing Company. ISBN 978-981-02-3778-3. Retrieved 2010-03-03.