Notes on CEP

Measurement is fundamental to virtually all human endeavors of practical importance. However, many measurements require accounting for error, and Circular Error Probable (CEP) is one method of quantifying deviations from a best estimate, commonly applied in the context of weapons systems and GPS. This article will derive CEP and explore its practical applications.

CEP is commonly defined as the radius of a circle in which 50% (i.e. the proportion \frac{1}{2} ) of observed data points are expected to fall. Our goal is to determine the radius q representing the circle enclosing the proportion Q \in [0,1] of the data.

To achieve this, lets summon the Rayleigh distribution, which will allow us to compute a good estimate of the radius. You can think of this distribution as a Gaussian distribution in which the random variable, in the case of the Gaussian would range from -\infty to \infty, would be compressed to the domain [0, \infty). The formula for this distribution is as follows:

\displaystyle P( R | \sigma^{2}) = \frac{R}{\sigma^{2}}e^{-\frac{R^{2}}{2\sigma^{2}}}

The density function itself resembles the form of the gaussian, in that they both belong to the exponential family of distributions. However, notice that this distribution is weighted by the variable R itself, imposing a sort of controlled “skewing” tendency on the graph. Examples of the probability density function are plotted below:

For a given Rayleigh distribution corresponding to the Gaussian variance \sigma^{2}, we can find the radius q which encircles Q percent of our data by adjusting the formula below:

\displaystyle \Huge Q = P(R \leq q) = \int \limits_{0}^{q} \frac{R}{\sigma^{2}}e^{-\frac{R^{2}}{2\sigma^{2}}}dR

\displaystyle \Huge = 1 - e^{- \frac{q^{2}}{2 \sigma^{2}}}

Therefore Q corresponds to the cumulative distribution function (CDF) of the Rayleigh distribution. By rearranging this CDF in terms of q, we obtain:

\displaystyle \Huge q = \sigma \sqrt{2 \ln{\frac{1}{1-Q}}}

In particular, for Q = 0.5, we obtain the CEP:

\displaystyle \Huge q = \sigma \sqrt{2 \ln{2}} = CEP

Consequently, the CEP is obtained by scaling the standard deviation by the factor \sqrt{2 \ln{2}}.

Leave a Reply

Your email address will not be published. Required fields are marked *