Maximum Correntropy Unscented Kalman Filter for Ballistic Missile Navigation System based on SINS/CNS Deeply Integrated Mode

Strap-down inertial navigation system/celestial navigation system (SINS/CNS) integrated navigation is a high precision navigation technique for ballistic missiles. The traditional navigation method has a divergence in the position error. A deeply integrated mode for SINS/CNS navigation system is proposed to improve the navigation accuracy of ballistic missile. The deeply integrated navigation principle is described and the observability of the navigation system is analyzed. The nonlinearity, as well as the large outliers and the Gaussian mixture noises, often exists during the actual navigation process, leading to the divergence phenomenon of the navigation filter. The new nonlinear Kalman filter on the basis of the maximum correntropy theory and unscented transformation, named the maximum correntropy unscented Kalman filter, is deduced, and the computational complexity is analyzed. The unscented transformation is used for restricting the nonlinearity of the system equation, and the maximum correntropy theory is used to deal with the non-Gaussian noises. Finally, numerical simulation illustrates the superiority of the proposed filter compared with the traditional unscented Kalman filter. The comparison results show that the large outliers and the influence of non-Gaussian noises for SINS/CNS deeply integrated navigation is significantly reduced through the proposed filter.


Introduction
With the development of missile technology, even higher accuracy of navigation is requested. Single navigation mode cannot satisfy the application requirement. Because of the uncertainty in the actual navigation system, a number of algorithms have been designed to improve the accuracy of state estimation [1]. To sum up, the navigation accuracy of ballistic missile has focused on two parts, navigation system construction and navigation filter algorithm designation.
For navigation system construction, a number of navigation modes have been designed to improve navigation accuracy. Among these modes, the strap-down inertial navigation system (SINS) has been playing an important role among navigation systems for several decades [2][3][4]. It can autonomously achieve a navigation mission completely relying on micro electro mechanical systems (MEMS) and has a strong anti-interference capability. Moreover, the attitude, the velocity, and the position of ballistic missiles can be easily obtained by its output information. However, the accumulated navigation errors caused by the inertial components cannot be avoided, especially in the long-term missile 1. A classical method involves creating filters directly for the systems under the conditions of non-Gaussian noises to deal with noises with different heavy-tailed distributions [23,24]. However, it is difficult to take effect in multidimensional situations, which limits its applicability [25]. 2. There are two efficient methods to handle non-Gaussian noises based on estimating the a posteriori probability density with a series of samples. The UKF approximates the mean and the covariance of the state estimation using unscented transformation (UT) through sigma points [26]. The ensemble Kalman filter (EnKF) is another algorithm to capture state estimation using samples set for the purpose of handling non-Gaussian noises [27]. 3. Another popular method is the Huber-based Kalman filter. It is a robust state estimator to handle the non-Gaussian noises based on the robust estimation theory [28]. Many navigation and target problems were often handled through it under the conditions of non-Gaussian noises in the past few years [29][30][31][32]. 4. Chang [33] proposed a new robust Kalman filter in recent years. He defined a judging index calculating the difference between the observation and the prediction based on Mahalanobis distance to solve the outliers, which can also successfully solve the observation noises that obey the thick-tailed distribution during the actual observation process.
Maximum correntropy criterion (MCC) is an information theoretic method proposed in recent years. It has been applied to the Bayesian estimation and the error analysis [34,35]. The maximum correntropy Kalman filter (MCKF), which is an improved Kalman filter algorithm using the MCC theory, is one of the latest designed estimators for the purpose of solving non-Gaussian noises. In addition, the maximum correntropy extended Kalman filter (MCEKF) and the maximum correntropy unscented Kalman filter (MCUKF) are proposed for the nonlinear system [36], and the MCUKF has been proven to show a better performance than MCEKF [37]. However, its robustness is not analyzed and it has not been applied to SINS/CNS deeply integrated navigation.
The key contributions of this paper are expressed as follows. First of all, the observability of the new deeply integrated navigation system is analyzed. The MCUKF is then deduced on the basis of maximum correntropy theory and unscented transformation to deal with the nonlinearity and non-Gaussian noises in the navigation system. The computational complexity of it is also analyzed. Furthermore, the robustness of the MCUKF is discussed under the condition of the non-Gaussian noises. Finally, the conditions of non-Gaussian noise, mainly including the large outliers and the Gaussian mixture noises, existing in SINS/CNS deeply integrated navigation are effectively handled by the MCUKF.
The organization of this paper proceeds as follows. SINS/CNS deeply integrated navigation is briefly reviewed and the observability is analyzed in Section 1. In Section 2, the algorithm steps are designed for the MCUKF and the robustness performance is analyzed compared with a Huber-based filter. In Section 3, simulations are demonstrated to compare the performances using several filter methods under the conditions of different kinds of noises. Finally, we draw conclusions in Section 4.

