Next Article in Journal
Structural Health Monitoring for a Z-Type Special Vehicle
Previous Article in Journal
Nanoporous Gold Films Prepared by a Combination of Sputtering and Dealloying for Trace Detection of Benzo[a]pyrene Based on Surface Plasmon Resonance Spectroscopy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Estimation of Multiple Fading Factors for GPS/INS Integrated Navigation Systems

1
School of Environment Science and Spatial Informatics, China University of Mining and Technology, Xuzhou 221116, China
2
Collaborative Innovation Center for Resource Utilization and Ecological Restoration of Old Industrial Base, China University of Mining and Technology, Xuzhou 221116, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(6), 1254; https://doi.org/10.3390/s17061254
Submission received: 12 April 2017 / Revised: 25 May 2017 / Accepted: 27 May 2017 / Published: 1 June 2017
(This article belongs to the Section Remote Sensors)

Abstract

:
The Kalman filter has been widely applied in the field of dynamic navigation and positioning. However, its performance will be degraded in the presence of significant model errors and uncertain interferences. In the literature, the fading filter was proposed to control the influences of the model errors, and the H-infinity filter can be adopted to address the uncertainties by minimizing the estimation error in the worst case. In this paper, a new multiple fading factor, suitable for the Global Positioning System (GPS) and the Inertial Navigation System (INS) integrated navigation system, is proposed based on the optimization of the filter, and a comprehensive filtering algorithm is constructed by integrating the advantages of the H-infinity filter and the proposed multiple fading filter. Measurement data of the GPS/INS integrated navigation system are collected under actual conditions. Stability and robustness of the proposed filtering algorithm are tested with various experiments and contrastive analysis are performed with the measurement data. Results demonstrate that both the filter divergence and the influences of outliers are restrained effectively with the proposed filtering algorithm, and precision of the filtering results are improved simultaneously.

1. Introduction

With a full solution for position, velocity, and attitude, the GPS and INS systems are both effective navigation techniques which have found broad applications. GPS, characterized by high precision, is able to provide accurate navigation and positioning information. INS is a self-navigation system which needs not to send and receive signals. The GPS/INS integrated navigation system manifests great superiority by integrating the advantages of GPS and INS, and performance of the integrated system is better than that of either the GPS or INS system alone [1]. Thus, the integration of GPS and INS gives a more complete and accurate navigation solution, and this integrated system has been widely adopted in the field of dynamic navigation and positioning in recent years. The Kalman filter is the most popular real-time optimal estimator [2], and it has been commonly applied in navigation applications. In the GPS/INS integrated navigation system, the Kalman filter is also the most commonly used data fusion method [3,4,5]. The parameter estimates of the Kalman filter are optimal when the statistical information of the noise is correctly specified. However, the required information is difficult to achieve when the model errors have time varying characteristics [6,7], and the performance will be degraded with improper statistical information and outliers, and sometimes may even lead to a filter divergence [7]. Moreover, it is unavoidable, especially when caused by a complex perturbation in a signal degraded environment [3]. As one of the basic procedures for the quality control of the Kalman filter, the Detection, Identification and Adaptation (DIA) methods were proposed to eliminate the effects of outliers [8]. In the DIA methods, the outliers are detected and identified firstly, then deviations of the state estimates caused by the outliers are eliminated, finally, the adaptation is implemented in which the model is recovered from the identified errors. However, the identification is quite difficult, especially when the measurements are not accurate [9]. For instance, the detection and identification process are difficult to implement when the measurements contain multiple outliers or the continuously-changing outliers. Therefore, in order to weaken the influence of the outliers, it is essential to research the filtering and estimation algorithms which are more robust in the GPS/INS integrated navigation system [9,10,11].
To achieve a better performance with the GPS/INS integrated navigation system, the modeling and estimation problem should be addressed properly [7]. The former is concerned with developing better error models and the latter requires proper use of the estimation algorithm and the measurement information. In the Kalman filter, filter divergence often happens when big errors of the system model exist, due to the effects of the improper a prior information [12]. With a constant fading factor, the fading Kalman filter was proposed to overcome the shortcomings of the Kalman filter [13]. In recent years, research about the fading Kalman filter focuses on the estimation of the fading factor. An optimal adaptive algorithm with a single fading factor for fading Kalman filter was proposed [14], nevertheless, a single fading factor can hardly guarantee the optimization of the whole system model, and the filter divergence will be likely to happen when handling complex system model [15]. Various adaptive Kalman filters have been developed to achieve a better estimator [9,16,17,18]. Then the multiple fading filter [15,19] was proposed to improve the performance of the fading filter. In the field of navigation and positioning, an algorithm based on the magnitude of the predicted residuals was proposed to reduce the dynamic model errors [20]. However, most of the above algorithms show few abilities to control the influences of the uncertain interferences. Moreover, the GPS/INS integrated navigation is very likely to be affected by uncertainties due to the exogenous disturbances, device damage, and inaccurate sensor noise statistical information [10]. By minimizing the estimation error in the worst case, the H-infinity filtering algorithm provides an effective way to improve the robustness of the system model [21] and it can be used to address the uncertainties in the measurement noise [18,22].
This paper focuses on a comprehensive filtering algorithm for the GPS/INS integrated navigation system by combining the advantages of the multiple fading filter and the H-infinity filter. Based on the concept of fading memory and the optimization of the filter, a new multiple fading factor is proposed. Experiments are implemented with the data collected under actual conditions. Performances of different filtering algorithms are evaluated and the contrastive analysis are implemented. Results demonstrate that the new multiple fading filter algorithm shows a higher accuracy whether the abnormal measurements ever exist. In addition, the proposed fading filtering algorithm is easy to implement with a relatively reasonable computation burden.
The rest of this paper is organized as follows: in Section 2, the theory of the fading filter and the H-infinity filter are introduced, and the nonlinear form of the H-infinity filter is listed. In Section 3, a new multiple fading factor is proposed and a comprehensive filtering algorithm is constructed. In Section 4, rules of the GPS/INS integrated navigation system are introduced, and equations of the system for the loosely-coupled integrated navigation system are listed. Section 5 presents the advantages of the proposed filtering algorithm through actual experiments, and the results are discussed and analyzed. The conclusions are provided in Section 6.

