In-Motion Alignment Method of SINS Based on Improved Kalman Filter under Geographic Latitude Uncertainty

: To realize the in-motion alignment of the strapdown inertial navigation system (SINS) under the geographic latitude uncertainty, we propose a latitude estimation and in-motion alignment method based on the integral dynamic window and polynomial ﬁtting (IDW-PF) and improved Kalman ﬁlter (IKF). First, the integral dynamic window (IDW) is designed to smooth out the high-frequency line motion interference and accelerometer noise. Second, the speciﬁc force integral is performed for a cubic polynomial ﬁtting (PF) with time as an independent variable to further suppress the line motion interference. Simultaneously, the latitude is estimated according to the geometric relationship between the angle of the gravitational acceleration vectors at different moments and the latitude. Finally, the IKF based on the multi-fading factor is designed for the in-motion alignment of SINS. A simulation experiment is conducted to verify the proposed latitude estimation and in-motion alignment method. The results indicate that the latitude can be estimated well by the method based on the IDW-PF; the mean and standard deviation of the estimated latitude can achieve − 0.016 ◦ and 0.013 ◦ within 300 s. The trapezoidal maneuvering path is optimal when IKF is used, the pitch error is 0.0002 ◦ , the roll error is 0.0009 ◦ and the heading error is − 0.0047 ◦ after the alignment ends at 900 s.


Introduction
The initial alignment is one of the key technologies for the strapdown inertial navigation system (SINS), and its accuracy and rapidity directly affect the performance of the navigation system [1].As reported in [2], some of the challenges associated with the stationary base alignment of SINS may have been resolved.However, this alignment has high requirements on the environment, which greatly reduces its agility and its anti-strike feature.Because the in-motion alignment method can effectively improve the carrier's maneuverability under the premise of ensuring the alignment accuracy, this technology has received more attention recently [3][4][5][6][7][8].The traditional initial alignment method of the SINS is generally carried out under the condition that the geographical location (mainly refers to latitude) is accurately known.However, in some environments, such as underwater, deep mountains and forests, deep tunnels, or when the carrier is parked in the hangar and cannot receive satellite positioning signals, there is also the need for initial alignment.Based on the above application scenarios, researching the initial alignment under latitude uncertainty has great theoretical significance and practical application value.• In recent years, G. Yan et al. [9] took the lead in proposing a latitude calculation method based on the output of inertial measurement components under the stationary base.The authors concluded that the initial alignment can be completed by using the estimation latitude information.Y. Wang et al. [10] used the gravitational acceleration vectors at two different moments to estimate the latitude value.It was shown that the initial alignment under dynamic conditions can be completed by effectively isolating the sloshing angular motion interference.The accelerometer will also be disturbed by the line vibration in addition to its instrument errors under the in-motion condition.As a result, the output error of the accelerometer is larger, which will have a greater impact on the latitude estimation.However, the above latitude estimation methods are applied to stationary bases and swaying bases.When the carrier is in motion, the above methods are not applicable.
For an in-motion environment, many researchers have studied the initial alignment with the latitude certainty.Unlike the stationary base alignment, external equipment is usually required to provide auxiliary information (such as the position or the velocity) to the in-motion alignment.The most popular aiding information for SINS is the global positioning system (GPS) [11][12][13][14][15][16][17].However, its availability and reliability are seriously affected by jamming, blocking and spoofing of the GPS signals.In response to this problem, odometers [18] or Doppler velocity log (DVL) [19][20][21][22][23] aided the in-motion alignment was proposed by many scholars.For these velocity-aided in-motion alignment conditions, if the geographic latitude is not available, the attitude cannot be determined.Hence, relying on the inertial component output of the SINS for the latitude self-estimation is very meaningful.It is of great significance to study how to conduct the in-motion alignment of SINS under geographic latitude uncertainty.It can improve the maneuverability and adaptability to the environment of the carrier.
The main contributions of this paper are summarized as follows: (1) In the latitude estimation, a method based on the integral dynamic window and polynomial fitting (IDW-PF) for SINS is proposed.The measurement value of the accelerometer is with integral transformation, and it is projected to the inertial system.Polynomial fitting (PF) is performed on the result of the specific force integration to achieve the second smoothing of the interference acceleration.The latitude can be estimated in real-time and accurately because of making full use of the inertial instrument data.(2) In the in-motion alignment, the improved Kalman filter (IKF) based on the multi-fading factor is designed to calculate the in-motion misalignment angles.Simulation experiments are performed to verify the performances of the proposed in-motion alignment method under latitude uncertainty.The results show that the trapezoidal path is with the best alignment effect based on the IKF.
The rest of this paper is organized as follows.The latitude estimation algorithm based on IDW-PF under the in-motion condition is proposed in Section 2. In Section 3, the IKF based on the multi-fading factor is designed.The in-motion alignment model is established in Section 4. The simulation results are given in Section 5. Finally, the conclusions are drawn in Section 6.

