Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Pilot Study of Individual Muscle Force Prediction during Elbow Flexion and Extension in the Neurorehabilitation Field
Next Article in Special Issue
A Novel Gravity Compensation Method for High Precision Free-INS Based on “Extreme Learning Machine”
Previous Article in Journal
One-Step Fabrication of Microchannels with Integrated Three Dimensional Features by Hot Intrusion Embossing
Previous Article in Special Issue
Anatomical Calibration through Post-Processing of Standard Motion Tests Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MEMS IMU Error Mitigation Using Rotation Modulation Technique

1
School of Aeronautics and Astronautics, University of Electronics Science and Technology of China, Chengdu 610000, China
2
School of Geomatics, Liaoning Technical University, Fuxin 123000, China
3
Department of Geomatics Engineering, The University of Calgary, Calgary, AB T2N 1N4, Canada
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(12), 2017; https://doi.org/10.3390/s16122017
Submission received: 6 August 2016 / Revised: 26 October 2016 / Accepted: 28 October 2016 / Published: 29 November 2016
(This article belongs to the Special Issue Inertial Sensors and Systems 2016)

Abstract

:
Micro-electro-mechanical-systems (MEMS) inertial measurement unit (IMU) outputs are corrupted by significant sensor errors. The navigation errors of a MEMS-based inertial navigation system will therefore accumulate very quickly over time. This requires aiding from other sensors such as Global Navigation Satellite Systems (GNSS). However, it will still remain a significant challenge in the presence of GNSS outages, which are typically in urban canopies. This paper proposed a rotary inertial navigation system (INS) to mitigate navigation errors caused by MEMS inertial sensor errors when external aiding information is not available. A rotary INS is an inertial navigator in which the IMU is installed on a rotation platform. Application of proper rotation schemes can effectively cancel and reduce sensor errors. A rotary INS has the potential to significantly increase the time period that INS can bridge GNSS outages and make MEMS IMU possible to maintain longer autonomous navigation performance when there is no external aiding. In this research, several IMU rotation schemes (rotation about X-, Y- and Z-axes) are analyzed to mitigate the navigation errors caused by MEMS IMU sensor errors. As the IMU rotation induces additional sensor errors, a calibration process is proposed to remove the induced errors. Tests are further conducted with two MEMS IMUs installed on a tri-axial rotation table to verify the error mitigation by IMU rotations.

1. Introduction

Tens of years ago, to maintain long term autonomous navigation performance, the inertial navigation system (INS) was designed based on high end inertial sensors such as ring laser gyro (RLG) and fiber optical gyro (FOG). However, the large size (usually very heavy as well) and high cost have limited their use to various applications. With the development of micro-electro-mechanical system (MEMS) technology, the MEMS-based inertial sensors became available to the commercial market in 1990s and they quickly found applications for indoor pedestrian and robotic navigation because of their low-cost, small-size, light-weight and low-power consumption [1,2]. MEMS IMU outputs, however, are corrupted with significant sensor errors, such as high frequency noise, bias, scale factors and installation errors. As a result, the navigation errors will accumulate quickly and deteriorate the navigation solution over a short time period [2,3,4,5]. For example, the position errors of a low-cost MEMS IMU will grow to kilometers within several minutes. To limit the accumulation of navigation errors in INS, aiding from other sensors is required. For example, GNSS position and velocity can be used to estimate the INS errors using an extended Kalman filter [6,7,8,9]. Various investigations have been conducted to employ external navigation information to limit the INS error accumulation. The neural networks are applied to mimic the behavior of navigation error accumulation of INS [10,11], and the fuzzy logic is also employed to adaptively adjust the parameters of the estimation algorithm [12]. The auto-regressive model was also proposed to interpret the randomness of the MEMS-based inertial sensor errors [3,13].
This paper investigates a different approach to limit the accumulation of navigation errors in MEMS-based INS. Instead of focusing on the use of external aiding information by existing methods, it uses a rotary INS to reduce the navigation error growth by rotating the IMU. A rotary INS is an inertial navigator in which the IMU is installed on a rotation platform. By rotating the IMU in certain pre-defined schemes, it can change the characteristics of the inertial sensor errors and consequently reduce the accumulation of navigation errors. For example, the rotation of IMU with a constant angular rate can modulate the constant inertial bias into periodic signals and an integration of the modulated inertial data over a complete rotation cycle can eliminate the bias impact on the navigation solutions. Rotary INS technique has been widely applied in marine navigation for submarines and warships [14,15,16]. However, only few researches have been conducted on rotary INS based on low-cost MEMS IMU, which has great potential applications in the near future and is the focus of this paper. A rotary MEMS-based INS has the potential to significantly increase the time period that it can bridge GNSS outages for GNSS/MEMS IMU integrated systems [17,18]. It also makes MEMS IMU possible to maintain much longer autonomous navigation performance as a self-contained navigation system when there is no external aiding available.
Although the requirement for a rotational platform would increase the complexity and cost of the use of the rotation modulation technique, it will be advantageous to attach the MEMS IMU to a rotating platform that is already built in the vehicle, such as the wheel for land vehicle, the airscrew for the helicopter, and the propeller for the boat, just to mention a few. In a recent research, gyros are attached to the wheel of a ground vehicle and the results have demonstrated that a low-cost MEMS IMU can provide a very accurate navigation solution using the rotation modulation technique [18]. In the near future, the advances in hardware technologies would also reduce the complexity and cost in developing rotational IMU devices, creating increased applications. In addition to the mitigation of navigation error accumulations in a non-aiding mode, the rotation of IMU will help improve the system observability as the observability of inertial sensor errors are related to IMU orientations [1]. Overall, the rotary MEMS-based inertial system is expected to significantly improve the navigation performance comparing to a non-rotary one.
Comparing to the low-cost MEMS IMU, FOG and RLG errors are orders of magnitude smaller, and thus the application of the rotation modulation technique to FOG and RLG based IMU can effectively mitigate the error accumulations. However for MEMS IMU, due to their significant bias instability, scale factor, installation errors, as well as noise, how can the application of the rotation modulation technique help reduce the navigation errors is not clear so far which requires investigation. Although the idea of MEMS-based rotary INS had already been proposed, it was verified only by simulations in previous research [17]. Real rotated tests must be conducted in order to fully investigate the feasibility of MEMS-based rotary systems for practical applications. This is the main objective of this contribution with the following specific research efforts:
(1)
Error analysis for IMU rotation about X- and Y-axes: As previous research indicated that the IMU rotation about the Z-axis cannot modulate the errors in the Z-axis, which results in accumulated azimuth and velocity errors, this research conducts error analysis with IMU rotations about the X- and Y-axes, and compares three different rotation schemes for navigation applications.
(2)
The calibration of gyro scale factor and installation errors: As the gyro scale factor and installation errors induce additional sensor errors when IMU rotates, which will deteriorate the error migration performance, a process to calibrate gyro scale factor and installation errors is proposed for the MEMS-based inertial system.
(3)
Experiments and tests with real rotated data: As MEMS IMU outputs contain significant bias instability error, the sensor biases will vary randomly over time, not behave as constants. Tests with real rotated data are necessary in order to study the feasibility of MEMS-based rotary inertial system in real scenario.
The remainder of this paper is organized as follows. Section 2 describes the rotary INS, including the mechanization algorithm and the error mitigations with different IMU rotation schemes, each rotating about X-, Y- and Z-axes, respectively. The related mathematical equations are also derived and provided in this section. As the IMU rotation induces additional sensor errors, a calibration process is proposed and introduced in Section 3. The tests of two different MEMS IMUs on a tri-axial rotation table are given to verify the proposed calibration algorithm and error mitigation performance using rotary INS in Section 4. The conclusions and future work are summarized in Section 5.

