Next Article in Journal
Toward the Ultimate-High-Speed Image Sensor: From 10 ns to 50 ps
Previous Article in Journal
Fast Visual Tracking Based on Convolutional Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Robust Unscented Kalman Filter via Fading Factor and Maximum Correntropy Criterion

1
The School of Automation, Beijing Institute of Technology, Beijing 100081, China
2
The Second Construction Limited Company of China Construction Eighth Engineering Division, Jinan 250014, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2018, 18(8), 2406; https://doi.org/10.3390/s18082406
Submission received: 1 June 2018 / Revised: 19 July 2018 / Accepted: 21 July 2018 / Published: 24 July 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
In most practical applications, the tracking process needs to update the data constantly. However, outliers may occur frequently in the process of sensors’ data collection and sending, which affects the performance of the system state estimate. In order to suppress the impact of observation outliers in the process of target tracking, a novel filtering algorithm, namely a robust adaptive unscented Kalman filter, is proposed. The cost function of the proposed filtering algorithm is derived based on fading factor and maximum correntropy criterion. In this paper, the derivations of cost function and fading factor are given in detail, which enables the proposed algorithm to be robust. Finally, the simulation results show that the presented algorithm has good performance, and it improves the robustness of a general unscented Kalman filter and solves the problem of outliers in system.

1. Introduction

In real-world applications, target tracking problems have attracted much attention such as maneuvering target tracking [1], ballistic target tracking [2] and multiple target tracking [3], etc. For getting better accuracy, efficiency and performance of tracking problems, the effect of noise needs to be reduced, especially of measurement noise. Therefore, the measurement noise needs to be restrained timely in the process of a tracking target. As is well known, filtering algorithms are powerful tools for suppressing the effect of noise. Kalman [4] first proposed the Kalman filter (KF) algorithm in 1960, and the filtering technique has developed quickly ever since. For the linear Gaussian process [4,5,6,7], the KF could get optimal recursive results based on the minimum mean square error (MMSE) estimation. It is well known that KF can employ the first two moments (mean and covariance) of state and measurement to obtain optimal estimates. It can be seen that the KF has been applied in numerous areas such as navigation, target tracking, communications [7], attitude determination [8], multiple sensors data [9] and many more. However, for amounts of data collected or sent by sensors, the KF experienced heavy computational load to get an optimal solution. In most practical applications, the systems are nonlinear. In this case, the performance of KF is unsatisfactory. Furthermore, when a dynamic model was contaminated by outliers [10], the KF degraded severely. Therefore, it is necessary to design an efficient algorithm that can suppress the outliers for the nonlinear problem.
In fact, for a nonlinear tracking problem, many nonlinear filtering algorithms were proposed to deal with the nonlinear systems—for example, extended Kalman filter (EKF), unscented Kalman filter (UKF), cubature Kalman filter (CKF), particle filter (PF), Gaussian sum filter (GSF), and so on. Particularly, it is well known that EKF deals with the nonlinear system by its first-order Taylor Series expansions and leads to a suboptimal solution. UKF, based on unscented transform (UT), is an improved filter compared to KF and EKF. It selects deliberately as sigma points to generalize the KF for both linear and nonlinear systems, and it can be used to propagate three order information for state estimation. Therefore, it can achieve better estimation accuracy and computational feasibility [11]. Zhan and Wan [12] developed the iterated unscented Kalman filter for passive target tracking. Based on the spherical-radial cubature rule, CKF is presented to approximate the Gaussian filter [13,14]. GSFs are based on the idea that a non-Gaussian probability density function (PDF) can be approximated by a sum of Gaussian PDFs [15]. Based on Bayesian estimation and a sequential Monte Carlo approach, Du et al. utilized PF to handle nonlinear and non-Gaussian problems, and the PF is applied in small target tracking in an optimal image sequence [16]. However, these nonlinear filters were susceptible to outliers and did not have robust property.
For measurement model contaminated by outliers, several filtering algorithms have been proposed [17,18,19,20,21,22,23]. Based on the Huber function, Wang et al. [17] presented the derivative-free filter to manage the measurement outliers, but it did not suppress the state outliers well. For measurement outliers, Durovic and Kovacevic [18] utilized the M-estimation to deal with measurement outliers in the presence of disturbance uncertainty. It can not deal with both state and measurement outliers. Their performance will diverge if the state and measurement models were contaminated by outliers simultaneously. To solve a nonlinear system with heavy-tailed noise, Wang et al. [19] studied a robust information filter to solve the measurements with a large error from the estimation process. However, it did not embed the fading factor into the framework of information filter. Under this case, Karlgaard [20] proposed an adaptive robust nonlinear filtering algorithm to resist the effects of outliers. For both the state and measurement outliers, Gandhi and Mili [21] introduced a generalized maximum likelihood type KF. However, the KF was limited to a linear system, and the evaluation of Jacobian matrices in EKF could be nontrivial and this leads to degraded performance. Chang et al. [22] investigated a robust filter to suppress the state and measurement outliers, and it utilized the robust property of the Huber function. The H filter has robustness to minimize the estimation error [23], and it can not accommodate the outliers well for outliers occurring randomly. However, these filters have robustness to some extent.
In order to enhance the robustness of the aforementioned filters, some fading factors were proposed to embed into the above filters to keep them stable and effective [24,25,26]. Wang et al. [24] introduced a modifiable fading factor to tackle the nonlinear estimation problem. After that, several researchers also investigated this problem. Yang et al. [25] investigated the adaptive robust filtering via the robust maximum likelihood estimation, but it can not control the dynamics model biases. Safarinejadian and Yousefi [26] proposed an adaptive fading memory KF to deal with static alignment of inertial navigation systems. Furthermore, Geng and Wang [27] utilized multiple fading factors in KF to handle the filtering divergence with inaccurate system noise. In [25,26,27], the fading factors were embedded into KF, but not nonlinear filters. Therefore, in order to overcome their shortcomings, Karlgaard [20] and Wang et al. [24] embedded the fading factors into nonlinear filters to handle nonlinear problems.
Motivated by the above discussion, this work proposes a new adaptive robust UKF scheme based on both fading factor and maximum correntropy criterion (MCC) to focus on the state estimation problems with measurement outliers. To the best knowledge of the authors, for the nonlinear tracking problem, the filtering algorithm based on MCC and fading factor had not been studied before. In our research, tracking a moving target by a sensor has been carried out to compare with UKF and an adaptive Huber unscented Kalman filter [24]. Furthermore, the prime dedications of this paper are as follows: (1) the cost function of the proposed filter is proposed, the fading factor is applied in the state model and the maximum correntropy criterion is used in the measurement model; (2) the proposed filter can suppress the effect of outliers effectively in the dynamics system (state or/and observation equation); and (3) the proposed filter is easy to perform and has better performance in the presence of outliers.
The structure of this paper is organized as follows. The Section 2 briefly introduces the fundamentals of the proposed filter. Section 3 proposes the adaptive robust unscented Kalman filter based on fading factor and MCC for a nonlinear system. Simulation results and comparisons are provided to confirm the feasibility and superiority of the proposed filter in Section 4. Finally, conclusions are shown in Section 5.

