Next Article in Journal
Baseball Player Behavior Classification System Using Long Short-Term Memory with Multimodal Features
Previous Article in Journal
Performance in Solar Orientation Determination for Regular Pyramid Sun Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Kalman Filters for Inertial Integrated Navigation

Department of Navigation Engineering, Naval University of Engineering, Wuhan 430000, China
*
Author to whom correspondence should be addressed.
Sensors 2019, 19(6), 1426; https://doi.org/10.3390/s19061426
Submission received: 25 January 2019 / Revised: 18 March 2019 / Accepted: 20 March 2019 / Published: 22 March 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
The current research on integrated navigation is mainly focused on filtering or integrated navigation equipment. Studies systematically comparing and analyzing how to choose appropriate integrated filtering methods under different circumstances are lacking. This paper focuses on integrated navigation filters that are used by different filters and attitude parameters for inertial integrated navigation. We researched integrated navigation filters, established algorithms, and examined the relative merits for practical integrated navigation. Some suggestions for the use of filtering algorithms are provided. We completed simulations and car-mounted experiments for low-cost strapdown inertial navigation system (SINS) to assess the performance of the integrated navigation filtering algorithms.

1. Introduction

An integrated navigation system is a system using two or more nonsimilar navigation systems to measure navigation information and calculate or correct the errors of the navigation system from measurements [1,2,3,4]. The integrated navigation system most commonly used is the inertial integrated navigation system because the navigation information is comprehensive and autonomous [5,6,7]. The main part of an inertial integrated navigation system is the inertial navigation system (INS). The INS is mainly divided into two types: platform inertial navigation system (PINS) and strapdown inertial navigation system (SINS). The PINS uses electromechanical control equipment to establish a physical platform, which places low requirements on the navigation computer. Its main disadvantages include its relatively complex structure, huge equipment volume, and heavy weight [8]. With the maturity of optical gyroscopic technology and the development of computers, SINS is gradually replacing PINS, becoming a more popular research topic [9]. SINS replaces the physical platform with a digital platform, which is simpler, smaller, and lighter. However, both PINS and SINS have disadvantages: the navigation accuracy diverges with time due to device errors and self-error-correction is unreliable [10,11,12]. As such, INS needs other navigation systems to form an integrated navigation system.
The inertial integrated navigation system has been used in many fields, such as land navigation, underwater navigation, and even aerospace and unmanned system navigation [13,14,15,16,17,18]. For example, these systems are usually integrated with a global positioning system (GPS), Roland C or an odometer for land navigation integration, and integrated with the Doppler Velocity Log (DVL) baseline system for underwater navigation integration. However, the technological key to achieving inertial integrated navigation is integrated filtering [19].
Filtering is an estimation method used for processing random data, and has been studied and applied as an estimator in related fields. Modern filtering techniques are represented by the Wiener filter and the Kalman filter [20]. The Wiener filter is a frequency domain filter, whereas the Kalman filter is designed in the time domain. The Wiener filter laid the foundation for modern control theory, but the filter experiences problems with real-time calculation and filtering design. Based on the least squares theory, Kalman designed a linear recursive Kalman filter (KF) algorithm. The KF algorithm is simple in structure and easy to execute and calculate on a computer. Although the KF is optimal in the estimation of linear systems, it cannot be applied to nonlinear systems. Farrell, Barth, and other scholars proposed the extended Kalman filter (EKF). It linearizes nonlinear system models using Taylor series and multivariate Jacobian matrices [21]. The EKF and its improved versions are the most widely used nonlinear Kalman filter algorithms in engineering. However, EKF has its drawbacks. If the degree of nonlinearity of the system is high, the estimate accuracy of EKF seriously decreases. In 1995, Julier and Uhlmann proposed the unscented Kalman filter (UKF) algorithm, which does not ignore the high-order terms of Taylor expansion by unscented transform (UT), and thus has high estimate accuracy in nonlinear systems [22]. Emerging filter algorithms, such as cubature Kalman filter (CKF) [23] and particle filter (PF) [24,25], have been theoretically verified in many fields. Zhou J. [24] and Jiao Z. [25] provided methods to reduce the computation complexity of the PF. The research on PF is still in the theoretical stage. Because of its huge calculation, PF is rarely used in integrated navigation filtering methods in engineering practice.
We can apply linear or nonlinear Kalman filters in inertial integrated navigation. These filters are called integrated navigation filters. However, this filtering application is not straightforward because integrated navigation filters have their own particularities in terms of attitude representations. Among the integrated navigation filters, the attitude representations vary. The common representations are rotation angles, the Euler angle, and the family of Rodrigues parameters, such as modified Rodrigues parameters (MRP), and quaternion [26,27]. The earliest integrated navigation filter was constructed from the Euler angle and inertial error equations, which is called Euler-KF in this paper. Since the inertial error equations are linear, the Euler-KF uses KF for filtering. Thus, the Euler-KF is simple and easy to apply. The application of Euler-KF can be traced back to the Apollo moon landing program. In the field of integrated navigation, the Euler-KF is traditionally called indirect integration due to its integrated structure.
Euler-KF can be used in many fields; however, it faces problems with model accuracy due to inertial error equations. To overcome the model accuracy problem with inertial error equations, inertial basic equations have been used in inertial integrated system models with the development of nonlinear filters. Compared with inertial error equations, the inertial basic equations are precision lossless models based on physical law. For inertial basic equations, the key point is how to represent the attitude transformation matrix and nonlinear filtering. The attitude transformation matrix describes a rotation of coordinates using attitude representations. It is the core of the SINS digital platform. For the attitude transformation matrix, different attitude representations can be used. Compared with three-dimensional (3D) attitude parameters like the Euler angle and the family of Rodrigues parameters, the constrained quaternion, as a four-dimensional attitude parameter, is more popular due to global nonsingularity. Then, the quaternion calculation rules can be divided into multiplicative quaternion and additive quaternion. Compared with additive quaternion, the multiplicative quaternion is more widely used because it is more consistent with the definition of quaternion. As the attitude transformation matrix uses quaternion, the filters need to use nonlinear Kalman filters such as EKF and UKF. Inertial integrated navigation based on EKF using multiplicative quaternion is called the multiplicative EKF (MEKF) algorithm.
L.J. Zhang [28] proposed a spacecraft attitude determination algorithm based on MEKF, and I.A. Ruicai et al. [29] applied the MEKF algorithm to low-precision microelectromechanical systems (MEMS) attitude estimation. Fangjun Qin et al. [30] proposed a sequential MEKF that updates the state covariance at each step of the sequential procedure. The essence of MEKF is EKF, so it inherits the disadvantages of EKF. In 2003, Crassidis first proposed a quaternion-based UKF for spacecraft attitude estimation, known as unscented quaternion estimator (USQUE) [31]. The USQUE is a layered filtering algorithm using a nonconstrained rotation vector that represents the attitude error to perform inner layer filtering recursion. In 2006, Crassidis introduced USQUE to land integrated navigation [32]. Zhou J., Edwan E., et al. [33] studied the application of USQUE for tightly integrated MEMS and GPS navigation system. Li Kailong et al. [34] modified USQUE to reduce the switch frequency between error-MRP and quaternion.
For attitude transformation matrixes, quaternion is not the only choice. Euler angles and the family of Rodrigues parameters can be also used for attitude transformation matrix. However, they have an attitude singularity problem. MRP can use its shadow MRP (SMRP) to avoid singularity. Karlgaard and Schaub proposed MRP-EKF and MRP-UKF using MRP and SMRP for inertial integrated navigation [35]. Recently, some studies focused on Euler angles for nonlinear inertial integrated filtering. With Euler angles, different rotation orders can form different Euler angles. Different Euler angles have different singularities. Thus, Euler angles can use a special rotation order to avoid singularity. Ran and Cheng [36] used double Euler angles under UKF for initial alignment, which is called DoEuler-UKF in this paper. Although still some problems exist with DoEuler-UKF, DoEuler-UKF has special advantages, like clear physical meaning and simpler algorithm structure. We think the DoEuler-UKF will be extensively researched.
In this paper, we focus on Euler-KF, MEKF, USQUE, MRP-UKF, and DoEuler-UKF. These five filters are currently common or research hotspots in integrated filtering methods. Although some studies proposed those five filters, systematic research comparing these filters in inertial integrated navigation is lacking. The published research on Kalman filters in the navigation field mainly focus on attitude estimation and filter improvement. To fill the literature gap about the characteristics and applicable settings of low-precision inertial navigation integrated navigation, this paper focuses on those filters, especially in terms of algorithm establishment, characteristics, the merits of the specific integrated navigation system, and the applicable environment. Through research and comparison of these five filters using simulation tests and a loosely coupled MEMS/GPS experiment, the performance of these filtering algorithms was systematically analyzed. As a result, we drew some conclusions about the various indicators of the five filters. Some suggestions are as a basis for the selection of integrated filtering methods in different situations.
The organization of this paper is as follows. In the review of existing theories, the basic SINS equations and the SINS error equations are introduced. In the filtering algorithm section, MEKF and USQUE are studied in detail. Finally, simulation testing and car-mounted experiments were conducted and are summarized in the experimental section. The results are provided and the performance of the five inertial integrated navigation filters are compared. A few meaningful conclusions are outlined in the summary.