2. Rotary Inertial Navigation System

The concept of rotary INS was initially proposed for the gimbaled inertial system by [19]. Afterwards, this technique was applied to a strapdown inertial system [19,20]. As certain rotations of IMU can mitigate the INS navigation errors without the requirement of external information, the rotary INS was extensively employed for the warship or submarine. The rotary INS normally contains two components, an IMU and a rotational platform, as shown in Figure 1. In the sequel, INS with non-rotating IMU is referred to as conventional INS to distinguish it from rotary INS.

2.1. Rotary INS Mechanization

As the IMU is rotating in the rotary INS, a new frame in which the inertial readings are collected is introduced, in addition to the traditional coordinate frames in conventional INS. This new coordinate frame can be referred as inertial sensor frame or IMU frame, and its axis are aligned with the sensitive axis of inertial sensors with the origin defined as the origin of IMU.
Generally speaking, the mechanization algorithm of rotary INS is very similar to the one of conventional INS. As the inertial readings are collected in sensor frame, a transformation process is required, as shown in Equations (1) and (2).
f i b b = C s b ( f i s s + f s b s )
ω i b b = C s b ( ω i s s + ω s b s )
where f i b b and ω i b b are the specific force and angular rate in body frame with respect to inertial frame, respectively; C s b is the transformation matrix from senor frame to body frame; f i s s and ω i s s are the specific force and angular rate in sensor frame with respect to inertial frame, respectively; and f s b s and ω s b s are the specific force and angular rate of body frame with respect to sensor frame; respectively. As the rotation of IMU does not introduce any linear movement, f s b s is a zero vector, while ω s b s is related to the IMU rotation schemes.
With the transformed specific force and angular rate in body frame, the mechanization of the conventional INS can be used to derive position, velocity and attitude solutions as shown in Figure 2, where C b n represents the transformation matrix from the body frame to the navigation frame, n represents the navigation frame, Ω represents the skew-symmetric matrix of angular rate ω, V represents the velocity, and φ , λ , h represent the latitude, longitude, and height, respectively. Apparently, the rotation angle between body frame and sensor frame is required in the transformation process, and this is usually measured by a device installed in rotation platform.

2.2. Error Mitigations with Different IMU Rotations

The inertial sensor errors are firstly introduced and the mitigation of navigation errors is then investigated using three different IMU rotation schemes, each rotating about X-, Y- and Z-axes, respectively. All analyses are conducted in static mode and the body frame is assumed to be aligned with the navigation frame to simplify the analysis.

2.2.1. Inertial Sensor Errors

The specific forces and angular rates collected by the accelerometer and gyro triads contain different types of errors, such as noise, turn on biases, scale factor errors and installation errors. The sensor errors are the major error source that causes navigation errors. Usually, the sensor error model for gyros and accelerometers can be described by Equations (3) and (4), respectively [21].
ω ˜ = ω + d + S g ω + N g ω + n w
f ˜ = f + b + S a f + N a f + n f
where ω ˜ is the actual gyro outputs, ω is the true gyro angular rate, d is the gyro bias, S g is the scale factor, N g is the installation errors, n ω is the gyro noise, f ˜ is the actual accelerometer outputs, f is the true specific force, b is the accelerometer bias, S a is the scale factor, N a is the installation errors, and n f is the gyro noise.
The gyro and accelerometer biases normally contain two parts, namely, deterministic part and random part. The deterministic part is constant but it differs each time when the sensor is turned on, known as turn on bias in some literatures. The random part varies with time and it is quantified by the bias instability [4,13,22]. The high-end IMU features very small gyro bias instability, e.g., 0.01°/h or even lower, but for MEMS IMU it is very significant, e.g., 10°/h–100°/h or even higher.
The scale factors of gyro and accelerometer cause sensor errors from the true value. Usually the scale factors for gyros and accelerometers are represented by the diagonal matrices in Equations (5) and (6), respectively, and they are quantified in unit of part per million (ppm). Normally the scale factor for MEMS IMU varies as a function of surrounding environments (such as temperature) [21,23].
S g = [ K g x 0 0 0 K g y 0 0 0 K g z ]
S a = [ K a x 0 0 0 K a y 0 0 0 K a z ]
where K g x , K g y , K g z are the scale factors for gyros along X-, Y- and Z-axes, respectively, and K a x , K a y , K a z are the scale factors for accelerometers along X-, Y- and Z-axes, respectively.
Due to the imperfection of inertial sensor assembling and installation, the three sensitive axes of the sensor triad are not perfectly orthogonal with each other, which cause the sensed inertial value of one axis to project into the two other axes. The installation error for gyros and accelerometers can be described by Equations (7) and (8), respectively, and they are quantified using the unit of angle. Similarly, the installation errors are also temperature-dependent [21].
N g = [ 0 K g x y K g x z K g y x 0 K g y z K g z x K g z y 0 ]
N a = [ 0 K a x y K a x z K a y x 0 K a y z K a z x K a z y 0 ]
where K g i j is the installation error between i-axis and j-axis, and K a i j is the accelerometer installation error between i-axis and j-axis (i, j = x, y, z).
Although the sensor biases, scale factors and installation errors vary in nature, they are considered to be constants in a short time period under stable temperature, e.g., during a complete rotation cycle, to simplify the analysis.

2.2.2. Error Mitigation by Rotating IMU about X Axis

