Next Article in Journal
A Survey of Vehicular VLC Methodologies
Previous Article in Journal
Enhancing Smart Building Surveillance Systems in Thin Walls: An Efficient Barrier Design
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Fast Self-Calibration Method for Dual-Axis Rotational Inertial Navigation Systems Based on Invariant Errors

Navigation Research Center, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(2), 597; https://doi.org/10.3390/s24020597
Submission received: 12 December 2023 / Revised: 12 January 2024 / Accepted: 15 January 2024 / Published: 17 January 2024
(This article belongs to the Topic Multi-Sensor Integrated Navigation Systems)

Abstract

:
In order to ensure that dual-axis rotational inertial navigation systems (RINSs) maintain a high level of accuracy over the long term, there is a demand for periodic calibration during their service life. Traditional calibration methods for inertial measurement units (IMUs) involve removing the IMU from the equipment, which is a laborious and time-consuming process. Reinstalling the IMU after calibration may introduce new installation errors. This paper focuses on dual-axis rotational inertial navigation systems and presents a system-level self-calibration method based on invariant errors, enabling high-precision automated calibration without the need for equipment disassembly. First, navigation parameter errors in the inertial frame are expressed as invariant errors. This allows the corresponding error models to estimate initial attitude even more rapidly and accurately in cases of extreme misalignment, eliminating the need for coarse alignment. Next, by utilizing the output of a gimbal mechanism, angular velocity constraint equations are established, and the backtracking navigation is introduced to reuse sensor data, thereby reducing the calibration time. Finally, a rotation scheme for the IMU is designed to ensure that all errors are observable. The observability of the system is analyzed based on a piecewise constant system method and singular value decomposition (SVD) observability analysis. The simulation and experimental results demonstrate that this method can effectively estimate IMU errors and installation errors related to the rotation axis within 12 min, and the estimated error is less than 4%. After using this method to compensate for the calibration error, the velocity and position accuracies of a RINS are significantly improved.

1. Introduction

Inertial navigation systems (INSs) are widely used in various automated unmanned vehicles for their advantages, such as good autonomy and strong concealment [1]. The precision of inertial navigation parameters directly affects the performance of an INS. In addition to using high-precision inertial devices, rotation modulation techniques can be introduced to improve the accuracy of an INS. The accuracy of rotational inertial navigation systems (RINSs) can be improved several times when using sensors of the same level. RINSs can self-calibrate by rotating their own gimbal mechanism to excite error parameters, which is more convenient and effective for error calibration [2,3]. Moreover, as the manufacturing cost and volume of inertial devices decrease, RINSs are gradually being applied to various systems and equipment [4,5,6,7,8].
Inertial measurement units (IMUs) typically require high-precision turntable calibration before leaving the factory. However, the service life of an RINS is generally more than ten years. During this period, IMUs can experience changes in the nominal error parameters due to component aging, the deformation of structural components that support the inertial devices, and other factors [9,10]. If the original calibration parameters are used over time, it will lead to a decrease in the accuracy of the inertial navigation system, which can have a negative impact on operational efficiency. Therefore, it is necessary to periodically conduct the high-precision autonomous calibration of RINSs to ensure the system’s performance when in use. However, the existing calibration methods take too long to calibrate and cannot quickly self-calibrate.
Currently, calibration methods are divided into discrete calibration and system-level calibration [11]. The former method involves removing the inertial navigation system, performing calibration using specialized high-precision turntables, and then reinstalling the inertial navigation system. The process is intricate and may introduce new installation errors [12,13]. The latter method uses navigation errors as observation variables to achieve the optimal estimation of inertial navigation system calibration parameters. This method eliminates the need for specific calibration sites and high-cost equipment. RINSs can be conveniently self-calibrated without an external turntable by exciting error parameters through the rotation of their own gimbal mechanisms. A 30-position calibration scheme was proposed for a dual-axis RINS in reference [14], which could calibrate 21 IMU errors by periodically rotating the axis.
Traditional system-level calibration methods for RINSs typically focus on calibrating IMU errors, including zero bias, scale factor errors, and installation errors. Due to the introduction of the gimbal mechanism, the installation errors between the IMU and the rotation axis also need to be accurately calibrated. Reference [15] utilized the difference in attitude between the two RINSs to estimate installation errors, but it could not achieve the self-calibration of the RINS. Reference [16] proposed a system-level self-calibration method for a dual-axis RINS to address two types of installation errors. However, this method requires the knowledge of an accurate initial attitude.
Currently, the primary method for system-level calibration uses the extended Kalman filter (EKF) to estimate errors [17,18]. In the traditional EKF framework, state error definitions, such as position and velocity, only consider magnitude differences and ignore directional differences. When the initial misalignment error is significant, it can lead to inconsistencies in the defined error coordinates [19]. Furthermore, due to the influence of specific force terms in the system matrix, the EKF defined with a linear error state encounters the problem of inconsistent variance estimation. This leads to the filtering estimate of errors that were initially unobservable, causing significant deviations and overly optimistic covariance estimates. To address these issues, the error equation can be projected and transformed to adapt to different application environments. Reference [20] eliminated the need for the high-speed integration of specific forces in the calculation of the Kalman filter transformation matrix by transforming the velocity error, which improved the accuracy and stability of the system. Reference [21] proposed the state transformation extended Kalman filter (ST-EKF), whereby the specific force vector in the system error model is replaced by the nearly constant gravity vector in a geographical coordinate system, ensuring consistent error state definitions and resolving the variance estimation inconsistency problem. In view of the strong non-linearity and high-dimensional problems of the RINS calibration error equation, the filtering method of sampling the probability space is also widely used in system-level calibration [22,23].
In recent years, the Lie group and manifold theories have been applied to some inertial-based applications [24,25]. By representing the states and their errors in double-direct space-isometric groups  S E 2 ( 3 ) , new group-based system equations and error equations have been designed [26]. If the system is a group affine, even when the attitude deviation is significant, the attitude can be obtained through a linear error model [27]. Based on the existing EKF framework, Barrau used the Lie group theory to construct navigation parameters, establishing a navigation parameter error equation based on invariant errors, and used an invariant extended Kalman filter (IEKF) to unify the initial alignment process, eliminating the need for coarse alignment [28]. Reference [29] studied the characteristics of  S E k ( 3 )  and the specific mathematical expressions of left-invariant and right-invariant errors under various coordinate frames. The right-invariant EKF projects the gravity vector into the inertial frame and has similar properties to the ST-EKF. The left-invariant EKF eliminates the influence of navigation parameters, and the filter performance does not depend on the selected coordinate system. Reference [30] investigated an attitude estimation method using the  S E 2 ( 3 ) EKF  with consideration of constant drift of the gyro. Therefore, this paper represents navigation parameter errors as invariant errors, transforming a non-linear system into a linear system, thereby improving the accuracy of calibration in cases of large misalignment angles.
To reduce the calibration time and improve calibration accuracy, this paper designs a rotation scheme and utilizes the output of a gimbal mechanism to establish extended measurement equations, rendering all errors observable. Subsequently, the observability of the system is analyzed based on the piecewise constant system (PWCS) observability method and singular value decomposition (SVD) observability analysis. Finally, the paper introduces backtracking navigation [31,32], enabling the sensor data to be reused to enhance the estimation accuracy within a short period.
The rest of this paper is organized as follows. Section 2 introduces the mathematical foundations for expressing navigation states using  S E k ( 3 ) . Section 3 presents the IMU error model, gimbal mechanism errors, and the navigation error state model based on  S E 2 ( 3 ) , demonstrating that this system is a group affine. Section 4 describes the designed rotation scheme, backtracking scheme, angular velocity constraints (AVCs) during rotation, the establishment of the measurement equations, and the observability analysis of the system. Section 5 validates the proposed method through simulations and experiments. Section 6 provides a summary and outlines future prospects.

2. Mathematical Preliminaries and the Framework of the Proposed Method