Integrated Navigation Principle
The traditional SINS/CNS integrated navigation system outputs the attitude, the velocity, and the position information of a missile through the IMU (including a gyro and an accelerometer). The star sensor of the CNS outputs the attitude information of the missile after star capturing and star identification. The attitude information calculated by the star sensor and the gyro is then integrated, and the integrated information is served as the measurement equation. The optimal estimation methods are then used to estimate the gyro drift to correct the gyro. Figure 1 shows the traditional SINS/CNS integrated navigation principle.
Although SINS/CNS integrated navigation system can correct the gyro drift to some extent, the error of the state variables still has a divergence [4,13,14,16]. SINS/CNS deeply integrated navigation can solve the problem. Based on the traditional integrated navigation system, the height, the pitching angle, and the azimuth angle are added as the measurement information according to the altimeter and the astronomical information. In fact, these three new measurements are the function of the position of the navigation target. Therefore, the estimation of the gyro drift and the accelerometer bias can be gained to correct the acceleration information of SINS. Figure 2 shows the navigation principle of the SINS/CNS deeply integrated mode.  The deeply integrated SINS/CNS integrated navigation system has not yet been realized in real scenarios. When it is applied on the missile, the star sensor measures not only the attitude information but also the partial position information of the missile. The attitude can be calculated through star identification. The position can be calculated by measuring the height angle and the azimuth angle of the star. An altimeter needs to be installed on the system to measure the height information of the missile. The altimeter can be a radar altimeter [38], a pressure altimeter [39], or a laser altimeter [40], decided according to the specific case. Combined with the position information from the star sensor, the position of the missile can be obtained. In this way, the deeply SINS/CNS integrated navigation system can be realized and applied to practical situations. The detailed principle of the system is introduced in the next section.

State Equation
Suppose that the navigation frame is the launch point frame (l-frame) and the SINS is strapdown-installed along the three axles of the missile. Then the state equation is identical to the traditional algorithm: where denotes the platform angles error (the navigation frame misalignment angle) in the l-frame; δv(t) = δv xl (t) δv yl (t) δv zl (t) denotes the velocity error in the l-frame; δr l (t) = δx l (t) δy l (t) δz l (t) denotes the position error in the l-frame; T denotes the process noise vector, including the random noises of the gyro and the accelerometer; G(t) denotes the process noise input matrix; F(t) denotes the process input matrix. The detailed information of F(t) and G(t) can be found in [41]-Equations (6-5) and (6-6), respectively. The solution of the differential Equation (1) for the discrete mode can be expressed as: Equation (2) can be obtained through the Laplace transformation which is introduced in Chapter 2 of Reference [41]. In simulation, the equation is solved through Runge-Kutta algorithm. And the linear version for Equation (2), which is used for MCEKF [36] can be expressed as: where k denotes the epoch, Φ (k − 1) and Γ (k − 1) are the discrete versions of f (X (k − 1)) and G (k − 1), respectively, which can be expressed as:

