Maneuvering Detection Using Multiple Parallel CUSUM Detector

The switching model tracking algorithm based on hard decisions is an important method to solve the maneuvering target tracking problem.The use of Doppler velocity not only helps shorten the delay time of maneuvering detection but also provides information about the target motion model. A novel target maneuvering detection method named Multiple Parallel Cumulative Sum (MCUSUM) for target multiple maneuvering models is proposed in this paper based on Doppler velocity. The main scheme of the proposed approach consists of the following: firstly, the problem framework of multiple model maneuvering detection is put forward; secondly, the statistic of acceleration is obtained through modeling the mapping relationship between Doppler velocity and the normal/tangential acceleration according to the geometry and kinematics; thirdly, the joint empirical distribution of the normal/tangential acceleration is obtained by the statistical experiment method and then the approximate joint probability distribution function of the normal/tangential acceleration is acquired by use of GaussianMixtureModel (GMM)with Expectation Maximization (EM) algorithm; fourthly, it is taken as the prior information of target maneuvering which is introduced to the likelihood ratio of prediction measurement residual by the marginalization method; finally, the standard Cumulative Sum (CUSUM) detector is extended as Multiple Parallel CUSUM detector. Simulation results show that M-CUSUM detector has a smaller maneuver onset detection delay time compared with similar detectors and has the ability of pattern recognition of target maneuvers.


Introduction
Tracking the manned maneuvering aircraft by use of radar is a maneuvering target tracking (MTT) problem [1][2][3].It is hard to accurately estimate the motion state of maneuvering target by using only single motion model, because the current radar technology cannot obtain the observation about accelerations directly supporting the selection of the target motion model.The switching model tracking (SMT) algorithm [4] based on (the) hard decision and the interacting multiple model (IMM) tracking algorithm [5][6][7] based on the soft decision are two kinds of basic methods to solve MTT.The SMT algorithms use a maneuvering detector to judge the target maneuvering onset.When the target maneuvering is successfully detected, the motion model of tracking algorithm is immediately switched.The timely and correct maneuvering detection is the key to the SMT algorithm.The most maneuver detectors [8] only use the position measurement conditioned on the linear Gauss hypothesis.As the position measurement cannot directly reflect the maneuver state, however, radar can obtain not only the target position measurement but also Doppler velocity measurement directly affected by the maneuvering [9].The maneuvering state of target can be obtained by use of Doppler velocity.Some scholars obtain the turn rate estimation from Doppler velocity [10][11][12].As the acceleration is directly related to the target maneuvering type, other scholars obtain the statistics of the target acceleration directly through Doppler velocity.Bizup and Brown [13] propose a method of maneuvering target detection using Doppler velocity.By assuming the constant turn rate motion, a statistic  min of acceleration is deduced from Doppler velocity.Ru et al. [14] combine the statistic  min with the tangential acceleration to obtain the statistic of total acceleration,  min2 , and  min3 .The main difference between  min2 and  min3 is whether Doppler velocity measurement is applied in the filtering.It turns out that the statistic  min2 or  min3 using Cumulative Sum (CUSUM) detector has a shorter detection onset delay time than  min [15].Lu et al. [16] use Mahalanobis distance and Euclidean distance optimization method to give two statistics of the total acceleration  min and  min , respectively, which overcome the problem that the acceleration has no solution when the noise is higher.Using the statistics in the above references to detect and track the target, there are the following problems.
Assuming a Constant Turn Rate Motion.The statistics  min ,  min2 , and  min3 are obtained by assuming a constant turn rate motion.However, in addition to the normal acceleration, the target maneuver usually has a tangential acceleration and the trajectory exhibits a curvilinear shape.There is an inherent system bias when the statistics  min ,  min2 , and  min3 are used in the maneuver target detection.

