Abstract

In the literature, the fading factor was constructed to overcome the shortage of model uncertainties in the Kalman filter. However, the a priori covariance matrix might be inflated abnormally by the fading factor once the measurement is unreliable. Thus, the fading factor may become invalid, and this problem is rarely discussed and tested. In this paper, squares of the Mahalanobis distance are introduced as the judging index, and the fading factor or the covariance inflation factor is adopted conditionally according to the hypothesis testing result. Therefore, an adaptive filtering scheme based on the Mahalanobis distance is put forward for the systems with model uncertainties. The proposed algorithm is implemented with the actual data collected by the integration of the global navigation satellite system (GNSS) and the inertial navigation system and INS (inertial navigation system) integrated systems (INS). For the systems with model uncertainties, experimental results demonstrate that the influences of both the outlying measurements and model errors are controlled effectively with the proposed scheme.

1. Introduction

As the linear estimator of the mean and covariance, the Kalman filter has become the classic fusion algorithm in many fields [1, 2], and it is implemented as the basic fusion method in the data processing with multiple sensors [3, 4]. It has been proved that the Kalman filter is optimal only when the assumptions of Gaussian-distributed process or measurement noise hold [5]. Unfortunately, the Kalman filter is susceptible to outlying measurements, and it performs inadequately with the model errors and uncertain statistical information. Aiming at the outlying measurements, many robust Kalman filtering algorithms were proposed. Outlier detection for all measurements is the most straightforward strategy, and the measurements with relatively large residuals will be rejected [6]. Nevertheless, all measurements should be tested, and this strategy may be inefficient and complex. The filters based on the median may be highly robust; however, many measurements are ignored, and the low efficiency limits its practical applications [7]. The H∞ filter was put forward aiming at the uncertain noises of measurement with the minimized estimation error for all possible disturbances [5]. However, the performance would be degraded significantly by the outliers [8]. As the generalized Kalman filter, the Bayesian estimator was derived robustly based on the M-estimation [9]. DIA (detection, identification, and adaptation) methods were developed to resist the influences of outliers, but the identification is valid only when the measurement is reliable [10]. Influences of the model errors can be weakened with many types of adaptive filters [1113]. In terms of adaptive filter, the MMAE (multiple-model-based adaptive estimation) and IAE (innovation-based adaptive estimation) [14] are two basic strategies. For the data fusion of GNSS/INS integrated systems, the IAE strategy performs better than the MMAE strategy [11]. Besides, another type of adaptive-robust filter where both the adaptation and robustness were considered simultaneously was developed using the adaptive factors and the M-estimation, including four adaptive factors and four detective statistics [15]. Based on this adaptive-robust filter, the influences on the adaptive factor from the outlying measurements were discussed [16], and an alternative strategy was constructed to control the influences of both model uncertainties and outliers. However, the covariance of the measurement noises remains constant at the normal epochs, and the relatively smaller outliers are more likely to be neglected which should be improved further.

In theory, the fading filter is another kind of adaptive filter, and it can control the influences of the dynamic model errors[17]. The most important problem for the fading filter is to construct a suitable factor . Then, the fading matrix was proposed to adjust the covariance matrix in different data channels [18, 19], and the fading matrix has been adopted in the data processing of GNSS/INS integrated systems [20]. The main superiority to the conventional Kalman filter lies on inflating the covariance matrix of the a priori state, namely, , and this indicates that the state estimation relies more on the current observation information. However, both the single factor and the fading matrix are constructed based on the residual vector. Therefore, the fading filter performs well only when the measurement is reliable, and the filter divergence may happen in the presence of outlying measurements, but this problem is seldom discussed. Consequently, aiming at the influences of outlying measurements, the robustness of the conventional fading filter should be improved further.

In this paper, a modified adaptive filtering scheme using the fading factors is proposed for the system with model uncertainty, and the Mahalanobis distance from the measurement to its prediction is applied to construct the judging index for the hypothesis test. Whether the fading filter or the robust estimation method is performed at each epoch depends on the result of hypothesis test. The proposed algorithm is tested using the actual data obtained through the self-developed GNSS/INS integrated systems in actual environment. Both the conventional fading filter and the proposed algorithm are implemented in the testing section of this paper. Results demonstrate that the proposed scheme is superior among all the other tested algorithms of this paper with or without the outlying measurements, and the filter divergence is restrained with the proposed algorithm.

