CEP Calculation Based on Weighted Bayesian Mixture Model

CEP is an important accuracy index in the test evaluation for guidance weapon systems. As for the samples origin from diverse populations, which is a general case in the practice, traditional method with one population is inaccurate to give estimate directly. A weighted method considering the credibility of the prior information is proposed for Bayesian estimation algorithm and the weighted estimation of normal distribution parameters is provided. And the statistical diversity between the weighted and classical methods is quantified, and upper bound of the calculation errors caused by prior distortion is deduced for mixture of finite normal distributions. Taking the estimation of the dispersion of the guided weapon impact points as an example, the conclusion is drawn that our method is credible.


Introduction
Circular Error Probable (CEP) is a common measure used in the evaluation test of guidance weapon systems, which can integrate precision with dispersion to assess the hit accuracy of projectiles [1].In traditional CEP assessment process, the impact point deviation is assumed to follow bivariate normal distribution, and based on the formula simplified by decorrelation coordinate transformation, the MLE (maximum likelihood estimator) of normal distribution parameters are calculated by specific samples, which are employed in the forthcoming point estimation, confidence bounds computation, and hypothesis testing of CEP [2].
There are two restrictions for above CEP evaluation procedure.One is the concrete computational difficulty.It is hard to solve the complex CEP equation precisely, so some approximative expressions are adopted in practical applications with different assumptions [2].In [3] 12 methods of calculating CEP are summarized and classified into three categories: univariate parameterization, binary computation, and numerical integration.
Another limitation is adaptability of the estimate process.The basis of bivariate normal distribution assumption is the central-limit theorem and engineering experience accumulated for ages [4].Compatibility and normality test are called for the hypothesis based on a certain confidence level.But the assumption is not suitable for complicated data, which demands expanding the ranges of CEP application to obtain more adaptable and reasonable results of hit accuracy [1,2,5,6].A bias-corrected estimator of CEP is provided based on Cornish-Fisher expansion in [2].And in [5,6], CEP is extended to ACEP (Area Circular Error Probability) to overcome the disadvantage that CEP is not combined with the self-characteristics of surface targets.Reference [7] considers the case that missile impact samples is not i.i.d. and provides a Bayesian procedure for estimating CEP.
For INS (Inertial Navigation System)/SAR (Synthetic Aperture Radar) integrated system, the system biases of impact point errors are quite distinct in diverse experiment conditions or scene features.Under the circumstances classical CEP evaluation is not applicable because it is tough to treat various samples as generated by the same normal population.In addition, the stochastic quality of scene matching also increases computational difficulty.In order to use the samples following diverse populations to assess accuracy of projectiles, this paper provides the CEP evaluation procedure under the circumstance that impact samples are not i.i.d.data.In Section 2, classical CEP is generalized to the CEP of multiconditional probability based on nonparametric mixture model.Then specific computational method for mixed CEP model is designed in Section 3, which includes simplification of general equation, point, and 2 Mathematical Problems in Engineering interval estimation.The algorithms of mixed and classical CEP are compared and analyzed in Section 4. In Section 5 some simulation experiments are implemented, the results of which indicate that the mixed CEP algorithm is more reasonable than traditional one under the circumstances that diversity of accuracy samples is obvious.

Modeling of Bayes Model of Scene Matching of Multiconditional Probability
2.1.Classical CEP Definition.The definition of CEP is distinct in kinds of literature, and a common one is adopted in this paper.CEP is the radius of a circle, where center is target and the probability of an impact point inside is 50% [1].
The corresponding mathematical description is as follows.Establish a rectangular coordinate system centered on target with (0, 0), and assume the deviation of downrange and crossrange noted as  and , respectively, both follows normal distribution; then the probability density function of  = (, )  is where   ,   are standard deviations of , ;   ,   are means of , ;  is correlation coefficient of  and , such that 0 ≤ || < 1.Then  in the following formula is the CEP of  = (, ) where (, ) is the same as (1).And (2) is called the general form of CEP equation.Suppose {u  | u  = (  ,   )  }  =1 are the impact deviation samples to be used.In previous CEP algorithms they are utilized to calculate   ,   ,   ,   ,  in (1), which are employed to estimate  by numerical integration or approximate approaches in (2).But samples collected in different test environments and influenced by complicate factors such as scene features are distinct, which is against the assumption of same population, and make the traditional CEP evaluation invalid.Under this condition we propose the mixed CEP model and algorithm of multiconditional probability.

