Analysis of Gyro Bias Depending on the Position of Inertial Measurement Unit in Rotational Inertial Navigation Systems

In this paper, a calibration method for gyro bias that changes depending on the position of the IMU (inertial measurement unit) is proposed to improve the navigation performance of RLG-based RINS (ring-laser-gyro-based rotational inertial navigation system). RINS is a navigation device that compensates for the inertial sensor errors by utilizing the rotation of the IMU. In previous studies, the rotation scheme of the IMU is designed assuming that inertial sensor errors are not affected by position of the IMU. However, changes in temperature distribution, direction of gravity, and dithering according to the rotation of the IMU affect the inertial sensor errors, such as gyro bias. These errors could degrade the long-term navigation performance of RLG-based RINS. To deal with this problem, this paper proposed a compensation method of the gyro bias that changes depending on the position of the IMU. First, RINS is reviewed using a dual-axis 16-position rotation scheme and RLG. Next, the attitude error of RLG-based RINS is derived utilizing navigation equations. The effect of the gyro bias change caused by the change in the IMU attitude for the navigation performance of RINS is analyzed based on navigation equations and simulations. Finally, system-level indirect calibrations for the Z–axis up position and Z–axis down position are performed to calculate the gyro bias change caused by the IMU attitude. The accuracy of the proposed calibration method is verified by long-term navigation test. The test results show that the proposed calibration method improves the navigation performance of RINS compared with the conventional calibration method.


Introduction
INS (inertial navigation system) is a device that calculates navigation information such as position, velocity, and attitude using measurements of the IMU (inertial measurement unit) [1][2][3][4]. INS is widely used in various fields such as weapons systems, robots, and aerospace because it can calculate accurate navigation information without the help of external sensors. However, the navigation information calculated by INS has a problem that sensor errors such as bias, scale factor, and misalignment are accumulated over time. These sensor errors could degrade the long-term navigation performance of the INS.
To deal with this problem, a rotational inertial navigation system (RINS) is developed. RINS, a type of INS, operates in a way that the IMU is continuously rotated according to a given rotation scheme [5][6][7][8]. As RINS can suppress sensor errors through rotation of the IMU, it has the advantage of being able to improve navigation performance without the need for constructing IMU with a high-accuracy gyro or accelerometer. Because of these depending on the rotation position of the IMU [31]. As the sensor error of RLG is affected by these characteristics, the sensor error of RLG also changes depending on the position of the IMU during the operation of RINS [32]. Therefore, in RLG-based RINS, there is a problem that the sensor error compensation accuracy of the rotation scheme suggested in previous research may be degraded [33][34][35].
To deal with this problem, this paper proposes a calibration method to improve the navigation performance of RLG-based RINS by pre-identifying and compensating for sensor errors in RLGs that change with the position of the IMU. The sensor error of the RLG affected by the position of the IMU was defined as a gyro bias according to the location and the system stage indirect calibration was introduced and identified.
In this paper, RINS is constructed using a three-axis RLG and a three-axis silicon accelerometer, and the rotation scheme of IMU is designed by importing a 16-position rotation scheme based on the dual-axis rotation table presented in [29]. The effect of gyro bias according to the position of the IMU on the navigation performance of RINS was analyzed through simulation. Finally, system-level indirect calibration was performed for various positions of the IMU to identify the gyro bias according to the position of the IMU. Navigation was performed using IMU set in the Z-axis up position and the Z-axis down position, and the gyro bias according to the position of the IMU was identified using the navigation information obtained during the navigation process. The gyro bias according to the location obtained through indirect calibration at the system stage was designed to be compensated in real time in the operation process of RINS. The navigation performance of the proposed calibration method was verified through long-term navigation test.
This paper is organized as follows. Section 2 describes the error model of RINS and the dual-axis 16 position rotation scheme for the IMU. Section 3 presents a procedure of identifying gyro bias that changes depending on the position of the IMU through systemlevel indirect calibration and conventional calibration methods. Section 4 verifies the performance of the proposed calibration method via long-term static state test. Finally, Section 5 presents the conclusion of this study.