As the IMU rotates about its X-axis with a rate of ω, the transformation matrix between the body frame and the sensor frame can be described by Equation (9).
C b s = [ 1 0 0 0 cos ω t sin ω t 0 sin ω t cos ω t ] = ( C s b ) T
By applying the above transformation matrix, the gyro and accelerometer biases in the navigation frame can be described by Equations (10) and (11), respectively [24,25]. These errors in north and vertical directions are modulated into periodic signals, and the attitude and velocity errors caused by such errors are self-eliminated after a complete rotation cycle as shown in Equations (12) and (13). As the gyro and accelerometer bias in the rotation axis cannot be modulated, the attitude and velocity error in the east direction propagates in the same way as in the conventional INS.
d n = C b n C s b d s = [ d x s d y s cos ω t + d z s sin ω t d y s sin ω t + d z s cos ω t ]
b n = C b n C s b b s = [ b x s b y s cos ω t + b z s sin ω t b y s sin ω t + b z s cos ω t ]
0 T d n d t = [ T d x s 0 0 ]
0 T b n d t = [ T b x s 0 0 ]
where d s = [ d x s d y s d z s ] T and b s = [ b x s b y s b z s ] T are the gyro drift and accelerometer bias in IMU frame, respectively; d n = [ d E d N d U ] T and b n = [ b E b N b U ] T are the gyro drift and accelerometer bias in navigation frame, respectively; and C b n is the transformation matrix from body frame to navigation frame.
Although IMU rotation can modulate the sensor constant biases, it also induces additional gyro biases proportional to the rotation rate due to the gyro scale factor of the X-axis as shown in Equation (14). After a complete rotation cycle, the rotation-induced gyro bias results in the accumulated attitude errors as shown in Equation (15). Similarly, the IMU rotation rate will be projected to the Y- and Z-axes due to gyro installation errors, and leads to additional gyro biases as shown in Equation (16). Although the rotation-induced gyro biases can be modulated, and the resulted attitude errors are removed after a complete rotation cycle as shown in Equation (17), it still degrades the navigation solutions within the rotation cycle [17,26].
δ ω S F n = C s n S g ω i s s = [ K g x ω K g y ( ω i e cos φ 1 + cos 2 ω t 2 + ω i e sin φ sin 2 ω t 2 ) K g z ( ω i e cos φ 1 cos 2 ω t 2 + ω i e sin φ sin 2 ω t 2 ) K g y ( ω i e cos φ sin 2 ω t 2 + ω i e sin φ 1 cos 2 ω t 2 ) + K g z ( ω i e cos φ sin 2 ω t 2 + ω i e sin φ 1 + cos 2 ω t 2 ) ]
0 T δ ω S F n d t = [ T K g x ω 1 2 T ( K g y + K g z ) ω i e cos φ 1 2 T ( K g y + K g z ) ω i e sin φ ]
δ ω N n = C s n N g ω i s s = [ K g x y ( ω i e cos φ cos ω t + ω i e sin φ sin ω t ) + K g x z ( ω i e cos φ sin ω t + ω i e sin φ cos ω t ) K g y x ω cos ω t K g z x ω sin ω t + K g y z ω i e ( sin φ 1 + cos 2 ω t 2 cos φ sin 2 ω t 2 ) K g z y ω i e ( cos φ sin 2 ω t 2 + sin φ 1 cos 2 ω t 2 ) K g y x ω sin ω t + K g z x ω cos ω t + K g y z ω i e ( sin φ sin 2 ω t 2 cos φ 1 cos 2 ω t 2 ) + K g z y ω i e ( cos φ 1 + cos 2 ω t 2 + sin φ sin 2 ω t 2 ) ]
0 T δ ω N n = [ 0 1 2 ( K g y z K g z y ) ω i e sin φ 1 2 ( K g y z + K g z y ) ω i e cos φ ]
where δ ω S F n is the sensor error caused by the scale factor, ω i e is the earth rotation rate in the navigation frame, and δ ω N n is the gyro error caused by the installation error.

2.2.3. Error Mitigation by Rotating IMU about Y and Z Axis

When IMU rotates about the Y-axis, both gyro and accelerometer biases in X- and Z-axes are modulated, so the attitude and velocity errors caused by those errors are automatically removed in east and vertical directions after a complete rotation cycle, while the errors in north direction propagate in the same way as in conventional INS, as shown in Equations (18) and (19). Similarly, the IMU rotation also induces the errors in the Y-axis due to the scale factor of the rotation axis, and the errors in the X- and Z-axes due to the gyro installation errors, as shown in Equations (20) and (21), respectively [26].
d n = C b n C s b d s = [ d x s cos ω t d z s sin ω t d y s d x s sin ω t + d z s cos ω t ]
b n = C b n C s b b s = [ b x s cos ω t b z s sin ω t b y s b x s sin ω t + b z s cos ω t ]
δ ω S F n = [ K g x ω i e sin φ sin 2 ω t 2 + K g z ω i e sin φ sin 2 ω t 2 K g y ( ω i e cos φ + ω ) K g x ω i e sin φ 1 cos 2 ω t 2 + K g z ω i e sin φ 1 + cos 2 ω t 2 ]
δ ω N n = [ K g x y ( ω i e cos φ + ω ) cos ω t + K g x z ω i e sin φ cos 2 ω t K g z x ω i e sin φ sin 2 ω t + K g z y ( ω i e cos φ + ω ) sin ω t K g y x ω i e sin φ sin ω t + K g y z ω i e sin φ cos ω t K g x y ( ω i e cos φ + ω ) sin ω t K g x z ω i e sin φ sin 2 ω t 2 K g z x ω i e sin φ sin 2 ω t 2 + K g z y ( ω i e cos φ + ω ) cos ω t ]
The gyro and accelerometer biases in X- and Y-axes are modulated when IMU rotates about Z-axis, therefore the attitude and velocity errors caused by such errors are self-eliminated in the east-north plane after a complete rotation cycle, though the errors in vertical direction still propagate in the same way as in conventional INS, as shown in Equations (22) and (23). Because of the gyro scale factor in the Z-axis and the installation errors, the IMU rotation also induces gyro biases, which can be described by Equations (24) and (25) [17,26].
d n = C b n C s b d s = [ d x s cos ω t + d y s sin ω t d x s sin ω t d y s cos ω t d z s ]
b n = C b n C s b b s = [ b x s cos ω t + b y s sin ω t b x s sin ω t b y s cos ω t b z s ]
δ ω S F n = [ ( K g x K g y ) ω i e cos φ sin ω t cos ω t ( K g x sin 2 ω t + K g y cos 2 ω t ) ω i e cos φ K g z ( ω i e sin φ + ω ) ]
δ ω N n = [ ω i e cos φ ( K g x y cos 2 ω t K g y x sin 2 ω t ) + ( ω i e sin φ + ω ) ( K g x z cos ω t K g y z sin ω t ) ω i e cos φ sin 2 ω t ( K g x y + K g y x ) / 2 + ( ω i e sin φ + ω ) ( K g x z sin ω t + K g y z cos ω t ) K g z x ω i e cos φ sin ω t + K g z y ω i e cos φ cos ω t ]
Table 1 summarizes the velocity errors caused by accelerometer biases and the attitude errors caused by gyro biases, gyro scale factors and gyro installation errors after a complete rotation cycle for IMU rotations about the X-, Y- and Z-axes.
For the IMU rotation about the X-axis, the unmodulated accelerometer bias causes the linear velocity errors and quadratic position errors over time in the east direction, as shown in Equation (26). The unmodulated gyro bias in the X-axis results in accumulated attitude error in the east direction, which causes projection error of the local gravity on the north direction. Eventually, it leads to quadratic velocity errors and cubic position errors over time as shown in Equation (27) [3,6]. Although the attitude errors caused by gyro scale factors can be found in the east, north and vertical directions, the east component is much more significant than other two components, as the IMU rotation rate is usually much greater than the earth rotation rate. Similar to the gyro biases, the rotation-induced gyro errors (due to the scale factor) in the east direction also leads to the quadratic velocity errors and cubic position errors in the north direction. As the rotation-induced gyro errors in the Y- and Z-axes due to the installation errors are modulated after a complete rotation cycle, the effect of attitude errors on the velocity and position errors are limited.
Similar to the IMU rotation about the X-axis, the unmodulated accelerometer bias in the Y-axis leads to the linear velocity errors and quadratic position errors in the north direction, while the unmodulated gyro bias and scale factor in the rotation axis result in the quadratic velocity errors and cubic position errors in the east direction, for the IMU rotation about the Y-axis.
For the IMU rotation about the Z-axis, the unmodulated accelerometer bias causes the linear velocity and quadratic position errors in the vertical direction, while the unmodulated gyro bias and scale factor of the rotation axis results in the azimuth errors, which lead to the horizontal velocity and position errors. For the IMU rotation about the X or Y-axis, as the velocity and position errors in the east or north directions accumulate quadratically and cubically over time, they are mostly much greater than the ones for the IMU rotation about the Z-axis, which is related to the coupling of the azimuth errors and the vehicle dynamics (usually can be expressed by a sine function times the vehicle velocity or displacement). Therefore, the IMU rotation about the Z-axis is more suitable than other two rotation schemes for navigation applications. The rest of the paper will concentrate on the use of this rotation scheme.
δ V E = b x s d t = b x s t , δ P E = δ V E d t = b x s t d t = 1 2 b x s t 2
δ V N = g ε E d t = g d x s t d t = 1 2 g d x s t 2 , δ P N = δ V N d t = 1 2 g d x s t 2 d t = 1 6 g d x s t 3
Other conclusions are summarized as follows:
(1)
The IMU rotation modulates the constant biases of inertial sensors that perpendicular to the rotation axis, and the attitude and velocity errors caused by such biases are self-eliminated after a complete rotation cycle;
(2)
The constant biases of inertial sensors in the rotation axis cannot be modulated, and the attitude and velocity errors caused by such errors propagate in the same way as in conventional INS;
(3)
The IMU rotation induces an extra error in the gyro of the rotation axis due to gyro scale factor, and this error results in accumulated attitude errors in the direction of the corresponding rotation axis;
(4)
The IMU rotation also induces extra errors in the gyros that are perpendicular to the rotation axis due to gyro installation errors, resulting in attitude and velocity errors.
(5)
As the rotation-induced errors are proportional to the rotation rate, which is usually much more significant than earth rotation rate, a calibration process is required to remove the rotation induced errors.