(1) Measurement of Attitude Error
Suppose that the attitude output of the gyro is θ g (k) = θ g (k) ϕ g (k) ψ g (k) T , and the attitude output of the star sensor is θ The three angles are equivalent to the three Euler angles, denoting roll angle, pitch angle, and yaw angle, respectively. The difference between the two sets of attitude angles can be denoted as δθ (k) = θ g (k) − θ s (k), which can be transformed into the error-angles of the platform and can be expressed as the function of the state Ψ(k). Denote Z 1 (k) = δθ(k) as the attitude error measurement, which can be expressed as: where M(k) = M 0 (k) O 3 × 12 , and M 0 (k) denotes the attitude error transition matrix, which can be expressed as: where ϕ(k) denotes the pitch angle, ψ(k) denotes the yaw angle, and V 1 (k) represents the measurement noises of the star sensor. The detailed derivation process of M 0 (k) can be obtained in [42].

(2) Measurement of Position Error
The measurement information relies on the azimuth angle, the pitch angle of the guide star, and the outputs of the altimeter. The measurement equation of the position error is derived in [4]. The measurement equation is expressed as: where δP(k) is the pitch angle difference between the measurement and the computation value; δA(k) is the azimuth angle difference between the measurement and the computation value; δh(k) is the difference between the height measurement of the altimeter and the computation value; C 1 i stands for the rotation matrix from the l-frame (the launching point coordinate frame which is also the navigation frame) to the i-frame (the geocentric inertial coordinate frame), which is a constant matrix and can be expressed as: where ζ denotes the launching angle of missile, φ 0 and γ 0 denote the longitude and the latitude of the launching point, respectively, C i e (k) denotes the rotation matrix from the i-frame to the e-frame (the earth-based frame), which can be expressed as: where ω e is the earth's rotation angular velocity, t(k) denotes the time in the k-th epoch. δr l (t) = δx l (k) δy l (k) δz l (k) denotes the position error in three directions; v P (k) and v A (k) are the observation noises of the pitch angle and the azimuth angle, respectively; v h (k) is the measurement noise of the altimeter; M 1 (k) and M 2 (k) are as follows: (11) where P(k) is the pitch angle of the guide star, A(k) is the azimuth angle of the guide star, h(k) is the height of the missile, R e is the earth radius, ϕ b (k) is the latitude of the missile, and λ b (k) is the longitude of the missile. Denote Z 2 (k) = δP(k) δA(k) δh(k) T , that is, where (5) and (12), the total measurement equation can be expressed:

