A Three-Stage Accelerometer Self-Calibration Technique for Space-Stable Inertial Navigation Systems.

As a specific force sensor, the tri-axis accelerometer is one of the core instruments in an inertial navigation system (INS). During navigation, its measurement error directly induces constant or alternating navigation errors of the same order of magnitude. Moreover, it also affects the estimation accuracy of gyro drift coefficients during the initial alignment and calibration, which will indirectly result in navigation errors accumulating over time. Calibration can effectively improve measurement accuracy of the accelerometer. Device-level calibration can identify all of the parameters in the error model, and the system-level calibration can accurately estimate part of these parameters. Combining the advantages of both the methods and making full use of the precise angulation of the space-stabilized platform, this paper proposes a three-stage accelerometer self-calibration technique that can be implemented directly in the space-stable INS. The device-level calibration is divided into two steps considering the large amount of parameters. The first step is coarse calibration, which identifies parameters except for the nonlinear terms, and the second step is fine calibration, which not only identifies the nonlinear parameters, but also improves the accuracy of the parameters identified in the first step. The follow-on system-level calibration is carried out on part of the parameters using specific force error and attitude error to further improve the calibration accuracy. Simulation result shows that by using the proposed three-stage calibration technique in the space-stable INS, the estimation accuracy of accelerometer error can reach 1 × 10 - 6   g order of magnitude. Experiment results show that after the three-stage calibration, the accuracy of latitude, longitude, and attitude angles has increased by over 45% and the accuracy of velocity has increased by over 22% during navigation.

1. Introduction

Inertial navigation systems (INSs) tend to be the first choice for marine navigation because they are able to output motion parameters completely, autonomously, and continuously [1,2,3]. In order to fulfill the requirement for long-term and high-precision marine navigation, it is essential for INSs to employ high-precision inertial instruments. The gyroscope should be able to work stably for a long time, of which the drift error needs to be small and the constant drift can be precisely compensated. And as for the accelerometer, high stability of the bias and the scale factor is required [4].

All of the inertial instruments must be precisely tested and calibrated before they can work properly to provide accurate measurement in the INS [5,6]. According to the observed quantities, accelerometer calibration can be divided into two types, device-level calibration and system-level calibration. In device-level calibration, the accelerometer is set to measure the gravity of Earth at different positions by using a dividing head, a centrifuge, or a turntable, and error parameters are identified according to the outputs of the accelerometer [7,8,9]. The accelerometer data measured during the multi-position rotation of the turntable are compared with the theoretical value of the local gravity, and then using the least square method, error parameters can be estimated from their difference. The system-level calibration is based on the propagation of error parameters of inertial instruments to the final navigation errors. Error parameters can be precisely identified with a proper rotation scheme that is designed considering the parameter observability [3].

Jiang et al. [9] proposed a device-level calibration method based on the static-base condition, and the precise attitude and attitude rate provided by a tri-axis precision turntable, where the calibration accuracy can reach order of magnitude. However, in practical applications, it is difficult to provide an ideal laboratory environment for the implementation of this method. Won et al. [10] and Ye et al. [11] calibrated the accelerometer on the basis that the amplitude of the output vector of the tri-axis accelerometer equals 1 g. Studies [12,13,14,15,16,17,18,19,20,21] presented several system-level calibration methods specific to the microelectromechanical system (MEMS) inertial measuring unit (IMU), but the calibration accuracy needs to be improved. The methods proposed in references [22,23,24,25,26,27,28,29,30] are appropriate for the rotational INSs, and several novel methods are presented for the hybrid INSs in references [31,32]. Gao et al. [22] used position error and velocity error as the observed quantities to estimate the nonlinear error coefficients. They also used velocity error and attitude error to estimate the error coefficients of the gyroscope and the accelerometer [23]. Pan et al. [33] calibrated the nonlinear term in the scale factor based on the change rate of velocity error. Gao et al. [27] and Liu et al. [31] proposed methods to compensate the size effect and the lever arm effect, respectively. As for the turntable rotation strategy, the 6-position [22], 9-position [24,32], 18-position [18,34], and 24-position [35] methods have been proposed. In existing works, calibration methods are rarely designed for the space-stable INSs. In addition, the complementary benefits of device-level and the system-level calibration have not been taken into consideration.