3. Calibration for the MEMS IMU

As the IMU rotation induces additional gyro biases, which eventually leads to navigation errors [17,26], a calibration process is proposed for rotary INS with IMU rotation about Z-axis. Although the full gyro scale factor and installation errors are represented by 9 parameters according to Equations (5) and (6), the proposed method only calibrate the gyro scale factor of Z-axis, K g z , and the installation errors, K g x z , K g y z , which motivate significant gyro errors when IMU rotates about Z-axis.

3.1. Gyro Error Model for Calibration

As the IMU rotation does not introduce any linear motion of IMU, the theoretical gyro outputs in body frame can be described by Equation (28), when body frame is static relative to local level frame.
ω i b b = C n b ω i b n = C n b ω i e n = C n b [ 0 ω i e cos φ ω i e sin φ ] = [ p ω i e q ω i e r ω i e ]
where and they satisfy the equality equation p 2 + q 2 + r 2 = 1 .
The theoretical and actual gyro outputs in the sensor frame can be described by Equations (29) and (30), respectively, and the gyro errors are given in Equation (31).
ω i s s = C b s ω i b b + ω b s s = [ ω i e ( p cos ω t + q sin ω t ) ω i e ( p sin ω t + q cos ω t ) r ω i e + ω ]
ω ˜ i s s = ω i s s + N ω i s s + d s
δ ω i s s = [ S g x ω i e ( p cos ω t + q sin ω t ) + K g x y ω i e ( p sin ω t + q cos ω t ) + K g x z ( r ω i e + ω ) + d x s K g y x ω i e ( p cos ω t + q sin ω t ) + S g y ω i e ( p sin ω t + q cos ω t ) + K g y z ( r ω i e + ω ) + d y s K g z x ω i e ( p cos ω t + q sin ω t ) + K g z y ω i e ( p sin ω t + q cos ω t ) + S g z ( r ω i e + ω ) + d z s ]
where N = [ S g x K g x y K g x z K g y x S g y K g y z K g z x K g z y S g z ] is a 3 × 3 matrix represents the combination of gyro scale factor and installation errors and d s = [ d x s d y s d z s ] T is the gyro biases in sensor frame.
The earth rotation rate can be ignored in the above error model because: (1) the earth rotation rate cannot be sensed in most MEMS IMU due to their significant gyro bias instability and noise; and (2) the IMU rotation rate is much more significant than earth rotation rate. By ignoring the earth rotation rate, the above gyro error model can be simplified as shown in Equation (32).
δ ω i s s = [ K g x z ω + d x s K g y z ω + d y s S g z ω + d z s ]
With different IMU rotation rates, the least square algorithm is applied to estimate the gyro scale factor, K g z , the gyro installation errors, K g x z , K g y z , as well as the gyro biases, [ d x s d y s d z s ] T , based on the error model given in Equation (32). The measurements are the gyro readings which can be described by Equation (33), and the measurement model for a single measurement is given by Equation (34).
Z j = ω ˜ i s s ω i s s [ ω ˜ i s , x s ω ˜ i s , y s ω ˜ i s , z s ω j ]
Z j = H j X = [ ω j 0 0 1 0 0 0 ω j 0 0 1 0 0 0 ω j 0 0 1 ] [ K g x z K g y z S g z d x s d y s d z s ]
where Z j represent the j th measurement, and ω j represents the IMU rotation rate with respect to the Z j .

3.2. Implementation of Calibration Process

The calibration process can be described as two steps: (1) the tested IMU remains still on the rotation platform for 60 s; and (2) the IMU rotates along with the rotation platform at a designated rotation rate about Z-axis for 60 s. The total calibration time is 2 min. According to the error model in Equation (32), the gyro biases can be estimated from the inertial data collected in static period, while the gyro scale factor, K g z , and installation errors, K g x z , K g y z , can be estimated from the collected data in rotation period. It is worthy to mention that, the MEMS IMU usually features significant non-linearity errors [27,28], so the designated rotation rate in calibration process should be the same as the one employed in rotation scheme. For example, for the MEMS-based rotary INS with IMU rotation about Z-axis at the rate of 10°/s, the calibration process should employ the same rotation rate to estimate the gyro scale factor and installation errors.
The error sources for the calibration error include the gyro sensor noise, the gyro bias instability (gyro biases are not constants but time-correlated variables), and the ignored earth rotation rate. The calibration errors of gyro biases caused by noise can be calculated using Equation (35), and the calibration errors of scale factor and installation errors caused by noise are given by Equation (36) [21]. Apparently, the longer calibration time can reduce the calibration errors, however it may increase the effect of gyro bias instability on calibration results. The reason we choose the 60 s for calibration time in each step will be further discussed in Section 4 along with the noise level and gyro bias instability of tested IMU.
σ b = A R W / T s t a t i c
σ s = A R W T r o t a t i o n / α r o t a t i o n
where ARW is the angular random walk of gyros, T s t a t i c is the time period that IMU remains still, T r o t a t i o n is the time period that IMU rotation, and α r o t a t i o n is the rotated angle of IMU.
It should be noted that the calibration method above is proposed to estimate gyro errors before each use of the rotary system. It requires the vehicle remains still during the calibration process, and this limits its use for on-line calibration of the rotary system. For on-line calibration, the external measurements are required, and an extended Kalman filter (EKF) is usually employed to estimate the sensor errors.

