Next Article in Journal
Enhanced Humidity Sensitivity with Silicon Nanopillar Array by UV Light
Previous Article in Journal
Methionine-Capped Gold Nanoclusters as a Fluorescence-Enhanced Probe for Cadmium(II) Sensing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Improved 5th-Cubature Kalman Filter in Initial Strapdown Inertial Navigation System Alignment for Large Misalignment Angles

Key Laboratory of Micro-Inertial Instrument and Advanced Navigation Technology, Ministry of Education, School of Instrument Science and Engineering, Southeast University, Nanjing 210096, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(2), 659; https://doi.org/10.3390/s18020659
Submission received: 22 December 2017 / Revised: 20 February 2018 / Accepted: 21 February 2018 / Published: 23 February 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
In view of the fact the accuracy of the third-degree Cubature Kalman Filter (CKF) used for initial alignment under large misalignment angle conditions is insufficient, an improved fifth-degree CKF algorithm is proposed in this paper. In order to make full use of the innovation on filtering, the innovation covariance matrix is calculated recursively by an innovative sequence with an exponent fading factor. Then a new adaptive error covariance matrix scaling algorithm is proposed. The Singular Value Decomposition (SVD) method is used for improving the numerical stability of the fifth-degree CKF in this paper. In order to avoid the overshoot caused by excessive scaling of error covariance matrix during the convergence stage, the scaling scheme is terminated when the gradient of azimuth reaches the maximum. The experimental results show that the improved algorithm has better alignment accuracy with large misalignment angles than the traditional algorithm.

1. Introduction

Initial alignment is the basis for obtaining accurate misalignment angles of Strapdown Inertial Navigation Systems (SINS), and its accuracy and efficiency will directly affect the SINS performance [1,2]. In an application environment with large misalignment angles and moving bases, the linear initial alignment algorithm is greatly restricted [3], while nonlinear algorithms have been increasingly used for initial alignment [4,5].
In order to apply the linear Kalman filter to a nonlinear system, the Extended Kalman Filter (EKF) [6] truncates the Taylor series expansion of the nonlinear function. The accuracy of EKF is low for strongly nonlinear systems, and there may be a problem in that the Jacobi matrix cannot be solved [7]. An Unscented Transformation (UT) is used in the Unscented Kalman Filter (UKF) [8] to replace the local linearization for obtaining the first two moments of the states. Although the precision of UKF is higher than EKF, there is a “dimension disaster” problem [9] in high dimensional applications of UKF, and parameter adjustment is needed when using UKF. Spherical Radial Cubature Rule (SRCR) is used in the third-degree CKF (hereinafter referred to as CKF3) to approximate the nonlinear integral and obtain the two moments of the states, which can achieve third-order accuracy [10]. CKF has a better numerical stability than UKF, and has been widely applied [11,12]. In order to further improve the accuracy of the CKF3, fifth-degree CKFs (hereinafter referred to as CKF5) [13,14] have been applied gradually [15]. CKF5 has problems of numerical stability and big calculation [16], and in this sense, CKF5 is more suitable for low dimensional systems. The adaptive Kalman filter is an effective method to restrain the divergence and improve the accuracy of filtering. The main idea of the adaptive Kalman filter is to scale the part of error covariance matrix or system noise matrix by calculating the ratio of matrix ranks of actual updating and filter updating [17,18]. The traditional method for estimating the updating covariance is the open window method, which uses the mean value of historical data and not takes into account the different weights of new and old updates. The accuracy and real-time performance of the traditional adaptive Kalman filter should be improved. Combining EKF and adaptive filtering [19], the update covariance matrix is directly introduced into the calculation of gain matrix, which improves the utilization rate of updates. However, EKF cannot effectively meet the accuracy requirement of the nonlinear initial alignment with large misalignment angles. Adaptive filtering is introduced into Unscented Particle Filter (UPF) [20], and the efficiency of convergence can be improved by partially scaling the error covariance matrix. Nevertheless, the precision of in-motion and nonlinear initial alignment with large misalignment angles should be further verified experimentally. The 5-th CKF is directly applied to the initial alignment [12], but experiments on initial alignment with large misalignment angles are lacking. The singular value decomposition can improve the robustness of the initial alignment algorithm [21], but there is no verification of actual data except for simulations.
In order to improve the efficiency and accuracy of the 5-th CKF in initial alignment situations, an entire scaling of the error covariance matrix method with the ratio of matrix ranks of actual updates and filtered updates is introduced in this paper. To solve the problem of oscillation or divergence of improved 5-th CKF caused by excessive scaling of the error covariance matrix, the scaling is terminated when the gradient of the azimuth error reaches the maximum. The improved algorithm needs to store the historical IMU data, GNSS data and the calculated azimuth misalignment angle data, and then recalculate it for alignment. To improve the stability of the 5-th CKF algorithm, the Singular Value Decomposition (SVD) algorithm is used instead of the Cholesky’s method in the matrix decomposition. The large misalignment angle model is simplified in this paper, which means that only the velocity errors and the misalignment angles are used as the states. The reduction of the states has the advantages of requiring less calculation and better calculation stability, which can reduce the chances of filter divergence.
In this paper, the nonlinear model of initial alignment with large misalignment angles is first analyzed and a simplified alignment model with 6-dimensional state variables is proposed. Secondly, the traditional CKF5 algorithm is briefly described in this paper, and then the improved algorithms based on fading factor for the updating covariance matrix, short-term scaling of the error covariance matrix and SVD decomposition are introduced. Finally, the effectiveness of the improved algorithms is verified by simulation experiments and the dynamic vehicle experiments.

2. Nonlinear Model of Initial Alignment with Large Misalignment Angle