Remaining of this paper is arranged as follows. The related theory of the fading filter is introduced in Section 2. In Section 3, a modified adaptive data fusion scheme is constructed based on the hypothesis test and the conventional fading filter. In Section 4, the dynamic model and measurement model for the GNSS/INS integrated systems are provided; then, a flowchart of the proposed algorithm is demonstrated. In Section 5, experiments with actual data are implemented, and performance the comparative algorithms together with the proposed scheme is tested. Section 6 provides the conclusions of this paper.

2. The Fading Filter

2.1. Basic Models of the Fading Filter

Kalman filter performs well when the assumptions of Gaussian distribution hold. However, the filter divergence may happen, resulting from big model errors. Aiming at this problem, the fading filter was proposed to limit the “memory length” of the Kalman filter [21]. Assume that and are the state vector at epoch and , respectively, denotes the state transition matrix, and is the state noise matrix; then, the dynamic equation is obtained and the a priori state is given by

Assume that and are the measurement matrix and measurement vector, respectively, is the measurement noise, and is the fading factor (in general, ). Then, the optimality criterion is defined by [21]:

Then, the iterative solution is listed below: where denotes the equivalent gain matrix, denotes the equivalent covariance matrix of , denotes the covariance matrix of , and and are the covariance matrices of measurement and state noises, respectively.

In the fading filter, the covariance matrix of is inflated for times when compared with that of the conventional Kalman filter, and this indicates that the current measurement information is treated with bigger weight [21]. Consequently, the model errors brought from the previous state are weakened with the fading factor.

2.2. Construction of the Fading Factor or Matrix

In the fading filter, the main work concentrates on the construction of the reasonable factor. In theory, the factor should be inflated if errors of is abnormal. A fading filter with the optimal fading factor was developed, and the simplified and applicable fading factor is given by [22] and , , , and , where means the trace of a matrix, is the residual vector, is the covariance matrix of , and .

Apart from the fading filter with a single factor, the fading matrix for the fading filter was developed [18]. Different from the single factor, the fading filter with the fading matrix is adjusted in multiple data channels according to the observabilities of the state vector elements. The fading matrix is derived by where is the dimension of measurement equation, is the th element of , is the th observable element of , is the th diagonal element of and , is the th diagonal element of , , and is the threshold determined by the Chi-square distribution under the predetermined significance level. In Equation (9), are estimated adaptively, and the other elements in should be 1.

3. Modified Adaptive Data Fusion Scheme

As depicted in Section 2, the current measurement information is treated with bigger weight in the fading filter, and the current measurement information must be reliable to achieve an ideal filtering performance. In the fading filter however, the circumstance, when the current measurement information is unreliable, is rarely discussed and tested with actual data. In fact, both a single fading factor and the fading matrix in the fading filter are constructed by the predicted residual vector obtained from the measurements. Therefore, the performance of the fading factor would be affected easily by unreliable measurements, and a mistaken adjustment towards the a priori covariance matrix may result in inferior performance or even the filter divergence. Thus, the fading factor or matrix should be implemented in a timely and more reasonable way, and the strategy on processing the unreliable measurements should be discussed further.

With the Gaussian distribution process and measurement noise, the discrete linear stochastic state space model is given by [5]

An optimal solution is obtained if the assumptions of Gaussian distribution process and measurement noise hold with the mean and covariance . Probability density function of the dimensional measurement is [1]

Nevertheless, would no longer hold once the outlying measurements exist or the measurement noises disobey the Gaussian distribution. Therefore, if does not hold, and this in turn indicates that an outlier exists at some epoch or the Gaussian distribution of measurement noise is contaminated. Accordingly, the hypothesis test is performed with the null hypothesis that the measurement corresponds to the assumptions. In fact, the square of the Mahalanobis distance from the measurement to its mean , namely, , is applied as the test statistic [23], and this test statistic has been adopted as to construct the robust Kalman filter [1]. According to the inertial assumptions, the test statistic should obey the Chi-square distinction with the degree of freedom . Then, the -quantile is obtained with the provided significance level . Obviously, the null hypothesis should be accepted if is less than ; under this circumstance, no outliers exist, and the key problem is to weaken the influences of the model errors. Otherwise, the null hypothesis should be rejected, and the outliers should be processed carefully.

In order to weaken the influences of unreliable measurements, is estimated and updated in an adaptive way. In general, components in the vector of the a priori state are correlated; the covariance matrix of the measurement noise is updated with the covariance inflation factor , namely, where denotes the standardized predicted residual and is fixed within 1.0 ~ 1.5. Consequently, is rewritten as