Rotational Inertial Navigation System (RINS)
Coordinate frames applied in this paper are defined as shown in Table 1.

Symbol Frames b
The body frame e The earth frame n The navigation frame i The orthogonal inertial frame The body frame is applied to describe the position, velocity, and orientation of the sensor platform. The body frame consists of an origin that is typically placed at the center of gravity and three orthogonal axes. Generally, the axes in the body frame are configured to the body such that the X-axis is pointing forward, the Y-axis is pointing to the right, and the Z-axis is pointing up.
The earth frame is a global reference frame with its origin at the center of the Earth and three orthogonal axes fixed to the Earth. In the earth frame, the x, y, and z axes are aligned to local east, north, and up axes. We should point out that the earth frame rotates with the Earth.
The navigation frame is a frame that directs east-north-up (ENU). The navigation frame is fixed to the platform and moves with the body frame. The navigation frame contains three orthogonal axes in which the axis points to true east, the axis points towards the interior of the Earth, and the axis completes the right-handed system pointing north.
The orthogonal inertial frame has its origin at the center of the Earth. In contrast to the earth frame, the orthogonal inertial frame does not rotate with Earth.

Systematic Error Model of RINS
This section describes systematic errors that occur during the operation of RINS. The systematic error of RINS may be analyzed using the position and acceleration error model. The position error model of RINS may be expressed as follows using the error model of INS presented in [36].
where φ and . φ mean the attitude error and rates of the attitude error, respectively. ω n in is the earth rate and δω n in is the earth rate error. n is the gyro error of the navigation frame. The acceleration error model of RINS may be expressed as follows in the same manner as the position error model.
where δv and δ . v are velocity error and acceleration error for RINS, respectively. f means specific force and δg means gravitational error. ∇ n is accelerometer error for the navigation frame.
The gyro error in Equation (1) and accelerometer error in Equation (2) consist of static error such as white noise and random walk and deterministic error such as bias, scale factor, and misalignment [37]. Stochastic error cannot be modeled and has little effect on gyro error. The deterministic error occurs according to a certain tendency and has a greater effect on navigation performance than the static error. Therefore, the gyro error and the accelerometer error of the RINS are generally expressed using a deterministic error.
Let ω b = ω x , ω y , ω z T be the input angular rate about body frame, the gyro error for navigation frame is defined as follows: where C n b is the coordinate transformation matrix for the body frame and δω b ib is the gyro error for the body frame. β x , β y , β z are the gyro bias; β xx , β yy , β zz are the gyro scale factor error; and β xy , β xz , β yx , β yz , β zx , β zy are the gyro misalignment.
Let f b = f x , f y , f z T be the specific force about the body frame, the accelerometer error for the navigation frame is defined using the deterministic error of the accelerometer: where δ f b is the accelerometer error for the body frame. α x , α y , α z are the accelerometer bias; α xx , α yy , α zz are the accelerometer scale factor error; and α xy , α xz , α yx , α yz , α zx , α zy are the accelerometer misalignment. The gyro error and accelerometer error defined in Equations (3) and (4) are compensated through the rotation of the IMU during the operation of RINS. The navigation performance of RINS can be effectively improved by rotating the IMU according to an elaborately designed rotation scheme.
In this paper, the dual-axis 16-position rotation scheme presented in [29] is applied to rotate the IMU of RINS. The IMU used in RINS was designed using a three-axis RLG and a three-axis silicon accelerometer.

Navigation Performance of RINS Based on the Dual-Axis 16-Position Rotation Scheme
This section analyzes the navigation performance of RINS designed in Section 2.1 based on system error. As system error is mainly affected by gyro, the navigation perfor- mance of RINS can be determined using position error. The dual-axis 16-position rotation scheme imported to rotate the IMU of RINS in this paper consists of rotation sequences 1 to 4, rotation sequences 5 to 8, rotation sequences 9 to 12, and rotation sequences 13 to 16. The dual-axis 16-position rotation scheme is designed to perform four rotations at position A and return to the initial position A. In the process of performing the rotation scheme, the IMU rotates at a constant speed and stops at positions A, B, C, and D for a certain period of time. The rotation order of the dual-axis 16-position rotation scheme is shown in Figure 1.
to rotate the IMU of RINS. The IMU used in RINS was designed using a three-axis RLG and a three-axis silicon accelerometer.