2. Fundamentals of the Proposed Filter

2.1. Maximum Correntropy Criterion

For any random variables X and Y, the correntropy is defined as
V ( X , Y ) = E X , Y [ ρ ( X , Y ) ] = ρ ( x , y ) d F X , Y ( x , y ) ,
where E [ · ] represents the mean value, and ρ ( · ) is a real-valued continuous, symmetric and nonnegative definite kernel function, respectively. F X , Y ( x , y ) is the joint probability density function of X and Y. However, F X , Y ( x , y ) is usually unknown, and numbers of data samples, ( x i , y i ) , ( i = 1 , 2 , , N ) can be obtained. Therefore, the correntropy can be computed as follows:
V ˜ ( X , Y ) = 1 N i = 1 N ρ ( x i , y i ) .
In theory, various kernel functions can be used. In general, the Gaussian kernel function is selected as follows:
κ G ( e i ) = κ G ( x i y i ) = 1 2 π σ exp { e i 2 2 σ 2 } ,
where e i = x i y i , and σ is the kernel width of correntropy, respectively. If e i = 0 , κ G ( e i ) can reach its maximum value. Therefore, the cost function of MCC is expressed by
J M C C = min ( i = 1 N ( κ G ( 0 ) κ G ( e i ) ) ) .

2.2. Cost Function of Adaptive Robust Kalman Filter

In this section, the cost function of the adaptive robust Kalman filter is derived by analyzing the following linear state model and measurement model:
x k = Φ k | k 1 x k 1 + w k 1 y k = H k x k + v k ,
where k denotes discrete time, x k R n represents n × 1 system state vector at time k, Φ k | k 1 is the state transition matrix, y k R m is m × 1 system measurement vector at time k, H k is observation matrix, w k 1 is the process noise with the covariance Q k 1 , v k is the measurement noise with the covariance R k , and they are uncorrelated zero-mean Gaussian white noises.
For simplicity, the problem of state and measurement outliers is focused on in this work. For the linear state space model (5), a combined cost function is used to perform two different criterions to the state model and measurement model. Utilizing the Bayesian maximum likelihood, the posterior mean estimate is derived by minimizing the following function
x ^ k | k = arg min ( x k x ^ k | k 1 P k | k 1 1 2 + H k x k y k R k 1 2 ) ,
where x A 2 = x T A x is the quadratic form with respect to A (A is nonnegative definite matrix ); x ^ k | k is the posteriori estimate, x ^ k | k 1 is the priori estimate, P k | k 1 is the covariance matrix of x ^ k | k 1 .
Denoting ξ k = R k 1 / 2 ( H k x k y k ) and utilizing MCC in Equation (4), the cost function of the filtering algorithm in Equation (6) is formulated as follows:
x ^ k | k = arg min ( x k x ^ k | k 1 P k | k 1 1 2 + j = 1 m ( k G ( 0 ) k G ( ξ k , j ) ) ) ,
where the term m is the dimension of the measurement model. Furthermore, the fading factor is embedded to strengthen the robustness property of the KF with model error. Then, the cost function in Equation (7) can be rearranged as
x ^ k | k = arg min ( x k x ^ k | k 1 ( α k P k | k 1 ) 1 2 + j = 1 m ( k G ( 0 ) k G ( ξ k , j ) ) ) .
Differencing Equation (8) with respect to x k , one has that
( α k P k | k 1 ) 1 ( x k x ^ k | k 1 ) i = 1 n k G ( ξ k , i ) ξ k , i ξ k , i x k = 0 .
Substituting Equation (3) into Equation (9), it can be obtained that
( α k P k | k 1 ) 1 ( x k x ^ k | k 1 ) + 1 2 π σ 3 i = 1 m exp ( ξ k , i 2 2 σ 2 ) ξ k , i ξ k , i x k = 0 .
Define the function
Ψ j = 1 2 π σ 3 i = 1 m exp ( ξ k , i 2 2 σ 2 )
and the matrix
Ψ = diag [ Ψ j ] j = 1 , , m .
Thus, we find out that Equation (10) can be redescribed as follows:
( α k P k | k 1 ) 1 ( x k x ^ k | k 1 ) + H k T R T / 2 Ψ ξ k = 0 .
Substituting ξ k into Equation (13), we arrive at
( α k P k | k 1 ) 1 ( x k x ^ k | k 1 ) + H k T R T / 2 Ψ R k 1 / 2 ( H k x k y k ) = 0 .
Equation (14) satisfies the minimization solution of the cost function as follows:
x ^ k | k = arg min ( x k x ^ k | k 1 P ˜ k | k 1 1 2 + H k x k y k R ˜ k 1 2 ) ,
where P ˜ k | k 1 = α k P k | k 1 and R ˜ k = R k T / 2 Ψ k 1 R k 1 / 2 . Comparing Equation (15) with Equation (6), it can be seen that covariance matrixes of them are different, and others are identical. According to the iterative equations of KF, the explicit solution of Equation (15) can be performed as
K k = P ˜ k | k 1 H k T ( H k P ˜ k | k 1 H k T + R ˜ k ) 1 ,
x ^ k | k = x ^ k | k 1 + K k ( y k H k x ^ k | k 1 ) ,
P ˜ k | k = ( I n K k H k ) P ˜ k | k 1 .
Remark 1.
From an engineering viewpoint, it can be seen that the real state vector x k can not be obtained. However, some methods are proposed to change the proportion of state vector, and they can not utilize the state vector—for example, maximum correntropy criterion, Huber function, and so on. Although the state vector was coped with in [22], the real state x k was replaced by estimation vector x ^ k | k 1 in the process.