4. Testing MEMS IMU on a Tri-Axial Rotation Table

By using a tri-axial rotation table as the rotation platform, two MEMS IMUs, namely MTi-G and NAV440 are tested. The rotation platform consists of a tri-axial rotation table and a console (computer) as shown in Figure 3. The tri-axial rotation table has three rotational frames, namely, outer frame, middle frame and inner frame. The console controls the position and rotation of these frames.
The MEMS IMUs are firmly installed on a piece of metal underneath the inner frame by screws as shown in Figure 4. Different rotation schemes, such as rotation about X-, Y- and Z-axes can be implemented by rotating the three frames. For example, with the IMU axis defined as shown in Figure 5, the rotation of outer frame rotates the IMU about its Z-axis when both the middle and inner frames remain at the level position (or angle position of 0°). An initialization process, after which both middle and inner frames are in the level position and the rotation axis of inner frame points to north direction, is required for the rotation table. Apparently, with IMU installations as shown in Figure 5, the IMU body frame is aligned to local level frame after the initialization process.
Table 2 and Table 3 summarize the technical parameters of the rotation platform and the characteristics of the tested IMU, respectively. All data are provided from the manufactures and available on their official website. Although the gyro white noise of NAV440 is not provided, the Allan variance method was utilized to obtain this value, which is about 0.4°/ h [13].
The frame rotation angle data (IMU rotation angle between body frame and sensor frame) and inertial data are recorded by the console with data rates of 50 Hz and 100 Hz, respectively. The frame rotation angle will be used to calculate transformation matrix between body frame and sensor frame, which are necessary to solve the attitude solutions.
The Allan variance method was also employed to study the random variations of the sensor biases as shown in Figure 6, which presents the calculated Allan variance of gyro data in the X-axis for both tested IMUs. The cluster time (or averaging time) ranges from 10 s to 300 s, and the most significant bias variance was observed at the cluster time of 10 s for MTi-G due to its high noise level. As the cluster time increases, the bias variances are reduced as the noises are eliminated by the averaging process. For both IMUs, the bias variances reach their minimums at the cluster time of around 120 s (about 24°/h and 11°/h for MTi-G and NAV440, respectively). For the longer cluster time, the bias variances are fluctuated due to the rate random walk, the exponentially correlated noise and the sinusoidal noise [29]. For the gyro data in the Y- and Z-axes, similar results are obtained.
The conventional INS static tests and the rotary INS static tests are conducted using both MEMS IMUs in a laboratory. The conventional INS static tests, in which the inertial data are collected while the MEMS IMU remains still on the rotation table for 6 min (the data from the 1st minute will be used to derive the gyro bias estimates), are conducted to investigate how fast the inertial errors accumulate with time in a non-rotary system and to provide a comparison to the rotary INS. The rotary INS static tests are conducted to study the feasibility of the MEMS-based rotary system, and to verify the inertial error mitigation performance through IMU rotations. Although higher IMU rotation rates can more effectively mitigate the navigation errors in theory [24,25], the effect of IMU rotation rate on navigation errors needs to be investigated for a MEMS-based rotary system, as the MEMS IMU features significant bias instability and noise level. Ten individual tests are conducted for the rotary INS static tests. In each individual test, the IMU rotates about Z-axis at a designated rate with a calibration process conducted in the beginning. The designated rotation rates for the ten individual tests are 10°/s, 20°/s, 30°/s, 40°/s, 50°/s, 60°/s, 70°/s, 80°/s, 90°/s and 100°/s, respectively. The time length of the rotated data for each individual test is 5 min (not including the calibration time), as after that the position errors are already accumulated to hundreds or thousands of meters. As mentioned in Section 3, the IMU rotation rate employed in the calibration process is the same as the one in the corresponding individual test.

4.1. Results of Conventional INS Static Tests

The means of the gyro data in three axes collected in the 1st minute provide the estimates of the gyro biases [30]. Based on the initial position, velocity and attitude from external information (body frame is aligned to local level frame after the initialization of rotation table), the navigation solutions are calculated based on the collected inertial data. As expected, the position, velocity and attitude errors accumulate fast with time. The Root mean square (RMS) values of the navigation errors for MTi-G and NAV440 are summarized in Table 4 and Table 5, respectively. After 5 min, the horizontal position errors accumulate to several kilometers for both MEMS IMUs.

4.2. Results of Rotary INS Static Tests

In each individual test, the gyro scale factor of Z-axis, K g z , the installation errors, K g x z , K g y z , and the gyro biases are estimated for MTi-G and NAV440 in the calibration process. Based on the gyro noise level of tested IMUs given in Table 3, the calibration errors caused by gyro noise for both IMUs are calculated and given in Table 6. It should be noted that the unit of installation errors are converted to ppm, and the calibration errors of scale factor and installation errors can also be converted to degree per hour. With the time length of 60 s, the calibration errors of gyro biases caused by the noises are 23°/h and 3°/h for MTi-G and NAV440, respectively, and similar calibration errors are also obtained for scale factor and installation errors.
The major calibration error sources contain the gyro noises and the gyro bias random variations. Although longer calibration time reduces calibration errors caused by the noise, it also increases the effect of gyro bias variations on calibration results. Therefore, the choosing of calibration time is a compromise of the two factors. For both tested IMUs, the calibration time is chosen as 120 s (60 s in each step) and the reason is that the calibration errors caused by the gyro noise are comparable to the minimum bias variations for the cluster time of 120 s as shown in Figure 6.
Table 7 summarizes the mean of the gyro data in both IMUs after the calibration process for each individual test. As we can see, the mean values of the gyro data in the X- and Y-axes are close to zero and the means of the gyro data in the Z-axis are close to the rotation rate after the calibration, which prove the rotation induced gyro biases and turn on biases are removed.
The rotary INS results are derived with two different data processing strategies, namely, partially calibrated data processing and fully calibrated data processing, to verify the effect of gyro scale factor and installation errors on navigation errors. For the partially calibrated data processing, the rotary navigation solutions are derived with only the estimates of gyro biases from the calibration results, so the gyro scale factor and installation errors remain in the rotary inertial data. For the fully calibrated data processing, the estimates of gyro installation errors, scale factor and biases are applied to derive the rotary INS solutions.