Navigation Performance of RINS Based on the Dual-Axis 16-Position Rotation Scheme
This section analyzes the navigation performance of RINS designed in Section 2.1 based on system error. As system error is mainly affected by gyro, the navigation performance of RINS can be determined using position error. The dual-axis 16-position rotation scheme imported to rotate the IMU of RINS in this paper consists of rotation sequences 1 to 4, rotation sequences 5 to 8, rotation sequences 9 to 12, and rotation sequences 13 to 16. The dual-axis 16-position rotation scheme is designed to perform four rotations at position A and return to the initial position A. In the process of performing the rotation scheme, the IMU rotates at a constant speed and stops at positions A, B, C, and D for a certain period of time. The rotation order of the dual-axis 16-position rotation scheme is shown in Figure 1.
where means the angular rates of the IMU in the process of conducting the rotation sequence. In the dual-axis 16-position rotation scheme, the coordinate transformation matrices for position A, B, C, and D C n b t (t = A, B, C, D) are shown in Equation (5).
During the process of the rotation sequence, the angular rate of the IMU ω b t (t = 1, . . . , 16) for the body frame could be organized into Equation (6).
where ω r means the angular rates of the IMU in the process of conducting the rotation sequence.
For the dual-axis 16-position rotation scheme, the coordinate transformation matrices for rotation sequence C n b t (t = 1, . . . , 16) are shown in Table 2. The position error in the operation phase of RINS with dual-axis 16-position rotation can be calculated by integrating the rate of change in position error in Equation (1). As the time to perform dual-axis 16-position rotation is relatively short, this paper assumes that the effect of the Earth's rotational angular velocity on the position error of RINS is negligibly small.
The position error occurring in the operation process of RINS is calculated as follows using the assumptions and Equations (1), (3), (5), and (6) and Table 2.
In Equation (7), T r is the rotation time of the IMU and T s is stationary time of the IMU.
In the process of conducting dual-axis 16-position rotation, the attitude error of RINS could be expressed as follows: As can be seen from Equation (8), gyro sensor errors such as bias, scale factor, and misalignment are compensated during the operation of RINS designed by applying dualaxis 16-position rotation. In the same way, accelerometer errors including bias, scale factor, and misalignment are compensated. Therefore, RINS in this paper can theoretically improve position errors on the east, north, and up axes.