2.3. Formation of the Fading Factor

The state estimate of Kalman filter depends on the ratio of new measurements and the ones which are based on predicted state vector, dynamics model, and all previous measurements. If the state model error is much larger than that of the measurement model, it is obtained that the information from new measurements will be ignored. Thus, the result of the filter will become poor. The fading factor in Equation (9) is adopted to ensure the performance of the filter in the presence of dynamics model error. For improving the utilization rate of the new measurements, the fading factor α k in Equation (8) becomes greater than 1. Subsequently, the variance matrix P k | k 1 is inflated. Therefore, it is obtained that the contribution of x ^ k | k 1 to x ^ k | k is reduced, and the impact of dynamics model error would be small.
Next, the fading factor α k in Equation (5) is analyzed and derived via the innovation sequence orthogonal principle [28]
E { v k + j v k T } = 0 j = 0 , 1 , 2 , k = 1 , 2 , ,
where E { . } and v k are expected value and innovation sequence, respectively.
For the state model and measurement model stated by Equation (5), the lemma in [28] is given as follows:
Lemma 1.
If x k x ^ k | k is significantly smaller than x k , then, for any j,
C j , k = E { v k + j v k } = Θ ( k + j , , k , x ^ k + j | k + j 1 , , x ^ k | k 1 ) ( P x k y k K k C 0 , k ) = 0 ,
where P x k y k is the cross covariance between state and measurement, K k represents the Kalman gain, C 0 , k equals to E { v k v k } , and Θ ( k + j , , k , x ^ k + j | k + j 1 , , x ^ k | k 1 ) can be written as follows:
Θ ( k + j , , k , x ^ k + j | k + j 1 , , x ^ k | k 1 ) = H k + j l = k + 1 k + j 1 Φ l ( I K l H l ) Φ k .
The proof of Lemma 1 is omitted here. Please refer to Zhou et al. [28] if needed.
It can be shown that Lemma 1 holds strictly for a linear system, and it is approximately true for a nonlinear system. A sufficient condition of Equation (20) is given as follows:
P x k y k K k C 0 , k = 0 .
Using the common criterion of UKF leads to K k = P x k y k P y k y k 1 , and
P x k y k ( I P y k y k 1 C 0 , k ) = 0 ,
where P y k y k is the error covariance matrix of measurement, and a sufficient condition of Equation (23) is given as follows:
P y k y k = C 0 , k .
According to equations of KF that P y k y k = H k P k | k 1 H k T + R k , Equation (24) can be rewritten as
H k P k | k 1 H k T + R k C 0 , k = 0 .
If there exists a fading factor such that Equation (25) holds, the innovations v k + j and v k become approximately orthogonal. Thus, Equation (25) can be redescribed as
H k α 0 P k | k 1 H k T = C 0 , k R k .
Noting that the terms H k , P k , C 0 , k and R k are full rank, by using the property of matrix trace, we have
t r [ H k α 0 P k | k 1 H k T ] = t r [ C 0 , k R k ] .
Because α 0 is a scalar, then it can be expressed as follows:
α 0 = t r [ C 0 , k R k ] t r [ H k P k | k 1 H k T ] .
It seems that Equation (28) can only be calculated if the state model (5) is a linear or linearized system. However, the linearization term H k P k | k 1 H k T can be readily approximated by the other nonlinear methods—for example, EKF and UKF, and so on. The focused UKF is utilized in this work. Therefore, the fading factor α 0 can be applied to the UKF framework, and it is redescribed as follows:
α 0 = t r [ C 0 , k R k ] t r [ P y k y k R k ] ,
where C 0 , k can be expressed as follows:
C 0 , k = v k v k T , k = 1 , λ C 0 , k 1 + v k v k T 1 + λ , k > 1 ,
where the forgetting factor λ is commonly set as 0.95 [19].
Remark 2.
In the framework of the proposed filter, the fading factor α 0 is a single fading factor. For the complicated multivariable systems in actual applications, it is not enough to use a fading factor. To overcome the shortcomings, Yang et al. [25] presented a multiple fading factor Kalman filter to deal with multivariable systems, and α 0 is replaced by a matrix of fading factors Λ = d i a g ( α 0 , α 1 , , α n ) .

3. The New Adaptive Robust Unscented Kalman Filter