Latitude Estimation Based on IDW-PF under the In-Motion Condition
As shown in Figure 1, Point O is the center of the latitude circle where the SINS is located.Point A and B are the positions in the inertial coordinate system at time t 1 and t 2 during the initial alignment, respectively.Point O is the intersection point of the vertical line of gravitational acceleration and the Earth's axis [9].g i b0 (t 1 ) and g i b0 (t 2 ) represent the gravitational acceleration vectors at t 1 and t 2 , respectively.θ denotes the angle between the vector g i b0 (t 1 ) and g i b0 (t 2 ).α indicates the angle of the Earth turning around from t 1 to t 2 .The angle between the latitude circle of point O and the vertical line of gravity is the latitude L. The relationship between the lengths of the geometric line segments can be obtained as: The latitude value can be estimated by ( 1) and ( 2): The angular velocity ie  of the Earth's rotation is known; hence, the angle  of rotation can be expressed as follows:

()
ie tt To reduce the influence of the actual inertial navigation specific force measurement noise and the interference movement of the base, the angle is generally calculated after smoothing with the specific force integral over a period of time.Then, it can be obtained as: where   The relationship between the lengths of the geometric line segments can be obtained as: The latitude value can be estimated by ( 1) and ( 2): The angular velocity ω ie of the Earth's rotation is known; hence, the angle α of rotation can be expressed as follows: To reduce the influence of the actual inertial navigation specific force measurement noise and the interference movement of the base, the angle is generally calculated after smoothing with the specific force integral over a period of time.Then, it can be obtained as: where g i b0 (t)dt, ∆t indicates the integration time.g i b0 (t) is calculated as: where g b (t) is the gravitational acceleration vector in b coordinate system at time t.
b can be updated by measuring the gyroscope information in real-time: .
where ω b ib is the output value of gyroscope and ω b ib × denotes the antisymmetric matrix of ω b ib .
Hence, the key to latitude estimation is to extract more accurate g i b0 in real-time.However, under the in-motion condition, besides its instrument error, the accelerometer is also disturbed by the line vibration.This leads to a large output error of the accelerometer, which has a large impact on the latitude estimation.In this paper, a method of latitude estimation based on IDW-PF under the in-motion condition is proposed.The algorithm diagram is shown in Figure 2 [24].
is also disturbed by the line vibration.This leads to a large output error of the accelerometer, which has a large impact on the latitude estimation.In this paper, a method of latitude estimation based on IDW-PF under the in-motion condition is proposed.The algorithm diagram is shown in Figure 2 [24].
That is, as long as there is the output of inertial instrument data, the current latitude estimation value can be obtained.The principle of integral dynamic window (IDW) is similar to that of dynamic window (DW).
In summary, the latitude estimation algorithm based on IDW-PF under the in-motion condition can be summarized as follows: (1) The integral value ) tt  by using (5).( 5) Calculate the angle α at time 2M t  by using (4).(6) The real-time latitude estimation value at time 2M t  can be obtained by using (3).

