Physics:Circular error probable

Short description: Ballistics measure of a weapon system's precision

CEP concept and hit probability. 0.2% outside the outmost circle.

Circular error probable (CEP), also circular error probability or circle of equal probability, is a measure of a weapon system's precision in the military science of ballistics. It is defined as the radius of a circle, centered on the aimpoint, that is expected to enclose the landing points of 50% of the rounds; said otherwise, it is the median error radius. That is, if a given munitions design has a CEP of 100 m, when 100 munitions are targeted at the same point, an average of 50 will fall within a circle with a radius of 100 m about that point.

There are associated concepts, such as the DRMS (distance root mean square), which is the square root of the average squared distance error, and R95, which is the radius of the circle where 95% of the values would fall in.

The concept of CEP also plays a role when measuring the accuracy of a position obtained by a navigation system, such as GPS or older systems such as LORAN and Loran-C.

Concept

circular bivariate normal distribution

20 hits distribution example

The original concept of CEP was based on a circular bivariate normal distribution (CBN) with CEP as a parameter of the CBN just as μ and σ are parameters of the normal distribution. Munitions with this distribution behavior tend to cluster around the mean impact point, with most reasonably close, progressively fewer and fewer further away, and very few at long distance. That is, if CEP is n metres, 50% of shots land within n metres of the mean impact, 43.7% between n and 2n, and 6.1% between 2n and 3n metres, and the proportion of shots that land farther than three times the CEP from the mean is only 0.2%.

CEP is not a good measure of accuracy when this distribution behavior is not met. Precision-guided munitions generally have more "close misses" and so are not normally distributed. Munitions may also have larger standard deviation of range errors than the standard deviation of azimuth (deflection) errors, resulting in an elliptical confidence region. Munition samples may not be exactly on target, that is, the mean vector will not be (0,0). This is referred to as bias.

To incorporate accuracy into the CEP concept in these conditions, CEP can be defined as the square root of the mean square error (MSE). The MSE will be the sum of the variance of the range error plus the variance of the azimuth error plus the covariance of the range error with the azimuth error plus the square of the bias. Thus the MSE results from pooling all these sources of error, geometrically corresponding to radius of a circle within which 50% of rounds will land.

Several methods have been introduced to estimate CEP from shot data. Included in these methods are the plug-in approach of Blischke and Halpin (1966), the Bayesian approach of Spall and Maryak (1992), and the maximum likelihood approach of Winkler and Bickert (2012). The Spall and Maryak approach applies when the shot data represent a mixture of different projectile characteristics (e.g., shots from multiple munitions types or from multiple locations directed at one target).

Conversion

While 50% is a very common definition for CEP, the circle dimension can be defined for percentages. Percentiles can be determined by recognizing that the horizontal position error is defined by a 2D vector which components are two orthogonal Gaussian random variables (one for each axis), assumed uncorrelated, each having a standard deviation [math]\displaystyle{ \sigma }[/math]. The distance error is the magnitude of that vector; it is a property of 2D Gaussian vectors that the magnitude follows the Rayleigh distribution, with a standard deviation [math]\displaystyle{ \sigma_d=\sqrt{2}\sigma }[/math], called the distance root mean square (DRMS). In turn, the properties of the Rayleigh distribution are that its percentile at level [math]\displaystyle{ F \in [0\%, 100\%] }[/math] is given by the following formula:

[math]\displaystyle{ Q(F, \sigma) = \sigma \sqrt{-2\ln(1 - F/100\%)} }[/math]

or, expressed in terms of the DRMS:

[math]\displaystyle{ Q(F, \sigma_d) = \sigma_d \frac{\sqrt{-2\ln(1 - F/100\%)}}{\sqrt{2}} }[/math]

The relation between [math]\displaystyle{ Q }[/math] and [math]\displaystyle{ F }[/math] are given by the following table, where the [math]\displaystyle{ F }[/math] values for DRMS and 2DRMS (twice the distance root mean square) are specific to the Rayleigh distribution and are found numerically, while the CEP, R95 (95% radius) and R99.7 (99.7% radius) values are defined based on the 68–95–99.7 rule

Measure of [math]\displaystyle{ Q }[/math]	Probability [math]\displaystyle{ F \, (\%) }[/math]
DRMS	63.213...
CEP	50
2DRMS	98.169...
R95	95
R99.7	99.7

We can then derive a conversion table to convert values expressed for one percentile level, to another. Said conversion table, giving the coefficients [math]\displaystyle{ \alpha }[/math] to convert [math]\displaystyle{ X }[/math] into [math]\displaystyle{ Y=\alpha.X }[/math], is given by:

From [math]\displaystyle{ X \downarrow }[/math] to [math]\displaystyle{ Y \rightarrow }[/math]	RMS ([math]\displaystyle{ \sigma }[/math])	CEP	DRMS	R95	2DRMS	R99.7
RMS ([math]\displaystyle{ \sigma }[/math])	1.00	1.18	1.41	2.45	2.83	3.41
CEP	0.849	1.00	1.20	2.08	2.40	2.90
DRMS	0.707	0.833	1.00	1.73	2.00	2.41
R95	0.409	0.481	0.578	1.00	1.16	1.39
2DRMS	0.354	0.416	0.500	0.865	1.00	1.21
R99.7	0.293	0.345	0.415	0.718	0.830	1.00

For example, a GPS receiver having a 1.25 m DRMS will have a 1.25 m × 1.73 = 2.16 m 95% radius.

Warning: often, sensor datasheets or other publications state "RMS" values which in general, but not always, stand for "DRMS" values. Also, be wary of habits coming from properties of a 1D normal distribution, such as the 68-95-99.7 rule, in essence trying to say that "R95 = 2DRMS". As shown above, these properties simply do not translate to the distance errors. Finally, mind that these values are obtained for a theoretical distribution; while generally being true for real data, these may be affected by other effects, which the model does not represent.

See also

• Probable error

Further reading

• Blischke, W. R.; Halpin, A. H. (1966). "Asymptotic Properties of Some Estimators of Quantiles of Circular Error". Journal of the American Statistical Association 61 (315): 618–632. doi:10.1080/01621459.1966.10480893. 
• Grubbs, F. E. (1964). "Statistical measures of accuracy for riflemen and missile engineers". Ann Arbor, ML: Edwards Brothers. Ballistipedia pdf
• MacKenzie, Donald A. (1990). Inventing Accuracy: A Historical Sociology of Nuclear Missile Guidance. Cambridge, Massachusetts: MIT Press. ISBN 978-0-262-13258-9. https://archive.org/details/inventingaccurac00dona. 
• Spall, James C.; Maryak, John L. (1992). "A Feasible Bayesian Estimator of Quantiles for Projectile Accuracy from Non-iid Data". Journal of the American Statistical Association 87 (419): 676–681. doi:10.1080/01621459.1992.10475269. 
• Winkler, V. and Bickert, B. (2012). "Estimation of the circular error probability for a Doppler-Beam-Sharpening-Radar-Mode," in EUSAR. 9th European Conference on Synthetic Aperture Radar, pp. 368–71, 23/26 April 2012. ieeexplore.ieee.org
• Wollschläger, Daniel (2014), "Analyzing shape, accuracy, and precision of shooting results with shotGroups". Reference manual for shotGroups

External links

• Circular Error Probable in Ballistipedia