Next Article in Journal
Sea–Sky Line and Its Nearby Ships Detection Based on the Motion Attitude of Visible Light Sensors
Previous Article in Journal
Dual-Channel Spectral Domain Optical Coherence Tomography Based on a Single Spectrometer Using Compressive Sensing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A System-Level Self-Calibration Method for Installation Errors in A Dual-Axis Rotational Inertial Navigation System

1
Navigation Research Center, Nanjing University of Aeronautics and Astronautics, Nanjing 210000, China
2
Overall Design Institute, Hubei Academy of Aerospace Technology, Wuhan 430000, China
3
Shanxi Baocheng Aviation Instrument Corporation, Baoji 721000, China
*
Author to whom correspondence should be addressed.
Submission received: 17 July 2019 / Revised: 5 September 2019 / Accepted: 11 September 2019 / Published: 16 September 2019
(This article belongs to the Section Physical Sensors)

Abstract

:
In a dual-axis rotational inertial navigation system (RINS), there are two kinds of installation errors, nonorthogonal installation errors of inertial sensors, and installation errors between the inertial measurement unit (IMU) and rotation axes. Traditionally, these two errors are not considered simultaneously. Thus, they are calibrated separately by different estimation algorithms and rotation schemes. In this paper, a system-level self-calibration method for installation errors of a dual-axis RINS is proposed. Based on the Kalman filter, the measurement model is reestablished to ensure that all installation errors can be estimated together. First, the relationship between the initial attitude and subsequent attitude of IMU during rotation is used as a constraint to estimate nonorthogonal installation errors of accelerometers, and installation errors between the IMU and rotation axes. Then, the angular rate of the rotation mechanism is used as a reference to estimate nonorthogonal installation errors of the gyros. The rotation scheme of the IMU is designed to make all installation errors observable, and the observability of the system is analyzed based on the piecewise constant system method. Simulation and laboratory experiment results suggest that installation errors can be effectively estimated by the proposed method, thereby avoiding the complex separating process.

1. Introduction

The errors of inertial sensors are the major cause which influences the accuracy of an inertial navigation system (INS). Such errors can be decreased by enhancing the accuracy of inertial sensors. However, the cost and volume of the system can be heavily increased. Instead of using costly inertial sensors, the performance of the INS can be improved through rotary modulation (RM). The constant errors of inertial sensors can be compensated by the rotating inertial measurement unit (IMU) periodically [1,2]. Hence, the rotational INS (RINS) has attracted increasing attention in recent years. For example, the RM technique was applied to enhance the performance of the INS in an autonomous underwater vehicle [3]. A fiber optic gyros (FOG) single-axis RINS was combined with an odometer to improve the positioning and orientation accuracy of a land vehicle [4]. To further improve the performance of the RINS, an improved 8-position rotation scheme and a novel 16-position rotation scheme were proposed for the optical gyro-based RINS [1]. Meanwhile, error propagation characteristics and rotation scheme designing were analyzed [5]. The effect of RINS’s main error source, including constant drift of gyro and zero bias of accelerometer under rotation conditions, was analyzed [6]. The compensation effects of gyros’ stochastic errors were analyzed, and the specific formulae were derived [7].
The accuracy of RINS is influenced by many error parameters [8], among which the installation errors are of key importance. In RINS, installation errors consist of two parts, one is caused by nonorthogonal mounting of inertial sensors. The other is caused by misalignments between IMU axes and rotation axes. Conventionally, researchers have only focused on the calibration of either of these installation errors and designed corresponding calibration algorithms and schemes to estimate one type of error. Following this method makes it complicated to obtain all installation errors because different estimation algorithms and rotation schemes need to be adopted in turn. For the first type of installation errors, one study used velocity measurement of the navigation system to implement calibration [9]. Velocity and position errors were used as a measurement to estimate installation errors and accelerometer nonlinearity errors by applying the Kalman filter [10]. Another study utilized attitude and velocity errors to identify installation errors by applying an extended Kalman filter [11]. The relationship between navigation errors and IMU errors was established using an error model of navigation solution and initial alignment, and then the installation errors were estimated through a total least squares method [12]. The above calibration methods do not consider the effect of installation errors between the IMU and rotation axes, which affects the accuracy of attitude. For the second type of installation error, one study utilized the sensitivity of the gyroscope to the residual rotation, which is caused by installation errors. The errors could be estimated by applying a thin-shell method [13]. Another study used specific force to estimate installation errors through singular value decomposition [14]. The attitude difference of two RINSs was utilized to estimate installation errors via the Kalman filter [15]. Although the above references discuss the calibration of installation errors between the IMU and rotation axes, their methods still cannot obtain all installation errors.
To ensure all error parameters are observable, a specific rotation scheme needs to be designed to implement error excitation in all states. The rotation scheme design principles for errors calibration in RINS was discussed [10]. A self-calibration scheme, which included three separate steps to reduce the correlations between error parameters in tri-axis RINS, was designed [16]. Another study designed a self-calibration scheme for the gimbaled IMU, in which the multi-position and continuous rotation steps were processed alternately, the dynamic errors were suppressed, and the excitation was augmented [17].
In this paper, a system-level self-calibration method for two types of installation errors of a dual-axis RINS is proposed to avoid the complicated separating calibration process. Existing models and algorithms are only applicable for the calibration of either of these two installation errors. Thus, different estimation algorithms must be adopted in turn to obtain all errors. The major breakthrough is that all installation errors in the dual-axis RINS can be calibrated using the proposed system-level method. The point of difference to previous researches is that these two installation errors are considered simultaneously. Based on the Kalman filter, the measurement model is reestablished to ensure that all errors can be estimated together. First, the relationship between initial attitude and sequent attitude calculated from specific forces is used to establish the observation equation, and nonorthogonal installation errors of accelerometers and installation errors between IMU axes and rotation axes can be estimated. Then, the estimated outputs of the first step are employed in the formulation of the relationship between gyroscope and angular rate of rotational mechanism, and nonorthogonal installation errors of gyros can be obtained. The rotation scheme of the IMU is designed to make all installation errors observable, and the observability of the system is analyzed based on the piecewise constant system (PWCS) method. Finally, the proposed self-calibration method is validated through simulation and experiment.

