Next Article in Journal
A Method of Detections’ Fusion for GNSS Anti-Spoofing
Previous Article in Journal
The Longitudinal Force Measurement of CWR Tracks with Hetero-Cladding FBG Sensors: A Proof of Concept
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Adaptive H-Infinity Filtering Algorithm for the GPS/INS Integrated Navigation

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.
Submission received: 13 October 2016 / Revised: 6 December 2016 / Accepted: 9 December 2016 / Published: 19 December 2016
(This article belongs to the Section Remote Sensors)

Abstract

:
The Kalman filter is an optimal estimator with numerous applications in technology, especially in systems with Gaussian distributed noise. Moreover, the adaptive Kalman filtering algorithms, based on the Kalman filter, can control the influence of dynamic model errors. In contrast to the adaptive Kalman filtering algorithms, the H-infinity filter is able to address the interference of the stochastic model by minimization of the worst-case estimation error. In this paper, a novel adaptive H-infinity filtering algorithm, which integrates the adaptive Kalman filter and the H-infinity filter in order to perform a comprehensive filtering algorithm, is presented. In the proposed algorithm, a robust estimation method is employed to control the influence of outliers. In order to verify the proposed algorithm, experiments with real data of the Global Positioning System (GPS) and Inertial Navigation System (INS) integrated navigation, were conducted. The experimental results have shown that the proposed algorithm has multiple advantages compared to the other filtering algorithms.

1. Introduction

Accurate Global Positioning System (GPS) measurement data can be used to control the Inertial Navigation System (INS) accumulative errors. On the other hand, the INS can be used in cases when the GPS fails. Thus, the integration of GPS and INS has been widely adopted in the field of dynamic navigation and positioning. The Kalman filter is one of the most celebrated real-time optimal estimators [1], and therefore it has become the basic data-fusion approach in navigation and positioning [2,3]. Unfortunately, the Kalman filter performance depends on many factors, thus it cannot meet the requirements of certain nonlinear systems and the performance is susceptible to outliers. Hence, many improved filtering algorithms were proposed, for instance, the unscented Kalman filter [4,5], the particle filter [6], the robust filter [7,8], the adaptive filter [9,10,11,12], and so forth. The nonlinear problem can be handled by the unscented Kalman filter or the particle filter, however, they both become unstable in a high-dimensional case. On the other hand, the cubature Kalman filter is a recently developed nonlinear filtering algorithm [13]. Compared with the unscented Kalman filter, the cubature Kalman filter shows better performance in stability; therefore, it has been adopted in the GPS/inertial measurement unit (IMU) integrated navigation system [14]. Many forms of robust filtering algorithms were proposed to control the influence of outliers, but the influence of the model errors was neglected. The multiple-model-based adaptive estimation (MMAE) and innovation-based adaptive estimation (IAE) represent the adaptive Kalman filtering approaches. In the MMAE, there is a set of Kalman filters which work in parallel under different models, while in the IAE the adaptation is applied directly to the statistical information. Moreover, the IAE is more applicable to the GPS/INS integrated navigation systems [9]. The DIA (detection, identification, and adaptation) methods represent the recursive testing procedures, and they were proposed to eliminate the effects of outliers [15]. In the DIA methods, firstly, the model errors are detected and identified, then deviations of the state estimates caused by the model errors are eliminated, and finally, the adaptation in which the model is recovered from the identified errors, is implemented. However, the identification is quite difficult, especially when the measurements are not accurate [16].
Recently, the adaptively robust filtering has been proposed for dynamic navigation and positioning [16,17]. The combination of adaptive filter and robust filter provides the required adaptivity and robustness, and the adaptive–robust filter has a better ability to mitigate the influences of dynamic model deviations and outlying measurements [16,18,19]. Nevertheless, the number of the measurements should be equal to or greater than the state vector dimension. Hence, the adaptively robust filtering can be applied in a limited number of cases [20]. It should be noticed that the aforementioned filtering algorithms do not consider the interference. The minimization of estimation error was used to develop the H-infinity filter without the information on statistical properties [21,22]. The H-infinity filter can be used to address the uncertainties in the measurement noise [20,23]. Hassibi has shown that the least mean squares (LMS)-based adaptive filtering algorithm represents the H-infinity optimal, and that the H-infinity algorithm is applicable in the uncertain environments [24]. However, the median-based filters might be highly robust, but not efficient [20]. Moreover, the H-infinity filter fails in the presence of outliers [8]. The performance of the H-infinity filter can be improved only if the effects of the outliers are controlled. The robust estimation method provides a way to control the influence of outliers [25,26,27]. Hence, the robust estimation method can be applied to improve the robustness and stability of the H-infinity filter. Accordingly, in order to control the outliers influence, a new filtering scheme based on combination of the adaptive filter and the H-infinity filter should be developed.
This paper focuses on a comprehensive filtering algorithm based on the combination of the adaptive filter and the H-infinity filter. According to the model errors and interference, a novel filtering algorithm intended for the GPS/INS integrated navigation is developed. The GPS/INS integrated dynamic navigation was tested using the cubature Kalman filter. The experimental results have shown that the proposed algorithm has better performance compared to the commonly used filtering algorithms.
The paper is organized as follows. In Section 2, a basic principle of the adaptive filter and the H-infinity filtering algorithms are introduced, and scaling factors are discussed. In Section 3, the adaptive H-infinity filtering algorithm and the robust estimation method are presented. The equations of dynamic model and measurement model intended for the GPS/INS integrated navigation are listed in Section 4. In Section 5, the proposed algorithm is verified by experiments and the experimental results are provided; in addition, a comparison to other filtering algorithms is performed and the obtained results are discussed. Lastly, a brief conclusion is given in Section 6.