4.2.1. MTi-G Results

By using two data processing methods, the obtained roll and pitch errors for the individual test with rotation rate of 10°/s are given in Figure 7. For the partially calibrated data processing, the attitude errors are modulated into oscillating signals because of the modulation of the gyro errors in X- and Y-axes through IMU rotation, and the oscillating period is the same as the rotation period. As the gyro biases are already removed, the remaining errors mostly consist of bias residuals and rotation-induced errors due to gyro installation errors. Similar characteristics can be observed in the horizontal velocity errors. For the fully calibrated data processing, the oscillating amplitudes are greatly reduced as shown in the figure. This is because the removal of gyro installation errors eliminates the effect of rotation-induced errors on navigation solutions.
For each individual test, the RMS values of the position, velocity and attitude errors are calculated using partially and fully data processing strategies, as shown in Figure 8 and Figure 9. Comparing to the errors in the conventional INS tests, the horizontal position, velocity, roll and pitch errors are reduced by modulation accelerometer and gyro errors in X- and Y-axes through IMU rotation in the partially calibrated case, although the gyro scale factor and installation errors are not removed yet. The major error sources of the navigation errors are gyro bias residuals (as the gyro biases vary with time) and the induced gyro biases by IMU rotation. As shown in Figure 9, the IMU rotation induces the gyro bias in Z-axis, which leads to azimuth error almost linearly increases over the rotation rate. The mean values of the RMS errors for horizontal position and velocity for all tested rotation rates are 1525.8 m and 14.6 m/s, respectively.
For the fully calibrated data processing, the removal of the rotation-induced gyro errors in the X- and Y-axes further reduce the roll and pitch errors, as well as the horizontal position and velocity errors, comparing to the results in the partially calibrated case. The azimuth errors are also significantly reduced because of the estimation of the gyro scale factor in the Z-axis. The mean value of RMS for the horizontal position errors and the velocity errors for all tests are further reduced to 1057.2 m and 10.4 m/s, respectively. The roll and pitch errors are found much smaller than the azimuth errors. This is because the gyro errors in the X- and Y-axes are modulated by IMU rotation about the Z-axis, while such errors in the Z-axis cannot be modulated. Although a higher IMU rotation rate can more effectively modulate sensor biases and mitigate error accumulations for FOG or RLG-based system [24,25], it can be noted that the navigation errors are not strongly related to the IMU rotation rate for MTi-G. This is because the MTi-G is a very low-cost MEMS IMU, which features significant bias instability (e.g., 3600°/h, as shown in Table 3) and high noise level (e.g., 3.0°/ h , as shown in Table 3). As a result, they affect the error mitigation performance in rotary system.

4.2.2. Nav440 Results

The rotated Nav440 data are also processed using the two data processing methods, and the obtained roll and pitch errors for the individual test with rotation rate of 30°/s are given in Figure 10. Similarly, the attitude errors are modulated into oscillating signals through IMU rotation, and the oscillating amplitude is greatly reduced in the fully calibrated case. Based on Figure 10, we can see that the greater oscillating amplitude in the partially calibrated case is mainly caused by rotation-induced errors due to gyro installation errors, while the much smaller oscillating amplitude in the fully calibrated case is caused by gyro bias residuals. Similar characteristics can be observed in horizontal velocity errors.
By using two different data processing methods, the obtained RMS values of navigation errors for each individual test are given in Figure 11 and Figure 12. For the partially calibrated data processing, the horizontal position and velocity errors as well as pitch and roll errors are reduced compared to the ones in conventional test, though the azimuth error still grows almost linearly over the IMU rotation rate due to the rotation-induced gyro errors in Z-axis. The mean values of the RMS errors for horizontal position and velocity are 544.8 m and 4.9 m/s, respectively. With the removal of gyro installation errors and scale factors in fully calibrated data processing, all navigation errors are further reduced, and the mean of RMS horizontal position and velocity errors are dropped to 213.8 m and 2.2 m/s, respectively.
It can be noted that the magnitude of the roll and pitch errors are nearly one order smaller than the azimuth errors due to the modulation of gyro errors in the X- and Y-axes. As the NAV440’s gyro bias instability and noise level are much lower than MTi-G, much smaller navigation errors are obtained for the rotary system based on NAV440.
Based on the results from the conventional INS static tests and the rotary INS static tests with two different MEMS IMUs, we can see that IMU rotations can modulate the sensor errors, and mitigate their effect on navigation solutions. The gyro scale factor and installation errors will cause additional navigation errors by IMU rotations, and their calibration is therefore necessary for a MEMS-based rotary system in order to achieve better navigation performance. The rotation rate is not the only factor that affects the error modulation, which will also be affected by the sensor errors, such as the bias instability and the noise. The results also indicate that the sensor errors can be more efficiently modulated for the MEMS IMUs that have smaller bias instability and lower noise level.

5. Conclusions

This paper investigates the mitigation of navigation errors due to significant MEMS inertial sensor errors by rotating IMU. Three IMU rotation schemes, each rotating about X-, Y- and Z-axes, respectively, are investigated. The related mathematical equations are derived and analyzed for each rotation scheme. A calibration method is proposed to remove the sensor biases, as well as gyro scale factors and installation errors for a MEMS-based rotary system. The tests with two MEMS IMUs on a tri-axial rotation table are conducted to verify the feasibility of MEMS-based rotary INS. The results indicate that with efficient data processing strategy, the IMU rotation can effectively modulate sensor errors and mitigate their effects on navigation solutions for a low-cost MEMS IMU. Kinematic field tests and other factors, such as vibration and wobble, that may affect the error mitigation in rotary INS, will be investigated in the future research.

Acknowledgments

This study is partly supported by Sichuan Province Science and Technology Project (Grant No. 2016GZ0062), National Natural Science Foundation of China (Grant No. 41604025), and Jiangxi Province Key Lab for Digital Land (Grant No. DLLJ201501).

Author Contributions