2. Analysis of the IMU Error Model

2.1. Coordinate System and Transformation

A few coordinates need to be established to describe spatial relationships among different parameter variables.

2.1.1. IMU Frame

The outputs of accelerometers and gyroscopes are measured in their respective coordinates, denoted by a-frame and g-frame. The three axes of a-frame are denoted as x a , y a , z a , and the three axes of g-frame are denoted by x g , y g , z g . Because a-frame and g-frame are both nonorthogonal frames, the outputs of inertial sensors need to be transformed to an orthogonal coordinate, called the IMU frame, denoted by s-frame. The three axes of s-frame are denoted by x s , y s , z s . The s-frame is defined as follows [10]: x s coincides with x a ,   y s is in the x a , y a plane, while z s forms a right-handed orthogonal frame with x s and y s . The spatial relationships among a-frame, g-frame, and s-frame are shown in Figure 1.
The transformation matrices from a-frame and g-frame to s-frame are expressed as follows:
C a s = 1 0 0 S a y z 1 0 S a z y S a z x 1 C g s = 1 S g x z S g x y S g y z 1 S g y x S g z y S g z x 1 ,
where S a y z , S a z x , and S a z y are nonorthogonal installation errors of accelerometers and S g x y , S g x z , S g y x , S g y z , S g z x , and S g z y are nonorthogonal installation errors of gyros.
In the process of rotation, s-frame rotates along with inertial sensors. A new frame needs to be introduced to describe the s-frame in the initial state, denoted by s0-frame.

2.1.2. Rotational Mechanism Frame

The rotational mechanism frame needs to be established to describe the spatial relationship between IMU axes and rotation axes. Aimed at the dual-axis RINS, in the initial state, the rotational mechanism frame is denoted by r0-frame; its three axes are defined by x r 0 , y r 0 , z r 0 . Assuming that the inner axis and outer axis are perpendicular, x r 0 coincides with the outer axis, z r 0 coincides with the inner axis, while y r 0 forms a right-handed orthogonal frame with x r 0 and z r 0 . In the initial state, r0-frame and s0-frame do not coincide with each other due to installation errors. There is a transformation between r0-frame and s0-frame, which is shown in Figure 2.
Considering that installation errors are small, the transformation matrix between r0-frame and s0-frame is represented as
C r 0 s 0 = 1 γ α γ 1 β α β 1 ,
where α , β , and γ are installation errors between IMU axes and rotation axes.
Meanwhile, a new frame, denoted by r-frame is used to describe the rotational mechanism frame in the process of rotation. The transformation matrix from r0-frame to r-frame can be represented as follows:
C r 0 r = R x A x R z A z = 1 0 0 0 c A x s A x 0 s A x c A x c A z s A z 0 s A z c A z 0 0 0 1 ,
where A x and A z are the rotation angle of the outer and inner axis, respectively, c and s represent the trigonometric operator cos and sin, respectively. Meanwhile, the transformation matrix between s-frame and r-frame satisfies Equation (2).

2.1.3. Body Frame

The body frame is fixed on the vehicle, denoted by b-frame. The three axes of b-frame are defined as x b , y b , z b . x b axis points rightward, y b axis points forward, and z b points upward. In this paper, b-frame coincides with s0-frame.

2.1.4. Navigation Frame

The navigation frame is denoted by n-frame, its three axes are defined as x n , y n , z n . x n axis points east, y n axis point north, and z n axis points up.

2.2. Error Modeling of Inertial Sensors