Combined with Equations
where the measurement matrix H(k) is: and O is the zero matrix. The measurement noises V (k) are expressed as: and it is independent of the state noises W (k) in Equation (3). Compared with the traditional mode, the deeply integrated mode adds the dimension of the measurements. In this way, the computational complexity will increase, which is analyzed in Section 3.3. In many cases, V (k) denotes the Gaussian white noises [4,13,14,[43][44][45]. However, in practice, due to the sensor faults or the complex environment, the measurement noises may be contaminated noises with different distributions or include large outliers. It can be shown as follows: where V (1) k denotes the Gaussian white noises, V k , and ∆V k is the outlier.
The non-Gaussian noises can influence the performance of the navigation system; therefore, it is necessary to use the robust filter to solve the problem. The maximum correntropy Kalman filter is newly designed to effectively tackle the non-Gaussian noise.

Observability Analysis
The observability analysis is needed to determine the efficiency of the Kalman filter designed for the purpose of estimating the system state variable [46]. If the system observability is poor, the estimation of the state will be inaccurate. Therefore, it is necessary to analyze the observability of the new integrated navigation system.
Considering the state Equation (3) and the measurement Equation (13), according to [47], the observability matrix U 1 of the deeply integrated navigation system is as follows: where H is the measurement matrix of the deeply integrated method, Φ denotes the discrete version of the state transition matrix, n x is the dimension of state variable, which is 15 in the deeply integrated navigation system, as described in Equation (1), so rank (U 1 ) = n x = 15. Hence, the SINS/CNS deeply integrated navigation mode is observable. For comparison, the observability matrix of the traditional navigation mode U 2 is analyzed as follows: where M is the measurement matrix of the traditional method. Thus, rank (U 2 ) = 6 < n x . Therefore, the traditional integrated navigation system is unobservable. In other words, the new navigation system, compared with the traditional one, will have a more precise estimation.

Maximum Correntropy Criterion
Correntropy [48] denotes the measure of an information entropy field, which defines the novel metric in the sample space, using the information available in the higher-order statistics of the signals. Given two random variables with the joint probability density function (PDF) p (a, b), the definition of correntropy is expressed as follows. A and B, the correntropy is defined as:

Definition 1. Given two random variables
where E denotes the expectation, κ λ ( , ) denotes a real-valued continuous, symmetric, and nonnegative definite kernel function λ > 0 denotes the kernel bandwidth.
The Gaussian kernel, which is expressed as follows, is used as the kernel function in Definition 1: where e = a − b. In many practical cases, there are a finite amount of samples to estimate the unknown joint PDF p (a, b). Hence, the sample mean estimator can calculate the correntropy as follows: where 5], and p (a, b) = 1 2 e −(a 2 + b 2 ) . The current correntropy value is depicted in Figure 3.
According to Figure 3, correntropy is considered a correlation scale in the joint space, which reaches maximum under the conditions that a (i) = b (i). Based on this, the cost function can be defined to represent the maximum correntropy criterion (MCC) as follows: From Figure 3, it is manifest that MCC represents a local similarity criterion that reflects the maximum error probability density from the view of probabilistic meaning [48]. Non-Gaussian noises, such as large outliers and Gaussian mixture noises, can be handled effectively on the basis of the property [49]. That is to say, MCC reaches the maximum value on the joint space's bisector because the value of Gaussian kernel maximizes on the line A = B, which is different from the mean square error (MSE) estimation [50]. Based on this, MCC has been used for parameter estimation in recent years [34,35,51].

The Maximum Correntropy Unscented Kalman Filter
The MCUKF is a novel algorithm combined with a UKF framework based on the MCC, and can perform better for nonlinear systems in non-Gaussian noise environments [37].
Consider the state equation and the measurement equation described in Equations (2) and (13) in Section 2: W (k − 1) and V (k) are mutually uncorrelated process noises and measurement noises, respectively. The mean of both noises are zero and the covariances of them are expressed as: Similar to the Kalman filter, the MCUKF includes the time update and the measurement update.

Time Update
The sigma points set are generated from the estimated stateX (k − 1|k − 1) and the covariance matrix P (k − 1|k − 1) at the previous step (k − 1).
where n x denotes the dimension number of the state, (n x + λ) P (k − 1|k − 1) i denotes the i-th column of (n x + λ) P (k − 1|k − 1) , and the composite scaling factor λ is expressed as: where κ is a parameter that is set to (3 − n x ) [37], and α controls the distribution conditions of the sigma points.
The unscented transformed points are given as: The predicted state and the covariance matrix are estimated as: where the weights of the sigma points are as follows: where α is expressed in Equation (28) and β is a weighted parameter that is non-negative. When the measurement data are not avaliable,X (k|k − 1) can replace the current state estimation. The sigma points of the prior state mean and the predicted covariance are set as: The sigma points are then transferred using the measurement equation as: The predicted measurement mean can be estimated as:

Measurement Update
Next, the measurement update is calculated based on the MCC. For the model described in the previous section, we have where I denotes the unit matrix, and u (k) = where B (k), B p (k|k − 1), and B r (k) is calculated on the basis of Cholesky decomposition of E u (k) u T (k) , the predicted covariance P (k|k − 1) , and the covariance of the measurement noises R (k), respectively. In fact, u denotes a random variable matrix, which consists of the prior estimation error of the state and the measurements prediction error. Equation (36) is multiplied by B −1 (k) on the left of both sides, and we get where , and e (k) = B −1 (k) u (k). In this way, e (k) denotes a new random variable matrix, which is just multiplied by a matrix B −1 . Now the MCC for Equation (24) based on the cost function can be rewritten as follows: where e i (k) denotes the i-th unit of e (k), d i (k) denotes the i-th unit of D (k), w i (k) denotes the i-th row of w (k), and N = n x + n z denotes the entire dimension of D (k), which represents the sum of the dimension of state variable X (k) and the measurement Z (k). In fact, d i (k) is equivalent to a i in Equation (23), and w i (k) X (k) is equivalent to b i in Equation (23). That is, e i (k) is equivalent to e i in Equations (22) and (23), and J MCC (X (k)) is equivalent to C λ (A, B) in Equation (22). In particular, according to Figure 3, when the error e i (k) reaches the minimum, the correntropy reaches its maximum value. In other words, the optimal estimation of the state can minimize the error and maximize the correntropy as well. In this way, the optimal estimate of X (k) is considered with respect to MCC.
Then, the optimal estimate under the MCC criterion of X (k) iŝ The optimal solution can thus be obtained by solving It follows that The optimal solution expressed by Equation (42) is actually a fixed-point equation of X (k) and can be rewritten as: with A fixed-point iterative algorithm can be readily obtained as [51]: whereX(k) t * denotes the solution at the fixed-point iteration t * .
Equation (45) can be further expressed as follows (the detailed derivation is in [52]): With the above derivations, we summarize the proposed MCUKF algorithm as follows: Step 1: Choose a proper kernel bandwidth λ and a small positive number ε ; the estimateX (0|0) and the covariance matrix P (0|0) are initialized; let k = 1.
Step 5: Use Equations (48)-(54) to computeX(k|k) t * : Step 6: Compare the estimation of the current step and the estimation of the last step. If Equation (55) holds, setX (k|k) =X(k|k) t * and continue to Step 6. Otherwise, go back to Step 4.
Step 7: Update the estimation covariance matrix, make k + 1 → k, and go back to Step 2.
For obtaining the optimal solution in Equation (42), the MCUKF updates the state estimation through a fixed-point algorithm, which is shown as Equation (55). Therefore, the small positive number ε provides a threshold for the fixed-point iteration. When the result of the iteration is in the threshold range, we call it convergence and the solution is optimal. The convergence of the solution will be fast because of the predicted stateX (k|k − 1) [52].

Computational Complexity
For the future implementation on real hardware, it is necessary to analyze the computational complexity of the proposed approach in the article. In terms of the floating points operations, the computation complexity of the basic equations of the MCUKF is analyzed in Table 1. Table 1. Computational complexities of some equations, where n denotes the dimension of the state variable, m denotes the dimension of the measurements. The traditional unscented Kalman filter is analyzed in [53]. The computation complexity of it is S UKF = 26 3 n 3 + 15n 2 + 10n 2 m + 5n + 8nm 2 + 6mn + m 3 + 3m 2 + 3m. (57) In fact, Equations (27)- (35) in the MCUKF are the same as that in the UKF. Equations (48)-(56) in the MCUKF mainly involves the different equations expressed by Equations (48)- (56). Equation (55) gives a fixed-point algorithm to update the posterior estimate of the state that can provide a stop condition in several steps. Assuming that the fixed-point iteration number is T, then the computation complexity of MCUKF is

Equation Addition/Subtraction and Multiplication Equation Addition/Subtraction and Multiplication
According to Equations (57) and (58), the computational complexity of the MCUKF is much higher than that of the UKF because of the fixed-point algorithm and more matrix inverse problems.
According to the complexity analysis above, the increasing dimension of the measurements of the deeply integrated mode has a higher computational complexity than the traditional mode.

