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.
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 , 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 and covariance . is quaternion attitude and is the other states . The discrete time state inertial integrated navigation process model can be given as .
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
The attitude error matrix is given by
For small angles the vector part of the quaternion is approximately equal to half angles, then
and
, where
is a small error-angle rotation vector. The linearized model error-kinematics is as follows:
where
and
.
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:
Then,
where
and
. The partials are shown as
where
and
. 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
. The state, state error, noise vector, and covariance used in MEKF are given by
where
is a 3 × 3 zeros matrix. After propagating filtering state, MEKF propagates covariance using an error equation. Then, the error equation is given by
where
The system noise is
and covariance is
and the noise coefficient is
with
Table 1.
The position partial differentials are
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,
and
are separately expressed in the measurement update. We introduce the measurement update by taking the loosely-coupled model as an example.
where
varies with the change in the measurement. When only using position information as the loosely-coupled measurement transfer matrix,
When only using velocity information as the loosely-coupled measurement transfer matrix:
When using both position and velocity information as the loosely-coupled measurement transfer matrix:
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 and . |
Propagate filtering state and covariance |
|
|
where is the discrete-time state transition matrix obtained via numerical solution. |
Measurement update: |
Compute filtering gain |
Compute state error with |
Compute quaternion states |
|
Compute other states |
|
Compute filtering covariance |
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)
, where the state vector is
. The corresponding filter covariance matrix is
. The USQUE algorithm is used to estimate
and then the corresponding quaternion attitude
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 |
Using a quaternion describe attitude and covariance |
Create sigma points using UT |
The corresponding quaternion error |
|
|
Compute the quaternion sigma points , then |
Propagate sigma points |
Compute error quaternion sigma points |
and |
Compute GRP sigma points |
and |
The predicted mean and covariance of state |
|
|
Measurement update: |
Compute filtering gain, |
|
Compute posterior mean and covariance, |
and |
with is an n-dimension unit matrix |
Attitude update: |
The corresponding quaternion is |
|
with |
and |
The update of quaternion is given by and , |
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:
where the MRP
when N = 2, and
is the attitude rotation vector. Thus, we know that the singularity appears at
where
n is a positive integer. Traditionally, the MRP can be transferred to SMRP
to avoid singularity. The relationship between MRP and SMRP is:
Since the SMRP is defined by
, 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:
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:
where
is the inverted Euler angle. Then, the corresponding attitude kinematic equation is given by:
The singularity of the forward Euler angle occurs when pitch and the singularity of the inverted Euler angle occurs when roll . 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.
ME2 is defined as:
where
is the reference state,
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.
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:
where
is the average estimation error,
is the number of total Monte-Carlo simulations,
is the total simulation time, and
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.
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:
where
is the average estimation error,
is the total experience time, and
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 , 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.