In this section, we will focus on the proposed unscented Kalman filter based on fading factor and maximum correntropy criterion. The cost function of the proposed filter in Equation (15) has been derived, and it is embedded into the structure of UKF. In addition, it is well known that UKF has some good properties: easy implementation, appropriate performance and computational feasibility. Therefore, it is very popular in the nonlinear system.
Suppose that the nonlinear discrete-time system can be modelled as follows:
x k = f ( x k 1 ) + ω k 1 y k = h ( x k ) + r k ,
where x k is the state vector with the covariance matrix P k , and the other terms have the same meaning as the above terms in Equation (5), respectively. The initial estimate state and its covariance matrix are given as x ^ 0 = E ( x 0 ) and P 0 = E [ ( x 0 x ^ 0 ) ( x 0 x ^ 0 ) T ] . Similarly, the estimate state and its covariance matrix at time k 1 are given as
x ^ k 1 = E ( x k 1 ) ,
P k 1 = E ( x k 1 x ^ k 1 ) ( x k 1 x ^ k 1 ) T .
The procedure of the proposed adaptive robust unscented Kalman filter depicted in Table 1 is described as follows:
Step 1: The Unscented Transform with 2n + 1 symmetric sigma points χ i , k 1 and weights update:
χ 0 , k 1 = x ^ k 1 , χ i , k 1 = x ^ k 1 + n + δ ( P k 1 ) i , i = 1 , 2 , , n , χ i , k 1 = x ^ k 1 n + δ ( P k 1 ) i , i = n + 1 , n + 2 , , 2 n ,
where δ = φ 2 ( n + κ ) n , φ , which ranges from 0 to 1, controls the distribution of sigma points. κ equals 3 n , and P k 1 is the Cholesky factor of P k 1 , respectively.
The corresponding weights ω 0 m and covariance ω 0 c are as follows:
ω 0 m = δ n + δ , ω 0 c = δ n + δ + ( 1 φ 2 + η ) , ω i m = ω i c = 1 2 ( n + δ ) i = 1 . , 2 n ,
where η is set as 2 generally for Gaussian distribution and relates with the prior distribution of state.
Step 2: Time update
χ i , k | k 1 = f ( χ i , k 1 ) x ^ k | k 1 = i = 0 2 n ω i m χ i , k | k 1 .
Step 3: Modified measurement covariance update
y i , k | k 1 = h ( χ i , k 1 ) y ^ k | k 1 = i = 0 2 n ω i m y i , k | k 1 ,
ξ k = R k 1 / 2 ( y ^ k | k 1 y k ) ,
Ψ = diag [ ψ k , j ( ξ k , j ) ] j = 1 , 2 , , m ,
R ˜ k = R k T / 2 Ψ 1 R k 1 / 2 .
Step 4: Fading factor update
v k = y k y ^ k | k 1 ,
C 0 , k = v k v k T , k = 1 , λ C 0 , k 1 + v k v k T 1 + λ , k > 1 ,
α 0 = t r [ C 0 , k R k ] t r i = 0 2 n ω i c ( y k , i y ^ k | k 1 ) ( y k y ^ k | k 1 ) T ,
α k = α 0 , α 0 > 1 , 1 , α 0 1 .
Step 5: Measurement update
P k | k 1 = α k i = 0 2 n ω i c ( χ i , k | k 1 x ^ k | k 1 ) ( χ i , k | k 1 x ^ k | k 1 ) T + Q k ,
P y k y k = α k i = 0 2 n ω i c ( y k , i y ^ k | k 1 ) ( y k , i y ^ k | k 1 ) T + R ˜ k ,
P x k y k = α k i = 0 2 n ω i c ( χ i , k | k 1 x ^ k | k 1 ) ( y i , k | k 1 y ^ k | k 1 ) T ,
P k | k = P k | k 1 K k P y k y k K k T ,
K k = P x k y k P y k y k 1 ,
x ^ k | k = x ^ k | k 1 + K k ( y k y ^ k | k 1 ) .
Remark 3.
In the process of nonlinear filter, it can be seen that both fading factor and maximum correntropy criterion are applied in the cost function of the new filter. We can obtain that the estimate state x ^ k is more accurate than estimate state x ^ k using other nonlinear filters such as EKF, UKF and other common nonlinear filters, etc. First, the proposed filter has better adaptive ability to balance the contribution between the process model information and measurement on the state vector. Second, the proposed filter can retain the effect of outliers. Particularly, when the state or/and measurement model are contaminated by outliers, the effectiveness of the proposed filter is much better than that of UKF.

4. Simulation and Comparison

To illustrate the practical applicability of the proposed nonlinear filtering algorithm, two classical filtering applications are employed in this section.

4.1. Radar Tracking System