In practice, the complex observation environment determines that is usually not a constant matrix, and the filter may perform better with adaptively updated . Therefore, at the epochs with normal measurements, is updated adaptively in a sequential way based on the Sage-Husa filter [24]; in this paper and the influences of the outlying measurements will be weakened further, namely, where and is the forgetting factor.

With the hypothesis test based on the Mahalanobis distance, the outlying measurement is identified. Meanwhile, the covariances of the outlying measurement noises and the normal measurements noises will be updated through the Equations (14) and (15), respectively. Both the judging index and are computed using the current measurements, and this indicates a strong sensitivity towards the abnormal measurements and a relatively light computational burden. Accordingly, a modified adaptive data fusion scheme is constructed. In this proposed scheme, the hypothesis test based on the Mahalanobis distance is implemented at every epoch. For each epoch, if the null hypothesis is accepted, both the fading filter and the sequential updating of are implemented. Otherwise, only the equivalent covariance matrix is performed, and this strategy is conducted adaptively. Therefore, a modified scheme is developed concerning the circumstance when the current measurement information is unreliable.

4. Models for GNSS/INS Integrated Systems

The GNSS/INS integrated systems have become a well applied technique in the field of dynamic navigation and positioning [25, 26]. GNSS and INS are integrated mainly in three types, and the loosely coupled type is applied in this paper. For the GNSS/INS and various other integrated systems, the nonlinear filter should be implemented aiming at the nonlinear problem. The cubature Kalman filter (CKF) was proposed to address the high-dimensional state estimation [27]. With a fifteen-dimension state vector , including the position error , the velocity error , the attitude error of the carrier and the bias of the gyroscope , and the accelerometer , the CKF is adopted in this paper, where means the earth frame and means the body frame. Elements of are listed below.

For the nonlinear discrete system where and are the nonlinear functions. In the CKF, the time updating process is derived by [27, 28]:

And the measurement updating process is derived by where denotes the cubature points generated from states and , denotes dimension of the state vector, denotes the square root of , , and and denote the propagated cubature points from states and measurements, respectively. Differences of the carrier position and velocity between GNSS and INS, namely, and , are adopted in the process of measurement update. Therefore, the vector of the measurement is constructed below.

Actually, the measurement equation is linear, and the measurement updating process should be modified as follows.

Obviously, the precision of the estimates relies largely on the quality of the measurements, and unreliable measurements may result in inestimable influences. In practice, however, unreliable measurements are inevitable. Therefore, the abnormal measurements must be deliberately processed. Both the fading factor or matrix and the covariance inflation factor have been adopted in the data processing of the GNSS/INS integrated systems. Aiming at the influences on the fading factor, the judging index constructed based on the Mahalanobis distance and the covariance inflation factor is implemented to improve the performance of the filtering process.

5. Performance Evaluation and Analysis

Data processing experiments using different filtering algorithms were designed and implemented. Equipped with the GNSS receiver (Trimble R8) and the inertial measurement unit (IMU, SPAN-CPT), the land vehicle was used as the testing carrier, and a three-axis open-loop gyroscope and three-axis MEMS accelerometers are included in the IMU. The average speed of the land vehicle was about 3.5 m/s. Meanwhile, another GNSS receiver called the base station was set on top of the building. Motion trajectory of the land vehicle (in the North and East directions) is demonstrated in Figure 1. The main technical parameters provided by the suppliers were listed in Table 1. Then, the data was collected with the land vehicle under real conditions. Position and velocity from GNSS were computed using double-difference pseudo-range measurement, and the position and velocity deviations were 0.25 m2 and (0.05 m/s)2, respectively. In the testing systems, the cut-off angle was set by 10°, and the sampling frequency was 1 Hz for GNSS and 100 Hz for INS. In this paper, the tightly coupled results provided by the commercial software IE using double-difference carrier phase measurements were regarded as the references.