Based on spatial relationship from Equation (1), the model of inertial sensors can be expressed as follows:
f s = C a s f a = 1 0 0 S a y z 1 0 S a z y S a z x 1 f x a f y a f z a = f x a S a y z f x a + f y a S a z y f x a + S a z x f y a + f z a ,
where f a is the output of accelerometers measured in a-frame, f s is the output of accelerometers expressed in s-frame.
ω s = C g s ω g = 1 S g x z S g x y S g y z 1 S g y x S g z y S g z x 1 ω x g ω y g ω z g = ω x g S g x z ω y g + S g x y ω z g S g y z ω x g + ω y g S g y x ω z g S g z y ω x g + S g z x ω y g + ω z g ,
where ω g is the output of gyros measured in g-frame, ω s is the output of gyros expressed in s-frame.

3. A Self-Calibration Method Based on the Kalman Filter

The structure of the self-calibration method is shown in Figure 3.
First, the vehicle attitude calculated from the rotation angle and IMU attitude should be the same as the IMU attitude in the initial state. This constraint can be used to establish the first observation equation. Specific forces are introduced to get IMU attitudes at each moment, and angles of the rotational mechanism are used to project IMU attitudes to b-frame. The nonorthogonal installation errors of accelerometers, and installation errors of IMU axes and rotation axes can be estimated through the first step.
Second, the angular rate of rotational mechanism can be utilized as a reference to calibrate nonorthogonal installation errors of gyros. However, the angular rates of the earth also exist in the output of gyros and should be compensated. The estimation results of the first step, along with the angle of the rotational mechanism, can be used to project the angular rate of the earth from n-frame to s-frame. The second observation equation is established based on the relationship among outputs of gyros, angular rate of the earth and rotational mechanism, then the nonorthogonal installation errors of gyros can be estimated.

3.1. Design of the Kalman Filter