We track a moving object by a radar which utilizes measurements of distances information. The two-dimensional system uses a single station for tracking targets, and the state vector includes position and velocity. The dynamics model moves within two dimensional plane according to the standard dynamics model [29]
X k = Φ X k 1 + Γ u k ,
where X k = [ x , x ˙ , y , y ˙ ] , x and y are Cartesian coordinates of the moving target, and x ˙ and y ˙ are correlative velocities of the moving target, respectively. u k is zero mean Gaussian white noise with the covariance Q = d i a g ( [ 10 4 , 10 4 ] ) . The transition matrix and noise excitation matrix are as follows:
Φ = 1 T 0 0 0 1 0 0 0 0 1 T 0 0 0 1 and Γ = T 2 / 2 0 T 0 0 T 2 / 2 0 T ,
where T = 1 is the time interval, and the measurement equation is
Z k = ( x k o x ) 2 + ( k k o y ) 2 + v k ,
where ( o x , o y ) = (0, 1000) is land station, the variance of r k is R = 5 , and v k is tuning parameter satisfying v 1 .
Next, we will analyze and compare estimation performance of the following nonlinear filters that are conventional unscented Kalman filter (UKF), adaptive Huber unscented Kalman filter (AHUKF) [24], and the proposed filtering algorithm (AMUKF). Those three filters are utilized to track position and velocity of the maneuvering vehicle.
The initial position and velocity of the maneuvering vehicle are set as [0 m, 1400 m] and [2 m/s, −10 m/s], respectively. The value of initial variance P 0 is d i a g [ 1 , 1 , 1 , 1 ] . Simulation time lasts 50 s, the parameter σ of Gaussian kernel function is set as 0.8 here. In this simulation, two cases are considered: (1) measurement outliers and (2) state and measurement outliers simultaneously.
For the first case, we set the measurement outliers Z ^ 20 = Z 20 + 2 (m) and Z ^ 35 = Z 35 1 (m). Under the same condition, Figure 1, Figure 2 and Figure 3 show the tracking performance of UKF, AHUKF and the proposed filter (AMUKF). Figure 1 shows the tracking results of the three filters above, and the effects of AMUKF are better than that of other filters. The performance of AHUKF and AMUKF is much better than that of UKF because the UKF is a nonlinear extension of the KF, and it is susceptible to outliers. In Figure 1, we can see that AHUKF has a robust property to some extent, and it can be seen that the impact of state outlier is to be suppressed to a certain extent. However, comparing with AMUKF, it is not better than that of AMUKF. In summary, we can draw the conclusion that AMUKF can suppress the outliers and has better robust properties. Thus, it has a better tracking result for true state. The tracking errors obtained by AHUKF are bigger than that obtained by AMUKF. Figure 2 and Figure 3 show the position errors and velocity errors about the three filters, respectively. In contrast, under the measurement outlier, the estimation position and velocity errors for AMUKF are always smaller than that of other filters.
For the second case, we also estimate position and velocity like the first case. The state outliers are set as X ^ 10 = X 10 + 0.2 × [5,1,5,1] T (m), X ^ 25 = X 25 − 0.4 × [5,1,5,1] T (m), and the measurement outliers are set as Z ^ 20 = Z 20 + 2 (m) , Z ^ 35 = Z 35 5 (m), respectively. From Figure 4, Figure 5 and Figure 6, we can obtain that the effect of the three different filters when the process and measurement have outliers simultaneously. From Figure 4, it can be seen that the tracking results of the three filters in case 2, and the performance of AMUKF is better than that of other filters. The results of AHUKF and AMUKF are much better than that of UKF because AHUKF has robust properties, the effectiveness of AHUKF is approaching that of AMUKF, but it is not better than that of AMUKF. The tracking errors obtained by AHUKF are bigger than that obtained by AMUKF. Figure 5 and Figure 6 show the position errors and velocity errors about the three filters, respectively. Finally, it can be concluded that the errors of AMUKF are smaller than that of AHUKF. From Figure 1 and Figure 4, we also see that the behaviours of three filters are convergent. Finally, the root mean square error (RMSE) of state is shown in Table 2, which shows that the presented filtering algorithm outperforms both UKF and AHUKF.

4.2. Mars Entry Model

In this subsection, the Mars entry model is considered. During Mars entry, the vehicle is affected by some uncertainty disturbances such as lift, drag, gravity, and measurement outliers, etc. Measurement outliers are investigated. The dynamics of the entry vehicle are described by [30]
r ˙ v ˙ = v D v v L cos σ ˜ v v × v × r v × r + L sin σ ˜ v × r v × r g r r ,
where r = [ r x , r y , r z ] T and v = [ v x , v y , v z ] T denote position vector and velocity vector, and σ ˜ is bank angle. The lift acceleration L, drag acceleration D and gravitational acceleration g are
L = 1 2 ρ ˜ v 2 s / C L M D = 1 2 ρ ˜ v 2 s / C D M g = g ˜ R 2 ,
where r and v are the scalar values of r and v , M is the mass of the vehicle, and g ˜ is the planet’s gravitational constant. Atmosphere density ρ ˜ is given by
ρ ˜ = ρ ˜ 0 exp ( r 0 r h s ) ,
where ρ ˜ 0 = 2 × 10 4 kg/m 3 is the nominal reference density, r 0 = 3 . 4372 × 10 6 m is the reference radius of Mars (40 km above Mars surface), and h s = 7500 m is the constant scale height.
For the measurement model, the two Mars orbiters–Mars Reconnaissance Orbiter (MRO) and Mars Express (MEX) used by Curiosity for relay communications are employed. The movement of the orbiter in polar coordinate can be constructed [31]:
r o b = a ( 1 e 2 ) 1 + e c o s ( θ ω ¯ ) r o b 2 θ ˙ = g ˜ a ( 1 e 2 ) ,
where r o b is the distance between the orbiter and the center of Mars, θ is the angle between the ascending node and the orbiter, a is the semi-major axis of the movement ellipse, e is the eccentricity ratio, ω ¯ is the argument of perigee, and g ˜ is the gravity constant of Mars. The ecliptic longitude and latitude of the orbiter are:
Ω o b = cos 1 [ cos θ cos ϕ o b ] + Ω , ϕ o b 0 , cos 1 [ cos θ cos ϕ o b ] + Ω , ϕ o b > 0 ,
where ϕ o b = sin 1 [ sin θ sin i sin k ] , Ω o b is ecliptic longitude, and ϕ o b is the ecliptic latitude, and i is orbital inclination (angle between longitude and latitude, k = π 2 ). Thus, the relative range measurement between the vehicle and the orbiter can be constructed:
R o b = ( r r o b ) T ( r r o b ) , R ˜ o b = R o b + ε o b ,
where r o b = ( cos ϕ o b cos Ω o b cos ϕ o b sin Ω o b sin ϕ o b ) T , and ε o b is the range measurement noise, whose prior covariance is set to be 100, meaning the range measurement error is 10 m.
The relative range measurements between Mars surface beacons (MSBs) and the vehicle can be obtained by:
R ˜ B i ^ = R B i ^ + ε R i ^ ,
where R B i ^ = ( x B i ^ r x ) 2 + ( y B i ^ r y ) 2 + ( z B i ^ r z ) 2 refers to the true range from the entry vehicle to the ith beacon, R B ˜ is the measurement range, and x B i ^ , y B i ^ and z B i ^ are the position coordinates of MSBs. ε R i ^ is the range measurement noise.
The observation of rate measurement can be given as:
V ˜ B i ^ = V i ^ + ε V i ^ V B i ^ = d R i ^ d t i ^ = 1 , , N ,
where V B i ^ is the true rate, V ˜ B i ^ is the measured rate, and ε V i ^ is the rate measurement noise, whose prior covariance is set to be 1, meaning that the measurement error rate is 1 m/s.
From the above measurement models, the overall observation model can be obtained as
y = R ˜ o b R ˜ B V ˜ B .
In this simulation, the measurement outliers are considered. We utilize UKF and AHUKF [24] to compare with the proposed filtering algorithm too. Their performance is to be analyzed as follows. First, the initial values are listed in Table 3. The parameters of the two orbiters are listed in Table 4. The simulation time is 500 s.
The initial position and velocity of maneuvering vehicle are also showed in Table 3. The initial errors matrix P 0 = d i a g ( [ 1000 ; 1000 ; 1000 ; 1000 ; 10 ; 10 ] ) . Q = 10 5 × d i a g ( [ 0.1 ; 0.1 ; 0.1 ; 0.1 ; 0.1 ; 0.1 ] ) ; R = d i a g ( [ 100 ; 100 ; 100 ; 1 ; 1 ; 1 ] ) . The parameter σ of Gaussian kernel function is set as 10 here. The measurement outliers are assumed as y ˜ 100 = y 100 + 0.3 × [ 800 , 800 , 800 , 800 , 30 , 30 ] T and y ˜ 200 = y 200 0.4 [ 800 , 800 , 800 , 800 , 30 , 30 ] T , respectively. From Figure 7, Figure 8 and Figure 9, obviously, it can be seen that the state errors of UKF become very large for measurement outliers. while others’ errors are very small. Therefore, we can determine that the performance of AHUKF and AMUKF is far better than that of UKF, especially AMUKF. Even when the measurement outliers occur, AMUKF can effectively track the movement of the vehicle. Finally, from Figure 7, Figure 8 and Figure 9, we also see that the behaviours of three filters are convergent.