2. Fading Filter and H-Infinity Filter

In this section, the basic rules of the fading filter are introduced and the recent fading factors are discussed, then the principles of the H-infinity filter and its nonlinear form are discussed.

2.1. Basic Rules of the Fading Filter

For the dynamic model equation and the measurement equation:
{ x k = Φ k , k 1 x k 1 + w k z k = H k x k + v k ,
the predicted state vector is presented as:
x k / k 1 = Φ k , k 1 x k 1 ,
where x k denotes the state vector at the epoch k , Φ k , k 1 denotes the state transition matrix, H k denotes the measurement matrix, z k denotes the measurement vector, w k and v k denote the state noise and measurement noise, respectively, and x k / k 1 denotes the predicted state vector.
The recursion approach of the fading filter is given by:
x k / k = x k / k 1 + K ¯ k ( z k H k x k / k 1 ) ,
K ¯ k = P ¯ k / k 1 H k T ( H k P ¯ k / k 1 H k T + R k ) 1 ,
P ¯ k / k 1 = S k Φ k , k 1 P k 1 Φ T k , k 1 + Q k ,
where K ¯ k denotes the equivalent gain matrix, P ¯ k / k 1 denotes the equivalent covariance matrix of the state vector, S k denotes the fading factor and S k 1 , P k 1 denotes the covariance matrix of the state vector at epoch k 1 , R k and Q k denote the covariance matrices of the measurement noise and the state noise, respectively.
In the standard Kalman filter, the covariance matrix of the state vector is given by:
P k / k 1 = Φ k , k 1 P k 1 Φ T k , k 1 + Q k .
In the fading filter, the covariance matrix of the a prior state is inflated for S k times which degrades the efficiency of the past state information, and the recent measurement information is regarded more seriously [12]. Obviously, compared with the standard Kalman filter, the state errors from the previous epoch are well controlled.

2.2. The Recent Fading Factors

The key problem in the fading filter is to find a suitable fading factor. With the optimization considered, a fading filter algorithm was proposed [14] where two forms of fading factors were derived, namely, the steepest descent method and the one-step method. It is quite difficult to apply the former method due to its complex computations. In practice, the latter is more applicable which is derived by:
S k = max { 1 , 1 n t r ( N k M k 1 ) } ,
where:
M k = H k Φ k , k 1 P k 1 Φ k , k 1 T H k T ,
N k = P V k H k Q k H k T R k ,
P V k = E ( V k V k T ) ,
V k = H k x k / k 1 z k ,
where t r ( ) denotes taking trace of a matrix, V k denotes the predicted residual vector, and P V k denotes the covariance matrix of V k , and P ^ V k = 1 k i = 1 k V i V i T [12]. Then Xia [14] proposed a simplified method, namely:
S k = max { 1 , t r ( N k ) / t r ( M k ) } ,
the fading factor will increase along with V k , and S k calculated in this way is optimal in theory.
A multiple fading filter algorithm was proposed by Geng [19], and the fading factor is given by:
S k = diag ( s 1 , s 2 , ... , s t , 1 , ... , 1 , ... , 1 ) ,
s i = max ( 1 , [ v i ( k ) ] 2 λ i 2 j i i ( k ) ε i b i i ( k ) j i i ( k ) ) ,   ( i = 1 , 2 , ... , t ) ,
where v i ( k ) is the i th element of V k , t is the dimension of the measurement equation, λ i is the i th observable element of H k , b i i ( k ) is the i th diagonal element of the matrix B k and B k = H k Q k 1 H k T + R k , j i i ( k ) is the i th diagonal element of the matrix J k and J k = Φ k P k 1 Φ k T , ε i is the threshold value fixed according to the Chi-square distribution at the given confident level. In this approach, only ( s 1 , s 2 , ... , s t ) can be estimated adaptively due to their observabilities, and the other elements should be fixed as one.

2.3. Principles of the H-Infinity Filter

Assume that the dynamic model equation and the measurement equation is given by:
{ x k = Φ k , k 1 x k 1 + w k y k = L k x k z k = H k x k + v k ,
where y k denotes the state vector to be estimated, L k denotes the given matrix, and L k is often fixed by an identity matrix. No statistical information of the noises is assumed before the filtering.
The cost function J is defined by [21]:
J = k = 1 N x k x ^ k 2 x 0 x ^ 0 P 0 1 2 + k = 1 N ( w k Q k 1 2 + v k R k 1 2 ) ,
where N denotes the total number of filtering time limit, x 0 denotes the initial value of the state vector x with the covariance matrix P 0 , x ^ 0 and x ^ k denote the estimated state vectors of x 0 and x k , respectively, and the expression x 0 x ^ 0 P 0 1 2 denotes ( x 0 x ^ 0 ) T P 0 1 ( x 0 x ^ 0 ) .
The estimate of the state vector should be achieved under the conditions of x ^ k = arg   min J . However, the closed-form approach of an optimal H-infinity filter can be achieved only under special conditions [23]. Accordingly, a suboptimal recursion algorithm was proposed through setting a threshold value γ which satisfied the following Riccati inequality [1]:
P k 1 + H k T H k γ 2 L k T L k > 0 ,
where P k is the covariance matrix of x k . Then, the recursion formulas of the H-infinity filter in the linear systems, which are different with the standard Kalman filter, are obtained [24]:
P k / k = P k / k 1 Φ k / k 1 P k / k 1 [ H k T   L k T ] R e , k 1 [ H k L k ] P k / k 1 Φ k / k 1 T ,
R e , k 1 = [ I 0 0 γ 2 I ] + [ H k L k ] P k / k 1 [ H k T   L k T ] .
It should be noticed that the robustness of the filter is related to the threshold value γ , and the filter becomes more robust when γ decreases. However, it may cause a filter divergence if γ is too close to zero [23].

3. Adaptive Estimation of the Multiple Fading Factor

In this section, the adaptive estimation of the new multiple fading factor is proposed, and a multiple fading filtering algorithm is constructed with the proposed multiple fading factor. Finally, the recursion formulas of the new algorithm are derived.
As mentioned above, the performance of the Kalman filter may degrade due to the past information. A constant fading factor was thus designed to address this problem [14]. However, a constant fading factor cannot guarantee optimal filtering absolutely. A multiple fading factor was proposed to improve the performance of a single fading factor [19], whereas it is unsatisfactory with respect to the robustness of the filtering. Accordingly, a new multiple fading filter is proposed to further improve the performance of the fading filter.
For the a prior covariance matrix in the multiple fading filter:
P k / k 1 = S k Φ k , k 1 P k 1 Φ T k , k 1 + Q k ,
S k = diag ( s 1 , s 2 , ... , s t , 1 , 1 , ... , 1 ) . To keep the symmetry of P k / k 1 , the Equation (20) can be rewritten as:
P k / k 1 = S k Φ k , k 1 P k 1 Φ T k , k 1 S k T + Q k .
Based on the observability and the optimal fading factor, a new multiple fading factor is constructed as:
s i = { max { 1 , [ N k ( i , i ) ] / [ M k ( i , i ) ] } , i = 1 , 2 , ... , t 1 ,   other ,
where N k ( i , i ) and M k ( i , i ) represent the i th diagonal elements of N k and M k , respectively. Similarly, only ( s 1 , s 2 , ... , s t ) can be estimated adaptively due to their observabilities, and the other elements should be fixed as one. The multiple fading factors in Equation (22) are constructed based on the optimal single fading factor in Equation (12) and they can guarantee the optimization of different data channels, thus the optimization of the whole filter is guaranteed.
Then, a new multiple fading H-infinity filtering algorithm is constructed. With the H-infinity filter adopted, the recursion formulas of the new algorithm for a linear system are listed below (the predicted state vector is calculated in the same way as Equation (2)):
P ¯ k / k 1 = S k Φ k , k 1 P k 1 / k 1 Φ k , k 1 T S k T + Q k ,
x k / k = x k / k 1 + K ¯ k ( z k H k x k / k 1 ) ,
K ¯ k = P ¯ k / k 1 H k T ( H k P ¯ k / k 1 H k T + R k ) 1 ,
P ¯ k / k = P ¯ k / k 1 Φ k / k 1 P ¯ k / k 1 [ H k T   L k T ] R ¯ e , k 1 [ H k L k ] P ¯ k / k 1 Φ k / k 1 T ,
R ¯ e , k 1 = [ I 0 0 γ 2 I ] + [ H k L k ] P ¯ k / k 1 [ H k T   L k T ] ,
where P ¯ k / k denotes the a posteriori covariance matrix of the state, and S k is fixed through (22).

4. The GPS/INS Integrated Navigation System

In this section, the dynamic model and the state estimator are listed. Then the loosely-coupled GPS/INS integrated navigation system is discussed in detail.
In recent years, three main forms of integration are developed: the loosely-coupled, the tightly-coupled and the deeply-coupled. Compared with the other two forms of integration, the loosely-coupled navigation system is easier to be implemented and the computation process is more concise [25,26].
A 15-dimension state vector is designed in the loosely-coupled GPS/INS integrated navigation system in this paper. And the parameters denote the deviations of position, velocities, attitudes, and the random bias of the gyroscope and accelerometer, respectively. The state vector X ^ is given by:
X ^ = [ δ x , δ y , δ z , δ v x , δ v y , δ v z , δ φ e , δ φ n , δ φ u , δ g x , δ g y , δ g z , δ a x , δ a y , δ a z ] .
The nonlinear differential error model for the INS system is defined by [6,27]:
{ Δ R ˙ e = Δ V e Δ V ˙ e = ( I 3 × 3 C e e ) f e + C b e b 2 Ω i e e Δ V e φ ˙ e = ( I 3 × 3 C e e ) ω i e e C b e ε b ˙ b = 0 ε ˙ b = 0 ,
where Δ R e and Δ V e are the position and velocity errors under the computer e frame, respectively, φ ˙ e is the attitude error between e frame and the platform e frame, I 3 × 3 is the identity matrix, C e e is the rotation matrix between e and e frame, C b e is the rotation matrix between the body frame and e frame, Ω i e e is the skew symmetric matrix of earth rotation rate ω i e e , b and ε b are the gyroscope and accelerometer drift errors under the body frame, respectively.
In this paper, the cubature Kalman filter is adopted to address the nonlinear problem of the integrated navigation system. For a discrete nonlinear system:
{ x k = f ( x k 1 ) + w k z k = h ( x k ) + v k ,
where f ( ) and h ( ) denote the known nonlinear functions. The cubature points are generated by:
ξ = m 2 [ 1 ] i ,   i = 1 , ... , m ,
where m denotes the number of the cubature points and m = 2 n , n denotes the dimension of the state vector, [ 1 ] denotes the integral and symmetrical set of points. Then equations of the cubature Kalman filter are given by [28].
(i) Time update
{ s k 1 / k 1 = S V D ( P k 1 / k 1 ) X k 1 , k 1 = s k 1 / k 1 ξ + x k 1 / k 1 X k / k 1 = f ( X k 1 , k 1 ) ,
x k / k 1 = 1 m i = 1 m X i , k / k 1 ,
P k / k 1 = 1 m i = 1 m X i , k / k 1 X i , k / k 1 T x k / k 1 x k / k 1 T + Q k .
(ii) Measurement update
{ s k / k 1 = S V D ( P k / k 1 ) X k / k 1 = s k / k 1 ξ + x k / k 1 Z k / k 1 = h ( X i , k / k 1 ) z k / k 1 = 1 m i = 1 m Z i , k / k 1 P z z , k / k 1 = 1 m i = 1 m Z i , k / k 1 Z i , k / k 1 T z k / k 1 z k / k 1 T + R k P x z , k / k 1 = 1 m i = 1 m X i , k / k 1 Z i , k / k 1 T x k / k 1 z k / k 1 T ,
then the final measurement update equations are obtained:
K k = P x z , k / k 1 P z z , k / k 1 1 ,
x k / k = x k / k 1 + K k ( z k z k / k 1 ) ,
P k / k = P k / k 1 K k P z z , k / k 1 K k T ,
where s denotes the square root of the covariance matrix P , X k 1 , k 1 denotes the cubature points of the states vector, X k / k 1 denotes the propagated cubature points, S V D denotes the singular value decomposition of a matrix; Z k / k 1 denotes the propagated cubature points of the measurement vector. Accordingly, the multiple fading filtering algorithm for the GPS/INS integrated navigation system is listed below:
(i) Time update
{ x k / k 1 = 1 m i = 1 m X i , k / k 1 P k / k 1 = S k ( 1 m i = 1 m X i , k / k 1 X i , k / k 1 T x k / k 1 x k / k 1 T ) S k T + Q k ,
(ii) Measurement update
{ x k / k = x k / k 1 + K k ( z k z k / k 1 ) K k = P x z , k / k 1 P z z , k / k 1 1 P k / k = P k / k 1 [ P x z , k / k 1 P k / k 1 ] [ P z z , k / k 1 R k + I P x z , k / k 1 T P x z , k / k 1 P k / k 1 γ 2 I ] 1 [ P x z , k / k 1 T P k / k 1 T ] .
In the loosely-coupled GPS/INS integrated navigation system, the differences of position and velocity between GPS and INS are regarded as the measurement inputs to the integration algorithm, namely:
Z ρ ( t ) = ρ GPS ρ INS ,
where ρ GPS and ρ INS are the position and velocity information output by the GPS and INS, respectively. Then the measurement vector is presented as:
Z k = [ r GPS r INS v GPS v INS ] ,
where Z k denotes the measurement vector of the integrated navigation system at epoch k , r GPS , r INS , v GPS , and v INS are the position and velocity information of GPS and INS, respectively. In the loosely-coupled navigation system, the measurement equation is linear, then the Equation (40) can be rewritten as:
{ x k / k = x k / k 1 + K k ( z k H k x k / k 1 ) K k = P k / k 1 H k T ( H k P k / k 1 H k T + R k ) 1 P k / k = P k / k 1 [ P k / k 1 H k T P k / k 1 ] [ H k P k / k 1 H k T + I H k P k / k 1 T P k / k 1 H k T P k / k 1 γ 2 I ] 1 [ H k P k / k 1 T P k / k 1 T ] .
Obviously, the observable variables are the position and velocity information of each direction, thus the corresponding fading factors can be estimated adaptively.

5. Testing

In this section, two cases of a loosely-coupled GPS/INS integrated navigation system are designed and implemented. Then the performance of the proposed filtering algorithm is discussed in comparison with the existing algorithms.

5.1. Experiment Schemes

To examine the performance of the proposed multiple fading filter, various actual experiments were performed. The data were obtained through a vehicle mounted GPS/INS integrated navigation system. In this system, two Trimble GPS receivers and a low cost inertial measurement unit (IMU) were adopted, and the main parameters of the IMU are listed in Table 1.
The update rate of GPS and IMU were set by 1 Hz and 100 Hz, respectively. The initial position and velocity errors are 1.0 m and 0.1 m/s, respectively. The initial position errors are set by 3.0 m, 3.0 m, 5.0 m, respectively, the initial velocity error is 0.5 m/s, and the initial alignment errors are set by 3°, 3°, 10°, respectively. Double difference GPS data were adopted with the position and velocity variances of 0.25 m2 and 0.0025 m2/s2, respectively. The precise results calculated by the double difference carrier phase were regarded as references. The motion trajectory of the land vehicle is displayed in Figure 1.
Then, two cases were designed to evaluate the stability and robustness of the new algorithm: the original measurements were adopted in case 1, and the measurements with perturbations in case 2. In each case, six schemes were performed:
  • Scheme 1: the cubature Kalman Filter (CKF);
  • Scheme 2: the H-infinity Filter (HF-CKF) ( γ was set as 2);
  • Scheme 3: the conventional multiple fading filter (MF-CKF);
  • Scheme 4: the new multiple fading filter (NMF-CKF);
  • Scheme 5: the multiple fading H-infinity filter (MHF-CKF) ( γ was set as 2);
  • Scheme 6: the new multiple fading H-infinity filter (NMHF-CKF) ( γ was set as 2).

5.2. Results of the Experiments

(1) Case 1:
In this case, each algorithm was performed with the original measurement data. The position errors in the X, Y, Z directions of the six schemes are demonstrated in Figure 2, Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7:
The position and attitude RMSEs of the each scheme are listed in Table 2 and Table 3, respectively.
Big disturbances are found in Figure 2 which may occur when the vehicle is passing through the speed hump, and this indicates that the robustness of the CKF algorithm should be strengthened. It is concluded from Figure 2 and Figure 3 that, with the effects of the uncertainties controlled, the HF-CKF algorithm performs much better than the CKF algorithm. Compared with the CKF algorithm, Figure 4 and Figure 5 demonstrate that both the conventional and the new multiple fading filters manifest certain ability of controlling the influence of the model deviations. By integrating the advantages of the multiple fading filter and the H-infinity filter, the MHF-CKF and NMHF-CKF algorithms show better performance than the other algorithms which can be learned from the Figure 6 and Figure 7, and this indicates that the influences of the model deviations and uncertain interferences are well controlled with the MHF-CKF and NMHF-CKF algorithms. In Table 2, RMSEs of the MF-CKF and NMF-CKF algorithms are much smaller than those of the CKF algorithm, and it can be concluded that the fading factors have reduced the weight of the past state information. Consequently, compared with the conventional Kalman filter, the impacts of the state model errors are controlled effectively with the multiple fading filter algorithm, and thus a higher precision is obtained. Figure 7 and Table 2 demonstrate that the new multiple fading filter performs the best among the six schemes, and both the amplitudes of errors and the RMSEs are further reduced. As displayed in Table 3, there exists a little differences of the attitude errors between each algorithm, and the RMSEs of the NMHF-CKF algorithm are slightly smaller than those of the other algorithms.
(2) Case 2:
In this case, ability of robustness for the above schemes is examined with both scattered and continuously-changing outliers. Thus the scattered outliers at the 160th, 260th, 360th, 460th, 560th epochs, respectively, and the continuously-changing outliers were added artificially to the GPS measurements at the 351th to the 380th epochs. Then the six schemes were implemented with the new data. Position errors of these schemes are displayed in Figure 8, Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13:
In this case, the position errors are mainly caused by outliers. It is clearly demonstrated from the Figure 8 and Figure 9 that the CKF and HF-CKF algorithms show little ability to resist the outliers. On the contrary, both the conventional and the new fading filters manifest a noticeable robustness which can be seen from the Figure 10, Figure 11, Figure 12 and Figure 13. By integrating the advantages of the H-infinity filter, filtering precision is furtherly improved with the MHF-CKF and the NMHF-CKF algorithms. In addition, error amplitudes of the new multiple fading filter are smaller than those of the other filtering algorithms, which indicates that effects of both scattered and continuously-changing outliers are weakened effectively. Then the position and attitude RMSEs of each scheme are listed in Table 4 and Table 5:
In Table 4, position RMSEs of the CKF and the HF-CKF algorithms are significantly affected, and the HF-CKF algorithm exerts certain robustness compared with the CKF algorithm. Position RMSEs of the conventional and the new multiple fading filtering algorithms change little in the presence of outliers, and the new multiple fading filter shows a higher accuracy than the conventional multiple fading filter whether the H-infinity filter is ever adopted. The Table 5 indicates that the proposed algorithm performs better than the other algorithms, and the influence of the attitude errors are weakened further. Due to the fading factors, multiple data channels are adjusted depending on their observabilities of the state vector, and the state model errors from the previous epoch are well controlled, thus the filter becomes more stable and the filter divergence caused by the model deviations is restrained. It should be pointed out however, in practice, N k may sometimes become negative definite and S k has to be fixed by 1 now, and this indicates that the state model error may not be controlled effectively at the current epoch [12].

6. Conclusions

In this study, a new multiple fading factor is proposed based on the thought of fading memory and the optimization of the filter, and a comprehensive filtering algorithm is constructed by integrating the proposed multiple fading filter and the H-infinity filter. The new filtering algorithm is implemented in the loosely-coupled GPS/INS integrated navigation system. The detail conclusions in this paper are summarized below:
(1)
The H-infinity filter based on the cubature Kalman filter performs better than the cubature Kalman filter, however, they are both affected significantly by the model errors and outliers. The conventional multiple fading filter shows certain robustness which can be furtherly enhanced.
(2)
With the influences of the model deviations and the uncertain interferences controlled simultaneously, performance of the conventional multiple fading filter is further improved, and the higher precision is obtained by integrating the advantages of the multiple fading filter and the H-infinity filter. The proposed multiple fading filter is implemented in the loosely-coupled GPS/INS integrated navigation system, and the stability and robustness are demonstrated with the original data and the data with both the scattered and the continuously-changing outliers.

Acknowledgments

This work was supported by the Fundamental Research Funds for the Central Universities (2017BSCXB40).

Author Contributions

Chen Jiang prepared the manuscript; Shu-Bi Zhang contributed to the theory studies and designed the experiments; Qiu-Zhao Zhang examined the quality of this work. All authors reviewed the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Godha, S.; Cannon, M.E. GPS/MEMS INS integrated system for navigation in urban areas. GPS Solut. 2007, 11, 193–203. [Google Scholar] [CrossRef]
  2. Kalman, R.E. A new approach to linear filtering and prediction problems. J. Basic Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef]
  3. Godha, S.; Lachapelle, G.; Cannon, M.E. Integrated GPS/INS system for pedestrian navigation in a signal degraded environment. In Proceedings of the 19th International Technical Meetings of the Satellite Division of the Institute of Navigation, Fort Worth, TX, USA, 26–29 September 2006. [Google Scholar]
  4. Gao, S.; Zhong, Y.; Zhang, X.; Shirinzadeh, B. Multi-sensor optimal data fusion for INS/GPS/SAR integrated navigation system. Aerosp. Sci. Technol. 2009, 13, 232–237. [Google Scholar] [CrossRef]
  5. Caron, F.; Duflos, E.; Pomorski, D.; Vanheeghe, P. GPS/IMU data fusion using multisensor Kalman filtering: introduction of contextual aspects. Inf. Fusion 2006, 7, 221–230. [Google Scholar] [CrossRef]
  6. Hide, C.; Moore, T.; Smith, M. Adaptive Kalman Filtering for Low-cost INS/GPS. J. Navig. 2003, 56, 143–152. [Google Scholar] [CrossRef]
  7. Mohamed, A.H.; Schwarz, K.P. Adaptive Kalman filtering for INS/GPS. J. Geod. 1999, 73, 193–203. [Google Scholar] [CrossRef]
  8. Teunissen, P.J.G. Quality control in integrated navigation systems. In Proceedings of the 1990’—A Decade of Excellence in the Navigation Sciences, IEEE PLANS’90 Position Location and Navigation Symposium, Las Vegas, NV, USA, 20–20 March 1990; pp. 158–165. [Google Scholar]
  9. Yang, Y.X.; He, H.B.; Xu, G.C. Adaptively robust filtering for kinematic geodetic positioning. J. Geod. 2001, 75, 109–116. [Google Scholar] [CrossRef]
  10. Zhao, L.; Qiu, H.Y.; Feng, Y.M. Analysis of a robust Kalman filter in loosely coupled GPS/INS navigation system. Measurement 2016, 80, 138–147. [Google Scholar] [CrossRef]
  11. Wu, F.M.; Yang, Y.X. A new two-step adaptive robust Kalman filtering in GPS/INS integrated navigation system. Acta Geod. Cartogr. Sin. 2010, 39, 522–527. [Google Scholar]
  12. Yang, Y.X. Adaptive Navigation and Dynamic Positioning; Surveying and Mapping Press: Beijing, China, 2006. [Google Scholar]
  13. Fagin, S.L. Recursive linear regression theory: Optimal filter theory and error analysis. IEEE Int. Conv. Rec. 1964, 12, 216–245. [Google Scholar]
  14. Xia, Q.J.; Sun, Y.X.; Zhou, C.H. An optimal adaptive algorithm for fading Kalman filter and its application. Acta Autom. Sin. 1990, 16, 210–216. [Google Scholar]
  15. Qian, H.M.; Lei, G.E.; Yu, P. Multiple fading factors Kalman filter and its application in SINS initial alignment. J. Chin. Inert. Technol. 2012, 20, 287–291. [Google Scholar]
  16. Lee, T.S. Theory and application of adaptive fading memory Kalman filters. IEEE Trans. Circuits Syst. 1988, 35, 474–477. [Google Scholar] [CrossRef]
  17. Yang, Y.X.; He, H.B.; Xu, T.H. Adaptive robust filtering for GPS positioning. Acta Geod. Cartogr. Sin. 2001, 30, 293–298. [Google Scholar]
  18. Chang, G.B. Kalman filter with both adaptivity and robustness. J. Process Control 2014, 24, 81–87. [Google Scholar] [CrossRef]
  19. Geng, Y.R.; Wang, J.L. Adaptive estimation of multiple fading factors in Kalman filter for navigation applications. GPS Solut. 2008, 12, 273–279. [Google Scholar] [CrossRef]
  20. Hu, C.W.; Chen, W.; Chen, Y.Q.; Liu, D.J. Adaptive Kalman filtering for vehicle navigation. J. Glob. Position. Syst. 2003, 2, 42–47. [Google Scholar] [CrossRef]
  21. Simon, D. Optimal State Estimation: Kalman, H and Nonlinear Approaches; John Wiley and Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  22. Duan, Z.S.; Zhang, J.X.; Zhang, C.S.; Mosca, E. Robust H 2 and H filtering for uncertain linear systems. Automatica 2006, 42, 1919–1926. [Google Scholar] [CrossRef]
  23. Hassibi, C.; Sayed, A.H.; Kailath, T. Linear estimation in Krein spaces-part II: applications. IEEE Trans. Autom. Control 1996, 41, 34–49. [Google Scholar] [CrossRef]
  24. Yue, X.K.; Yuan, J.P. H-infinity filtering algorithm and its application in GPS/SINS integrated navigation system. Acta Aeron Autica Astronaut. Sin. 2001, 22, 366–368. [Google Scholar]
  25. Gao, W.G.; Yang, Y.X.; Zhang, S.C. Adaptive robust Kalman filtering based on the current statistical model. Acta Geod. Cartogr. Sin. 2006, 35, 15–18. [Google Scholar]
  26. Wang, X.M.; Ni, W.B. An improved particle filter and its application to an INS/GPS integrated navigation system in a serious noisy scenario. Meas. Sci. Technol. 2016, 27, 095005. [Google Scholar] [CrossRef]
  27. Zhang, Q.Z.; Meng, X.L.; Zhang, S.B.; Wang, Y.J. Singular value decomposition-based robust cubature Kalman filtering for an integrated GPS/SINS navigation system. J. Navig. 2015, 68, 549–562. [Google Scholar] [CrossRef]
  28. Zhao, Y.W. Performance evaluation of cubature Kalman filter in a GPS/IMU tightly-coupled navigation system. Signal Proc. 2016, 119, 67–79. [Google Scholar] [CrossRef]
Figure 1. The motion trajectory of the land vehicle.
Figure 1. The motion trajectory of the land vehicle.
Sensors 17 01254 g001
Figure 2. Position errors of the CKF algorithm.
Figure 2. Position errors of the CKF algorithm.
Sensors 17 01254 g002
Figure 3. Position errors of the HF-CKF algorithm.
Figure 3. Position errors of the HF-CKF algorithm.
Sensors 17 01254 g003
Figure 4. Position errors of the MF-CKF algorithm.
Figure 4. Position errors of the MF-CKF algorithm.
Sensors 17 01254 g004
Figure 5. Position errors of the NMF-CKF algorithm.
Figure 5. Position errors of the NMF-CKF algorithm.
Sensors 17 01254 g005
Figure 6. Position errors of the MHF-CKF algorithm.
Figure 6. Position errors of the MHF-CKF algorithm.
Sensors 17 01254 g006
Figure 7. Position errors of the NMHF-CKF algorithm.
Figure 7. Position errors of the NMHF-CKF algorithm.
Sensors 17 01254 g007
Figure 8. Position errors of the CKF algorithm.
Figure 8. Position errors of the CKF algorithm.
Sensors 17 01254 g008
Figure 9. Position errors of the HF-CKF algorithm.
Figure 9. Position errors of the HF-CKF algorithm.
Sensors 17 01254 g009
Figure 10. Position errors of the MF-CKF algorithm.
Figure 10. Position errors of the MF-CKF algorithm.
Sensors 17 01254 g010
Figure 11. Position errors of the NMF-CKF algorithm.
Figure 11. Position errors of the NMF-CKF algorithm.
Sensors 17 01254 g011
Figure 12. Position errors of the MHF-CKF algorithm.
Figure 12. Position errors of the MHF-CKF algorithm.
Sensors 17 01254 g012
Figure 13. Position errors of the NMHF-CKF algorithm.
Figure 13. Position errors of the NMHF-CKF algorithm.
Sensors 17 01254 g013
Table 1. Main technological parameters of IMU and GPS.
Table 1. Main technological parameters of IMU and GPS.
SensorsRandom BiasRandom Constant NoiseRMSEs
Gyroscope20 (°)/h0.067 (°)/h1/2-
Accelerometer5 mg50 μg/h1/2-
GPS receiver--Position: 0.5 m; Velocity: 0.05 m/s
Table 2. Position RMSEs of each scheme (meter).
Table 2. Position RMSEs of each scheme (meter).
AxisCKFHF-CKFMF-CKFNMF-CKFMHF-CKFNMHF-CKF
X 0.1290.0960.1080.0930.0950.088
Y 0.2840.2040.2070.1810.1950.167
Z 0.1880.1280.1270.1210.1260.112
Table 3. Attitude RMSEs of each scheme (degree).
Table 3. Attitude RMSEs of each scheme (degree).
AxisCKFHF-CKFMF-CKFNMF-CKFMHF-CKFNMHF-CKF
X 0.0430.0420.0420.0380.0400.036
Y 0.0510.0490.0520.0420.0480.040
Z 0.3240.2780.2750.2330.2430.215
Table 4. Position RMSEs of each scheme (meter).
Table 4. Position RMSEs of each scheme (meter).
AxisCKFHF-CKFMF-CKFNMF-CKFMHF-CKFNMHF-CKF
X 0.4820.2570.1240.1090.1140.092
Y 0.6720.4090.2240.1910.1980.178
Z 0.4880.2730.1400.1270.1300.121
Table 5. Attitude RMSEs of each scheme (degree).
Table 5. Attitude RMSEs of each scheme (degree).
AxisCKFHF-CKFMF-CKFNMF-CKFMHF-CKFNMHF-CKF
X 0.0460.0440.0430.0410.0420.038
Y 0.0520.0500.0560.0480.0500.046
Z 0.3260.3120.3010.2540.2650.239

Share and Cite

MDPI and ACS Style

Jiang, C.; Zhang, S.-B.; Zhang, Q.-Z. Adaptive Estimation of Multiple Fading Factors for GPS/INS Integrated Navigation Systems. Sensors 2017, 17, 1254. https://doi.org/10.3390/s17061254

AMA Style

Jiang C, Zhang S-B, Zhang Q-Z. Adaptive Estimation of Multiple Fading Factors for GPS/INS Integrated Navigation Systems. Sensors. 2017; 17(6):1254. https://doi.org/10.3390/s17061254

Chicago/Turabian Style

Jiang, Chen, Shu-Bi Zhang, and Qiu-Zhao Zhang. 2017. "Adaptive Estimation of Multiple Fading Factors for GPS/INS Integrated Navigation Systems" Sensors 17, no. 6: 1254. https://doi.org/10.3390/s17061254

APA Style

Jiang, C., Zhang, S. -B., & Zhang, Q. -Z. (2017). Adaptive Estimation of Multiple Fading Factors for GPS/INS Integrated Navigation Systems. Sensors, 17(6), 1254. https://doi.org/10.3390/s17061254

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