The SINS updating algorithm includes attitude-updating, speed-updating and position-updating, and the corresponding differential equations are expressed with Equations (1) to (3) [22]:
C ˙ b n = C b n ( ω n b b × ) ,
v ˙ n = C b n f s f b ( 2 ω i e n + ω e n n ) × v n + g n ,
[ L ˙ λ ˙ h ˙ ] = [ 0 1 / R M h 0 sec L / R N h 0 0 0 0 1 ] v n .
where C b n is the direction cosine matrix of the attitude converting from the measurement reference frame (abbreviated as b frame) to the navigation reference frame (abbreviated as n frame), ω n b b = ω i b b ( ω i e b + ω e n b ) is the rotational angular rate of b frame relative to n frame, i is the inertial reference frame, ( ω n b b × ) is the antisymmetric matrix of ω n b b , the vector ν n = [ ν E   ν N   ν U ] T expresses the eastern velocity, northern velocity and vertical velocity, respectively, f s b b is the output of the accelerometer, L is the latitude of the location of carrier, λ is the longitude of the location of carrier, h is the height of the location of carrier, RMh is the radius of meridian of the location of carrier and RNh is the radius of prime vertical circle of the location of carrier. Assuming n′ and n are the calculation reference frame and ideal navigation reference frame, respectively, the misalignment (expressed as Euler angle) between the two navigation reference frames can be expressed as:
ϕ ˙ = C ω 1 ω n n n ,
ω n n n = ( I C n n ) ω ^ i n n + C n n δ ω i n n C b n ( ε b + w g b ) ,
where ω n n n is the relative rotational angular rate between the n reference frame and n′ reference frame, εb is the random constant drift of Fiber Optic Gyroscope (FOG) and w g b is the random drift of FOG. C ω 1 and C n n can be expressed as:
C ω 1 = [ cos ϕ y 0 sin ϕ y sin ϕ x sin ϕ y / cos ϕ x 1 sin ϕ x cos ϕ y / cos ϕ x sin ϕ y / cos ϕ x 0 cos ϕ y / cos ϕ x ]
C n n [ cos ϕ y cos ϕ z sin ϕ x sin ϕ y sin ϕ z cos ϕ y sin ϕ z + sin ϕ x sin ϕ y sin ϕ z cos ϕ x sin ϕ y cos ϕ x sin ϕ z cos ϕ x cos ϕ z sin ϕ x sin ϕ y cos ϕ z + sin ϕ x cos ϕ y sin ϕ z sin ϕ y sin ϕ z sin ϕ x cos ϕ y cos ϕ z cos ϕ x cos ϕ y ] .
Under the reference frame, Equation (2) can be written as:
v ^ ˙ n = C b n f ^ s f b ( 2 ω ^ i e n + ω ^ e n n ) × v ^ n + g ^ n .
Subtracting Equation (2) from Equation (8) and neglecting the calculation error of gravitational acceleration and higher-order terms, we obtain:
δ v ˙ n = [ I ( C n n ) T ] C b n f ^ s f b + ( C n n ) T C b n b ( 2 ω ^ i e n + ω ^ e n n ) × δ v n ,
where ν ^ n = ν n + δ ν n ,   f ^ s f b = f s f b + δ f s f b ,   b = δ f s f b , and we have:
ε ˙ b = 0 ,
˙ b = 0 .
The purpose of initial alignment is to quickly identify the misalignment angles. It is difficult to estimate the errors of inertial sensors accurately, and that is not the purpose of initial alignment, then the errors of εb and ∇b can be regarded as system noises. Therefore, the states used in initial alignment are misalignment angles φ and the velocity error δνn, which can reduce the dimensions of the filter, and thus reduce the computation and improve the numerical stability of CKF5. In order to simplify the filtering error equation, a SINS nonlinear error model which can reflect the characteristics of SINS errors has important practical value in the initial alignment [20]. When the SINS and GNSS are integrated for initial alignment under moving base conditions, GNSS can provide relatively accurate velocity and location information, then Equations (4) and (9) can be simplified as:
ϕ ˙ = C ω 1 [ ( I C n n ) ω ^ i n , g n ε n ] + C ω 1 C b n w g b ,
δ v ˙ n = [ I ( C n n ) T ] C b n f ^ s f b + ( C n n ) T n + ( C n n ) T C b n w a b ,
where w a b is the random drift of accelerometer, ω ^ i n , g n is the angular velocity of the Earth determined by the information from GNSS. With the difference between the velocity of SINS measured by GNSS and the calculated velocity working as the measurement of z, the nonlinear filtering equations of initial alignment with large misalignment angles can be expressed as:
{ x ˙ = f ( x ) + g ( x ) w , z = H x + v ,
where x = [φ (δνn)T], w = [ ( w g b ) T ( w a b ) T ] , H = [03×3 I3×3], ν is the error vector of measurement and the function of f can be obtained from Equation (12). The model of initial alignment (Equations (12) and (13)) takes into account of the moving base situation, so the estimation of the initial alignment angles will theoretically not be affected by the line movement and angle movement of SINS.

3. CKF5 and Improved CKF5

3.1. CKF5 Algorithm

For the integral of the global nonlinear function I(f) and variable x = ry in the spherical-radial coordinate system, I(f) can be written as:
I ( f ) = 0 U n f ( r y ) r n 1 e r 2 d σ ( y ) d r ,
where δTδ = 1, r ∈ [0, ∞), Un is n-dimensional unit hypersphere and σ(g) is the element of Un. I(f) can be divided into two parts, spherical integral S(r) and radial integral R:
{ S ( r ) = U n f ( r y ) d σ ( y ) , R = 0 S ( r ) r n 1 e r 2 d r .
The numerical integration is approximated with the summation of weighted sigma points in the algorithm of CKF, then Equation (16) can be expressed as:
{ S ( r ) j = 1 N s ω s , j f ( r y j ) , R i = 1 N r ω r , i f ( r i ) .
Then:
I ( f ) i = 1 N r j = 1 N s ω r , i ω s , j f ( r i y j ) .
The m-order approximation of S(r) can be written as:
S = { 2 i = 1 n Γ [ ( d i + 1 ) / 2 ] Γ [ ( | d | + n ) / 2 ] ,   { d i }   are   even   numbers , 0 ,   at   least   one   odd   number   in   { d i } ,
where d = [p1,p2,…,pn], | d = i = 1 n d i | , | d | m , and Γ(g) is the gamma function. Note that when m = 5, the value of S is calculated with f 1 ( y ) = 1 ,   f 2 ( y 1 ) = y 1 2 ,   f 3 ( y 1 ) = y 1 4 and f 4 ( y 1 , y 2 ) = y 1 2 y 2 2 . The points set on the Un is defined as:
{ [ e ] n = { ± { e i } i = 1 n } , [ s ] n = { ± { s i } i = 1 n ( n 1 ) } ,
where { e i } i = 1 n and { s i } i = 1 n ( n 1 ) are defined as:
{ { e i } i = 1 n = { e 1 , e 2 , , e n } , { s i } i = 1 n ( n 1 ) = { ( e k ± e l ) / 2 :   k , l = 1 , 2 , , n .   and   k < l } ,
where e i is the ith column of identity matrix I n × n . Then the Equation (19) can be expressed as:
S ( r ) ω 1 i = 1 2 n f ( r [ e ] i ) + ω 2 i = 1 2 n ( n + 1 ) f ( r [ s ] i ) ,
where ω 1 = ( 4 n ) A n 2 n ( n + 2 ) , ω 2 = A n n ( n + 2 ) and A n = { [ a ] n , [ s ] n } . For the symmetry, S 1 ( r ) = 1 ,   S 2 ( r ) = r 2 and S 3 ( r ) = r 4 are needed to calculated in the solution of 5-order integral of S ( r ) . Then
{ ω r , 1 r 1 0 + ω r , 2 r 2 0 = 1 2 Γ ( n / 2 ) , ω r , 1 r 1 2 + ω r , 2 r 2 2 = 1 2 Γ ( n / 2 + 1 ) = n 4 Γ ( n / 2 ) , ω r , 1 r 1 4 + ω r , 2 r 2 4 = 1 2 Γ ( n / 2 + 2 ) = 1 2 ( 1 2 n + 1 ) ( n / 2 ) Γ ( n / 2 ) .
Let r 1 = 0 , then we have:
{ r 2 = 1 2 n + 1 , ω r , 1 = 1 ( n + 2 ) Γ ( n / 2 ) , ω r , 2 = n 2 ( n + 2 ) Γ ( n / 2 ) .
Then it follows from Equations (19) and (23) that the integrated result of Equation (18) can be obtained. The algorithm flow of CKF5 is shown as follows:
Prediction phase:
  • Let x k 1 ~ N ( x k 1 ;   x ^ k 1 | k 1 ,   P k 1 | k 1 ) be the estimated state at time k 1 and decompose P k 1 | k 1 with Cholesky’s method:
    P k 1 | k 1 = S k 1 | k 1 S k 1 | k 1 T .
  • Compute cubature points:
    { x 1 = S k 1 | k 1 0 + x ^ k 1 | k 1 , x 2 = S k 1 | k 1 n + 2 s j + + x ^ k 1 | k 1 , x 3 = S k 1 | k 1 n + 2 s j + + x ^ k 1 | k 1 , x 4 = S k 1 | k 1 n + 2 s j + x ^ k 1 | k 1 , x 5 = S k 1 | k 1 n + 2 s j + x ^ k 1 | k 1 , x 6 = S k 1 | k 1 n + 2 e j + x ^ k 1 | k 1 , x 7 = S k 1 | k 1 n + 2 e j + x ^ k 1 | k 1 ,
    where { s j + } { 1 2 ( e k + e l ) :   k < l ,   k , l = 1 , 2 , , n } and { s j } { 1 2 ( e k e l ) :   k < l ,   k , l = 1 , 2 , , n } .
  • Compute the predicted mean x ^ k | k 1 and predicted covariance P k | k 1 :
    x ^ k | k 1 = w 1 f ( x 1 ) + w 2 j = 1 n ( n 1 ) / 2 ( f ( x 2 ) + f ( x 3 ) ) + w 2 j = 1 n ( n 1 ) / 2 ( f ( x 4 ) + f ( x 5 ) ) + w 3 j = 1 n ( f ( x 6 ) + f ( x 7 ) ) ,
    P k | k 1 = w 1 f ( x 1 ) f ( x 1 ) T + w 2 j = 1 n ( n 1 ) / 2 ( f ( x 2 ) f ( x 2 ) T + f ( x 3 ) f ( x 3 ) T ) + w 2 j = 1 n ( n 1 ) / 2 ( f ( x 4 ) f ( x 4 ) T + f ( x 5 ) f ( x 5 ) T ) + w 3 j = 1 n ( f ( x 6 ) f ( x 6 ) T + f ( x 7 ) f ( x 7 ) T ) x ^ k | k 1 x ^ k | k 1 T + Q k 1 ,
    where w 1 = 2 n + 2 , w 2 = 1 ( n + 2 ) 2 , w 3 = 4 n 2 ( n + 2 ) 2 and Q k 1 is system noise covariance matrix.
Update phase:
4.
Decompose P k | k 1 with Cholesky’s method:
P k | k 1 = S k | k 1 S k | k 1 T .
5.
Compute cubature points:
{ x 1 = S k | k 1 0 + x ^ k | k 1 , x 2 = S k | k 1 n + 2 s j + + x ^ k | k 1 , x 3 = S k | k 1 n + 2 s j + + x ^ k | k 1 , x 4 = S k | k 1 n + 2 s j + x ^ k | k 1 , x 5 = S k | k 1 n + 2 s j + x ^ k | k 1 , x 6 = S k | k 1 n + 2 e j + x ^ k | k 1 , x 7 = S k | k 1 n + 2 e j + x ^ k | k 1 .
6.
Compute the estimated measurement z ^ k | k 1 , the gain matrix K k and the covariance matrixes of P k | k 1 z z , P k | k 1 x z :
z ^ k | k 1 = w 1 h ( x 1 ) + w 2 j = 1 n ( n 1 ) / 2 ( h ( x 2 ) + h ( x 3 ) ) + w 2 j = 1 n ( n 1 ) / 2 ( h ( x 4 ) + h ( x 5 ) ) + w 3 j = 1 n ( h ( x 6 ) + h ( x 7 ) ) ,
P k | k 1 z z = w 1 h ( x 1 ) h ( x 1 ) T + w 2 j = 1 n ( n 1 ) / 2 ( h ( x 2 ) h ( x 2 ) T + h ( x 3 ) h ( x 3 ) T ) + w 2 j = 1 n ( n 1 ) / 2 ( h ( x 4 ) h ( x 4 ) T + h ( x 5 ) h ( x 5 ) T ) + w 3 j = 1 n ( h ( x 6 ) h ( x 6 ) T + h ( x 7 ) h ( x 7 ) T ) z ^ k | k 1 z ^ k | k 1 T + R k ,
P k | k 1 x z = w 1 x 1 h ( x 1 ) T + w 2 j = 1 n ( n 1 ) / 2 ( x 2 h ( x 2 ) T + x 3 h ( x 3 ) T ) + w 2 j = 1 n ( n 1 ) / 2 ( x 4 h ( x 4 ) T + x 5 h ( x 5 ) T ) + w 3 j = 1 n ( x 6 h ( x 6 ) T + x 7 h ( x 7 ) T ) x ^ k | k 1 z ^ k | k 1 T ,
K k = P k | k 1 x z P k | k 1 z z 1 .
7.
Update x ^ k | k and P k | k :
P k | k = P k | k 1 K k P k | k 1 z z K k T ,
x ^ k | k = x ^ k | k 1 + K k ( z k z ^ k | k 1 ) ,
where z k is the measurement at time k .

3.2. Improved CKF5 Algorithm

The filtering precision of CKF5 can reach the fourth order, and the accuracy of CKF3 can reach the second order. The amount of sigma points of CKF5 is 2 n 2 + 1 , and that of CKF3 is 2 n , so the cost of computing of CKF3 is less than that of CKF5. In order to estimate the numerical stability of nonlinear filter, the stability factor can be defined as [23]:
I = i = 1 N | ω i | ,
where ω i is the weight of sigma point and N is the amount of sigma points. If I = 1 , the numerical integration has absolute stability theoretically, on the contrary, non-absolute stability for I 1 . When I 1 , the cumulative truncation error may leads to a non-positive error covariance matrix, even further to the error of calculation or divergence. For CKF3, I = 1 and for CKF5, I 1 , which means that CKF3 has better numerical stability than CKF5. For CKF5, we have:
I = | 2 n + 2 | + | n ( n 1 ) ( n + 2 ) 2 | + | n ( 4 n ) 2 ( n + 2 ) 2 | = 2 n + 2 + n ( n 1 ) ( n + 2 ) 2 + n ( n 4 ) 2 ( n + 2 ) 2 ,   n > 4 , = 3 n 2 2 n + 8 2 ( n + 2 ) 2 ,   n > 4 .
According to Equation (14), n = 6 and I = 13 / 16 , the robustness of CKF is poor in this situation. The non-positive covariance matrix caused by strong nonlinear or data errors may bring many problems in Kalman filters, such as filtering error or divergence, which also can lead to the error of Cholesky decomposition in CKF. In order to improve the numerical stability of CKF5, and considering that the singular value decomposition (SVD) has better numerical robustness than Cholesky [24], SVD is introduced into CKF5.
For s linear optimal Kalman filter, when the gain matrix K k is optimal, the innovation sequence k should be orthogonal. This conclusion can be extended to nonlinear filtering. The adaptive Kalman (hereinafter referred to as AKF) scales the error covariance matrix or system noise matrix by calculating the ratio of matrix ranks of actual innovation and filtered innovation [17,18]. Essentially, the traditional AKF only partly scale the error covariance matrix in the next filtering period, and the utilization rate of innovation is low. For the innovation sequence, the contribution of the new innovation to the state estimation is large than the old one. Therefore, ICKF5 uses the exponent fading factor to weigh the innovation sequence, and then scale the entire error covariance matrix in the next filtering period. Assuming the current filter period is k and the weight set of innovation sequence is { α i } , we have:
i = 1 k α i = 1 ,   α i 1 = α i ξ ,
where ξ is fading factor and ξ ( 0 ,   1 ) , empirically, ξ ( 0.6 ,   1 ) . To illustrate the impact of ξ on α i , ξ was taken as 0.7, 0.8, and 0.9, respectively, for the convenience of observation, the innovations are weighted by the exponent fading factor from k = 100 . The distribution of α i is shown in Figure 1.
Figure 1 shows that if the ξ is big, the weight of the current innovation is small, correspondingly, the weights of old innovations are big, which means that the old innovation has bigger influence than the new innovation on calculating the covariance matrix of innovation sequence. The opposite conclusion can be obtained if ξ is small. A recursive formula for the covariance matrix of the innovation sequence k is given below. When the forgetting factor ξ is determined, then:
ρ k = i = 1 k 1 ξ i .
It follows from Equations (39) and (40) that:
α i = ξ k i ρ k ,   i = 1 , 2 , , k .
The open window method is used in traditional AKF, the covariance matrix P ^ , k of innovation k at time k can be expressed as:
P ^ , k = 1 k 1 i = 1 k i i T .
The weight of 1 / ( k 1 ) in Equation (42) can be replaced by the weights of α i in Equation (41), then a recursive formula between P ^ , k and P ^ , k 1 can be expressed as:
P ^ , k = i = 1 k ξ k i ρ k Z i Z i T = Z k Z k T ρ k + i = 1 k 1 ξ k i Z i Z i T ρ k = Z k Z k T ρ k + ξ ρ k 1 ρ k ρ k 1 i = 1 k 1 ξ k i 1 Z i Z i T = Z k Z k T ρ k + ξ ρ k 1 ρ k ρ k 1 i = 1 k 1 ξ k i 1 Z i Z i T = Z k Z k T ρ k + ρ k 1 ρ k ( 1 ρ k 1 i = 1 k 1 ξ k i 1 Z i Z i T ) = Z k Z k T ρ k + ρ k 1 ρ k P ^ , k 1 .
In order to improve the feedback effect of innovation, we scale the error covariance matrix of the next filtering period β k times entirely. The increase of the error covariance matrix makes the filter gain K k increase and the error correction ability of the innovation in the next filtering period is strengthened. This method can reduce the dependence of the filter on the model, and then improve the filtering accuracy and reduce the risk of divergence. β k is defined as:
β k = M a x { t r a c e ( P ^ , k R k ) t r a c e ( H k P k | k 1 H k T ) ,   1 } ,
where t r a c e ( · ) is the operator of obtaining the rank of a matrix, R k is the measurement error covariance matrix at time k , H k is the measurement matrix at time k . The error covariance matrix of next filtering period is scaled as:
P k + 1 | k = β k P k + 1 | k .
The algorithm of ICKF5 is summarized as follows:
Prediction phase:
(1)
Input the initial values for filtering: x k 1 | k 1 ,   P k 1 | k 1 ,   Q k 1 ,   β k 1 = 1 ,   R k .
(2)
Compute cubature points x i ,   ( i = 1 , , 7 ) according to Equation (26).
(3)
Compute the predicted mean x ^ k | k 1 and predicted covariance P k | k 1 according to Equations (27) and (28), respectively.
(4)
Scaling the entire P k | k 1 : P k | k 1 = β k 1 P k | k 1 .
Update phase:
(5)
Decompose P k | k 1 with the SVD method.
(6)
Compute cubature points x i ,   ( i = 1 , , 7 ) according to Equation (30).
(7)
Compute the estimated measurement z ^ k | k 1 , the gain matrix K k and the covariance matrixes of P k | k 1 z z , P k | k 1 x z according to Equations (31)–(34), respectively.
(8)
Update x ^ k | k , P k | k and β k according to Equations (35), (36) and (44), respectively.
The scaling of the error covariance matrix by ICKF5 can effectively improve the utilization of innovation and speed up the convergence of initial alignment. ICKF5 is a kind of suboptimal algorithm, if the scaling value of β k is too large, the system will be in a state of overshoot, which may lead to oscillation or divergence. For example, when getting the initial alignment in the convergence stage, the gain K k and the error covariance matrix P k | k 1 should ideally remain unchanged or little changed. If the feedback for error covariance matrix is big at the stage of alignment convergence, it may cause the filtering convergence to be unstable or lead to the divergence. In order to reduce the possibility of divergence of initial alignment filtering and ensure the stationarity after convergence, we use a short-term feedback strategy that we stop the scaling of β k in ICKF5 before the filtering occurs oscillating or diverging, which can ensure the convergence of initial alignment. For convenience, we call the short-term feedback algorithm IICKF5. The first step of IICKF5 is that the ICKF5 is used to complete the initial alignment firstly, then the gradient sequence ψ i = diff ( ψ i ) of the set of azimuth misalignment angle { ψ i } outputted by ICKF5 can be obtained, where d i f f ( · ) is gradient operator. The second step, the time T max corresponding to the maximum value of the gradient is chosen as the time to terminate the feedback to the error covariance matrix P k | k 1 . T max can be expressed as:
T max = M a x { a b s { ψ } } .
To be brief, ICKF5 is used in the period of [ 1 , T max ) and CKF5 is used in [ T max , T e n d ] , where T e n d is the end time of initial alignment filtering. The flow chart of IICKF5 is shown in Figure 2.
Compared to ICKF5, IICKF5 needs to recalculate the stored IMU data and the misalignment angle set {xk} output by ICKF5. IICKF5 has a certain delay compared with ICKF5, and some storage space is needed.

3.3. Simulation and Its Analysis

In order to verify the effectiveness of ICKF5 and IICKF5, a simulation of initial alignment under static base conditions is performed. The simulation conditions are set as follows: the gyro constant drift is ε C = 0.01   ( / h ) , random drift is ε N = 0.001   ( / h ) ; the accelerometer zero bias is C = 50   μ g , random deviation is N = 5   ( μ g / h ) ; the local geographic latitude is 34 ° , longitude is 108 ° ; the misalignment angles are set as ϕ ( 0 ) = [ 20 ° ,   40 ° ,   170 ° ] and the simulation period is 800   s . The general parameters of the initial alignment algorithms are set with the same values. The initial value of X 0 is set as X 0 = [ 0 ; 0 ; 0 ; 0 ; 0 ; 0 ] . The matrix of Q is determined according to the random drift of IMU, and Q k = d i a g { ( 0.001 ) ;   ( 0.001 ) ;   ( 0.001 ) ;   ( 5 u g ) ;   ( 5 u g ) ;   ( 5 u g ) } 2 . The initial error covariance matrix is set as P 0 = d i a g ( [ [ 10 ; 10 ; 180 ] π / 180 ;   [ 1 ; 1 ; 1 ] ] ) 2 . Where d i a g ( · ) is the operator of obtaining the diagonal matrix with diagonal elements. For UKF, α = 1 e 3 , β = 2 and λ = 0 . EKF, UKF [25], CKF3, CKF5, ICKF5 and IICKF5 algorithms are applied to the nonlinear initial alignment. In addition to IICKF5, the calculated misalignment angles are as shown in Figure 3.
Figure 3 shows that the pitching misalignment angles converge to 20° and the rolling misalignment angles converge to −40° with the exception of EKF, and only the azimuth misalignment angles of CKF5 and ICKF5 converge to 170°.
The gradient values of the azimuth angle errors of ICKF5 are shown in Figure 4.
Figure 4 shows that there is an obvious mutation in the azimuth angle error at t = 4   s , which is the overshoot caused by ICKF5, so T max = 4   s . The estimated azimuth angles of CKF5, ICKF5 and IICKF5 are shown in Figure 5.
As seen in Figure 5, IICKF5 enters the convergence stage at about 400 s, while CKF5 and ICKF5 converge slower than IICKF5. Figure 3 also shows that the convergences of horizontal misalignment angles are faster than that of the azimuth misalignment angles, which is due to the fact that the observability of the horizontal misalignment angles is high under the velocity matching mode. There is a coupling relationship between the horizontal misalignment angles and the bias of accelerometer. When the variation of heading angle is small, the accuracies of the horizontal misalignment angles depend on the bias of the accelerometer. The simulation experiment in this paper is based on static base, so the variation of heading angle is very small, and the horizontal misalignment has a slow drift due to the bias of the accelerometer. From this point of view, the estimation period of the horizontal misalignment angles is not the longer the better. For the sake of accuracy, we select the estimations of horizontal misalignment angles at 50 s and select the estimation value of azimuth misalignment angle at the last moment (800 s) as the outputs of alignment, and the alignment errors are shown in Table 1.
Table 1 shows that the pitch and the azimuth misalignment errors of IICKF5 are the smallest, and the azimuth misalignment error is close to the limitd alignment precision of 0.05°, which meets the initial alignment requirement.

4. Dynamic Vehicle Experiment

The experimental equipment is shown in Figure 6. The GNSS real-time differential system consisted of two model ProPak6 receivers (NovAtel, Calgary, AB, Canada). The attitude reference system is a NovAtel SPAN system (NovAtel, Calgary, AB, Canada), which can reach centimeter level precision. An LN-200 IMU named from the Novatel Company is used for the initial alignment experiment. The bias of the accelerometer is 0.3 mg and its scale factor is 300 ppm. The constant drift, the random drift and the scale factor of gyro are ε C = 1   ( / h ) , ε N = 0.07 ( / h ) and 100 ppm, respectively. The output rates of the LN-200, SPAN system and GNSS system are 200 Hz, 100 Hz and 5 Hz, respectively.
The period of initial alignment is 700 s. In addition to IICKF5, the calculated azimuth angle errors are as shown in Figure 7.
For the azimuth angle errors, Figure 7 shows that EKF, UKF, CKF3 did not converge to the accurate azimuth angle, ICKF5 has faster convergence rate than CKF5, and ICKF5 has much more oscillations in the convergence process than CKF5 because of the overshoot. For the pitch angle and the roll angle errors, we can see that CKF5 and ICKF5 have greater oscillations than other algorithms, which is caused by the inaccuracy of model noises, such as the GPS’s measurement noise, IMU noise and so on. The nonlinear filter of high order precision amplifies the noises with high frequencies, which easily results in larger estimation errors if the estimated state values are small. When obtaining the gradient of the azimuth angles of ICKF5 in Figure 7c, the gradient values are as shown in Figure 8.
Figure 8 shows that there is a mutation in the azimuth angle error at t = 50   s , so T max = 50   s . The estimated misalignment angles of all the algorithms are shown in Figure 9.
From Figure 9 we can find that IICKF5 has a faster convergence rate than CKF5 and ICKF5 except for the pitching misalignment angle, and the azimuth misalignment angle of IICKF5 has better convergence stability than ICKF5. The alignment errors at the end of alignment (at 700 s) are shown in Table 2.
Table 2 shows that IICKF5 has better accuravy than CKF5 in estimating the rolling misalignment angle and azimuth misalignment angle. When the alignment is finished, the alignment errors of ICKF5 and IICKF5 did not converge to small values (e.g., 1°), so more time is needed for fine alignment. In order to further compare the convergence abilities of different algorithms, the mean and standard deviation of azimuth misalignment errors are calculated in the range of [100 s, 700 s], which can be regarded as a relatively complete convergence period for ICKF5 and IICKF5, and the results are shown in Table 3.
From Table 3 we can see that IICKF5 has the smallest errors of mean and standard deviation. In this sense, the improved algorithm is suitable for quickly identifying a large azimuth misalignment angle in a complex alignment environment. ICKF5 is a kind of real-time algorithm that represents an improvement on CKF5, which does not need to process historical IMU data and no storage space for the IMU is needed. IICKF5 needs to obtain the gradient value of the historical azimuth angle, so the historical IMU data need to be stored for the second initial alignment filtering. For the 700 s period of in-motion initial alignment, the length of IMU data is 70,000 (in order to synchronize with the reference system), and the storage space for IMU data is 2869 KB. The data stored are double and floating types with 64 bits. The update rate of GNSS is 5 Hz, so the storage space it needs is small, and the size is 10 KB. The capacity of the azimuth misalignment angle data is 6 KB. The total storage space needed for IICKF5 is 2885 KB, and the additional burden on the navigation computer is small.

5. Conclusions

It is not only convenient to calculate the innovation variance matrix in a recursive and exponent fading factor way, but also to effectively use the current and old updates. The speed and precision of convergence can be effectively improved by scaling error covariance matrix with the ratio of matrix ranks of actual updates and filter updates. In order to solve the problem of oscillation or divergence caused by excessive scaling of the error covariance matrix, the scaling scheme is terminated when the azimuth gradient reaches the maximum. The improved CKF can improve the precision of initial alignment better than the traditional CKF. The length of filtering time is also dependent on the experience, and the robustness of the improved algorithm is also worthy of further study.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (No. 51375087), the Scientific Research Foundation of Graduate School of Southeast University (No. YBJJ1737) and the Postgraduate Research & Practice Innovation Program of Jiangsu Province (No. KYCX17_0065).

Author Contributions

Wei Wang and Xiyuan Chen conceived and designed this study. Wei Wang performed the experiments and wrote the paper. Xiyuan Chen reviewed and edited the manuscript. All authors read and approved this manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dmitriyev, S.P.; Stepanov, O.A.; Shepel, S.V. Nonlinear filtering methods application in INS alignment. IEEE Trans. Aerosp. Electron. Syst. 1997, 33, 260–271. [Google Scholar] [CrossRef]
  2. Chang, L.B.; Hu, B.Q.; Li, A.; Qin, F.J. Strapdown inertial navigation system alignment based on marginalized unscented Kalman filter. IET Sci. Meas. Technol. 2013, 7, 128–138. [Google Scholar] [CrossRef]
  3. Li, H.; Pan, Q.; Wang, X.; Jiang, X.; Deng, L. Kalman Filter Design for Initial Precision Alignment of a Strapdown Inertial Navigation System on a Rocking Base. J. Navig. 2015, 68, 184–195. [Google Scholar] [CrossRef]
  4. Hao, Y.L.; Mu, H.W.; Jia, H.M. Application of ICDKF in initial alignment of large azimuth misalignment in SINS. Syst. Eng. Electron. 2013, 35, 152–155. [Google Scholar]
  5. Chatzi, E.N.; Smyth, A.W. The unscented Kalman filter and particle filter methods for nonlinear structural system identification with non-collocated heterogeneous sensing. Struct. Control Health Monit. 2009, 16, 99–123. [Google Scholar] [CrossRef]
  6. Cox, H. On the estimation of state variables and parameters for noisy dynamic systems. IEEE Trans. Autom. Control 1964, 9, 5–12. [Google Scholar] [CrossRef]
  7. Crassidis, J.L.; Markley, F.L. Unscented filtering for spacecraft attitude estimation. J. Guid. Control Dyn. 2003, 26, 536–542. [Google Scholar] [CrossRef]
  8. Julier, S.J.; Uhlman, J.K.; Durrant-Whyte, H.F. A new approach for filtering nonlinear systems. In Proceedings of the IEEE American Control Conference, Seattle, WA, USA, 21–23 June 1995; pp. 1628–1632. [Google Scholar]
  9. Uchiyama, S.; Takemoto, K.; Satoh, K.; Yamamoto, H.; Tamura, H. MR Platform: A basic body on which mixed reality applications are built. In Proceedings of the International Symposium on Mixed and Augmented Reality (ISMAR 2002), Washington, DC, USA, 30 September–1 October 2002; pp. 248–256. [Google Scholar]
  10. Arasaratnam, I.; Haykin, S. Cubature Kalman filters. IEEE Trans. Autom. Control 2009, 54, 1254–1269. [Google Scholar] [CrossRef]
  11. Sun, F.; Tang, L.J. Initial alignment of large azimuth misalignment angle in SINS based on CKF. Chin. J. Sci. Instrum. 2012, 33, 327–333. [Google Scholar]
  12. Zhang, Y.G.; Huang, Y.L.; Wu, Z.M.; Li, N. Moving state marine SINS initial alignment based on high degree CKF. Math. Probl. Eng. 2014, 2014, 546107. [Google Scholar] [CrossRef]
  13. Jia, B.; Xin, M.; Cheng, Y. High-degree cubature Kalman filter. Automatica 2013, 49, 510–518. [Google Scholar] [CrossRef]
  14. Jia, B.; Xin, M. Multiple sensor estimation using a new fifth-degree cubature information filter. Trans. Inst. Meas. Control 2015, 37, 15–24. [Google Scholar] [CrossRef]
  15. Zhang, Y.G.; Huang, Y.L.; Li, N. SINS Initial Alignment Based on Fifth-degree Cubature Kalman Filter. In Proceedings of the 2013 IEEE International Conference on Mechatronics and Automation (ICMA), Takamatsu, Japan, 4–7 August 2013; pp. 401–406. [Google Scholar]
  16. Xu, B.; Zhang, P.; Wen, H.; Wu, X. Stochastic stability and performance analysis of cubature Kalman filter. Neurocomputing 2016, 186, 218–227. [Google Scholar] [CrossRef]
  17. Ding, W.; Wang, J.; Rizos, C.; Kinlyside, D. Improving adaptive Kalman estimation in GPS/INS integration. J. Navigat. 2007, 60, 517–529. [Google Scholar] [CrossRef]
  18. Olfati-Saber, R. Kalman-Consensus filter: Optimality, stability, and performance. In Proceedings of the 48th IEEE Conference on Decision and Control held jointly with 2009 28th Chinese Control Conference, Shanghai, China, 15–18 December 2009; pp. 7036–7042. [Google Scholar]
  19. Fang, J.; Yang, S. Study on Innovation Adaptive EKF for In-Flight Alignment of Airborne POS. IEEE Trans. Instrum. Meas. 2011, 60, 1378–1388. [Google Scholar]
  20. Sun, J.; Xu, X.S.; Liu, Y.T.; Zhang, T.; Li, Y. Initial Alignment of Large Azimuth Misalignment Angles in SINS Based on Adaptive UPF. Sensors 2015, 15, 21807–21823. [Google Scholar] [CrossRef] [PubMed]
  21. Huang, X.Y.; Tang, X.Q.; Wu, M. Application of 5th-degree CKF in SINS nonlinear initial alignment. Syst. Eng. Electron. 2015, 37, 633–638. [Google Scholar]
  22. Yan, G.M.; Yan, W.S.; Xu, D.M. A SINS nonlinear error model reflecting better characteristics of SINS errors. J. Northwest. Polytech. Univ. 2009, 27, 511–516. [Google Scholar]
  23. Wu, Y.X.; Hu, D.W.; Wu, M.P.; Hu, X.P. A numerical-integration perspective on Gaussian filters. IEEE Trans. Signal Process. 2006, 54, 2910–2921. [Google Scholar] [CrossRef]
  24. Long, R.; Qin, Y.Y.; Jia, J.C. Observable degree analysis of SINS initial alignment based on singular value decomposition. In Proceedings of the 2008 IEEE International Symposium on Knowledge Acquisition and Modeling Workshop, Wuhan, China, 21–22 December 2008; pp. 444–448. [Google Scholar]
  25. Yan, G.M.; Yan, W.S. Application of simplified UKF in SINS initial alignment for large misalignment angles. J. Chin. Inert. Technol. 2008, 16, 253–264. [Google Scholar]
Figure 1. The distribution of weights of innovations under different ξ .
Figure 1. The distribution of weights of innovations under different ξ .
Sensors 18 00659 g001
Figure 2. The flow chart of IICKF5.
Figure 2. The flow chart of IICKF5.
Sensors 18 00659 g002
Figure 3. The estimated misalignment angles ( ϕ ( 0 ) = [ 20 ° ,   40 ° ,   170 ° ] ): (a) Pitching misalignment angle; (b) rolling misalignment angle; (c) azimuth misalignment angle.
Figure 3. The estimated misalignment angles ( ϕ ( 0 ) = [ 20 ° ,   40 ° ,   170 ° ] ): (a) Pitching misalignment angle; (b) rolling misalignment angle; (c) azimuth misalignment angle.
Sensors 18 00659 g003
Figure 4. The gradient values of azimuth angle error.
Figure 4. The gradient values of azimuth angle error.
Sensors 18 00659 g004
Figure 5. The estimated azimuth misalignment angles of CKF5, ICKF5 and IICKF5.
Figure 5. The estimated azimuth misalignment angles of CKF5, ICKF5 and IICKF5.
Sensors 18 00659 g005
Figure 6. The field test setup.
Figure 6. The field test setup.
Sensors 18 00659 g006
Figure 7. The alignment errors of large azimuth misalignment angle of dynamic vehicle experiment: (a) Pitch errors; (b) roll errors; (c) heading errors.
Figure 7. The alignment errors of large azimuth misalignment angle of dynamic vehicle experiment: (a) Pitch errors; (b) roll errors; (c) heading errors.
Sensors 18 00659 g007
Figure 8. The gradient values of azimuth angle error (dynamic vehicle experiment).
Figure 8. The gradient values of azimuth angle error (dynamic vehicle experiment).
Sensors 18 00659 g008
Figure 9. The alignment errors of different algorithms: (a) Pitch errors; (b) roll errors; (c) heading errors.
Figure 9. The alignment errors of different algorithms: (a) Pitch errors; (b) roll errors; (c) heading errors.
Sensors 18 00659 g009
Table 1. Comparison of alignment errors (°).
Table 1. Comparison of alignment errors (°).
EKFUKFCKF3CKF5ICKF5IICKF5
error of pitch (°)−7.02050.08990.12700.15900.08290.0291
error of roll (°)7.8234−0.1464−0.1634−0.2078−0.2949−0.2738
error of heading (°)62.7783157.736161.290213.6403−1.2458−0.0878
Table 2. Comparison of alignment errors with in-motion base (°).
Table 2. Comparison of alignment errors with in-motion base (°).
EKFUKFCKF3CKF5ICKF5IICKF5
error of pitch (°)−3.2814−1.2868−2.6771−3.71893.02705.9149
error of roll (°)−1.2960−0.28673.37515.4245−2.90301.8357
error of heading (°)−139.5622−266.3212−204.34577.998642.5004−1.9540
Table 3. Mean and Standard Deviation (SD) of azimuth misalignment errors (°).
Table 3. Mean and Standard Deviation (SD) of azimuth misalignment errors (°).
EKFUKFCKF3CKF5ICKF5IICKF5
Mean (°)−155.2618−264.1352−201.19621.6920−6.1625−0.6852
SD (°)38.89353.806211.535910.230741.20564.6229

Share and Cite

MDPI and ACS Style

Wang, W.; Chen, X. Application of Improved 5th-Cubature Kalman Filter in Initial Strapdown Inertial Navigation System Alignment for Large Misalignment Angles. Sensors 2018, 18, 659. https://doi.org/10.3390/s18020659

AMA Style

Wang W, Chen X. Application of Improved 5th-Cubature Kalman Filter in Initial Strapdown Inertial Navigation System Alignment for Large Misalignment Angles. Sensors. 2018; 18(2):659. https://doi.org/10.3390/s18020659

Chicago/Turabian Style

Wang, Wei, and Xiyuan Chen. 2018. "Application of Improved 5th-Cubature Kalman Filter in Initial Strapdown Inertial Navigation System Alignment for Large Misalignment Angles" Sensors 18, no. 2: 659. https://doi.org/10.3390/s18020659

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