2. The Adaptive Kalman Filter and the H-Infinity Filter

2.1. The Adaptive Kalman Filter

Various filtering algorithms have been developed in order to mitigate the influence of model errors, such as the adaptive Kalman filter. The adaptive algorithms intended for the GPS/INS integrated navigation can be divided into multiple-model adaptive estimation, adaptive stochastic modeling, and covariance scaling [28]. In the multiple-model adaptive estimation, a set of Kalman filters works in parallel, under different models of the filter’s statistical information [9]. A single or multiple factor(s) of the state covariance matrix can be used to improve the stability and convergence performance, thus, it is important to select a suitable scaling factor. In general, the scaling factor is determined empirically, based on experience or innovation and residual information [23]. In reference [12], two adaptive factors were derived, and the filtering performance of the algorithm based on the optimal adaptive factor was superior to the one based on the nonoptimal adaptive factors.
The dynamic model and measurement model are given by:
{ x k = Φ k , k 1 x k 1 + w k z k = H k x k + v k ,
where x k is the state vector, Φ k , k 1 presents the state transition matrix, H k is the measurement matrix, z k denotes the measurement vector, and w k and v k are the system noise and measurement noise, respectively. The predicted state vector and its covariance matrix are defined by:
x k / k 1 = Φ k , k 1 x k 1 ,
P k / k 1 = Φ k , k 1 P k 1 Φ T k , k 1 + Q k ,
where x k / k 1 is the predicted state vector, and P k / k 1 and Q k denote the predicted covariance matrices of x k / k 1 and w k , respectively.
According to the rule of the least squares estimation, the following can be written:
V k T R k 1 V k + V x k / k 1 T P k / k 1 1 V x k / k 1 = min ,
where V k and V x k / k 1 denote the residual vectors of z k and x k / k 1 , respectively.
The solution of the standard Kalman filter is obtained by:
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 ( z k H k x k / k 1 ) ,
where R k is the covariance matrix of the measurement noise. If the risk function is defined by:
V k T R k 1 V k + α k V x k / k 1 T P k / k 1 1 V x k / k 1 = min ,
then the solution of the adaptive Kalman filter is obtained by:
x k / k = ( H k T R k 1 H k + α k P k / k 1 1 ) 1 ( H k T R k 1 z k + α k P k / k 1 1 x k / k 1 ) ,
where α k denotes the adaptive factor. The equivalent expression of Equation (8) is in reference [27]:
x k / k = x k / k 1 + 1 α k P k / k 1 H k T ( 1 α k H k P k / k 1 H k T + R k ) ( z k H k x k / k 1 ) ,
Four types of the adaptive factors and error-learning statistics were summarized in reference [17]. The predicted residual statistic was adopted if no redundant information exists, and the adaptive factor based on the three-segment function was [16,29]:
α k = { 1    | Δ V ¯ k | c 0 c 0 | Δ V ¯ k | ( c 1 | Δ V ¯ k | c 1 c 0 ) 2 c 0 < | Δ V ¯ k | c 1 0    | Δ V ¯ k | > c 1 ,
where Δ V ¯ k denotes the learning statistic based on the predicted residual, Δ V ¯ k = ( V ¯ k T V ¯ k tr ( P V ¯ k ) ) 1 2 , P V ¯ k = H k P k / k 1 H k T + R k , V ¯ k denotes the innovation vector, further, 1.0 c 0 1.5 , and 3.0 c 1 8.5 . In Equation (9), α k 0 , therefore, when Equation (9) is used, the adaptive factor based on the two-segment function [17], defined by Equation (11), should be adopted:
α k = { 1 | Δ V ¯ k | c c | Δ V ¯ k | | Δ V ¯ k | > c ,
where 1.0 c 2.5 , and the optimal value of c is 1.0 [12,17].

2.2. Basic Principle of the H-Infinity Filter

The H-infinity filter is a special form of the Kalman filter, and the principle of the H-infinity filter is based on the H-infinity optimal estimation that guarantees the smallest estimation energy error for all possible disturbances of the fixed energy [30]. The cost function of a nonlinear filter is defined by:
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 is the number of filtering epochs, x 0 is the initial value of the state vector x with the covariance matrix P 0 , x ^ 0 and x ^ k are 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 ) . In general, the estimate, x ^ k , which satisfies x ^ k = arg   min J , should be found in order to minimizes J . Actually, an analytical solution for the optimal H-infinity filter is difficult to achieve. Thus, a suboptimal recursion algorithm is developed by setting a threshold value, γ , which satisfies the following Riccati inequality [31]:
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 , L k is the coefficient matrix of the constraint equation, and L k is generally defined by the unit matrix I . The recursion equations of the H-infinity filter are listed below [31]:
x k / k 1 = Φ k , k 1 x k 1 ,
P k / k 1 = Φ k , k 1 P k 1 / k 1 Φ k , k 1 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 ] ,
The H-infinity filter minimizes the estimation error in the worst case, which makes it more robust than the standard Kalman filter. In addition, the H-infinity filter becomes more robust when the constraint parameter γ decreases. However, the value of γ should not be in the vicinity of zero, because that might cause the divergence of the H-infinity filter [32]. In general, γ is fixed to certain value, which is chosen by experience. In addition, the H-infinity filter represents a rigorous method for the system with an unreliable model [30].

3. A Novel Adaptive H-Infinity Filtering Algorithm

In the adaptive filtering algorithms, it is assumed that the measurements are reliable and that the predicted residual can reflect the deviations of the model information. However, some of the adaptive filtering algorithms cannot control the effects of the outliers. As mentioned before, the H-infinity filter can overcome the uncertainties in the measurements, but the effects of the outliers cannot be controlled just by the H-infinity filter, thus, some robust Kalman filtering algorithms were developed [7].
In this paper, an integrated adaptive H-infinity filtering algorithm is proposed. An adaptive factor was constructed in order to control the influence of the dynamic model errors, and the H-infinity filter was adopted to address the uncertain interference. The recursion equations of the adaptive H-infinity filter can be expressed by:
x k / k 1 = Φ k , k 1 x k 1 ,
P k / k 1 = Φ k , k 1 P k 1 / k 1 Φ k , k 1 T + Q k ,
x k / k = x k / k 1 + K ¯ k ( z k H k x k / k 1 ) ,
K ¯ k = 1 α k P k / k 1 H k T ( 1 α k H k P k / k 1 H k T + R k ) 1 ,
P k / k = 1 α k P k / k 1 1 α k 2 Φ 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 ] + 1 α k [ H k L k ] P k / k 1 [ H k T L k T ] ,
where K ¯ k denotes the equivalent gain matrix of the H-infinity filter. Since, there were no redundant measurements, the adaptive factor based on the two-segment function was adopted, and the predicted residual was selected as the statistic value.
In the integrated adaptive H-infinity filtering algorithm, the robust estimation method was employed to decrease the effects of outlying measurements. The selection of the equivalent weight and equivalent covariance matrix should be considered in the robust estimation. In general, the components of the predicted state vector are correlated, so the equivalent covariance matrix constructed by scaling factor, λ i j , is applied. The scaling factor is defined by reference [33]:
λ i j = λ i i λ j j ,
λ i i = { 1 | V ¯ X ¯ k i | c | V ¯ X ¯ k i | c | V ¯ X ¯ k i | > c ,
where | V ¯ X ¯ k i | denotes the component of the predicted residual vector, 1.0 c 1.5 , and λ j j can be expressed the same as λ i i . Thus, the equivalent covariance matrix R ¯ k is obtained by
R ¯ k = λ i j R k ,
The entire process of the proposed filtering algorithm is shown in Figure 1.

4. The GPS/INS Integrated Navigation

Recently, in the GPS/INS integrated navigation, three types of integration have been developed: loosely coupled, tightly coupled, and deeply coupled [34]. In the loosely coupled integrated navigation, a 15-dimension state vector is adopted, which includes the deviations of position, velocities, attitudes, and the noises of gyroscope and accelerometer. Thus, the state vector x ^ is defined 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 ] ,
Since, the GPS/INS integrated navigation system represents a nonlinear system, the cubature Kalman filter is applied to address the nonlinear problem. For a discrete nonlinear system, it can be written as follows:
{ 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, and [ 1 ] denotes the point ( 1 0 ) in this paper. Then, the equations of the cubature Kalman filter are given by reference [14]:
(a) 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 ,
(b) 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 Y k / k 1 = h ( X i , k / k 1 ) y k / k 1 = 1 m i = 1 m Y i , k / k 1 P z z , k / k 1 = 1 m i = 1 m Y i , k / k 1 Y i , k / k 1 T y k / k 1 y k / k 1 T + R k P x z , k / k 1 = 1 m i = 1 m X i , k / k 1 Y i , k / k 1 T x k / k 1 y k / k 1 T ,
then the final measurement update equations are obtained by:
{ 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 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; and Z k / k 1 denotes the propagated cubature points of the measurement vector. Accordingly, the adaptive H-infinity filtering algorithm intended for the GPS/INS integrated navigation system is summarized as follows:
(a) Time update
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 ,
(b) 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 = 1 α k P k / k 1 [ P x z , k 1 α k P k / k 1 ] [ P z z , k R k + I P x z , k T P x z , k 1 α k P k / k 1 γ 2 I ] 1 [ P x z , k T ( 1 α k P k / k 1 ) T ] ,
In addition, R k should be replaced with the equivalent covariance matrix R ¯ k in order to make the algorithm more robust.
In the loosely coupled integrated navigation, the differences in position and velocity, Z ρ , between GPS and INS are considered as the external measurements, namely:
Z ρ ( t ) = ρ GPS ρ INS ,
where ρ GPS , namely, x GPS , y GPS , z GPS , v XGPS , v YGPS , and v ZGPS , presents the output information of the GPS in the WGS-84 coordinate system, while the output information of the INS is ρ INS , namely, x INS , y INS , z INS , v XINS , v YINS , and v ZINS . Then, the measurement equation is as follows:
Z k = [ r GPS r INS v GPS v INS ] ,
where Z k represents the measurement vector of the integrated system, and r GPS , r INS and v GPS , v INS are position and velocity output information of GPS and INS, respectively. Compared to the loosely coupled, the tightly coupled integration has better performance in terms of precision. However, the loosely coupled navigation can be implemented easily, and the computation process is more concise [35].

5. Test Cases and Data Analysis

5.1. Test Case 1

In this case, the GPS and INS data were obtained by a small aerial vehicle equipped with inertial sensors of automotive-grade quality and a GPS receiver [36]. Four schemes were implemented to examine the performance of the proposed filtering algorithm:
Scheme 1: Kalman Filter (KF);
Scheme 2: Adaptive Kalman Filter (AKF) (two-segment function was adopted and c = 1.0 );
Scheme 3: H-infinity Filter (HF);
Scheme 4: Adaptive H-infinity Filter (AHF) (two-segment function was adopted and c = 1.0 , and the robust estimation method was adopted with the double-factors and c = 1.5 , where c denotes the criterion of the double-factors equivalent covariance matrix);
Position errors of all mentioned schemes are demonstrated in the following:
As it can be seen in previous figures, there is little difference between Figure 2 and Figure 3, which indicates that there are limited effects of the model errors. Due to the minimization of the worst-case estimation error, the H-infinity filter performs slightly better than the adaptive Kalman filter, which can be seen from Figure 4 and Table 1. Nonetheless, Figure 5 indicates that a noticeable improvement was obtained by the adaptive H-infinity filtering algorithm in terms of dynamic model errors and the uncertain interferences.
The Root Mean Squares Errors (RMSEs) of the schemes are presented in Table 1.

5.2. Test Case 2

In this experiment, data were collected by a vehicle equipped with a GPS/INS integrated navigation system, which was composed of two GPS receivers and an inertial measurement unit (IMU). One of the GPS receivers was set as a reference station, and another receiver as well as the IMU were mounted on the vehicle. GPS data were calculated by the double difference pseudorange with variances of 0.25 m2 and 0.0025 m2/s2, respectively. The sampling frequencies of GPS and IMU were 1 Hz and 100 Hz, respectively. The precise results achieved by the double difference carrier phase were regarded as references. Afterwards, four schemes were performed in the integrated navigation system:
Scheme 1: Kalman Filter (KF);
Scheme 2: Adaptive Kalman Filter (AKF) (two-segment function was adopted and c = 1.0 );
Scheme 3: H-infinity Filter (HF);
Scheme 4: Adaptive H-infinity Filter (AHF) (two-segment function was adopted and c = 1.0 , and the robust estimation method was adopted with the double-factors where c = 1.5 ).

5.2.1. Experiments without Outliers

Position errors of used schemes are demonstrated in the following:
Due to the vehicle’s movement over the bumps, disturbances can be found in both Figure 6 and Figure 7, which indicates that the robustness of the KF and AKF algorithms should be improved. Because of decrease of the dynamic model errors, the AKF algorithm shows a slightly better performance than the KF algorithm. Apparently, the uncertain interference is controlled, thus, the HF and AHF algorithms (Figure 8 and Figure 9) perform much better than the first two algorithms. The robustness of the AHF algorithm is improved by the robust estimation method, therefore, the error amplitude is reduced. The RMSE of each scheme is given in Table 2.
As mentioned before, RMSEs of the AKF and HF algorithms, presented in Table 2, are both smaller than RMSE of the KF algorithm. The integration of the AKF and HF algorithms into the AHF algorithm provides better performance compared to the other algorithms.

5.2.2. Experiments with Outliers

In these experiments, we used data from the Section 5.2.1, but the single gross errors were added artificially to the GPS measurement data at the 160th, 360th, 460th, and 560th epochs, respectively, and the constantly changed outliers were added to all epochs between the 251th and 280th epochs. Then, the previously used schemes were implemented again, and the obtained position errors of all schemes are displayed in Figure 10, Figure 11, Figure 12 and Figure 13.
According to the presented results, it can be concluded that in the presence of outliers, filtering results were mostly dependent on outliers. Moreover, Figure 10 and Figure 11 show that the KF and AKF filtering results have similar trends, while the latter figure manifests a reduced error. As it can be seen in Figure 12, the filtering results of the HF algorithm are stable except in the epochs that contain outliers. Furthermore, Figure 13 shows that the proposed filtering algorithm has better performance than other algorithms, so the influence of single and constantly changed outliers is reduced using the AHF algorithm. The RMSE values of all used scheme are presented in Table 3.
Compared to the Kalman filter, the filtering precision of the AKF and HF algorithms is improved. However, both AKF and HF algorithms are significantly affected by the outliers. On the other hand, the AHF algorithm shows better fault tolerance and robustness compared to other schemes. The adaptive H-infinity filter is more robust because of the robust estimation method, based on the control of dynamic model errors and uncertain interference. In all presented cases, RMSEs of the AHF algorithm are the smallest for all coordinates, which means that the positions calculated by the AHF algorithm are in good agreement with the actual positions.

6. Conclusions

An integrated adaptive H-infinity filtering algorithm for mitigation of positioning errors caused by dynamic model errors and uncertain interference is presented. The proposed filtering algorithm is improved by a robust estimation method, wherein both single and constantly changed outliers are considered. In addition, the proposed algorithm was verified by the GPS/INS integrated navigation. The results have shown that the proposed algorithm has better performance than other algorithms.
The conclusions can be summarized as follows:
(1)
The adaptive Kalman filtering algorithms were developed in order to reduce the positioning errors, and the proper adaptive factors were selected. The H-infinity filtering algorithm performed well in the GPS/INS integrated navigation system that contained the uncertainties. However, the performance was greatly affected by the outliers.
(2)
The integration of the adaptive Kalman filter, the H-infinity filter, and the robust estimation method provided the AHF algorithm, which can address the deviations caused by dynamic model errors and interference. Since the proposed algorithm was verified by real data, a wide application of the proposed AHF algorithm in dynamic navigation and positioning can be expected.

Acknowledgments

The authors are grateful for the support of the National Natural Science Foundation of China (NO.41504032), the Natural Science Foundation of Jiangsu Province (NO.BK20150175), the Specialized Research Fund for the Doctoral Program of Higher Education of China (NO.20130095110022) and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Author Contributions

Chen Jiang originated this work and prepared the manuscript; Shu-bi Zhang contributed to the theory studies; Qiu-zhao Zhang designed the experiments and improved the quality of this work. All authors reviewed the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kalman, R.E. A new approach to linear filtering and prediction problems. J. Basic Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef]
  2. Sang, S.B.; Zhai, R.Y.; Zhang, W.D.; Sun, Q.R.; Zhou, Z.Y. A self-developed indoor three-dimensional pedestrian localization platform based on MEMS sensors. Sens. Rev. 2015, 35, 157–167. [Google Scholar] [CrossRef]
  3. Neto, P.; Mendes, N.; Moreira, A.P. Kalman filter-based yaw angle estimation by fusing inertial and magnetic sensing: A case study using low cost sensors. Sens. Rev. 2015, 35, 244–250. [Google Scholar] [CrossRef]
  4. Wan, E.A.; van der Merwe, R. The Unscented Kalman Filter for Nonlinear Estimation. In Proceedings of the IEEE 2000, Adaptive Systems for Signal Processing, Communications and Control Symposium (AS-SPCC), Lake Louise, AB, Canada, 1–4 October 2000; pp. 153–158.
  5. Abd Rabbou, M.; El-Rabbany, A. Integration of GPS precise point positioning and MEMS-based INS using unscented particle filter. Sensors 2015, 15, 7228–7245. [Google Scholar] [CrossRef] [PubMed]
  6. Gustafsson, F. Particle filter theory and practice with positioning applications. IEEE Aerosp. Electron. Syst. Mag. 2010, 25, 53–82. [Google Scholar] [CrossRef]
  7. Koch, K.R.; Yang, Y.X. Robust Kalman filter for rank deficient observation model. J. Geod. 1998, 72, 436–441. [Google Scholar] [CrossRef]
  8. Gandhi, M.A.; Mili, L. Robust Kalman filter based on a generalized maximum-likelihood-type estimator. IEEE Trans. Signal Process. 2010, 58, 2509–2520. [Google Scholar] [CrossRef]
  9. Mohamed, A.H.; Schwarz, K.P. Adaptive Kalman filtering for INS/GPS. J. Geod. 1999, 73, 193–203. [Google Scholar] [CrossRef]
  10. Ding, W.D.; Wang, J.L.; Rizos, C.; Kinlyside, D. Improving adaptive Kalman estimation in GPS/INS integration. J. Navig. 2007, 60, 517–529. [Google Scholar] [CrossRef]
  11. Yang, Y.X.; Xu, T.H. An adaptive Kalman filter based on Sage windowing weights and variance components. J. Navig. 2003, 56, 231–240. [Google Scholar] [CrossRef]
  12. Yang, Y.X.; Gao, W.G. An optimal adaptive Kalman filter. J. Geod. 2006, 80, 177–183. [Google Scholar] [CrossRef]
  13. Arasaratnam, I.; Haykin, S. Cubature kalman filters. IEEE Trans. Autom. Control 2009, 54, 1254–1269. [Google Scholar] [CrossRef]
  14. Zhao, Y.W. Performance evaluation of cubature Kalman filter in a GPS/IMU tightly-coupled navigation system. Signal Process. 2016, 119, 67–79. [Google Scholar] [CrossRef]
  15. Teunissen, P.J.G. Quality Control in Integrated Navigation Systems. In Proceedings of the IEEE Symposium on Position Location and Navigation. A Decade of Excellence in the Navigation Sciences (IEEE PLANS’90), Las Vegas, NV, USA, 20–23 March 1990; pp. 158–165.
  16. 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]
  17. Yang, Y.X.; Ren, X.; Xu, Y. Main progress of adaptively robust filter with applications in navigation. J. Navig. Position. 2013, 1, 9–15. [Google Scholar]
  18. Yang, Y.X.; Gao, W.G. Comparison of two fading filters and adaptively robust filter. Geosp. Inf. Sci. 2007, 10, 200–203. [Google Scholar] [CrossRef]
  19. Yang, Y.X.; Cui, X.Q. Adaptively robust filter with multi adaptive factors. Surv. Rev. 2008, 40, 260–270. [Google Scholar] [CrossRef]
  20. Chang, G.B. Robust Kalman filtering based on Mahalanobis distance as outlier judging criterion. J. Geod. 2014, 88, 391–401. [Google Scholar] [CrossRef]
  21. 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]
  22. 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]
  23. 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]
  24. Hassibi, B.; Kailath, T. H∞ Adaptive Filtering. In Proceedings of the 1995 International Conference on IEEE Acoustics, Speech, and Signal Processing (ICASSP-95), Detroit, MI, USA, 9–12 May 1995; Volume 2, pp. 949–952.
  25. Zhou, J.W. Classical theory of errors and robust estimation. Acta Geod. Cartogr. Sin. 1989, 18, 115–120. [Google Scholar]
  26. Yang, Y.X. Robust estimation of geodetic datum transformation. J. Geod. 1999, 73, 268–274. [Google Scholar] [CrossRef]
  27. Yang, Y.X.; Xu, T.H. GNSS receiver autonomous integrity monitoring (RAIM) algorithm based on robust estimation. Geod. Geodyn. 2016, 7, 117–123. [Google Scholar] [CrossRef]
  28. Almagbile, A.; Wang, J.L.; Ding, W.D. Evaluating the performances of adaptive Kalman filter methods in GPS/INS integration. J. Glob. Position. Syst. 2010, 9, 33–40. [Google Scholar] [CrossRef]
  29. Yang, Y.X. Adaptive Navigation and Dynamic Positioning; Surveying and Mapping Press: Beijing, China, 2006. [Google Scholar]
  30. Simon, D. Optimal State Estimation: Kalman, H and Nonlinear Approaches; John Wiley and Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  31. Chen, Y.R.; Yuan, J.P. An improved robust H multiple fading fault-tolerant algorithm for INS/GPS integrated navigation. J. Astronaut. 2009, 30, 930–936. [Google Scholar]
  32. Hassibi, C.; Sayed, A.H.; Kailath, T. Linear estimation in Krein spaces-part II: Applications. IEEE Trans. Automat. Control 1996, 41, 34–49. [Google Scholar] [CrossRef]
  33. Yang, Y.X.; Song, L.J.; Xu, T.H. Robust parameter estimation for geodetic correlated observations. Acta Geod. Cartogr. Sin. 2002, 31, 95–99. [Google Scholar]
  34. Zhang, Q.Z.; Stephenson, S.; Meng, X.L.; Zhang, S.B.; Wang, Y.J. A new robust filtering for a GPS/INS loosely coupled integration system. Surv. Rev. 2016, 48, 181–187. [Google Scholar] [CrossRef]
  35. 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]
  36. Gleason, S.; Gebre-Egziabher, D. GNSS Applications and Methods; Artech House: London, UK, 2009. [Google Scholar]