Theorical Analysis of the Gyro Bias Change Caused by Position of the IMU
To improve the navigation performance of RINS by applying dual-axis 16-position rotation, the gyro sensor error should not be affected by the position of the IMU. However, the gyro error is affected by the gyro internal temperature distribution, gravitational direction, and dithering, which change with the position of the IMU. The gyro error affected by the position of the IMU can degrade the compensation accuracy of the gyro error such as bias, scale factor error, and non-alignment error, as well as the navigation performance of RINS. Therefore, in order to precisely compensate for the gyro error by performing dual-axis 16-position rotation, it is necessary to identify and compensate for the gyro error that varies with the location of the IMU.
In this section, the effect of the gyro error according to position of the IMU described in Section 2.2 is analyzed. This paper defined the gyro bias that changes according to the position of the IMU as δβ = δβ x , δβ y , δβ z T . The gyro bias that changes according to the position of the IMU is affected by the position of the IMU and the input angular velocity is not affected. The gyro error of Equation (3) may be recalculated as follows by importing a gyro bias that changes according to the position of the IMU.
The effect of the gyro bias, which changes according to the position of the IMU, on the position error of RINS is recalculated as follows using Equations (1), (5), and (6) and Table 2.
According to Equation (10), the position error occurring during the operation of RINS is calculated by adding the integration of the gyro bias that changes according to the position of the IMU and the integration of the remaining gyro sensor error. As dual-axis 16position rotation applied in this paper can compensate for the remaining gyro sensor error, the navigation performance of RINS can be improved by identifying and compensating for the gyro bias that changes according to the position of the IMU.
The position error occurring in the operation process of RINS is summarized as follows by substituting Equations (5) and (6) and Table 2 into Equation (10).
As shown in Equation (11), the gyro bias that changes according to the position of the IMU during the RINS operation phase is compensated in the x-axis and y-axis. However, the z-axis gyro bias δβ z , which varies with the position of the IMU, causes an up-axis position error proportional to T s and δβ z . In the case of gyro errors including scale factor and misalignment and accelerometer errors including bias, scale factor, and misalignment, these errors are compensated for through the dual-axis 16-position rotation scheme.
For example, if T s is 180 s and δβ z is 0.005 deg/h, the up-axis attitude error during the operation of RINS is approximately 0.002 deg. Because the up-axis attitude error leads to the latitude error and longitude error, it is necessary to compensate for the gyro bias that changes according to the position of the IMU to improve navigation performance.
In this paper, numerical simulation is performed to evaluate the effect of the gyro bias that changes according to the position of the IMU. The changes in navigation performance were compared by performing navigation simulations on conditions that compensate δβ x , δβ y , and δβ z and non-compensating conditions. The position where the navigation simulation is performed is 36.6 • latitude and 123.6 • longitude and the navigation time is 1 week. The angular velocity ω s of the IMU rotating according to the dual-axis 16-position rotation is 2 deg/s and the time T s of the IMU stopping at positions A, B, C, and D is 180 s. The sensor specification of RINS used in the simulation is shown in Table 3. The position error obtained by performing a navigation simulation according to the conditions in Table 3 is shown in Figure 2. The position error obtained by performing a navigation simulation according to the conditions in Table 3 is shown in Figure 2. The navigation simulation results show that the up-axis average attitude error when not compensating for the gyro bias that changes according to position of the IMU is 0.0338 and the up-axis average attitude error when compensating for the gyro bias that changes according to position of the IMU is 0.0302 . Consequently, the up-axis attitude error can be improved by 10.65% by compensating for the gyro bias that changes according to the position of the IMU.
The latitude error and longitude error obtained by performing the navigation simulation are shown in Figure 3 and the horizontal axis position error at the end of the navigation is summarized in Table 4.   Table 4, by compensating for the gyro bias that changes according to position of the IMU, the latitude error was improved by 0.53% and the longitude error was improved by 59.44%. Therefore, the position error of RINS could be effectively improved by compensating for the gyro bias that changes according to position of the IMU. The navigation simulation results show that the up-axis average attitude error when not compensating for the gyro bias that changes according to position of the IMU is 0.0338 deg and the up-axis average attitude error when compensating for the gyro bias that changes according to position of the IMU is 0.0302 deg. Consequently, the up-axis attitude error can be improved by 10.65% by compensating for the gyro bias that changes according to the position of the IMU.
The latitude error and longitude error obtained by performing the navigation simulation are shown in Figure 3 and the horizontal axis position error at the end of the navigation is summarized in Table 4. The position error obtained by performing a navigation simulation according to the conditions in Table 3 is shown in Figure 2. The navigation simulation results show that the up-axis average attitude error when not compensating for the gyro bias that changes according to position of the IMU is 0.0338 and the up-axis average attitude error when compensating for the gyro bias that changes according to position of the IMU is 0.0302 . Consequently, the up-axis attitude error can be improved by 10.65% by compensating for the gyro bias that changes according to the position of the IMU.
The latitude error and longitude error obtained by performing the navigation simulation are shown in Figure 3 and the horizontal axis position error at the end of the navigation is summarized in Table 4.   Table 4, by compensating for the gyro bias that changes according to position of the IMU, the latitude error was improved by 0.53% and the longitude error was improved by 59.44%. Therefore, the position error of RINS could be effectively improved by compensating for the gyro bias that changes according to position of the IMU.  As shown in Table 4, by compensating for the gyro bias that changes according to position of the IMU, the latitude error was improved by 0.53% and the longitude error was improved by 59.44%. Therefore, the position error of RINS could be effectively improved by compensating for the gyro bias that changes according to position of the IMU.