5. Conclusions

In this paper, a new adaptive robust unscented Kalman filter is obtained by combining both the robust property of maximum correntropy criterion and the adaptive property of the fading factor with the superiority of UKF. The cost function of the proposed filter includes both the fading factor to the state model and the maximum correntropy criterion to the measurement model. The maximum correntropy criterion can enhance the ability of being insensitive to outliers, and the fading factor can improve the ability of strongly tracking the state estimate. Therefore, the proposed filtering algorithm can track the state strongly and is insensitive to outliers. After that, the process of fading factor is derived in detail. Two numerical examples are given to illustrate the effectiveness of the proposed filter. UKF and AHUKF are selected for comparison with the proposed filtering algorithm in simulation. When state or/and measurement outliers occur, the proposed filter can track the target effectively and suppress the effect of outliers. Therefore, from the theoretical and numerical simulation results, it can be concluded that the proposed filter has not only the robustness of the fading factor but also the accuracy and flexibility of nonlinear systems. In the future, we will study the stability of the adaptive robust unscented Kalman filter based on the maximum correntropy criterion.

Author Contributions

L.Y. contributed with the concept and design of the manuscript, simulations, drafting and revising of the manuscript. Z.D. contributed with the design of the manuscript, drafting and revising the manuscript. B.H. and Y.X. contributed with revising of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China under Grant 61673067.

Acknowledgments