One-to-Many Mapping Relationship.
There is a one-to-many mapping relationship between Doppler velocity and the normal acceleration.It is assumed that the target velocity direction change is less than 2 in the radar scanning period, so the one-to-many mapping is restricted to a pair of four mappings in the statistics  min ,  min2 , and  min3 .In practice, the tracking radar scanning period is less than 1 s and the maximum instantaneous turn rate of the manned maneuvering aircraft is not more than /4 [17].Therefore, if one-to-four mapping is applied to track such target, it will cause the dispersion of probability weight and affect the sensitivity of detector as some unlikely turns are considered.
Unknown Motion Model.The detector of total acceleration based on Doppler velocity can detect target maneuver rapidly and distinguish the maneuvering and nonmaneuvering state.After the maneuvering is detected, the motion model with larger acceleration noise is usually chosen to filter.However, the motion model mismatch problem cannot be completely solved by increasing the acceleration noise.If a maneuver detector can not only detect whether the target is making a maneuvering or not but also recognize the target maneuvering model, we could select the matched motion model to improve the tracking accuracy with a lower acceleration noise.
To solve the above problems, the mapping relationship between Doppler velocity and the normal/tangential acceleration is deduced without assuming a constant turn rate motion.Meanwhile, it is assumed that the target velocity direction change is less than  during the radar scanning period, and one-to-four mapping relationship is restricted to one-to-two which reduces the dispersion degree of probability weight.In this paper, two parameters, namely, normal/tangential acceleration, are used to describe the target maneuver, and the single total acceleration is not adopted.The joint empirical distribution of the normal/tangential acceleration is obtained by the statistical experiment method [15] and approximated by use of Gaussian Mixture Model (GMM) with Expectation Maximization (EM) algorithm [18].Then the approximate joint probability distribution function (PDF) of the normal/tangential accelerations is obtained.The approximate joint PDF is taken as the a priori information of target maneuvering which is introduced to calculate the likelihood ratio of prediction measurement residual by the marginalization method.A Multiple Parallel Cumulative Sum (M-CUSUM) detector is proposed by extending the standard CUSUM detector.In addition, it has been proposed that the use of Doppler velocity measurement can improve tracking accuracy [19][20][21].In this paper, Doppler velocity is used separately in the maneuvering detector and filtering.As the measurement equation with Doppler velocity is highly nonlinear, the traditional EKF algorithm [22] does not work well.We use the measurement conversion algorithm proposed by Jifeng and Huimin to construct the linear KF filter.
The content of this paper is designed as follows.We describe the tracking and detection problem and propose a concept framework of multiple maneuver model detection and recognition in Section 2. We present the method of using Doppler velocity to deduce the normal/tangential acceleration in Section 3. We present how to obtain the empirical distribution of the normal/tangential acceleration and the approximated GMM with EM algorithm in Section 4. We describe M-CUSUM detector for multiple maneuver model detection and recognition in Section 5.In Section 6, firstly the results of the empirical distribution of the normal/tangential acceleration and GMM fitting are shown in five scenes; secondly, M-CUSUM proposed in this paper is compared with other maneuver detection methods; finally, it is applied to the SMT algorithm and is compared with IMM algorithm and CSAF algorithm.We summarize the research content of this paper in Section 7. In Appendix A, the modified measurement conversion algorithm and filter algorithm are given.In Appendix B, The EM algorithm is given to learn the parameters in GMM.In Appendix C, the state space equation of CSAF algorithm is given in two-dimensional space.

Target Maneuver Model and Measurement
Model.An aircraft is moving in a two-dimensional plane whose acceleration can be decomposed into a normal acceleration perpendicular to the velocity direction and a tangential acceleration along the velocity direction.The normal/tangential accelerations are used as the input vector of the motion model and then the target maneuver model [18] is expressed as where the state vector is F and Γ(X  ) are as follows: where  is the radar scanning period;   is the target velocity direction.
Assume that the radar is located at the origin of coordinate .The measurement equation of radar which includes Doppler velocity is where Z  = [  ,   , V , ]  is the observation vector;   is the range;   is the azimuth; V , is Doppler velocity; V  is a white Gaussian noise and its covariance matrix is where "±" is "+" if   ≥ 0 and "±" is "−" if   < 0. ( 0 ) The target does not maneuver; A  ∈ {A | A = 0} for  = 1, . . ., , where  is the radar sample time.

