Abstract

With the rapid development of inertial technology, the rotational inertial system has been widely used for high-accuracy navigation. In addition, the single-axis rotational inertial navigation system is one of the most popular navigation systems due to its low cost. However, the fact which cannot be ignored is that although single axis rotational inertial navigation system can restrain divergence of attitude error, it cannot eliminate velocity error, which is mainly caused by the gyro misalignment and scale factor error related to rotational axis. A novel calibration method established on filter has been proposed. Position, velocity, and attitude are chosen to be the state variables. In order to estimate the gyro errors related to rotational axis, these errors must be included in the state equation. Correspondingly, zero velocity and rate of turntable are chosen to be measured. Simulation has been carried out to verify the correctness of the theory, and the real test has been performed to further demonstrate the validity of the method. By compensating the main error estimated by filter, it can be found that the accumulated errors in velocity and attitude are decreased. So, the precision of navigation is greatly improved with the proposed method.

1. Introduction

Inertial navigation system (INS) can obtain position, attitude, and velocity of the craft by resolving the data sampled by its inertial measurement unit, which is IMU and contains three orthogonal gyros and accelerometers [1]. Thus, the INS is self-contained and widely used in military field such as airplanes, submarines, and ships. INS can be divided into two categories: platform and strap-down inertial navigation system. Compared with platform INS, strap-down INS has the characteristics of higher reliability, smaller volume, better maintenance, and lower cost [24].

Device error is affected by the precision of the navigation. Meanwhile, velocity and attitude errors will accumulate with time quickly [5]. There are two major approaches to enhance the precision of inertial navigation. One is to improve the performance of gyro and accelerometer. However, this method is expensive and restricted by technique. The other one is error compensation. Compensation can be performed based on the external information such as Global Positioning System (GPS) and turntable [6]. Self-compensation is another way to eliminate the device error. Rotation modulation technology has been one of the most widely used self-compensation methods, especially in the field of aerospace for its low cost and reliability.

In rotational inertial navigation system (RINS), IMU is placed on the indexing mechanization which is fixed on the craft [79]. Thus, the self-compensation can be performed by designing appropriate rotating sequence. RINS can be divided into single-axis, dual-axis, and triaxial-axis or above [1012].

The advantage of rotation modulation is that it can realize self-compensation of sensor errors during navigation process [13, 14]. However, errors still exist and affect the precision of the navigation.

In triaxis RINS, instrument errors, such as accelerometer nonlinearity and inner lever-arm, can be accurately compensated through optimal estimation with velocity and position error measurements. Meanwhile, proper rotation scheme design can help modulate and estimate error parameters of IMU. Taking commercial cost into consideration, triaxis and dual-axis RINSs are more expensive than single-axis [1518]. In consequence, single-axis RINS will be mainly discussed in this paper.

In single-axis RINS, gyro drift is the main concern. Shang et al. [19] pointed out that single-axis rotation can modulate the constant drift perpendicular to the rotational axis, while gyro drift on the spin axis still diverge with time. In order to deal with the constant drift of spin axis, identification methods of gyro drift have been proposed. Yu et al. [20] calibrated the single-axis gyro drift by using radial basis function network. A least square method was established to calibrate the axial gyro drift precisely by Hu et al. [21]. Similar methods of calibration of gyro drifts are shown in [2224].

Besides the gyro drift, error propagation of the single axis rotational INS showed that scale factor error and misalignment errors related to the rotational axis cannot be modulated [25, 26]. The scale factor error along the rotational axis was the main error that affected the single-axis rotational INS in [27]. However, no paper can give answer to this question. Therefore, a novel method which can calibrate the instrument errors is proposed in this article.

Compared with other methods, the method proposed in this paper is suitable for single-axis RINS to calibrate main device errors by using zero velocity and attitude errors. Thus, compensation can be made to reduce the error accumulation during the navigation process. Based on this, the precision of the single-axis RINS can be promoted.

The advantages of the proposed method are shown as follows: according to single-axis rotational inertial navigation system, this paper provides the deduction of error propagation in terms of misalignment and scale factor error related to rotational axis. Besides, it points out that in single-axis RINS, attitude error caused by gyro misalignment related to rotational axis can be modulated, but velocity error still exists. Simulation has been carried out to verify this thesis. By observing the velocity and rate of turntable, gyro misalignment and scale factor errors of the rotational axis can be estimated under stationary base.