Identifyning Gyro Bias Change Caused by Position of the IMU
This section describes techniques for identifying and compensating for gyro biases that vary with the location of the IMU. In this paper, a system-level indirect calibration technique was adopted to identify the gyro bias that changes according to the position of the IMU and compensated according to the IMU position condition.
In system-level indirect calibration techniques, RINS performs navigation for a short period of time in a static state [33]. Under static state and short-term navigation conditions, it can be assumed that the effect of the transport rate and Coriolis effects on the navigation performance of RINS is negligibly small. Based on the above assumption, the position error of RINS in the process of performing the system-level indirect calibration technique is obtained again as follows using Equation (1).
By assuming that the effect of the earth rate on the position error of RINS is very small and residual gyro/accelerometer errors except gyro bias are sufficiently small, the position error of Equation (12) is summarized as shown in Equation (13).
In a similar manner, the acceleration error of RINS in the process of performing the system-level indirect calibration technique is obtained again as follows: By assuming that the acceleration error of RINS does not change over time, the first derivative of the acceleration error over time in Equation (14) could be calculated as shown in Equation (15).
The gyro bias that changes according to the position of the IMU may be calculated using the up-axis position error of Equation (13) and the horizontal axis acceleration error change rate of Equation (15). In this paper, system-level indirect calibration was repeatedly performed in the Z-axis up position and Z-axis down position to obtain the gyro bias that changes according to the position of the IMU.
The rotation sequence for calculating the gyro bias at the Z-axis up position consists of self-alignment at the initial position, navigation at the initial position, and navigation at the post-rotation position. The navigation time for the initial position and the post-rotation position is 10 min and the position of RINS is shown in Figure 4.  In this paper, analysis of the rotation sequence applies sum frames including body frame and navigation frame. The detailed descriptions of these frames are prese in Section 2. In Figure 4, the azimuth angle of the initial position is 0 degrees and azimuth angle of the post-rotation position is 180 degrees. The coordinate transform matrix for the initial position and the post-rotation position is as shown in Equation , could be determined using the second derivative of the north-axis velocity and , could be determined using the second derivative of the east-axis velocity e , could be calculated using the rates of up-axis attitude error before and after rota The gyro bias for the Z−axis down position may be obtained in the same mann the gyro bias for the Z−axis up position. The rotation sequence for obtaining the gyro at the Z-axis down position is shown in Figure 5. In this paper, analysis of the rotation sequence applies sum frames including the body frame and navigation frame. The detailed descriptions of these frames are presented in Section 2. In Figure 4, the azimuth angle of the initial position is 0 degrees and the azimuth angle of the post-rotation position is 180 degrees. The coordinate transformation matrix for the initial position and the post-rotation position is as shown in Equation (16).
The gyro bias for the Z-axis up position is obtained by substituting Equation 16 into Equations (13) and (15).
β x,u could be determined using the second derivative of the north-axis velocity error and β y,u could be determined using the second derivative of the east-axis velocity error. β z,u could be calculated using the rates of up-axis attitude error before and after rotation.
The gyro bias for the Z−axis down position may be obtained in the same manner as the gyro bias for the Z−axis up position. The rotation sequence for obtaining the gyro bias at the Z-axis down position is shown in Figure 5.
In Figure 5, the coordinate transformation matrix for the initial position and the post-rotation position is as shown in Equation (18). In Figure 5, the coordinate transformation matrix for the initial position and the rotation position is as shown in Equation (18). (13) and (15). The gyro bias that changes according to the position of the IMU may be obt using Equations (17) and (19). The gyro bias that changes according to the position IMU is repeatedly identified for various temperature conditions and positional cond of the IMU.