Robustness Analysis
The robustness is a very important feature for the navigation filter algorithm, especially in the complex space environment. One of the most significant advantages of MCC is robustness that can tackle non-Gaussian noises. In fact, it uses a method similar to the Huber-based filter (HF) to improve the robustness [54]. Both of them use the weighted function to tackle the non-Gaussian noises. HF is a classical method to tackle the non-Gaussian noises, and it has been proved in detail [55]. Additionally, some articles have demonstrated the ability of MCC in accommodating non-Gaussian noises, especially outliers and Gaussian mixture noises [56,57]. Next, the algorithm robustness is compared with HF through the weighting function.
HF is created based on the concept of the robustness through the maximum likelihood method. Assume that the maximum likelihood estimator (MLE) stems from the data observation Z = {z 1 , ..., z N } and the state variable X = {x 1 , ..., x N }. Then e = Z − HX is drawn from the distribution p (e; X), which is known apart from the parameters X. In this way, the MLE of X can be defined as: Assuming that e 1 , ..., e N are i.i.d. observations, then it can be rewritten as: Let ρ (e i ; X) = a (− ln p (e i ; X)) + b, where a > 0 and b is uncorrelated with X. Thus, X MLE obeys In fact, the errors in the target track model are always considered as pure Gaussian white noises with zero mean and constant variance σ 2 in some parametric models. In these cases, the maximum likelihood estimation can be solved. However, in practice, the external environment influences the measurement equipment and causes the filter algorithm to break down. Hence, the non-Gaussian noises in the measurements cannot be avoided. A typical non-Gaussian noises are Gaussian mixture noises, which obey the Gaussian distributions with different mean values and variances with most errors. For instance, assume that the distribution is expressed as: where F ( ), F 0 ( ), and F c ( ) denote the Gaussian mixture distribution, the known normal Gaussian distribution, and an unknown contaminated distribution, respectively, and υ is the contamination ratio. It is much smaller than 1. For Equation (62), it is of necessity to solve the maximum likelihood estimation through Equation (61).
For the purpose of solving the problem, Huber proposed that ρ can change freely within limits, leading to the birth of M-estimation, which is an effective way of solving Equation (61).
M-estimator is defined as a generalization of MLE as: where ρ (e) denotes the cost or penalty function. Let where λ denotes the kernel bandwidth. Then Therefore, MCC is equivalent to Equation (64): According to the M-estimation problem, ψ MCC (e i ), as well as the weighting function of MCC, can be rewritten as: where ρ MCC (e i ) = dρ MCC (e i ) de i . Therefore, On the basis of analysis of the weighting function above, we will compare it with the Huber filter to prove the superiority of MCC.
The weighted functions of them are shown in Table 2. Table 2. The weighted functions of two methods.

Item Penalty Function Weighted Function
Huber Some researchers have performed robustness analysis of MCC in theory. The robustness of two methods are shown through Figure 4 to intuitively analyze the robustness of MCC.
In Figure 4, the weight of HF decreases rapidly as the estimate residual grows. However, the weight of MCC approaches zero faster compared with that of the Huber method. In addition, with the growth of the error, the weight of the Huber method is not zero. In other words, Huber-based methods can select more measurements containing large errors than MCC-based methods, which may cause large errors to the filter algorithm, although with small weights. Therefore, the filter based on the MCC is more robust than the Huber-based filter.

Simulation Results
In this section, the superiority of SINS/CNS deeply integrated navigation is proved firstly. The performances of the MCUKF are then illustrated under the conditions of the Gaussian noises, the large outliers, and the Gaussian mixture noises, respectively. The estimation performance are compared on the basis of the following indices: where m is the simulation time steps and L represents the total simulation numbers. In addition, the time cost of every method is given to compare the computational complexity. All methods are implemented in the same precision (64-point floating point) in MATLAB running on a PC with processor Intel(R) Core(TM) i7-7500U CPU 2.70GHz and with 8 GB of installed memory (RAM). The simulation parameters are as Table 3. The missile trajectory is shown in Figure 5. To prove the high performance of the new integrated navigation system, a simulation is done only under the condition of the Gaussian white noises. Compared with the traditional integrated mode, the position and the attitude angle residual are shown in Figures 6 and 7.  Table 4. From the Figures and the table above, the navigation position error has evidently improved when constructing the measurement Z 2 in Equation (12). In addition, the SINS/CNS deeply integrated navigation mode for the ballistic missile can tackle the divergence of the position error effectively, which obeys the observability analysis results. However, the newly proposed navigation system needs more time to gain the results.
Next, we will compare the MCUKF performance with other algorithms under the conditions of the Gaussian noises, the large outliers, and the Gaussian mixture noises for SINS/CNS deeply integrated navigation.