The paper is organized as follows. Section 2 provides the error mechanism of single axis rotational INS. Section 3 introduces the design of filter for single-axis rotational system. Section 4 describes the process of the simulation, which is carried out to verify the error mechanism and the proposed method. In order to further demonstrate the proposed method, experiment has been performed and results are shown in Section 5. Conclusion is drawn in Section 6.

2. Error Mechanism

The denotation of frames is shown in Figure 1, where denotes earth-centered inertial frame, denotes earth-centered frame, denotes navigation frame, and denotes body frame.

2.1. Wander Azimuth Frame

The wander frame is widely used in navigation system, and it can be obtained by rotating with respect to the local geodetic frame z axis by the wander frame angle α.

In this paper, considering the fact calibration is performed on stationary base, α is chosen to be zero. The error model in local geodetic frame can be considered the same as wander azimuth frame in a short time.

2.2. Error Analysis

To further demonstrate the error mechanism in this paper, simplification and assumption are given.

The IMU can be rotated about the vertical body axis with constant angular velocity :

The gyro angular velocity error equation is shown above, where denotes the measured angular velocity in the body frame, denotes the true angular velocity, denotes the angular velocity error, denotes gyro bias, denotes the scale factor error, and denotes the misalignment:where denotes the angular velocity provided by turntable, which rotates around the body spin axis with constant rate and denotes the earth angular velocity in the body frame.

Gyro error model is given by equation (1). Considering the angular velocity provided by turntable is much larger than the earth angular velocity and gyro bias, the approximation of gyro model is shown as follows:

The simplified gyro error model shows that the most important error sources are , , . Correspondingly, the velocity error accumulated by inertial device error will directly affect the precision of navigation. The rotation sequence from navigation frame to body frame is yaw, pitch, and roll, respectively:where is the direction cosine matrix of body frame, is used to transfer the initial body frame to the wander azimuth navigation frame, denotes the initial yaw angle, denotes the initial roll angle,, denotes the initial pitch angle, and denotes rotation angle around spin axis from initial frame to body frame .

The attitude error can be defined as (concrete deduction is shown in Appendix)where denotes the attitude error at the beginning and denotes the accumulated attitude error. From equation (6), it can be seen that the accumulated leveling attitude error can be modulated by rotation. However, the initial attitude error still exists and has great effect on the velocity error.

Besides, scale factor error of the rotational axis gyro will affect yaw angle. The velocity error dynamic equations demonstrate the relationship between attitude error and velocity error as follows:

Velocity errors can be obtained by the integral of equation (7). To research the factors that affect velocity, expansion of integral is shown as follows (concrete deduction is given in Appendix):

From the final expansion of equation (8), it can be seen that terms in the first bracket are converged because trigonometric function is bounded. Terms in the second bracket are diverging with time. Thus, the drift of velocity is mainly caused by the items in the second bracket that contains initial attitude error, gyro misalignment, and scale factor error related to rotational axis.

Similarly, the error velocity can be expressed as

Obviously, the main angular velocity error above, which is induced by scale factor error and misalignment error related to body spin axis, cannot be modulated during the rotation. The navigation attitude error will accumulate quickly because of this. Thus, the precision of the INS navigation system can be improved by using the self-calibration method to estimate the errors related to body spin axis.

3. Design of Filter for Single-Axis Rotational System

According to the error mechanism, the proposed method should estimate the instrument errors including scale factor errors and misalignment of the related rotation axis. Besides, the navigation information, position, velocity, and attitude should also be contained in the measurement matrix. Correspondingly, the measurement model of the proposed method should be established by parameters that can be obtained precisely.

3.1. State Equation

The state model chosen in this paper is listed as follows:where , , and denote position errors, , , and denote velocity errors, , , and denote attitude errors, and , , and denote deterministic unknown errors.(1)Instrument error equationsThe gyro errors are modelled as constant during the filtering process:(2)Attitude error equations:where denotes dynamic error, denotes gyro error under navigation system, denotes the navigation to earth frame angular rotation rate in the navigation frame, and denotes the Earth’s rotation rate.(3)Velocity error equations:where denotes the sensed output of i-axis accelerometer, and denote the angular position errors, denotes the altitude error, and denotes the error of accelerometer.(4)Position error equations:where denotes the earth radius.(5)Gyro error model:where denotes the angular velocity error caused by gyro, denotes the bias of i-axis gyro, denotes the scale factor error of i-axis gyro, and denotes the misalignment error of gyro in i-axis related to j-axis.

