A Cost-Effective Tracking Algorithm for Hypersonic Glide Vehicle Maneuver Based on Modified Aerodynamic Model

In order to defend the hypersonic glide vehicle (HGV), a cost-effective single-model tracking algorithm using Cubature Kalman filter (CKF) is proposed in this paper based on modified aerodynamic model (MAM) as process equation and radar measurement model as measurement equation. In the existing aerodynamic model, the two control variables attack angle and bank angle cannot be measured by the existing radar equipment and their control laws cannot be known by defenders. To establish the process equation, the MAM for HGV tracking is proposed by using additive white noise to model the rates of change of the two control variables. For the ease of comparison several multiple model algorithms based on CKF are presented, including interacting multiple model (IMM) algorithm, adaptive grid interacting multiple model (AGIMM) algorithm and hybrid grid multiple model (HGMM) algorithm. The performances of these algorithms are compared and analyzed according to the simulation results. The simulation results indicate that the proposed tracking algorithm based on modified aerodynamic model has the best tracking performance with the best accuracy and least computational cost among all tracking algorithms in this paper. The proposed algorithm is cost-effective for HGV tracking.


Introduction
Hypersonic glide vehicle (HGV) [1] is a rapid strike weapon which can make long-range glide in near space by aerodynamic force at greater than Mach 5.It is a great threat to world peace that HGV can even strike anyplace in the world within two hours.Therefore, some papers about HGV tracking are published for the research of defending this weapon.
The wide range of maneuvers in three-dimension makes it really difficult to track the HGV and researchers prefer multiple model algorithms to single model algorithms for tracking the so-called strong maneuvering targets.Though the study of HGV tracking is in primary stage, there are still some representative papers about tracking some simple but not actual HGV trajectory.Zhang et al. [2] study the HGV's Sanger trajectory [3] (without lateral maneuver) tracking using interacting multiple model (IMM) [4] algorithm based on constant acceleration (CA) model, constant turning (CT) [5] model and Singer model.Qin et al. [6] study the near space target tracking with adaptive grid interacting multiple model (AGIMM) [7] algorithm and the target trajectory is a simple combination of constant-velocity motion, constant-turning motion and constant-acceleration motion.Furthermore, the AGIMM with three different CT models is used to track the near space target, which is just called near space target but does not consist of any character of HGV.Zhai et al. [8] use adaptive IMM algorithm based on the dynamic model to study the tracking and predicting the HGV trajectory with the fixed known attack angle and the varying unknown bank angle.The attack angle and the bank angle are the two control variables.However, in real HGV trajectory tracking, the attack angle and the bank angle are unknown, which is different from all mentioned above.Thus, tracking the maneuvering trajectory with unknown attack angle and unknown bank angle will be studied in this paper using several algorithms.
If we track HGV just using some kinematic models, such as CA model, constant velocity (CV) model, CT model and so on, we may miss much important HGV information.To develop new HGV tracking algorithms using as much information as possible, in every model filter, we establish the process equation based on HGV aerodynamic model and the measurement equation according to radar measurement principle.However, the process equations of the two control variables are unknown by defenders.Hence, an effective method is needed to be proposed for modeling the rates of change of the two control variables.Furthermore, if we do not know how the two control variables change from the discrete time k-1 to k, we can also use multiple model tracking algorithms, including IMM algorithm, AGIMM algorithm and HGMM [9] algorithm.We will use different algorithms mentioned above to track the maneuvering trajectory with unknown attack angle and bank angle, and compare the accuracy and computational cost in the whole tracking stage.
Cubature Kalman filter (CKF) [10] is proposed by Arasaratnam based on the Cubature rules with rigorous mathematical derivation.CKF is more accurate than extended Kalman filter (EKF) [11] in nonlinear filter and more accurate than unscented Kalman filter (UKF) [12] in high-degree nonlinear system filter [13].Furthermore, the system of tracking HGV is a high-degree nonlinear system.According to the above considerations, every model filter in algorithms in this paper can be constructed using CKF.
Radar is measurement equipment widely used in object detection and target tracking.In this paper, we conduct HGV tracking using only one sensor namely radar and in the future we could conduct HGV tracking using multiple sensors because of more measurement information and higher accuracy.If we used multiple sensors to measure HGV, the fuzzy multiple-sensor data fusion Kalman filter approach [14] proposed by James A. Rodger can improve the tracking performance effectively, which will be studied in the future.
This paper is organized as follows.In Section 2, the tracking algorithm based on modified aerodynamic model (MAM) is proposed, including the HGV aerodynamic model, radar measurement model and the MAM tracking algorithm.In Section 3, the multiple model tracking algorithms are proposed, including IMM tracking algorithm, AGIMM tracking algorithm and HGMM tracking algorithm.In Section 4, two typical HGV trajectories are generated and tracked using different algorithms by simulation, and the tracking performances are compared and discussed systematically.Finally, some valuable conclusions are drawn in Section 5.