The gyro bias for the Z-axis down position is obtained by substituting Equatio into Equations
The gyro bias that changes according to the position of the IMU obtained by per ing the above process is compensated in real time in the operation process of the The process of realizing the proposed calibration method and verification tests a scribed in Section 4.

Experiment Results
In this section, the accuracy of the proposed calibration method in Section 3 is ve by performing stationary state tests. A 1 Nm/h class of RINS consisting of a three-axi and a three-axis silicon accelerometer is used for the stationary state test. This pape structed the IMU by adopting the sensor specification in Table 3. The RLG arrange in the IMU for the Z-axis up position and Z-axis down position is described in Fig  During the operation of RINS, the IMU rotates according to the dual-axis 16-positi tation scheme presented in [29]. Gyro bias, which varies with the location of the IMU identified by performing system-level indirect calibration under various temperatur ditions. The gyro bias for the Z-axis down position is obtained by substituting Equation (18) into Equations (13) and (15).
The gyro bias that changes according to the position of the IMU may be obtained using Equations (17) and (19). The gyro bias that changes according to the position of the IMU is repeatedly identified for various temperature conditions and positional conditions of the IMU.
The gyro bias that changes according to the position of the IMU obtained by performing the above process is compensated in real time in the operation process of the RINS. The process of realizing the proposed calibration method and verification tests are described in Section 4.

Experiment Results
In this section, the accuracy of the proposed calibration method in Section 3 is verified by performing stationary state tests. A 1 Nm/h class of RINS consisting of a three-axis RLG and a three-axis silicon accelerometer is used for the stationary state test. This paper constructed the IMU by adopting the sensor specification in Table 3. The RLG arrangement in the IMU for the Z-axis up position and Z-axis down position is described in Figure 6. During the operation of RINS, the IMU rotates according to the dual-axis 16-position rotation scheme presented in [29]. Gyro bias, which varies with the location of the IMU, was identified by performing system-level indirect calibration under various temperature conditions.