In order to apply the Lie group theory [33] in RINS system-level self-calibration, this section introduces the mathematical preliminaries of matrix Lie group and Lie algebra to facilitate subsequent error modeling and formula derivation. The matrix Lie group  S E 2 ( 3 )  is used to represent the extended pose as [34].
S E 2 ( 3 ) = { T = [ C 0 2 × 3 | v   p I 2 × 2 ] R 5 × 5 | C S O ( 3 ) v , p R 3 × 1 }
where  S O ( 3 )  is the special orthogonal group and  C , v , p  are the attitude rotation matrix, speed, and position, respectively. The inverse of  T  is provided by
T 1 = [ C T 0 2 × 3 | C T v C T p I 2 × 2 ] S E 2 ( 3 )
The Lie algebra associated with  S E 2 ( 3 )  is provided by
s e 2 ( 3 ) = { Ξ = ζ R 5 × 5 | ζ R 9 × 1 }
where
ζ = [ ρ ν φ ] = [ φ × 0 2 × 3 | ν   ρ 0 2 × 2 ] R 5 × 5 , φ , ν , ρ R 3 × 1
The Lie algebra and Lie group can be transformed into each other via exponential mapping  Exp ( )  and logarithmic mapping  Log ( )  as follows:
T = Exp ( ζ ) = [ exp ( φ ) 0 2 × 3 | J l ( φ ) ν   J l ( φ ) ρ I 2 × 2 ]
ζ = Log ( T ) = [ φ a J l 1 ( φ ) v J l 1 ( φ ) p ]
where
φ = cos 1 ( ( tr ( C ) 1 ) / 2 )
exp ( φ ) = cos φ I 3 × 3 + ( 1 cos φ ) a a T + sin φ ( a × )
J l ( φ ) = ( sin φ / φ ) I 3 × 3 + ( 1 ( sin φ / φ ) ) a a T + ( ( 1 cos φ ) / φ ) ( a × )
J l 1 ( φ ) = ( φ / 2 ) ( cot φ / 2 ) I 3 × 3 + ( 1 ( φ / 2 ) ( cot φ / 2 ) ) a a T ( φ / 2 ) ( a × )
φ = φ a  is the rotation vector,  φ  denotes the rotation angle, and  a  is the unit-length axis of rotation.  ( × )  denotes the skew-symmetric matrix operation. When  φ  is assumed to be small, then  exp ( φ ) I 3 × 3 + ( φ × ) J l ( φ ) I 3 × 3 + ( φ × ) / 2 , and  J l 1 ( φ ) I 3 × 3 ( φ × ) / 2 .
In addition to the traditional IMU error, this article considers the installation error of the gimbal mechanism and stores the gyro, accelerometer, and gimbal mechanism signals outputs in real time. In order to quickly and accurately estimate the error parameters, the navigation parameters are established in the form of invariant errors, the linear Kalman filter is used to accurately estimate the calibration errors, and the stored data is reused through the backtracking algorithm to shorten the filter convergence time. Finally, a fast, autonomous, and accurate dual-axis rotation inertial navigation system-level self-calibration algorithm is achieved. The structure of the Kalman filter based on the invariant error during self-calibration is shown in Figure 1.

3. Gimbal Mechanism Installation Error Model, IMU Error Model, and Navigation Error Model

To better analyze the mechanism of the calibration process, several coordinate systems are defined in Table 1.

3.1. Gimbal Mechanism Installation Error Model

For a dual-axis RINS, the IMU is installed on the dual-axis gimbal mechanism. Due to gimbal mechanism installation errors, the r-frame needs to be established to describe the spatial relationship between the r-frame and the s-frame. There is a transformation between the r-frame and the s-frame in the initial state, which is shown in Figure 2.
Generally, the gimbal mechanism installation error is small, and the transformation matrix between the r-frame and the s-frame is represented as:
C r s = [ 1 γ β γ 1 α β α 1 ] = I 3 × 3 + ( μ × )
where  μ = [ α , β , γ ] T  is the installation angle error between the s-frame and the r-frame.

3.2. IMU Error Model

The coordinate system formed with the output of the three-axis accelerometer is the a-frame. The coordinate system formed with the output of the three-axis gyroscope is the g-frame. Due to the installation errors, the a-frame and g-frame are both non-orthogonal coordinate systems. Under normal circumstances, the a-frame and the g-frame are not coincident. Therefore, it is necessary to convert the three-axis acceleration output and the three-axis gyroscope output into the same sensor coordinate system, which is the s-frame. The definition of the s-frame is as follows [16]:  y s  is in the  x a y a  plane,  x s  coincides with  x a , while  z s  forms a right-handed orthogonal frame with  x s  and  y s . The spatial relationships among the g-frame, a-frame, and s-frame are shown in Figure 3.
The IMU error model is defined in the s-frame. The gyros and accelerometers with errors are described as follows:
ω i s s = ( I 3 × 3 C g s K g ) ω i g g ε f i s s = ( I 3 × 3 C a s K a ) f i a a
where  ε = [ ε x , ε y , ε z ] T  and  = [ x , y , z ] T  are the constant gyro drifts and accelerometer biases.  ω i g g = [ ω i g x g , ω i g y g , ω i g z g ] T  and  f i a a = [ f i a x a , f i a y a , f i a z a ] T  are the input of the gyro and accelerometer.  K g  and  K a  are the scale factor errors of the gyros and accelerometers.  k g i ( i = x , y , z )  is the gyro-scale factor errors corresponding to the Xg-axis, Yg-axis, and Zg-axis.  k a i ( i = x , y , z )  is the accelerometer-scale factor errors corresponding to the Xg-axis, Yg-axis, and Zg-axis. They are expressed as follows:
K g = [ k g x 0 0 0 k g y 0 0 0 k g z ] , K a = [ k a x 0 0 0 k a y 0 0 0 k a z ]
C g s C a s  are the installation errors of the gyros and accelerometers. Through Figure 3 C g s  and  C a s  can be expressed as follows:
C g s = [ 0 S g x z S g x y S g y z 0 S g y x S g z y S g z x 0 ] = ( S g )                   C a s = [ 0 0 0 S a y z 0 0 S a z y S a z x 0 ] = ( S a )
where  S g = [ S g x z , S g x y , S g y z , S g y x , S g z y , S g z x ] T  is the installation angle error between the g-frame and s-frame.  S a = [ S a y z , S a z y , S a z x ] T  is the installation angle error between the a-frame and s-frame.  ( )  is the operation that constructs  S g  into  C g s ( )  is the operation that constructs  S a  into  C a s .
Through the above analysis, it can be confirmed that the gyro and accelerometer errors in the r-frame are expressed as:
δ ω i r r = [ ( S g ) + K g + ( μ × ) ] ω i r r + ε + w g δ f i r r = [ ( S a ) + K a + ( μ × ) ] f i r r + + w a
where  ε , , K g , K a , S g , S a , μ  are random constant variables, that is, their derivatives are all zero.  w g  is the measurement noise of the gyros.  w a  is the measurement noise of the accelerometer.

3.3. Navigation Error State Model in the Inertial Frame

According to the SINS mechanization in the in0-frame in [24], the transformed dual-axis RINS mechanization in the in0-frame is provided by:
[ C ˙ r i n 0 v ¯ ˙ e n i n 0 r ˙ e n i n 0 ] = [ C r i n 0 ( ω i r r × ) C r i n 0 f i r r + C n i n 0 G n v ¯ e n i n 0 ]
where  C r i n 0  is the direction of the cosine matrices from the r-frame to the in0-frame.  v ¯ e n i n 0  is the auxiliary velocity vector.  r e n i n 0  is the position in the in0-frame.  ω i r r  is the angular velocity expressed in the r-frame measured with the gyros.  f i r r  is the specific force in the r-frame measured with the accelerometers.  G n  is the global gravitational vector, and the auxiliary velocity vector is provided by:
v ¯ e n i n 0 = v e n i n 0 + ( ω i n 0 e i n 0 × ) r e n i n 0
where  v e n i n 0  is the velocity in the in0-frame.  ω i n 0 e i n 0  is the earth’s rotation in the in0-frame.
The relationship between  r e n i n 0  and geographical  p e n n = [ λ , L , h ] T  is provided by:
r e n i n 0 = C e i n 0 r e n e = C e i n 0 Γ ( p e n n )
where  λ L , and  h  are the longitude, latitude, and height, respectively, obtained from the GPS.
r e n e = Γ ( p e n n ) = [ ( R N + h ) cos L cos λ ( R N + h ) cos L sin λ [ R N ( 1 e 2 ) + h ] sin L ]
The attitude transformation matrix  C e i n 0  is provided by:
C e i n 0 = C e 0 i n 0 C e e 0
where  C e 0 i n 0  is determined by the initial location and  p e n n ( 0 ) C e e 0  is determined by the earth rotation and calibration time.
C e 0 e = [ cos ( ω i e t ) sin ( ω i e t ) 0 sin ( ω i e t ) cos ( ω i e t ) 0 0 0 1 ] C e 0 i n 0 = [ sin λ 0 cos λ 0 0 sin L 0 cos λ 0 sin L 0 sin λ 0 cos L 0 cos L 0 cos λ 0 cos L 0 sin λ 0 sin L 0 ]
The relationship between  v ¯ e n i n 0  and ground velocity  v e n n
v ¯ e n i n 0 = C e i n 0 ( ω i n 0 e e × ) r e n e + C n i n 0 v e n n
It can be used to initialize  v ¯ e n i n 0  for the velocity calculation in (22) and can also be used as the measurement for the dual-axis RINS calibration.  v e n n  can be obtained using a speed sensor, such as an odometer.
Formulating  C r i n 0 v ¯ e n i n 0 , and  r e n i n 0  as elements of the matrix Lie group  S E 2 ( 3 ) :
χ = [ C r i n 0 0 2 × 3 | v ¯ e n i n 0 r e n i n 0 I 2 × 2 ]
According to Equation (23), the differential equation of the  χ  can be calculated as:
χ ˙ = Γ u ( χ ) = χ W 1 + W 2 χ = [ 0 G i n 0 0 0 0 1 0 0 0 ] χ + χ [ ( ω i r r × ) f i r r 0 0 0 1 0 0 0 ]
It is obvious to verify that any solution  χ 1 , χ 2  of Equation (23) satisfies Equation (24). Hence, the dynamical equation  Γ u ( χ )  is group affine. The group-affine system owns the log-linear property of the corresponding error propagation [23]. Such striking property is just the fundamental of the linear KF based an initial alignment with an arbitrary large misalignment [24].
Γ u ( χ 1 ) χ 2 + χ 1 Γ u ( χ 2 ) χ 1 Γ u ( I 5 ) χ 2 = Γ u ( χ 1 χ 2 )
If  χ  represents the ground truth value and  χ ˜  represents the estimated value of  χ , the invariant error  δ χ  is defined in Table 2. In error modeling based on  S O ( 3 ) , only the attitude is established in the Lie group space, and the velocity error and position error are established in the Euclidean space, which is  δ χ S O ( 3 ) + R 6 . Different from the error modeling based on  S O ( 3 ) , the invariant error considers the direction errors of velocity and position and mathematically solves the problem of inconsistent coordinate systems defined by velocity error and position error, which is  δ χ S E 2 ( 3 ) .