2. Theoretical Review

In this section, the reference frames commonly used to derive the strapdown inertial integrated navigation system, the properties of basic SINS equations, and the error SINS equations are provided.
The coordinate frame is the research foundation of SINS. The SINS attitude transformation matrix is usually related to two important coordinate frames: body frame (b-frame, Right-Front-Up, RFU) and navigation frame (n-frame, East-North-Up, ENU). The b-frame is fixed onto the vehicle body and rotates with it. The outputs of IMU are expressed with respect to the b-frame. The n-frame is the reference frame that performs navigation calculation. Other coordinates used in this paper are the Earthframe (e-frame) and the inertial frame (i-frame). The outputs of IMU of SINS should be transformed from b-frame to n-frame for navigation calculation using the attitude transformation matrix. The vector rb with respect to b-frame can be transformed to rn with respect to n-frame:
r n = C b n r b
C b n = C n b T
where C b n is the attitude transformation matrix used to describe a rotation of coordinates using different attitude representations.

2.1. Basic SINS Equations

Basic SINS equations are the core of the system model of integrated filtering methods, which contain an attitude part, a velocity part, a position part, a gyro measurement part, and an accelerometer measurement part. For basic SINS equations, the key parts are the attitude transformation matrix and attitude kinematic equation because they form the essence of the SINS digital platform. The attitude transformation can be represented by the Euler angle, MRP, and quaternion. When using the Euler angle, the attitude kinematic equation is given by
θ ˙ γ ˙ ψ ˙ = 1 cos θ cos θ cos γ 0 cos θ sin γ sin θ sin γ cos θ cos γ sin θ sin γ 0 cos γ ω n b b
where φ = θ ; γ ; ψ , and θ is pitch, γ is roll, ψ is yaw, ω n b b is angular velocity of the b-frame relative to the n-frame, ω i b b is the gyro output, and ω i n n is the angular velocity of the n-frame relative to the i-frame. The attitude transformation matrix C n b ( φ ) using the Euler angle is given by
C n b ( φ ) = cos γ cos ψ + sin γ sin ψ sin θ sin ψ cos θ sin γ cos ψ cos γ sin ψ sin θ cos γ sin ψ + sin γ cos ψ sin θ cos ψ cos θ sin γ sin ψ cos γ cos ψ sin θ sin γ cos θ sin θ cos γ cos θ
When using MRP, the attitude kinematic equation is
σ ˙ = 1 4 1 σ 2 I 3 × 3 + 2 σ × + 2 σ σ T ω n b b
where σ is MRP and I 3 × 3 is a unit matrix. The attitude transformation matrix C n b ( σ ) using MRPis given by
C n b ( σ ) = I 3 × 3 σ × I 3 × 3 + σ × 2 = I 3 + 8 σ × 2 4 1 σ 2 σ × 1 + σ 2 2
Then, the attitude kinematic equation is:
q ˙ = 1 2 Ξ ( q ) ω n b b = 1 2 Ω ( ω n b b ) q
where q = ρ ; q 4 is a quaternion with ρ = q 1 ; q 2 ; q 3 vector part and q 4 scalar part and Ω ( ω n b b ) = [ ω n b b × ] ω n b b ( ω n b b ) T 0 . Then, the × cross product matrix is defined by
x × = 0 x 3 x 2 x 3 0 x 1 x 2 x 1 0   with   x = x 1 ; x 2 ; x 3
The attitude transformation matrix C b n ( q ) using a quaternion is
C n b ( q ) = Ψ ( q ) T Ξ ( q ) = q 1 2 q 2 2 q 3 2 + q 4 2 2 ( q 1 q 2 q 3 q 4 ) 2 ( q 1 q 3 + q 2 q 4 ) 2 ( q 1 q 2 + q 3 q 4 ) q 1 2 + q 2 2 q 3 2 + q 4 2 2 ( q 2 q 3 q 1 q 4 ) 2 ( q 1 q 3 q 2 q 4 ) 2 ( q 2 q 3 + q 1 q 4 ) q 1 2 q 2 2 + q 3 2 + q 4 2
where
Ξ ( q ) = q 4 I 3 × 3 + ρ × ρ T   and   Ψ ( q ) = q 4 I 3 × 3 ρ × ρ T
where I3×3 is a 3 × 3 unit matrix. The attitude kinematic equation is different using different attitude representations. This is the core of basic SINS equations. Then, other parts of the basic SINS equations are [37]
λ ˙ = v N R M + h
ϕ ˙ = v E ( R N + h ) cos λ
h ˙ = v U
v ˙ E = v E ( R N + h ) cos λ + 2 ω i e e v N sin λ v E v U R N + h 2 ω i e e v U cos λ + f E
v ˙ N = v E ( R N + h ) cos λ + 2 ω i e e v E sin λ v E v U R M + h + f N
v ˙ U = v E 2 R N + h + v N 2 R M + h + 2 ω i e e v E cos λ g f U
where (9)–(11)are position kinematic equations and (12)–(14) are velocity kinematic equations; is the position, λ is the latitude, ϕ is the longitude, and h is the height; v = v E ; v N ; v U is the velocity; g is the gravity vector; is the Earth’s rotation rate, which is 7.292115 × 10−5 rad/s (WGS-84); is the specific force vector in n-frame, expressed in the b-frame by; ω i b b is the gyro measurements; and R M and R M are the radius of Earth:
R M = R e ( 1 e 2 ) / ( 1 e 2 sin 2 λ ) 1.5
R N = R e / ( 1 e 2 sin 2 λ ) 0.5
where R e = 6,378,137 m and e = 0.0818. The gyro measurement equation is
ω ˜ i b b = ω i b b + ε + η g v
ε ˙ = η g u
where ω ˜ i b b is the outputs of the gyro with the gyro bias; η g v and η g u are zero-mean Gaussian white-noise process with spectral densities given by σ g v 2 and σ g u 2 , respectively. The accelerometer measurement equation is
f ˜ b = f b + + η a v
˙ = η a u
where f ˜ b is the outputs of the accelerometer with accelerometer bias , and η a v and η a u are zero-mean Gaussian white-noise process with spectral densities given by σ a v 2 and σ a u 2 , respectively.

2.2. SINS Error Equations