The Improved Kalman Filtering Algorithm
The state equation and observation equation of the stochastic linear discrete system are given as [25]: In Figure 2, −f i b0 (t N ) indicates the calculated gravitational acceleration projected on the inertial coordinate system at time t N .It can be calculated by updating the accelerometer output f i b and the gyroscope output ω b ib at the current time.According to ( 3)-( 7), take f i b0 (t 1 ) obtained from the calculation value at time t 1 as a reference.The real-time latitude estimation value at time t N can be obtained from the calculation value That is, as long as there is the output of inertial instrument data, the current latitude estimation value can be obtained.The principle of integral dynamic window (IDW) is similar to that of dynamic window (DW).In summary, the latitude estimation algorithm based on IDW-PF under the in-motion condition can be summarized as follows: (1) The integral value v dt, which is used as the reference of the IDW.
(2) Calculate the integral value of v (3) Perform polynomial fitting with time t as the independent variable to the calculation The angle θ at time t 2M can be obtained from the integral value v (5) Calculate the angle α at time t 2M by using ( 4). ( 6) The real-time latitude estimation value at time t 2M can be obtained by using (3).

The Improved Kalman Filtering Algorithm
The state equation and observation equation of the stochastic linear discrete system are given as [25]: where X k is the state matrix at t k , Φ k,k−1 is the one-step transition matrix from t k−1 to t k , Γ k,k−1 is the system noise driving matrix, Z k is the system measurement sequence, H k is the measurement matrix, V k is the measurement noise sequence and W k is the system incentive noise sequence.The standard KF recursive equations are as follows [26]: Remote Sens. 2022, 14, 2581 5 of 15 (3) Filter gain (4) One-step prediction mean square error (5) Estimated mean square error where Xk,k−1 is the state one-step prediction matrix, Xk−1 is the estimated state matrix at t k−1 , Xk is the estimated state matrix at t k , K k is the filter gain matrix, P k,k−1 is the one-step prediction mean square error matrix, R k is the measurement noise variance matrix, P k−1 is the error covariance matrix of the optimal filter value at t k−1 , Q k−1 is system noise variance matrix at t k−1 , P k is the estimated mean square error matrix and I is the unit matrix.
During the calculation process of KF, filtering anomalies or even divergences often occur.The main reasons are the inaccurate system models, the inaccurate noise statistical models, the accumulation of the rounding errors, etc.They cause the current measurement value to reduce the correction effect of the estimation value and the old measurement value to increase the correction effect of the estimation value [27].The proposed IKF introduces a diagonal matrix fading factor into (12) so that the residual sequences at different times remain orthogonal everywhere, namely where D k+1 is the diagonal matrix with multiple fading factors, and it can be determined by the following method [27]: where are the coefficients predetermined by the prior knowledge.It can increase the corresponding α k of the components that are prone to sudden changes.If the system has no prior knowledge, take α k = 1.In (16), C k can be expressed as where tr[•] denotes the trace of the matrix "•".
and it is the weakening factor.V 0,k is the mean square error.
The multi-fading factor KF is used to fade each data channel at different rates.It can better adjust the gain matrix.Even when the system reaches a steady state, the gain matrix can be adjusted adaptively to enhance the robustness of the model mismatch and the system destabilization.More accurate and stable estimation results will be obtained.