Figure 1. A flow chart of the adaptive H-infinity algorithm.
Figure 1. A flow chart of the adaptive H-infinity algorithm.
Sensors 16 02127 g001
Figure 2. Position errors of the Kalman filter (KF) algorithm.
Figure 2. Position errors of the Kalman filter (KF) algorithm.
Sensors 16 02127 g002
Figure 3. Position errors of the adaptive Kalman filter (AKF) algorithm.
Figure 3. Position errors of the adaptive Kalman filter (AKF) algorithm.
Sensors 16 02127 g003
Figure 4. Position errors of the H-infinity filter (HF) algorithm.
Figure 4. Position errors of the H-infinity filter (HF) algorithm.
Sensors 16 02127 g004
Figure 5. Position errors of the adaptive H-infinity filter (AHF) algorithm.
Figure 5. Position errors of the adaptive H-infinity filter (AHF) algorithm.
Sensors 16 02127 g005
Figure 6. Position errors of the KF algorithm.
Figure 6. Position errors of the KF algorithm.
Sensors 16 02127 g006
Figure 7. Position errors of the AKF algorithm.
Figure 7. Position errors of the AKF algorithm.
Sensors 16 02127 g007
Figure 8. Position errors of the HF algorithm.
Figure 8. Position errors of the HF algorithm.
Sensors 16 02127 g008
Figure 9. Position errors of the AHF algorithm.
Figure 9. Position errors of the AHF algorithm.
Sensors 16 02127 g009
Figure 10. Position errors of the KF algorithm.
Figure 10. Position errors of the KF algorithm.
Sensors 16 02127 g010
Figure 11. Position errors of the AKF algorithm.
Figure 11. Position errors of the AKF algorithm.
Sensors 16 02127 g011
Figure 12. Position errors of the HF algorithm.
Figure 12. Position errors of the HF algorithm.
Sensors 16 02127 g012
Figure 13. Position errors of the AHF algorithm.
Figure 13. Position errors of the AHF algorithm.
Sensors 16 02127 g013
Table 1. Root mean square error (RMSE) values of the schemes (m).
Table 1. Root mean square error (RMSE) values of the schemes (m).
AxisKFAKFHFAHF
X 0.0460.0430.0390.035
Y 0.0470.0450.0380.034
Z 0.0580.0490.0430.042
Table 2. RMSE values of the schemes (m).
Table 2. RMSE values of the schemes (m).
AxisKFAKFHFAHF
X 0.1290.1050.0960.078
Y 0.2840.2510.2040.121
Z 0.1880.1600.1280.094
Table 3. RMSE values of the schemes (m).
Table 3. RMSE values of the schemes (m).
AxisKFAKFHFAHF
X 0.4460.2820.2320.083
Y 0.4940.3380.2980.128
Z 0.4130.2660.2220.099

Share and Cite

MDPI and ACS Style

Jiang, C.; Zhang, S.-B.; Zhang, Q.-Z. A New Adaptive H-Infinity Filtering Algorithm for the GPS/INS Integrated Navigation. Sensors 2016, 16, 2127. https://doi.org/10.3390/s16122127

AMA Style

Jiang C, Zhang S-B, Zhang Q-Z. A New Adaptive H-Infinity Filtering Algorithm for the GPS/INS Integrated Navigation. Sensors. 2016; 16(12):2127. https://doi.org/10.3390/s16122127

Chicago/Turabian Style

Jiang, Chen, Shu-Bi Zhang, and Qiu-Zhao Zhang. 2016. "A New Adaptive H-Infinity Filtering Algorithm for the GPS/INS Integrated Navigation" Sensors 16, no. 12: 2127. https://doi.org/10.3390/s16122127

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