3.3.1. Error State Model Based on  S O ( 3 )

Derive both sides of Equation (27) simultaneously and substitute it into Equation (16). It can be deduced that the error model in the inertial frame with  φ i n 0 δ v ¯ e n i n 0 , and  δ r e n i n 0  as navigation parameter errors is:
φ ˙ i n 0 = C ˜ r i n 0 δ ω i r r δ v ¯ ˙ e n i n 0 = ( C ˜ r i n 0 f ˜ i r r ) × φ i n 0 + C ˜ r i n 0 δ f i r r δ r ˙ e n i n 0 = δ v ¯ e n i n 0
Derive both sides of Equation (29) simultaneously and substitute it into Equation (16). It can be deduced that the error model in the r-frame with  φ r δ v ¯ e n i n 0 , and  δ r e n i n 0  as navigation parameter errors is:
φ ˙ r = ( ω ˜ i r r × ) φ r δ ω i r r δ v ¯ ˙ e n i n 0 = C ˜ r i n 0 ( f ˜ i r r × ) φ r + C ˜ r i n 0 δ f i r r δ r ˙ e n i n 0 = δ v ¯ e n i n 0
It can be seen from Equations (30) and (31) that the error state model contains  C ˜ r i n 0 , which is dependent on the trajectory and is affected by the application environment of the calibration model.
For the above two error models, their measurement models are consistent and can be expressed as:
z s o = [ z s o v z s o r ] = [ δ v ¯ e n i n 0 δ r e n i n 0 ] = [ v ¯ ˜ e n i n 0 v ¯ e n i n 0 r ˜ e n i n 0 r e n i n 0 ] + V s o
where  v ¯ ˜ e n i n 0  and  r ˜ e n i n 0  are the velocity and position solved by the dual-axis RINS.  v ¯ ˜ e n i n 0  and  r ˜ e n i n 0  are the velocity and position solved using Equations (18) and (22).
For the error modeling in Equations (31) and (32), the feedback correction for the navigation parameters is provided by Equation (33).
{ C ^ r i n 0 = exp ( φ ^ i n 0 ) C ˜ r i n 0 , C ^ r i n 0 = C ˜ r i n 0 exp ( φ ^ r ) v ¯ ^ e n i n 0 = v ¯ ˜ e n i n 0 δ v ¯ e n i n 0 r ^ e n i n 0 = r ˜ e n i n 0 δ r e n i n 0

3.3.2. Right-Invariant Error State Model Based on  S E 2 ( 3 )

Denote the Lie algebra  d x ¯ R = [ φ i n 0 T , d v ¯ i n 0 T , d r i n 0 T ] T  corresponding to the invariant error  δ χ R s e . According to Equation (5), we have:
C r i n 0 C ˜ i n 0 r = exp ( φ i n 0 ) I 3 × 3 + ( φ i n 0 × )
d v ¯ i n 0 = J l 1 ( φ i n 0 ) ( v ¯ e n i n 0 exp ( φ i n 0 × ) v ¯ ˜ e n i n 0 ) = δ v ¯ e n i n 0 + v ¯ ˜ e n i n 0 × φ i n 0
d r i n 0 = J l 1 ( φ i n 0 ) ( r e n i n 0 exp ( φ i n 0 × ) r ˜ e n i n 0 ) = δ r e n i n 0 + r ˜ e n i n 0 × φ i n 0
The differential equation of  φ i n 0 d v ¯ i n 0 , and  d r ¯ i n 0  are provided by:
( φ ˙ i n 0 × ) = C ˙ r i n 0 C ˜ i n 0 r + C r i n 0 C ˜ ˙ i n 0 r = ( C r i n 0 δ ω i r r ) ×
d v ¯ ˙ i n 0 = δ v ¯ ˙ e n i n 0 + v ¯ ˜ ˙ e n i n 0 × φ i n 0 + v ¯ ˜ e n i n 0 × φ ˙ i n 0 = G i n 0 × φ i n 0 C r i n 0 δ f i r r ( v ¯ ˜ e n i n 0 × ) C r i n 0 δ ω i r r
d r ˙ i n 0 = δ r ˙ e n i n 0 + r ˜ ˙ e n i n 0 × φ i n 0 + r ˜ e n i n 0 × φ ˙ i n 0 = d v ¯ i n 0 ( r ˜ e n i n 0 × ) C r i n 0 δ ω i r r
Add the calibration errors to the state vector. The corresponding state-space model is provided by:
[ φ ˙ i n 0 d v ¯ ˙ i n 0 d r ˙ i n 0 x ˙ c ] = [ 0 3 × 3 0 3 × 3 0 3 × 3 ( G i n 0 × ) 0 3 × 3 0 3 × 3 0 3 × 3 I 3 × 3 0 3 × 3 | F R c _ 0 24 × 33 ] [ φ i n 0 d v ¯ i n 0 d r i n 0 x c ] + [ C r i n 0 0 3 × 3 ( v ¯ ˜ e n i n 0 × ) C r i n 0 C r i n 0 ( r ˜ e n i n 0 × ) C r i n 0 0 3 × 3 _ 0 24 × 6 ] [ w g w a ]
where  x c = [ ε T , T , K g T , K a T , S g T , S a T , μ T ] T  is the state vector of the parameters to be calibrated.  F R c  is provided by:
F R c = [ C r i n 0 0 3 × 3 C r i n 0 D ( ω i r r ) 0 3 × 3 C r i n 0 Γ g ( ω i r r ) 0 3 × 3 C r i n 0 ( ω i r r × ) ( v ¯ ˜ e n i n 0 × ) C r i n 0 C r i n 0 ( v ¯ ˜ e n i n 0 × ) C r i n 0 D ( ω i r r ) C r i n 0 D ( f i r r ) ( v ¯ ˜ e n i n 0 × ) C r i n 0 Γ g ( ω i r r ) C r i n 0 Γ a ( ω i r r ) ( v ¯ ˜ e n i n 0 × ) C r i n 0 ( ω i r r × ) + C r i n 0 ( f i r r × ) ( r ˜ e n i n 0 × ) C r i n 0 0 3 × 3 ( r ˜ e n i n 0 × ) C r i n 0 D ( ω i r r ) 0 3 × 3 ( r ˜ e n i n 0 × ) C r i n 0 Γ g ( ω i r r ) 0 3 × 3 ( r ˜ e n i n 0 × ) C r i n 0 ( ω i r r × ) ]
where  D ( ) Γ g ( )  and  Γ a ( )  are provided by:
D ( a ) = [ a x 0 0 0 a y 0 0 0 a z ] , Γ g ( a ) = [ a y a z 0 0 0 0 0 0 a x a z 0 0 0 0 0 0 a x a y ] , Γ a ( a ) = [ 0 0 0 a x 0 0 0 a x a y ] , a = [ a x a y a z ]
The right-invariant error measurement equation is provided by:
z R s e = [ z R s e v z R s e r ] = [ δ v ¯ e n i n 0 δ r e n i n 0 ] = [ d v ¯ i n 0 + v ¯ ˜ e n i n 0 × φ i n 0 d r i n 0 + r ˜ e n i n 0 × φ i n 0 ] + V R s e
It can be seen from Equation (40) that the navigation parameter system equations in the right-invariant model are only related to the known quantity  G i n 0  and is trajectory-independent. Although the system is group affine and can use the linear-KF to solve non-linear problems, the calibration parameters are strongly coupled with the navigation parameters, which affects the calibration accuracy.
For the right-invariant error modeling in Equation (40), the feedback correction for the navigation parameters is provided by:
C ^ r i n 0 = exp ( φ ^ i n 0 ) C ˜ r i n 0 v ¯ ^ e n i n 0 = exp ( φ ^ i n 0 ) v ¯ ˜ e n i n 0 + J l ( φ ^ i n 0 ) d v ^ i n 0 r ^ e n i n 0 = exp ( φ ^ i n 0 ) r ˜ e n i n 0 + J l ( φ ^ i n 0 ) d r ^ i n 0

3.3.3. Left-Invariant Error State Model Based on  S E 2 ( 3 )

Denote the Lie algebra  d x ¯ L = [ φ r T , d v ¯ r T , d r r T ] T  corresponding to the invariant error  δ χ L s e . According to Equation (5), we have:
C ˜ i n 0 r C r i n 0 = exp ( φ r ) = I 3 × 3 + ( φ r × )
d v ¯ r = J l 1 ( φ r ) ( C ˜ r i n 0 T δ v ¯ e n i n 0 ) = C ˜ r i n 0 T δ v ¯ e n i n 0
d r ¯ r = J l 1 ( φ i n 0 ) ( C ˜ r i n 0 T δ r e n i n 0 ) = C ˜ r i n 0 T δ r e n i n 0
The differential equations of  φ r d v ¯ r , and  d r ¯ r  are provided by:
( φ ˙ r × ) = C ˜ ˙ i n 0 r C r i n 0 + C ˜ i n 0 r C ˙ r i n 0 = [ ( ω i r r × ) φ r ] × ( δ ω i r r × )
d v ¯ ˙ r = C ˜ ˙ r i n 0 T δ v ¯ e n i n 0 C ˜ r i n 0 T δ v ¯ ˙ e n i n 0 = ( ω ˜ i n 0 r r × ) d v ¯ r ( f i r r × ) φ r δ f i r r
d r ¯ ˙ r = C ˜ ˙ r i n 0 T δ r e n i n 0 C ˜ r i n 0 T δ r ˙ e n i n 0 = ( ω ˜ i n 0 r r × ) d r ¯ r + d v ¯ r
Add the calibration errors to the state vector. The corresponding state-space model is provided by:
[ φ ˙ r d v ¯ ˙ r d r ˙ r x ˙ c ] = [ ( ω ˜ i r r × ) 0 3 × 3 0 3 × 3 ( f ˜ i r r × ) ( ω ˜ i r r × ) 0 3 × 3 0 3 × 3 I 3 × 3 ( ω ˜ i r r × ) | F L c 0 3 × 24 _ 0 24 × 33 ] [ φ r d v ¯ r d r r x c ] + [ I 3 × 3 0 3 × 3 0 3 × 3 I 3 × 3 0 3 × 3 0 3 × 3 _ 0 24 × 6 ] [ w g w a ]
where  F L c  is provided by:
F L c = [ I 3 × 3 0 3 × 3 D ( ω i r r ) 0 3 × 3 Γ g ( ω i r r ) 0 3 × 3 ( ω i r r × ) 0 3 × 3 I 3 × 3 0 3 × 3 D ( f i r r ) 0 3 × 6 Γ a ( ω i r r ) ( f i r r × ) ]
The left-invariant error measurement equation is provided by:
z L s e = [ z L s e v z L s e r ] = [ δ v ¯ e n i n 0 δ r e n i n 0 ] = [ C ˜ r i n 0 d v r C ˜ r i n 0 d r r ] + V L s e
Compared with the right-invariant model,  C ˜ r i n 0  is only in the measurement equation of the left-invariant model, and the system equation is only related to  ω i r r  and  f i r r , which is less dependent on the trajectory and its error propagation would be autonomous.
For the left-invariant error modeling in Equation (51), the feedback correction for the navigation parameters is provided by:
C ^ r i n 0 = C ˜ r i n 0 exp ( φ ^ r ) v ¯ ^ e n i n 0 = v ¯ ˜ e n i n 0 + C ˜ r i n 0 J l ( φ ^ r ) d v ^ r r ^ e n i n 0 = r ˜ e n i n 0 + C ˜ r i n 0 J l ( φ ^ r ) d r ^ r
In summary, the left-invariant error model is selected as the filtering algorithm in this paper. The state variables of the proposed method are chosen as follows:
X = [ φ r T , d v ¯ r T , ε T , T , K g T , K a T , S g T , S a T , μ T ] T
where  φ r = [ φ x r , φ y r , φ z r ]  is left-invariant attitude error in the Euler angles.  d v ¯ r = [ d v x r , d v y r , d v z r ]  is left-invariant velocity error.  ε = [ ε x , ε y , ε z ]  is the gyro drift error.  = [ x , y , z ]  is the accelerometer bias error.  K g = [ K g x , K g y , K g z ]  is the gyro-scale factor error.  K a = [ K a x , K a y , K a z ]  is the accelerometer-scale factor error.  S a = [ S a y x , S a z x , S a z y ]  is the accelerometer installation angle error.  S g = [ S g x y , S g x z , S g y x , S g y z , S g z x , S g z y ]  is the gyro installation angle error.  μ = [ μ x , μ y , μ z ]  is installation angle error between the s-frame and r-frame.
The state equation and the measurement equation can be designed as follows:
X ˙ = F X + G W Z L s e = δ v ¯ e n i n 0 = H L s e X + V L s e
where  W  and  V L s e  are the process noise and measurement noise, respectively.  G  is the process noise matrix. The matrix  F  and the matrix  H L s e  are expressed as:
F = [ ( ω ˜ i r r × ) 0 3 × 3 ( f ˜ i r r × ) ( ω ˜ i r r × ) | F L c _ 0 24 × 30 ] , H L s e = [ 0 3 × 3 C ˜ r i n 0 0 3 × 24 ]
The algorithm flow for discretizing the continuous Kalman filter is shown in Algorithm 1.
Algorithm 1 Discrete Kalman filter algorithm flow
(1) State and variance matrix initialization
X ^ 0 = X ( t 0 ) , P ^ 0 = P ( t 0 )
(2) State propagation process
X ^ k / k 1 = Φ k / k 1 X ^ k 1 P ^ k / k 1 = Φ k / k 1 P ^ k 1 Φ k / k 1 T + Γ k 1 Q Γ k 1 T Φ k / k 1 I + F ( t k 1 ) Δ T + F ( t k 1 ) F ( t k 1 ) Δ T 2 / 2 Γ k 1 ( I + F ( t k 1 ) Δ T / 2 + F ( t k 1 ) F ( t k 1 ) Δ T 2 / 6 ) G
(3) Update filter state
K k = P k / k 1 H k T ( H k P k / k 1 H k T + R k ) 1 P k = ( I K k H k ) P k / k 1 X ^ k = X ^ k / k 1 + K k ( z H k X ^ k / k 1 )
(4) Update navigation state
where  Δ T  is the discrete time.  Q  and  R  are device random noise.

4. The Fast Self-Calibration Method and Observability Analysis of the Dual-Axis RINS

4.1. Rotation Scheme Design and Observability Analysis

In order to accurately estimate all the state variables in the self-calibration procedure, the rotation scheme should be properly designed to ensure that all the state variables are observable. Reference [35] proposed design principles for a rotation scheme, which included pointing upward and downward alternately for each accelerometer and a bidirectional rotation on each gyro. According to the above principles and the limitations of the gimbal mechanism, the rotation scheme proposed in this article is shown in Figure 4. In a rotation period, first, the outer axis rotates 270 degrees and −270 degrees (corresponding to A and B in Figure 4). Then, the inner axis rotates to −90 degrees (corresponding to C in Figure 4), and the outer axis rotates bidirectionally (corresponding to D and E in Figure 4). Finally, the inner axis rotates bidirectionally (corresponding to F and G in Figure 4) and returns to position A. During the rotation, the angular velocity is  ω , and it stops for  T s  seconds after each rotation of 90 degrees. According to Equation (16), except for  ε  and  , the remaining calibration errors require  ω i r r  and  f i r r  to excite. In order to excite the error of the gyros, the gyros should be continuously sensitive to the angular velocity of the gimbal mechanism. In the same way, to excite the error of the accelerometers, the accelerometer should be in the vertical direction of the horizontal plane for a long time. Therefore, in order to fully excite the IMU error,  ω  is set to 18°/s and the stop time  T s  is set to 30 s. The complete rotation period is 600 s.
This article uses the PWCS method to verify the feasibility of the designed rotation scheme. In each short time period, the state matrix and measurement matrix are assumed to be constant matrices. At this time, the linear time-varying system transforms into a linear time-invariant system, and the observable matrix of each time period can be obtained. In the i-th time period, the observability matrix of the dynamic system can be expressed as:
Q i = [ H i T , ( H i F i ) T , ( H i F i 2 ) T , , ( H i F i n 1 ) T ] T
where  H i  and  F i  are the measurement matrix and the state matrix in the i-th time interval and  n  is the dimension of the state. The stripped observability matrix (SOM) is:
Q S O M ( r ) = [ Q 1 T , Q 2 T , Q 3 T , , Q r T ] T
When the rank of  Q S O M ( r )  is n, the system is fully observable, while if the rank of  Q S O M ( r )  is less than n, the system is not fully observable. Performing a singular value decomposition on  Q S O M ( r ) , Equation (60) can be obtained:
Q S O M ( r ) = U Σ V T = i = 1 n u i σ i v i T
where  U = [ u 1 , u 2 u r m n ]  and  V = [ v 1 , v 2 v n ]  are  r m n × r m n -dimensional and  n × n -dimensional unitary matrices, respectively.  Σ = [ Σ 0 , 0 n × ( r m n n ) ] T Σ 0 = diag ( σ 1 , σ 2 , σ n ) , and  σ i ( i = 1 , , m )  are the singular values of  Q S O M ( r ) .
Define the singular value  σ k  corresponding to any element  X k  of the state vector  X  to be the singular value  σ i  that maximizes  { ( u i T Y / σ i ) v i } k ( i = 1 , 2 , , m ) . The observability degree of  X k  is defined as:
η k = σ k / σ 0
where  Y = [ Y 1 T , Y 2 T , Y r T ] T Y i = [ Z i T , Z ˙ i T , Z ( n 1 ) T ] T  is the measurement in each time period, and its derivatives of each order.  σ 0  is the singular value corresponding to the state component that can be directly obtained from external information.
Corresponding to the left-invariant error model of Equation (57), the rank of the system is 27.  Q S O M ( r )  is less than the state dimension, indicating that the system is not completely observable. Calculate the observability degrees of each variable at the end of calibration. The observable degrees of  ε x ε y ε z x y z k a x k a y k a z k g x k g y k g z S g x y S g x z S g y x S g y z S g z x S g z y S a y x S a z x S a z y α β , and  γ  are 2.23 × 10−4, 2.25 × 10−4, 1.48 × 10−4, 0.062, 0.063, 0.062, 8.56 × 10−4, 7.51 × 10−5, 2.3 × 10−8, 0.29, 0.33, 0.61, 5.52 × 10−4, 7.71 × 10−6, 6.31 × 10−4, 7.71 × 10−4, 5.78 × 10−4, 5.07 × 10−4, 0.21, 0.26, 0.29, 2.29 × 10−14, 2.29 × 10−14, and 3.64 × 10−16 respectively. The observability degrees of  α β , and  γ  are the worst and is unobservable, which is the reason why the rank of  Q S O M ( r )  and  n  differ by 3. The system error of the gyros and the installation error angle between the r-frame and the s-frame are poorly observable degrees. In response to this situation, the angular velocity constraint equation is added to the measurement equation to render all the parameters observable.

4.2. Angular Velocity Constraint Equation

During the calibration process, the gimbal mechanism drives the IMU rotation. The angular rate of the gimbal mechanism can be utilized as a reference to calibrate the systematic errors of the gyros and the installation error between the s-frame and the r-frame. The outputs of the gyros in the s-frame can be expressed as follows:
ω i g g = C s g C r s ( C i n 0 r ω i e i n 0 + ω b r r + ε ) + w g
where  ω i e i n 0 = [ 0 , ω i e cos L 0 , ω i e sin L 0 ] T  is the angular velocity of the earth expressed in the in0-frame.  L 0  is the latitude in the initial location.  ω b r r  is the angular velocity of the gimbal mechanism.
By substituting Equations (11) and (12) into Equation (62), the second-order small quantity can be ignored and the IMU measurement constraint equation can be expressed as:
ω i g g = ( I 3 × 3 + S g + K g + μ × ) ω 1 + ε + w g ω 1 = C i n 0 r ω i e i n 0 + ω b r r
Equation (63) can be converted into the measurement equation of the calibration parameters:
Z A V C = H A V C X + V A V C = ω i g g ω 1 H A V C = [ 0 3 × 6 , I 3 × 3 , 0 3 × 3 , D ( ω 1 ) , 0 3 × 3 , Γ g ( ω 1 ) , 0 3 × 3 , ( ω 1 ) × ]
Combine Equations (32), (57) and (64) to obtain the new extended measurement equation:
Z L s e A V C = H L s e A V C X + V L s e A V C = [ Z L s e T , Z A V C T ] T , H = [ H L s e T , H A V C T ] T
Z s o A V C = H s o A V C X + V s o A V C = [ z z o v T , Z A V C T ] T , H s o A V C = [ H s o v T , H A V C T ] T
where  H s o v = [ 0 3 × 3 , I 3 × 3 , 0 3 × 24 ] .
For the left-invariant error model with  F  and  H  as the system matrix and the measurement matrix, its rank of  Q S O M ( r )  is 30. The observability degrees of  ε x ε y ε z x y z k a x k a y k a z k g x k g y k g z S g x y S g x z S g y x S g y z S g z x S g z y S a y x S a z x S a z y α β , and  γ  are 0.98, 0.98, 0.98, 2.71 × 10−4, 2.72 × 10−4, 2.71 × 10−4, 0.2, 0.2, 0.2, 2.41 × 10−3, 2.21 × 10−3, 2.36× 10−3, 0.23, 0.21, 0.21, 0.21, 0.22, 0.21, 1.3 × 10−3, 1.1 × 10−3, 8.9 × 10−4, 0.38, 0.35, and 0.23, respectively. The angular velocity constraint equation is added to the measurement equations, which improves the system’s observability, and all calibration parameters can be observed.
From the perspective of the analytical analysis of the variable observability, the accelerometer calibration error has a linear relationship with the first derivative of the velocity error. The gyro error is linearly related to the derivative of the second-order velocity error. It can be seen from the angular velocity constraint equation that there is a linear relationship between the gyros systematic error parameter and  Z g . Adding it to the measurement equation can more directly observe the gyro system error and the installation error between the s-frame and r-frame, effectively shortening the calibration time.

4.3. Backtracking Scheme Design in the Inertial Frame

Backtracking navigation reduces the calibration time by utilizing stored data multiple times. The backtracking scheme is shown in Figure 5. During the calibration process, the SINS update algorithm in the inertial frame is applied for forward navigation. During the first forward navigation process, the data of the IMU and gimbal mechanism are stored. Subsequently, using the stored data multiple times, the calibration parameter error can be accurately estimated in a short time.

5. Simulation Test and Experimental Analysis

This section is devoted to numerically evaluating the calibration performance based on different error state models. The calibration method based on the error state model (30) and measurement model (32) is denoted as RSO-KF. The calibration method based on the error state model (31) and measurement model (32) is denoted as LSO-KF. The calibration method based on the error state model (30) and measurement model (65) is denoted as RSOAVC-KF. The calibration method based on the error state model (31) and measurement model (65) is denoted as LSOAVC-KF. The calibration method based on the error state model (40) and measurement model (66) is denoted as RSEAVC-KF. The calibration method based on the error state model (51) and measurement model (66) is denoted as LSEAVC-KF (the proposed method). The above methods all apply the backtracking method to unify the algorithm process and verify the impact of different error modeling methods on the calibration algorithm. During the simulation and experiment, the dual-axis RINS remained stationary. The location was obtained from a GPS, and the speed was zero.

5.1. Simulation Test

To verify the effectiveness of the proposed algorithm, the simulations are demonstrated in this section. The initial location was 31.94° N, 118.79° E, 15 m. The initial attitude error was 0.1°(roll), 0.1°(pitch), 0.5°(heading). The random gyro drifts were 0.005 deg/h, and the random accelerometer bias was 50 μg. The error of the IMU-scale factor was 100 ppm. The error of the IMU installation angle was 20″. The error of installation angle between the s-frame and r-frame was 200″. The error of the angle of the gimbal mechanism was 0.5″. The error of the angular rate of the gimbal mechanism was 0.7″/s. The rotation scheme is consistent with that described in Section 4.1. The number of backtracking navigations was four.
The simulation was divided into two parts. First, the validity of the angular velocity constraints was verified. Then, the advantages of accuracy were verified and compared with the error modeling method based on  S O ( 3 ) .

5.1.1. Verification of the Observability Analysis

The observability was assessed with a covariance analysis. The convergence process of the root mean square (RMS) of the calibration error covariance is shown in Figure 6. It can be seen from Figure 6 that the RMSs of the  k g z S g y z S a y x α , and  γ  error covariance are large. The observability is consistent with the analysis in Section 4. After introducing AVC constraints into the measurement equation, all the calibration error parameters could be estimated well. Moreover, the RMS of the error covariance became smaller faster.

5.1.2. Comparison of the Calibration Accuracy

In order to eliminate the influence of unobservable installation errors in traditional calibration algorithms, AVC constraints are added as a comparison algorithm to verify the superiority of the calibration algorithm based on invariant error modeling. The four algorithms are evaluated for comparison.
Figure 7, Figure 8, Figure 9 and Figure 10 describe the estimation results of the calibration errors; the red dotted lines represent real values, and the solid lines represent the estimation result of the proposed method (LSEAVC-KF). It took 600 s in the first stage, and then the time spent on backtracking navigation was calculated using the Matlab platform to be 106.1 s. We can see all errors present good convergence with the proposed method within 12 min. Except for the gyro drift error, the rest of the calibration errors converged to near the true value in the first stage. The estimation error of the gyro drift was less than 4%, and the rest of the calibration errors were less than 2%.
In order to further verify the superiority of the proposed method, the estimation accuracy of the remaining three methods was compared, and 100 Monte Carlo simulations were performed for the calibration. The statistical results of the mean error (ME) and root mean square error (RMSE) of each error parameter are shown in Table 3. Among the four methods, LSEAVC-KF had the best effect. The error of the invariant error modeling method was significantly smaller than the error of the error modeling method based on  S O ( 3 ) , which proves that the invariant error can express the error characteristics more accurately.
In order to prove the calibration performance of LSEAVC-KF in the case of large misalignment angles, the calibration errors under different heading angle errors are shown in Table 4. With an increasing heading angle error, it can be observed that the estimation error experiences a slight growth, but it remains within an acceptable range. The proposed algorithm demonstrates greater applicability and can accurately calibrate errors even in harsh environmental conditions. In comparison to traditional calibration methods, it is less reliant on complex instrumentation and offers a faster convergence rate.

5.2. Experiment

In order to verify the feasibility of the proposed method in the actual environment, this paper uses a dual-axis SINS to conduct the experiments, as shown in Figure 11. The dual-axis SINS consists of three fiber optic gyroscopes and three quartz accelerometers. The specific performance is shown in Table 5.
Since the true error of the IMU cannot be known, the accuracy of the error parameter estimation was verified through navigation performance. The higher the accuracy of the IMU parameter, the better the navigation performance. Hence, a navigation experiment with different calibration results was adopted. LSOAVC-KF was chosen as the traditional algorithm for comparison. The velocity errors and position errors under different calibration results are shown in Figure 12 and Figure 13.
The statistical results of the navigation errors are listed in Table 6. In Table 6, Ve and Vn denote the east and north velocities, respectively, in the navigation frame. Pe and Pn denote the east and north positions, respectively, in the navigation frame. It can be seen in Table 6 and in Figure 12 and Figure 13 that the navigation performance of the inertial navigation system was improved with LSEAVC-KF. The maximum error and root mean square error of LSEAVC-KF are smaller than those of LSOAVC-KF.
Therefore, it can be seen from the simulation and experimental results that the proposed algorithm can quickly and accurately estimate the IMU error and installation error between the r-frame and s-frame within 12 min, effectively improving the navigation performance of a dual-axis rotational inertial navigation system.

6. Conclusions

Inertial navigation systems require rapid and accurate error calibration to ensure their performance. A system-level self-calibration method for a dual-axis rotational inertial navigation system based on invariant errors is proposed. Navigation errors are represented as an element of  S E 2 ( 3 ) , allowing the system to converge quickly even during large misalignment angles. Compared to error modeling methods based on  S O ( 3 ) , the proposed method is not affected by the calibration environment and does not require a coarse alignment. The method utilizes the information from the gimbal mechanism to establish angular velocity constraint equations, rendering all errors simultaneously estimated. The rotation scheme for the dual-axis rotational inertial navigation system is designed and discussed, and the observability of the system is analyzed using the PWCS method and an SVD-based observability analysis method. A backtracking navigation algorithm is used to shorten the calibration time. Finally, the superiority of the proposed algorithm is verified through simulations and laboratory experiments, demonstrating that the algorithm can rapidly and accurately estimate IMU systematic errors and installation errors between the r-frame and s-frame. Therefore, the proposed calibration algorithm has significant practical value and can significantly improve navigation accuracy.

Author Contributions

Conceptualization, X.S. and P.L.; methodology, X.S.; software, X.S.; validation, X.S. and W.G.; formal analysis, X.S. and R.L.; investigation, X.S.; resources, X.S. and J.L.; data curation, X.S.; writing—original draft preparation, X.S.; writing—review and editing, X.S.; visualization, X.S.; supervision, J.L.; project administration, J.L.; funding acquisition, J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is sponsored by the National Key Research and Development Program (2022YFB3904300), Key Program of the National Natural Science Foundation of China (U2233215), National Natural Science Youth Foundation of China (62203216), and Natural Science Youth Foundation of the Jiangsu Province (BK20220886).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Groves, P.D. Principles of GNSS, inertial, and multisensor integrated navigation systems, [Book review]. IEEE Aerosp. Electron. Syst. Mag. 2015, 30, 26–27. [Google Scholar] [CrossRef]
  2. Wang, B.; Ren, Q.; Deng, Z.; Fu, M. A self-calibration method for nonorthogonal angles between gimbals of rotational inertial navigation system. IEEE Trans. Ind. Electron. 2014, 62, 2353–2362. [Google Scholar] [CrossRef]
  3. Ren, Q.; Wang, B.; Deng, Z.; Fu, M. A multi-position self-calibration method for dual-axis rotational inertial navigation system. Sens. Actuators A Phys. 2014, 219, 24–31. [Google Scholar] [CrossRef]
  4. Zhang, Q.; Wang, L.; Liu, Z.; Zhang, Y. Innovative self-calibration method for accelerometer scale factor of the missile-borne RINS with fiber optic gyro. Opt. Express 2016, 24, 21228–21243. [Google Scholar] [CrossRef] [PubMed]
  5. Wei, Q.; Zha, F.; Chang, L. Novel rotation scheme for dual-axis rotational inertial navigation system based on body diagonal rotation of inertial measurement unit. Meas. Sci. Technol. 2022, 33, 095105. [Google Scholar] [CrossRef]
  6. Li, Q.; Li, K.; Liang, W. A Dual-Axis Rotation Scheme for Long-Endurance Inertial Navigation Systems. IEEE Trans. Instrum. Meas. 2022, 71, 1–10. [Google Scholar] [CrossRef]
  7. Ishibashi, S.; Tsukioka, S.; Yoshida, H.; Hyakudome, T.; Sawa, T.; Tahara, J.; Aoki, T.; Ishikawa, A. Accuracy Improvement of an Inertial Navigation System brought about by the Rotational Motion. In Proceedings of the OCEANS 2007—Europe, Aberdeen, Scotland, 18–21 June 2007; pp. 1–5. [Google Scholar]
  8. Gao, J.; Li, K.; Chen, Y. Study on Integration of FOG Single-Axis Rotational INS and Odometer for Land Vehicle. IEEE Sens. J. 2018, 18, 752–763. [Google Scholar] [CrossRef]
  9. Zhang, Q.; Niu, X.; Shi, C. Impact Assessment of Various IMU Error Sources on the Relative Accuracy of the GNSS/INS Systems. IEEE Sens. J. 2020, 20, 5026–5038. [Google Scholar] [CrossRef]
  10. Zhang, H.; Zheng, W.; Tang, G. Stellar/inertial integrated guidance for responsive launch vehicles. Aerosp. Sci. Technol. 2012, 18, 35–41. [Google Scholar] [CrossRef]
  11. Bo, X.; Feng, S. A FOG Online Calibration Research Based on High-precision three-axis Turntable. In Proceedings of the 2009 International Asia Conference on Informatics in Control, Automation and Robotics, Bangkok, Thailand, 1–2 February 2009; pp. 454–458. [Google Scholar]
  12. Zheng, Z.; Han, S.; Zheng, K. An eight-position self-calibration method for a dual-axis rotational Inertial Navigation System. Sens. Actuators A-Phys. 2015, 232, 39–48. [Google Scholar] [CrossRef]
  13. Sun, F.; Xia, J.; Ben, Y.; Zhang, X.; IEEE. Four-position Drift Measurement of SINS Based on Single-axis Rotation. In Proceedings of the 2012 IEEE/ION Position, Location and Navigation Symposium (PLANS), Myrtle Beach, SC, USA, 23–26 April 2012; pp. 818–823. [Google Scholar]
  14. Wei, Q.; Zha, F.; He, H.; Li, B. An Improved System-Level Calibration Scheme for Rotational Inertial Navigation Systems. Sensors 2022, 22, 7610. [Google Scholar] [CrossRef]
  15. Hu, P.; Xu, P.; Chen, B.; Wu, Q. A Self-Calibration Method for the Installation Errors of Rotation Axes Based on the Asynchronous Rotation of Rotational Inertial Navigation Systems. IEEE Trans. Ind. Electron. 2018, 65, 3550–3558. [Google Scholar] [CrossRef]
  16. Bai, S.; Lai, J.; Lyu, P.; Xu, X.; Liu, M.; Huang, K. A System-Level Self-Calibration Method for Installation Errors in A Dual-Axis Rotational Inertial Navigation System. Sensors 2019, 19, 4005. [Google Scholar] [CrossRef]
  17. Huang, Z.; Du, P.; Kosterev, D.; Yang, B. Application of extended Kalman filter techniques for dynamic model parameter calibration. In Proceedings of the 2009 IEEE Power & Energy Society General Meeting, Calgary, AB, Canada, 26–30 July 2009; pp. 1–8. [Google Scholar] [CrossRef]
  18. Chen, D.; Yuan, P.; Wang, T.; Cai, Y.; Tang, H. A normal sensor calibration method based on an extended Kalman filter for robotic drilling. Sensors 2018, 18, 3485. [Google Scholar] [CrossRef]
  19. Li, K.; Chang, L.; Chen, Y. Common Frame Based Unscented Quaternion Estimator for Inertial-Integrated Navigation. IEEE ASME Trans. Mechatron. 2018, 23, 2413–2423. [Google Scholar] [CrossRef]
  20. Scherzinger, B.M.; Reid, D.B. Modified strapdown inertial navigator error models. In Proceedings of the 1994 IEEE Position, Location and Navigation Symposium (PLANS’94), Las Vegas, NV, USA, 11–15 April 1994; pp. 426–430. [Google Scholar] [CrossRef]
  21. Wang, M.; Wu, W.; Zhou, P.; He, X. State transformation extended Kalman filter for GPS/SINS tightly coupled integration. GPS Solut. 2018, 22, 112. [Google Scholar] [CrossRef]
  22. He, H.; Tang, H.; Xu, J.; Liang, Y.; Li, F. A SINS/USBL System-Level Installation Parameter Calibration With Improved RDPSO. IEEE Sens. J. 2023, 23, 17214–17223. [Google Scholar] [CrossRef]
  23. Cheng, J.; Liu, P.; Wei, Z.; Luo, G. A novel MEMS-RIMU self-calibration method based on gravity vector observation. Meas. Sci. Technol. 2021, 32, 055108. [Google Scholar] [CrossRef]
  24. Zhang, T.; Wu, K.; Song, J.; Huang, S.; Dissanayake, G. Convergence and Consistency Analysis for a 3-D Invariant-EKF SLAM. IEEE Robot. Autom. Lett. 2017, 2, 733–740. [Google Scholar] [CrossRef]
  25. Potokar, E.; Norman, K.; Mangelson, J. Invariant Extended Kalman Filtering for Underwater Navigation. IEEE Robot. Autom. Lett. 2021, 6, 5792–5799. [Google Scholar] [CrossRef]
  26. Barrau, A.; Bonnabel, S. The Invariant Extended Kalman Filter as a Stable Observer. IEEE Trans. Autom. Control. 2017, 62, 1797–1812. [Google Scholar] [CrossRef]
  27. Chang, L.; Qin, F.; Xu, J. Strapdown Inertial Navigation System Initial Alignment Based on Group of Double Direct Spatial Isometries. IEEE Sens. J. 2022, 22, 803–818. [Google Scholar] [CrossRef]
  28. Barrau, A.; Bonnabel, S. Invariant kalman filtering. Annu. Rev. Control Robot. Auton. Syst. 2018, 1, 237–257. [Google Scholar] [CrossRef]
  29. Luo, Y.; Guo, C.; You, S.; Hu, J.; Liu, J. SE2 (3) based Extended Kalman Filter and Smoothing for Inertial-Integrated Navigation. arXiv 2021, arXiv:2102.12897. [Google Scholar]
  30. Chang, L.B. SE(3) based extended Kalman filter for attitude estimation. J. Chin. Inert. Technol. 2020, 28, 499–504, 550. [Google Scholar]
  31. Wang, J.; Chen, X.; Shao, X. An Adaptive Multiple Backtracking UKF Method Based on Krein Space Theory for Marine Vehicles Alignment Process. IEEE Trans. Veh. Technol. 2023, 72, 3214–3226. [Google Scholar] [CrossRef]
  32. Lin, Y.; Miao, L.; Zhou, Z. A high-accuracy initial alignment method based on backtracking process for strapdown inertial navigation system. Measurement 2022, 201, 111712. [Google Scholar] [CrossRef]
  33. Sola, J.; Deray, J.; Atchuthan, D. A micro Lie theory for state estimation in robotics. arXiv 2018, arXiv:1812.01537. [Google Scholar]
  34. Tang, H.; Xu, J.; Chang, L.; Shi, W.; He, H. Invariant Error-Based Integrated Solution for SINS/DVL in Earth Frame: Extension and Comparison. IEEE Trans. Instrum. Meas. 2023, 72, 1–17. [Google Scholar] [CrossRef]
  35. Gao, P.; Li, K.; Wang, L.; Liu, Z. A Self-Calibration Method for Accelerometer Nonlinearity Errors in Triaxis Rotational Inertial Navigation System. IEEE Trans. Instrum. Meas. 2017, 66, 243–253. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the proposed method in this paper.
Figure 1. Schematic diagram of the proposed method in this paper.
Sensors 24 00597 g001
Figure 2. Definition of the s-frame and r-frame in the initial state.
Figure 2. Definition of the s-frame and r-frame in the initial state.
Sensors 24 00597 g002
Figure 3. Definition of the sensor frame, accelerometer frame, and gyro frame. (a) Sensor frame and accelerometer frame; (b) sensor frame and gyro frame.
Figure 3. Definition of the sensor frame, accelerometer frame, and gyro frame. (a) Sensor frame and accelerometer frame; (b) sensor frame and gyro frame.
Sensors 24 00597 g003
Figure 4. One period of rotation scheme.
Figure 4. One period of rotation scheme.
Sensors 24 00597 g004
Figure 5. Backtracking scheme.
Figure 5. Backtracking scheme.
Sensors 24 00597 g005
Figure 6. RMSs of the calibration error covariance.
Figure 6. RMSs of the calibration error covariance.
Sensors 24 00597 g006aSensors 24 00597 g006b
Figure 7. Estimation of IMU bias errors. (a ε x ε y ε z ; (b x y z .
Figure 7. Estimation of IMU bias errors. (a ε x ε y ε z ; (b x y z .
Sensors 24 00597 g007
Figure 8. Estimation of IMU-scale factor errors. (a k a x k a y k a z ; (b k g x k g y k g z .
Figure 8. Estimation of IMU-scale factor errors. (a k a x k a y k a z ; (b k g x k g y k g z .
Sensors 24 00597 g008
Figure 9. Estimation of gyro installation angle errors. (a S g x y S g x z S g y x  (b S g y z S g z x S g z y .
Figure 9. Estimation of gyro installation angle errors. (a S g x y S g x z S g y x  (b S g y z S g z x S g z y .
Sensors 24 00597 g009
Figure 10. Estimation of accelerometer installation angle errors. (a S a y x S a z x S a z y ; estimation of installation angle errors between the s-frame and r-frame. (b α β , and  γ .
Figure 10. Estimation of accelerometer installation angle errors. (a S a y x S a z x S a z y ; estimation of installation angle errors between the s-frame and r-frame. (b α β , and  γ .
Sensors 24 00597 g010
Figure 11. Internal structure of the dual-axis rotational inertial navigation system (RINS).
Figure 11. Internal structure of the dual-axis rotational inertial navigation system (RINS).
Sensors 24 00597 g011
Figure 12. (a) Velocity errors along east. (b) Velocity errors along north.
Figure 12. (a) Velocity errors along east. (b) Velocity errors along north.
Sensors 24 00597 g012
Figure 13. (a) Position errors along east. (b) Position errors along north.
Figure 13. (a) Position errors along east. (b) Position errors along north.
Sensors 24 00597 g013
Table 1. Definitions of the reference frames.
Table 1. Definitions of the reference frames.
Reference FrameDefinition
a-frameAccelerometer installation frame
g-frameGyro installation frame
s-frameSensor frame
s0-frameSensor frame in the initial state
b-frameBody frame, which is the orthogonal reference frame of a dual-axis RINS. In this paper, the b-frame coincides with the s0-frame
n-frameThe navigation frame, which is the orthogonal reference frame aligned with the east-north-up geodetic axes
in0-frameThe navigation inertial frame, which is the fixed axes of the n-frame in the inertial space at the beginning of the calibration process
e-frameEarth-centered earth-fixed (ECEF) frame
e0-frameECEF inertial frame in the initial state
r-framer-frame, which is the orthogonal reference frame of the gimbal mechanism. Its three axes are defined by  x r , y r , z r x r  coincides with the outer axis,  z r  coincides with the inner axis, while  y r  forms a right-handed orthogonal frame with  x r  and  z r  
r0-framer-frame in the initial state
Table 2. Definitions of invariant error.
Table 2. Definitions of invariant error.
Invariant Error Model Based on  S E 2 ( 3 ) NameCorresponds to Error Model Based on  S O ( 3 )
δ χ R s e = χ χ ˜ 1 = [ C r i n 0 C ˜ r i n 0 T 0 2 × 3 | v ¯ e n i n 0 C r i n 0 C ˜ r i n 0 T v ¯ ˜ e n i n 0 r e n i n 0 C r i n 0 C ˜ r i n 0 T r ˜ e n i n 0 I 2 × 2 ]  (26) Right
Invariant
I 3 × 3 + ( φ i n 0 × ) = C r i n 0 C ˜ r i n 0 T δ v ¯ e n i n 0 = v ¯ ˜ e n i n 0 v ¯ e n i n 0 δ r e n i n 0 = r ˜ e n i n 0 r e n i n 0  (27)
δ χ L s e = χ ˜ 1 χ = [ C ˜ r i n 0 T C r i n 0 0 2 × 3 | C ˜ r i n 0 T ( v ¯ e n i n 0 v ¯ ˜ e n i n 0 ) C ˜ r i n 0 T ( r e n i n 0 r ˜ e n i n 0 ) I 2 × 2 ]  (28) Left
Invariant
I 3 × 3 + ( φ r × ) = C ˜ r i n 0 T C r i n 0 δ v ¯ e n i n 0 = v ¯ ˜ e n i n 0 v ¯ e n i n 0 δ r e n i n 0 = r ˜ e n i n 0 r e n i n 0  (29)
Table 3. Simulation results from the four methods.
Table 3. Simulation results from the four methods.
ParameterRSEAVC-KFRSOAVC-KFLSEAVC-KFLSOAVC-KF
MERMSEMERMSEMERMSEMERMSE
ε x  (°/h)0.00030.00420.00050.00580.00020.00350.00040.0052
ε y  (°/h)0.00060.00410.00030.00420.00050.00390.00020.0042
ε z  (°/h)−0.00030.0034−0.00040.0039−0.00010.0020−0.00050.0039
x  (μg)−0.21310.8136−0.27940.8935−0.14340.7252−0.23270.868
y  (μg)0.15950.740.36890.81360.07960.73750.20490.744
z  (μg)0.0530.39530.12450.44520.04360.37710.09830.44
k a x  (ppm)0.14470.16680.18080.19380.14290.16440.16440.169
k a y  (ppm)−0.14930.1658−0.1520.169−0.14640.1646−0.15170.1674
k a z  (ppm)−0.14490.1584−0.15610.1739−0.13960.1535−0.15530.1585
k g x  (ppm)0.03210.7175−0.20590.7305−0.02220.6592−0.18090.7255
k g y  (ppm)0.42860.7040.52411.9110.30870.67760.50950.9855
k g z  (ppm)0.02670.32160.06540.40940.00550.2980.0610.389
S g x y  (″)0.03760.30520.19860.35490.0150.25790.18280.3166
S g x z  (″)0.05940.26300.08110.28880.05230.22640.06970.2843
S g y x  (″)0.23520.36580.44520.51740.22150.34410.39350.51
S g y z  (″)0.05320.16630.07620.17590.05260.13750.070.1702
S g z x  (″)0.35090.44420.49050.54140.33380.41420.47500.4712
S g z y  (″)0.45190.48010.47410.49910.33290.35470.45230.4833
S a y x  (″)0.14720.35180.42570.53640.12960.34540.38660.4828
S a z x  (″)0.37520.59530.5420.65520.3650.46480.52140.6058
S a z y  (″)0.39490.47420.51940.58360.33630.44680.51330.5067
α  (″)0.15220.29780.27480.31390.13010.17540.25260.3022
β  (″)−0.1150.3132−0.25360.3381−0.1120.2757−0.21360.3158
γ  (″)0.00360.2950.18520.35630.00250.24610.16370.3353
Table 4. Simulation results of different heading angles.
Table 4. Simulation results of different heading angles.
Parameter30°60°90°120°150°180°
ε x  (°/h)0.00010.00020.00030.00040.00020.00030.0004
ε y  (°/h)0.00020.00020.00030.00040.000050.00040.0005
ε z  (°/h)−0.00020.0003−0.00030.0003−0.00030.0002−0.0003
x  (μg)−0.82820.8047−0.76210.678−0.67230.9058−1.1413
y  (μg)−0.44690.4493−0.49530.4189−0.49320.5359−0.7016
z  (μg)0.28360.0030.44080.54250.06340.35090.8046
k a x  (ppm)0.23100.21000.15280.35610.10590.47950.2095
k a y  (ppm)−0.1640−0.1014−0.1108−0.1751−0.21350.1915−0.1813
k a z  (ppm)−0.1199−0.1785−0.1568−0.1210−0.1812−0.1964−0.2170
k g x  (ppm)0.08340.5543−0.79751.10180.9184−0.1006−0.9222
k g y  (ppm)0.56520.46800.25680.47500.5810−0.72521.1168
k g z  (ppm)0.06180.32180.09110.66600.20120.3106−0.6793
S g x y  (″)−0.16690.1989−0.12520.2270−0.13890.15850.5027
S g x z  (″)0.28900.36140.14230.51800.03310.30290.5712
S g y x  (″)−0.3021−0.3992−0.2673−0.5448−0.3001−0.5239−0.7781
S g y z  (″)0.16650.18870.35530.72220.36280.64990.9026
S g z x  (″)0.15060.25190.26450.40420.19840.23530.5648
S g z y  (″)−0.2413−0.1235−0.3786−0.4283−0.4807−0.4395−0.5220
S a y x  (″)0.34480.86150.40110.70380.31260.57320.9560
S a z x  (″)0.11600.39930.33790.38200.53920.32550.9529
S a z y  (″)0.48300.38930.45170.44360.32570.37200.7392
α  (″)0.05630.06110.20520.21740.29380.24040.3120
β  (″)−0.10740.5557−0.18080.3342−0.35890.5153−0.7725
γ  (″)0.03070.05730.09520.04950.38310.63510.6807
Table 5. Specifications of the dual-axis rotational inertial navigation system.
Table 5. Specifications of the dual-axis rotational inertial navigation system.
CharacteristicsDescription
Sampling frequency200 Hz
Gyro-constant drift 0.005 deg/h
Gyro-stochastic drift 0.001 deg/sqrt(h)
Gyro-scale factor repeatability 50 ppm  ( 1 σ )
Gyro-constant drift repeatability 0.001 deg/h  ( 1 σ )
Gyro installation error 0.001 rad
Accelerometer constant bias 10 μg
Accelerometer-scale factor repeatability 50 ppm  ( 1 σ )
Accelerometer constant bias repeatability 10 μg  ( 1 σ )  
Accelerometer installation error 0.001 rad
Rotation angle accuracy0.5″
Rotation velocity accuracy0.7″/s
Table 6. Comparison of the velocity and position errors.
Table 6. Comparison of the velocity and position errors.
Parameter LSEAVC-KFLSOAVC-KF
Ve error (m/s)RMSE
Maximum
0.0595
0.0870
0.1002
0.1515
Vn error (m/s)RMSE
Maximum
0.0403
0.0661
0.0706
0.1092
Pe error (m)RMSE
Maximum
110.0472
175.9277
130.3365
216.0348
Pn error (m)RMSE
Maximum
40.4820
100.6326
82.8674
185.4080
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, X.; Lai, J.; Lyu, P.; Liu, R.; Gao, W. A Fast Self-Calibration Method for Dual-Axis Rotational Inertial Navigation Systems Based on Invariant Errors. Sensors 2024, 24, 597. https://doi.org/10.3390/s24020597

AMA Style

Sun X, Lai J, Lyu P, Liu R, Gao W. A Fast Self-Calibration Method for Dual-Axis Rotational Inertial Navigation Systems Based on Invariant Errors. Sensors. 2024; 24(2):597. https://doi.org/10.3390/s24020597

Chicago/Turabian Style

Sun, Xin, Jizhou Lai, Pin Lyu, Rui Liu, and Wentao Gao. 2024. "A Fast Self-Calibration Method for Dual-Axis Rotational Inertial Navigation Systems Based on Invariant Errors" Sensors 24, no. 2: 597. https://doi.org/10.3390/s24020597

APA Style

Sun, X., Lai, J., Lyu, P., Liu, R., & Gao, W. (2024). A Fast Self-Calibration Method for Dual-Axis Rotational Inertial Navigation Systems Based on Invariant Errors. Sensors, 24(2), 597. https://doi.org/10.3390/s24020597

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