The Proposed In-Motion Alignment Method
The IKF-based alignment technology for SINS requires auxiliary navigation equipment to complete the alignment.Considering that the external velocity of the vehicle is easy to obtain, this paper focuses on the use of external reference velocity-aided alignment methods.The difference between the velocity output by SINS and the velocity provided by the auxiliary navigation equipment is used as the observation.Attitude misalignment angle and inertial device error are estimated by IKF so as to complete the initial alignment.
The instability of the vertical channel is considered; the upward velocity and the z-axis accelerometer constant bias are not regarded as the state variables in this paper.Hence, the state vectors selected in the study are as follows: where δV E is the eastward velocity error, δV N is the northward velocity error, φ E is the eastward misalignment angle, φ N is the northward misalignment angle, φ U is the upward misalignment angle, δL is the latitude error, δλ is the longitude error, ∇ x is the x-axis accelerometer bias, ∇ y is the y-axis accelerometer bias, ε x is the x-axis gyroscope constant drift, ε y is the y-axis gyroscope constant drift and ε z is the z-axis gyroscope constant drift.
The state equation of the system is expressed as [28,29]: .
where X(t) is the system state vector, F(t) is the system matrix, G(t) is the system noise matrix and w(t) is the system noise vector.The system matrix F(t) is expressed as: where where L is the latitude, f E , f N and f U are the eastward acceleration, northward acceleration and upward acceleration, respectively.V E , V N and V U are the eastward velocity, northward velocity and upward velocity, respectively.ω ie is the angular velocity of the Earth's rotation.R e is the Earth's semi-major axis, and R n is the radius of principal curvature.C n b is the attitude matrix, and all nonzero elements of the matrix T(t) are expressed as: The horizontal velocity of the SINS is assumed to be expressed as [30]: where V ES and V NS are the eastward velocity and the northward velocity of the SINS, respectively.V ET and V NT are the true eastward velocity and northward velocity of the carrier, respectively.δV ES and δV NS are the eastward velocity error and the northward velocity error of the SINS, respectively.
The external reference velocity is assumed to be expressed as [30]: where V ED and V ND are the eastward velocity and the northward velocity of the external equipment, respectively.δV ED and δV ND are the eastward velocity error and the northward velocity error of the external equipment, respectively.Hence, the measurement equation of the system can be expressed as: where So far, the model of velocity-aided in-motion alignment has been established.Then, the IKF can be used for filtering, and the process is introduced as follows: (1) State prediction (2) Calculate the diagonal matrix with multiple suboptimal fading factors, where where 0 ≤ ρ ≤ 1, and it is the forgetting factor.(3) Calculate the prediction state covariance matrix using ( 14).(4) Calculate the gain matrix using (11).
(5) Calculate the state estimation and its covariance matrix at t.
At this point, the in-motion alignment method based on IKF under geographic latitude uncertainty has been implemented.Then the initial attitude can be obtained.

Simulation Results and Discussions
In order to verify the latitude estimation and in-motion alignment method proposed in this paper, firstly, the simulation experiments of the latitude estimation based on the IDW-PF are carried out.Then, the simulation experiments of in-motion alignment based on the IKF are performed for various motion paths.The external reference velocity error is not considered when studying the influence of different maneuver tracks on the in-motion alignment; here, the external reference velocity error is set to zero.The latitude used in the simulation experiment is the latitude estimation value based on the IDW-PF.
The specifications of the inertial measurement unit (IMU) are listed in Table 1.The instrument output frequency is 200 Hz.At the time of the experiment, the system has completed the coarse alignment, and the pitch, roll and heading attitude errors at the beginning of fine alignment are 0.1 • , 0.1 • and 1 • , respectively.The linear velocity of the carrier is 2 m/s; the turn angle rate is 0.9 • /s.The simulation time of latitude estimation and in-motion are set as 300 s and 900 s, respectively.The time parameters are set as: t N = 8 s and t M = 4 s.In this subsection, three kinds of latitude estimation algorithms are performed, which are based on the DW, IDW and IDW-PF, respectively.Figure 3 shows the curves of the latitude estimation error, δL denotes the latitude error, and the statistical results of the latitude estimation error are listed in Table 2.