Experiment Setup
This section verifies the performance of the calibration technique proposed in this paper by performing stationary state tests. A 1 Nm/h class of RINS consisting of a three-axis RLG and a three-axis silicon accelerometer is used for the stationary state test. During the operation of RINS, the IMU rotates according to the dual-axis 16-position rotation scheme presented in [29]. Gyro bias, which varies with the location of the IMU, was identified by performing system-level indirect calibration under various temperature conditions. The system-level indirect calibration performed to identify the gyro bias depending on the position of IMU is performed in the same order as shown in Figures 6 and 7. The data processing PC is responsible for collecting navigation data and calculating calibration coefficients. The DC power supply supplies the power required for the system-level indirect calibration process. RINS is installed inside the temperature chamber of the rate table and rotates according to the input of the rate table control console.
This section verifies the performance of the calibration technique proposed in this paper by performing stationary state tests. A 1 Nm/h class of RINS consisting of a threeaxis RLG and a three-axis silicon accelerometer is used for the stationary state test. During the operation of RINS, the IMU rotates according to the dual-axis 16-position rotation scheme presented in [29]. Gyro bias, which varies with the location of the IMU, was identified by performing system-level indirect calibration under various temperature conditions.
The system-level indirect calibration performed to identify the gyro bias depending on the position of IMU is performed in the same order as shown in Figures 6 and 7. The data processing PC is responsible for collecting navigation data and calculating calibration coefficients. The DC power supply supplies the power required for the system-level indirect calibration process. RINS is installed inside the temperature chamber of the rate table and rotates according to the input of the rate table control console. The system step indirect calibration test is performed on the Z-axis up position condition of the IMU and the Z-axis down position condition. The calibration test is repeated 37 times in a temperature range between 6.54 °C and 58.14 °C. The rotational angular velocity of the IMU is 2 deg/s. During the process of the calibration test, navigation was conducted for about 1500 s. Gyro bias depending on the position of IMU is calculated using Equations (17) to (19). To reflect the temperature effect on the gyro bias depending on the position of the IMU, the calculation process is repeated for various temperature conditions. Figure 8 shows that the gyro bias changes according to the position of the IMU identified by performing the calibration test. The system step indirect calibration test is performed on the Z-axis up position condition of the IMU and the Z-axis down position condition. The calibration test is repeated 37 times in a temperature range between 6.54 • C and 58.14 • C. The rotational angular velocity of the IMU is 2 deg/s. During the process of the calibration test, navigation was conducted for about 1500 s.
Gyro bias depending on the position of IMU is calculated using Equations (17) to (19). To reflect the temperature effect on the gyro bias depending on the position of the IMU, the calculation process is repeated for various temperature conditions. Figure 8 shows that the gyro bias changes according to the position of the IMU identified by performing the calibration test. Gyro bias depending on the position of IMU is calculated using Equations (17) to (19). To reflect the temperature effect on the gyro bias depending on the position of the IMU, the calculation process is repeated for various temperature conditions. Figure 8 shows that the gyro bias changes according to the position of the IMU identified by performing the calibration test. As shown in Figure 8, in the X-axis and Z-axis, the effect of the gyro bias according to the position of the IMU is large. On the other hand, in the Y-axis, the effect of the gyro bias according to the position of the IMU is hardly shown.
These identified gyro bias sets are modeled as compensation functions based on third-order polynomial. The compensation functions are utilized for real-time compensation. Therefore, navigation performance of the RINS can be improved using an appropriate compensation function depending on the position of the IMU.
For RINS used in this paper, the X-axis gyro and the Y-axis gyro are arranged similarly in the Z-axis up position and Z-axis down position. However, the identified gyro bias changed by 0.015 deg/h on the X-axis. In the case of the Z-axis gyro, heat transfer by gravity from the electric board to the mono block can be expected to have an effect, but As shown in Figure 8, in the X-axis and Z-axis, the effect of the gyro bias according to the position of the IMU is large. On the other hand, in the Y-axis, the effect of the gyro bias according to the position of the IMU is hardly shown.
These identified gyro bias sets are modeled as compensation functions based on thirdorder polynomial. The compensation functions are utilized for real-time compensation. Therefore, navigation performance of the RINS can be improved using an appropriate compensation function depending on the position of the IMU.
For RINS used in this paper, the X-axis gyro and the Y-axis gyro are arranged similarly in the Z-axis up position and Z-axis down position. However, the identified gyro bias changed by 0.015 deg/h on the X-axis. In the case of the Z-axis gyro, heat transfer by gravity from the electric board to the mono block can be expected to have an effect, but only 0.004 deg/h has changed. This difference in gyro bias appears differently depending on each RINS. Therefore, it can be interpreted that the calibration method proposed in this paper can improve the navigation performance of RINS.
The gyro bias identified under the Z-axis up position condition and the Z-axis down position condition varies by about 0.015 deg/h depending on the position of the IMU on the X-axis and by about 0.004 deg/h depending on the position of the IMU on the Z-axis. As analyzed in Section 3.1, the gyro bias that changes according to the position of the IMU may degrade the navigation performance of RINS. Therefore, it is possible to improve navigation performance by compensating for the gyro bias that changes depending on the position of the IMU identified through calibration tests during the operation of the RINS.

Results of the Static State Navigation Test
In this section, we validate the improved navigation performance by compensating for the gyro bias that changes according to the position of the IMU. In this paper, longterm stationary navigation was performed to determine the navigation performance of RINS. Long-term stationary navigation tests are designed to store the gyro bias obtained in Section 4.1 on the navigation computer of RINS, compensate for the gyro bias that changes according to the position of the IMU in real time, and perform the navigation. For the static-state navigation test, RINS is mounted on a leveled surface and performs the navigation for about 7 days. In this paper, the navigation performances of the proposed calibration method and conventional calibration method were evaluated based on the horizontal axis position error.
The results of long stationary navigation tests are used to assess the effect of gyro biases depending on the position of the IMU on navigation performance. In this paper, navigation performance was determined using the latitude error, hardness error, position error, and time root mean square (tRMS) position error at the end of the stationary navigation test [37]. The tRMS position error is calculated as in Equation (20).
where t n means navigation time and E i means the horizontal-axis position error for t i . Latitude error and longitude error obtained by performing the stationary state test for a long time are shown in Figure 9 and horizontal axis position error and tRMS error are shown in Figure 10. The navigation performance obtained by the stationary navigation test is summarized in Table 5.