Case 1: Gaussian Noises
In this section, we will apply the MCUKF on SINS/CNS deeply integrated navigation under the conditions of the Gaussian noises to compare performances with other filter algorithms including the HKF [55], the IEKF, the NIF, the PLF, the UKF, and the MCEKF [36]. The simulation parameters are set as Table 1. The simulation results are shown in Figures 8 and 9, and Table 5.  According to Table 5, PLF demonstrates the highest accuracy among the seven filters in this condition because of the a posteriori estimation and the iteration. It can be observed that robust filters cannot always perform as well as the PLF, the NIF, or the UKF under the condition of the Gaussian noises.

Case 2: Large Outliers
In this section, the case in which the measurement value including large outliers is considered. The attitude outliers are added artificially to the attitude angle data at the 5000 th and 9000 th epochs, respectively. The position outliers are also added artificially to the position data at the 5000 th epochs. The constant outliers are added to all epochs between the 4001 th to the 4005 th to the attitude angle data. The constant outliers are added to all epochs between 4001 th and the 4005 th to the position data. The obtained position errors and the attitude angle errors are displayed in Figures 10 and 11, and Table 6.
According to the presented results, it can be concluded that, in the presence of outliers, the filter results mostly depend on the outliers. In addition, Figures 10 and 11 show that the HKF, the MCEKF, and the MCUKF can resist the influence of the outliers to some extent. Furthermore, the two figures show that the MCUKF has better performance than the other algorithms, so the influence of the single and the constant outliers is reduced using the MCUKF algorithm. The RMSE values and time cost of all used filters are presented in Table 6.
Compared with other filters, the filtering algorithm of the MCEKF and the MCUKF is improved to a certain degree, especially for the outliers. The MCUKF has better performance than other filters and lower time cost than the iterated filters.

Case 3: Gaussian Mixture Noises
In this section, the case in which the measurement noises are the Gaussian mixture noises, whose distributions are as follows: Star Sensor Noise: v s ∼ 0.9N 0, (3 ) 2 + 0.1N 0, (30 ) 2 .
Altimeter Noise: v h ∼ 0.9N 0, (50m) 2 + 0.1N 0, (500m) 2 . The obtained position errors and the attitude angle errors are displayed in Figures 12 and 13, and Table 7.   Table 7 illustrates the RMSE of the position, the attitude angle, and the time cost of every method. It can be observed that the MCUKF achieves better performance than other filters.

Conclusions
SINS/CNS deeply integrated navigation, which consists in modified SINS/CNS navigation, is introduced, and the observability of the navigation system was analyzed. However, the simulation results demonstrate that the new system needs more time. Measurement noises are often the research focus with regard to navigation systems. The non-Gaussian noises, including the large outliers and the Gaussian mixture noises, always exist during the measurement process. On the basis of the maximum correntropy theory, the maximum correntropy unscented Kalman filter (MCUKF) is a newly designed robust unscented Kalman filter (UKF). The robustness for state estimation, which can resist the effects of non-Gaussian noises effectively, is also analyzed and compared with a Huber-based filter. Simulated experimental results demonstrate that the MCC-based filter provides much better performance when tackling the non-Gaussian noises of SINS/CNS deeply integrated navigation, but with a higher computational complexity than the traditional UKF.
In practical applications, the added star sensors and the altimeter of the deeply integrated system will increase the load of the missile. Future work will focus on the installation method of the measurement facilities and applications in real scenarios. The coupled problem of multiple sensors may also need to be solved.