HGV Tracking with Modified Aerodynamic Model
From above consideration, we know CKF is a good choice for high-degree nonlinear system, just like HGV tracking system.Based on CKF theory, the tracking problem of a nonlinear dynamic system can be defined by the following state-space model with additive noise in discrete time [15].
Process equation: Measurement equation: where X k is the state variable at discrete time k; and {w k } and ν R,k are independent process and measurement Gaussian noise sequence with zero means and covariance matrixes Q k and R k , respectively.Furthermore, the process equation will be established based on HGV aerodynamic model and the measurement equation will be obtained based on radar principle.

HGV Aerodynamic Model
Assuming the earth is a homogeneous sphere with rotation [16] and sideslip angle is zero, the three-dimensional point-mass dynamics of HGV in semi-speed coordinate system are described by the following equations of motion [17]: where V, θ, σ, r, λ and φ are the state variables, which stands for speed, flight path angle, velocity azimuth angle measured from the North in a clockwise direction, radial distance from the center of the earth to the vehicle, longitude and latitude, respectively.Bank angle ν and attack angle α are the two unknown control variables, which will be estimated in the following sections for tracking HGV.g is gravity acceleration.t is the time.ω e is a coefficient related to the earth's rotation.L and D are lift acceleration and drag acceleration, respectively, which are defined as following: where q = 1 2 ρV 2 is dynamic pressure.S M is the HGV reference area.m is the HGV mass.ρ is mass density of air.HGV is designed based on the high-lift common aerodynamic vehicles (CAV-H) [18] so we can study the HGV aerodynamic model using the CAV-H aerodynamic parameters: m = 908 kg, S M = 0.484 m 2 .Based on the U.S. Standard Atmosphere, 1976(USSA1976) [19], the mass density of air denoted as ρ and the sound speed a can be functions of height as ρ = 1.225 × 10 −0.1133h (5) where h is height in kilometers.
The drag and lift coefficient C D and C L are functions of attack angle α and Mach number Ma, which can be expressed as: where C D0 , C D1 , C D2 , C D3 , C L0 , C L1 , C L2 and C L3 are some fixed numerical coefficients.Thus, we have almost established aerodynamic model.Furthermore, there are some important parameter constraints that must be considered, such as dynamic pressure q, aerodynamic load n and heat flux .Q.These constraints can be expressed as where q max , n max and .
Q max are the maximum dynamic pressure, maximum aerodynamic load and maximum heat flux.ρ 0 is the mass density of air at sea level.K is a coefficient associated with the head of HGV.
When the bank angle ν is not zero, the trajectory will maneuver laterally.Usually, in the HGV controlling, it is common to keep the attack angle fixed and keep the bank angle varying to make maneuvering laterally and skipping longitudinally trajectory.The bank angle is the main control variable and the attack angle is the auxiliary control variable for HGV controlling.Furthermore, in HGV trajectory design, the nominal attack angle is set based on engineering experience and the bank angle is used as the main control variable within the constraints of air heating, loading and dynamic pressure [20][21][22][23].There are three main reasons for using above control approach [24]: First, the attack angle determines the flight distance.Furthermore, once the flight mission is determined, the attack angle is determined correspondingly.
Second, the attack angle has a direct effect on aerodynamic force.The adjustments to attack angle frequently make aerodynamic force complicated, resulting in aggravating the burden of control system.
Third, the proper thermodynamic environment is provided by using big angle for HGV flight.It is harmful to the thermal protection system when adjusting attack angle frequently.
However, in some papers [17,25] the attack angle is not a constant and the two control variables can change at the same time to a large degree.
However, we will study HGV tracking using the typical trajectory with the control variables changing at the same time in this paper, which will be more convincing to validate the tracking algorithm.
Through analysis of Equations ( 3)-( 7), we can know that the bank angle and attack angle have an important effect on flight distance.Furthermore, for an appropriate flight distance the variation ranges are usually set as follows [8]: attack angle, 6 • to 12 • ; bank angle, −20 • to 20 • .
Taking the variation ranges of attack angle and bank angle into consideration, the control variables do not change at a big degree.From the aspect of kinematic models to study the trajectory HGV maneuvers strongly, but from the aspect of aerodynamic model to study the trajectory HGV does not maneuver that strongly in the long flight distance.