The states of the filter are installation errors, which are modeled as random constants, their dynamic models are expressed as follows:
X = [ 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 z S a z x S a z y α β γ ] X k = F X k 1 + W k = I 12 X k 1 + W k ,
where I 12 is a 12-dimension identity matrix, W k is system noise.
Attitude calculated by projecting IMU attitude to b-frame equals to the IMU attitude in the initial state. Thus, this can be regarded as a constraint to construct an observation equation. Based on the attitude calculation process in RINS, the following equation can be obtained
C n b = C s b C n s = C s s 0 C n s = C n s 0 ,
where C n b is the attitude matrix of the vehicle, C n s is the attitude matrix of the IMU during the rotation process, C n s 0 is the attitude matrix of the IMU in the initial state, C s b and C s s 0 are the attitude matrices which can project the IMU attitude to b-frame.
Since there are installation errors between IMU axes and rotation axes, Equation (7) can be rewritten as
C n s = C s 0 s C n s 0 = C r s C r 0 r C s 0 r 0 C n s 0 C n s = C 11 C 12 C 13 C 21 C 22 C 23 C 31 C 32 C 33 C r s = 1 γ α γ 1 β α β 1 C r 0 r = A 11 A 12 A 13 A 21 A 22 A 23 A 31 A 32 A 33 C s 0 r 0 = 1 γ α γ 1 β α β 1 C n s 0 = C 11 0 C 12 0 C 13 0 C 21 0 C 22 0 C 23 0 C 31 0 C 32 0 C 33 0 ,
where C r s and C s 0 r 0 are transfer matrices between the IMU and rotational mechanism, C r 0 r is the transfer matrix calculated from the angle of the rotational mechanism. The third column of C n s and C n s 0 can be represented as follows by specific forces:
C 13 = sin r cos θ = f x s g C 23 = sin θ = f y s g C 33 = cos r cos θ = f z s g C 13 0 = sin r 0 cos θ 0 = f x s 0 g C 23 0 = sin θ 0 = f y s 0 g C 33 0 = cos r 0 cos θ 0 = f z s 0 g
where C 13 , C 23 , C 33 are entries of C n s . C 13 0 , C 23 0 , C 33 0 are entries of C n s 0 . r and θ are roll and pitch of the IMU during the rotation process, r 0 and θ 0 are roll and pitch of the IMU in the initial state, g is the gravity vector, f s and f s 0 are outputs of accelerometers expressed in the rotation process and initial state, respectively.
From Equations (8) and (9), Equation (10) is obtained as follows:
f x s f y s f z s = C r s C r 0 r C s 0 r 0 f x s 0 f y s 0 f z s 0
By substituting f s and f s 0 into Equation (10) with Equation (4), the second-order small quantity can be ignored, then the first observation equation at k can be represented as
Z 1 k = H 1 X k + V 1 k = 0 3 × 6 H ¯ X k + V 1 k Z 1 k = f x a k + A 11 f x a 0 + A 12 f y a 0 + A 13 f z a 0 f y a k + A 21 f x a 0 + A 22 f y a 0 + A 23 f z a 0 f z a k + A 31 f x a 0 + A 32 f y a 0 + A 33 f z a 0 X k = S g x y k S g x z k S g y x k S g y z k S g z x k S g z y k S a y z k S a z x k S a z y k α k β k γ k T H ¯ 11 = A 12 f x a 0 H ¯ 12 = A 13 f y a 0 H ¯ 13 = A 13 f x a 0 H ¯ 14 = A 31 f x a 0 + A 13 f x a 0 + A 32 f y a 0 A 11 f z a 0 + A 33 f z a 0 H ¯ 15 = A 12 f z a 0 A 13 f y a 0 H ¯ 16 = A 21 f x a 0 A 12 f x a 0 + A 11 f y a 0 A 22 f z a 0 A 23 f z a 0 H ¯ 21 = A 22 f x a 0 f x a H ¯ 22 = A 23 f y a 0 H ¯ 23 = A 23 f x a 0 H ¯ 24 = A 23 f x a 0 A 21 f z a 0 H ¯ 25 = A 31 f x a 0 A 32 f y a 0 A 23 f y a 0 + A 22 f z a 0 A 33 f z a 0 H ¯ 26 = A 11 f x a 0 A 22 f x a 0 + A 21 f y a 0 + A 12 f y a 0 + A 13 f z a 0 H ¯ 31 = A 32 f x a 0 H ¯ 32 = A 33 f y a 0 f y a H ¯ 33 = A 33 f x a 0 + f x a H ¯ 34 = A 33 f x a 0 A 11 f x a 0 A 12 f y a 0 A 31 f z a 0 A 13 f z a 0 H ¯ 35 = A 21 f x a 0 + A 22 f y a 0 A 33 f y a 0 + A 32 f z a 0 + A 23 f z a 0 H ¯ 36 = A 32 f x a 0 + A 31 f y a 0
Equation (11) is based on the invariance of the vehicle’s attitude in a static situation. During the calibration process, the IMU rotates with the rotation mechanism while the vehicle is static. Therefore, the initial attitude of the IMU should equal to the vehicle attitude that can be obtained by projecting the subsequent attitude of the IMU during rotation to b-frame. Considering that the vehicle is in a static situation while the IMU rotates with the axes and rotation needs to stop at certain points during the calibration, the effect of a disturbance, such as moving or centripetal acceleration, can be ignored. Specific force is used to calculate the attitude of the IMU. Thus, the measurement z in Equation (11) includes specific force from accelerometers and angular data from the rotation mechanism.
In the rotation process, the outputs of gyros expressed in s-frame can be written as
ω s = ω i e s + ω r s = C n s ω i e n + C r s ω r = C 11 C 12 C 13 C 21 C 22 C 23 C 31 C 32 C 33 0 ω i e c L ω i e s L + 1 γ α γ 1 β α β 1 ω ¯ r o t a ,
where ω s is the output of gyros in s-frame, ω i e s is the angular rate of the earth expressed in s-frame, ω r s and ω ¯ r o t a are angular rates of the rotational mechanism expressed in s-frame and r-frame, respectively. L is the latitude. The estimation results of the first step are utilized in Equation (12). First, the outputs of accelerometers in the initial state are recalculated via Equation (4) with the estimation of the nonorthogonal installation errors. A modified C n s 0 is obtained. Then an estimation of installation errors between IMU axes and rotation axes along with angle of the rotational mechanism is used to calculate C n s through Equation (8). C n s is substituted into Equation (12) and the angular rate of the earth measured in s-frame can be obtained. Estimation of installation errors between IMU axes and rotation axes combined with the angular rate of the rotational mechanism is substituted into Equation (12) to get the residual rotation ω r s . By substituting ω s in Equation (12) with Equation (5), the second observation equation at k can be represented as
Z 2 k = H 2 X k + V 2 k = H ¯ 0 3 × 6 X k + V 2 k Z 2 k = C 12 ω i e n c L + C 13 ω i e n s L + ω ¯ x r o t a k γ ω ¯ y r o t a k + α ω ¯ z r o t a k ω i s x s k C 22 ω i e n c L + C 23 ω i e n s L + γ ω ¯ x r o t a k + γ ω ¯ y r o t a k β ω ¯ z r o t a k ω i s y s k C 32 ω i e n c L + C 33 ω i e n s L α ω ¯ x r o t a k + β ω ¯ y r o t a k + β ω ¯ z r o t a k ω i s z s k X k = S g x y k S g x z k S g y x k S g y z k S g z x k S g z y k S a y z k S a z x k S a z y k α k β k γ k T H ¯ = ω z g ω y g 0 0 0 0 0 0 ω z g ω x g 0 0 0 0 0 0 ω y g ω x g
Equation (13) is based on that angular rate of the rotation mechanism and can be used as a high-precision reference. During the rotation, the output of gyros includes the angular rate of the earth and rotation mechanism and nonorthogonal installation error terms. The angular velocity of each axes can be directly obtained from the rotation mechanism. Thus, the measurement z in Equation (13) includes the angular rate from gyros, earth, and rotation mechanism.

3.2. Rotation Scheme Design and Observability Analysis