In the space-stable INS, the platform is stabilized in inertial space by the gimbal servo system. With no applied torque, the moment of momentum of the rotor in the free-rotor gyroscope maintains fixed orientation in inertial space, which can ensure the high precision of the system [1]. As there is a stabilized platform in a space-stable INS, self-calibration based on such an own condition becomes the current development trend. This paper attempts to combine the advantages of full-parameter estimation in the device-level calibration and high-precision estimation of partial parameters in the system-level calibration. A three-stage self-calibration technique is proposed for the high-precision accelerometer calibration in practical applications. Firstly, device-level self-calibration is conducted directly in the system. It is divided into two steps considering the large amount of parameters. The first step is coarse calibration which identifies parameters except for the nonlinear terms, and the second step is fine calibration which identifies the nonlinear terms and simultaneously improves the accuracy of parameters identified in the first step. Secondly, the follow-on system-level calibration is carried out on part of the parameters using the errors of specific force and attitude angles to further improve the calibration accuracy. The proposed method is specific to the space-stable INS and has three contributions. First, taking advantage of the high-precision angulation ability of the gimbaled platform, the two-stage device-level calibration can be directly implemented in the system. Second, based on the error propagation of the residual accelerometer model error to navigation errors, the system-level calibration uses specific force error and attitude error to form the observation model. Third, the proposed three-stage calibration method merges the complementary advantages of both the device-level calibration and the system-level calibration to achieve full-parameter calibration of the accelerometer measurement model as well as high-precision estimation of part of these parameters. This can help further improve the performance of the system.

The rest of this paper is arranged as follows: In Section 2, the accelerometer measurement model is introduced. The device-level coarse calibration and device-level fine calibration are presented in Section 3 and Section 4, respectively. The system-level calibration method based on the errors of specific force and attitude angles is proposed in Section 5. Simulation and experiment results are discussed in Section 6. And finally conclusions are drawn in Section 7.

2. Accelerometer Measurement Model

For the convenience of following discussion, several coordinate frames are defined first, and the superscript of a vector will indicate which frame it is in.

—the local geographic navigation frame (n-frame), where the axis and the axis are in the local-level plane pointing north and east, respectively, and the axis is angled vertically downward. The gravity vector has the same orientation with the axis in this coordinate system.

—the body-fixed coordinate frame (b-frame). The origin is usually located at the center of the mass of the vehicle. The axis points forward along the longitudinal axis, the axis points to the bottom, and the axis points rightward to form a right-handed coordinate system.

—the platform-fixed coordinate system (p-frame). The axis is along the axis of revolution of the platform. And ideally, and are along the revolution axes of the inner gimbal and the middle gimbal, respectively.

—the accelerometer coordinate system (a-frame), which is fixed on the tri-axis accelerometer. The , and axes are along the three sensitive axes, respectively. This coordinate system is non-orthogonal due to the installation error.

For a static system, the specific force measured in the n-frame should be: where g is the scalar gravity. Transforming it to the p-frame yields: where the direction cosine matrix from the n-frame to the b-frame is: and the direction cosine matrix from the b-frame to the p-frame is:

In Equation (3), are yaw, pitch, and roll, respectively. In Equation (4), is the rotation angle of the platform’s revolution axis relative to the revolution axis of the inner gimbal, h is the rotation angle of the revolution axis of the middle gimbal relative to that of the outer gimbal, and q is the rotation angle of the revolution axis of the outer gimbal relative to the base. The rotation angle of the revolution axis of the inner gimbal relative to that of the middle gimbal is zero, for the outer gimbal is tracking the rotation angle of the inner gimbal, and for the revolution axes of the platform, the inner gimbal and the middle gimbal are orthogonal.

The geometric relationship between the orthogonal p-frame and the non-orthogonal a-frame is shown in Figure 1. It is assumed that installation error angles are all small quantities. denote the three installation error angles in the x-direction, y-direction, and z-direction, respectively. Each of these three small angles can be decomposed into two rotation angles around the axes of the p-frame, namely , and .

Figure 1

The geometric relationship between the p-frame and the a-frame.

According to the geometric relationship shown in Figure 1, the direction cosine matrix from the p-frame to the a-frame can be calculated as:

Thus the specific force measured by the tri-axis accelerometer satisfies an equation, as: where is the matrix of scale factors, is the bias vector, and is the matrix of nonlinear coefficients of the scale factors. In this paper, a three-stage calibration technique is proposed, to estimate and compensate all the parameters shown in Equation (6). The first two stages will lay a foundation of full-parameter and full-value estimations of model coefficients. The third stage will then correct part of the coefficients and further improve the calibration accuracy. Although the final calibration accuracy will mainly depend on the novel system-level calibration method, it must work on the basis of the device-level calibration. Therefore, the first two stages of the device-level calibration designed for the space-stable INS will also be introduced to form the integrated three-stage calibration technique.