As to data fusion of the GNSS/INS integrated systems of the paper, the CKF was adopted as the basic filter, and time interval of the data fusion was one second. In the experiments, both the initial data and the data with artificially added outliers were adopted to test the performances of different algorithms more clearly. Consequently, two cases were designed and performed. Each case consists of four algorithms, and the differences between the references and the results of each algorithm are thought to be errors. Comparing the fading filters with a fading factor and fading matrix, the fading matrix performs better by adjusting the covariance matrix in multiple data channels simultaneously, and the ability to resist the model mismatch is enhanced. Therefore, the fading matrix was adopted in this paper. For all algorithms of this paper, the initial value of the covariance are fixed by experience. Thus, the statistical deviations of the models exist, and it is suitable to test the proposed estimation algorithm. Four different algorithms of each case are designed below.

Algorithm 1. The conventional CKF.

Algorithm 2. The fading filter with fading matrix (MF-CKF).

Algorithm 3. The fading filter with fading matrix is implemented at the epochs when the null hypothesis is accepted (Partial-MF-CKF).

Algorithm 4. The fading filter with fading matrix is implemented while the null hypothesis is accepted; otherwise, the covariance inflation factor is adopted for outlying measurements and sequential updating for normal measurements (RMF-CKF).

Case 5. In this case, the initial data collected through the self-developed GNSS/INS integrated navigation systems is processed with the above four algorithms, respectively. Positioning errors of the land vehicle of each algorithm are demonstrated in Figures 25.

As mentioned above, each algorithm was implemented based on the initial measurements collected under ideal observation conditions, and little outliers exist in the measurements. Accordingly, positioning errors were mainly brought by the model errors. Comparing Figures 25, apparently, error amplitudes in , and directions of the CKF algorithm are much bigger than the other algorithms, and this indicates that the performance of the conventional CKF algorithm is improved significantly. Moreover, stability of the filter was improved, and the filter divergence resulted from the model uncertainties was restrained effectively with the last three algorithms. In theory, main differences between the Algorithms 2 and 3 lie on the circumstance when the fading matrix is adopted. With the initial measurements, the null hypothesis was accepted at most epochs; thus, the results of the Algorithms 2 and 3 are similar. In the MF-CKF algorithm, the fading matrix was implemented at all epochs, and the model errors were well-considered. Except for the epochs when the null hypothesis was accepted, the covariance inflation factor in the RMF algorithm was adopted at other epochs, and the influences of abnormal measurements were weakened. Consequently, error amplitudes of the processed algorithm displayed in Figure 4 are smaller than those of other three algorithms.

The root mean square error (RMSE) of these positioning errors was computed to examine the performances of each algorithm quantificationally. Figure 6 displays RMSEs in , and directions of each algorithm. Meanwhile, the detailed RMSEs (position , velocity , and attitude (Yaw) are taken as examples) of these four algorithms are listed in Table 2.

In terms of RMSE values in Figure 6 and Table 2, CKF is much bigger than other three algorithms, and this indicates that the robustness and stability of the conventional CKF should be enhanced further. Comparing the RMSE values of MF-CKF and Partial-MF-CKF algorithms, the former performed better since the influences of model errors were more significant than those of abnormal measurements. Based on the Partial-MF-CKF, the covariance inflation factor was implemented at the epochs when the null hypothesis was rejected. Thus, superiorities of the Mahalanobis distance, the fading filter, and the robust estimation are combined in the proposed algorithm. Meanwhile, RMSE values of the other schemes are all bigger than those of the RMF-CKF algorithm, and this indicates that better performance is achieved with the processed algorithm. In Table 2, obviously, RMSEs of all the fading filtering algorithms are smaller than the CKF algorithm since effects of the model errors are weakened, and the RMF-CKF algorithm outperforms the other algorithms in terms of accuracy. Compared with the conventional MF algorithm, the accuracy of positioning is improved for 17.7% and 13.0% in and directions, respectively, using the RMF-CKF algorithm.

Case 6. In this case, the continuously changing and the separated positioning outliers were introduced artificially into the measurements to test the robustness and stability of each algorithm. Thus, the perturbation data based on the initial data was constructed. All the above four algorithms were implemented based on the perturbation data, and positioning errors of each algorithm are illustrated in Figures 710.