The Problem Framework of Multiple Maneuvering
( 1 ) The target starts to maneuver at unknown time ; For the above maneuvering detection problem, the traditional maneuver detector can accomplish the judgment between the maneuver and nonmaneuver, that is to say, to solve a simple binary hypothesis [23] test task.However, decision-makers often want to get the information of the target maneuver type.Then hypothesis  1 can be decomposed into multiple branches and the maneuver detection problem can be described as a new form: ( 0 ) The target does not maneuver; A  ∈ {A | A = 0} for  = 1, . . ., .
A  represents the th maneuvering model;  is the total number of maneuvering target models.In particular, the nonmaneuver

The Normal/Tangential Acceleration
The key of reasonable multiple maneuvering detector is to obtain the information of normal/tangential acceleration.In the aspect of radar measurements, Doppler velocity contains the target maneuvering information.And [13] establishes the mapping relationship between Doppler velocity and the normal acceleration on the assumption of a constant turn rate motion.This paper renews to build a mapping relationship between Doppler velocity and the normal/tangential acceleration without the assumption.
We assume that the target moves in 2D plane of coordinate , as shown in Figure 1.In Figure 1,  is the angle between the target velocity direction and range extending direction; V is the size of target velocity.The basic relationship can be obtained: By formula (6), formula (7) can be obtained: As  ∈ [−, ], Assume that the target velocity direction denoted as  −1 is known at time  − 1.After the radar measurement (  ,   , V , ) coming at time ,   may be as follows: The variable Δ  ∈ [−,] is denoted as the change of target velocity direction during  in the following formula: If the direction from  −1 to   is the clockwise direction, denote Δ  < 0; otherwise, denote Δ  > 0. In the radar scanning period , there are an infinite number of ways from  −1 to   .It is assumed that the change of target velocity direction during  is less than 2 in [13].Under this assumption, there are four turn ways from  −1 to   .The maximum instantaneous turn rate of the manned maneuvering aircraft is less than /4, such as F-22.And the scanning period of the common Doppler radar is less than 1 s under the tracking pattern.Then, we assume that the change of target velocity direction during  is less than .Under this assumption, there are two turn ways from  −1 to   , which are the turn ways from  −1 to  (1)   with Δ  ( (1)   ,  −1 ) ≤  and from  −1 to  (2)   with Δ  ( (2)   ,  −1 ) ≤  as shown in Figure 2. Thus, the possible turn ways may only be "1" and "2," excluding "3" and "4" in Figure 2.
The turn rate denoted as   meets the left hand system: Then,   is obtained by formulas ( 9), (10), and (11): The normal acceleration can be calculated by This implies an assumption that   and   are constant in the period .As  is very short, this assumption is satisfied in practice.Then,  , is denoted as the following formula: The variable  , is unknown at time .Generally, the tangential acceleration of aircraft is provided by its engine thrust and the instantaneous change of thrust is slow.In a short period 2, we can assume  , =  ,−1 .Then by formula (7), we can calculate If the exact values of and  −1 are known, we could calculate  , and  , by formulas ( 13) and ( 14).Due to the random noise, it is assumed that the radar measurements V , , V ,−1 , V ,−2 ,   ,  −1 ,  −2 are random variables with the Gauss distribution.In addition to the above radar measurements, other variables are unknown which are considered as random variables obeying some probability distributions.All variables in formulas ( 13) and ( 14) are treated as random variables.Then  , and  , are functions of random variables.The joint PDF ( , ,  , ) describes the target maneuvering.However, there is a high degree of nonlinear relationship between  , ( , ) and its arguments; thus PDFs of its arguments cannot directly pass to  , ( , ).Here, we use the statistical simulation experiment to obtain the empirical distribution of ( , ,  , ).