3. Device-Level Coarse Calibration

For the tri-axis accelerometer, there are totally 18 error coefficients that need to be calibrated: three scale factors, ; three biases, ; three nonlinear coefficients of the scale factors, ; six installation error angles , ; and three zero-point offsets of the gimbals’ angular sensors, of the platform axis, of the inner gimbal axis, and of the middle gimbal axis. The platform needs to be rotated to 18 positions. Table 1 lists the order of rotation, and the corresponding rotation angles of the gimbals and components of the p-frame gravity vector [3], where the symbol denotes that the angle remains unchanged. Figure 2 shows the directions of the accelerometer’s sensitive axes in detail, where N, W, and U denote the local northern, western and up directions, respectively, and x, y, and z denote the orientations of the tri-axis accelerometer.

Table 1

The order of rotation for the 18-position calibration, and the corresponding rotation angles of the gimbals and components of the p-frame gravity vector.

The Order of Rotation	Rotation Angles of Gimbals			Components of the p-Frame Gravity Vector
	θ_1	θ_2	θ_3	gxp	gyp	gzp
1	θ_10	θ_20	θ_30	−gθ_30	gθ_20	g
2	↓	θ_20+π4	↓	(22θ_10−θ_30)g	22g	22g
3	↓	θ_20+π2	↓	g(θ_10−θ_30)	g	−gθ_20
4	θ_10+π4	↓	↓	22g	22g	−gθ_20
5	θ_10+π2	↓	↓	g	g(θ_30−θ_10)	−gθ_20
6	↓	θ_20+π4	↓	22g	(θ_30−22θ_10)g	22g
7	↓	θ_20+π	↓	−gθ_20	gθ_30	−g
8	↓	↓	θ_30−π4	−22g(θ_10+θ_20)	−22g	−22g
9	↓	↓	θ_30−π2	−gθ_10	−g	−gθ_30
10	θ_10+3π4	↓	↓	−22g	−22g	−gθ_30
11	θ_10+π	↓	↓	−g	gθ_10	−gθ_30
12	↓	↓	θ_30−π4	−22g	22g(θ_10+θ_20)	−22g
13	↓	θ_20+5π4	θ_30	(22θ_10+θ_30)g	22g	−22g
14	θ_10+3π4	θ_20+3π2	↓	−22g	22g	22gθ_20
15	θ_10+π2	θ_20+7π4	↓	−22g	22g(θ_30+θ_20)	22g
16	↓	θ_20	θ_30−π4	22g(−θ_10+θ_20)	−22g	22g
17	θ_10+π4	↓	θ_30−π2	22g	−22g	gθ_30
18	θ_10	↓	θ_30−3π4	22g	−22g(θ_10+θ_20)	−22g

Figure 2

The 18 positions chosen for accelerometer calibration.

During coarse calibration, the nonlinearity of scale factors is ignored. Thus the accelerometer measurement equation shown in Equation (6) can be simplified as: where is the gravity vector sensed by the tri-axis accelerometer.

Each position takes 20 s for the rapid measurement. Based on the outputs of the tri-axis accelerometer at different positions and Table 1, 18 equations can be obtained according to Equation (7) to preliminarily estimate 15 parameters except the nonlinear coefficients. Scale factors and biases can be coarsely estimated by solving the equation set corresponding to the six rotation positions (No. 1, 3, 5, 7, 9, and 11) where each sensitive axis of the tri-axis accelerometer points upward or downward along the axis in turn. The other equations corresponding to the rest of the rotation positions (No. 2, 4, 6, 8, 10, and 12~18), where one of the sensitive axes is parallel to the axis or the axis, and other two sensitive axes have included angles of 45° with the vertical direction, are used to estimate zero-point offsets of the angular sensors and installation error angles.

Estimated values of zero-point offsets, and , are used to coarsely correct the orientation of the p-frame, and to make the directions of the middle gimbal axis and the inner gimbal axis approach the directions of axes of the p-frame. Concretely, the inner gimbal axis and the middle gimbal axis are inversely rotated and , respectively, and then coarse calibration is re-executed. This procedure will repeat until the re-estimated and decrease to negligibly small quantities.