The authors will express their gratitude to the editors and reviewers for their valuable and constructive comments that can help to enhance the quality of this work.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Yu, Y. Consensus-Based Distributed Mixture Kalman Filter for Maneuvering Target Tracking in Wireless Sensor Networks. IEEE Trans. Veh. Technol. 2016, 65, 8669–8681. [Google Scholar] [CrossRef]
  2. Kumru, F.; Aka, T.; Altunta, E. Ballistic Target Tracking Using Multiple Model Kalman Filter with a Priori Ballistic Information. In Proceedings of the 25th Signal Processing and Communications Applications Conference, Antalya, Turkey, 15–18 May 2017. [Google Scholar]
  3. Yang, T.; Prashant, G. Probabilistic Data AssociationFeedback Particle Filter for Multiple Target Tracking Applications. J. Dyn. Syst. Meas. Control 2018, 140. [Google Scholar] [CrossRef]
  4. Kalman, R.E. A New Approach to Linear Filtering and Prediction Problems. Trans. ASME J. Basic Eng. 1960, 81, 35–45. [Google Scholar] [CrossRef]
  5. Lim, C.K.; Mba, D. Switching Kalman Filter for Failure Prognostic. Mech. Syst. Signal Process. 2015, 52–53, 1–6. [Google Scholar] [CrossRef]
  6. Feng, Z.G.; Teo, K.L.; Ahmed, N.U.; Zhao, Y.; Yan, W.Y. Optimal Fusion of Sensor Data for Kalman Filtering. Discret. Contin. Dyn. Syst. Ser. A 2017, 14, 483–503. [Google Scholar]
  7. Psiaki, M. Kalman Filtering and Smoothing to Estimate Real Valued States and Integer Constants. J. Guidance Control. Dyn. 2010, 33, 1404–1417. [Google Scholar] [CrossRef]
  8. Baroni, L. Kalman Filter for Attitude Determination of a CubeSat Using Low-cost Sensors. Comput. Appl. Math. 2017, 7–8, 1–12. [Google Scholar] [CrossRef]
  9. Sun, S.; Deng, Z. Multi-sensor Optimal Information Fusion Kalman Filter. Automatica 2004, 40, 1017–1023. [Google Scholar] [CrossRef]
  10. Dall’ Anese, E.; Bazerque, J.A.; Giannakis, G.B. Group Sparse Lasso For Cognitive Network Sensing Robust to Model Uncertainties and Outliers. Phys. Commun. 2012, 5, 161–172. [Google Scholar] [CrossRef]
  11. Julier, S.; Uhlmann, J. Unscented Filtering and Nonlinear Estimation. Proc. IEEE 2004, 92, 401–422. [Google Scholar] [CrossRef] [Green Version]
  12. Zhan, R.; Wan, J. Iterated Unscented Kalman Filter for Passive Target Tracking. IEEE Trans. Aerosp. Electron. Syst. 2007, 43, 1155–1163. [Google Scholar] [CrossRef] [Green Version]
  13. Arasaratnam, I.; Haykin, S. Cubature Kalman Filters. IEEE Trans. Autom. Control 2009, 54, 1254–1269. [Google Scholar] [CrossRef] [Green Version]
  14. Liu, H.; Wu, W. Strong Tracking Spherical Simplex-Radial Cubature Kalman Filter for Maneuvering Target Tracking. Sensors 2017, 17, 741. [Google Scholar] [CrossRef] [PubMed]
  15. Terejanu, G.; Singla, P.; Singh, T.; Scott, P. Adaptive Gaussian Sum Filter for Nonlinear Bayesian Estimation. IEEE Trans. Autom. Control 2011, 56, 2151–2156. [Google Scholar] [CrossRef]
  16. Du, H.; Wang, W.; Bai, L. Observation Noise Modeling Based Particle Filter: An Efficient Algorithm for Target Tracking in Glint Noise Environment. Neurocomputing 2015, 56, 155–166. [Google Scholar] [CrossRef]
  17. Wang, X.; Cui, N.; Guo, J. Huber-Based Unscented Filtering and Its Application to Vision-Based Relative Navigation. IET Radar Sonar Navig. 2010, 4, 134–141. [Google Scholar] [CrossRef]
  18. Durovic, Z.; Kovacevic, B. Robust Estimation with Unknown Noise Statistics. IEEE Trans. Autom. Control 1999, 44, 1292–1296. [Google Scholar] [CrossRef]
  19. Wang, Y.; Zheng, W.; Sun, S. Robust Information Filter Based on Maximum Correntropy Criterion. J. Guidance Control Dyn. 2016, 39, 1–6. [Google Scholar] [CrossRef]
  20. Karlgaard, C.; Schaub, H. Adaptive Nonlinear Huber-Based Navigation for Rendezvous in Elliptical Orbit. J. Guidance Control Dyn. 2011, 34, 388–402. [Google Scholar] [CrossRef]
  21. Gandhi, M.; Mili, L. Robust Kalman Filter Based on a Generalized Maximum-Likelihood-Type Estimator. IEEE Trans. Signal Process. 2010, 58, 2509–2520. [Google Scholar] [CrossRef]
  22. Chang, L.; Hu, B.; Chang, G.; Li, A. Multiple Outliers Suppression Derivative-Free Filter Based on Unscented Transformation. J. Guidance Control Dyn. 2015, 35, 1902–1906. [Google Scholar] [CrossRef]
  23. Coutinho, D.; de Souza, C.; Barbosa, K.; Trofino, A. Robust Linear H Filter Design for a Class of Uncertain Nonlinear Systems: An LMI approach. SIAM J. Control. Optim. 2009, 48, 1452–1472. [Google Scholar] [CrossRef]
  24. Wang, Y.D.; Zheng, W.; Sun, S. Adaptively Robust Unscented Kalman Filter for Tracking a Maneuvering Vehicle. J. Guidance Control Dyn. 2014, 37, 1696–1701. [Google Scholar] [CrossRef]
  25. Yang, Y.; He, H.; Xu, G. Adaptively Robust Filtering for Kinematic Geodetic Positioning. J. Geod. 2001, 75, 109–116. [Google Scholar] [CrossRef]
  26. Safarinejadian, B.; Yousefi, M. Static Alignment of Inertial Navigation Systems Using an Adaptive Multiple Fading Factors Kalman Filter. Syst. Sci. Control. Eng. 2015, 3, 351–359. [Google Scholar] [CrossRef]
  27. Geng, Y.; Wang, J. Adaptive Estimation of Multiple Fading Factors in Kalman Filter for Navigation Applications. GPS Solut. 2008, 12, 273–279. [Google Scholar] [CrossRef]
  28. Zhou, D.; Frank, P. Strong Tracking Kalman Filtering of Nonlinear Time-Varying Stochastic Systems with Coloured Noise: Application to Parameter Estimation and Empirical Robustness Analysis. Int. J. Control 1996, 65, 295–307. [Google Scholar] [CrossRef]
  29. Wu, Y.; Hu, D.; Wu, M.; Hu, X. A Numerical Integration Perspective on Gaussian Filters. IEEE Trans. Signal Process. 2006, 54, 2910–2921. [Google Scholar] [CrossRef]
  30. Yu, Z.S.; Cui, P.Y.; Zhu, S.Y. On the observability of Mars entry navigation using radiometric measurements. Adv. Space Res. 2014, 54, 1513–1524. [Google Scholar] [CrossRef]
  31. Wang, X.C.; Xia, Y.Q. Navigation Strategy with the Spacecraft Communications Blackout for Mars Entry. Adv. Space Res. 2014, 55, 1264–1277. [Google Scholar] [CrossRef]