Radar Measurement Model
Assuming that, in the spherical coordinate system, the coordinates of the origin of the radar measurement coordinate system are [ r 0 λ 0 φ 0 ] T , the coordinates of HGV centroid are T in the spherical coordinate system and x y z T in the radar measurement coordinate system.According to coordinate transforming relations, the transforming equations can be denoted as where F is the transform matrix of the radar measurement coordinate system to the geocentric rectangular coordinate system and can be denoted as where R Z , R X and R Y is the rotation matrixes around the Z-axis, X-axis and Y-axis, respectively.
Assuming the HGV coordinate in radar measurement coordinate system is x y z T , based on the radar principle [26] RAE can be expressed as: In the HGV tracking, we assume the radar measurement noises are known independent white noises, so the radar measurement equations can be expressed as T , then the radar measurement equation can be expressed as The statistical characteristics of radar measurement noises v R can be expressed as where σ 2 R , σ 2 A and σ 2 E are measurement error deviations of range, azimuth and elevate, respectively.

Modified Aerodynamic Model Tracking Algorithm
It is known by us that the existing radar equipment cannot directly measure the attack angle and bank angle of HGV.Furthermore, the change law of the two control variables is unknown by defenders as an important secret of a strike weapon.Thus, we need to find an effective method for modeling the change law of the two control variables.Because HGV tracking is still in primary stage, there is not an existing method for modeling the change law.
To develop the MAM algorithm, we analyze the characters of HGV in near space first.HGV maneuvers by aerodynamic force with a Mach number bigger than 5 in near space and HGV can bear damage to the vehicle structure when bank angle and attack angle change strongly.Furthermore, in the practical control of HGV, the control variables usually change continuously in a certain degree for control stability.Thus, we think the rates of change of control variables in a short time are not very large and unknown by defenders.
Then, the two control variables can be modeled as constant Gaussian Markov process and we can use additive white noise to model the rates of change of attack angle and bank angle [27] as where w t,α and w t,ν are white noise with an certain effect on attack angle and bank angle that accounts for the unknown control input.
Thus, the complete process equation can be established.The augmented state variables at discrete time k can be defined as The basic process equations can be written by making Equation (3) discretization as Furthermore, the augmented equations about the control variables can be written by make Equation ( 15) discretization as where T is the time interval.w V,k , w θ,k , w σ,k , w φ,k , w λ,k , w r,k , w α,k and w ν,k are the components of w k .We consider that the changes of the attack angle and the bank angle at the discrete time k and (k−1) are caused by the time interval T multiplied by white noise w k−1,α and w k−1,ν , respectively.The method to model the rates of change of the two control variables is similar to the method in CV model to model the rate of change of velocity [28].In CV model the ideal equation is modified as .
where the white noise w (t) has a "small" effect on x that accounts for unpredictable modeling errors due to turbulence, etc.After getting the process and measurement equations, according to CKF, we can get the predicated state X k at discrete time k using the state X k−1 at discrete time k−1 and the new measurement Y k .CKF operates only on means and covariances of the conditional densities encountered in time and measurement updates [10].To develop this filter, a third-degree spherical-radial cubature rule is used to compute the various multi-dimensional Gaussian-weighted moment integrals.CKF provides an efficient solution even for high-dimensional nonlinear filtering problems.

HGV Tracking with Multiple Model Algorithms
To establish the process equation based on unmodified aerodynamic model with unknown attack angle and unknown bank angle, we can use multiple model algorithms.In the following section the three-model multiple model algorithms are established based on aerodynamic model with attack angle unmodified and bank angle modified, and the nine-model IMM algorithm is established based on unmodified aerodynamic model with both attack angle and bank angle unmodified.
Multiple model algorithms for target tracking include fixed structure multiple model (FSMM) algorithm and variable structure multiple model (VSMM) [29] algorithm.When we apply the FSMM to a complicated real system estimation, sometimes we encounter two problems [9]: First, the chosen model set may not cover the full range of the mode and the truth may lie between adjacent models.Second, even if the chosen model set is large enough to cover the full range, use of all those models does not guarantee performance improvement, not to mention the prohibitively large computational cost.To overcome those defects of FSMM estimation and increase cost-effectiveness, VSMM is proposed in [29][30][31].IMM algorithm is one of the most popular and most cost-effective FSMM algorithms [32].AGIMM algorithm and HGMM algorithm are two typical VSMM algorithms.Tracking algorithms using IMM algorithm, AGIMM algorithm and HGMM algorithm will be briefly showed in this section as follows.

IMM Tracking Algorithm
In this section we establish the HGV IMM tracking algorithms with three-model IMM algorithm and nine-model IMM algorithm based on aerodynamic equations and radar measurement equations.In the three-model IMM algorithm the distinguishing aerodynamic model parameters are combination of three known fixed attack angles and corresponding unknown varying bank angles.The three fixed attack angles are 6 • , 9 • and 12 • .The process equation of corresponding unknown varying bank angle can be set as In the nine-model IMM algorithm, the distinguishing aerodynamic model parameters are combination of three known fixed attack angles and three known fixed bank angles.The three known fixed bank angles are −20 • , 0 • and 20 • .The transition between the N models is governed by a Markov chain, named as the transition probability ρ ij , which is selected at the beginning of the algorithm.
The IMM algorithm can be divided into four parts namely interacting, Kalman filtering (model updating), model probability calculation and estimate combination.For ease of comparison with other algorithms and considering above discussion about different nonlinear Kalman filters, we use CKF in model updating.Then, the IMM algorithm can be demonstrated as follows:

Interaction of the Estimates
The mixed state estimate for the CKF model M j at discrete time k−1 is calculated using outputs of the models X i k−1/k−1 , the corresponding model probability ω i k−1 , the transition probability ρ ij and is given by where N is the model number.The error covariance of the interactive state for the filter model M j at discrete time k−1 is where P i k−1/k−1 is the estimate covariance for the i-th filter at discrete time k−1.

Model Updating
We set X

Model Probability Calculation
The likelihood of the CKF model is given by the innovation Y j k and the covariance of the innovation where M is the dimension of the radar measurement; Y j k is the residual of actual measurement and prediction measurement in CKF; and The model probabilities ω j k are given by Bayes rule from

Estimate Combination
The final output is generated by combining the state estimates from each model as