4. Device-Level Fine Calibration

Fine calibration will be executed on the basis of coarse calibration. The orientation of the platform has been corrected and new data will be measured in more accurate rotation positions. In order to improve the data collection accuracy, the measurement time of each position is extended to 10 min. Moreover, to improve the data processing accuracy, the estimation error of coarse calibration is estimated instead of the full values of error parameters.

According to Equation (6), the ideal specific force in the p-frame satisfies:

If the result of coarse calibration and outputs of the tri-axis accelerometer during fine calibration () have been obtained, the coarse estimation of can be calculated as: where “hatted” symbols denote the estimated values.

In order to reduce the requirement for position control accuracy of the platform during fine calibration, the observed variable is not the direct difference between the coarsely estimated value and the ideal value of the specific force at one certain position, but the difference between the squared amplitudes of them. Subtracting the squared amplitude of Equation (9) from that of Equation (8), and then taking half of the result yields the observed variable, as: where:

According to Equations (10)–(15), the observation equation at position j can be written as: where the unknown error vector is:

Based on the data collected at the 18 positions, an -dimension observation equation can be formed based on Equation (16), and the unknown variables shown in Equation (17) can be solved by the least squares method. Finally, using scale factors to convert the solved and adding the result to the coarse estimation will yield a full-value fine estimation of 12 error coefficients. It should be noted that , are linearly correlated in pairs and cannot be decomposed from each other. Therefore the fine calibration can only improve the accuracy of , but it cannot calibrate the rotation characteristics of the a-frame.

5. System-Level Calibration Based on Specific Force Error and Attitude Error

Errors of inertial instruments will induce the navigational errors of INSs, including the positional error, the velocity error and the attitude error. Inversely, the independently observed navigation errors such as the specific force error, the attitude error and the position error can be used to calibrate the main error parameters in the system. As the position error is affected by both the gyroscope error and accelerometer error, it is generally difficult to directly extract the accelerometer-induced term from the position error. On the other hand, the specific force error and the horizontal attitude error are mainly caused by the accelerometer error. Therefore this paper employs the specific force and horizontal attitude angles to correct part of the accelerometer’s error coefficients. The prerequisite for accurate observation in such a system-level calibration is that the system should be installed on a static and horizontal base.

5.1. Calibration Using Specific Force Error

Under the ideal condition, that there is no installation error between the base and the stabilized platform, and that the measurement of gimbals’ rotation angle is accurate, there is an approximate relationship: where L denotes the latitude. On the static and horizontal base, the local scalar gravity g can be accurately calculated. Then, according to Equations (1) and (18), the ideal specific force vector in the p-frame can be written as:

Similar to the device-level calibration method presented in the previous section, the difference of the squared amplitudes of the measured and the ideal specific force vectors is chosen as the observed variable, as:

A 5-dimensional state vector is constructed as , where:

Thus the observation equation is:

Considering Equation (19), the observation matrix is:

The five components of are actually the mean value, the sine and the cosine fundamental waves, and the sine and the cosine second harmonics of the Earth’s rotation rate in the specific force error. Assuming that the observation noise is a zero-mean Gaussian white noise, can be estimated by a Kalman filter.

5.2. Calibration Using Attitude Error

On the static and horizontal base, the attitude error of the INS with both vertical damping and horizontal velocity damping mainly consists of a constant value and a 24 h long period component, which is induced by the gyroscope error, the accelerometer error, and the rotation angle errors of platform’s gimbals. Using the space-stable mechanization, the azimuth error will accumulate over time due to the gyroscope error, while the horizontal attitude errors are mainly caused by the accelerometer error and have no such accumulation. Therefore, the horizontal attitude errors are chosen as the observed quantities for accelerometer calibration.

Neglecting the short-period misalignment angle of gyro case rotation, the attitude error on the static and horizontal base has an expression as: where , , and are the transverse tilt error, the longitudinal tilt error, and the azimuth error of the base, respectively, and the superscript e denotes the earth-fixed frame (e-frame).

Substituting the e-frame specific force error into Equation (28) yields the expressions of pitch error and roll error. Both of them consist of the mean value, the fundamental wave, and the second harmonic of the Earth’s rotation rate. The five Fourier coefficients of the pitch error are:

The other five Fourier coefficients of the roll error are:

The subscripts indicate which component each expression is (1 for the fundamental waves, 2 for the second harmonics, s for the sine functions, and c for the cosine functions). After the error curves of pitch and roll are obtained, these coefficients can be extracted by Fourier expansion.

5.3. Parameter Separation Algorithm for Combined Error Coefficients of Accerelometer

In Section 5.1 and Section 5.2, 15 combined coefficients (Equations (21)–(25), (29)–(33), and (34)–(38)) are estimated by static Kalman filter and Fourier expansion. Twelve of them are selected and divided into three groups to extract the residual error parameters of the tri-axis accelerometer. Since only the fundamental waves and second harmonics of pitch and roll errors are concerned, the requirement for the horizontality of the base need not be very strict.

The first group consists of three equations, which are the coefficients of the sine and cosine second harmonics of Earth’s rotation rate about . They are the linearly correlated Equations (24), (32), and (37). Thus can be solved using their equally weighted average, as:

Similarly, the second group consists of Equations (25), (33), and (38), and can be solved as:

Another six parameters (, , , , , ) can be solved using the least square method based on the remaining six combined coefficients, according to:

The six parameters in Equation (41) can be directly compensated by substituting them into the accelerometer’s error model. Also, and will be indirectly compensated by deducting the error terms about them from the n-frame specific force expression.

6. Simulation and Experiment

6.1. Simulation Result

As the final accuracy of the three-stage calibration mainly depends on the performance of the last stage, the simulation of the system-level calibration is conducted first. In the simulation, the location of the static base is set at 40° N. The non-orthogonal installation error of the tri-axis accelerometer is quite a stable mechanical error, and after the 18-position two-stage device-level calibration, the residuals are generally less than 10″. Thus, the residuals of accelerometer’s error parameters are assigned as Table 2 and Table 3. In addition, a zero-mean random noise with an amplitude of is added to each sensitive axis of the tri-axis accelerometer.

Table 2

Residuals of scale factor errors and biases.

ΔSF_x/SF_x	ΔSF_y/SF_y	ΔSF_z/SF_z	∇_x	∇_y	∇_z
10^−4	2×10^−4	−10^−4	−10^−3 m/s^2	5×10^−4 m/s^2	10^−3 m/s^2

Table 3

Residuals of installation error angles.

Δβ_xz	Δβ_yz	Δβ_zy	Δβ_xy	Δβ_yx	Δβ_zx
5″	10″	−5″	−10″	5″	−5″

Truth values of the system-level combined coefficients are calculated according to the error parameters shown in Table 2 and Table 3, and are listed in Table 4 as reference values. The combined coefficients are estimated by system-level calibration based on specific force error and attitude error, and are also listed in Table 4. Comparison shows that the estimation accuracy of each combined coefficient is better than , which is equivalent to an accelerometer calibration accuracy that is better than g. It indicates that the method can be used for the accelerometer calibration of high-precision INSs.

Table 4

Simulation result of system-level calibration of the combined accelerometer error coefficients (unit: rad).

Combined Coefficients	Reference Values	Estimated Values	Errors
x_1	1.1228×10^−4	1.1256×10^−4	2.8×10^−7
x_2	−8.645×10^−5	−8.612×10^−5	3.3×10^−7
x_3	8.217×10^−5	8.206×10^−5	−1.1×10^−7
x_4	−5.000×10^−5	−4.957×10^−5	4.3×10^−7
x_5	1.212×10^−5	1.185×10^−5	−2.7×10^−7
δϕ_1s	−1.3319×10^−4	−1.3343×10^−4	−2.4×10^−7
δϕ_1c	−3.542×10^−5	−3.575×10^−5	−3.3×10^−7
δϕ_2s	3.830×10^−5	3.848×10^−5	1.8×10^−7
δϕ_2c	−9.28×10^−6	−9.47×10^−6	−1.9×10^−7
δθ_1s	−5.521×10^−5	−5.536×10^−5	−1.5×10^−7
δθ_1c	1.3629×10^−4	1.3695×10^−4	6.6×10^−7
δθ_2s	−5.97×10^−6	−5.76×10^−6	2.1×10^−7
δθ_2c	−2.462×10^−5	−2.484×10^−5	−2.2×10^−7

6.2. Experiment Results