Bayes Nonparametric Mixture Model.
The Bayes mixture model combined with random density function is as follows: where (; ) is density function with different parameters .
And () is random density function as the following discrete form usually: where (  ,   ) ∞ =1 is prior information and    is characteristic function with various .

Mixed CEP Model for Scene Matching of Multiconditional
Probability.Define mixed CEP as the radius of a circle centered on the target, such that the probability of an impact point inside the circle is 50%, while impact deviation does not follow the same bivariate normal population.For INS/SAR integrated navigation system, the impact biases are influenced by scene matching, of which the error sources such as scene features, types, and matching number are quite distinct and follow diverse populations [8].In order to assess hit accuracy of the system by mixed CEP, modeling for scene matching is essential beforehand.Since the types and occurring probability of scene in one test are unknown before the experiment implement, they can be regarded as random variables while modeling based on Bayes nonparametric model frame.
Suppose synthetical test conditions of scene matching are classified into  types denoted as ( Consider the conditional probability density function of impact errors on ℎ  is normal as follows: where    ,    are standard deviations of ,  on condition Then the joint probability density function of impact errors is The discrete form is Substitute (8) for (, ) in ( 2), and  in the following equation is the mixed CEP of multiconditional probability

Mixed CEP Algorithm of Multiconditional Probability
3.1.Simplification of General CEP Equation.For condition ℎ  , take orthogonal transformation for (6) as follows: where Then (, ) is transformed to (, ), of which  and  are independent of condition ℎ  .And probability density function of (, ) on condition ℎ  is where    ,  V  are standard deviations of (, ) on condition ℎ  ;    ,  V  are means of (, ) on condition ℎ  .Then ( 9) is turned into Take polar coordinate transformation as follows: and then (13) is simplified from the form of double integral to the following quadratic integral: where where Suppose prior information of text conditions is explicit; that is, where   are all constants, such that Substitute μ  , μV  , σ  , σV  ,   for    ,  V  ,    ,  V  , [(ℎ  )] in (15); then mixed CEP can be calculated approximately with numerical integration and bisection algorithm.

Mixed CEP Interval Estimation of Multiconditional Probability.
Bootstrap method [9,10] is called for calculating mixed CEP confidence interval estimation of multiconditional probability here.Suppose R is the plug-in point estimator of .  denote  quantile of ; that is, ( ≤   ) = , and corresponding estimator is Ĉ .Then the calculation formulas of confidence interval [  ,   ] and confidence upper bound are given as follows: where Ĉ1−/2 , Ĉ/2 , Ĉ are approximately computed by the following steps.

Performance Analysis of the Algorithm
Thus the errors of single normal population approximated to the bivariate mixed population of two normal distributions are where By (25), it indicates that the means of two algorithms are the same when μ = 0 or   =   /, ∀ = 1, 2, . . ., .But as is shown in (26), variances of two algorithms are distinct under these circumstances.
As displayed in Figure 1, the mixed population is still unimodal, of which the theoretical, classical, and mixed CEP are 1.414 m, 1.215 m, and 1.420 m.Take sample from normal populations 1 and 2 with sizes 20 and 30, respectively, of which the dispersion is shown in Figure 2. In Figure 2   (16%) impact in the circle of classical CEP, which does not accord with the definition of CEP.Meanwhile there are 28 samples (56%) in the circle of mixed CEP, which illustrates that the mixed CEP calculation algorithm in this paper is more adaptable and accurate contrasted with theoretical and classical CEP under the circumstance.
In order to reduce the influence of random factors to CEP calculation result, resample 10 and 100 times under the circumstance that weights, sample size, and parameters of two normal populations are fixed, and the averages of corresponding computing results are shown in Table 1.
In Table 1, it indicates that the average precision of mixed CEP algorithm is larger than classical one on the condition

Mean of real probability
Prior mean of probability of resampling, and the probability when mixed CEP is more precise than typical one is large (about 77% in resampling 100 times) under the circumstance that samples are fixed.In addition, viewed from the perspective of containing samples, the proportion of samples in the circle of mixed CEP much more comply with the CEP definition, of which the average is about 44% in resampling times.Meanwhile the corresponding value of classical CEP circle is far from the ideal 50%, of which the average is about 21%.

Example 2.
The influence of prior distortion to mixed CEP errors for two normal mixture population is discussed in this subsection, of which the related parameters are the same as Section 5.1.In this case the input of prior parameter is only p1 ∈ [0, 1], while p2 can be replaced by p2 = 1 − p1 .So the relationship between p1 and mixed CEP can be expressed as the following equation: where  1 and  2 denote density functions of normal distribution 1 and 2. Suppose p1 = 0.4 is the real prior information, and the exact mixed CEP is 1.414 m.Then CEP and corresponding errors are calculated on the conditions of p1 = 0.1 ( = 0, 1, . . ., 10), which are shown in Table 2.
As is shown in Table 2, the CEP error Δ R increases as the prior error Δp 1 increases.And we can see the maximum of CEP under prior distortion is achieved by the boundary point p1 = 0.4, which is in accordance with Proposition 1.
In addition, we can conclude that prior distortion can deviate the mixed CEP from theoretic value indeed, but the upper bound of errors is estimable, which is quantified in Proposition 1.

Example 3.
The assessment of CEP in the case of small samples is discussed in this subsection.Assume 9 impact error samples are obtained after a group of flight tests and divided into 3 classes under corresponding conditions, which is shown in Table 3.In addition, real and prior mean of probability of test conditions are given in Table 4. Then the CEP point estimators and confidence intervals can be calculated after resampling 1000 times under the different assumptions, which is shown in Table 5.
The CEP confidence upper bounds under different assumptions are compared and analyzed here.As is shown in Table 5, the assessment results by mixed CEP method are closer to the mixed CEP of accurate prior information compared with classical single algorithm, even with the distorted prior.It illustrates that the adaptive ability of mixed CEP algorithm is promising in the case of small samples.

Conclusions
We have investigated the problem that classical CEP method may be unreasonable and may give inaccurate results for the impact samples following different populations.For purpose of addressing it, the mixed CEP is modeled in the context that the errors caused by scene matching follows diverse populations, and corresponding algorithms of point estimator and confidence interval estimation are designed.Taking mixed population of two normal distributions as an example, the statistical properties of traditional and mixed CEP algorithm are compared, and it concludes that prior distortion can deviate the mixed CEP from theoretic value indeed, but the upper bound of errors is estimable.
In numerical experiments, CEP calculated by mixed algorithm is close to theoretical value in the case of the population mixing two normal distributions, while traditional CEP method is unreasonable and unprecise, which illustrates that the mixed CEP algorithm designed in this paper is adaptable under the circumstance.In addition, the response relationship of prior parameters and CEP is analyzed in the case, which verifies Proposition 1 from the perspective of statistical simulation.And on the condition of small samples with mixture population of three normal distributions, the confidence upper bounds of mixed CEP with real and biased prior parameters are calculated, which indicates that the mixed CEP algorithm is credible for small samples.
In the case of impact samples following the mixed population of finite normal distributions, the mixed CEP algorithm proposed in this paper is credible.However, the effect of the method to other types of mixture models is not explicit, which is what we aim to investigate furthermore.
the prior probabilities of test experiments are fixed as [(ℎ  )] =   (0 ⩽   ⩽ 1) and populations (  , Σ  ) are independent, then the mean and variance estimates of mixed population (, ) will be

Figure 1 :
Figure 1: Probability density function of two normal mixture populations.

Figure 2 :
Figure 2: Sample dispersion and comparison of CEP by different algorithms.

Table 1 :
Contrast average CEP mixed with classical algorithms by resampling.

Table 2 :
CEP and corresponding errors with different prior information.

Table 3 :
Value of impact errors.

Table 4 :
Actual and prior mean of probability of test conditions.

Table 5 :
CEP under different assumptions of the population.