Difference between revisions of "Elliptic Error Probable"

From ShotStat
Jump to: navigation, search
 
Line 1: Line 1:
 
<p style="text-align:right"><B>Previous:</B> [[Precision Models]]</p>
 
<p style="text-align:right"><B>Previous:</B> [[Precision Models]]</p>
  
See section 2.3 of http://www.arl.army.mil/arlreports/2013/ARL-TR-6494.pdf for rectangular and elliptical statistics on asymmetric bivariate normal distributions.
+
= Preserving Asymmetry =
 +
There are two significant sources of dispersion asymmetry in ballistics:
 +
# Crosswind variance will cause excess dispersion only in the horizontal axis.
 +
# Muzzle velocity variance will cause excess dispersion that appears only in the vertical axis, and which can become quite significant at longer shooting distances.
  
If we are going to make separate estimates of variance in each axis, <math>s_x^2, s_y^2</math> note that the confidence interval in one dimension uses <math>\chi^2</math> values with (''n'' - 1) degrees of freedom, and:
+
Given (''x, y'') coordinate data of impacts on a target, we can preserve asymmetry in the precision models and their application by treating each axis separately.  [[Closed_Form_Precision#Variance_Estimates|As described in Closed Form Precision]], for variance in one dimension:
:&nbsp; <math>\sigma_x^2 \in \left[\, \frac{(n-1)s_x^2}{\chi_2^2}, \ \frac{(n-1)s_x^2}{\chi_1^2} \,\right]</math>
+
 
 +
:&nbsp; <math>\hat{\sigma}_x^2 = s_x^2 = \frac{1}{n-1}\sum{(x_i - \bar{x})^2}</math>
 +
 
 +
The confidence interval for variance in one dimension uses <math>\chi^2</math> values with (''n'' - 1) degrees of freedom:
 +
:&nbsp; <math>\hat{\sigma}_x^2 \in \left[\, \frac{(n-1)s_x^2}{\chi_{\frac{\alpha}{2},n-1}^2}, \ \frac{(n-1)s_x^2}{\chi_{1-\frac{\alpha}{2},n-1}^2} \,\right]</math>
 +
 
 +
As before, the estimates can be brought into unbiased standard deviation terms by multiplying the square root of estimated variance by the [[Closed_Form_Precision#Gaussian_correction_factor|Gaussian correction term]]: <math>\hat{\sigma}_x=c_G(n) \sqrt{s_x^2}</math>
 +
 
 +
= Covering Ellipse =
 +
The radii (''a, b'') of an ellipse that will cover proportion ''p'' of shots are given by the Rayleigh quantile function:
 +
* <math>a=\sigma_x \sqrt{-2 \ln(1-p)}</math>
 +
* <math>b=\sigma_y \sqrt{-2 \ln(1-p)}</math>
 +
 
 +
For example, here is a simulation of shots with <math>\sigma_x=3, \sigma_y=1</math> overlaid with the 50% covering ellipse:
 +
 
 +
[[File:Ellipse50CEP.png|470px]]
 +
 
 +
= Other target shapes =
 +
As detailed in Section 2.3 of [https://apps.dtic.mil/sti/tr/pdf/ADA588846.pdf ARL-TR-6494]<ref>Strohm, Luke (2013) ''An Introduction to the Sources of Delivery Error for Direct-Fire Ballistic Projectiles''</ref>, numerical integration can calculate the hit probabilities for arbitrary target shapes and sizes.

Latest revision as of 16:37, 20 November 2023

Previous: Precision Models

Preserving Asymmetry

There are two significant sources of dispersion asymmetry in ballistics:

  1. Crosswind variance will cause excess dispersion only in the horizontal axis.
  2. Muzzle velocity variance will cause excess dispersion that appears only in the vertical axis, and which can become quite significant at longer shooting distances.

Given (x, y) coordinate data of impacts on a target, we can preserve asymmetry in the precision models and their application by treating each axis separately. As described in Closed Form Precision, for variance in one dimension:

  \(\hat{\sigma}_x^2 = s_x^2 = \frac{1}{n-1}\sum{(x_i - \bar{x})^2}\)

The confidence interval for variance in one dimension uses \(\chi^2\) values with (n - 1) degrees of freedom:

  \(\hat{\sigma}_x^2 \in \left[\, \frac{(n-1)s_x^2}{\chi_{\frac{\alpha}{2},n-1}^2}, \ \frac{(n-1)s_x^2}{\chi_{1-\frac{\alpha}{2},n-1}^2} \,\right]\)

As before, the estimates can be brought into unbiased standard deviation terms by multiplying the square root of estimated variance by the Gaussian correction term\[\hat{\sigma}_x=c_G(n) \sqrt{s_x^2}\]

Covering Ellipse

The radii (a, b) of an ellipse that will cover proportion p of shots are given by the Rayleigh quantile function:

  • \(a=\sigma_x \sqrt{-2 \ln(1-p)}\)
  • \(b=\sigma_y \sqrt{-2 \ln(1-p)}\)

For example, here is a simulation of shots with \(\sigma_x=3, \sigma_y=1\) overlaid with the 50% covering ellipse:

Ellipse50CEP.png

Other target shapes

As detailed in Section 2.3 of ARL-TR-6494[1], numerical integration can calculate the hit probabilities for arbitrary target shapes and sizes.

  1. Strohm, Luke (2013) An Introduction to the Sources of Delivery Error for Direct-Fire Ballistic Projectiles