An appropriate rotation scheme should be designed to ensure all state variables observable.
The rotation scheme refers to the design principles in study [10], which involve pointing upward and downward alternately for each accelerometer, and bidirectional rotation on each gyro. The specific rotation scheme goes as follows. The outer axis rotates bidirectionally first and stops every time the rotation angle is on a multiple of −90 degree. Second, the inner axis rotates to −90 degree. The outer axis rotates bidirectionally as in the last step. Finally, the inner axis rotates bidirectionally. Then, the first period of rotation is over, and the second repeats the above process. The angular velocity is set as 20 deg/s, the duration for every stopping point is 200 s. Equation (11) is used when the IMU stops during the first period of rotation, which is designed for avoiding the influence of centripetal acceleration. Equation (13) is used when the IMU rotates during the second period of rotation so that the gyro installation error can be activated.
One period of rotation for the installation error calibration is shown in Figure 4, and the angle changes of the inner and outer axis during the whole calibration scheme are shown in Figure 5.

3.2.1. The Installation Errors’ Excitation of the First Period

1. Rotation along x s of the IMU
In the initial state, x s of the IMU nearly coincides with the outer axis of rotational mechanism and the IMU rotates along its xs axis when outer axis rotates. Then, Equation (13), as the rotation angle is on a multiple of −90 degree, can be approximated as
f x a = g γ + g α , f y a = g S a z x A x = π 2 f x a = 2 g γ A x = π f x a = g γ g α , f y a = g S a z x A x = 3 π 2 .
It can be seen that α , γ , and S a z x are activated when the IMU rotates along x s axis.
2. Rotation along y s of IMU
After the inner axis of rotational mechanism rotates to −90 degree, y s of the IMU nearly coincides with the outer axis, and the output of accelerometers at this moment approximates g 0 0 T . When outer axis rotates, the IMU rotates along its y s axis, and Equation (13) can be approximated as
f x a = g S a z y , f y a = g S a y z g β + g γ A x = π 2 f y a = 2 g S a y z + 2 g γ , f z a = 0 A x = π f x a = g S a z y , f y a = g S a y z + g β + g γ A x = 3 π 2
It can be seen that β , γ , S a y z , and S a z y are activated.
3. Rotation along z s of the IMU
After rotation along z s of the IMU, the output of the accelerometers still approximates g 0 0 T . Due to the rotation of the inner axis in the second step, the inner axis cannot rotate −270 degrees. Thus, Equation (13) can be approximated as
f x a = g S a y z , f z a = g S a z x + g S a z y g α g β A z = π 2 f z a = 2 g S a z y 2 g α A z = π
It can be seen that α , β , S a y z , S a z x , and S a z y are activated.

3.2.2. The Installation Errors’ Excitation of the Second Period

When the IMU rotates along the z s axis, the components of the angular rate caused by nonorthogonal installation errors can be expressed as
ω g x = S g x y ω s z ω g y = S g y x ω s z .
It can be seen that S g x y and S g y x are activated. In the same way, S g x z , S g z x , S g y z , and S g z y can be activated when the IMU rotates along the other two axes.
In this paper, the PWCS method is used to validate the observability of the system. When the time interval is small enough, the state matrix and observation matrix can be regarded as constant matrices in every time interval. The linear time-varying system becomes a linear time-invariant system in every time interval, which is the PWCS. In the i t h time interval, the observable matrix of the dynamic system is
Q i = H i T H i A i T H i A 2 i T H i A n 1 i T T ,
where A i and H i are the state matrix and the observation matrix in the i t h time interval, n is the dimension of the state. The total observability matrix (TOM) of the dynamic system is
Q T O r = Q 1 Q 2 e A 1 Δ 1 Q 3 e A 2 Δ 2 e A 1 Δ 1 Q r e A r 1 Δ r 1 e A 1 Δ 1 ,
where i is the i t h time interval, and the stripped observability matrix (SOM) of the dynamic system is
Q S O r = Q 1 T Q 2 T Q 3 T Q r T T .
When the rank of Q S O r is n , the system if fully observable, while the rank of Q S O r is less than n , the system is not fully observable. In our self-calibration scheme, i is set 0.005s, and the rank change of the system is shown in Figure 6.
From Figure 6, we can see that the dynamic system reaches full rank after 0.2 s, which is not the convergence time but indicates the system can be fully observable. Therefore, all state variables are observable in the proposed self-calibration scheme.

4. Simulation and Experiment

4.1. Simulation and Analysis