Shuang Du and Wei Sun conceived the methodology and designed the prototype system, as well as designed and performed the experiments; Shuang Du processed the data; Shuang Du, Wei Sun and Yang Gao analyzed the data; and Shuang Du and Yang Gao wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hide, C.; Moore, T.; Hill, C.; Abdularhim, K. Investigating the use of rotating foot mounted inertial sensor for positioning. In Proceedings of the 25th ION GNSS, Nashville, TN, USA, 17–21 September 2012.
  2. Zhou, Z.; Li, Y.; Zhang, J.; Rizos, C. Integrated Navigation System for a Low-Cost Quadrotor Aerial Vehicle in the Presence of Rotor Influences. J. Surv. Eng. 2016, 4, 1–13. [Google Scholar] [CrossRef]
  3. Nassar, S. Improving the Inertial Navigation System (INS) Error Model for INS and INS/DGPS Applications. Ph.D. Thesis, The University of Calgary, Calgary, AB, Canada, 2003. [Google Scholar]
  4. Nassar, S.; Schwarz, K.P.; El-Sheimy, N. INS and INS/GPS Accuracy Improvement Using Autoregressive Modeling of INS Sensor Errors. In Proceedings of the ION NTM, San Diego, CA, USA, 26–28 January 2003.
  5. Shin, E.; Niu, X.; El-Sheimy, N. Performance Comparison of the Extended and the Unscented Kalman Filter for Integrated GPS and MEMS-Based Inertial Systems. In Proceedings of the ION NTM, San Diego, CA, USA, 24–26 January 2005.
  6. Shin, E. Estimation Techniques for Low-Cost Inertial Navigation. Ph.D. Thesis, University of Calgary, Calgary, AB, Canada, 2005. [Google Scholar]
  7. Antonio, A. GNSS/INS Integration Methods. Ph.D. Thesis, The University of Calgary, Calgary, AB, Canada, 2010. [Google Scholar]
  8. Fang, J.; Going, X. Predictive iterated Kalman Filter for INS/GPS integration and its application to SARmotion compensation. IEEE Trans. Instrum. Meas. 2010, 6, 909–915. [Google Scholar] [CrossRef]
  9. Du, S.; Gao, Y. Inertial aided cycle slip detection and identification for integrated PPP GPS and INS. Sensors 2012, 12, 14344–14362. [Google Scholar] [CrossRef] [PubMed]
  10. Kaygisiz, B.H.; Erkmen, A.M.; Erkmen, I. GPS/INS enhancement using neural networks for autonomous ground vehicle applications. In Proceedings of the 2003 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS 2003), Las Vegas, NV, USA, 27–31 October 2003.
  11. Noureldin, A.; El-Shafie, A.; Bayoumi, M. GPS/INS integration utilizing dynamic neural networks for vehicular navigation. Inf. Fusion 2011, 1, 48–57. [Google Scholar] [CrossRef]
  12. Noureldin, A.; El-Shafie, A.; El-Sheimy, N. Adaptive neurofuzzy module for inertial navigation system/global positioning system integration utilising position and velocity updates with realtime cross-validation. Radar Sonar Navig. IET 2007, 1, 388–396. [Google Scholar] [CrossRef]
  13. Du, S.; Gao, Y. Integration of Precise Point Positioning and the Latest MEMS IMU for Precise Applications. In Proceedings of the ION GNSS 2012, Nashville, TN, USA, 18–21 September 2012.
  14. Levinson, D.E.; Willcocks, M. The next generation marine inertial navigation is here now. In Proceedings of the IEEE Position Location and Navigation Symposium, Las Vegas, NV, USA, 11–15 April 1994.
  15. Levinson, D.E.; Majure, R. Accuracy enhancement techniques applied to the Marine Ring Laser Inertial Navigator (MARLIN). J. Inst. Navig. 1987, 34, 71–80. [Google Scholar] [CrossRef]
  16. Heckman, D.W.; Baretela, L.M. Improved Affordability of High Precision Submarine Inertial Navigation by Insertion of Rapidly Developing Fiber Optic Gyro Technology. In Proceedings of the IEEE PLANS, San Diego, CA, USA, 13–16 March 2000.
  17. Sun, W.; Wang, D.X.; Xu, L.W.; Xu, L.L. MEMS-based rotary strapdown inertial navigation system. Measurement 2013, 46, 2585–2596. [Google Scholar] [CrossRef]
  18. Collin, J. MEMS IMU Carouseling for ground vehicles. IEEE Trans. Veh. Technol. 2015, 64, 2242–2251. [Google Scholar] [CrossRef]
  19. Geller, E.S. Inertial system platform rotation. IEEE Trans. Aerosp. Electr. Syst. 1968, aes-4, 557–568. [Google Scholar] [CrossRef]
  20. Giovanni, S.C.; Levinson, D.E. Performance of a Ring Laser Strapdown Marine Gyrocompass. Navigation 1981, 28, 311–341. [Google Scholar] [CrossRef]
  21. Niu, X.; Li, Y.; Zhang, H.P.; Wang, Q.J.; Ban, Y.L. Fast Thermal Calibration of Low-Grade Inertial Sensors and Inertial Measurement Units. Sensors 2013, 13, 12192–12217. [Google Scholar] [CrossRef] [PubMed]
  22. Hou, H. Modeling Inertial Sensor Errors Using Allan Variance. Master’s Thesis, The University of Calgary, Calgary, AB, Canada, 2004. [Google Scholar]
  23. Guo, J.; Zhong, M.Y. Calibration and Compensation of the Scale Factor Errors in DTG POS. IEEE Trans. Instrum. Meas. 2013, 62, 2784–2794. [Google Scholar] [CrossRef]
  24. Ben, Y.Y.; Chai, Y.L.; Gao, W.; Sun, F. Analysis of Error for a Rotating Strap-down Inertial Navigation System with Fiber Gyro. J. Mar. Sci. Appl. 2010, 9, 419–424. [Google Scholar] [CrossRef]
  25. Sun, F.; Sun, W.; Gao, W.; Ben, Y.Y. Research on the Technology of Rotational Motion for FOG Strapdown Inertial Navigation System. In Proceedings of the 2009 IEEE International Conference on Mechatronics and Automation, Changchun, China, 9–12 August 2009.
  26. Du, S.; Sun, W.; Gao, Y. An Investigation on MEMS IMU Error Mitigation Using Rotation Modulation Technique. In Proceedings of the ION GNSS, Tempa, FL, USA, 8–12 September 2014.
  27. Artese, G.; Trecroci, A. Calibration of a low cost MEMS INS sensor for an integrated navigation system. In Proceedings of the International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Beijing, China, 3–4 September 2008.
  28. He, K.P. Parametric Identification and Error Compensation of MEMS Inertial Sensors and IMU. Ph.D. Thesis, Harbin Engineering University, Harbin, China, 2008. [Google Scholar]
  29. Hou, H.; El-Sheimy, N. Inertial Sensor Errors Modeling Using Allan Variance. In Proceedings of the ION GNSS 2003, Portland, OR, USA, 9–12 September 2003.
  30. Skalund, J. Optimizing Georeferencing of Airborne Survey Systems by INS/DGPS. Ph.D. Thesis, The University of Calgary, Calgary, AB, Canada, 1999. [Google Scholar]