AGIMM Tracking Algorithm
In order to solve real-world problems the use of only a small number of models is not good enough but the use of more models does not necessarily improve the performance and increases the computational burden considerably.Thus, the concept of variable structure multiple model estimation is proposed and several VSMM algorithms are developed by researchers.An effective VSMM algorithm is to make adaptive the grid of the parameters that characterize the possible modes.
In AGIMM tracking algorithm a coarse grid relating about attack angle is set up initially, and then the grid is adjusted recursively.HGV tracking using AGIMM algorithm can be demonstrated as follows: Because the real attack angle (unknown) is assumed to be within the continuous range [6,12], we can consider a three-model IMM algorithm with time-varying supporting digraph D k and the model set α 1,k , α 2,k , α 3,k for grid values (vertices of Initially, the AGIMM is started with the coarse grid At each recursive time-step ( k → k + 1 ), the grid is adjusted as following two-step adaptation logic.

Grid Center Readjustment
In the AGIMM algorithm, α 2,k is assumed to be the grid center and can be obtained by

Grid Distances Readjustment
In this paper, we divide the grid distances readjustment situation into three cases no-jump, up-jump and down-jump.Case 1. No-jump: where 35 is a threshold for detecting an unlikely model; and ξ = 1 is a model separation distance (design parameters).Case 2. Up-jump: where t 2 = 0.65 is a threshold for detecting the significant mode.Case 3. Down-jump: Thus, the AGIMM algorithm for HGV tracking is defined.Through the adaptive grid method the attack angle α 1,k , α 2,k and α 3,k are closer to the real unknown attack angle which is good to improve the tracking accuracy.

HGMM Tracking Algorithm
Stimulated by the idea that adding a new model set which is close to the true model will improve the performance of multiple-model estimation [33], the HGMM estimator is developed using the hybrid grid to cover a large continuous mode space with a relatively small number of models at a given accuracy level.
If the true mode is known and added to the model set, the estimator will converge to the true mode with probability 1, and the optimal state estimate can be obtained [34].However, the true model is unknown actually, and a common way is to add the optimal mode estimate, which is statistically close to the true mode [35].
In HGMM algorithm the model set about attack angle is a combination of a fixed coarse grid and an adaptive fine grid.Denote by A k the adaptive fine grid and by M the fixed coarse subset.One cycle of HGMM algorithm consists of three steps: first, obtain the state estimate X based on M; second, obtain the fine subset using the state estimate X; and, third, run VSMM.
Then the HGMM tracking algorithm can be demonstrated as follows:

Determination of the Coarse Subset
The fixed coarse grid is set as the same as that in AGIMM algorithm and can be written as

Determination of the Fine Subset
Two classes of approaches are presented in [36], which are simpler and easier to implement.The two approaches are: fine grid with fixed grid distance, and fine grid with adaptive grid distance.We use the first approach in this paper.
We need to get an attack angle estimate αk based on M using IMM algorithm and then we can form the fine grid by quantizing the region whose center is αk as following: where d is the fixed grid distance, and in this paper we set d = 1 • .

Simulation and Discussion
In this section, the HGV maneuvering trajectory is generated in semi-speed coordinate system and is to be measured and tracked in radar measurement coordinate system.The simulation configuration is defined in Section 4.1.The simulation results using MAM algorithm, three-model IMM algorithm, nine-model IMM algorithm, AGIMM algorithm and HGMM algorithm are compared and analyzed in Section 4.2.

Simulation Configuration
We set the initial augmented state of the HGV as follows: V = 6000 m/s; θ = −0.1 rad; σ = 0 rad; r = 6.378 × 10 6 + 50 × 10 3 m; λ = 0 • ; φ = 0 • .To design a trajectory in a global strike mission, some constraints of real flight must be considered such as waypoints for reconnaissance, no-fly zones for threat avoidance, maximum dynamic pressure, maximum aerodynamic load and maximum heat flux.However, in this paper, we mainly study the HGV tracking algorithms and focus on the tracking performance of different algorithms.Thus, it is proper to design some simple but typical trajectories to test the tracking algorithms regardless of above constraints.
To validate the effectiveness of the HGV maneuvering trajectory tracking algorithms proposed in this paper, two maneuvering methods are given from where V 1 = 5000 m/s, V 2 = 4500 m/s, α max = 12 • and α min = 6 • ; t denotes the flight time.Equation (37) shows the HGV maneuvers laterally at the biggest degree.Furthermore, in the real world, controlling the bank angle cannot jump 20 • immediately.The change law of Equation ( 37) is to challenge the tracking algorithm in simulation.Equation (38) shows the HGV maneuvers laterally by the bank angle changing in the sine law with period 80 s and amplitude 20 • .Then, two HGV trajectory scenarios can be generated as Figure 1.
where 1 V = 5000 m/s, 2 V = 4500 m/s, max α = 12  and min α = 6  ; t denotes the flight time.Equation (37) shows the HGV maneuvers laterally at the biggest degree.Furthermore, in the real world, controlling the bank angle cannot jump 20° immediately.The change law of Equation ( 37) is to challenge the tracking algorithm in simulation.Equation (38) shows the HGV maneuvers laterally by the bank angle changing in the sine law with period 80 s and amplitude 20°.Then, two HGV trajectory scenarios can be generated as Figure 1.In Figure 1, Scenario 1 is generated based on Equation ( 37) and Scenario 2 is generated based on Equation (38).From Figure 1 we can know that HGV maneuvers laterally more strongly when bank angle changes according to Equation (37).Assume the radar is located at the position: In Figure 1, Scenario 1 is generated based on Equation (37) and Scenario 2 is generated based on Equation (38).From Figure 1 we can know that HGV maneuvers laterally more strongly when bank angle changes according to Equation (37).Assume the radar is located at the position: r = 6.378 × 10 6 + 200 m; λ = −π/32 rad; φ = π/16 rad.The standard deviation of the radar measurement errors is set as σ R σ A σ E = 10 m 0.5 mrad 0.5 mrad (39) Every process noise covariance matrix in all tracking models is set as In the algorithm, simulation the time interval T is 0.1 s.The model switching probability ρ ij for the multiple model algorithm is a constant matrix with a large value on the diagonal elements.The three-model switching probability can be set the same for three-model IMM, AGIMM and HGMM algorithm as The nine-model switching probability for nine-model IMM algorithm can be set as The initial model probability for the three-model tracking algorithm is taken assuming the second model is near to the real model of HGV; that is, The initial model probability for the nine-model tracking algorithm is taken assuming the first model is near to the real model of HGV; that is,

Simulation Results Comparison and Discussion
We using several different algorithms to track HGV maneuvering trajectory and run 100 Monte-Carlo simulations for every algorithm on each trajectory.We use root mean square error (RMSE) of position and velocity to compare the performances of different tracking algorithms.The RMSE, average RMSE (ARMSE) and peak RMSE (PRMSE) can be defined as where n MC = 100 is the Monte-Carlo simulation number; k/k is the predicted state at discrete time k at i-th Monte-Carlo simulation; X k is the real state at discrete time k.
As defenders, velocity tracking error and position tracking error are the main indicators for HGV tracking performance [2,6,8] and more general target tracking performance [37][38][39].Thus, in following tracking performance comparison, we mainly use velocity RMSE and position RMSE.V RMSE and r RMSE can represent the RMSE of velocity and position, respectively.Moreover, the attack angle error and bank angle error at discrete time k have a certain effect on the velocity and position at discrete time k+1.Thus, V RMSE, r RMSE, α RMSE and ν RMSE using different algorithms will be compared and analyzed as follows.
To compare the tracking performance of MAM algorithm with that of FSMM (three-model IMM and nine-model IMM) algorithm and compare the tracking performance of MAM algorithm with that of VSMM (AGIMM and HGMM) algorithm, we choose a point in every 2 s to plot the former and the later as Figure 2 for tracking Scenario 1 and as Figure 3 for tracking Scenario 2. Furthermore, the average RMSE and peak RMSE using different algorithms are listed as Table 1 for tracking Scenario 1 and as Table 2     V RMSE using AGIMM for tracking scenario 1 V RMSE using HGMM for tracking scenario 1 V RMSE using MAM for tracking scenario 1 r RMSE using AGIMM for tracking scenario 1 r RMSE using HGMM for tracking scenario 1 r RMSE using MAM for tracking scenario 1 αRMSE using AGIMM for tracking scenario 1 αRMSE using HGMM for tracking scenario 1 αRMSE using MAM for tracking scenario 1 νRMSE using AGIMM for tracking scenario 1 νRMSE using HGMM for tracking scenario 1 νRMSE using MAM for tracking scenario 1  From Tables 1 and 2, we can know that the MAM algorithm outperforms multiple model algorithms in almost every ARMSE and PRMSE of velocity and position.Furthermore, VSMM algorithm performs better than FSMM algorithm for tracking Scenario 1 and Scenario 2. In VSMM algorithms, HGMM algorithm performs better than AGIMM algorithm.Although nine-model IMM algorithm uses more models, it does not perform better than three-model IMM algorithm.The tracking results using different algorithms will be analyzed and discussed in detail as follows.
Using MAM algorithm has reduced the average V RMSE by 9.36% than using three-model IMM algorithm for tracking Scenario 1, by 9.52% than using AGIMM algorithm for tracking Scenario 1, by 8.45% than using HGMM algorithm for tracking Scenario 1, by 12.18% than using three-model IMM algorithm for tracking Scenario 2, by 12.20% than using AGIMM algorithm for tracking Scenario 2, by 9.50% than using HGMM algorithm for tracking Scenario 2, respectively.Using MAM algorithm has reduced the average r RMSE by a certain degree than using other algorithms.The results that MAM (single model) algorithm outperforms other algorithms mentioned in this paper for tracking HGV trajectories contradicts the usual opinion that multiple model tracking algorithm outperforms single model (MAM) tracking algorithm for this so-called strong maneuvering target.
The HGV can make maneuver at big degree in long flight distance and maneuver strongly in three dimension from the aspect of kinematic models as Figure 1.However, the maneuver ability is not that strong as we think from the aspect of aerodynamic model because the maneuver is caused by the two control variables and the rates of change of the two control variables are not large.First, by comparing nine-model IMM algorithm with three-model IMM algorithm, we know the modeling V RMSE using AGIMM for tracking scenario 2 V RMSE using HGMM for tracking scenario 2 V RMSE using CKF for tracking scenario 2 r RMSE using AGIMM for tracking scenario 2 r RMSE using HGMM for tracking scenario 2 r RMSE using MAM for tracking scenario 2 αRMSE using AGIMM for tracking scenario 2 αRMSE using HGMM for tracking scenario 2 αRMSE using MAM for tracking scenario 2 νRMSE using AGIMM for tracking scenario 2 νRMSE using HGMM for tracking scenario 2 νRMSE using MAM for tracking scenario 2   From Tables 1 and 2, we can know that the MAM algorithm outperforms multiple model algorithms in almost every ARMSE and PRMSE of velocity and position.Furthermore, VSMM algorithm performs better than FSMM algorithm for tracking Scenario 1 and Scenario 2. In VSMM algorithms, HGMM algorithm performs better than AGIMM algorithm.Although nine-model IMM k−1 as the input of Cubature Kalman filter j and after once recurrence the output of Cubature Kalman filter j can be X j k/k and P j k/k .
using 9-model IMM for tracking scenario 1 V RMSE using 3-model IMM for tracking scenario 1 V RMSE using MAM for tracking scenario 1 r RMSE using 9-model IMM for tracking scenario 1 r RMSE using 3-model IMM for tracking scenario 1 r RMSE using MAM for tracking scenario 1 αRMSE using 9-model IMM for tracking scenario 1 αRMSE using 3-model IMM for tracking scenario 1 αRMSE using MAM for tracking scenario 1 νRMSE using 9-model IMM for tracking scenario 1 νRMSE using 3-model IMM for tracking scenario 1 νRMSE using MAM for tracking scenario
9-model IMM for tracking scenario 2 V RMSE using 3-model IMM for tracking scenario 2 V RMSE using MAM for tracking scenario 2 r RMSE using 9-model IMM for tracking scenario 2 r RMSE using 3-model IMM for tracking scenario 2 r RMSE using MAM for tracking scenario 2 αRMSE using 9-model IMM for tracking scenario 2 αRMSE using 3-model IMM for tracking scenario 2 αRMSE using MAM for tracking scenario 2 νRMSE using 9-model IMM for tracking scenario 2 νRMSE using 3-model IMM for tracking scenario 2 νRMSE using MAM for tracking scenario

Table 1 .
for tracking Scenario 2. The performances for stably tracking Scenario 1.

Table 2 .
The performances for stably tracking Scenario 2.

Table 1 .
The performances for stably tracking Scenario 1.

Table 2 .
The performances for stably tracking Scenario 2.