The simulations are first implemented to validate the effectiveness of the proposed self-calibration method. The random error of the gyro was set at 0.005 deg/h, and the random error of the accelerometers was 30 μ g , the measurement noise from gyros and accelerometers had the magnitude of 0.2 d e g / h and 0.3 μ g , respectively. The IMU was in a dynamic process during the calibration, it rotated at an angular rate up to 20 degrees per second, which was a much greater magnitude than the measurement noise. Thus, the measurement noise of IMU can be ignored. The error of the angle of the rotation mechanism was 0.5 , while the error of the angular rate of rotation mechanism was 0.7 / s . The initial location was 114.19659 ° E , 30.93041 ° N . The angular rate of rotation mechanism was set as 20   d e g / s   (72,000   d e g / h ). The covariance Q was set as d i a g w 2 w 2 , w = 1 , the dimension of Q was 12, and R in the two Kalman filters were set as d i a g v 1 2 v 1 2 v 1 2 ,   v 1 = 1 μ g and d i a g v 2 2 v 2 2 v 2 2 ,   v 2 = 1 d e g / h , respectively. Measurement data of the IMU in simulation is shown in Figure 7 and Figure 8. The real values of installation errors and their estimation error are shown in Table 1. Figure 9, Figure 10 and Figure 11 describe the comparisons of estimation results of all installation errors, the blue lines represent the estimation of installation errors, and the red dotted lines represent real values.
The time in Figure 9, Figure 10 and Figure 11 was relative, which indicated the convergent tendency of installation errors when the errors were fully activated and could be calibrated. From Table 1 and Figure 9, Figure 10 and Figure 11, we can see all installation errors present good convergence with the proposed method. The estimation errors are all less than 4%.

4.2. Experiment

To validate the effectiveness of the proposed method, laboratory experiments were carried out with an actual dual-axis RINS which consists of fiber gyros and quartz accelerometers, whose detailed specifications of dual-axis RINS are shown in Table 2. The covariance Q and R were set to the same as in the simulations. The installation accuracy between the inner axis and outer axis of the dual-axis RINS used in this paper was better than ± 4 arc-seconds, thus the inner axis is approximately regarded as perpendicular to the outer axis.
The dual-axis RINS was statically placed on a turntable to isolate external disturbance while the IMU rotates along with the rotation mechanism inside the dual-axis RINS to conduct the calibration.
Considering that the error parameter’s reference cannot be acquired in real situations, estimation results of installation errors using two separating calibration methods mentioned in [11] and [13] were regarded as a reference to prove that the proposed self-calibration method can get comparable accuracy, thereby avoiding the complex separating process. The internal structure of the dual-axis RINS is shown in Figure 12.
The experimental procedure is shown in Figure 13. The first period of rotation of the IMU was conducted at first, and Equation (11) was used to estimate nonorthogonal errors of the accelerometers and installation errors between the IMU and rotation axes. Then the results obtained from the first step were sent to the next calibration process. The second period of rotation of the IMU was conducted, and Equation (13) was used to get nonorthogonal errors of gyros. In the end, estimations from the proposed method were compared with counterparts from the separating methods [11,13].
The reference of installation errors and their estimation errors are shown in Table 3. Measurement data of the IMU in the experiment is shown in Figure 14 and Figure 15. Figure 16, Figure 17 and Figure 18 describe the comparisons of the estimation results of all installation errors, the blue lines represent the estimation results of installation errors, and the red line represents the reference.
The time in Figure 16, Figure 17 and Figure 18 was relative, which indicated the convergent tendency of installation errors when the errors were fully activated and could be calibrated. From Table 3 and Figure 16, Figure 17 and Figure 18, we can see that all estimation results present good convergence with the proposed method. All installation errors can be estimated together by the proposed method, which acquires comparable results compared with the separating methods. The estimation errors were all less than 6%.

5. Conclusions

A system-level self-calibration method for installation errors of a dual-axis RINS was proposed in this paper. The measurement model was reestablished to ensure that all errors can be estimated together, rather than calibrated separately by utilizing different estimation algorithms and rotation schemes. The relationship between initial IMU attitude and subsequent IMU attitude and angular rate of the rotational mechanism were used to estimate two types of installation errors of RINS, that is, the nonorthogonal installation errors of inertial sensors and installation errors between the IMU axes and rotation axes. The rotation scheme of the IMU was discussed and designed, and the PWCS method was used to analyze the observability of the system. Finally, simulation and laboratory experiments were conducted to validate the proposed method. The results showed that all installation errors of the dual-axis RINS can be accurately estimated by the proposed system-level self-calibration method, thereby avoiding complex separating calibration process.

Author Contributions

Conceptualization, S.B. and J.L.; methodology, S.B. and P.L.; software, S.B. and X.X.; validation, S.B. and M.L.; formal analysis, S.B. and K.H.; investigation, S.B.; writing—original draft preparation, S.B.; writing—review and editing, J.L. and P.L.

Funding