The SINS error equations can be derived from the basic SINS equations. As with basic SINS equations, SINS error equations have an attitude error part, velocity error part, and position error part. The key process derived from basic SINS equations is how to present the attitude error equation. The attitude error equation reflects the transition from quaternion to Euler angle error δ φ . The velocity error equation and position error equation are mainly formed by assuming the approximation condition, simplification, and ignoring the infinitesimal term. As such, the SINS error equations can be obtained. By analyzing the derivation principle of the SINS error equations, the SINS error equations affect the integrated navigation estimate effect to a large extent in some cases when the motion is severe or under large misalignment angles. For SINS error equations, the importance is the attitude equations using the Euler angle error. The attitude error equations are
δ φ ˙ E = δ φ N ω i e e sin λ + v E tan λ R N + h δ φ U ω i e e cos λ + v E R N + h δ v N R M + h + δ h v N ( R M + h ) 2 ε E b
δ φ ˙ N = δ φ E ω i e e sin λ + v E tan λ R N + h δ φ U v N R M + h δ λ ω i e e sin λ + δ v E R N + h δ h v E ( R N + h ) 2 ε N b
δ φ ˙ U = δ φ E ω i e e cos λ + v E R N + h + δ φ N v N R M + h + δ λ ω i e e cos λ + v E R N + h sec 2 λ    + δ v E tan λ R N + h δ h v E tan λ ( R N + h ) 2 ε U b
where δ ( ) is the position or velocity error. The other velocity or position parts can be obtained from corresponding references. Then, the SINS error equations are linear and the KF will be optimal for solving the filtering problem. Four aspects of connections and differences between basic SINS equations and error equations are summarized as follows:
(1)
The SINS error equations are formed by the equivalent transformation and approximation processing of the basic SINS equations
(2)
The SINS error equations are linear equations, and the basic SINS equations are nonlinear equations.
(3)
Compared with the SINS error equations, the basic SINS equations directly reflect the change in navigation parameters, and therefore can more accurately reflect the actual motion of the carrier.
(4)
For the attitude error equation, the linearization premise is that the attitude error angle is a small angle, so for large misalignment angles, the equation loses the linear assumption and becomes a nonlinear equation.

3. Kalman Filters for Inertial Integrated Navigation

For SINS, the first step is to determine the attitude representation. The second step involves selecting the corresponding Kalman filter according to the integrated navigation system model. The relationships between attitude representation, model, and integrated filtering methods are shown in Figure 1.
For the SINS error equations, the Euler-KF is used for estimation. For the basic SINS equations, nonlinear Kalman filters can be used for estimation. Then, the nonlinear Kalman filters, like EKF and UKF with different attitude representations, can form different integrated navigation filters. EKF includes additive EKF (AEKF) and MEKF using quaternion. Compared with MEKF, AEKF is rarely used in engineering. Thus, we mainly studied MEKF. UKF includes additive UKF (AUKF) and USQUE using a quaternion. Although AUKF may exist in theory, we were unable to find any published literature on this filter. Thus, we mainly studied USQUE. Finally, some 3Dattitude representations, like MRP and Euler angle, were used in UKF for integrated navigation filters: MRP-UKF and DoEuler-UKF. For the 3D integrated navigation filters, the main problem is how to avoid singularity. Thus, we discussed the singularity avoidance of MRP-UKF and DoEuler-UKF.

3.1. Multiplicative Extended Kalman Filter

The attitude quaternion parameter normalization constraint is easily satisfied for the calculation of a single quaternion. However, in the EKF nonlinear filter calculation, problems may be encountered. In addition, the covariance matrix in EKF is 3 × 3, whereas a quaternion is a four-dimensional parameter, thus EKF cannot be directly applied to the strapdown inertial integrated navigation system with a quaternion due to the dimension mismatch. The basic idea of MEKF is to reduce the dimension by a small error quaternion. The unconstrained 3D attitude error parameters were estimated using multiplicative quaternions to provide global nonsingular attitude descriptions. For the nonlinear model using nonlinear function local linear features, we applied the first-order Taylor expansion to the model, and obtained the Jacobian matrix of the model, so the Jacobian matrix could be applied to the propagation of the covariance matrix. Since the difference between the MEKF algorithm and the EKF algorithm is caused by attitude, we emphasized the attitude part of the algorithm. The attitude part mainly influences the propagation of covariance.
The quaternion has no physical meaning. The MEKF algorithm reduces the dimensionality of the error quaternion by approximating the quaternion attitude and the real quaternion attitude δ q , which is reduced to the three-dimensional δ α , so that the Jacobian matrix can be obtained. The covariance matrix can be propagated. In the measurement update, due to the quaternion attitude, the quaternion attitude cannot be expressed in a unified form with the position, velocity, or other information about the 3D vector. Therefore, the measurement update of the attitude is expressed by the quaternion. After the measurement update, the state update is added to complete the estimation of the states.
In the following part, we derive the MEKF process combined with the strapdown inertial integrated navigation system. The state chooses quaternion attitude, position, velocity, gyro bias, and accelerometer bias. Due to the dimension problem, the state vectors are divided into an attitude part and an ‘other’ part, so we define the state x ^ k 1 = x ^ k 1 q , x ^ k 1 e T and covariance P k 1 . x ^ k 1 q is quaternion attitude and x ^ k 1 e is the other states x ^ k 1 e = p ^ k 1 , v ^ k 1 , ε ^ k 1 , ^ k 1 T . The discrete time state inertial integrated navigation process model can be given as x k = f ( x k 1 ) + w k 1 .
The basic goal of MEKF is to reduce the dimension with an error quaternion. MEKF uses the error quaternion as the state vector in the time update covariance propagation using the multiplicative rule. MEKF starts by the definition of the error quaternion. The error quaternion is defined by the multiplicative error as
δ q = q q ^ 1   with   δ q = δ ρ ; δ q 4 .
The attitude error matrix is given by
C b n ( δ q ) = C b n ( q ) C b n ( q ^ ) T .
For small angles the vector part of the quaternion is approximately equal to half angles, then δ ρ δ α / 2 and δ q 4 1 , where δ α is a small error-angle rotation vector. The linearized model error-kinematics is as follows:
δ α ˙ = [ ω ^ n b b × ] δ α + δ ω i b b C b n ( q ^ ) δ ω i n n
δ q ˙ 4 = 0
where δ ω i b b = ω i b b ω ^ i b b and δ ω i n n = ω i n n ω ^ i n n . δ ω i n n can be computed by first-order Taylor series expansion. Usually, the true quaternion is close to the estimated quaternion. Therefore, we can reduce the quaternion to a three-dimensional vector:
δ ω i b b = [ I 3 × 3 Δ ε + I 3 × 3 η g v ]   with   Δ ε = ε ε ^
Then,
δ α ˙ = ( ω ˜ i b b ε ^ ) × δ α Δ ε η g v C n b ( q ) ω i n n p p ^ , v ^ Δ p C n b ( q ) ω i n n v p ^ Δ v
where Δ p = p p ^ and Δ v = v v ^ . The partials are shown as
ω i n n p = v N ( R M + h ) 2 R M λ 0 v N ( R M + h ) 2 ω i e e sin λ v E ( R N + h ) 2 R N λ 0 v E ( R N + h ) 2 ω i e e cos λ + v E sec 2 λ R N + h v E tan λ ( R N + h ) 2 R N λ 0 v E tan λ ( R N + h ) 2
ω i n n v = 0 1 R M + h 0 1 R N + h 0 0 tan λ R N + h 0 0
where R N λ = R e e 2 sin λ cos λ ( 1 e 2 sin n 2 λ ) 3 / 2 and R M λ = 3 R e ( 1 e 2 ) e 2 sin λ cos λ ( 1 e 2 sin 2 λ ) 5 / 2 . After the dimension reduction is completed, the calculation of the Jacobian matrix can be performed to propagate the covariance matrix. Then, we determine the state error Δ x = δ α q ^ k 1 , Δ p ^ k 1 , Δ v ^ k 1 , Δ ε ^ k 1 , Δ ^ k 1 T . The state, state error, noise vector, and covariance used in MEKF are given by
x = q ; p ; v ; ε ; , Δ x = δ α ; Δ p ; Δ v ; Δ ε ; Δ , w = η g v ; η g u ; η a v ; η a u
Q = σ g v 2 I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 σ g u 2 I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 σ a v 2 I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 σ a u 2 I 3 × 3
where 0 3 × 3 is a 3 × 3 zeros matrix. After propagating filtering state, MEKF propagates covariance using an error equation. Then, the error equation is given by
Δ x ˙ = F Δ x + G w
where
F = F 11 F 12 F 13 F 14 0 3 × 3 0 3 × 3 F 22 F 23 0 3 × 3 0 3 × 3 F 31 F 32 F 33 0 3 × 3 F 35 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 .
The system noise is w = η g v ; η g u ; η a v ; η a u and covariance is
Q = σ g v 2 I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 σ g u 2 I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 σ a v 2 I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 σ a u 2 I 3 × 3
and the noise coefficient is
G = I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 C b n ( q ^ ) 0 3 × 3 0 3 × 3 I 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 0 3 × 3 I 3 × 3
with Table 1.
The position partial differentials are
p ˙ p = v N ( R M + h ) 2 R M λ 0 v N ( R M + h ) 2 v E sec λ ( R N + h ) 2 R N λ + v E sec λ tan λ R N + h 0 v E sec λ ( R N + h ) 2 0 0 0
p ˙ v = 0 1 R M + h 0 sec λ R N + h 0 0 0 0 1
The velocity partial differentials are
v ˙ p = Y 11 0 Y 13 Y 21 0 Y 23 Y 31 0 Y 33
v ˙ v = Z 11 Z 12 Z 13 Z 21 Z 22 Z 23 Z 31 Z 32 0
with Table 2 and Table 3.
g λ = 9.780327 1.06048 × 10 2 sin λ cos λ 4.64 × 10 5 ( sin λ c o s 3 λ sin 3 λ cos λ )
g h = 3.0877 × 10 6 + 4.4 × 10 9 sin 2 λ + 1.44 × 10 13 h
In the measurement update, either the method of tightly coupled or super tightly coupled is adopted, then the measurement equation is nonlinear. For the MEKF algorithm, the filtering process of the measurement update is basically the same as the time update. Using the loosely-coupled model, the measurement equation is a linear equation. The difference from Kalman filtering is that, since the quaternion attitude dimension does not match other state quantities, δ α ^ k and δ x ^ k e are separately expressed in the measurement update. We introduce the measurement update by taking the loosely-coupled model as an example.
y k = H k x k + η k
where H k varies with the change in the measurement. When only using position information as the loosely-coupled measurement transfer matrix,
H p , k = [ 0 3 × 3 I 3 × 3 0 3 × 9 ]
When only using velocity information as the loosely-coupled measurement transfer matrix:
H v , k = [ 0 3 × 6 I 3 × 3 0 3 × 6 ] .
When using both position and velocity information as the loosely-coupled measurement transfer matrix:
H p , v , k = 0 3 × 3 I 3 × 3 0 3 × 3 0 6 × 3 0 3 × 3 0 3 × 3 I 3 × 3 0 6 × 3 .
In the state update, the entire filtering process is completed by updating the state vectors by using the errors. The MEKF algorithm is summarized in Algorithm 1.
Algorithm 1: The multiplicative extended Kalman filter (MEKF) algorithm.
Time Propagation:
Define the filtering state x ^ k 1 q ^ and P k 1 .
Propagate filtering state and covariance
x ^ k / k 1 ( q ^ ) = f x ^ k 1 ( q ^ )
P k / k 1 = Φ k 1 P k 1 Φ k 1 T + W k 1
where Φ k 1 is the discrete-time state transition matrix obtained via numerical solution.
Measurement update:
Compute filtering gain K k = P k / k 1 H k T ( H k P k / k 1 H k T + R k ) 1
Compute state error δ x ^ k = K k y k h ( x ^ k / k 1 ) = δ α ^ k δ x ^ k e with δ x ^ k e = δ p ^ k ; δ v ^ k ; δ ε ^ k ; δ ^ k
Compute quaternion states
q ^ k = q ^ k / k 1 + 1 2 Ξ ( q ^ k / k 1 ) δ α ^ k
Compute other states
x ^ k e = x ^ k / k 1 e + δ x ^ k e
Compute filtering covariance P k = ( I n K k H k ) P k / k 1

3.2. Unscented Quaternion Estimator

Because UKF cannot be directly applied to the strapdown inertial integrated navigation system with quaternion presentation attitude, guaranteeing solutions to the attitude quaternion constraint problem of the attitude quaternion parameter and the matching constraint problem of the quaternion filter variance matrix is difficult. USQUE solves the problem of UKF’s attitude estimation in the equation of state is “layered filtering”. The filtering process is divided into two layers: the inner layer filtering uses the generalized error Rodrigues parameter to update the attitude state, and the outer layer filtering uses the quaternion to perform the attitude. The quaternion and the generalized error Rodrigues parameter use the attitude transformation relationship of the multiplicative error quaternion as a bridge for switching. The basic idea of USQUE is shown in Figure 2.
According to Figure 2, the USQUE filtering process can be summarized. The equation of state and the measurement equation have been described in detail above. Similar to MEKF, the state is divided into two parts: attitude vector and nonattitude vector. The nonattitude part is like MEKF but the attitude part is a generalized Rodrigues parameter (GRP) x ^ k 1 δ , where the state vector is x ^ k 1 = x ^ k 1 δ ; x ^ k 1 e . The corresponding filter covariance matrix is P k 1 . The USQUE algorithm is used to estimate x ^ k δ and then the corresponding quaternion attitude x ^ k q in the state estimation. For the loose-coupled method, the measurement equation is a linear equation and the measurement update process is the same as for the linear Kalman filter. After the measurement update is completed, the attitude GRP is converted to a quaternion, which is the purpose of the attitude update. The algorithm is summarized in Algorithm 2.
Algorithm 2: The unscented quaternion estimator (USQUE) algorithm.
Time Propagation:
Define the filtering state
Using error MRP describe attitude x ^ k 1 ( δ ) = x ^ k 1 δ ; x ^ k 1 e
Using a quaternion describe attitude x ^ k 1 ( q ) = x ^ k 1 q ; x ^ k 1 e and covariance P k 1
Create sigma points using UT χ k 1 i = χ k 1 δ i χ k 1 e i = s i g m a x ^ k 1 δ , P k 1
The χ k 1 δ i corresponding quaternion error χ k 1 δ q i = δ ρ k 1 i ; δ q 4 , k 1 i
δ q 4 , k 1 i = a χ k 1 δ i 2 + f f 2 + 1 a 2 χ k 1 δ i 2 f 2 + χ k 1 δ i 2
δ ρ k 1 i = f 1 a + δ q 4 , k 1 i χ k 1 δ i
Compute the quaternion sigma points χ k 1 q i = χ k 1 δ q i q ^ k 1 , then χ k 1 i = χ k 1 q i ; χ k 1 e i
Propagate sigma points χ k / k 1 i = f χ k 1 i = χ k / k 1 q i ; χ k / k 1 e i
Compute error quaternion sigma points
χ k / k 1 δ q i = χ k / k 1 q i x k / k 1 q ¯ 1 and χ k / k 1 δ q i = δ ρ k / k 1 i ; δ q 4 , k / k 1 i
Compute GRP sigma points
χ k / k 1 δ i = f δ ρ k / k 1 i a + δ q 4 , k / k 1 i and χ k / k 1 i = χ k / k 1 δ i ; χ k / k 1 e i
The predicted mean and covariance of state
x ^ k / k 1 = i = 0 2 n w i χ k / k 1 i
P k / k 1 = i = 0 2 n w i χ k / k 1 i x ^ k / k 1 χ k / k 1 i x ^ k / k 1 T + W k 1
Measurement update:
Compute filtering gain,
K k = P k / k 1 H k T H k P k / k 1 H k T + R k 1
Compute posterior mean and covariance,
x ^ k = x ^ k / k 1 + K k y k H k x ^ k / k 1 and x ^ k = x ^ k 1 δ ; x ^ k e
P k = ( I n × n K k H k ) P k / k 1 with I n × n is an n-dimension unit matrix
Attitude update:
The x ^ k δ corresponding quaternion is
x ^ k δ q i = δ ρ k i , δ q k , 4 i T
with
δ q 4 , k = a x ^ k δ 2 + f f 2 + 1 a 2 x ^ k δ 2 f 2 + x ^ k δ 2 and δ ρ k = f 1 a + δ q 4 , k x ^ k δ
The update of quaternion is given by x ^ k q = x ^ k δ q x k / k 1 q ¯ and δ k 0 , x ^ k ( q ) = x ^ k q ; x ^ k e
Then enter the next filtering cycle.

3.3. Three-Dimensional Integrated Navigation Filters

For 3D attitude representations, like MRP and Euler angle, their integrated navigation filters help avoid singularity. In this section, we briefly show this problem.
According to definition of the family of Rodrigues parameters:
R P = e tan ( ϑ 2 N )
where the MRP σ = e tan ( ϑ 4 ) when N = 2, and ϑ is the attitude rotation vector. Thus, we know that the singularity appears at 2 π ± 4 n π where n is a positive integer. Traditionally, the MRP can be transferred to SMRP σ s to avoid singularity. The relationship between MRP and SMRP is:
σ s = σ σ 2
Since the SMRP is defined by σ s = e cot ( ϑ 4 ) , the singularities of MRP and SMRP are different. Thus, the MRP-UKF exploits this feature to avoid singularity. For MRP-UKF, we not only consider the switch in attitude representation, but also the switch in filtering covariance. The switch in attitude representation affects the stability of filtering covariance. To stabilize the filtering covariance, the switch of filtering covariance in MRP-UKF can be shown by:
A = 1 ( σ ) T σ I 3 × 3 + 2 σ × + 2 σ ( σ ) T
T = d i a g A ; I n × n
A s = 1 ( σ s ) T σ s I 3 × 3 + 2 σ s × + 2 σ s ( σ s ) T
T s = d i a g A s ; I n × n
P s = ( T s T T ) P ( T s T T ) T 1 + σ T σ 4
For the Euler angle in UKF compared with MRP, the singularity of the Euler angle is complex because the singularity is related to the rotate order of Euler angle for the attitude transformation matrix. Different rotate orders of the Euler angle in the attitude transformation matrix will have different singularities. The DoEuler-UKF exploits this feature to avoid singularity. DoEuler-UKF uses two different rotate orders of Euler angles to avoid singularity. One of the Euler angles is called forward Euler angle (Figure 3a) and the attitude transformation matrix is outlined by Equation (4) and the corresponding attitude kinematic equation is outlined in Equation (3).
The other Euler angle is called the inverted Euler angle (Figure 3b) and the attitude transformation matrix is given by:
C n b ( φ r ) = cos γ r cos ψ r cos γ r sin ψ r sin γ r cos ψ r sin γ r sin θ r cos θ r sin ψ r cos ψ r cos θ r + sin ψ r sin γ r sin θ r cos γ r sin θ r sin ψ r sin θ r + cos ψ r cos θ r sin γ r cos θ r sin ψ r sin γ r cos ψ r sin θ r cos γ r cos θ r
where φ r = θ r ; γ r ; ψ r is the inverted Euler angle. Then, the corresponding attitude kinematic equation is given by:
θ ˙ r γ ˙ r ψ ˙ r = 1 cos γ r cos γ r sin γ r sin θ r sin γ r cos θ r 0 cos θ r cos γ r sin θ r cos γ r 0 sin γ cos θ r ω n b b
The singularity of the forward Euler angle occurs when pitch θ ± π / 2 and the singularity of the inverted Euler angle occurs when roll γ r ± π / 2 . Thus, the DoEuler-UKF uses this different singularity to switch the attitude to avoid singularity. Compared with MRP-UKF, the DoEuler-UKF do not seem to need switching of the filtering covariance due to the lack of change in the attitude representation.

3.4. Discussionof Integrated Navigation Filters

To highlight the characteristics of each algorithm, some remarks are provided and summarized.
Firstly, the states of the Euler-KF are state errors. The linear model of the Euler-KF system model involves accuracy lossy equations. Therefore, the estimation accuracy is low. However, the calculation amount is small. The integrated structure of navigation devices is more flexible and easier to implement and may be used to form a multi model system.
Secondly, the states of the MEKF are navigation parameters. The system model of MEKF includes basic SINS equations, which are accuracy lossless equations. The research core of MEKF involves determining how to acquire the attitude error to propagate filtering covariance. This is mainly reflected in Equations (13)–(15). MEKF has the advantages of EKF, but also has the disadvantages of EKF. Although MEKF does not have model accuracy loss, filtering accuracy loss occurs. Because MEKF is also easy to implement, it is widely used in engineering.
Thirdly, the states of USQUE are navigation parameters. The system model of USQUE includes basic SINS equations, which are also accuracy lossless equations. UKF has higher accuracy compared to EKF. However, the quaternion cannot be directly applied to UKF. The USQUE is centered on the idea of layered filtering. The filtering process is divided into two layers: the inner layer filtering uses GRP to update the inner attitude state, and the outer layer filtering uses the quaternion to propagate the global attitude state. This attitude transformation is depicted in Figure 2. Compared with other algorithms, USQUE has a unique attitude update that transforms GRP to a quaternion.
Finally, MRP-UKF and DoEuler-UKF are 3D attitude integrated navigation filters. The states of those filters are navigation parameters. The system models are also basic SINS equations. Compared with MEKF and USQUE using a global nonsingular quaternion, MRP-UKF and DoEuler-UKF avoid the attitude singularity problem. Some other studies focused on how to realize the switching of filtering covariance or attain stable filtering.

4. Simulation Test and Experiments

A simulation test and a car-mounted experiment were conducted to comprehensively compare the performance of the five filtering methods. The car-mounted experiment used MEMS/GPS integration. From the analysis of the results, we determined the merits of the different integrated navigation filters under different conditions. Thus, we are able to provide some useful suggestions for the selection of integrated navigation filters.

4.1. Simulation

We compared the estimation performance for attitude, position, velocity, and computing. We conducted a loosely coupled SINS/GPS test with velocity as the measurement. The simulation test trajectory will be shown in Figure 4. The total simulation time was N = 1300 s. The total Monte-Carlo time was M = 50. In this simulation test, two metrics based on the mean error (ME) were used for the filtering accuracy check: ME1 is the performance of a single Monte-Carlo run and ME2 is the performance of the whole Monte-Carlo runs expressed in simulation time.
ME1 is defined as:
M E 1 ( m ) = 1 N k = 1 N x r e f , k x ^ k k m     m = 1 , 2 , , M
ME2 is defined as:
M E 2 ( k ) = 1 M m = 1 M x r e f , k x ^ k m k     k = 1 , 2 , , N
where x r e f , k is the reference state, x ^ k is the estimate, k is the number of simulation, and m is the number of Monte-Carlo runs. In this test, ME1 was used for the position and velocity estimates and ME2 was used for the attitude estimate. The initialization conditions are shown in Table 4.
The states are defined by 15-dimensional vector include attitude, velocity, position, and inertial device errors.
P 0 = d i a g P q 0 ; P v 0 ; P p 0 ; P ε 0 ; P 0 2   with   P q 0 = [ 3.0462 e 4 , 3.0462 e 4 , 3.0462 e 4 ] T
P v 0 = [ 0.01 , 0.01 , 0.01 ] T ,     P p 0 = [ 2.4582 e 14 , 2.4582 e 14 , 1 ] T ,
P ε 0 = [ 2.3504 e 11 , 2.3504 e 11 , 2.3504 e 11 ] T ,   P 0 = [ 9.5655 e 5 , 9.5655 e 5 , 9.5655 e 5 ] T
Q 0 = d i a g 0.01 ; 0.001 ; 100 ; 10 2
R 0 = d i a g R v 0 2 R v 0 = [ 0.01 , 0.01 , 0.01 ] T
In the loosely-coupled strapdown inertial integrated navigation system with velocity as the measurement, the estimation performance for attitude, position, and velocity of the five methods are shown in Figure 5, Figure 6 and Figure 7 and Table 5. The average estimation errors equations in Table 5 are given by:
δ x ^ = i = 0 T m = 1 M x ^ i m M T i = 0 , 1 , , T m = 1 , 2 , , M
where δ x ^ is the average estimation error, M is the number of total Monte-Carlo simulations, T is the total simulation time, and x ^ i m is the absolute value of the error.
According to the simulation results, regardless of attitude position and velocity, the estimation accuracies of UKF series filters are better than EKF series filters. The nonlinear model filter is more accurate than the linear model filter. USQUE, DoEuler-UKF, and MRP-UKF are almost equal in estimation accuracy, whereas Euler-KF is the worst. In attitude estimation, MEKF and USQUE are more stable than Euler-KF, DoEuler-UKF, and MRP-UKF. In terms of computation time, Euler-KF and MEKF have a considerable advantage; DoEuler-UKF is also far better than USQUE and MRP-UKF. When comparing MRP-UKF and USQUE, MRP-UKF is slightly better. The equations in USQUE, MEKF, MRP-UKF, and DoEuler-UKF are nonlinear basic SINS equations, whereas the Euler-KF equations are linear SINS error equations. To meet the linear model requirements, SINS error equations usually include partial ellipsis and linearization approximation, but SINS equations are more precise. As a result, USQUE, MEKF MRP-UKF, and DoEuler-UKF need more computation time. MEKF is a special form of EKF that is worse than UKF in estimation accuracy. Although MEKF is linearized, the Jacobian matrix of MEKF is created to propagate covariance and is unrelated to the propagation of the state. MEKF is worse than UKF series filters, but better than Euler-KF in terms of estimation accuracy. For UKF series filters researched in this paper, MRP-UKF, and DoEuler-UKF, due to the switching problem, the estimation accuracy is also slightly inadequate. As the model used in USQUE has the highest degree of nonlinearity, it has the longest computation time.
Fortightly-coupled or super tightly-coupled systems, the measurement model becomes nonlinear. USQUE MRP-UKF and DoEuler-UKF require more computation time and the estimation accuracy of Euler-KF worsens. MEKF is the first choice both in accuracy or computation time. Through many iterations, MEKF estimation accuracy is almost equal to USQUE and it needs less computation time. This explains why EKF series filters are widely used in engineering.

4.2. Car-Mounted Experiment

In the car-mounted experiments, MEMS was used for low-precision SINS to examine the performance of the five filters. The car-mounted experimental platform included MEMS strapdown inertial measurement equipment XW-IMU5220, a navigation-grade ring laser SINS for attitude reference, and a GPS receiver.
For the MEMS/GPS car-mounted experiment, we chose attitude, velocity, position, and gyro bias for comparison. In this car-mounted experiment, we focused on the estimated performance of attitude and gyro bias with the velocity and position measurements obtained from the GPS. The specifications of the gyroscopes and accelerometers of the MEMS IMU are listed in Table 6. The accuracy of the GPS receiver was 0.1 m/s for velocity and less than 2 m for position. Since the performance of high-precision laser SINS is much better than the low-precision MEMS, a highly-precise laser SINS was chosen for the attitude reference system. The specifications of MEMS and some initialization conditions are provided in Table 7.
P 0 = d i a g P q 0 ; P v 0 ; P p 0 ; P ε 0 2
with   P q 0 = [ 3.0462 e 4 , 3.0462 e 4 , 3.0462 e 4 ] T ,   P v 0 = [ 0.01 , 0.01 ] T ,   P p 0 = [ 2.4582 e 14 , 2.4582 e 14 ] T
P ε 0 = [ 2.3504 e 11 , 2.3504 e 11 , 2.3504 e 11 ] T
Q 0 = d i a g 0.01 ; 0.001 ; 100 ; 10 2
R 0 = d i a g R v 0 ; R p 0 2
R v 0 = [ 0.01 , 0.01 ] T ;   R p 0 = [ 2.4582 e 14 , 2.4582 e 14 ] T
The experiment test trajectory involved driving the car around on the ground. The experiment time was 250 s.
For the MEMS-based SINS, the estimation performance for attitude and gyro bias of the five filters are displayed in Figure 8 and Figure 9 and Table 8. The average estimation errors equation in Table 8 are given by:
δ x ^ = i = 0 T x ^ i T i = 0 , 1 , , T
where δ x ^ is the average estimation error, T is the total experience time, and x ^ i is the absolute value of the error.
According to the experiment, both in attitude and gyro bias, it is found that the estimation performance of USQUE is better than MEKF. This is because the MEKF algorithm uses the same model as USQUE. However, in the attitude estimation part, in order to calculate the Jacobian matrix, the δ α approximation is used instead of δ q , with some precision loss. In addition, the use of the Jacobian matrix also has loss of precision, and Euler-KF uses the SINS error equation. The model makes some approximations in the linearization process, so that although the Kalman filter is optimal in the linear model estimation, because of the model’s precision loss, the estimation accuracy is the worst. The accuracy of the model and calculation method is a double-edged sword. The high precision comes at the cost of calculation. DoEuler-UKF has a high estimation performance both in accuracy and computation time. When comparing with USQUE, the accuracy is almost equal and the computation time is half, because of the Euler-angle having less computation than quaternion. When using the MRP-UKF estimate of pitch and roll, the accuracy is pretty good, but when estimating gyro bias, the accuracy of MRP-UKF is a little poor.
According to the results and analysis of the above integrated navigation experiments, we drew some conclusions about the various indicators of the five filters. Table 9 summarizes the five filtering algorithm properties, characteristics, and estimation performances.
Table 9 shows that Euler-KF requires the least amount of calculation, but the estimation accuracy is the worst due to its linear model. Euler-KF is a mature integrated filtering method. The amount of calculation of MEKF is fairly small and the estimation accuracy is moderate. Although MEKF does not lose model accuracy, it loses filtering accuracy. MEKF is widely used in engineering. USQUE requires a huge amount of calculation, but has the best estimation accuracy. At present, USQUE is a research hotspot. MRP-UKF requires a large amount of calculation and moderate estimation accuracy. Notably, DoEuler-UKF requires a moderate amount of calculation, but its estimation accuracy is almost the same as USQUE. DoEuler-UKF produces good integrated navigation filtering performance. However, scholars have not paid enough attention to DoEuler-UKF.

5. Conclusions

We provide the following principles and recommendations for the use of integrated navigation filtering methods according to our research on the five filters.
First, USQUE has high estimation accuracy and high computational cost. For situations with high-accuracy estimate requirements, especially for aeronautics and astronautics, the calculation cost is secondary because high accuracy is required to meet the safety and stability requirements. Thus, it is generally recommended to use USQUE. MEKF has slightly lower accuracy than USQUE with a reasonable calculation cost. For nonlinear measurement functions, MEKF has higher accuracy through iterations with minimal calculation cost. Thus, MEKF is widely used for practical engineering. DoEuler-UKF has almost the same accuracy as USQUE, while having a moderate calculation cost, but it is not as stable as USQUE or MEKF. In some relatively stable systems, DoEuler-UKF will perform the best. MRP-UKF has moderate accuracy and high computational cost. Finally, in situations with slow motion and relatively low accuracy requirements, Euler-KF would be suitable for application because the filtering function of Euler-KF is linear and easy to implement and understand.

Author Contributions

Data curation, K.L.; Resources, M.Z.; Software, M.Z., K.L., B.H. and C.M.; Writing—original draft, M.Z.; Writing—review & editing, M.Z. and K.L.

Funding

This research was funded by the National Natural Science Foundation of China (No. 61703419, No.61873275, No.41474061).

Acknowledgments

The author would like to thank. Crassidis and Junkins for providing the MATLAB programs for the examples of Fangjun Qin [30] in http://www.buffalo.edu/johnc/estim_book.htm. These programs have given the authors a more thorough understanding of integrated navigation filter algorithms and accelerated the corresponding investigation in this paper.

Conflicts of Interest

The authors declare that there is no conflict of interest.

References

  1. Alandry, B.; Latorre, L.; Mailly, F.; Nouet, P. A fully integrated inertial measurement unit: Application to attitude and heading determination. IEEE Sens. J. 2011, 11, 2852–2860. [Google Scholar] [CrossRef]
  2. Zhang, L.; Chen, M.; He, H. A method research on robust fault diagnosis of integrated navigation systems. In Proceedings of the 8th International Conference on Electronic Measurement and Instruments (IEEE), Xi’an, China, 16–18 August 2007; pp. 3-326–3-330. [Google Scholar]
  3. Li, K.; Chang, L.; Hu, B. Unscented attitude estimator based on dual attitude representations. IEEE Trans. Instrum. Meas. 2015, 64, 3564–3576. [Google Scholar] [CrossRef]
  4. Wu, Y.X.; Pan, X.F. Velocity/position integration formula, Part I: Application to in-flight coarse alignment. IEEE Trans. Aerosp. Electron. Syst. 2013, 49, 1006–1023. [Google Scholar] [CrossRef]
  5. Wang, D.; Xu, X.; Zhu, Y. A Novel Hybrid of a Fading Filter and an Extreme Learning Machine for GPS/INS during GPS Outages. Sensors 2018, 18, 3863. [Google Scholar] [CrossRef]
  6. Jalving, B.; Gade, K.; Svartveit, K.; Willumsen, A.B.; Sorhagen, R. DVL Velocity Aiding in the HUGIN 1000 Integrated Navigation System; Norwegian Defense Research Establishment: Kjeller, Norway, 2004. [Google Scholar]
  7. Qin, Y.; Zhang, H. Kalman Filter and Integrated Navigation; Northwestern Polytechnical University Press: Xi’an, China, 2015. [Google Scholar]
  8. Gu, D.Q.; El-Sheimy, N.; Hassan, T.; Syed, Z. Coarse Alignment for Marine SINS Using Gravity in the Inertial Frame as a Reference. In Proceedings of the IEEE/ION Position, Location and Navigation Symposium (PLANS), Monterey, CA, USA, 5–8 May 2008; pp. 961–965. [Google Scholar]
  9. Ren, H.; Kazanzides, P. Investigation of attitude tracking using an integrated inertial and magnetic navigation system for hand-held surgical instruments. IEEE/ASME Trans. Mechatron. 2012, 17, 210–217. [Google Scholar] [CrossRef]
  10. Levinson, E.; Horst, J. The next generation marine inertial navigator is here now. In Proceedings of the IEEE Position Location and Navigation Symposium, Las Vegas, NV, USA, 11–15 April 1994; pp. 121–127. [Google Scholar]
  11. Qin, Y. Inertial Navigation; Science Press: Beijing, China, 2005. [Google Scholar]
  12. Chen, Y. Principle of Inertial Navigation; National Defense Industry Press: Beijing, China, 2007. [Google Scholar]
  13. Xu, R.; Ding, M.; Qi, Y.; Yue, S. Performance Analysis of GNSS/INS Loosely Coupled Integration Systems under Spoofing Attacks. Sensors 2018, 18, 4108. [Google Scholar] [CrossRef] [PubMed]
  14. Crassidis, J.L.; Markley, F.; Cheng, Y. Survey of Nonlinear Attitude Estimation Methods. J. Guid. Control Dyn. 2007, 30, 12–28. [Google Scholar] [CrossRef]
  15. Dissanayake, G.; Sukkarieh, S.; Nebot, E. The aiding of a low-cost strapdown inertial measurement unit using vehicle model constraints for land vehicle applications. IEEE Trans. Robot. Autom. 2001, 17, 731–747. [Google Scholar] [CrossRef]
  16. Sun, F.; Xia, J.Z.; Lan, H.Y.; Zhang, Y. DVL-aided Paralled Algorithm for Marine Attitude and Heading Reference System. J. Comput. Inf. Syst. 2014, 10, 1019–1027. [Google Scholar]
  17. Liu, Y.; Liu, F.; Gao, Y.; Zhao, L. Implementation and Analysis of Tightly Coupled Global Navigation Satellite System Precise Point Position/Inertial Navigation System with Insufficient Satellites for Land Vehicle Navigation. Sensors 2018, 18, 4305. [Google Scholar] [CrossRef]
  18. Chang, L.; Li, Y.; Xue, B. Initial Alignment for Doppler Velocity Log aided Strapdown Inertial Navigation System with Limited Information. IEEE/ASME Trans. Mechatron. 2016. [Google Scholar] [CrossRef]
  19. Groves, P.D. Principles of GNSS, Inertial, and Multisensor Integrated Navigation Systems; Artech House: Norwood, MA, USA, 2008. [Google Scholar]
  20. Wiener, N. Extrapolation, Interpolation, and Smoothing of Stationary Time Series; MIT Press: Cambridge, MA, USA, 1949. [Google Scholar]
  21. Farrell, J.; Barth, M. The Global Positioning System and Inertial Navigation; McGraw-Hill: New York, NY, USA, 1998. [Google Scholar]
  22. Wang, Y.; Sun, S.; Li, L. Adaptively robust unscented Kalman filter for tracking a maneuvering vehicle. J. Guid. Control Dyn. 2014, 37, 1697–1701. [Google Scholar] [CrossRef]
  23. Feng, K.; Li, J.; Zhang, X.; Zhang, X.; Shen, C.; Cao, H. An Improved Strong Tracking Cubature Kalman Filter for GPS/INS Integrated Navigation Systems. Sensors 2018, 18, 1919. [Google Scholar] [CrossRef] [PubMed]
  24. Zhou, J.; Knedlik, S.; Loffeld, O. INS/GPS Tightly-coupled Integration using Adaptive Unscented Particle Filter. J. Navig. 2010, 63, 491–511. [Google Scholar] [CrossRef]
  25. Jiao, Z.; Zhang, R. Kalman/Particle Filter for Integrated Navigation System. Adv. Mater. Res. 2013, 756, 2142–2146. [Google Scholar] [CrossRef]
  26. Hao, Y.; Xu, A. A Modified Extended Kalman Filter for a Two-Antenna GPS/INS Vehicular Navigation System. Sensors 2018, 18, 3809. [Google Scholar] [CrossRef] [PubMed]
  27. Costanzi, R.; Fanelli, F.; Monni, N.; Ridolfi, A.; Allotta, B. An attitude estimation algorithm for mobile robots under unknown magnetic disturbances. IEEE/ASME Trans. Mechatron. 2016, 21, 1900–1911. [Google Scholar] [CrossRef]
  28. Zhang, L.J. Multiplicative filtering for spacecraft attitude determination. J. Natl. Univ. Def. Technol. 2013, 35, 41–48. [Google Scholar]
  29. Jia, R. Attitude Estimation Algorithm for Low Cost MEMS Based on Quaternion EKF. Chin. J. Sens. Actuators 2014, 27, 90–95. [Google Scholar]
  30. Qin, F.; Chang, L.; Jiang, S.; Zha, F. A Sequential Multiplicative Extended Kalman Filter for Attitude Estimation Using Vector Observations. Sensors 2018, 18, 1414. [Google Scholar] [CrossRef]
  31. Crassidis, J.L.; Markley, F. Unscented Filtering for Spacecraft Attitude Estimation. J. Guid. Control Dyn. 2003, 26, 536–542. [Google Scholar] [CrossRef]
  32. 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]
  33. Zhou, J.; Edwan, E.; Kinedlik, S.; Loffeld, O. Low-cost INS/GPS with nonlinear filtering methods. In Proceedings of the 13th IEEE Conference on Information Fusion, Edinburgh, UK, 26–29 July 2010; pp. 1–8. [Google Scholar]
  34. Li, K.; Hu, B.; Chang, L. Modified Quaternion Unscented Kalman filter based on Modified Rodriguez Parameters. J. Navig. 2015. submitted. [Google Scholar]
  35. Karlgaard, C.D.; Schaub, H. Nonsingluar attitude filtering using modified Rodrigues parameters. J. Astronaut. Sci. 2010, 4, 777–791. [Google Scholar]
  36. Ran, C.; Cheng, X. A Direct and Non-singular UKF Approach using Euler angle kinematics for integrated navigation system. Sensors 2016, 16, 1415. [Google Scholar] [CrossRef] [PubMed]
  37. Markley, F.L.; Crassidis, J.L. Fundamentals of Spacecraft Attitude Determination and Control; Springer: New York, NY, USA, 2014. [Google Scholar]
Figure 1. The relationship between attitude, model, and filtering methods.
Figure 1. The relationship between attitude, model, and filtering methods.
Sensors 19 01426 g001
Figure 2. USQUE algorithm.
Figure 2. USQUE algorithm.
Sensors 19 01426 g002
Figure 3. (a)The forward rotate order of Euler angle; (b) the inverted rotate order of Euler angle.
Figure 3. (a)The forward rotate order of Euler angle; (b) the inverted rotate order of Euler angle.
Sensors 19 01426 g003
Figure 4. The simulation test trajectory.
Figure 4. The simulation test trajectory.
Sensors 19 01426 g004
Figure 5. (a) Pitch errors of each filter in ME2; (b) Roll errors of each filter in ME2; (c) Yaw errors of each filter in ME2.
Figure 5. (a) Pitch errors of each filter in ME2; (b) Roll errors of each filter in ME2; (c) Yaw errors of each filter in ME2.
Sensors 19 01426 g005
Figure 6. (a) Latitude errors of each filter; (b) Longitude errors of each filter; (c) Height errors of each filter.
Figure 6. (a) Latitude errors of each filter; (b) Longitude errors of each filter; (c) Height errors of each filter.
Sensors 19 01426 g006
Figure 7. (a) East velocity errors of each filter; (b) North velocity errors of each filter; (c) Up velocity errors of each filter.
Figure 7. (a) East velocity errors of each filter; (b) North velocity errors of each filter; (c) Up velocity errors of each filter.
Sensors 19 01426 g007
Figure 8. (a) Pitch estimate of each filter; (b) Roll estimate of each filter; (c) Yaw estimate of each filter.
Figure 8. (a) Pitch estimate of each filter; (b) Roll estimate of each filter; (c) Yaw estimate of each filter.
Sensors 19 01426 g008aSensors 19 01426 g008b
Figure 9. (a) East gyro bias estimate of each filter; (b) North gyro bias estimate of each filter; (c) Up gyro bias estimate of each filter.
Figure 9. (a) East gyro bias estimate of each filter; (b) North gyro bias estimate of each filter; (c) Up gyro bias estimate of each filter.
Sensors 19 01426 g009aSensors 19 01426 g009b
Table 1. Calculation equations of Jacobian matrix.
Table 1. Calculation equations of Jacobian matrix.
F 11 = ( ω ˜ i b b ε ^ ) × F 12 = C n b ( q ^ ) ω i n n p p ^ , v ^ F 13 = C n b ( q ^ ) ω i n n v p ^ F 14 = I 3 × 3 0 3 × 3
0 3 × 3 F 22 = p ˙ p p ^ , v ^ F 23 = p ˙ v p ^ 0 3 × 3 0 3 × 3
F 31 = C b n ( q ^ ) f b × F 32 = v ˙ p p ^ , v ^ F 33 = v ˙ v p ^ , v ^ 0 3 × 3 F 35 = C b n ( q ^ )
Table 2. Calculation equations of velocity partial differentials.
Table 2. Calculation equations of velocity partial differentials.
Y 11 = v E v N sec 2 λ R N + h v E v N tan λ ( R N + h ) 2 R N λ + 2 ω i e e v N cos λ + v E v U ( R N + h ) 2 R N λ + 2 ω i e e v U sin λ
Y 13 = v E v U v E v N tan λ ( R N + h ) 2
Y 21 = v E 2 sec 2 λ R N + h + v E 2 tan λ ( R N + h ) 2 R N λ 2 ω i e e v E cos λ + v E v U ( R M + h ) 2 R M λ
Y 23 = v E 2 tan λ ( R N + h ) 2 + v N v U ( R M + h ) 2
Y 31 = v E 2 ( R N + h ) 2 R N λ v N 2 ( R M + h ) 2 R M λ 2 ω i e e v E sin λ g λ
Y 33 = v E 2 ( R N + h ) 2 v N 2 ( R M + h ) 2 g h
Table 3. Calculation equations of velocity partial differentials.
Table 3. Calculation equations of velocity partial differentials.
Z 11 = v U + v N tan λ R N + h Z 12 = v E tan λ R N + h + 2 ω i e e sin λ Z 13 = v E R N + h 2 ω i e e cos λ
Z 21 = 2 v E tan λ R N + h + 2 ω i e sin λ Z 22 = v U R M + h Z 23 = v N R M + h
Z 31 = 2 v E R N + h + 2 ω i e e cos λ Z 32 = 2 v N R M + h 0
Table 4. Initialization conditions.
Table 4. Initialization conditions.
FilterInitialization StateFiltering CovarianceSystem Noise CovarianceMeasurement Noise Covariance
USQUE/MEKF x ^ 0 ( q 0 ) P 0 Q 0 R 0
MRP-UKF x ^ 0 ( σ 0 ) , σ 0 = ρ / ( 1 + q 4 )
DoEuler-UKF x ^ 0 ( φ 0 ) , φ 0 q 0
Euler-KF x ^ 0 ( 0 )
Table 5. Average estimation error and computation time of each filter.
Table 5. Average estimation error and computation time of each filter.
Average Estimation ErrorEuler-KFMEKFUSQUEMRP-UKFDoEuler-UKF
Pitch error (arcmin)1.2310.5190.2780.4910.354
Roll error (arcmin)1.0690.5290.3780.5640.457
Yaw error (arcmin)24.06418.5639.62511.6889.969
Latitude error (m)0.9820.1560.0980.1430.132
Longitude error (m)0.2450.1770.1540.1670.162
Height error (m)11.1670.4750.3980.4010.399
East velocity error (m/s)0.004620.006520.005180.005310.00523
North velocity error (m/s)0.02540.00720.00580.00920.0123
Up velocity error (m/s)0.1670.0120.0080.0080.009
Computation time (s)3.4910.2633.6526.7418.65
Table 6. Specifications of the XW-IMU5220.
Table 6. Specifications of the XW-IMU5220.
Technical IndexesGyroscopeAccelerometer
Dynamic range ± 150   ° / s ± 10   g
Bias 0.5   ° / s 0.005   g
Bias stability 0.02   ° / s 0.001   g
Bias repeatability 0.02   ° / s 0.002   g
Update rate100 Hz100 Hz
Table 7. The initialization conditions.
Table 7. The initialization conditions.
FilterStateInitialization
MEKF q , v N , v E , λ , ϕ , ε q 0
v N , v E and λ , ϕ provided by laser SINS
ε = 0 3 × 1
USQUE q , v N , v E , λ , ϕ , ε q 0 with x δ = 0 3 × 1
v N , v E and λ , ϕ provided by laser SINS
ε = 0 3 × 1
MRP-UKF σ , v N , v E , λ , ϕ , ε σ 0
v N , v E and λ , ϕ provided by laser SINS
ε = 0 3 × 1
DoEuler-UKF θ , γ , ψ , θ r , γ r , ψ r , v N , v E , λ , ϕ , ε θ 0 , γ 0 , ψ 0 / θ r 0 , γ r 0 , ψ r 0
v N , v E and λ , ϕ provided by laser SINS
ε = 0 3 × 1
Euler-KF φ , δ v N , δ v E , δ λ , δ ϕ , δ ε 0 10 × 1
Table 8. Average estimation error and computation time of each filter.
Table 8. Average estimation error and computation time of each filter.
Average Estimation ErrorEuler-KFMEKFUSQUEMRP-UKFDo-Euler UKF
Pitch error (°)4.8710.3210.2391.1790.284
Roll error (°)1.8430.5920.2150.4850.424
Yaw error (°)7.4255.6962.2143.3594.218
East gyro bias error (°/s)0.00120.00140.00130.00140.0012
North gyro bias error (°/s)8.69 × 10−48.82 × 10−47.36 × 10−48.66 × 10−48.51 × 10−4
Up gyro bias error (°/s)0.00190.00220.00180.00140.0024
Computation time (s)5.50719.85961.35555.21932.687
Table 9. Indicators of integrated navigation filters.
Table 9. Indicators of integrated navigation filters.
Integrated Filtering MethodCore Equation of Filtering MethodAmount of CalculationEstimation AccuracyCharacteristic
Euler-KFLinearEquation (10)LeastWorstMature
MEKFNonlinear Table 1Fairly smallModerateWidely used in engineering
USQUENonlinear Table 2HugeBestResearch hotspots
MRP-UKFNonlinear Equations (29)–(31)largeModerateNeed improvement
DoEuler-UKFNonlinear Equations (32)–(33)ModerateBetterPotential

Share and Cite

MDPI and ACS Style

Zhang, M.; Li, K.; Hu, B.; Meng, C. Comparison of Kalman Filters for Inertial Integrated Navigation. Sensors 2019, 19, 1426. https://doi.org/10.3390/s19061426

AMA Style

Zhang M, Li K, Hu B, Meng C. Comparison of Kalman Filters for Inertial Integrated Navigation. Sensors. 2019; 19(6):1426. https://doi.org/10.3390/s19061426

Chicago/Turabian Style

Zhang, Mengde, Kailong Li, Baiqing Hu, and Chunjian Meng. 2019. "Comparison of Kalman Filters for Inertial Integrated Navigation" Sensors 19, no. 6: 1426. https://doi.org/10.3390/s19061426

APA Style

Zhang, M., Li, K., Hu, B., & Meng, C. (2019). Comparison of Kalman Filters for Inertial Integrated Navigation. Sensors, 19(6), 1426. https://doi.org/10.3390/s19061426

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