Since positioning outliers were introduced into the measurements, the outliers, instead of the model errors, became the principal factor in affecting the filtering performance. It is demonstrated in Figures 710 that amplitudes of all algorithms were much bigger than those of Case 1, and this denotes that all algorithms were affected significantly by the outlying measurements. Meanwhile, as to error amplitudes, all the algorithms with fading matrix are better than the CKF algorithm. Since the outlying measurements were not addressed effectively in the CKF and Partial-MF-CKF algorithms, Figures 7 and 9 illustrate that effects of the outlying measurements are more apparent. As shown in Figures 7 and 9, the fading matrix was implemented when the null hypothesis was accepted and the model errors were well controlled; thus, the Partial-MF-CKF algorithm performed superior than the CKF. In the MF-CKF algorithm, the fading matrix is implemented at all epochs; the fading matrix performs better when the null hypothesis is accepted. However, the fading filter becomes instable when the null hypothesis is rejected. Therefore, once the fading matrix is applied mistakenly, the MF-CKF algorithm may become inferior to the CKF algorithm. Comparing Figures 7 and 8, it is illustrated that amplitudes at some epochs of the MF-CKF algorithm are even bigger than those of the CKF algorithm, and this is the influences brought from the outlying measurements. In other words, implementing the fading matrix mistakenly may affect the filtering performance at the current together with the future epochs. In the RMF-CKF algorithm, both the model errors and the outlying measurements were considered. Meanwhile, the fading matrix was implemented at the relatively reasonable epochs. In Figure 10, amplitudes in three directions are much smaller than the other algorithms, and a better performance is achieved with the RMF-CKF algorithm. However, as to the constantly changing outliers, the values may small enough to go through the hypothesis test, and the fading matrix will be implemented mistakenly. Therefore, at some epochs in Figure 10, performance of the RMF-CKF algorithm was still apparently affected by the constantly changing outliers.

Figure 11 depicts RMSEs of each algorithm, and the detailed RMSEs of these four algorithms are provided in Table 3.

Comparing the RMSEs in Figures 6 and 11, it is concluded that all algorithms were affected by the outlying measurements. Under the effects of outlying measurements, the model errors were weakened effectively; thus, the MF-CKF algorithm performed better than the CKF algorithm. Although both the MF and Partial-MF algorithms were affected significantly by the outlying measurements, the fading matrix was implemented only when the null hypothesis was accepted, and the negative effects from mistakenly used fading matrix were avoided in the Partial-MF-CKF algorithm. Therefore, the Partial-MF-CKF algorithm performed better than the MF-CKF algorithm. Corresponding to the qualitative conclusions, RMSEs of the RMF-CKF algorithm are much smaller than the other algorithms; apparently, better performance is achieved with the proposed algorithm. Comparing Tables 2 and 3, it is demonstrated that the performances of all algorithms were degraded by the outlying measurements, and the filters with fading matrix are more stable and robust than CKF algorithm. By integrating the superiorities of the multiple fading filter and the robust estimation method, the smaller RMSE values and a relatively stable performance were obtained with the proposed algorithm. In this circumstance, the accuracy of positioning are improved for 60.1% and 58.5% in and directions, respectively, using the RMF-CKF algorithm, and this in turn demonstrates that the MF algorithm is affected by the outlying measurements seriously.

6. Conclusions

In this paper, a rarely discussed problem about how to apply the fading factor or matrix with unreliable measurements adaptively is researched, and a modified adaptive filter and an alternative scheme are developed. In the proposed scheme, the hypothesis test was performed using the statistics derived from the Mahalanobis distance, and the actual measurements were collected with self-developed GNSS/INS integrated systems; then, the contrastive experiments and analysis were implemented to test the performance of the proposed algorithm. Detailed conclusions of this paper are summarized below: (1)Compared with the conventional CKF, the adaptive filter with the fading matrix can control the negative effects of model uncertainties; however, the performance may be degraded significantly when there exist outlying measurements(2)In the fading filter, unreliable measurements will construct unreliable fading matrix, and performance of the filter may be even inferior to the conventional Kalman filter if the fading matrix is applied mistakenly. Therefore, it is suggested that the fading matrix be applied in a timely and more reasonable way(3)Statistics derived from the Mahalanobis distance can be adopted to test the quality of measurements, and stability of the fading filter is improved with the hypothesis test. By integrating the advantages of the fading matrix and the covariance inflation factor, effects of both the model errors and the outlying measurements are controlled using the proposed scheme for the systems with model uncertainties. Also, the proposed adaptive fusion scheme is more suitable when there exist outliers in the measurements

Data Availability

Some or all data, models, or code that support the findings of this study are available from the corresponding author upon reasonable request (email: [email protected]).

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This research project was supported by the Natural Science Foundation of Henan Province (212300410198), the Nation Key Research Program of China (2017YFE0131400), and the State Key Laboratory of Geodesy and Earth’s Dynamics (SKLGED 2020-3-1-E).