Next Article in Journal
Fault Diagnosis Method for a Mine Hoist in the Internet of Things Environment
Previous Article in Journal
Using Temporal Covariance of Motion and Geometric Features via Boosting for Human Fall Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Strong Tracking Cubature Kalman Filter for GPS/INS Integrated Navigation Systems

1
Key Laboratory of Instrumentation Science & Dynamic Measurement, Ministry of Education, North University of China, Taiyuan 030051, China
2
National Key Laboratory for Electronic Measurement Technology, North University of China, Taiyuan 030051, China
*
Author to whom correspondence should be addressed.
Submission received: 8 May 2018 / Revised: 8 June 2018 / Accepted: 8 June 2018 / Published: 12 June 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
The cubature Kalman filter (CKF) is widely used in the application of GPS/INS integrated navigation systems. However, its performance may decline in accuracy and even diverge in the presence of process uncertainties. To solve the problem, a new algorithm named improved strong tracking seventh-degree spherical simplex-radial cubature Kalman filter (IST-7thSSRCKF) is proposed in this paper. In the proposed algorithm, the effect of process uncertainty is mitigated by using the improved strong tracking Kalman filter technique, in which the hypothesis testing method is adopted to identify the process uncertainty and the prior state estimate covariance in the CKF is further modified online according to the change in vehicle dynamics. In addition, a new seventh-degree spherical simplex-radial rule is employed to further improve the estimation accuracy of the strong tracking cubature Kalman filter. In this way, the proposed comprehensive algorithm integrates the advantage of 7thSSRCKF’s high accuracy and strong tracking filter’s strong robustness against process uncertainties. The GPS/INS integrated navigation problem with significant dynamic model errors is utilized to validate the performance of proposed IST-7thSSRCKF. Results demonstrate that the improved strong tracking cubature Kalman filter can achieve higher accuracy than the existing CKF and ST-CKF, and is more robust for the GPS/INS integrated navigation system.

1. Introduction

GPS/INS integration navigation systems have drawn great attention due to their widespread applications in the areas of dynamic navigation and positioning [1,2,3]. Depending on the information fusion level, three main forms of integration structures have been developed: loosely coupled [4], tightly coupled [5] and deeply coupled [6]. In this paper, only the loosely coupled integrated navigation system is considered due to its easier implementation and appropriate performance. For GPS/INS integration, the most celebrated mechanism is the Kalman type filter and it has been widely used in the practical systems. In the Kalman type filter, the optimal estimation can be achieved when the distributions of both process noise and measurement noise are assumed to be Gaussian and the corresponding statistics are accurately known [7]. However, due to the influence of vehicle’s severe maneuvers in the actual applications, the Gaussian assumption may be violated and the statistics may be uncertain, which will result in a degraded performance of the conventional KF-based GPS/INS integration system [8,9]. Therefore, in the presence of process uncertainties, there is a great demand to develop an effective filter for better suppressing them, which is the main focus of this work.
To address the filtering problems with process uncertainty effectively, a large number of fruitful algorithms have been investigated, such as the particle filter (PF) [10], the variational Bayesian-based adaptive Kalman filter (VB-AKF) [11,12,13,14], the maximum correntropy-based Kalman filter (MC-KF) [15,16,17,18,19] and the Huber’s M-Estimation-based Kalman filter (HKF) [8]. The particle filter is based on the sequence Monte Carlo sampling technique and can solve the nonlinear non-Gaussian system state estimation problem effectively [10]. It is, however, difficult to use the particle filter to cope with the high dimensional filter problem since the computational complexity increases exponentially with the dimensions of the state. The variational Bayesian-based adaptive Kalman filter can estimate the statistics of the time-varying process noise variance online by approximating the joint posterior distribution of the latent state and noise variance by a factorized free form distribution. However, the dynamic model designed for variance parameters used in the VB-AKF is heuristic and relatively rough, which may degrade its estimation accuracy when the process noise covariance changes frequently [14]. The maximum correntropy-based Kalman filter is capable of suppressing the larger process uncertainty by maximizing the correntropy of the predicted error and residual. However, it is very difficult to determine the appropriate value of kernel bandwidth size in actual applications, which may degrade the estimation performance of the MC-KF [16]. The Huber’s M-Estimation-based Kalman filter handles the process uncertainty by employing the cost function of M-estimator to construct weighting matrix, which can rescale the prior state estimate covariance as a result the process uncertainty are suppressed. Unfortunately, the influence function of the HKF doesn’t redescend, which results in limited estimation accuracy [20].
Another way of solving the problems of process uncertainty is to use adaptive filters where the online estimates of process noise statistics are obtained together with the dynamic state. The adaptive filters fall into two categories, i.e., the Multiple Model Adaptive Kalman filter (MM-AKF) [21] and the Innovation-based Adaptive Kalman filter (IAKF) [22]. The multiple model adaptive Kalman filter handles the process uncertainty by running a bank of Kalman filters with different stochastic models in parallel. However, it is only suitable for systems with known dynamics since it operates under the assumption that one of the models in the model bank is the correct one [22,23]. The Innovation-based Adaptive Kalman filter can estimate the appropriate covariance matrix of process noise by forcing the innovation sequence of the Kalman filter to be white Gaussian noise sequence with zero mean. However, this algorithm suffers from the limitation that the convergence to the right process noise covariance matrix can’t be guaranteed and rather large windows of data are required to achieve a reliable estimation of process noise covariance matrix, making it only suitable for slowly changing systems [22].
A feasible approach to deal with the process uncertainty is to employ the strong tracking filter [24,25,26]. The principle of strong tracking filter is that the strong tracking factor can maintain the residual sequence orthogonal by introducing the time-variant suboptimal fading factor matrix to the state prediction covariance matrix. As a result, the strong tracking filter provides a strong robustness against process uncertainty. However, the conventional strong tracking filter may suffer from the following problems: (1) The suboptimal fading factor is incorporated in the whole filtering process, which may result in the loss of precision when there is no process uncertainty existed; (2) The symmetry of the prediction covariance matrix cannot be guaranteed when the suboptimal fading factor with different diagonal elements is carried out on the error covariance matrix, which may degrade the filtering performance and even lead to the divergence.
Thus, in order to solve the issues mentioned above, an improved strong tracking cubature Kalman filter that could be applied to solve the high dimensional GPS/INS integrated navigation problem is proposed in our paper. First, the hypothesis testing method is adopted to identify process model uncertainty, which avoids the loss of precision when there is no process uncertainty in the dynamic system. In addition, a defined suboptimal fading factor, which is based on the idea of Cholesky decomposition, is proposed to ensure the symmetry of error covariance matrix when the fading factor is introduced. Second, the seventh-degree spherical simplex-radial cubature rule, which is used to compute the more exact Gaussian type integral, is integrated into the strong tracking cubature Kalman filter to improve the estimation accuracy and numerical stability of the GPS/INS integrated navigation system. Finally, the car-mounted experiment is conducted to evaluate the effectiveness of the proposed approach. The experiment results indicate that the proposed IST-7thSSRCKF method has better robustness for the suppression of process uncertainties as compared with CKF and STCKF for the GPS/INS integrated navigation problem.
The reminder of this paper is organized as follows: the strong tracking filter and cubature Kalman filter are briefly reviewed in Section 2. Section 3 presents the improved STCKF for SINS/GPS integration. In Section 4, experiments are carried out to evaluate the performance of the proposed approach. Concluding remarks are given in Section 5.

2. The Strong Tracking Filter and Cubature Kalman Filter