Results of the Static State Navigation Test
In this section, we validate the improved navigation performance by compens for the gyro bias that changes according to the position of the IMU. In this paper, term stationary navigation was performed to determine the navigation performan RINS. Long-term stationary navigation tests are designed to store the gyro bias obta in Section 4.1 on the navigation computer of RINS, compensate for the gyro bias changes according to the position of the IMU in real time, and perform the navigation the static-state navigation test, RINS is mounted on a leveled surface and perform navigation for about 7 days. In this paper, the navigation performances of the prop calibration method and conventional calibration method were evaluated based on the izontal axis position error.
The results of long stationary navigation tests are used to assess the effect of biases depending on the position of the IMU on navigation performance. In this p navigation performance was determined using the latitude error, hardness error, pos error, and time root mean square (tRMS) position error at the end of the stationary gation test [37]. The tRMS position error is calculated as in Equation (20).
where means navigation time and means the horizontal-axis position error fo Latitude error and longitude error obtained by performing the stationary stat for a long time are shown in Figure 9 and horizontal axis position error and tRMS are shown in Figure 10. The navigation performance obtained by the stationary navig test is summarized in Table 5.     The long-term stationary state test results show that the calibration technique proposed in this paper improves the navigation performance of RINS. Based on the navigation performance at the end of the long-term stationary state test, the tRMS position error of 1.2921 Nm was measured by applying the calibration technique proposed in this paper. On the other hand, when the gyro bias that changes according to the position of the IMU was not compensated for, a tRMS position error of 2.4603 Nm was measured.
On average, the navigation performance was improved by 69.83% by identifying and compensating for the gyro bias, which changes according to the position of the IMU. Applying the calibration technique proposed in this paper, the horizontal-axis position error was improved by 79.86% and the tRMS position error was improved by 47.48%. In particular, it was confirmed that the longitude error was improved by 82.12% through compensation for the gyro bias, which changes according to the position error of the IMU.

Conclusions
In this paper, to improve the navigation performance of RINS, we propose a calibration technique that identifies and compensates for the gyro bias that changes according to the position of the IMU. First, RINS was designed using RLG and IMU was rotated by importing the dual-axis 16-position rotation scheme. The effect of the gyro bias that changes according to the position of the IMU on the navigation performance of RINS designed in this paper was theoretically analyzed and verified through navigation simulation.
Gyro bias, which changes according to the location of the IMU, was identified by applying a system-level indirect calibration technique. By identifying the gyro bias for the Z-axis up position condition and the Z-axis down position condition of the IMU, the gyro bias that changes according to the position of the IMU was obtained. The gyro bias that changes according to the position of the IMU identified by repeating system-level indirect calibration under various temperature conditions is compensated for in real time during the operation of RINS.
The performance of the proposed technique was verified by performing long-term stationary navigation tests. A long-term stationary navigation test was performed according to the compensation condition of the gyro bias that changes depending on the position of the IMU, and the longitude error, latitude error, horizontal-axis position error, and tRMS position error were measured. As a result of the test, it was confirmed that the navigation performance of RINS was improved by up to 82.12% by applying the calibration technique proposed in this paper. Thus, it is concluded that the proposed calibration technique sufficiently improves the navigation performance of RINS.
Finally, the calibration technique proposed in this paper performed system-level indirect calibration only for the Z-axis up position condition and the Z-axis down position condition to identify the gyro bias that changes according to the position of the IMU. Therefore, future studies will improve the navigation performance of RINS in a way that identifies and compensates for the gyro bias that changes according to the position of the IMU for more diverse positional conditions.