Simulation Experiment of the Latitude Estimation Based on the IDW-PF
In this subsection, three kinds of latitude estimation algorithms are performed, which are based on the DW, IDW and IDW-PF, respectively.Figure 3 shows the curves of the latitude estimation error, L  denotes the latitude error, and the statistical results of the latitude estimation error are listed in Table 2.It can be seen from Figure 3 and Table 2 that all three methods can realize self-estimation of latitude.However, due to the influence of accelerometer noise, the standard deviation (SD) of the DW method is relatively large.The IDW method effectively suppresses the influence of the accelerometer noise on the estimation accuracy through the integral effect.The IDW-PF method achieves the second smoothing of the interference acceleration by performing the PF on the result of the specific force integration.Hence, the IDW-PF method has the best performance with the fastest convergence speed and the highest estimation accuracy.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Uniform Linear Motion
The initial attitude angles are assumed to be as follows: the pitch angle is 0 , the roll angle is 0 , the heading angle is 0 , 90 , 180 and 270 .Figure 4 shows the attitude error of the in-motion alignment of uniform linear motion when the heading angle is 90 .Figure 5 shows the estimation values of the gyroscope and accelerometer constant error when the heading angle is 90 . ,  and  are pitch error, roll error  It can be seen from Figure 3 and Table 2 that all three methods can realize selfestimation of latitude.However, due to the influence of accelerometer noise, the standard deviation (SD) of the DW method is relatively large.The IDW method effectively suppresses the influence of the accelerometer noise on the estimation accuracy through the integral effect.The IDW-PF method achieves the second smoothing of the interference acceleration by performing the PF on the result of the specific force integration.Hence, the IDW-PF method has the best performance with the fastest convergence speed and the highest estimation accuracy.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Uniform Linear Motion
The initial attitude angles are assumed to be as follows: the pitch angle is 0 • , the roll angle is 0 • , the heading angle is 0 • , 90 • , 180 • and 270 • .Figure 4 shows the attitude error of the in-motion alignment of uniform linear motion when the heading angle is 90 • .Figure 5 shows the estimation values of the gyroscope and accelerometer constant error when the heading angle is 90 • .δθ, δγ and δφ are pitch error, roll error and heading error, respectively.∆ x and ∆ y are accelerometer bias of x-axis and y-axis, respectively.ε x and ε y are gyroscope drift of x-axis and y-axis, respectively.The system is open-loop during the alignment process, the alignment ends at 900 s, and then the attitude is compensated.Table 3 shows the statistical results of the attitude error, the estimation values of the gyroscope and accelerometer constant error when the heading angle is 0 • , 90 attitude is compensated.Table 3 shows the statistical results of the attitude error, the estimation values of the gyroscope and accelerometer constant error when the heading angle is 0 , 90 , 180 and 270 .attitude is compensated.Table 3 shows the statistical results of the attitude error, the estimation values of the gyroscope and accelerometer constant error when the heading angle is 0 , 90 , 180 and 270 .It can be seen from Figures 4 and 5 and Table 3 that when the initial heading angle is 90 , the pitch, roll and heading errors measured at 900 s are sufficiently low and are  It can be seen from Figures 4 and 5 and Table 3 that when the initial heading angle is 90 • , the pitch, roll and heading errors measured at 900 s are sufficiently low and are 0.0011 • , −0.0022 • and 0.0345 • , respectively.The bias of the x-axis and y-axis accelerometer is 1.036 µg and −1.550 µg, respectively.The constant drift of x-axis and y-axis gyroscopes is 0.073 • /h and 0.011 • /h, respectively.Although the alignment accuracy can meet the alignment requirements, there are large errors in estimating the accelerometer constant bias and gyroscope constant drift.It can be seen from Table 3 that the above conclusions are true for multiple heading angles.Moreover, the estimation accuracy of the misalignment angle is stable under the condition of uniform linear motion maneuvering.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (Circular Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 • , the roll angle is 0 • and the heading angle is 0 • .When the influence of the circular maneuver mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 6.The results of the attitude error and the constant error estimation for the gyroscope and accelerometer are shown in Figures 7 and 8, respectively.
Remote Sens. 2022, 14, x FOR PEER REVIEW 11 of 16 meet the alignment requirements, there are large errors in estimating the accelerometer constant bias and gyroscope constant drift.It can be seen from Table 3 that the above conclusions are true for multiple heading angles.Moreover, the estimation accuracy of the misalignment angle is stable under the condition of uniform linear motion maneuvering.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (Circular Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 , the roll angle is 0 and the heading angle is 0 .When the influence of the circular maneuver mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 6.The results of the attitude error and the constant error estimation for the gyroscope and accelerometer are shown in Figure 7 and Figure 8, respectively.meet the alignment requirements, there are large errors in estimating the accelerometer constant bias and gyroscope constant drift.It can be seen from Table 3 that the above conclusions are true for multiple heading angles.Moreover, the estimation accuracy of the misalignment angle is stable under the condition of uniform linear motion maneuvering.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (Circular Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 , the roll angle is 0 and the heading angle is 0 .When the influence of the circular maneuver mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 6.The results of the attitude error and the constant error estimation for the gyroscope and accelerometer are shown in Figure 7 and Figure 8, respectively.When the alignment ends at 900 s, the attitude error, accelerometer bias and gyroscope drift are measured.The pitch error is 0.0002 , the roll error is −0.0010 and the heading error is 0.0016 .The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 45.19 µg and 43.84 µ g, respectively.Moreover, they are close to the bias value of 50 µ g set by the simulation, which is explained that the accelerometer bias has been effectively estimated.The estimated constant drift of the x-axis and y-axis gyroscope is 0.043 /h and 0.061 /h, respectively.There is a certain gap between the constant drift of 0.006 /h set by the simulation, which is explained that the gyroscope constant drift is not well estimated.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (Trapezoidal Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 , the roll angle is 0 and the heading angle is 0 .When the influence of the trapezoidal maneuvering mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 9.The results of the attitude error and the constant error estimation for the gyroscope and accelerometer are shown in Figure 10 and Figure 11, respectively.When the alignment ends at 900 s, the attitude error, accelerometer bias and gyroscope drift are measured.The pitch error is 0.0002 • , the roll error is −0.0010 • and the heading error is 0.0016 • .The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 45.19 µg and 43.84 µg, respectively.Moreover, they are close to the bias value of 50 µg set by the simulation, which is explained that the accelerometer bias has been effectively estimated.The estimated constant drift of the x-axis and y-axis gyroscope is 0.043 • /h and 0.061 • /h, respectively.There is a certain gap between the constant drift of 0.006 • /h set by the simulation, which is explained that the gyroscope constant drift is not well estimated.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (Trapezoidal Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 • , the roll angle is 0 • and the heading angle is 0 • .When the influence of the trapezoidal maneuvering mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 9.The results of the attitude error and the constant error estimation for the gyroscope and accelerometer are shown in Figures 10 and 11 When the alignment ends at 900 s, the attitude error, accelerometer bias and gyroscope drift are measured.The pitch error is 0.0002 , the roll error is −0.0010 and the heading error is 0.0016 .The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 45.19 µg and 43.84 µ g, respectively.Moreover, they are close to the bias value of 50 µ g set by the simulation, which is explained that the accelerometer bias has been effectively estimated.The estimated constant drift of the x-axis and y-axis gyroscope is 0.043 /h and 0.061 /h, respectively.There is a certain gap between the constant drift of 0.006 /h set by the simulation, which is explained that the gyroscope constant drift is not well estimated.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (Trapezoidal Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 , the roll angle is 0 and the heading angle is 0 .When the influence of the trapezoidal maneuvering mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 9.The results of the attitude error and the constant error estimation for the gyroscope and accelerometer are shown in Figure 10 and Figure 11, respectively.The attitude error, accelerometer bias and gyroscope drift were measured at 900 s.The pitch error is 0.0002 , the roll error is 0.0009 and the heading error is −0.0047 .The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 51.09 µg and 48.51 µ g, respectively.Moreover, they are close to the bias value 50 µ g set by the simulation, which is explained that the accelerometer bias was effectively estimated.The estimated constant drift of x-axis and y-axis gyroscopes is 0.0095 /h and 0.0096 /h, respectively.There is a certain gap between the constant drift of 0.006 /h set by the simulation, which shows that there is a certain error in the constant drift of the gyroscope estimated in the maneuvering mode.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (S-Shaped Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 , the roll angle is 0 and the heading angle is 0 .When the influence of the s-shaped maneuvering mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 12.The results of attitude error and constant error estimation for the gyroscope and accelerometer are shown in Figures 13 and 14, respectively.The attitude error, accelerometer bias and gyroscope drift were measured at 900 s.The pitch error is 0.0002 , the roll error is 0.0009 and the heading error is −0.0047 .The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 51.09 µg and 48.51 µ g, respectively.Moreover, they are close to the bias value 50 µ g set by the simulation, which is explained that the accelerometer bias was effectively estimated.The estimated constant drift of x-axis and y-axis gyroscopes is 0.0095 /h and 0.0096 /h, respectively.There is a certain gap between the constant drift of 0.006 /h set by the simulation, which shows that there is a certain error in the constant drift of the gyroscope estimated in the maneuvering mode.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (S-Shaped Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 , the roll angle is 0 and the heading angle is 0 .When the influence of the s-shaped maneuvering mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 12.The results of attitude error and constant error estimation for the gyroscope and accelerometer are shown in Figures 13 and 14, respectively.The attitude error, accelerometer bias and gyroscope drift were measured at 900 s.The pitch error is 0.0002 • , the roll error is 0.0009 • and the heading error is −0.0047 • .The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 51.09 µg and 48.51 µg, respectively.Moreover, they are close to the bias value 50 µg set by the simulation, which is explained that the accelerometer bias was effectively estimated.The estimated constant drift of x-axis and y-axis gyroscopes is 0.0095 • /h and 0.0096 • /h, respectively.There is a certain gap between the constant drift of 0.006 • /h set by the simulation, which shows that there is a certain error in the constant drift of the gyroscope estimated in the maneuvering mode.

Simulation Experiment of Alignment Based on IKF When the Carrier Performs Turning Maneuver (S-Shaped Path)
The initial attitude angles are assumed to be as follows: the pitch angle is 0 • , the roll angle is 0 • and the heading angle is 0 • .When the influence of the s-shaped maneuvering mode on the in-motion alignment is considered, the maneuvering trajectory is shown in Figure 12.The results of attitude error and constant error estimation for the gyroscope and accelerometer are shown in Figures 13 and 14, respectively.The pitch error is 0.0003 , the roll error is 0.0011 and the heading error is 0.0039 when the alignment ends at 900 s.The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 49.02 µg and 46.68 µ g, respectively.Moreover, they are close to the bias value of 50 µ g set by the simulation, which is explained that the accelerometer bias has been effectively estimated.The estimated constant drift of the x-axis and y-axis gyroscope is 0.0088 /h and 0.0022 /h, respectively.There is a certain gap between the constant drift of 0.006 /h set by the sim-   The pitch error is 0.0003 , the roll error is 0.0011 and the heading error is 0.0039 when the alignment ends at 900 s.The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 49.02 µg and 46.68 µ g, respectively.Moreover, they are close to the bias value of 50 µ g set by the simulation, which is explained that the accelerometer bias has been effectively estimated.The estimated constant of the x-axis and y-axis gyroscope is 0.0088 /h and 0.0022 /h, respectively.There is a certain gap between the constant drift of 0.006 /h set by the sim-   The pitch error is 0.0003 , the roll error is 0.0011 and the heading error is 0.0039 when the alignment ends at 900 s.The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 49.02 µg and 46.68 µ g, respectively.Moreover, they are close to the bias value of 50 µ g set by the simulation, which is explained that the accelerometer bias has been effectively estimated.The estimated constant drift of the x-axis and y-axis gyroscope is 0.0088 /h and 0.0022 /h, respectively.There is a certain gap between the constant drift of 0.006 /h set by the sim- The pitch error is 0.0003 • , the roll error is 0.0011 • and the heading error is 0.0039 • when the alignment ends at 900 s.The misalignment angle is estimated well in the maneuvering mode.The bias of the x-axis and y-axis accelerometer is 49.02 µg and 46.68 µg, respectively.Moreover, they are close to the bias value of 50 µg set by the simulation, which is explained that the accelerometer bias has been effectively estimated.The estimated constant drift of the x-axis and y-axis gyroscope is 0.0088 • /h and 0.0022 • /h, respectively.There is a certain gap between the constant drift of 0.006 • /h set by the simulation, which shows that there is a certain error in the constant drift of the gyroscope estimated in the maneuvering mode.
The results of the alignment experiments between the three types of turning maneuvers are compared; it can be seen that, although the turning maneuver can have a positive effect on the estimation of the attitude misalignment angle and the IMU error, there are some differences in the effects of the motion path of the carrier, so the in-motion alignment effect is not only affected by the maneuver of the carrier but also by the path selection.Among the above three kinds of paths, the trapezoidal path has the best alignment effect, and the effect of the s-shaped path is better than that of the circular path.Hence, it is better to choose a trapezoidal path or an s-shaped path when practical conditions permit.

Conclusions
In-motion alignment method of SINS under geographic latitude uncertainty was proposed.The proposed method contains a method of latitude estimation and in-motion alignment method.Firstly, the latitude estimation algorithm based on the IDW-PT was designed.This IDW-PT can suppress and compensate for the interference of the carrier's line vibration.Secondly, the in-motion alignment method based on the IKF was devised.The multi-fading factor is introduced into the KF to fade each data channel at different rates.The gain matrix can be adjusted adaptively to enhance the robustness of the model mismatch and system destabilization and accordingly achieve more accurate and stable estimation results.Finally, the simulation experiment was shown to prove that the proposed method can estimate the geographic latitude accurately and applies to the in-motion initial alignment.The mean and standard deviation of the estimated latitude can achieve −0.016 • and 0.013 • within 300 s.Under turning maneuver conditions, the convergence accuracy and speed of the attitude misalignment angle and the inertial device error can be improved.Moreover, the trapezoidal path was shown to be the optimal maneuvering path.The pitch error is 0.0002 • , the roll error is 0.0009 • and the heading error is −0.0047 • after the alignment ends at 900 s.Based on the simulation performance of the proposed method, it is expected to be applied in real-time application.The method can enhance the navigation performance and the dynamic performance of the carrier.In the future, further investigations should be developed for the more serious situation such as the large misalignment angle, nonhorizontal path and multi-path combination.

Figure 1 .
Figure 1.Geometric relationship diagram for the latitude estimation.

ω
acceleration vector in b coordinate system at time t. is the output value of gyroscope and ( ) b ib  ω denotes the antisymmetric matrix of b ib ω .Hence, the key to latitude estimation is to extract more accurate 0 b i g in real-time.However, under the in-motion condition, besides its instrument error, the accelerometer

Figure 1 .
Figure 1.Geometric relationship diagram for the latitude estimation.

Figure 4 .
Figure 4. Attitude error of in-motion alignment of uniform linear motion.

Figure 4 .
Figure 4. Attitude error of in-motion alignment of uniform linear motion.

Figure 4 .
Figure 4. Attitude error of in-motion alignment of uniform linear motion.

Figure 7 .
Figure 7. Attitude error of in-motion alignment of circular maneuver.

Figure 7 .
Figure 7. Attitude error of in-motion alignment of circular maneuver.Figure 7. Attitude error of in-motion alignment of circular maneuver.

Figure 7 .
Figure 7. Attitude error of in-motion alignment of circular maneuver.Figure 7. Attitude error of in-motion alignment of circular maneuver.

Figure 10 .
Figure 10.Attitude error of in-motion alignment of trapezoidal maneuver.

Figure 13 .
Figure 13.Attitude error of in-motion alignment of s-shaped maneuver.

Figure 13 .
Figure 13.Attitude error of in-motion alignment of s-shaped maneuver.

Figure 13 . 16 Figure 12 .
Figure 13.Attitude error of in-motion alignment of s-shaped maneuver.

Figure 13 .
Figure 13.Attitude error of in-motion alignment of s-shaped maneuver.

Table 1 .
Parameters settings of IMU.
5.1.Simulation Experiment of the Latitude Estimation Based on the IDW-PF

Table 2 .
Statistical results of latitude estimation errors.

Table 2 .
Statistical results of latitude estimation errors.

Table 3 .
Statistical results of uniform linear motion.
yIt can be seen from Figures4 and 5and Table3that when the initial heading angle

Table 3 .
Statistical results of uniform linear motion.
y 

Table 3 .
Statistical results of uniform linear motion.