Figure 1. Target tracking performance in case 1.
Figure 1. Target tracking performance in case 1.
Sensors 18 02406 g001
Figure 2. Position tracking performance in case 1.
Figure 2. Position tracking performance in case 1.
Sensors 18 02406 g002
Figure 3. Velocity tracking performance in case 1.
Figure 3. Velocity tracking performance in case 1.
Sensors 18 02406 g003
Figure 4. Target tracking performance in case 2.
Figure 4. Target tracking performance in case 2.
Sensors 18 02406 g004
Figure 5. Position tracking performance in case 2.
Figure 5. Position tracking performance in case 2.
Sensors 18 02406 g005
Figure 6. Velocity tracking performance in case 2.
Figure 6. Velocity tracking performance in case 2.
Sensors 18 02406 g006
Figure 7. Position and velocity tracking errors on the x-axis.
Figure 7. Position and velocity tracking errors on the x-axis.
Sensors 18 02406 g007
Figure 8. Position and velocity tracking errors on the y-axis.
Figure 8. Position and velocity tracking errors on the y-axis.
Sensors 18 02406 g008
Figure 9. Position and velocity tracking errors on the z-axis.
Figure 9. Position and velocity tracking errors on the z-axis.
Sensors 18 02406 g009
Table 1. AMUKF. AMUKF is short for Adaptive Robust Unscented Kalman Filter.
Table 1. AMUKF. AMUKF is short for Adaptive Robust Unscented Kalman Filter.
First set x 0 , k | k = [0 m, 1400 m, 2 m/s, −10 m/s] T .
P 0 , k | k = d i a g [ 1 , 1 , 1 , 1 ] .
Then iterate the follow,
for i = 0 : N 1
Reformulate the augmented covariance
time update
χ i , k | k 1 = f ( χ i , k 1 )
x ^ k | k 1 = i = 0 2 n ω i m χ i , k | k 1
Measurement update:
x ^ i , k | k 1 = [ χ i , k | k 1 χ i , k | k 1 + μ ( P i , k | k 1 ) χ i , k | k 1 μ ( P i , k | k 1 ) ]
y i , k | k 1 = h ( x ^ i , k | k 1 )
y ^ k | k 1 = Σ i = 0 2 n ω k y i , k | k 1
Calculate α k
v k = y k y ^ k | k 1
C 0 , k = v k v k T k = 1 λ C 0 , k 1 + v k v k T 1 + λ k > 1
α k = α 0 α 0 > 1 1 α 0 1
P k | k 1 = α k tr i = 0 2 n ω i c ( E i ) ( E i ) T + Q k
P y k y k = α k i = 0 2 n ω i c ( y i , k | k 1 y ^ k | k 1 ) ( y i , k | k 1 y ^ k | k 1 ) T + R ˜ k
P x k y k = α k i = 0 2 n ω i c ( χ i , k | k 1 y ^ k | k 1 ) ( y i , k | k 1 y ^ k | k 1 ) T
K k = P x k y k P y k y k 1 x ^ k | k = x ^ k | k 1 + K k ( y k y ^ k | k 1 )
P k | k = P k | k 1 K k P y k y k K k T
end
where E i = χ i , k | k 1 x ^ k | k 1 , μ = n + κ ,
W 0 = κ / ( n + κ ) , W i = κ / ( n + κ ) ,
i = 1 , 2 , , 2 n and κ is a tune parameter.
More details about the selection of κ can be seen in [11,22].
Table 2. RMSE of State. RMSE is short for root mean square error.
Table 2. RMSE of State. RMSE is short for root mean square error.
FilterRMSE of xRMSE of x ˙ RMSE of yRMSE of y ˙
UKF28.10.3151450.117
AHUKF26.70.0171390.037
AMUKF25.80.0111350.029
Table 3. Parameters setting of initial conditions. MSBs is short for Mars surface beacons.
Table 3. Parameters setting of initial conditions. MSBs is short for Mars surface beacons.
Initial SettingNotationValues
Initial position ( r x 0 , r y 0 , r z 0 ) (−3.92 km, −3099.09 km, −1663.11 km)
Initial velocity ( v x 0 , v y 0 , v z 0 ) (463.25 m/s, −1528.75 m/s, 5268.14 m/s)
MSBs’ locations (1) ( x B 1 , y B 1 , z B 1 ) (875.35 km, −2914.43 km, −1509.77 km)
MSBs’ locations (2) ( x B 2 , y B 2 , z B 2 ) (410.25 km, −2955.32 km, −1624.04 km)
Vehicle massM2804 kg
Vehicle cross-sections15.9 m 2
Table 4. Parameters of the orbiters. MRO and MEX are short for Mars Reconnaissance Orbiter and Mars Express, respectively.
Table 4. Parameters of the orbiters. MRO and MEX are short for Mars Reconnaissance Orbiter and Mars Express, respectively.
MROMEX
semi-major axis a 3663.7 km 8572.2 km
eccentricity ratio e 0.0089 rad 0.5770 rad
argument of perigee ω ¯ 4.7124 rad 2.7911 rad
orbital inclination i 1.6154 rad 1.5055 rad

Share and Cite

MDPI and ACS Style

Deng, Z.; Yin, L.; Huo, B.; Xia, Y. Adaptive Robust Unscented Kalman Filter via Fading Factor and Maximum Correntropy Criterion. Sensors 2018, 18, 2406. https://doi.org/10.3390/s18082406

AMA Style

Deng Z, Yin L, Huo B, Xia Y. Adaptive Robust Unscented Kalman Filter via Fading Factor and Maximum Correntropy Criterion. Sensors. 2018; 18(8):2406. https://doi.org/10.3390/s18082406

Chicago/Turabian Style

Deng, Zhihong, Lijian Yin, Baoyu Huo, and Yuanqing Xia. 2018. "Adaptive Robust Unscented Kalman Filter via Fading Factor and Maximum Correntropy Criterion" Sensors 18, no. 8: 2406. https://doi.org/10.3390/s18082406

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop