Next Article in Journal
Motion Capture Technology in Sports Scenarios: A Survey
Next Article in Special Issue
A Machine Learning-Based Tropospheric Prediction Approach for High-Precision Real-Time GNSS Positioning
Previous Article in Journal
Gait Pattern Analysis: Integration of a Highly Sensitive Flexible Pressure Sensor on a Wireless Instrumented Insole
Previous Article in Special Issue
Modified RTK-GNSS for Challenging Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Improved Initial Alignment Method Based on SE2(3)/EKF for SINS/GNSS Integrated Navigation System with Large Misalignment Angles

1
College of Internet of Things, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
2
State Key Laboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai 200240, China
3
Portland Institute, Nanjing University of Posts and Telecommunications, Nanjing 210003, China
4
Key Laboratory of Modern Agricultural Equipment and Technology, Ministry of Education and Jiangsu Province, Jiangsu University, Zhenjiang 212013, China
5
School of Agricultural Engineering, Jiangsu University, Zhenjiang 212013, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(9), 2945; https://doi.org/10.3390/s24092945
Submission received: 25 March 2024 / Revised: 17 April 2024 / Accepted: 2 May 2024 / Published: 6 May 2024
(This article belongs to the Special Issue GNSS Signals and Precise Point Positioning)

Abstract

:
This paper proposes an improved initial alignment method for a strap-down inertial navigation system/global navigation satellite system (SINS/GNSS) integrated navigation system with large misalignment angles. Its methodology is based on the three-dimensional special Euclidean group and extended Kalman filter (SE2(3)/EKF) and aims to overcome the challenges of achieving fast alignment under large misalignment angles using traditional methods. To accurately characterize the state errors of attitude, velocity, and position, these elements are constructed as elements of a Lie group. The nonlinear error on the Lie group can then be well quantified. Additionally, a group vector mixed error model is developed, taking into account the zero bias errors of gyroscopes and accelerometers. Using this new error definition, a GNSS-assisted SINS dynamic initial alignment algorithm is derived, which is based on the invariance of velocity and position measurements. Simulation experiments demonstrate that the alignment method based on SE2(3)/EKF can achieve a higher accuracy in various scenarios with large misalignment angles, while the attitude error can be rapidly reduced to a lower level.

1. Introduction

The initial alignment is a crucial component of the strap-down inertial navigation system (SINS) and is assessed based on its speed and accuracy [1]. Currently, there are two main categories of traditional alignment methods based on misalignment angles—small misalignment angles linear alignment and large misalignment angles nonlinear alignment [2]. Linear alignment models and linear filtering algorithms for small misalignment angles are well established, while research on the alignment problems for large misalignment angles has mainly focused on nonlinear models and nonlinear filtering algorithms [3]. However, these approaches can lead to errors in model linearization, increased computational complexity, and reduced filtering accuracy for large misalignment angles. In the 21st century, modern warfare has demonstrated the significance of high-tech conditions, where precision is a key consideration, alongside speed. Consequently, it is essential to find ways to achieve fast and highly accurate initial alignment under large misalignment angles, for applications like weapon launchers and guided weapons that require emergency mobile transfer. Addressing this urgent problem is crucial from both a system performance perspective in modern warfare and from practical alignment environments, data utilization, and engineering software development [4,5,6]. Hence, studying the linearized uniform initial alignment algorithm without coarse alignment at any misalignment angles is of paramount importance and offers significant theoretical and engineering benefits.
The current most popular mode of navigation is the SINS/Global navigation satellite system (GNSS) integrated navigation model [7,8]. Satellite navigation offers both high accuracy and low cost. In the SINS/GNSS integrated navigation system, GNSS not only improves navigation accuracy by correcting inertial navigation information, but also provides initial position and velocity information for SINS [9]. However, the GNSS cannot directly provide attitude information. In recent years, nonlinear system state estimation has gained attention in the field of inertial navigation for combined navigation and initial alignment. This includes methods such as the extended Kalman filter (EKF) [10], unscented Kalman filter (UKF) [11,12], particle filter (PF) [13], cubature Kalman filter (CKF) [14], unscented particle filter (UPF) [15], distributed Kalman filter (DKF) [16], or a combination of multiple nonlinear filters [17,18]. Among these, the traditional EKF has two drawbacks. Firstly, it requires high accuracy in the initial value of the system state. If the initial value is based on the actual situation, the filter may struggle to converge. Secondly, the EKF can lead to inconsistency. When a new observation is received, the EKF calculates the covariance matrix of the current state based on the linearization of the previous state. However, the actual value of the covariance matrix may not align with this calculated value. Compared to the EKF, the UKF is more practical. However, both the traditional UKF and EKF require prior statistics on known system noise and measurement noise. In practical applications, the accuracy of the filter is inevitably affected by environmental limitations and algorithm problems, leading to a decrease in accuracy or even divergence.
In the traditional EKF framework used for the GNSS/SINS integrated navigation system, state errors such as position and velocity are defined only by considering differences in size, completely ignoring differences in direction. This oversight can lead to inconsistent definitions of state error coordinate systems. In recent years, a new filter has been developed based on Lie groups and manifold theory. Its core idea is to define the states and errors such as attitude, velocity, and position in the group space, rather than in the traditional Euclidean space. According to the multiplicative closure property of Lie groups and the affine property between Lie algebras, a more compact dynamic equation of position and coordinate state errors, considering attitude errors, has been redesigned to meet the requirement that the error transfer matrix or measurement matrix remain unchanged or slowly change and to achieve the independence of F or H matrix and state estimation. Therefore, it can be unchanged, also known as invariant EKF. The error definition corresponding to the invariant EKF has good self-consistency, which effectively overcomes the defect that the traditional EKF is too dependent on the initial value, and has better convergence and consistency of filtering [19]. It is more rigorous, in theory, and has been well applied in inertial navigation fields such as robot attitude estimation, simultaneous localization and mapping (SLAM), and visual odometers (VOs) [20].To address this issue, this paper adopts the invariant EKF concept. It defines states and errors, such as attitude, velocity, and position, in group space. Based on this approach, new inertial navigation error equations and measurement update methods are deduced. The main contributions of this paper can be summarized as follows: (1) Constructing the attitude, velocity, and position states as three-dimensional special Euclidean group (SE2(3))/EKF elements, which takes into account the zero bias errors of gyroscopes and accelerometers. This facilitates the formation of a mixed group vector error. The alignment method based on SE2(3)/EKF filtering has demonstrated faster convergence speed, higher accuracy, and greater computational efficiency than the traditional EKF. (2) Attaining the accurate position of the vehicle during the alignment process is crucial for improving the alignment accuracy of the system. Moreover, the position alignment aids in directly transitioning the vehicle to the autonomous navigation stage after the GNSS-aided alignment. Consequently, it enhances the vehicle’s applicability, by eliminating the need to obtain the vehicle’s position again at the end of the alignment. Accurate and real-time position alignment can improve the attitude alignment accuracy of the system.
The paper follows this outline—Section 2 establishes the inertial navigation error equation under the Lie group framework. Section 3 presents the filtering model for the SINS/GNSS integrated navigation system based on SE2(3)/EKF. The simulation experiment was conducted in Section 4. Finally, Section 5 concludes.

2. Inertial Navigation Error Equation under the Lie Group Framework

SE2(3)/EKF does not simply use the difference between the estimated state and the real state as the error, but more carefully considers the consistency of the coordinate frame of the state definition and provides a more rigorous and compact mathematical form by redefining the error in Lie group space. In order to describe the state error, a lie group state error model including attitude, velocity, and position is designed and the gyro bias error and accelerometer bias error are still defined in the traditional vector space.
The related Cartesian reference coordinate frames used in this study are defined as follows [21]:
(1)
b-frame: Body coordinate system, with its three axes pointing to the right–front–up (R-F-U) of the carrier, respectively, denoted as x b y b z b ;
(2)
n-frame: It indicates the navigation frame and it is the frame used by the SINS to calculate the navigation parameters, denoted as Eastward–Northward–Upward (E-N-U);
(3)
e-frame: Earth coordinate system, with its origin at the geocenter. The x-axis is the intersection of the geocenter pointing to the prime meridian and the equator, the z-axis is the geocenter pointing to the north pole, and the y-axis forms a right-handed coordinate system with the x-axis and z-axis, denoted as x e y e z e ;
(4)
i-frame: Inertial coordinate system. It is a non-rotating coordinate system in inertial space, denoted as x i y i z i .
The schematic diagram of the coordinate system mentioned above is shown in Figure 1. ω i e is the angular velocity of the Earth’s rotation and L is the latitude of the SINS.

2.1. SINS Navigation Differential and Error Equations

Selecting the E-N-U geographic coordinate system as the navigation reference coordinate system for the SINS, the attitude differential equation using the n-frame as the reference frame is as follows [22]:
C ˙ b n = C b n ( ω i b b × ) [ ( ω i e n + ω e n n ) × ] C b n
where C b n denotes the attitude matrix of the b-frame relative to the n-frame, ω i b b is the angular velocity of the b-frame relative to the i-frame, ω i e n is the rotational angular velocity of the Earth, ω i e n = [ 0 ω i e cos L ω i e s i n L ] T , ω e n n is the angular velocity generated by the relative motion of the n-frame to the e-frame, ω e n n = [ v N R M + h v E R N + h v E R N + h tan L ] T , and h is the geographical altitude. R M is the principal radius of curvature along the meridional section, R N is the principal radius of curvature along the prime-vertical normal section, v n is the velocity in the n-frame, and v n = [ v E v N v U ] T . v E , v N , and v U are the eastern, northern, and upward velocity, respectively. × denotes the conversion of vectors into oblique symmetric matrices.
The differential equation for attitude error can be derived as follows [22]:
ϕ ˙ n = ϕ n × ω i n n + δ ω i n n C b n δ ω i b b
where ϕ is attitude error, ω i n n = ω i e n + ω e n n , δ ω i n n is the calculation error, δ ω i n n = δ ω i e n + δ ω e n n , δ ω i e n = [ 0 ω i e sin L δ L ω i e cos L δ L ] T = M 1 δ p n , M 1 = [ 0 0 0 ω i e sin L 0 0 ω i e cos L 0 0 ] , and δ ω i b b is the measurement error of the gyroscope.
δ ω e n n = [ δ v N R M + h + v N δ h ( R M + h ) 2 δ v E R N + h v E δ h ( R N + h ) 2 tan L δ v E R N + h + v E sec 2 L δ L R N + h v E tan L δ h ( R N + h ) 2 ] = M a v δ v n + M 2 δ p n
where M a v = [ 0 1 R M + h 0 1 R N + h 0 0 tan L R N + h 0 0 ] , M 2 = [ 0 0 v N ( R M + h ) 2 0 0 v E ( R N + h ) 2 v E sec 2 L R N + h 0 v E tan L ( R N + h ) 2 ] , p n = [ L λ h ] T , λ is the longitude, δ p n = [ δ L δ λ δ h ] T is the position error, and δ L , δ λ , and δ h are the latitude error, longitude error, and height errors, respectively. δ v n is the velocity error in the n-frame, δ v n = [ δ v E δ v N δ v U ] T , and δ v E , δ v N , and δ v U are the eastern, northern, and upward velocity errors, respectively.
The velocity differential equation defined under the local navigation system is as follows [18]:
v ˙ n = C b n f b ( 2 ω i e n + ω e n n ) × v n + g n
where f b is the specific force measured using the accelerometer, 2 ω i e n × v n is the Coriolis acceleration caused by the motion of the carrier and the rotation of the Earth, ω e n n × v n is the centripetal acceleration caused by the movement of the carrier towards the ground, g n is gravitational acceleration, and ( 2 ω i e n + ω e n n ) × v n + g n is collectively referred to as harmful acceleration.
The corresponding differential equation for velocity error is as follows [23]:
δ v ˙ n = f n × ϕ n + v n × ( 2 δ ω i e n + δ ω e n n ) ( 2 ω i e n + ω e n n ) × δ v n + C b n δ f b + δ g n
where f n is the projection of the specific force vector in the n-frame, δ f b is the measurement error of the accelerometer, and δ g n is the gravitational acceleration error.
The position differential equation defined under the local navigation system is as follows [23]:
p ˙ n = [ L ˙ λ ˙ h ˙ ] T = [ 0 1 R M + h 0 sec L R N + h 0 0 0 0 1 ] v n = M p v v n
where M p v = [ 0 1 R M + h 0 sec L R N + h 0 0 0 0 1 ] .
The corresponding differential equation for positional error is as follows [23]:
δ p ˙ n = [ δ L ˙ δ λ ˙ δ h ˙ ] T = [ 0 1 R M + h 0 sec L R N + h 0 0 0 0 1 ] δ v n + [ 0 0 v N ( R M + h ) 2 v E sec L tan L R N + h 0 v E sec L ( R N + h ) 2 0 0 0 ] δ p n = M p v δ v n + M p p δ p n
The error model for gyroscopes and accelerometers is as follows:
{ ω ^ i b b = ω i b b + ε b + w g f ^ b = f b + b + w a
where ω ^ i b b and f ^ b are actual measured values of gyroscopes and accelerometers, respectively. ε b is the random constant drift of the gyroscopes, b is the random constant bias of the accelerometers, and ε ˙ b = 0 , ˙ b = 0 . w g , and w a are the Gaussian white noises of the accelerometers and gyroscopes, respectively.

2.2. Left Invariant Error Equation of Inertial Navigation in the Lie Group Framework

In the Lie group framework, the n-frame is chosen as the projection coordinate system for the SINS/GNSS integrated system. The attitude, velocity, and position of the n-frame are defined as a group, χ , [24] as follows:
χ = [ C b n v n p n 0 1 × 3 1 0 0 1 × 3 0 1 ]
where C b n S O ( 3 ) , S O ( 3 ) is a three-dimensional special orthogonal group, 0 1 × 3 is a zero matrix with one row and three columns, and χ S E ( 3 ) is a more concise group that includes elements of attitude, velocity, and position.
By using the properties of Lie groups [25], we can obtain the following:
χ 1 = [ C n b v b p b 0 1 × 3 1 0 0 1 × 3 0 1 ]
where χ 1 denotes the inverse matrix of χ and, similarly, χ 1 S E ( 3 ) .
If χ represents the true state and χ ^ represents the nominal state, the error of the attitude, velocity, and position states is defined as η S E ( 3 ) . The error in Lie group space can be classified into two types—left invariance, η = χ ^ 1 χ , and right invariance, η = χ χ ^ 1 . Research has shown that left invariance can achieve the invariance or gradual change of the measurement matrix, while the observations of GNSS theoretically belong to the category of left invariant observations [25]; therefore, this paper focuses on the left invariant error for processing.
Based on Equations (1) and (2), the specific expression for the left invariant error is given by:
η = χ ^ 1 χ = [ C ^ n b v ^ b p ^ b 0 1 × 3 1 0 0 1 × 3 0 1 ] [ C b n v n p n 0 1 × 3 1 0 0 1 × 3 0 1 ] = [ C ^ n b C b n C ^ n b v n v ^ b C ^ n b p n p ^ b 0 1 × 3 1 0 0 1 × 3 0 1 ]
The errors corresponding to the attitude, velocity, and position under the defined Lie group are as follows:
C ^ n b C b n = E x p ( ϕ b ) = exp ( ϕ b × )
J ρ v b = C ^ n b v n v ^ b = C ^ n b v n C ^ n b v ^ n = C ^ n b δ v n
J ρ p b = C ^ n b p n p ^ b = C ^ n b p n C ^ n b p ^ n = C ^ n b δ p n
where ϕ b is the attitude error angles, E x p ( ) represents the mapping between Lie algebras and Lie groups, and exp ( ) represents exponentiation. exp ( ϕ b × ) = I 3 + sin θ θ ( ϕ b × ) + 1 cos θ θ 2 ( ϕ b × ) 2 and θ = | ϕ | . J is the Jacobian matrix of the Rodriguez formula and J = n = 0 1 ( n + 1 ) ! ( ϕ b × ) = I 3 + 1 cos θ θ 2 ( ϕ b × ) + θ sin θ θ 3 ( ϕ b × ) 2 . J ρ v b is the new definition of velocity error and J ρ p b is the new definition of positional error.
From the error form defined in Equation (11), we can observe that the velocity and position errors defined in the Lie group framework include attitude terms. This takes into account the differences in numerical magnitude and direction between true values and estimates. Therefore, the error definition is more reasonable and concise compared to the traditional vector space method of differencing, which solves the problem of inconsistent benchmarks for defining velocity error states.
According to the characteristics of Lie groups and Lie algebras, the left invariant error satisfies the following [26,27]:
η = [ exp ( ϕ n × ) J ρ v b J ρ p b 0 1 × 3 1 0 0 1 × 3 0 1 ]
The new attitude error equation is derived as follows:
d d t ( C ^ n b C b n ) = C ^ ˙ n b C b n + C ^ n b C ˙ b n = [ C ^ n b ( ω ^ i n n × ) ( ω ^ i b b × ) C ^ n b ] C b n + C ^ n b [ C b n ( ω i b b × ) ( ω i n n × ) C b n ] = C ^ n b ( ω ^ i n n × ) C b n ( ω ^ i b b × ) C ^ n b C b n + C ^ n b C b n ( ω i b b × ) C ^ n b ( ω i n n × ) C b n C ^ n b [ ( ω ^ i n n ω i n n ) × ] C b n [ ( ω i b b + δ ω i b b ) × ] [ I 3 + ( ϕ b × ) ] + [ I 3 + ( ϕ b × ) ] ( ω i b b × ) δ ω i n n × [ I 3 + ( ϕ b × ) ] + ( ϕ b × ω i b b ) × δ ω i b b × δ ω i n n × + ( ϕ b × ω i b b ) × δ ω i b b ×
where C ^ n b C b n I 3 + ( ϕ b × ) , δ ω i n n = ω ^ i n n ω i n n , δ ω i b b = ω ^ i b b ω i b b , and the second-order small quantities ( δ ω i n n × ) ( ϕ b × ) and ( δ ω i b b × ) ( ϕ b × ) are neglected in the process of formula derivation.
The new velocity error equation is derived as follows:
d d t ( J ρ v b ) = C ^ ˙ n b δ v n C ^ n b δ v ˙ n = [ C ^ n b ( ω ^ i n n × ) ( ω ^ i b b × ) C ^ n b ] δ v n C ^ n b [ f n × ϕ n + v n × ( 2 δ ω i e n + δ ω e n n ) ( 2 ω i e n + ω e n n ) × δ v n + C b n δ f b + δ g n ] ( ω ^ i b b × ) J ρ v b δ f b + ϕ b × f b C n b δ g n + ( C n b ω i e n ) × J ρ v b C n b ( v n × ) ( δ ω i e n + δ ω i n n )
where δ f b = f ^ b f b , δ ω i e n = ω ^ i e n ω i e n , and δ g n = g ^ n g n .
The new position error equation is derived as follows:
d d t ( J ρ p b ) = C ^ ˙ n b δ p n C ^ n b δ p ˙ n = [ C ^ n b ( ω ^ i n n × ) ( ω ^ i b b × ) C ^ n b ] δ p n C ^ n b ( M p v δ v n + M p p δ p n ) = ( ω ^ i b b × ) C ^ n b δ p n C ^ n b ( ω ^ i n n × ) δ p n C ^ n b M p v δ v n C ^ n b M p p δ p n ( ω ^ i b b × ) J ρ p b + J ρ v b C n b ( p n × ) δ ω e n n + ( C n b ω i e n ) × J ρ p b
By substituting Equations (13) and (14) into δ ω i e n and δ ω i n n , we can obtain the following equations:
δ ω i e n = M 1 δ p n = M 1 C ^ b n J ρ p b
δ ω i n n = ( M 1 + M 2 ) δ p n + M a v δ v n = ( M 1 + M 2 ) C ^ b n J ρ p b M a v C ^ b n J ρ v b
Finally, the error equations of SINS can be written as follows:
d d t ( C ^ n b C b n ) = C n b ( M 1 + M 2 ) C ^ b n J ρ p b C n b M a v C ^ b n J ρ v b ω i b b × ϕ b ( ε b + w g )
d d t ( J ρ v b ) = C n b ( v n × ) ( 2 M 1 + M 2 ) C b n J ρ p b + [ C n b ( v n × ) M a v C b n ( ω i b b × ) ( C n b ω i e n ) × ] J ρ v b f b × ϕ b C n b δ g n ( b + w a )
d d t ( J ρ p b ) = [ C n b ( p n × ) M 2 C b n ( ω i b b × ) + ( C n b ω i e n ) × ] J ρ p b + [ I 3 + C n b ( p n × ) M a v C b n ] J ρ v b

3. Filtering Model Based on SE2(3)/EKF

3.1. State Equation Based on SE2(3)

If the state vector is X and the system noise vector is W, then the state equation based on SE2(3) is as follows [28]:
X ˙ = F X + G W
where F is the state transition matrix, G is the system noise allocation matrix, and E [ W W T ] = Q .
X = [ ( ϕ b ) T ( J ρ v b ) T ( J ρ p b ) T ( ε b ) T ( b ) T ] T
F = [ ω ^ i b b × C ^ n b M a v C ^ b n C ^ n b ( M 1 + M 2 ) I 3 0 3 f ^ b × C ^ n b ( v ^ n × ) M a v C ^ b n ω ^ i b b × C ^ n b ( ω ^ i e n × ) C ^ b n C ^ n b ( v ^ n × ) ( 2 M 1 + M 2 ) 0 3 I 3 0 3 M p v C ^ b n M p v 0 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 ]
G = [ I 3 0 3 0 3 I 3 0 3 0 3 0 3 0 3 0 3 0 3 ]
W = [ w g w a ]
where I 3 is a 3 × 3 unit vector and 0 3 is a 3 × 3 zero matrix.

3.2. Measurement Equation Based on SE2(3)

The GNSS can provide position and velocity information in integrated navigation as measurements for Kalman filtering, to suppress the divergence of inertial navigation calculation results. Therefore, the measurement equation based on the SE2(3) can be written as follows [29]:
Z = [ v S I N S n v G N S S n p S I N S n p G N S S n ] = H X + V
where v S I N S n and v G N S S n are the velocity of the SINS and the GNSS in the n-frame, respectively. p S I N S n and p G N S S n are the position of the SINS and the GNSS in the n-frame, respectively. H is the measurement matrix and V is the white noise for velocity measurement and position measurement of satellite receivers.
H = [ v ^ n I 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 0 3 ]
V = [ V v V p ]
where V v and V p are the white noise for satellite receiver velocity measurement and position measurement, respectively.

3.3. SE2(3)/EKF Algorithm

Assuming the sampling interval of IMU is Δ t and the discretized state transition matrix is Φ , then the SE2(3)/EKF algorithm is as follows [30]:
(1)
One-step state prediction
X ^ k , k 1 = Φ k , k 1 X ^ k 1
(2)
State estimation
X ^ k = X ^ k , k 1 + δ X ^ k
(3)
State estimation error
δ X ^ k = K k ( Z k H k X ^ k , k 1 )
(4)
Filter gain
K k = P k , k 1 H k T [ H k P k , k 1 H k T + R k ] 1
(5)
One-step prediction mean square error
P k , k 1 = Φ k , k 1 P k 1 Φ k , k 1 T + Q k 1
(6)
Estimating mean square error
P k = ( I K k H k ) P k , k 1 ( I K k H k ) T + K k R k K k T
where X ^ k , k 1 is the state one-step prediction matrix, Φ k , k 1 is the one-step transition matrix from time t k 1 to time t k , and Φ k , k 1 = I + F Δ t . X ^ k 1 is the state estimation matrix at time t k 1 , δ X ^ k is the state estimation error matrix at time t k , K k is the filtering gain matrix at time t k , Z k is the measurement matrix at time t k , P k , k 1 is the one step prediction mean square error matrix, R k is the measurement noise variance matrix, Q k 1 is the system noise variance matrix at time t k 1 , Q k 1 = G k 1 ( q k 1 Δ t ) G k 1 T , G k 1 is the discretized system noise allocation matrix at time t k 1 , and q is the variance intensity matrix corresponding to the system noise matrix W . P k is the estimated mean square error matrix.

4. Simulation Results

This section presents the simulation experiments conducted to verify the validity and feasibility of the proposed method. The specific parameters of the inertial sensor used in the simulation experiments are as follows: the constant drift and random walk of the gyroscopes are represented by values 0.03 / h and 0.001 / h , respectively. The constant bias and random bias of the accelerometers are represented by values 100   μ g and 10   μ g , respectively. The data output frequency is 200   H z and the calculation period of SINS is 5   m s . The performance parameters of the GNSS are as follows: the velocity measurement accuracy is represented by a value of 0.1   m / s , the position measurement accuracy is represented by a value of 1   m , and the data output frequency is represented by a value of 1   H z . The initial position is set as follows: 118.786365   E , 32.057313   N , and the height is 0   m . The initial velocity is represented by a value of [ 0 0 0 ] T   m / s . The simulation time is 900 s. The state and trajectory of the vehicle are illustrated in Figure 2 and Figure 3, respectively.
To fully demonstrate the performance of the proposed algorithm, contrast experiments of SE2(3)/EKF and the EKF algorithm, as well as position alignment experiments, will be designed. These experiments will compare the traditional EKF with the SE2(3)/EKF algorithm proposed in this paper, using three different initial misalignment angle scenarios. The performance of various algorithms will be evaluated to verify the advantages of the SE2(3)/EKF algorithm.

4.1. Experiment 1

In this experiment, the misalignment angles in the three directions are sequentially set [ 1 1 40 ] T . The mean and variance of the attitude error and position error are recorded in Table 1 for the 50 s leading up to the end of alignment. The attitude angle error are as shown in Figure 4, Figure 5 and Figure 6, and the position alignment error are as shown in Figure 7, Figure 8 and Figure 9.
Based on the comprehensive information from the chart, both methods can converge the attitude angles to a relatively stable range in this scenario. The figure shows that the time taken to converge the horizontal attitude angles estimated by the two methods is not significantly different; both methods converge quickly. However, the convergence time for the heading angles is longer. The data provided in Table 1 indicate that the attitude angle accuracy of the two methods after alignment is relatively close. Overall, under the misalignment angles set in Experiment 1, both methods perform equally. This is because the initial misalignment angles are small at this time and the nonlinearity of the traditional EKF model is not severe. Additionally, defining the velocity error as the direct difference between the true value and the estimated value is considered reasonable, allowing both alignment methods to converge quickly.

4.2. Experiment 2

In this experiment, the misalignment angles in the three directions are sequentially set [ 1 1 95 ] T . The mean and variance of the attitude error and position error are recorded in Table 2 for the 50 s leading up to the end of alignment. The attitude angle error are as shown in Figure 10, Figure 11 and Figure 12, and the position alignment error are as shown in Figure 13, Figure 14 and Figure 15.
From the figures, it can be observed that, due to further increasing the initial heading angle error, the nonlinearity of the traditional model is enhanced. At this point, the EKF method still uses a linear error differential equation error model, disregarding the influence of high-order system terms. As a result, significant model errors occur, leading to a slower convergence speed in attitude compared to the SE2(3)/EKF method. According to Table 2, the convergence accuracy of the EKF method is lower than that of SE2(3)/EKF. The SE2(3)/EKF method redefines errors more rigorously within the Lie group framework, resulting in a better performance for larger misalignment angles.

4.3. Experiment 3

In this experiment, the misalignment angles in the three directions are sequentially set [ 20 20 165 ] T . The mean and variance of the attitude error and position error are recorded in Table 3 for the 50 s leading up to the end of alignment. The attitude angle error are as shown in Figure 16, Figure 17 and Figure 18, and the position alignment error are as shown in Figure 19, Figure 20 and Figure 21.
The figure shows that, when the misalignment angles in three directions reach a severe level, the nonlinearity of the model becomes extremely severe. Additionally, defining the velocity error as a direct subtraction becomes even more unreasonable. Consequently, the traditional EKF algorithm is no longer able to converge, while the SE2(3)/EKF method still converges relatively quickly. According to Table 3, it is evident that, in this scenario, the attitude angle error calculated by SE2(3)/EKF can still maintain a level close to Experiment 1 and Experiment 2, further demonstrating the stable characteristics of the SE2(3)/EKF alignment method.

5. Conclusions

This paper proposed a novel initial alignment method for the SINS/GNSS integrated navigation system with large misalignment angles based on SE2(3)/EKF. The left invariant SE2(3)/EKF adopts a group-vector mixed error model, which allows the linear state error based on Lie algebra to effectively capture the nonlinear error on the Lie group. This compatibility with the linear assumption of EKF filtering results in a higher precision dynamic model, improved measurement update accuracy, and better performance in dynamic initial alignment, especially when considering attitude error. Simulation experiments demonstrated that the EKF method can achieve attitude convergence in scenarios with small initial misalignment angles. On the other hand, the alignment method based on SE2(3)/EKF can converge to higher accuracy in the different large misalignment angle scenarios simulated in this paper, quickly reducing attitude errors to a lower level.
In future research, we focus on studying a more reasonable difference resistance adaptive filtering method within the framework of SE2(3). This will enable us to adapt to complex GNSS observation environments or high dynamic scenes of carriers, expanding the application scenario of this algorithm. Additionally, we aim to conduct real-time vehicle experiments to validate the proposed methodology for practical use.

Author Contributions

Conceptualization, J.S.; methodology, J.S. and Y.C.; software, J.S. and B.C.; validation, J.S., Y.C. and B.C.; formal analysis, J.S.; investigation, J.S.; writing—original draft preparation, J.S.; writing—review and editing, J.S.; funding acquisition, J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by the National Nature Science Foundation of China under Grant 62203231, the State Key Laboratory of Ocean Engineering (Shanghai Jiao Tong University) under Grant GKZD010084, the China Postdoctoral Science Foundation under Grant 2020M681685, the Postdoctoral Research Funding Project of Jiangsu Province under Grant 2021K161B, the Jiangsu Province and Education Ministry Co-sponsored Synergistic Innovation Center of Modern Agricultural Equipment, grant number XTCX2009, Open Funding from the Key Laboratory of Modern Agricultural Equipment and Technology(Jiangsu University), Ministry of Education, grant number MAET202301.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The abbreviations used in this manuscript are as follows:
SINSStrap-down Inertial Navigation System
GNSSGlobal Navigation Satellite System
SESet Euclidean
EKFExtended Kalman Filter
UKFUnscented Kalman Filter
PFParticle Filter
CKFCubature Kalman Filter
UPFUnscented Particle Filter

References

  1. Fang, J.C.; Wan, D.J. A Fast Initial Alignment Method for Strap-down Inertial Navigation System on Stationary Base. IEEE Trans. Aerosp. Electron. Syst. 1996, 32, 1501–1504. [Google Scholar] [CrossRef]
  2. Li, J.C.; Gao, W.; Zhang, Y.; Wang, Z.C. Gradient Descent Optimization-based Self-alignment Method for Stationary SINS. IEEE Trans. Instrum. Meas. 2019, 68, 3278–3286. [Google Scholar] [CrossRef]
  3. Zhang, B.; Wang, X.D.; Lu, H.; Hao, Z.J.; Gu, C.C. Application of Adaptive Robust CKF in SINS/GPS Initial Alignment with Large Azimuth Misalignment Angles. Math. Prob. Eng. 2021, 2021, 1–6. [Google Scholar]
  4. Tang, C.Y.; Chen, L.; Chen, J.F. Improvement to classical coning algorithms in maneuver performance based on a match correction structure. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2019, 233, 2478–2488. [Google Scholar] [CrossRef]
  5. Tang, C.Y. Non-initial attitude error analysis of SINS by error closed-loop and energy-level constructing. Aerosp. Sci. Technol. 2020, 106, 106149. [Google Scholar] [CrossRef]
  6. Li, Y.C.; Cai, Y.F.; Li, Z.; Feng, S.Z.; Wang, H.; Sotelo, M.A. Map-based localization for intelligent vehicles from bi-sensor data fusion. Expert Syst. Appl. 2022, 203, 117586. [Google Scholar] [CrossRef]
  7. Zhu, C.Y.; Liu, X.Q.; Chen, H.L.; Tian, X. Automatic cruise system for water quality monitoring. Int. J. Agric. Biol. Eng. 2018, 11, 244–250. [Google Scholar] [CrossRef]
  8. Ruan, C.Z.; Zhao, D.; Sun, Y.P.; Hong, J.Q.; Ding, S.H.; Luo, J. Design and testing of a control system associated with the automatic feeding boat for farming Chinese river crabs. Comput. Electron. Agric. 2018, 150, 14–25. [Google Scholar] [CrossRef]
  9. Cui, B.; Chen, W.; Weng, D.; Wei, X.; Sun, Z.; Zhao, Y.; Liu, Y. Observability-Constrained Resampling-Free Cubature Kalman Filter for GNSS/INS with Measurement Outliers. Remote Sens. 2023, 15, 4591. [Google Scholar] [CrossRef]
  10. Zhang, J.L.; Qin, Y.Y.; Mei, C.B. Application of EKF and SUKF in SINS Alignment on Great-Swing Base. Meas. Control Technol. 2013, 32, 8–11. [Google Scholar]
  11. Zhang, T.; Xu, X.S. Nonlinear Initial Alignment of INS based on Simplified UKF. J. Chin. Inerti. Technol. 2011, 19, 537–542. [Google Scholar]
  12. Cao, S.Y.; Gao, H.L.; You, J. In-Flight Alignment of Integrated SINS/GPS/Polarization/Geomagnetic Navigation System Based on Federal UKF. Sensors 2022, 22, 5985. [Google Scholar] [CrossRef] [PubMed]
  13. Xu, X.S.; Liu, X.Y. Improved Particle Filter Algorithm in SINS Initial Alignment. J. Chin. Inerti. Technol. 2016, 24, 299–305. [Google Scholar]
  14. Sun, F.; Tang, L.J. Initial Alignment of Large Azimuth Misalignment Angles in SINS Based on CKF. Chin. J. Sci. Instrum. 2012, 33, 327–333. [Google Scholar]
  15. Sun, J.; Xu, X.S.; Liu, Y.T.; Zhang, T.; Li, Y. Initial Alignment of Large Azimuth Misalignment Angles in SINS Based on Adaptive UPF. Sensors 2015, 15, 21807–21823. [Google Scholar] [CrossRef] [PubMed]
  16. Xu, Y.; Cao, J.; Shmaliy, Y.S.; Zhuang, Y. Distributed Kalman Filter for UWB/INS Integrated Pedestrian Localization under Colored Measurement Noise. Satell. Navig. 2021, 2, 22. [Google Scholar] [CrossRef]
  17. Gao, H.L.; You, J.; Cao, S.Y. In-Flight Alignment Method of Integrated SINS/GPS Navigation System Based on Combined PF-UKF Filter. J. Shanghai Jiao Tong Univ. 2022, 56, 1447–1452. [Google Scholar]
  18. Wang, J.; Zhang, T.; Tong, J.W.; Yan, Y.X. Application of KF/UPF in Initial Alignment of Large Azimuth Misalignment of SINS. Syst. Eng. Electron. 2018, 40, 2775–2781. [Google Scholar]
  19. Laan, N.; Cohen, M.; Arsenault, J.; Forbes, R. The Invariant Rauch-Tung-Striebel Smoother. IEEE Robot. Automat. Lett. 2020, 5, 5067–5074. [Google Scholar] [CrossRef]
  20. Potokar, E.; Norman, K.; Mangelson, J. Invariant extended kalman filtering for underwater navigation. IEEE Robot. Automat. Lett. 2021, 6, 5792–5799. [Google Scholar] [CrossRef]
  21. Titterton, D.; Weston, J. Strap-down inertial navigation technology-2nd edition—[book review]. IEEE Aerosp. Electron. Syst. Mag. 2005, 20, 33–34. [Google Scholar] [CrossRef]
  22. Xu, Y.; Wang, D.; Shmaliy, Y.S.; Chen, X.Y.; Shen, T.; Bi, S.H. Dual Free-Size LS-SVM Assisted Maximum Correntropy Kalman Filtering for Seamless INS-Based Integrated Drone Localization. IEEE Trans. Ind. Electron. 2023, 71, 9845–9854. [Google Scholar] [CrossRef]
  23. Huang, Y.; Liu, X.; Shao, Q.; Wang, Z. Virtual Metrology Filter-Based Algorithms for Estimating Constant Ocean Current Velocity. Remote Sens. 2023, 15, 4097. [Google Scholar] [CrossRef]
  24. Rossmann, W. Lie Groups: An Introduction through Linear Groups; Oxford University Press: Oxford, UK, 2006. [Google Scholar]
  25. Hall, B.C. Lie Groups, Lie Algebras, and Representations; Springer: New York, NY, USA, 2007. [Google Scholar]
  26. Schwarz, G.W. Smooth Functions Invariant under the Action of a Compact Lie Group. Topology 1975, 14, 63–68. [Google Scholar] [CrossRef]
  27. Lageman, C.; Trumpf, J.; Mahony, R. Gradient-like Observers for Invariant Dynamics on a Lie Group. IEEE Trans. Autom. Control 2008, 55, 367–377. [Google Scholar] [CrossRef]
  28. Cui, B.B.; Wei, X.H.; Chen, X.Y.; Li, J.Y.; Li, L. On Sigma-Point Update of Cubature Kalman Filter for GNSS/INS under GNSS-Challenged Environment. IEEE Trans. Veh. Technol. 2019, 68, 8671–8682. [Google Scholar] [CrossRef]
  29. Wang, G.C.; Xu, X.S.; Yao, Y.Q.; Tong, J.W. A Novel BPNN-based Method to Overcome the GPS Outages for INS/GPS System. IEEE Access 2019, 7, 82134–82143. [Google Scholar] [CrossRef]
  30. Simo, S. Bayesian Filtering and Smoothing; Cambridge University Press: Cambridge, UK, 2013. [Google Scholar]
Figure 1. Schematic diagram of the coordinate system.
Figure 1. Schematic diagram of the coordinate system.
Sensors 24 02945 g001
Figure 2. The dynamic characteristics of the vehicle during the experiment.
Figure 2. The dynamic characteristics of the vehicle during the experiment.
Sensors 24 02945 g002
Figure 3. The trajectory of the vehicle during the experiment.
Figure 3. The trajectory of the vehicle during the experiment.
Sensors 24 02945 g003
Figure 4. Eastward attitude angle error of misalignment angles ( [ 1 1 40 ] T ).
Figure 4. Eastward attitude angle error of misalignment angles ( [ 1 1 40 ] T ).
Sensors 24 02945 g004
Figure 5. Northward attitude angle error of misalignment angles ( [ 1 1 40 ] T ).
Figure 5. Northward attitude angle error of misalignment angles ( [ 1 1 40 ] T ).
Sensors 24 02945 g005
Figure 6. Upward attitude angle error of misalignment angles ( [ 1 1 40 ] T ).
Figure 6. Upward attitude angle error of misalignment angles ( [ 1 1 40 ] T ).
Sensors 24 02945 g006
Figure 7. Latitude alignment error of misalignment angles ( [ 1 1 40 ] T ).
Figure 7. Latitude alignment error of misalignment angles ( [ 1 1 40 ] T ).
Sensors 24 02945 g007
Figure 8. Longitude alignment error of misalignment angles ( [ 1 1 40 ] T ).
Figure 8. Longitude alignment error of misalignment angles ( [ 1 1 40 ] T ).
Sensors 24 02945 g008
Figure 9. Height alignment error of misalignment angles ( [ 1 1 40 ] T ).
Figure 9. Height alignment error of misalignment angles ( [ 1 1 40 ] T ).
Sensors 24 02945 g009
Figure 10. Eastward attitude angle error of misalignment angles ( [ 1 1 95 ] T ).
Figure 10. Eastward attitude angle error of misalignment angles ( [ 1 1 95 ] T ).
Sensors 24 02945 g010
Figure 11. Northward attitude angle error of misalignment angles ( [ 1 1 95 ] T ).
Figure 11. Northward attitude angle error of misalignment angles ( [ 1 1 95 ] T ).
Sensors 24 02945 g011
Figure 12. Upward attitude angle error of misalignment angles ( [ 1 1 95 ] T ).
Figure 12. Upward attitude angle error of misalignment angles ( [ 1 1 95 ] T ).
Sensors 24 02945 g012
Figure 13. Latitude alignment error of misalignment angles ( [ 1 1 95 ] T ).
Figure 13. Latitude alignment error of misalignment angles ( [ 1 1 95 ] T ).
Sensors 24 02945 g013
Figure 14. Longitude alignment error of misalignment angles ( [ 1 1 95 ] T ).
Figure 14. Longitude alignment error of misalignment angles ( [ 1 1 95 ] T ).
Sensors 24 02945 g014
Figure 15. Height alignment error of misalignment angles ( [ 1 1 95 ] T ).
Figure 15. Height alignment error of misalignment angles ( [ 1 1 95 ] T ).
Sensors 24 02945 g015
Figure 16. Eastward attitude angle error of misalignment angles ( [ 20 20 165 ] T ).
Figure 16. Eastward attitude angle error of misalignment angles ( [ 20 20 165 ] T ).
Sensors 24 02945 g016
Figure 17. Northward attitude angle error of misalignment angles ( [ 20 20 165 ] T ).
Figure 17. Northward attitude angle error of misalignment angles ( [ 20 20 165 ] T ).
Sensors 24 02945 g017
Figure 18. Upward attitude angle error of misalignment angles ( [ 20 20 165 ] T ).
Figure 18. Upward attitude angle error of misalignment angles ( [ 20 20 165 ] T ).
Sensors 24 02945 g018
Figure 19. Latitude alignment error of misalignment angles ( [ 20 20 165 ] T ).
Figure 19. Latitude alignment error of misalignment angles ( [ 20 20 165 ] T ).
Sensors 24 02945 g019
Figure 20. Longitude alignment error of misalignment angles ( [ 20 20 165 ] T ).
Figure 20. Longitude alignment error of misalignment angles ( [ 20 20 165 ] T ).
Sensors 24 02945 g020
Figure 21. Height alignment error of misalignment angles ( [ 20 20 165 ] T ).
Figure 21. Height alignment error of misalignment angles ( [ 20 20 165 ] T ).
Sensors 24 02945 g021
Table 1. Mean and variance of attitude error and position error ( [ 1 1 40 ] T ).
Table 1. Mean and variance of attitude error and position error ( [ 1 1 40 ] T ).
AlgorithmAttitude and Position ErrorMeanVariance
SE2(3)/EKFPitch Error−0.7871″0.1214
Roll Error1.2537″0.0926
Heading Error0.0595′0.0087
Latitude Error0.4381 m0.0412
Longitude Error−0.2522 m0.0201
Height Error−0.3418 m0.0603
EKFPitch Error0.9947″0.1216
Roll Error1.9856″0.0945
Heading Error−0.3037′0.0107
Latitude Error0.4882 m0.0415
Longitude Error−0.3065 m0.0203
Height Error−0.5243 m0.0616
Table 2. Mean and variance of attitude error and position error ( [ 1 1 95 ] T ).
Table 2. Mean and variance of attitude error and position error ( [ 1 1 95 ] T ).
AlgorithmAttitude and Position ErrorMeanVariance
SE2(3)/EKFPitch Error11.0881″0.2950
Roll Error−7.7115″0.1303
Heading Error−1.0848′0.0259
Latitude Error0.2507 m0.1201
Longitude Error0.8113 m0.0552
Height Error1.6058 m0.1191
EKFPitch Error57.2824″0.3147
Roll Error−17.7575″0.1366
Heading Error−2.5052′0.0445
Latitude Error0.5569 m0.1238
Longitude Error2.3491 m0.0612
Height Error3.1618 m0.1307
Table 3. Mean and variance of attitude error and position error ( [ 20 20 165 ] T ).
Table 3. Mean and variance of attitude error and position error ( [ 20 20 165 ] T ).
AlgorithmAttitude and Position ErrorMeanVariance
SE2(3)/EKFPitch Error 13.641 0.4387
Roll Error 17.0 1.2976
Heading Error 4.5116 0.2122
Latitude Error−1.8307 m0.0204
Longitude Error2.2376 m0.1855
Height Error4.1076 m0.1372
EKFPitch Error 5.4181 × 10 2 0.4905
Roll Error 2.5211 × 10 2 2.2026
Heading Error 38.9656 × 10 2 0.1851
Latitude Error5.3826 m0.1123
Longitude Error22.3550 m0.6989
Height Error16.1874 m0.1525
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sun, J.; Chen, Y.; Cui, B. An Improved Initial Alignment Method Based on SE2(3)/EKF for SINS/GNSS Integrated Navigation System with Large Misalignment Angles. Sensors 2024, 24, 2945. https://doi.org/10.3390/s24092945

AMA Style

Sun J, Chen Y, Cui B. An Improved Initial Alignment Method Based on SE2(3)/EKF for SINS/GNSS Integrated Navigation System with Large Misalignment Angles. Sensors. 2024; 24(9):2945. https://doi.org/10.3390/s24092945

Chicago/Turabian Style

Sun, Jin, Yuxin Chen, and Bingbo Cui. 2024. "An Improved Initial Alignment Method Based on SE2(3)/EKF for SINS/GNSS Integrated Navigation System with Large Misalignment Angles" Sensors 24, no. 9: 2945. https://doi.org/10.3390/s24092945

APA Style

Sun, J., Chen, Y., & Cui, B. (2024). An Improved Initial Alignment Method Based on SE2(3)/EKF for SINS/GNSS Integrated Navigation System with Large Misalignment Angles. Sensors, 24(9), 2945. https://doi.org/10.3390/s24092945

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