This research was funded by National Natural Science Foundation of China (61703207, 61973160); Jiangsu Provincial Natural Science Foundation of China (BK20170801); Aeronautical Science Foundation of China (2018ZC52037).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yuan, B.; Liao, D.; Han, S. Error Compensation of an Optical Gyro INS by Multi-Axis Rotation. Meas. Sci. Technol. 2012, 23, 025102. [Google Scholar] [CrossRef]
  2. Zhang, Q.; Wang, L.; Liu, Z.; Feng, P. An Accurate Calibration Method Based on Velocity in a Rotational Inertial Navigation System. Sensors 2015, 15, 18443–18458. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. 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, Aberdeen, UK, 18–21 June 2007. [Google Scholar]
  4. Gao, J.; Li, K.; Chen, Y. Study on Integration of Fog Single-Axis Rotational INS and Odometer for Land Vehicle. IEEE Sens. J. 2017, 18, 752–763. [Google Scholar] [CrossRef]
  5. Yin, H.; Yang, G.; Song, N.; Jiang, R.; Wang, Y. Error Modulation Scheme Analyzing for Dual-Axis Rotating Fiber-Optic Gyro Inertial Navigation System. Sens. Lett. 2012, 10, 1359–1363. [Google Scholar] [CrossRef]
  6. Ben, Y.; Chai, Y.; Gao, W.; Sun, F. Analysis of Error for a Rotating Strap-Down Inertial Navigation System with Fibro Gyro. J. Mar. Sci. Appl. 2010, 9, 419–424. [Google Scholar] [CrossRef]
  7. Lv, P.; Lai, J.; Liu, J.; Nie, M. The Compensation Effects of Gyros’ Stochastic Errors in a Rotational Inertial Navigation System. J. Navig. 2014, 67, 1069–1088. [Google Scholar] [CrossRef]
  8. Dai, M.; Lu, J. A Full-Parameter Self-Calibration Method Based on Inertial Frame Filtering for Triaxis RINS Under Swaying Base. IEEE Sens. J. 2019, 19, 2170–2180. [Google Scholar] [CrossRef]
  9. Wang, L.; Wang, Y.; Zhang, Q.; Gao, P. Self-Calibration Method Based on Navigation in High-Precision Inertial Navigation System with Fiber Optic Gyro. Opt. Eng. 2014, 53, 064103. [Google Scholar] [CrossRef]
  10. Gao, P.; Li, K.; Wang, L.; Liu, J. A Self-Calibration Method for Accelerometer Nonlinearity Errors in Triaxis Rotational Inertial Navigation System. IEEE Trans. Instrum. Meas. Sci. Technol. 2017, 66, 243–253. [Google Scholar]
  11. 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. 2015, 62, 2353–2362. [Google Scholar] [CrossRef]
  12. 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]
  13. Song, N.; Cai, Q.; Yang, G.; Yang, G. Analysis and Calibration of the Mounting Errors between Inertial Measurement Unit and Turntable in Dual-Axis Rotational Inertial Navigation System. Meas. Sci. Technol. 2013, 24, 115002. [Google Scholar] [CrossRef]
  14. Kang, L.; Ye, L.; Song, K.; Zhou, Y. Attitude Heading Reference System Using Mems Inertial Sensors with Dual-Axis Rotation. Sensors 2014, 14, 18075–18095. [Google Scholar] [CrossRef] [PubMed]
  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. Gao, P.; Li, K.; Wang, L.; Liu, Z. A Self-Calibration Method for Tri-Axis Rotational Inertial Navigation System. Meas. Sci. Technol. 2016, 27, 115009. [Google Scholar] [CrossRef]
  17. Cao, Y.; Cai, H.; Zhang, S.; Li, A. A New Continuous Self-Calibration Scheme for a Gimbaled Inertial Measurement Unit. Meas. Sci. Technol. 2011, 23, 015103. [Google Scholar] [CrossRef]