3.2. Measurement Model

The angular velocity provided by rate turntable is much larger than that of the earth, so the errors are mainly caused by scale factor errors and misalignment errors related to rotational axis. Thus, the error of gyro output and the velocity during navigation process can be utilized as observations to estimate the corresponding errors.(1)Given that the IMU is tested on stationary base, the velocity of IMU should be zero. The output velocity of IMU can be defined as the measurement :The measurement equation of sensed velocity of tested IMU can be described aswhere the measurement matrix is(2)The attitude error always exists and the computed attitude matrix translating from navigation frame to body frame is represented aswhere denotes tilt error and denotes attitude direction cosine matrix without attitude error. The input rate of turntable and the rate of earth rotation can be considered as references, and the output rate of IMU can be measured by gyros. Thus, the measurement is defined to be the difference between the measured angular velocity and the true angular velocity:and the measurement matrix of can be expressed as(3)Flow chart of self-calibration filtering is given in Figure 2.The principle of the proposed calibration method in this paper to estimate the unknown errors is shown in Figure 2. Figure 2 uses angular velocity and zero velocity as the observed quantity. A 12-state Kalman filter is established by the INS dynamic equations and simplified gyro error model.

4. Simulation

To support the theory discussed above and evaluate the performance of the proposed method, simulations are carried out in Section 4.2.

4.1. Condition

The IMU errors in this simulation are defined as follows: gyro bias is 0.1°/h with a white noise of 0.01°/h/√Hz, accelerometers bias is 100 μg with a white noise of 10 μg/√Hz, misalignment of accelerometers is 20″, scale factor error of accelerometer is 100 ppm, the misalignment of gyro about rotational axis (z-axis of body frame) is 60″ and the others are 20″, and the scale factor error of gyro around rotation axis is 300 ppm and the others are 100 ppm. The main error parameters are listed in Table 1. Data sampling frequency is 50 Hz.

In order to validate the theory above, the initial attitude angles of roll, pitch, and yaw are 20°, 30°, and 40°, respectively. The tested IMU is under stationary base. The tested IMU rotates with the angular velocity of 10°/s, and the length of the test period t is 36 s.

4.2. Simulation Results

Assume that there is no attitude error at the beginning of navigation and the whole navigation process lasts 360 s, including 10 whole periods t. Based on the error sources listed in Tables 1 and 2, the simulation results are evaluated by 30 Monte Carlos runs. The state variables of instrument errors are , , and , and the values are 59.96″,59.96″, and 300 ppm, respectively. In Table 3, it is obvious that the estimated parameters obtained by 12-state filter are very close to the real values.

Figure 3 describes the error covariance root-mean-square curve of the state variables during the filtering process. The curve shows that the algorithm has better convergence speed, and the state variables , , and have good observability.

Figure 4 describes the estimated gyro scale factor errors and misalignment errors around rotational axis in 30 Monte Carlos runs. Similar to error covariance root-mean-square curve, error estimated curve also converges quickly after filtering. The steady values of estimated gyro scale factor errors and misalignment errors can be compensated during the navigation process. Navigation results after compensation are listed in Table 4.

Figure 5 describes the velocity error during navigation process without filter. The curve shows that velocity error accumulates with time quickly, which is the same as the unknown attitude error and the inertial device error related to rotational axis. Thus, the error mechanism deduction has been verified.

Table 4 provides a comparison of the navigation results under different conditions. Through 360 s navigation without compensation, the velocity errors Vx and Vy are 2.982 m/s and −17.548 m/s, and the attitude errors of roll, pitch, and yaw angles are 0.467°, 0.303°, and 0.911°. Meanwhile, Table 4 provides navigation results under the same condition after compensation. The velocity errors Vx and Vy are 0.065 and −0.024, and the attitude errors of roll pitch and yaw are −0.044°, 0.034°, and −0.092°, respectively.

Thus, simulation results show that velocity error caused by scale factor errors and misalignment errors related to rotation axis cannot be modulated by the single-axis rotational inertial navigation system. With the proposed method, the corresponding errors can be calibrated precisely. Simulation results show that navigation performance is greatly enhanced by calibration and compensation.