In this section, the strong tracking filter and cubature Kalman filter are briefly reviewed and discussed respectively. Consider the following nonlinear discrete-time dynamic system with additive process and measurement noise:
{ x k = f ( x k 1 ) + w k 1 z k = h ( x k ) + v k
where x k n × 1 , z k m × 1 are the state vector of the dynamic system and measurement vector at discrete time k ; f ( ) denotes the nonlinear dynamic model and h ( ) represents the measurement model. w k 1 and v k are assumed to be independent process and measurement Gaussian noise sequences with zero means and variances Q k 1 and R k , respectively.

2.1. Strong Tracking Filter

The strong tracking filter, proposed by Zhou et al. [27], is essentially a kind of Kalman filter which is based on the orthogonal principle of the output residual error. In the strong tracking filter, the time-varing suboptimal fading factor λ k is introduced to the predicted error covariance matrix P k | k 1 , which adjusts the gain matrix K k online in real time, making sure that the innovations remain orthogonal and the available information in the residuals ε k can be extracted completely. In this way, the algorithm has strong robustness against process uncertainty. The general structure of the strong tracking filter algorithm can be expressed as follows:
ε k = z k h ( x ^ k | k 1 )
x ^ k | k = x ^ k | k 1 + K k ε k
where x ^ k | k is the estimate value at discrete time k , ε k denotes the residual sequence of measurement. The question of strong tracking filter is how to adjust the time-varing gain matrix online according to the measurement innovation. As we known, if there is no process uncertainty in the system model, the strong tracking filter boils down to the standard Kalman filter and the time-varing gain matrix K k can be determined directly as follows:
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 = F k | k 1 P k 1 F k | k 1 T + Q k 1
where F k | k 1 = f x | x = x ^ k 1 and H k = h x | x = x ^ k | k 1 are the system transition matrix and observation matrix, respectively. However, the process uncertainty in actual application results in a deviation between the state estimation and actual state, leading to the nonorthogonality among the residuals. Actually, the nonorthogonality indicates that the available information in the measurement is not fully used in the filtering estimation. Thus, the suboptimal fading factor λ k is introduced to the state predicted error covariance matrix, which will deliberately decrease the effect of historical information in the process estimation and keep the residual sequence mutually orthogonal. The predicted error covariance matrix after incorporating the fading factor can be expressed as follows:
P k | k 1 = λ k F k | k 1 P k 1 F k | k 1 T + Q k 1
where the suboptimal fading factor λ k is determined by the orthogonal principle and a popular choice of this can be given by [28]:
λ k = { α c k ,       α c k 1     1 ,          α c k < 1 , c k = t r [ N k ] t r [ M k ]
where α ( α 1 ) denotes the proportional coefficient of fading factor λ k , which is a constant and can be predetermined by a prior knowledge. t r [ N k ] and t r [ M k ] represent the trace of the matrix N k and M k , respectively. N k and M k are defined as follows:
N k = V k H k Q k 1 Η k T β R k
M k = H k F k | k 1 P k 1 F k | k 1 T H k T
V k = { ε 1 ε 1 T ,                                        k = 0 ρ V k 1 + v k v k T 1 + ρ ,     k 1
where V k denotes the residual error sequence covariance matrix. β denotes the softening factor of the residual sequence, which is used to improve the smoothness and accuracy of the state estimation, generally β = 4.5 [28]. ρ is the fading factor satisfying 0.95 ρ < 0.995 , usually, we set ρ to be 0.95 [29].

2.2. Cubature Kalman Filter

The cubature Kalman filter was first presented by Haykin to solve high-dimensional nonlinear filtering problems [30]. The heart of cubature Kalman filter is to use the third-degree cubature rule to approximate the integral of form (nonlinear function × Gaussian density) instead of approximating the nonlinear function. Based on the third-degree cubature rule, a set of 2 n cubature points with the same weight are chosen to propagate the state and covariance matrix at each update cycle. The integral with respect to the general Gaussian weighted integral can be approximated by the cubature rule as follows:
I ( f ) = n f ( x ) N ( x ; x ¯ , P ) d x = n f ( P x + x ¯ ) N ( x ; 0 , I ) d x i = 1 m ω i f ( P ξ i + x ¯ )
where P = P P T ; P can be calculated by the Cholesky decomposition or the singular value decomposition. f ( x ) is the arbitrary nonlinear function, n is the integral domain, x ¯ is the mean of x , m is the total number of cubature points. { ξ i , ω i } are the i th cubature point and the corresponding weight, respectively, which are given by:
ξ i = m 2 [ 1 ] i
ω i = 1 m , i = 1 , 2 , m = 2 n
where [ 1 ] i denotes the i th column vector of the points set [ 1 ] . For example, [ 1 ] 2 represents the following set of points:
{ ( 1 0 ) , ( 0 1 ) , ( 1 0 ) , ( 0 1 ) }
Under the assumption that the initial state x ^ 0 | 0 and initial error covariance matrix P ^ 0 | 0 are known, the computational steps of the cubature Kalman filter can be summarized as follows [30]:

2.2.1. Prediction

(1)
Factorize the covariance and evaluate the cubature points ( i = 1 , 2 , 2 n ):
P k 1 | k 1 = S k 1 | k 1 S k 1 | k 1 T
X i , k 1 | k 1 = S k 1 | k 1 ξ i + x ^ k 1 | k 1
(2)
Evaluate the propagated cubature points through the process model:
X i , k | k 1 = f ( X i , k 1 | k 1 )
(3)
Estimate the predicted state and the corresponding error covariance:
x ^ k | k 1 = 1 2 n i = 1 2 n X i , k | k 1
P k | k 1 = 1 2 n i = 1 2 n X i , k | k 1 X i , k | k 1 T x ^ k | k 1 x ^ k | k 1 T + Q k 1

2.2.2. Update

(1)
Factorize the covariance and evaluate the cubature points ( i = 1 , 2 , 2 n ):
P k | k 1 = S k | k 1 S k | k 1 T
X i , k | k 1 = S k | k 1 ξ i + x ^ k | k 1
(2)
Evaluate the propagated cubature points through the observation model:
Z i , k | k 1 = h ( X i , k | k 1 )
(3)
Estimate the predicted measurement:
z ^ k | k 1 = 1 2 n i = 1 2 n Z i , k | k 1
(4)
Estimate the covariance and Kalman gain:
P z z , k | k 1 = 1 2 n i = 1 2 n Z i , k | k 1 Z i , k | k 1 z ^ k | k 1 z ^ k | k 1 T + R k
P x z , k | k 1 = 1 2 n i = 1 2 n X i , k | k 1 Z i , k | k 1 T x ^ k | k 1 z ^ k | k 1 T
K k = P x z , k | k 1 / P z z , k | k 1
(5)
Estimate the updated state and the corresponding error covariance:
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

3. An Improved Strong Tracking 7thSSRCKF Algorithm

3.1. The Improvement of Strong Tracking Kalman Filter

In order to improve the tracking performance of the traditional strong tracking filtering algorithm for multi-variable nonlinear system, an improved time-varying fading factor is proposed in this section. In the proposed improved strong tracking Kalman filter, the hypothesis testing method is adopted to identify process model uncertainty, which avoids the loss of precision when there is no process uncertainty in the dynamic system. In addition, a defined suboptimal fading factor, which is based on the idea of least squares and Cholesky decomposition, is proposed to ensure the symmetry of error covariance matrix when the fading factor is introduced.

3.1.1. Process Uncertainty Identification

For the dynamic system without process uncertainty, the innovation ε k should be zero-mean Gaussian-distribution with covariance H k P k | k 1 H k T + R k when the filtering is stable [31]. So the square of the Mahalanobis distance of the innovation should obey a χ 2 distribution [32], then we define the statistical information of predicted residual vector as follows:
γ k = ε k T [ H k ( F k | k 1 P k 1 F k | k 1 T + Q k 1 ) H k T + R k ] 1 ε k χ 2 ( m )
where γ k obeys a Chi-square distribution with m degrees of freedom, m is the number of state variable that can be observed directly. According to the hypothesis testing theory, for a chosen significant level α , we have:
P ( χ 2 < χ α , m 2 ) = 1 α
where the α -quantile χ α , m 2 of the Chi-square distribution is predetermined. In this paper, the value is set α = 0.05 ,   χ α , m 2 = 12.592 . If the actual γ k is less than this α -quantile, i.e., γ k < χ α , m 2 , the standard cubature Kalman filter is carried out. Otherwise, if γ k χ α , m 2 , it can be conclude with high probability ( 1 α ) that there exists process uncertainty in the dynamic system. In this case, the improved strong tracking cubature Kalman filter is used.

3.1.2. Improved Strategy for Fading Factor

In order to make different channel have different fading rate and respective filter adjustment capability in the multi-dimensional GPS/INS integrated navigation systems, a defined multiple time-varing suboptimal fading factor matrix Λ k to the prediction error covariance matrix is introduced:
P k | k 1 = Λ k F k | k 1 P k 1 F k | k 1 T + Q k 1
where Λ k = d i a g { 1 , 1 , 1 , λ 4 , k , λ 5 , k , , λ 9 , k , 1 , 1 } is the matrix of multiple time-varying suboptimal fading factor. Specifically, { λ 4 , k , λ 5 , k , , λ 9 , k } is used to adjust the state variable that can be directly observed in the GPS/INS integrated navigation system which includes the velocity error and position error, whereas the other elements of Λ k is set as 1 since they cannot be estimated adaptively due to their unobservablities. The defined multiple time-varing fading factor Λ k can be determined as follows:
According to the orthogonality principle [31], in order to satisfy E [ ε k + j ε k T ] = 0 , we need to make the following equation hold:
P k | k 1 H k T K k V k 0
Substituting (4) into (32) yields:
P k | k 1 H k T P k | k 1 H k T ( H k P k | k 1 H k T + R k ) 1 V k 0
Simplifying Equation (33), we have:
H k P k | k 1 H k T = V k R k
Substituting (31) and the fading factor Λ k into (34) yields:
H k Λ k F k | k 1 P k 1 F k | k 1 T H k T = V k H k Q k 1 H k T R k
To be specific, for the GPS/INS integrated navigation system with the integration mode of velocity and position, the observation matrix can be expressed as:
H k = [ 0 6 × 3 I 6 × 6 0 6 × 6 ]
Here we define:
Δ k = H k Λ k = [ 0 6 × 3 λ 4 , k 0 0 0 0 0 0 λ 5 , k 0 0 0 0 0 0 λ 6 , k 0 0 0 0 0 0 λ 7 , k 0 0 0 0 0 0 λ 8 , k 0 0 0 0 0 0 λ 9 , k 0 6 × 6 ]
G k = F k | k 1 P k 1 F k | k 1 T H k T
Substituting (37) and (38) into (35) yields:
Δ k G k = N k
i.e.,:
[ 0 6 × 3 λ 4 , k 0 0 0 0 0 0 λ 5 , k 0 0 0 0 0 0 λ 6 , k 0 0 0 0 0 0 λ 7 , k 0 0 0 0 0 0 λ 8 , k 0 0 0 0 0 0 λ 9 , k 0 6 × 6 ] G k = N k
Using the method of least squares we have:
λ 4 , k = j = 1 6 G k [ 4 , j ] N k [ 1 , j ] j = 1 6 G k 2 [ 4 , j ]
Similarly, we obtain the λ 5 , k λ 9 , k . Then the multiple suboptimal fading factor Λ k for the low-cost GPS/INS integration navigation system can be estimated adaptively as:
Λ i , k = { j = 1 6 G k [ i , j ] N k [ i 3 , j ] j = 1 6 G k 2 [ i , j ] , i = 4 , 5 , , 9 1   ,   i = 1 , 2 , 3 , 10 , , 21
Unfortunately, the symmetry of the prediction covariance matrix P k | k 1 in (31) cannot be guaranteed when the suboptimal fading factor Λ k with different diagonal elements, i.e., λ 4 , k λ 5 , k λ 9 , k , is carried out on the matrix F k | k 1 P k F k | k 1 T . To solve the problem, the square root filter method is considered. Based on the idea of Cholesky decomposition, the fading factor matrix Λ k can be decomposed into:
Λ k = Λ ^ k · Λ ^ k T
where Λ ^ k = d i a g { 1 , 1 , 1 , λ 4 , k , λ 5 , k , , λ 9 , k , 1 , 1 } . Hence, by introducing the improved two multiple fading factor Λ ^ k , the prediction error covariance matrix can be rewritten as:
P k | k 1 = Λ ^ k F k | k 1 P k F k | k 1 T Λ ^ k T + Q k 1
By using Equation (44), the properties of the prediction error covariance matrix such as symmetric positive definite can be ensured and the stability of the filter can be therefore enhanced.

3.2. The Seventh-Degree Spherical Simplex-Radial Cubature Rule for Cubature Kalman Filter

The cubature Kalman filter based on the third-degree spherical-radial cubature rule has been successfully applied in many practical applications. However, it still has limited estimation accuracy. For example, the cubature Kalman filter can’t calculate exactly the Gaussian weighted integral of such simple polynomial functions as x 1 2 x 2 2 , where x 1 and x 2 are the two components of the Gaussian random vector function [33]. In order to overcome the shortcoming and further improve the accuracy and efficiency of the traditional CKF, a high-degree cubature Kalman filter based on the seventh-degree spherical simplex-radial rule is developed in this section.
In the seventh-degree spherical simplex-radial cubature rules, the following integral is considered:
I ( f ) = n f ( x ) exp ( x T x ) d x
where f ( x ) is some arbitrary nonlinear function, n is the domain of integration, and exp ( x T x ) is the weighting function. Let x = r s with s T s = 1 and r = x T x , Equation (45) can be rewritten in the spherical-radial coordinate system as:
I ( f ) = 0 U n f ( r s ) r n 1 exp ( r 2 ) d σ ( s ) d r
where U n = { s n | s T s = 1 } is the spherical surface and σ ( ) is the spherical surface measure or the area element on U n . Then, the Equation (46) can be decomposed into the spherical integral S ( r ) = U n f ( r s ) d σ ( s ) with the weighting function ω f ( s ) = 1 , and the radial integral 0 S ( r ) r n 1 exp ( r 2 ) d r with the weighing function ω f ( r ) = r n 1 exp ( r 2 ) .

3.2.1. Seventh-Degree Spherical Simplex Rule

The spherical integral S ( r ) = U n f ( r s ) d σ ( s ) can be derived based on the transformation group of the regular simplex, with vertices a j = [ a j , 1 , a j , 2 , , a j , n ] T , j = 1 , 2 , , n + 1 as follows [34]:
a j , i { 0 , i > j ( n + 1 ) ( n j + 1 ) n ( n j + 2 ) , i = j n + 1 n ( n i + 2 ) ( n i + 1 ) , i < j
For example, when the dimension n = 2, we can obtain that { a j } = { [ 1 , 0 ] T , [ 1 / 2 , 3 / 2 ] T , [ 1 / 2 , 3 / 2 ] T } . Based on the projection from the midpoints of the vector a j , we can obtain the following point sets [34]:
{ b j } { n 2 ( n 1 ) ( a i + a j ) : i < j }
{ c j } { n 3 ( n 2 ) ( a i + a j + a m ) : i < j < m , m = 1 , 2 , , n + 1 }
{ d j } { n 10 n 6 ( a i + 3 a j ) : i j }
where { b j } are the set of midpoints of the vertices a j projected onto the surface of the sphere U n , { c j } are projections of the centroid of vertices a j faces onto U n , and { d j } are projections of the selected edge points of the vertices a j onto the surface of U n .
Taking the sets { a j } , { b j } , { c j } and { d j } as cubature points and further employing the central symmetry of the cubature formula, we can obtain the seventh-degree spherical simplex rule with ( n + 1 ) ( n 2 + 8 n + 6 ) / 3 cubature points as follows [34,35]:
S 7 ( r ) = | A n | 36 n ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) { n 3 ( 9 n 2 793 n + 1800 ) j = 1 n + 1 [ f ( r a j ) + f ( r a j ) ] + 144 ( n 1 ) 3 ( 4 n ) j = 1 n ( n + 1 ) / 2 [ f ( r b j ) + f ( r b j ) ] + 486 ( n 2 ) 3 j = 1 ( n 1 ) n ( n + 1 ) / 6 [ f ( r c j ) + f ( r c j ) ] + ( 10 n 6 ) 3 j = 1 n ( n + 1 ) [ f ( r d j ) + f ( r d j ) ] }
where A n = 2 Γ ( 1 2 ) / Γ ( n 2 ) = 2 π n / Γ ( n 2 ) is the surface area of the unit sphere, the Gamma function is defined as Γ ( n ) = 0 x n 1 e x d x with the properties of Γ ( 1 / 2 ) = π and Γ ( n + 1 ) = n Γ ( n ) .

3.2.2. Seven-Degree Radial Rule

Based on the idea of moment matching, the radial integral can be calculated by the following moment equation:
0 S ( r ) r n 1 e r 2 d r i = 1 N r w r , i S ( r i )
where r i and w r , i denote the quadrature point and the weights, respectively. N r is the total number of points; S ( r ) = r l is a monomial in r , with l an even integer. Then, the left hand side of the Equation (52) is reduced to 1 2 Γ ( n + l 2 ) with Γ ( n ) = 0 x n 1 e x d x , i.e.,:
0 r l r n 1 e r 2 d r = 1 2 Γ ( n + l 2 )
Because the spherical rule and the resultant spherical-radial cubature rule are fully symmetric, we only need to match the even-degree monomials. The quadrature point r i and w r , i can therefore be obtained by solving the different even-degree monomials matching. For the seventh-degree radial rule, we require that the points and weights satisfy the following four equations:
{ w r , 1 r 1 0 + w r , 2 r 2 0 = 1 2 Γ ( n 2 ) w r , 1 r 1 2 + w r , 2 r 2 2 = 1 2 Γ ( n 2 + 1 ) = n 4 Γ ( n 2 ) w r , 1 r 1 4 + w r , 2 r 2 4 = 1 2 Γ ( n 2 + 2 ) = 1 2 ( n 2 + 1 ) ( n 2 ) Γ ( n 2 ) w r , 1 r 1 6 + w r , 2 r 2 6 = 1 2 Γ ( n 2 + 3 ) = 1 2 ( n 2 + 2 ) ( n 2 + 1 ) ( n 2 ) Γ ( n 2 )
By solving the above four equations of zeroth, second, fourth, and sixth order moments, the points and weights for the seventh-degree radial rule with 2 quadrature points can be obtained as follows:
{ r 1 = 1 2 4 + 2 n + 2 4 + 2 n r 2 = 1 2 4 + 2 n 2 4 + 2 n
{ w r , 1 = ( 1 4 1 2 4 + 2 n ) Γ ( n 2 ) w r , 2 = ( 1 4 + 1 2 4 + 2 n ) Γ ( n 2 )
A detailed calculation of the points and weights can be found in [33]. Then, the seventh-degree radial rule can be obtained as follows:
0 S ( r ) r n 1 e r 2 d r = w r , 1 S ( r 1 ) + w r , 2 S ( r 2 )
where w r , 1 , r 1 , w r , 2 , r 2 are given by Equations (55) and (56).

3.2.3. Seventh-Degree Spherical Simplex-Radial Rule

Combining Equations (51) and (57), the seventh-degree Spherical Simplex-Radial (SSR) rule can be formulated as:
n f ( x ) N ( x ; x ¯ , P x ) d x = ( 1 2 1 4 + 2 n ) n 2 ( 9 n 2 793 n + 1800 ) 36 ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) j = 1 n + 1 [ f ( x ¯ ( n + 2 + 4 + 2 n ) P x a j ) + f ( x ¯ + ( n + 2 + 4 + 2 n ) P x a j ) ] + ( 1 2 1 4 + 2 n ) 4 ( n 1 ) 3 ( 4 n ) n ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) k = 1 n ( n + 1 ) / 2 [ f ( x ¯ ( n + 2 + 4 + 2 n ) P x b k ) + f ( x ¯ + ( n + 2 + 4 + 2 n ) P x b k ) ] + ( 1 2 1 4 + 2 n ) 486 ( n 2 ) 3 36 n ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) k = 1 ( n 1 ) n ( n + 1 ) / 6 [ f ( x ¯ ( n + 2 + 4 + 2 n ) P x c k ) + f ( x ¯ + ( n + 2 + 4 + 2 n ) P x c k ) ] + ( 1 2 1 4 + 2 n ) ( 10 n 6 ) 3 36 n ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) k = 1 n ( n + 1 ) [ f ( x ¯ ( n + 2 + 4 + 2 n ) P x d k ) + f ( x ¯ + ( n + 2 + 4 + 2 n ) P x d k ) ]

3.3. Steps of the IST-7thSSRCKF

Based on the basic framework of spherical simplex-radial CKF, the proposed IST-7thSSRCKF algorithm is constructed by incorporating the strong tracking filter (STF) with an improved fading factor into the time update equation of seventh-degree SSRCKF. Hence, both the robustness against the process uncertainty and high accuracy can be achieved. The proposed IST-7thSSRCKF involves four steps for an iteration: (1) The state prediction; (2) Process uncertainty identification and calculation of the improved fading factor; (3) Covariance prediction using the improved strong tracking technique; (4) The measurement update.
Initialization: initialize state vector x ^ 0 | 0 and state covariance matrix P 0 | 0
  • Step 1. State prediction
Generate cubature points by employing the seventh-degree spherical simplex-radial cubature rule as follows:
P k 1 | k 1 = S k 1 | k 1 S k 1 | k 1 T
X i , k 1 | k 1 = S k 1 | k 1 ξ i + x ^ k 1 | k 1
ξ i = { x ¯ + ( n + 2 + 4 + 2 n ) P x a i i = 1 , , n + 1 , x ¯ ( n + 2 + 4 + 2 n ) P x a i n 1 i = n + 2 , , 2 n + 2 , x ¯ + ( n + 2 + 4 + 2 n ) P x b i 2 n 2 i = 2 n + 3 , , ( n 2 + 5 n + 4 ) / 2 , x ¯ ( n + 2 + 4 + 2 n ) P x b i ( n 2 + 5 n + 4 ) / 2 i = ( n 2 + 5 n + 6 ) / 2 , , n 2 + 3 n + 2 , x ¯ + ( n + 2 + 4 + 2 n ) P x c i n 2 3 n 2 i = n 2 + 3 n + 3 , , ( n 3 + 6 n 2 + 17 n + 12 ) / 6 , x ¯ ( n + 2 + 4 + 2 n ) P x c i ( n 3 + 6 n 2 + 17 n + 12 ) / 6 i = ( n 3 + 5 n 2 + 18 n + 18 ) / 6 , , ( n 3 + 3 n 2 + 8 n + 6 ) / 3 , x ¯ + ( n + 2 + 4 + 2 n ) P x d i ( n 3 + 3 n 2 + 8 n + 6 ) / 3 i = ( n 3 + 3 n 2 + 8 n + 9 ) / 3 , , ( n 3 + 6 n 2 + 11 n + 6 ) / 3 , x ¯ ( n + 2 + 4 + 2 n ) P x d i ( n 3 + 6 n 2 + 11 n + 6 ) / 3 i = ( n 3 + 6 n 2 + 11 n + 9 ) / 3 , , , ( n 3 + 9 n 2 + 14 n + 6 ) / 3 .
where n denotes the dimension of state vector x k ; the points sets { a i } , { b i } , { c i } and { d i } can be obtained by Equations (47)–(50).
Propagate the cubature points through the nonlinear state equation, the predicted state x ^ k | k 1 and predicted error covariance without the fading factor P k | k 1 l can be given by:
χ i , k | k 1 * = f ( χ i , k 1 )
x ^ k | k 1 = i = 1 ( n 3 + 9 n 2 + 14 n + 6 ) / 3 ω i χ i , k | k 1 *
P k | k 1 l = i = 1 ( n 3 + 9 n 2 + 14 n + 6 ) / 3 ω i ( χ i , k | k 1 * x ^ k | k 1 ) ( χ i , k | k 1 * x ^ k | k 1 ) T + Q k 1
where X i , k | k 1 denote the propagated cubature points, Q k 1 denotes the covariance matrix of process noise, ω i represent the corresponding weights of the cubature points and can be obtained by:
ω i = { ( 1 2 1 4 + 2 n ) n 2 ( 9 n 2 793 n + 1800 ) 36 ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) i = 1 , , 2 n + 2 , ( 1 2 1 4 + 2 n ) 4 ( n 1 ) 3 ( 4 n ) n ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) i = 2 n + 3 , , n 2 + 3 n + 2 , ( 1 2 1 4 + 2 n ) 486 ( n 2 ) 3 36 n ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) i = n 2 + 3 n + 3 , , ( n 3 + 3 n 2 + 8 n + 6 ) / 3 , ( 1 2 1 4 + 2 n ) ( 10 n 6 ) 3 36 n ( n + 1 ) 3 ( n + 2 ) ( n + 4 ) i = ( n 3 + 3 n 2 + 8 n + 12 ) / 6 , , ( n 3 + 9 n 2 + 14 n + 6 ) / 3 .
  • Step 2. Process uncertainty identification and calculation of the improved fading factor
Calculate the innovation ε k according to Equation (29) and identify the process uncertainty by using the hypothesis testing method described in:
(a)
If ε k < χ α , m 2 , Λ k = d i a g { 1 , 1 , 1 , 1 , 1 } , the strong tracking cubature Kalman filter therefore reduces to the standard cubature Kalman filter.
(b)
If ε k χ α , m 2 , perform the improved strong tracking cubature Kalman filter. The improved multiple fading factor used in the IST-7thSSRCKF is calculated as follows:
Λ i , k = { j = 1 6 G k [ i , j ] N k [ i 3 , j ] j = 1 6 G k 2 [ i , j ] , i = 4 , 5 , , 9 1 ,   i = 1 , 2 , 3 , 10 , , 21
  • Step 3. Covariance prediction using the improved strong tracking technique
Then, by introducing the improved multiple fading factor according to Equation (64), the modified prediction covariance can be calculated as:
P k | k 1 = Λ i , k ( P k | k 1 l Q k 1 ) Λ i , k T + Q k 1
  • Step 4. Measurement update
By utilizing the predicted state x ^ k | k 1 and the modified predicted covariance P k | k 1 in (65), we can obtain the modified predicted measurement z ^ k | k 1 , the innovation covariance matrix P z z , k | k 1 and the cross covariance matrix P x z , k | k 1 as follows:
X i , k | k 1 = c h o l ( P k | k 1 ) ξ i + x ^ k | k 1
Z i , k | k 1 = h ( X i , k | k 1 )
z ^ k | k 1 = i = 1 ( n 3 + 9 n 2 + 14 n + 6 ) / 3 ω i Z i , k | k 1
P z z , k | k 1 = i = 1 ( n 3 + 9 n 2 + 14 n + 6 ) / 3 ω i ( Z i , k | k 1 z ^ k | k 1 ) ( Z i , k | k 1 z ^ k | k 1 ) T + R k
P x z , k | k 1 = i = 1 ( n 3 + 9 n 2 + 14 n + 6 ) / 3 ω i ( X i , k | k 1 x ^ k | k 1 ) ( Z i , k | k 1 z ^ k | k 1 ) T
Finally, the state estimate x ^ k | k and the corresponding covariance P k | k at time k are calculated by Equations (26)–(28).

4. Performance Evaluation

In this section, the simulation and experimental results are presented to validate the proposed method. First, a new nonlinear filtering model that is suitable for the low-cost GPS/INS integrated navigation system is designed. Second, simulations and car-mounted experiments are carried out to evaluate the performance of different filters.

4.1. Design of SINS/GPS Filtering Model

In this paper, a 21-dimension state vector is adopted in the GPS/INS loosely coupled navigation system, which is given by:
x = [ ψ δ v δ p b g b f δ b g δ b f ] T
where ψ = [ ψ N ψ U ψ E ] , δ v = [ δ V N δ V U δ V E ] and δ p = [ δ λ δ L δ h ] are the attitude error vector, velocity error vector and position error vector respectively. b g denotes the static bias and δ b g is the dynamic bias noise of gyroscope expressed in body frame. Similarity, b f and δ b f are the static and dynamic bias of three axis accelerometer under the body frame respectively. The nonlinear system error model of INS is defined by:
{ ψ ˙ = C w [ ( I C n p ) w i n n + δ w i n n C b p δ w i b b ] δ v ˙ = ( I C p n ) C b p f i b b + C b p δ f i b b + δ v n × ( 2 w i e n + w e n n ) + v n × ( 2 δ w i e n + δ w e n n ) δ λ ˙ = δ V E R N + h sec L + δ L V E secL R N + h tan L δ L ˙ = δ V N R N + h δ h ˙ = δ V U b ˙ g = 0 b ˙ f = 0 δ b ˙ g = 1 τ g δ b g + η g δ b ˙ f = 1 τ f δ b f + η f
where C w is the transformation matrix from angle rate to Euler angles, C n p is the rotation matrix from n-frame(ideal navigation frame) to p-frame(actual navigation frame), C b p is the rotation matrix from b-frame(body frame) to p-frame. f i b b denotes the specific force expressed in body frame, δ f i b b is the corresponding error vector. ω b a c represents the rotation velocity of a-frame with respected to b-frame expressed in c-frame and δ ω b a c are the corresponding error vector. ψ is the misalignment angle and ψ N ψ U ψ E are the component in north, up and east direction. V E , V N , V U are the velocity component in east, north and up direction, L , λ , h are the latitude, longitude, and height, and δ V E , δ V N , δ V U , δ L , δ λ , δ h are the corresponding error. R N is the normal radius. τ g and τ f represent the correlation time of 1st-order Markovian process for gyroscope and accelerometer. η g and η f are the zero-mean Gaussian white noise process.

4.2. Simulation for SINS/GPS Integration

In this section, the proposed method is compared with the CKF [30], STCKF [25] and STSSRCKF [26] in the GPS/INS integrated navigation application. This application has been widely used as a benchmark to validate the performance of nonlinear filter due to its practical application value. The nonlinear model described in (72) will be employed as the dynamic process model, which can be given by:
x k = f ( x k 1 ) + v k 1
where f ( . ) is the nonlinear dynamic model function in Equation (72), v k 1 is Gaussian noise with zero mean and covariance Q k . Q k is determined from the inertial sensor’s stochastic error and set as Q k = diag([0.2 deg / h , 0.2 deg / h , 0.2 deg / h , 0.003 m / s 2 / H z , 0.003 m / s 2 / H z , 0.003 m / s 2 / H z , 0.2 deg / h , 0.2 deg / h , 0.2 deg / h , 0.003 m / s 2 / H z , 0.003 m / s 2 / H z , 0.003 m / s 2 / H z ])2. In order to show how the filters perform under process uncertainty, the process noise are generated according to [36]:
v k { N ( 0 ,   Q k ) w . p . 0.8 N ( 0 ,   100 Q k ) w . p . 0.2
where w . p . denotes “with probability”, i.e., eighty percent of process noise are generated from Gaussian distribution with nominal covariance Q k , and twenty percent of process noise are drawn from Gaussian distribution with severely increased covariance.
The measurement model can be formulated as follows:
z k = H k x k + ω k 1
where z k = ( V g p s N V i n s N V g p s U V i n s U V g p s E V i n s E L g p s L i n s h g p s h i n s λ g p s λ i n s ) , the subscripts (GPS and ins) represent the velocity and geographical position obtained from GPS and INS, respectively. H k = [ 0 6 × 3 , I 6 × 6 , 0 6 × 12 ] is the observation matrix. ω k 1 is the white Gaussian measurement noise with zero mean and covariance R k = diag([0.05 m/s, 0.05 m/s, 0.05 m/s, 1 m, 3 m, 1 m])2, which is determined from the GPS’s stochastic measurement error.
The simulation scenario is described in Table 1. For the first 30 s, the vehicle accelerates from 0 m/s to 30 m/s at the constant acceleration of 1 m / s 2 , and then head up and decelerates from 30 m/s to 24 m/s at the constant acceleration of −1 m / s 2 . Finally, the vehicle keeps this constant velocity and conducts the 8-driving. The three dimensional trajectory of the vehicle is shown in Figure 1.
In this simulation, the sample rate of IMU and GPS are 1000 Hz and 10 Hz, respectively. The output rate of GPS/SINS integration navigation system is 10 Hz. The initial state vector and the associated covariance are set as x ^ 0 | 0 = 0 1 × 21 , and P0|0 = diag([0.005 rad, 0.015 rad, 0.005 rad, 0.5 m/s, 0.5 m/s, 0.5 m/s, 10 m, 30 m, 10 m, 0.2 deg/s, 0.2 deg/s, 0.2 deg/s, 10 mg, 10 mg, 10 mg, 0.2 deg/s, 0.2 deg/s, 0.2 deg/s, 10 mg, 10 mg, 10 mg])2, respectively.
Similar to [30], we conduct 50 independent Monte Carlo runs for fair comparison and use the root mean square error(RMSE) as the metrics to compare the performance of four filters. The RMSEs of the attitude and velocity obtained by these filters are shown in Figure 2, respectively.
As can be seen from Figure 2, both the RMSEs from the proposed IST-7thSSRCKF, STSSRCKF and STCKF are smaller than that of the existing CKF. It means that the strong tracking-based filter can suppress the process uncertainty effectively. Moreover, we can also see from Figure 2 that the proposed IST-7thSSRCKF using the seventh-degree spherical simplex-radial cubature rule outperforms the STCKF and STSSRCKF as expected, which means that the estimation accuracy of the proposed IST-7thSSRCKF is further improved by the new Gaussian integral rule. Thus, the proposed method is more robust against process uncertainty for the GPS/INS integration navigation system as compared with the existing filters.

4.3. Car-Mounted Experiment for SINS/GPS Integration

In order to verify the effectiveness and superiority of the proposed algorithm, a car-mounted experiment was carried out by the advanced navigation system research group of the North University of China (Taiyuan, China). As shown in Figure 3, the experimental platform is composed of a LCI-1 tactical grade inertial measurement unit (IMU) and a Propak satellite receiver, which uses Novatel synchronized position attitude navigation (SPAN) as the reference solution [37]. The attitude, velocity and position accuracy of the Novatel SPAN reference system are 0.01 deg, 0.05 m/s and 0.1 m, respectively. The sensors specification of our self-made MINS/GPS integration navigation system is listed in Table 2. The experiment is carried out in North University of China and the total test time is 152 s. In the experiment test, the car firstly moved smoothly from 0 s to 79 s and then maneuvered severely from 79 s to 152 s. The test trajectory is shown in Figure 4. The reference attitude, velocity and position provided by the high-precision integrated navigation system during the whole test are shown in Figure 5. The original output of the tri-axis gyroscope and accelerometer are shown in Figure 6. It can be clearly seen from Figure 6 that there are much stochastic process uncertainties during time interval 80–152 s.

4.4. Performance Comparison with Different Filtering Algorithms

In this section, the proposed filter is tested and compared with the existing CKF [30], strong tracking cubature Kalman filter (ST-CKF) [25] and strong tracing spherical simplex-radial cubature Kalman filter (STSSRCKF) [26] in the car-mounted experiment for the low-cost GPS/INS integration navigation. For all the four methods, the initial state and corresponding error covariance are chosen as x ^ 0 | 0 = [ 0 1 × 9 , 0.2 deg/s, 0.2 deg/s, 0.2 deg/s, 16 mg, 16 mg, 16 mg, 0.0018 deg/s, 0.0018 deg/s, 0.0018 deg/s, 0.1 mg, 0.1 mg, 0.1 mg] and P0|0 = diag([0.005 rad, 0.015 rad, 0.005 rad, 0.1 m/s, 0.1 m/s, 0.1 m/s, 5 m, 10 m, 5 m, 0.2 deg/s, 0.2 deg/s, 0.2 deg/s, 16 mg, 16 mg, 16 mg, 0.0018 deg/s, 0.0018 deg/s, 0.0018 deg/s, 0.1 mg, 0.1 mg, 0.1 mg])2, respectively. The iSnitial value of the nominal process noise covariance matrix and measurement noise covariance matrix are respectively set as Q k = diag([0.3 deg / h , 0.3 deg / h , 0.3 deg / h , 0.3 m / s 2 / H z , 0.3 m / s 2 / H z , 0.3 m / s 2 / H z , 0.3 deg / h , 0.3 deg / h , 0.3 deg / h , 0.3 m / s 2 / H z , 0.3 m / s 2 / H z , 0.3 m / s 2 / H z ])2 and R k = diag([0.1 m/s, 0.1 m/s, 0.1 m/s, 5 m, 10 m, 5 m])2. In the ST-CKF, the forgetting factor ρ and softening factor β are empirically selected as 0.95 and 4.5 for optimal purpose. In the proposed improved strong tracking seventh-degree cubature Kalman filter (IST-7thSSRCKF), the parameters are set as: α = 0.05 , χ α , m 2 = 12.592 . The proposed algorithm and existing algorithms are all carried out by Matlab 2016 on a computer with an Intel core i5-3320 CPU at 2.60 GHz, 4 GB memory, and Windows 7 operating system. For a fair comparison, we conduct M = 20 independent runs. All the filters use the same parameter in each run.
The velocity and attitude error obtained from the four filtering algorithms when M = 10 are respectively shown in Figure 7 and Figure 8. Owing to the fact that the position error estimation is mainly determined by the measurement noise matrix R and is not sensitive to the time-varying process noise, it is unable to reflect the true performance of the four different filters in the maneuvering stage. For this reason, we do not present the result of the position error here. It can be seen from Figure 7 and Figure 8 that, in the smooth moving stage during intervals 0–79 s, the proposed method performs quite similar with the existing filters. The results are not surprising, this is because that in this case the Gaussian assumption for the statistical value of process noise can be satisfied. While for the maneuvering stage in the intervals of 80–152 s, the proposed IST-7thSSRCKF and the existing strong tracking based filters such as ST-CKF and ST-SSRCKF, exhibit smaller navigation error than the traditional CKF. This is because that in the maneuvering stage the Gaussian assumption of process noise is violated and the statistical value of process noise is changed due to the influence of the vehicle’s severe maneuvers, which results in a much degraded performance of CKF. In contrast, the severe vehicle maneuverability can be tracked appropriately by the ST-CKF, ST-SSRCKF and IST-7thSSRCKF as shown in Figure 7 and Figure 8. The performance improvement is mainly due to the strong tracking technique used in these filters, since the time-varying process uncertainty can be better suppressed by introducing the time-variant suboptimal fading factor matrix to the state prediction covariance matrix to maintain the residual sequence orthogonal. Furthermore, we can also see that the IST-7thSSRCKF outperforms the ST-CKF and ST-SSRCKF, which indicates that the 7th-degree SSR rule is more accurate in computing the multi-dimensional integrals than the SR and SSR rules. However, the accuracy improvement is at the cost of slightly larger computational burden as shown in Table 3. The reason is that the computation time is approximately proportional to the number of points used in these filters. The number of points required in the IST-7thSSRCKF is 4510, which is greater than that of CKF, STCKF and ST-SSRCKF, which are being 42, 42 and 44, respectively. In this respect, the implementation time of the proposed method is almost 100 times of CKF, STCKF and ST-SSRCKF. Actually, for the GPS/INS integrated navigation system with high performance navigation computer, the large computational burden can be afforded. Thereby, the proposed IST-7thSSRCKF is a better choice to deal with process uncertainty in terms of estimation accuracy when comparing with the investigated other three methods.
In order to compare the overall performance of the four methods and intuitively illustrate the effectiveness of the proposed filter for low-cost GPS/INS integration, the root mean square error (RMSE) of the velocity and attitude that averaged across all time instances is used as a comparison metric, which are defined as follows:
R M S E = 1 N k = 1 N ( x k x ^ k )
where N is the total sample number, x k and x ^ k denote the true value and the estimated value of velocity and attitude at time k, respectively. The RMSEs of the four methods in estimating the velocity and attitude of the low-cost GPS/INS integrated navigation system when M = 10 are listed in Table 4. The RMSEs of attitude error and velocity error when M = 1:20 are shown in Figure 9, respectively. It can be clearly seen from Table 4 and Figure 9 that the CKF exhibits large error in the maneuvering stage than the STCKF, STSSRCKF and IST-7thSSRCKF.
The STCKF and STSSRCKF maintain a good performance but are less accurate than the IST-7thSSRCKF. Therefore, taking both the estimation accuracy and computational efficiency into account, we can conclude that the proposed IST-7thSSRCKF can achieve high estimation accuracy with slightly higher computational complexity, and has better robustness for the suppression of process uncertainties.

5. Conclusions

In this paper, an improved strong tracking cubature Kalman filter is proposed to suppress the process uncertainty induced by the severe maneuvering for the low-cost GPS/INS integrated navigation systems. Based on the improved strong tracking technique, the process uncertainty can be detected and suppressed by modifying the prior state estimate covariance online according to the change in vehicle dynamics. Moreover, the seventh-degree spherical simplex-radial cubature rule, which can be used to compute the more exact Gaussian type integral, is integrated into the strong tracking cubature Kalman filter to further enhance the estimation accuracy. Thus, the proposed IST-7thSSRCKF possesses the advantage of 7thSSRCKF’s high accuracy and strong tracking filter’s strong robustness against process uncertainties. The car-mounted experiments are utilized to demonstrate that the proposed IST-7thSSRCKF can achieve high estimation accuracy and has better robustness for the suppression of process uncertainties.

Author Contributions

All authors were involved in the study in this manuscript. K.F. proposed the idea and wrote the paper; J.L. (Jie Li) and J.L. (Jun Liu) provided the theory analysis and guidance; X.Z. (Xiaoming Zhang) conceived and deigned the experiments. Y.Y. and X.Z. (Xi Zhang) carried out the experiments. C.S. and H.C. wrote the code of the proposed algorithm and analyzed the experimental results. All authors reviewed and approved the manuscript.

Funding

This work was supported in part by the National Natural Science Funds for Distinguished Young Scholars (51225504), the National Natural Science Foundation of China (51575500, 51705477), the Fund for Shanxi ‘1331 project’ Key Subject Construction and the Foundation for Middle-Aged and Young Talents in Higher Education Institutions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wu, Z.; Yao, M.; Ma, H.; Jia, W.; Tian, F. Low-cost antenna attitude estimation by fusing inertial sensing and two-antenna GPS for vehicle-mounted satcom-on-the-move. IEEE Trans. Veh. Technol. 2013, 62, 1084–1096. [Google Scholar] [CrossRef]
  2. Sun, D.; Petovello, M.G.; Cannon, M.E. Ultratight GPS/reduced-imu integration for land vehicle navigation. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 1781–1791. [Google Scholar] [CrossRef]
  3. Crassidis, J.L. Sigma-point Kalman filtering for integrated GPS and inertial navigation. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 750–756. [Google Scholar] [CrossRef]
  4. Zhao, L.; Qiu, H.; Feng, Y. Analysis of a robust Kalman filter in loosely coupled GPS/INS navigation system. Measurement 2016, 80, 138–147. [Google Scholar] [CrossRef]
  5. Falco, G.; Pini, M.; Marucco, G. Loose and tight GNSS/INS integrations: Comparison of performance assessed in real urban scenarios. Sensors 2017, 17, 255. [Google Scholar] [CrossRef] [PubMed]
  6. Oh, S.H.; Hwang, D.-H. Low-cost and high performance ultra-tightly coupled GPS/INS integrated navigation method. Adv. Space Res. 2017, 60, 2691–2706. [Google Scholar] [CrossRef]
  7. Dan, S. Optimal State Estimation: Kalman, H Infinity, and Nonlinear Approaches; Wiley-Interscience: Hoboken, NJ, USA, 2006. [Google Scholar]
  8. Chang, L.; Li, K.; Hu, B. Huber’s m-estimation-based process uncertainty robust filter for integrated INS/GPS. IEEE Sens. J. 2015, 15, 3367–3374. [Google Scholar] [CrossRef]
  9. Zhong, M.; Guo, J.; Cao, Q. On designing PMI Kalman filter for INS/GPS integrated systems with unknown sensor errors. IEEE Sens. J. 2015, 15, 535–544. [Google Scholar] [CrossRef]
  10. Chorin, A.J.; Tu, X. A tutorial on particle filters for online nonlinear/nongaussian bayesia tracking. Esaim Math. Model. Numer. Anal. 2012, 46, 535–543. [Google Scholar] [CrossRef]
  11. Sarkka, S.; Nummenmaa, A. Recursive noise adaptive Kalman filtering by variational bayesian approximations. IEEE Trans. Autom. Control 2009, 54, 596–600. [Google Scholar] [CrossRef]
  12. Dong, P.; Jing, Z.; Leung, H.; Shen, K. Variational bayesian adaptive cubature information filter based on wishart distribution. IEEE Trans. Autom. Control 2017, 62, 6051–6057. [Google Scholar] [CrossRef]
  13. Li, K.; Chang, L.; Hu, B. A variational bayesian-based unscented Kalman filter with both adaptivity and robustness. IEEE Sens. J. 2016, 16, 6966–6976. [Google Scholar] [CrossRef]
  14. Wang, S.-Y.; Yin, C.; Duan, S.-K.; Wang, L.-D. A modified variational bayesian noise adaptive Kalman filter. Circuits Syst. Signal Process. 2017, 36, 4260–4277. [Google Scholar] [CrossRef]
  15. Chen, B.; Principe, J.C. Maximum correntropy estimation is a smoothed map estimation. IEEE Signal Process. Lett. 2012, 19, 491–494. [Google Scholar] [CrossRef]
  16. Chen, B.; Liu, X.; Zhao, H.; Principe, J.C. Maximum correntropy Kalman filter. Automatica 2017, 76, 70–77. [Google Scholar] [CrossRef] [Green Version]
  17. Kulikova, M.V. Square-root algorithms for maximum correntropy estimation of linear discrete-time systems in presence of non-gaussian noise. Syst. Control Lett. 2017, 108, 8–15. [Google Scholar] [CrossRef]
  18. Liu, X.; Chen, B.; Xu, B.; Wu, Z.; Honeine, P. Maximum correntropy unscented filter. Int. J. Syst. Sci. 2017, 48, 1607–1615. [Google Scholar] [CrossRef]
  19. Wang, G.; Li, N.; Zhang, Y. Maximum correntropy unscented Kalman and information filters for non-gaussian measurement noise. J. Frankl. Inst. 2017, 354, 8659–8677. [Google Scholar] [CrossRef]
  20. Huang, Y.; Zhang, Y. A new process uncertainty robust student’s t based Kalman filter for SINS/GPS integration. IEEE Access 2017, 5, 14391–14404. [Google Scholar] [CrossRef]
  21. Hide, C.; Moore, T.; Smith, M. Adaptive Kalman filtering for low-cost INS/GPS. J. Navig. 2003, 56, 143–152. [Google Scholar] [CrossRef]
  22. Karasalo, M.; Hu, X. An optimization approach to adaptive Kalman filtering. In Proceedings of the 48th IEEE Conference on Decision and Control and 28th Chinese Control Conference, Shanghai, China, 16–18 December 2011; pp. 2333–2338. [Google Scholar]
  23. Bilik, I.; Tabrikian, J. Mmse-based filtering in presence of non-gaussian system and measurement noise. IEEE Trans. Aerosp. Electron. Syst. 2010, 46, 1153–1170. [Google Scholar] [CrossRef]
  24. Hu, G.; Gao, S.; Zhong, Y.; Gao, B.; Aleksandar, S. Modified strong tracking unscented Kalman filter for nonlinear state estimation with process model uncertainty. Int. J. Adapt. Control Signal Process. 2015, 29, 1561–1577. [Google Scholar] [CrossRef]
  25. Huang, W.; Xie, H.; Shen, C.; Li, J. A robust strong tracking cubature Kalman filter for spacecraft attitude estimation with quaternion constraint. Acta Astronaut. 2016, 121, 153–163. [Google Scholar] [CrossRef]
  26. Liu, H.; Wu, W. Strong tracking spherical simplex-radial cubature Kalman filter for maneuvering target tracking. Sensors 2017, 17, 741. [Google Scholar] [CrossRef] [PubMed]
  27. Zhou, D.H.; Frank, P.M. Strong tracking filtering of nonlinear time-varying stochastic systems with coloured noise: Application to parameter estimation and empirical robustness analysis. Int. J. Control 1996, 65, 295–307. [Google Scholar] [CrossRef]
  28. Jwo, D.J.; Wang, S.H. Adaptive fuzzy strong tracking extended Kalman filtering for GPS navigation. IEEE Sens. J. 2007, 7, 778–789. [Google Scholar] [CrossRef]
  29. Lin, C.L.; Chang, Y.M.; Hung, C.C.; Tu, C.D.; Chuang, C.Y. Position estimation and smooth tracking with a fuzzy-logic-based adaptive strong tracking Kalman filter for capacitive touch panels. IEEE Trans. Ind. Electron. 2015, 62, 5097–5108. [Google Scholar] [CrossRef]
  30. Arasaratnam, I.; Haykin, S. Cubature Kalman filters. IEEE Trans. Autom. Control 2009, 54, 1254–1269. [Google Scholar] [CrossRef]
  31. Xiong, H.; Tang, J.; Xu, H.; Zhang, W.; Du, Z. A robust single GPS navigation and positioning algorithm based on strong tracking filtering. IEEE Sens. J. 2018, 18, 290–298. [Google Scholar] [CrossRef]
  32. Chang, G. Robust Kalman filtering based on mahalanobis distance as outlier judging criterion. J. Geod. 2014, 88, 391–401. [Google Scholar] [CrossRef]
  33. Jia, B.; Xin, M.; Cheng, Y. High-degree cubature Kalman filter. Automatica 2013, 49, 510–518. [Google Scholar] [CrossRef]
  34. Lu, J.; Darmofal, D.L. Higher-dimensional integration with gaussian weight for applications in probabilistic design. SIAM J. Sci. Comput. 2008, 26, 613–624. [Google Scholar] [CrossRef]
  35. Genz, A.; Monahan, J. Stochastic integration rules for infinite regions. SIAM J. Sci. Comput. 1998, 19, 426–439. [Google Scholar] [CrossRef]
  36. Huang, Y.; Zhang, Y.; Li, N.; Chambers, J. A robust gaussian approximate fixed-interval smoother for nonlinear systems with heavy-tailed process and measurement noises. IEEE Signal Process. Lett. 2016, 23, 468–472. [Google Scholar] [CrossRef]
  37. Feng, K.; Li, J.; Zhang, X.; Shen, C.; Bi, Y.; Zheng, T.; Liu, J. A new quaternion-based Kalman filter for real-time attitude estimation using the two-step geometrically-intuitive correction algorithm. Sensors 2017, 17, 2146. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The three dimensional trajectory of the vehicle.
Figure 1. The three dimensional trajectory of the vehicle.
Sensors 18 01919 g001
Figure 2. RMSEs of the attitude (left) and velocity (right).
Figure 2. RMSEs of the attitude (left) and velocity (right).
Sensors 18 01919 g002
Figure 3. Setup of experimental platform.
Figure 3. Setup of experimental platform.
Sensors 18 01919 g003
Figure 4. Car-mounted test trajectory.
Figure 4. Car-mounted test trajectory.
Sensors 18 01919 g004
Figure 5. The reference attitude (left) and velocity (right) during the whole test.
Figure 5. The reference attitude (left) and velocity (right) during the whole test.
Sensors 18 01919 g005
Figure 6. The original outputs of the tri-axis gyroscope (left) and accelerometer (right).
Figure 6. The original outputs of the tri-axis gyroscope (left) and accelerometer (right).
Sensors 18 01919 g006
Figure 7. The attitude error in smooth stage (left) and maneuvering stage (right) from different filters.
Figure 7. The attitude error in smooth stage (left) and maneuvering stage (right) from different filters.
Sensors 18 01919 g007
Figure 8. The velocity error in smooth stage (left) and maneuvering stage (right) from different filters.
Figure 8. The velocity error in smooth stage (left) and maneuvering stage (right) from different filters.
Sensors 18 01919 g008
Figure 9. The RMSEs of attitude error (left) and velocity error (right) when M = 1:20.
Figure 9. The RMSEs of attitude error (left) and velocity error (right) when M = 1:20.
Sensors 18 01919 g009
Table 1. Description of vehicle motion.
Table 1. Description of vehicle motion.
Time (s)Motion
0–30Accelerate
30–36Head up and decelerate
37–47Uniform
48–878-driving
Table 2. Specifications of the SINS.
Table 2. Specifications of the SINS.
QuantityGyroscopeAccelerometer
Range±300°/s±10 g
Bias12°/h5 mg
Random walk0.28 deg / h 90 ug / H z
Table 3. Number of points and computation complexity of different filter for each run.
Table 3. Number of points and computation complexity of different filter for each run.
FilersPoints NumberTime (s)
CKF420.011
ST-CKF420.016
ST-SSRCKF440.018
IST-7thSSRCKF45100.968
Table 4. RMSEs of different filters in the maneuvering stage when M = 10.
Table 4. RMSEs of different filters in the maneuvering stage when M = 10.
FilersCKFST-CKFST-SSRCKFIST-7thSSRCKF
Azimuth (deg)1.701.191.110.96
Pitch (deg)0.740.490.240.19
Roll (deg)0.300.220.210.18
North Velocity (m/s)0.390.220.170.12
Up Velocity (m/s)0.500.270.250.16
East Velocity (m/s)0.490.230.180.14

Share and Cite

MDPI and ACS Style

Feng, K.; Li, J.; Zhang, X.; Zhang, X.; Shen, C.; Cao, H.; Yang, Y.; Liu, J. An Improved Strong Tracking Cubature Kalman Filter for GPS/INS Integrated Navigation Systems. Sensors 2018, 18, 1919. https://doi.org/10.3390/s18061919

AMA Style

Feng K, Li J, Zhang X, Zhang X, Shen C, Cao H, Yang Y, Liu J. An Improved Strong Tracking Cubature Kalman Filter for GPS/INS Integrated Navigation Systems. Sensors. 2018; 18(6):1919. https://doi.org/10.3390/s18061919

Chicago/Turabian Style

Feng, Kaiqiang, Jie Li, Xi Zhang, Xiaoming Zhang, Chong Shen, Huiliang Cao, Yanyu Yang, and Jun Liu. 2018. "An Improved Strong Tracking Cubature Kalman Filter for GPS/INS Integrated Navigation Systems" Sensors 18, no. 6: 1919. https://doi.org/10.3390/s18061919

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