Figure 1. Structure of rotary inertial navigation system (INS).
Figure 1. Structure of rotary inertial navigation system (INS).
Sensors 16 02017 g001
Figure 2. Flowchart of rotary INS mechanization.
Figure 2. Flowchart of rotary INS mechanization.
Sensors 16 02017 g002
Figure 3. Tri-axial rotation platform.
Figure 3. Tri-axial rotation platform.
Sensors 16 02017 g003
Figure 4. Installation of Micro-Electron-Mechanical System (MEMS) IMU on tri-axial rotation table.
Figure 4. Installation of Micro-Electron-Mechanical System (MEMS) IMU on tri-axial rotation table.
Sensors 16 02017 g004
Figure 5. Rotation table set-up for IMU rotation about Z-axis.
Figure 5. Rotation table set-up for IMU rotation about Z-axis.
Sensors 16 02017 g005
Figure 6. Allan variances of gyro data in X-axis for MTi-G and NAV440.
Figure 6. Allan variances of gyro data in X-axis for MTi-G and NAV440.
Sensors 16 02017 g006
Figure 7. Roll and pitch errors for rotary INS with MTi-G.
Figure 7. Roll and pitch errors for rotary INS with MTi-G.
Sensors 16 02017 g007
Figure 8. Root mean square (RMS) values of horizontal position and velocity errors for rotary INS with different IMU rotation rates using MTi-G.
Figure 8. Root mean square (RMS) values of horizontal position and velocity errors for rotary INS with different IMU rotation rates using MTi-G.
Sensors 16 02017 g008
Figure 9. RMS values of attitude errors for rotary INS with different IMU rotation rates using MTi-G.
Figure 9. RMS values of attitude errors for rotary INS with different IMU rotation rates using MTi-G.
Sensors 16 02017 g009
Figure 10. Roll and pitch errors for rotary INS with NAV440.
Figure 10. Roll and pitch errors for rotary INS with NAV440.
Sensors 16 02017 g010
Figure 11. RMS values of horizontal position and velocity errors for rotary INS with different IMU rotation rates using NAV440.
Figure 11. RMS values of horizontal position and velocity errors for rotary INS with different IMU rotation rates using NAV440.
Sensors 16 02017 g011
Figure 12. RMS values of attitude errors for rotary INS with different IMU rotation rates using NAV440.
Figure 12. RMS values of attitude errors for rotary INS with different IMU rotation rates using NAV440.
Sensors 16 02017 g012
Table 1. Errors after a complete rotation cycle for IMU rotation about X-, Y- and Z-axes.
Table 1. Errors after a complete rotation cycle for IMU rotation about X-, Y- and Z-axes.
Errors After a Complete Rotation CycleRotation About X AxisRotation About Y AxisRotation About Z Axis
Accelerometer biases [ T b x s 0 0 ] [ 0 T b y s 0 ] [ 0 0 T b z s ]
Gyro biases [ T d x s 0 0 ] [ 0 T d y s 0 ] [ 0 0 T d z s ]
Gyro scale factors [ T K g x ω 1 2 T ( K g y + K g z ) ω i e cos φ 1 2 T ( K g y + K g z ) ω i e sin φ ] [ 0 T K g y ( ω i e cos φ + ω ) T 2 ( K g x + K g z ) ω i e sin φ ] [ 0 T 2 ( K g x + K g y ) ω i e cos φ K g z ( ω i e sin φ + ω ) T ]
Gyro installation errors [ 0 1 2 ( K g y z K g z y ) ω i e sin φ 1 2 ( K g y z + K g z y ) ω i e cos φ ] [ T 2 ( K g x z K g z x ) ω i e sin φ 0 0 ] [ T 2 ( K g x y K g y x ) ω i e cos φ 0 0 ]
Table 2. Technical parameters of tri-axial rotation platform.
Table 2. Technical parameters of tri-axial rotation platform.
Position Accuracy (°)Rotation Rate Accuracy (°/s)Maximum Rotation Rate (°/s)
1 × 10−51 × 10−5±100
Table 3. Characteristics of tested IMU.
Table 3. Characteristics of tested IMU.
CharacteristicsMTi-GNAV440
Range (°/s)−300~300−400~400
Gyro bias instability (°/h)360020
Gyro white noise (°/ h )3.0-
Table 4. RMS of navigation errors for conventional INS with MTi-G in static mode.
Table 4. RMS of navigation errors for conventional INS with MTi-G in static mode.
Latitude (m)Longitude (m)Height (m)
2412.31526.8173.8
Vn (m/s)Ve (m/s)Vu (m/s)
33.220.51.2
Pitch (°)Roll (°)Azimuth (°)
1.831.322.31
Table 5. RMS of navigation errors for conventional INS with NAV440 in static mode.
Table 5. RMS of navigation errors for conventional INS with NAV440 in static mode.
Latitude (m)Longitude (m)Height (m)
829.12820.11033.0
Vn (m/s)Ve (m/s)Vu (m/s)
4.826.79.0
Pitch (°)Roll (°)Azimuth (°)
0.480.600.73
Table 6. Calibration errors caused by gyro noise for tested IMUs.
Table 6. Calibration errors caused by gyro noise for tested IMUs.
Rotation Rate (°/s)MTi-GNAV440
σ b (°/h) σ s (ppm) σ b (°/h) σ s (ppm)
1023.2645.53.186.1
20322.743.0
30215.228.7
40161.421.5
50129.117.2
60107.614.2
7092.212.3
8080.710.8
Table 7. Mean of gyro data after calibration for each individual rotary INS test.
Table 7. Mean of gyro data after calibration for each individual rotary INS test.
Rotation Rate (°/s)MTi-GNAV440
X (°/s) Y (°/s)Z (°/s)X (°/s)Y (°/s)Z (°/s)
10−2.0 × 10−31.0 × 10−210.022−2.0 × 10−3−1.1 × 10−29.990
20−7.1 × 10−31.4 × 10−320.0209.8 × 10−3−1.1 × 10−319.990
30−1.4 × 10−31.9 × 10−430.0247.9 × 10−39.6 × 10−429.996
40−5.9 × 10−4−7.8 × 10−340.0088.7 × 10−32.2 × 10−339.999
50−4.8 × 10−3−2.5 × 10−350.0012.0 × 10−3−2.3 × 10−350.011
60−1.5 × 10−31.7 × 10−360.0084.5 × 10−4−2.5 × 10−360.005
701.1 × 10−23.8 × 10−369.992−2.9 × 10−4−2.7 × 10−370.004
803.3 × 10−3−1.6 × 10−379.988−1.1 × 10−34.7 × 10−380.002
90−8.3 × 10−3−1.6 × 10−389.998−3.3 × 10-35.2 × 10−390.005
1009.4 × 10−33.1 × 10−399.993−5.6 × 10−33.7 × 10−3100.008

Share and Cite

MDPI and ACS Style

Du, S.; Sun, W.; Gao, Y. MEMS IMU Error Mitigation Using Rotation Modulation Technique. Sensors 2016, 16, 2017. https://doi.org/10.3390/s16122017

AMA Style

Du S, Sun W, Gao Y. MEMS IMU Error Mitigation Using Rotation Modulation Technique. Sensors. 2016; 16(12):2017. https://doi.org/10.3390/s16122017

Chicago/Turabian Style

Du, Shuang, Wei Sun, and Yang Gao. 2016. "MEMS IMU Error Mitigation Using Rotation Modulation Technique" Sensors 16, no. 12: 2017. https://doi.org/10.3390/s16122017

APA Style

Du, S., Sun, W., & Gao, Y. (2016). MEMS IMU Error Mitigation Using Rotation Modulation Technique. Sensors, 16(12), 2017. https://doi.org/10.3390/s16122017

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