5. Real Test

5.1. Experiment Set up

To further validate the proposed method, the real calibration test is established. The test IMU consists of three orthogonal gyros whose accuracy is 0.1°/h and accelerometers whose accuracy is 100 μg. The test is carried out with personal computer made by DELL, industrial personal computer, the self-designed computer monitoring software, and navigation computing software.

Data are sampled by the frequency of 100 Hz. The turntable shown in Figure 6 is utilized to simulate single-axis rotational mechanization, whose rotation rate accuracy is 0.0001°/s. Besides, its positional accuracy is 2 arc-second. Initial angles of pitch, roll, and yaw are 0.0456°, 0.0279°, and 0.0197°. The rotation rate of turntable is 10°/s corresponding to the simulation.

In the beginning, the turntable remains stationary for 140 s to collect data for initial alignment. Then, the IMU rotates at the angular velocity of 10°/s around its z-axis at least 900 s for navigation. The turntable should be rotated to its original position after navigation.

5.2. Experiment Result

Gyro errors of rotation axis are estimated by the designed self-calibration filter. Estimated values of , , and are shown in Figure 3. It is obvious that estimation curve converges very quickly. The steady values of estimated scale factor errors and misalignment errors , , and related to rotational axis are utilized as calibration results. Calibration results estimated by filter are listed in Table 5, where is 154.2″, is −107.7″, and is 84.0 ppm.

Figure 6 shows the navigation results of accumulated velocity errors before compensation. It is evident that velocity accumulates quickly, and it can be proved that single-axis rotational INS cannot modulate scale factor error and misalignment errors related to rotation axis. Figure 7 displays the velocity estimated curve after compensating for the corresponding errors. Compared with the estimated curve before compensation, the error accumulates more slowly. It proves that the proposed method can decrease or eliminate the effects caused by gyro in rotation axis. Figure 8 shows the Vx, Vy curve after compensation.

Navigation results about velocity and attitude are listed in Table 6. Before compensation, the velocity errors in level drift are 4.9222 m/s and −1.9773 m/s, respectively, while after compensation are 0.2950 m/s and 0.5888 m/s. In terms of attitude errors, the improvement after compensation is also obvious. Compared with the initial attitude, the attitude errors accumulated during navigation and resolved by single-axis rotational INS in pitch, roll, and yaw are 0.0426 degree, 0.1037 degree, and 0.8712 degree. Similarly, attitude errors in pitch, roll, and yaw angle after compensation are 0.0011 degree, 0.0834 degree, and 0.0771 degree. Data express that the precision of navigation is enhanced at least 20%.

Thus, the proposed method with filter can estimate the device error related to rotational axis effectively. The precision of the single axis rotational INS is improved by compensating the estimated steady value of corresponding errors.

6. Conclusions

In the single-axis rotational inertial navigation system, instrument errors related to rotational axis cannot be modulated and these errors will have great impact on the navigation precision. Error mechanism and equations of velocity and attitude have been shown in this paper, and a calibration method with filter has been proposed. The state model should contain device errors which cannot be modulated by rotation. In the stationary base, zero velocity and the known turntable angular velocity can be chosen as the measurement. IMU errors can be estimated and compensated online to reduce the navigation error. Simulation has been carried out to verify the correctness of error mechanism and the feasibility of the proposed method, the experiment results further illustrate the validation of the proposed method, and the precision of the single-axis rotational INS can be improved.

It is shown that the precision of navigation is enhanced at least 20%. The proposed method using filter can estimate the device error related to rotational axis effectively. The precision of the single axis rotational INS has been improved by compensating the estimated steady value of corresponding errors.

Appendix

Deduction of equation (7):

In the equation, can be approximated aswhere , then inserting equation (5) into , we get

Deduction of equation (8): the velocity error accumulated with time can be written as the integral of equation (7):

Normally, accelerometer error can be modulated by rotation under stationary base. Thus, the simplified velocity model can be approximately written as

Combining with equation (7), equation (A.5) can be expressed as

Similarly, the error velocity can be expressed as

Data Availability

Due to project team constraints, other data materials cannot be provided for the time being.

Conflicts of Interest

The authors declare that they have no conflicts of interest.