In order to validate the proposed three-stage accelerometer calibration method, a static-base navigation test has been conducted using a high-precision space-stable INS. The INS is in the navigation grade class and has been introduced in reference [36]. Two dual-axis free-rotor gyros and a tri-axis quartz accelerometer are mounted on a gimbaled platform. According to the outputs of the free-rotor gyros, the platform is servo-controlled to be stabilized in the inertial space. The gyros have an ultra-low drift of less than 0.001°/h. The measurement error of the accelerometer is of a g order of magnitude, and its stability is better than g per year. The static measurement accuracy of gimbals’ rotation angles can reach the level of 1″. The position error of this navigation system can reach the order of magnitude of 1 n mile per day. The data used in the experiment are collected when the system is installed on a static and horizontal base.

Navigation errors before calibration, after the original two-stage 18-position coarse and fine self-calibration, and after the entire three-stage calibration are compared in Figure 3, Figure 4 and Figure 5 and Table 5. Since the system employs vertical damping, height error and vertical velocity error can be neglected and are not presented. The output position of the system is compared with the reference position of the base to obtain the position error. The output velocity and horizontal attitude angles are actually the velocity error and horizontal attitude error, for the base is static and horizontal. All experiment results here have been normalized.

Figure 3

Horizontal position errors. (a) Longitude errors; (b) latitude errors.

Figure 4

Horizontal velocity errors. (a) Northern velocity errors; (b) eastern velocity errors.

Figure 5

Attitude errors. (a) Roll errors; (b) pitch errors; (c) yaw errors.

Table 5

Maximums of navigation errors.

Navigation Errors	Before Calibration	Two-Stage Calibration	Improvement afterTwo-Stage Calibration	Three-Stage Calibration	Improvement afterThree-Stage Calibration
Longitude	0.84	0.58	30.95%	0.33	60.71%
Latitude	0.86	0.60	30.23%	0.26	69.77%
Northern Velocity	0.85	0.70	17.65%	0.56	34.12%
Eastern Velocity	0.87	0.78	10.34%	0.67	22.99%
Roll	0.89	0.50	43.82%	0.38	57.30%
Pitch	0.95	0.58	38.95%	0.43	54.74%
Yaw	0.80	0.54	32.50%	0.44	45.00%

The system has employed gyro case rotation techniques to decrease the gyro drift. Although its influence on the platform’s motion has been deducted during navigation calculating, there are still 8-minute residuals in the velocity and attitude error curves. As the cycle of these residuals are much shorter than the test time, they are severely squeezed in the error curves shown in Figure 4 and Figure 5.

Accelerometer parameters will change after the high-precision INS works for a period of time, which will consequently lead to the increase of navigation errors. Therefore, on-site and online calibration is necessary to improve the navigation precision. Latitude error, longitude error, and attitude errors are reduced by over 30% after the two-stage calibration, and by over 45% after the three-stage calibration. Since the velocity error mainly consists of Schuler oscillation, and the error term is related to the gyro case rotation residual, the improvement of velocity accuracy is not conspicuous. The velocity error has decreased by over 10% after the two-stage calibration, and by over 22% after the three-stage calibration.

7. Conclusions

A three-stage tri-axis accelerometer self-calibration technique is proposed with a thorough theoretical and mathematical analysis. Compared with existing calibration methods, the proposed three-stage calibration method has the following features and advantages: (1) Taking full advantage of high-precision angulation of the gimbaled platform in the space-stable INS, the two-stage device-level calibration is implemented directly on the platform of the space-stable INS. Thus the entire three-stage self-calibration can be implemented in the INS after the tri-axis accelerometer has been mounted on the platform. Besides, during the device-level calibration, zero-point offsets of the gimbals’ angular sensors are simultaneously estimated and used to correct the platform’s orientation. This will help improve the system’s navigation performance. (2) Different from the existing calibration methods which usually estimate error parameters of both accelerometer and gyroscope simultaneously, the proposed system-level calibration is designed based on error propagation to choose observed quantities only resulting from accelerometer measurement error. This can get rid of the influence of gyro drifts. (3) Combining device-level calibration and system-level calibration, the proposed method can achieve both full-parameter estimation of the accelerometer measurement model and high-precision correction for part of the parameters, while existing calibration methods normally concentrate on one of the two aspects. Both simulation and experiment results validate the proposed method. The simulation result shows that using the proposed three-stage calibration technique, the estimation accuracy of accelerometer error is better than g. Experiment results show that after implementing the three-stage calibration, during navigation, the accuracy of latitude, longitude, and attitude angles has increased by over 45%, and the accuracy of velocity has increased by over 22%.