Empirical Distribution and Its Approximated PDF
The vector (  ,   ) has uncountable combinations that correspond to maneuvering models.Obviously, it is convenient for us to study maneuvering detection by determining some typical maneuvering models through some combinations.If we select  typical maneuvering models, we should build  joint PDFs denoted as {  (  ,   )}  =1 whose empirical distribution functions are obtained through statistical simulation experiments.In the experiment, the standard 2D curve motion model is adopted to obtain the true trajectory of the target as the following formula: where   ,   ,  V , and   are the independent Gauss white noise;    ,    ,   V , and    are the corresponding standard deviation.In the experiment, multiple real trajectories given   and   could be obtained by formula (15) and by initializing different , , V, .The measurement vector (, , V  ) from Doppler radar is obtained by the following formula: where "±" is "+" if  ≥ 0 and "±" is "−" if  < 0; "∓" is opposite to "±";   ,   , and  V  are the independent Gauss white noise;    ,    , and   V  are the corresponding standard deviation.The radar measurement sequence can be obtained through inputting a real trajectory into formula (16).In addition to radar measurement, the value of V −1 ,  −1 and  −2 should be known to calculate  , and  , by formulas ( 13) and (14).We use a filter to estimate them.In order to solve the problem of the high nonlinearity of measurement equation ( 4), [14] proposes a measurement conversion algorithm that converts the measurement in Polar coordinate system to the pseudo-Gauss noise measurement in Cartesian coordinates system and then used KF for state estimation.For the better application in this paper, the filter has been modified as shown in Appendix A. X−1 = [x −1 , V,−1 , ŷ−1 , V,−1 ]  can be obtained by filtering and V−1 , α−1 , and α−2 can be calculated by the following formulas: The variables â, and â, can be calculated with those values by formulas ( 13) and (14).For the variable â, , there exist two possible values, which are â(1) , and â(2) , .In [13,14], the statistic  min = min{ (1)  , ,  (2)  , } is employed.However, when | − | is small, â(1) and â(2)  are similar.They will be taken approximately with the same probability as an estimate of   .Then  min will introduce a statistical error, since the smaller value between â(1) and â(2)  is selected.A new statistic   is proposed in this paper based on the statistic  min ; that is, where rand{â (1)  , , â(2) , } is to select a value between â(1) , and â(2) , randomly;  ≈ /18 is the experiential threshold value.We apply the Monte Carlo simulation method under every initial state given a maneuvering model (  ,   ).In each simulation, the target firstly moves at constant velocity in accordance with the initial state and then moves in accordance with the maneuvering model (  ,   ) after time   .Record   and   at time   and the empirical distribution denoted as (  ,   ) could be acquired.We denote the approximate joint PDF   (  ,   ) of (  ,   ) which can be obtained by fitting (  ,   ).This is a fitting problem of multidimensional nonlinear PDF.GMM can be used to approximate an arbitrary joint PDF.Therefore, GMM is used to approximate (  ,   ) and EM algorithm is applied to learn the parameters of GMM as shown in Appendix B.