Figure 1. Definition of the inertial measurement unit (IMU) frame, accelerometers frame, and gyro frame. (a) IMU frame and accelerometer frame; (b) IMU frame and gyro frame.
Figure 1. Definition of the inertial measurement unit (IMU) frame, accelerometers frame, and gyro frame. (a) IMU frame and accelerometer frame; (b) IMU frame and gyro frame.
Sensors 19 04005 g001
Figure 2. Definition of IMU frame and rotational mechanism frame in the initial state.
Figure 2. Definition of IMU frame and rotational mechanism frame in the initial state.
Sensors 19 04005 g002
Figure 3. Structure of system-level self-calibration method.
Figure 3. Structure of system-level self-calibration method.
Sensors 19 04005 g003
Figure 4. One period of rotation in the self-calibration scheme.
Figure 4. One period of rotation in the self-calibration scheme.
Sensors 19 04005 g004
Figure 5. Angle changes of the inner and outer axis during the whole calibration. (a) Inner axis angle; (b) Outer axis angle.
Figure 5. Angle changes of the inner and outer axis during the whole calibration. (a) Inner axis angle; (b) Outer axis angle.
Sensors 19 04005 g005
Figure 6. Rank change of the system.
Figure 6. Rank change of the system.
Sensors 19 04005 g006
Figure 7. Measurement data of gyros in simulation (a) x-axis; (b) y-axis; (c) z-axis.
Figure 7. Measurement data of gyros in simulation (a) x-axis; (b) y-axis; (c) z-axis.
Sensors 19 04005 g007
Figure 8. Measurement data of accelerometers in simulation (a) x-axis; (b) y-axis; (c) z-axis.
Figure 8. Measurement data of accelerometers in simulation (a) x-axis; (b) y-axis; (c) z-axis.
Sensors 19 04005 g008
Figure 9. Estimation of accelerometers nonorthogonal installation errors. (a) Sayz; (b) Sazx; (c) Sazy.
Figure 9. Estimation of accelerometers nonorthogonal installation errors. (a) Sayz; (b) Sazx; (c) Sazy.
Sensors 19 04005 g009
Figure 10. Estimation of installation errors between rotation axes and IMU axes. (a) α ; (b) β ; (c) γ .
Figure 10. Estimation of installation errors between rotation axes and IMU axes. (a) α ; (b) β ; (c) γ .
Sensors 19 04005 g010
Figure 11. Estimation of gyros nonorthogonal installation errors. (a) Sgxy; (b) Sgxz; (c) Sgyx; (d) Sgyz; (e) Sgxz; (f) Sgyx.
Figure 11. Estimation of gyros nonorthogonal installation errors. (a) Sgxy; (b) Sgxz; (c) Sgyx; (d) Sgyz; (e) Sgxz; (f) Sgyx.
Sensors 19 04005 g011
Figure 12. Internal structure of the dual-axis rotational inertial navigation system (RINS).
Figure 12. Internal structure of the dual-axis rotational inertial navigation system (RINS).
Sensors 19 04005 g012
Figure 13. The flow diagram of the experiment.
Figure 13. The flow diagram of the experiment.
Sensors 19 04005 g013
Figure 14. Measurement data of gyros in experiment (a) x-axis; (b) y-axis; (c) z-axis.
Figure 14. Measurement data of gyros in experiment (a) x-axis; (b) y-axis; (c) z-axis.
Sensors 19 04005 g014
Figure 15. Measurement data of accelerometers in experiment (a) x-axis; (b) y-axis; (c) z-axis.
Figure 15. Measurement data of accelerometers in experiment (a) x-axis; (b) y-axis; (c) z-axis.
Sensors 19 04005 g015
Figure 16. Estimation of accelerometers nonorthogonal installation errors. (a) Sayz; (b) Sazx; (c) Sazy.
Figure 16. Estimation of accelerometers nonorthogonal installation errors. (a) Sayz; (b) Sazx; (c) Sazy.
Sensors 19 04005 g016
Figure 17. Estimation of installation errors between rotation axes and IMU axes. (a) α ; (b) β ; (c) γ .
Figure 17. Estimation of installation errors between rotation axes and IMU axes. (a) α ; (b) β ; (c) γ .
Sensors 19 04005 g017
Figure 18. Estimation of gyros nonorthogonal installation errors. (a) Sgxy; (b) Sgxz; (c) Sgyx; (d) Sgyz; (e) Sgxz; (f) Sgyx.
Figure 18. Estimation of gyros nonorthogonal installation errors. (a) Sgxy; (b) Sgxz; (c) Sgyx; (d) Sgyz; (e) Sgxz; (f) Sgyx.
Sensors 19 04005 g018
Table 1. Estimation results.
Table 1. Estimation results.
ParametersReal Values (’’)Estimation Errors (’’)
S g x y 160.389
S g x z 160.575
S g y x 160.474
S g y z 160.488
S g z x 160.544
S g z y 160.406
S a y z 100.255
S a z x 100.164
S a z y 100.119
α 1000.145
β 1000.042
γ 1000.067
Table 2. Specifications of the dual-axis rotational inertial navigation system (RINS).
Table 2. Specifications of the dual-axis rotational inertial navigation system (RINS).
CharacteristicsDescription
Sampling frequency200 Hz
Gyro bias stability0.005 ° / h
Accelerometers bias stability30   μ g
Rotation angle accuracy0.5
Rotation velocity accuracy0.7 / s
Table 3. Estimation results.
Table 3. Estimation results.
ParametersReal Values (’’)Estimation Errors (’’)
S g x y −233.6150.540
S g x z −195.1801.054
S g y x −79.4670.573
S g y z 39.6030.698
S g z x −421.0752.294
S g z y −84.6600.559
S a y z 31.3730.538
S a z x −36.9111.610
S a z y −31.4901.884
α −13.8280.299
β 111.3831.422
γ −247.5181.502

Share and Cite

MDPI and ACS Style

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. https://0-doi-org.brum.beds.ac.uk/10.3390/s19184005

AMA Style

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(18):4005. https://0-doi-org.brum.beds.ac.uk/10.3390/s19184005

Chicago/Turabian Style

Bai, Shiyu, Jizhou Lai, Pin Lyu, Xiaowei Xu, Ming Liu, and Kai Huang. 2019. "A System-Level Self-Calibration Method for Installation Errors in A Dual-Axis Rotational Inertial Navigation System" Sensors 19, no. 18: 4005. https://0-doi-org.brum.beds.ac.uk/10.3390/s19184005

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