St · Statistics & Probability

What Is the Difference Between the Rician and Gaussian Distributions?

Published By the 9to5Science editorial team

At a glance

  1. 1. Identify the variable

    A Gaussian variable can take any real value; a Rician variable is a nonnegative magnitude.

  2. 2. Model Gaussian components

    Take two independent Gaussian components with equal variance and a possibly nonzero mean vector.

  3. 3. Compute the magnitude

    The square root of the sum of the components’ squares has a Rician distribution.

  4. 4. Check special cases

    A zero mean vector gives a Rayleigh distribution; a strong mean can make the magnitude approximately Gaussian.

The key difference is what each distribution models. A Gaussian (normal) distribution describes a real-valued quantity that can lie anywhere on the real line. A Rician distribution describes a magnitude: it is the length of a two-dimensional Gaussian vector, so it cannot be negative.

This connection explains why a Rician model appears in applications such as radio signal amplitudes. The signal’s two components can each contain Gaussian noise, even though the measured amplitude itself is not Gaussian.

How the Rician distribution arises

Let XX and YY be independent Gaussian random variables with the same variance. Their means can be viewed as coordinates of a vector with length ν\nu:

X∼N(ν,σ2),Y∼N(0,σ2).X \sim N(\nu,\sigma^2), \qquad Y \sim N(0,\sigma^2).

The magnitude

R=X2+Y2R=\sqrt{X^2+Y^2}

has a Rician distribution with parameters ν\nu and σ\sigma. The direction of the mean vector does not affect the distribution of the magnitude, so it is convenient to place that mean along the XX-axis as above. Here, ν≥0\nu\geq 0 is the length of the mean vector, and σ\sigma is the standard deviation of either Gaussian component.

For r≥0r\geq 0, its probability density function is

fR(r)=rσ2exp⁡(−r2+ν22σ2)I0(rνσ2),f_R(r)=\frac{r}{\sigma^2}\exp\left(-\frac{r^2+\nu^2}{2\sigma^2}\right) I_0\left(\frac{r\nu}{\sigma^2}\right),

and fR(r)=0f_R(r)=0 for r<0r<0. The function I0I_0 is the modified Bessel function of the first kind of order zero. It appears because the joint Gaussian density is being converted from Cartesian coordinates to a radial magnitude.

Main differences at a glance

FeatureGaussian variableRician magnitude
Values it can takeAny real numberNonnegative numbers only
Typical constructionOne real-valued measurementLength of a two-component Gaussian vector
Parameters in this modelMean and standard deviationMean-vector length ν\nu and component standard deviation σ\sigma
ShapeSymmetric about its meanGenerally asymmetric; shape depends on ν/σ\nu/\sigma

A Gaussian distribution may describe either component, XX or YY, in the construction. It does not generally describe their magnitude RR.

Special cases and useful relationships

When ν=0\nu=0, both components are centered at zero. The Rician distribution then reduces to the Rayleigh distribution, with density

fR(r)=rσ2exp⁡(−r22σ2),r≥0.f_R(r)=\frac{r}{\sigma^2}\exp\left(-\frac{r^2}{2\sigma^2}\right), \qquad r\geq 0.

More generally, the Rician distribution is a noncentral chi distribution with two degrees of freedom. This provides a useful link to the chi and chi-squared distributions.

The squared magnitude also has a simple expectation:

E[R2]=E[X2+Y2]=ν2+2σ2.\mathbb{E}[R^2]=\mathbb{E}[X^2+Y^2]=\nu^2+2\sigma^2.

This is the mean square magnitude, not the square of the mean magnitude: in general, E[R2]≠(E[R])2\mathbb{E}[R^2]\ne (\mathbb{E}[R])^2.

Worked example

Suppose X∼N(3,1)X\sim N(3,1) and Y∼N(0,1)Y\sim N(0,1) independently, and measure R=X2+Y2R=\sqrt{X^2+Y^2}. The component standard deviation is σ=1\sigma=1, and the mean vector has length ν=3\nu=3. Therefore, RR is Rician with parameters (3,1)(3,1).

At magnitude r=3r=3, its density is

fR(3)=3exp⁡(−32+322)I0(9)=3e−9I0(9)≈0.405.\begin{aligned} f_R(3) &=3\exp\left(-\frac{3^2+3^2}{2}\right)I_0(9)\\ &=3e^{-9}I_0(9)\\ &\approx 0.405. \end{aligned}

The density is zero at r=0r=0 for these parameters, and it is zero for every negative value. By contrast, the Gaussian component XX can be negative: P(X<0)=P(Z<−3)≈0.00135P(X<0)=P(Z<-3)\approx 0.00135 for a standard normal variable ZZ. These are different random variables, so their probability distributions need not match.

Also, the mean square magnitude is E[R2]=32+2(1)2=11\mathbb{E}[R^2]=3^2+2(1)^2=11, giving a root-mean-square magnitude of 11≈3.32\sqrt{11}\approx 3.32.

When can a Rician variable be approximated by a Gaussian?

If ν\nu is much larger than σ\sigma, most of the vector’s length comes from its nonzero mean rather than the noise. In that regime, the magnitude is often approximately normal, centered near ν\nu with spread near σ\sigma. The approximation is less suitable when ν\nu is comparable to σ\sigma, and a Gaussian approximation can still assign probability to negative values even though a magnitude cannot be negative.

Common mistakes

  • Calling the magnitude Gaussian because its components are Gaussian. A nonlinear transformation such as X2+Y2\sqrt{X^2+Y^2} changes the distribution.
  • Treating a Rician variable as able to be negative. It represents a length, so its support is [0,∞)[0,\infty).
  • Assuming every pair of Gaussian components gives a Rician magnitude. The standard Rician model requires independent components with equal variance. Unequal variances or dependence generally lead to a different model.
  • Mixing parameter conventions. Here, σ\sigma is the standard deviation of each component and ν\nu is the mean-vector length. Check definitions when using a source or software package, since notation and parameterizations can vary.

References

  1. Continuous Univariate Distributions, Volume 1 — N. L. Johnson, S. Kotz, and N. Balakrishnan
  2. Probability, Random Variables, and Stochastic Processes — Athanasios Papoulis and S. Unnikrishna Pillai

Spotted an error? Send a correction. Read how we write and check articles in our editorial policy.