M-CUSUM Detector
The purpose of this paper is to design a multiple maneuvering model detector that can deal with the compound hypothesis test problem described in Section 2.2.The maneuvering detector can be divided into two groups: one is a batch sliding window detector and the other is a sequential detector.Since the radar measurement is coming in order, the sequential detection is more suitable for the maneuvering detection.CUSUM as a kind of sequential detector can ensure the fastest detection at a certain decision error rate under the condition that the input vector A is known.For a simple binary hypothesis testing of the conventional maneuvering detection described in Section 2.2, the standard CUSUM detector is Decision rules are as follows: (1) Accept  1 , if   ≥ , and the maneuvering time m = min{ :   ≥ }.(2) Continue to detect and  =  + 1, if   < .Z 1: ≜ (Z 1 , . . ., Z  );  is the detection threshold (decision error rate); Z is the measurement prediction residual at time .The Multiple Parallel CUSUM detector, called M-CUSUM detector, is proposed in the paper for the detection and recognition of  maneuver patterns.The specific form is and  detection thresholds are denoted as { 1 ,  By formula (20), it is the key of the detection and recognition of  maneuver patterns to calculate ( Z |  0 , Z 1:−1 ) and ( Z |   , Z 1:−1 ).Assume that ( Z |  0 , Z 1:−1 ) is a zero mean Gaussian distribution in general: where S  is the measurement prediction residual error covariance at time .( Z |   , Z −1 ) is the measurement prediction residual likelihood function and is the Gaussian distribution with the mean H(Γ( X−1 )A  ): The measurement conversion algorithm (see Appendix A) can approximately represent H(Γ( X−1 )A  ) as a linear form H  Γ( X−1 )A  .Although the real value of A  is unknown, the GMM approximation joint PDF   (A  ) of A  can be obtained by the statistical simulation experiment.Then   (A  ) can be used as a priori information of target maneuvering and we adopt the marginal method to calculate ( Z |   , Z −1 ): where   (A  ) is the matrix form of GMM;    is the mean vector; Σ   is the covariance matrix;    is the weight value.

Simulation Experiment
The main purpose of the simulation experiment is to verify M-CUSUM detector and compare it with other maneuvering detectors.The SMT algorithms using these detectors are also compared with state-of-the-art IMM algorithm and current statistical model adaptive filter (CSAF) algorithm [24,25].   = 2 m,    = 2 m,   V = 0.1 m/s, and    = 0.005 rad.The standard deviations in formula ( 16) are set as    = 50 m,    = 0.01 rad, and   V  = 1 m/s.The initial state of target ( 0 ,  0 , V 0 ,  0 ) is set as (40 km, 60 km, 300 m/s, /6).Doppler radar is located at the origin of coordinate  and its scanning period is set as  = 0.5 s.Design five scenarios for performance evaluation.The target moves at a constant velocity (that is the nonmaneuver model A 0 = (0, 0)) during 0 s∼20 s in each scenario.The target starts to maneuver at   = 21 s and keeps 20 s.In five scenarios, the maneuver models of target are designed as A 1 = (0, 10), A 2 = (0, −20), A 3 = (10, 0), A 4 = (10, 10), and A 5 = (5, 20), as shown in Figure 3.

Tracking Algorithm Settings.
The filters are divided into two types according to the different state equations adopted.The first filter is one whose state equation [26] does not include the input vector U  : Formula ( 24) is often taken as the state equation by the standard tracking algorithm where the common SMT algorithms could be compared.The state transform matrix F in formula (24) is the constant velocity model (CV) and we design two level noise covariance matrices Q low and Q high , respectively, as the nonmaneuvering model and the maneuvering model.The SMT algorithm switches the motion model between Q low and Q high based on the maneuvering detector.The model set of IMM algorithm consists of two CV motion models with Q low and Q high , respectively.In the simulation, we especially set Q low ( = 10) and Q high ( = 1000).
The second filter is one whose state equation includes the input vector U  , as described in the following formula: For the STM algorithm with M-CUSUM detector and IMM algorithm, the input vector U  is denoted as Γ(X  )A  .Then the state equation ( 25) is the same as the state equation (1).Given F and Q  , the different A  represents the different maneuvering model.If A  = A 0 , the state equation ( 25) describes the nonmaneuver model, that is, CV; if A  ∈ {A 1 , A 2 , A 3 , A 4 , A 5 }, the state equation ( 25) describes the maneuver models corresponding to the five scenarios.In the simulation, the switching models of M-CUSUM detector are denoted as {A  } 5 =1 and the model set of IMM algorithm is denoted as {A  } 5 =0 .The state transform matrix F in formula (25) also is CV and the noise covariance matrixes are set as Q  0 ( = 10) and Q   ( = 600),  = 1, 2, 3, 4, 5.
Reference [25] proposes a direct method of estimating the tangential and normal accelerations of maneuvering targets in three-dimensional space.The method takes the modified Rayleigh distribution as the statistical model of the tangential and normal accelerations that directly are introduced into the state vector X and updates the state noise covariance matrix Q  online to adapt to the current maneuver of target.The method is called the current statistical model and adaptive Kalman filter (CSAF) algorithm in the paper.In order to compare with M-CUSUM, CSAF algorithm should be deduced in two-dimensional space and modified with the measurement from Doppler radar.The discrete U  , F  , and Q  in the state equation ( 25) can be deduced through (C.5) given in Appendix C for the current statistical model, and the adaptive Kalman filter algorithm can refer to [25] and Appendix A.
Two types of filter are set up to provide a comparative platform for different tracking algorithms.The maneuver detector based on the binary hypothesis can only be used in the detection of nonmaneuver model and maneuver model but cannot give the specific information of maneuver model.Thus it can only be used in the first filter but not in the second filter.M-CUSUM is a multiple maneuver model detector, which can realize the hard decision among multiple maneuver models, and the IMM algorithm can realize the soft decision among multiple maneuver models.Both M-CUSUM and IMM can be used in the tracking problems with the first filter and second filter.CSAF algorithm can realize the adaptive modification of U  and Q  which can be compared with M-CUSUM and IMM with the second filter only.For two types of filters, the measurement conversion algorithm in Appendix A is used.The difference is that the different state equation is used.
The average onset detection delay time (n − ) is used to measure the performance of the maneuver detector.The root mean square error (RMSE) of state estimation is used to measure the performance of the tracking algorithm: where  is Monte Carlo simulation time;   is the real state vector; x is the state estimation vector.=1 are acquired by the statistical experiment.We should assign many values to the initial state ( 0 ,  0 , V 0 ,  0 ) for experiments, because (  ,   ) is affected by the velocity magnitude, velocity direction, and position.In each dimension, a number of typical levels are selected for uniform experiments, and the specific values of parameters are shown in Table 1.
In Table 1, there are 162 initial state values.We execute 30 Monte Carlo simulations for each maneuvering model in {A  } 5 =1 and each initial state value in Table 1.Record the value of (  ,   ) at time   = 21 s, and we can get 4860 groups of data for each maneuver model.Then the empirical distribution function   (  ,   ) corresponding to one maneuvering model A  can be obtained.Two-dimensional GMM by EM algorithm in Appendix B is used to fit the empirical distribution, and then the approximate PDF    (  ,   ) can be acquired.Specifically, three two-dimensional Gauss distribution functions are used to fit the distributions.The scatter diagrams and the approximate GMM diagrams of each maneuver model in {A  } 5 =1 are shown in Figure 4.And the parameters of GMM are recorded in Table 2. From Figure 4, different maneuvering models correspond to different distribution patterns, which proves that the probability distribution of variables in formulas ( 13) and ( 14) cannot be directly transmitted to   or   .It is also claimed that the normal acceleration is related to the tangential acceleration without the assumption of a constant turn rate motion.Thus, it will introduce an error if the joint PDF (  ,   ) is decomposed into the product of marginal PDF (  )(  ).

Detection Performance Comparison.
In the simulation experiment, the M-CUSUM detector is compared with the measurement residual (MR), Input Estimate (IE), CUSUM, and  min and  min2 detectors.In the calculation of the average delay time (n−), the false alarm probability is set as   = 0.01 and the miss probability is set as   = 0.1.For MR and IE detectors, the window size is set as 5.For  min detector, the threshold value of fading memory average (FMA) is set as 5000.For  min2 and CUSUM detectors, the threshold value can be calculated by  = log((1 −   )/  ) [27], that is,  = 1.9542.For comparison, set the threshold values of M-CUSUM detector as  1 =  2 = ⋅ ⋅ ⋅ =  5 = .Table 3 summarizes detection delay time of all above algorithms and the detectors have longer detection delay time for relatively small maneuver scenarios as A 1 and A 3 .
min ,  min2 , and M-CUSUM detectors taking advantage of Doppler velocity measurement are more sensitive to the target maneuver compared with MR, IE, and CUSUM detectors.But when maneuver models (A 4 and A 5 ) have both the normal acceleration and the tangential acceleration, M-CUSUM detector has a shorter detection delay time.It is the main reason why M-CUSUM detector uses a more accurate prior knowledge of the normal/tangential acceleration.(28)

Estimation Performance
We execute 100 Monte Carlo simulations for the five scenarios in the first filter.Table 4 records the position and velocity of the average RMSE with different tracking algorithms for each scene.The velocity and position of RMSE of different tracking algorithms for each scene are shown at the sampling time in Figure 5, respectively.For the maneuvering target, the sooner it switches the motion model, the more conducive it is to the state estimation for the SMT algorithm.It is concluded that the smaller delay time leads to the lower RMSE, which means the better tracking accuracy according to Tables 3 and 4. Figure 5 shows that M-CUSUM has the smallest RMSE in target maneuver phase compared with MR, IE, CUSUM,  min ,  min2 , and CUSUM.The difference among the average position RMSE of each detector is small.However, the large difference is reflected in the average velocity RMSE in Table 4.When the target maneuvering is stronger, the average velocity RMSE is higher.This shows that the degree of mismatch between the real maneuvering model and the filter motion model is more directly reflected in the velocity than in the position, which is the reason why the efficiency of the detectors based on Doppler velocity measurement is higher than that based on the position measurement.
The average position and velocity RMSE of IMM algorithm are higher than the SMT algorithms.The SMT algorithms using the matched target motion model Q low have lower RMSE during the nonmaneuver phase 0∼20 s as shown in Figure 5.As the model noise of IMM algorithm is larger than Q low , it has a higher RMSE during the nonmaneuver phase 0∼20 s.However, as IMM algorithm could adaptively adjust the weights of two models Q low and Q high to match the target maneuver model, it has a lower RMSE during the maneuver phase 21 s∼40 s.The difference in average RMSE between IMM algorithm and SMT algorithms mainly comes from the accumulation of nonmaneuver phase as shown in Table 4.The tracking algorithms based on M-CUSUM,  min , and  min2 detectors have the comparable estimation performance.Moreover, since M-CUSUM detector has the normal/tangential acceleration in the target maneuver, the corresponding algorithm obtains the better estimation accuracy.The CV motion model is more mismatched to the maneuver model with larger normal accelerations as A 2 and A 5 compared with other maneuver models as A 1 and A 3 as shown in Table 4 and Figure 5.
In the second filter, IMM algorithm adopts six models: {A  } 5 =0 .The initial state vector  0 and the initial estimation error covariance matrix P 0 are the same as the first filter.The initial model probability denoted as  0 and model transition probability denoted as   of IMM algorithm are set as  0 = [0.90.02 0.02 0.02 0.02 0.02] , .9 0.02 0.02 0.02 0.02 0.02 0.02 0.9 0.02 0.02 0.02 0.02 0.02 0.02 0.9 0.02 0.02 0.02 0.02 0.02 0.02 0.9 0.02 0.02 0.02 0.02 0.02 0.02 0.9 0.02 0.02 0.02 0.02 0.02 0.02 0.9 The parameters of CSAF algorithm are set as where  max  and  max − are the positive and negative maximum of tangential acceleration;  max  and  max − are the positive and negative maximum of normal acceleration;   and   are the maneuvering frequency of the tangential and normal acceleration, respectively.The initial state vector X 0 and its standard deviation of CSAF are set as (31) The initial estimation error covariance matrix of CSAF denoted as P 0 is set as (32) If the SMT algorithm with the M-CUSUM detector detects the target maneuver, it will adopt the maneuver model A  ,  = 1, 2, 3, 4, 5, which firstly reaches the detection threshold to replace the term A  = A 0 in state equation (1).From the perspective of position and velocity, IMM algorithm and SMT algorithm in the second filter have higher estimation accuracy than in the first filter according to Figure 6.The more important reason is that the motion models {A  } 5 =0 of state equation ( 1) are more suitable to the maneuver states of the target in the five scenes than the motion models Q low and Q high .The SMT algorithm with M-CUSUM detector is significantly better than the IMM algorithm in two strong maneuver scenarios as A 2 and A 5 as shown in Figure 6.
First of all, this is because the IMM algorithm has the competition problem among multiple motion models, which can cause the probability weights to be dispersed by the unmatched models, and the system error is introduced into the state fusion.
Secondly, M-CUSUM detector can achieve higher recognition rates in strong maneuver scenes compared with other scenes, as shown in Table 5.The recognition rate of A 2 and A 5 maneuver models based on M-CUSUM detector is 93% and 95%, respectively.

Recognition results
A Then, the target state correction effect of the input vector is obvious in the strong maneuver scene, because the magnitude of input is higher than the noise level.However, as the magnitudes of input vector and noise are similar for weak maneuvering scenes, the target state correction of input vector will be interfered with the noise, which can lead to the performance degradation of the SMT algorithm with M-CUSUM detector.
The variables U  and Q  are modified online in CSAF algorithm as the current statistical model is used.Then the performance of CSAF is similar to IMM.CSAF algorithm has bigger position and velocity RMSE than SMT algorithm with M-CUSUM detector in the nonmaneuver stage, as shown in Figure 6.This is because the normal acceleration estimation obtained by the filter directly is more sensitive to the position measurement error than the target maneuvering in the nonmaneuvering stage.The tangential and normal accelerations in the iterative process are mutually independent with the modified Rayleigh distribution in CSAF algorithm.However, it has been proven in Section 6.2.1 that they are correlative.The error due to the independence assumption is constantly accumulated and then CSAF algorithm has significant tracking error around 30 s∼40 s of the maneuvering scenarios A 4  and A 5 that have both the tangential and normal acceleration, as shown in Figures 6(g)-6(j).However, CSAF algorithm has the comparability with STM algorithm and IMM algorithm in the maneuvering scenarios A 1 , A 2 , and A 3 that have only the normal acceleration or the tangential acceleration.

Conclusions
A theoretical framework to support the detection and recognition of multiple maneuvering models is proposed in the paper.Multiple Parallel CUSUM detector named as M-CUSUM is proposed by extending the standard CUSUM detector.The normal and tangential accelerations based on Doppler velocity measurement are deduced without the assumption of a constant turn rate motion.A new statistic of normal acceleration is proposed to reduce the error based on the statistic  min .The joint PDF of normal and tangential accelerations is used to describe the target maneuvering model.Its joint empirical PDF is obtained by the statistical experiment method and approximated by use of GMM with EM algorithm.The computer simulation results show that M-CUSUM detector not only has a good detection performance but also can recognize the target maneuver model to increase the matching degree of motion model.According to V  = V  cos  + V  sin , we can obtain the measurement equation: where  =  − where The measurement conversion algorithm needs (V  , V  ) at the current time.However, as the measurement conversion algorithm is executed before the filter, we use the velocity estimate at the previous time in practice.Then, we can directly estimate the target state using KF on the linear Gauss hypothesis in this way.If state equation ( 1) is used, the specific filtering algorithm is described as follows: One step state prediction and error covariance prediction:

B. GMM with EM Algorithm
The approximate joint PDF   (  ,   ) of (  ,   ) is modeled as GMM which is two-dimensional Gauss distribution function:  (B.7) The Lagrange multiplier method is used to estimate parameter ψ+1 .ψ+1 is taken as the initial value for the next step.Repeat the above steps until convergence.Finally, we can obtain the parameters of the approximation distribution function   (  ,   ) by EM algorithm.

C. Acceleration Estimation with CS Motion Model
In order to estimate the normal acceleration   and tangential acceleration   ,   and   are introduced into the state vector X:

Figure 3 :
Figure 3: The real trajectories of the five scenarios.
Comparison.In the first filter, the estimation performances of SMT algorithms with MR, IE, CUSUM,  min ,  min2 , and M-CUSUM detectors and IMM algorithm are compared.These tracking algorithms all use Doppler velocity measurement in the state estimation.The initial state vector  0 and its standard deviation are set as  0 = (40 km, 245 m/s, 60 km, 150 m/s) ,  0 = 30 m,  V = 5 m/s,  0 = 30 m,  V = 5 m/s.(27)The initial estimation error covariance matrix denoted as P 0 , the initial model probability denoted as  0 , and model transition probability of IMM algorithm denoted as   are set as

Figure 4 :
Figure 4: Scatter and GMM diagrams corresponding to five maneuver models.
Velocity RMSE of A4

Figure 5 :
Figure 5: RMSE of the first filter.
Velocity RMSE of A4

Figure 6 :
Figure 6: RMSE of the second filter.

Table 2 :
The parameters of GMM.

Table 3 :
The average delay time of maneuver detection (s).

Table 4 :
The average position RMSE (m) and average velocity RMSE (m/s).

Table 5 :
Recognition results of 100 Monte Carlo simulations based on M-CUSUM.