An Improved Strapdown Inertial Navigation System Initial Alignment Algorithm for Unmanned Vehicles

Along with the development of computer technology and informatization, the unmanned vehicle has become an important equipment in military, civil and some other fields. The navigation system is the basis and core of realizing the autonomous control and completing the task for unmanned vehicles, and the Strapdown Inertial Navigation System (SINS) is the preferred due to its autonomy and independence. The initial alignment technique is the premise and the foundation of the SINS, whose performance is susceptible to system nonlinearity and uncertainty. To improving system performance for SINS, an improved initial alignment algorithm is proposed in this manuscript. In the procedure of this presented initial alignment algorithm, the original signal of inertial sensors is denoised by utilizing the improved signal denoising method based on the Empirical Mode Decomposition (EMD) and the Extreme Learning Machine (ELM) firstly to suppress the high-frequency noise on coarse alignment. Afterwards, the accuracy and reliability of initial alignment is further enhanced by utilizing an improved Robust Huber Cubarure Kalman Filer (RHCKF) method to minimize the influence of system nonlinearity and uncertainty on the fine alignment. In addition, real tests are used to verify the availability and superiority of this proposed initial alignment algorithm.


Introduction
The unmanned vehicle has been used as a platform for aerial photogrammetry, marine monitoring, geodetic surveying, hazard state investigation and security protection based on different sensors equipped on it. The navigation information of the unmanned vehicle provided by its navigation system, including the position, the velocity, the heading and the horizontal attitude, is the premise and guarantee of the normal working. The precision of navigation system is directly related to its track accuracy and safety.
The Global Navigation Satellite System (GNSS) can provide three-dimensional navigation information including position and velocity of everywhere on the earth with high-accuracy by utilizing satellite signals. Caused by signal blocking and multipath effect, GNSS has poor location accuracy in urban areas and forests normally [1]. By utilizing the acceleration and angular rate measured by gyroscopes and accelerometers, the Strapdown Inertial Navigation System (SINS) can calculate the vehicle's velocity, position and attitude simultaneously [2]. Thus, compared with GNSS, SINS is an autonomous navigation system, which does not depend on external information, such as radio signal or radiates electromagnetic waves. Since SINS can track and reflect the vehicle's maneuvering in time, the generated navigation data have the characteristics of high accuracy in a short-term, good stability and high data update rate. However, SINS is a time integration system and the navigation error is accumulated with time rapidly because of inertial sensors' errors. Since the satellite navigation integration-based suppression method requires that individual subnetworks used for integration have large differences and the whole process takes a longer time compared with other algorithms. As a result, it is urgent to propose an effective signal denoising method to preprocess the inertial sensor data, enhancing the performance of coarse alignment.
In the process of fine alignment, the inertial sensors' error is estimated and compensated using the optimal estimation algorithm to enhance the precision of initial attitude matrix [20]. The most commonly used estimations are based on a Kalman filter (KF) that can only handle linear systems and require that the noise statistics should be accurately known [21,22]. However, in the unmanned vehicle's SINS, the performance of inertial sensors makes it difficult to determine the system and measurement noise, especially when the unmanned vehicle is maneuvering [23]. In addition, the performance of inertial sensors also leads to unsatisfactory results of the coarse alignment, resulting in the initial misalignment angle not satisfying the small-angle assumption [24]. Therefore, nonlinear and robust state estimation algorithms, such as Extended Kalman filter (EKF), Unscented Kalman filter (UKF), Cubature Kalman filter (CKF), H ∞ filter, and Huber filter, have been studied to solve the problems of error model nonlinearity, noise uncertainty, and bad external interference [6,25].
In order to enhance the initial alignment accuracy of low-precision SINS for unmanned vehicles, a novel initial alignment algorithm is proposed in this paper. In this novel initial alignment algorithm, the inertial sensor signal is preprocessed by utilizing an improved EMD denoising method based on Extreme Learning Machine (ELM) [26] and Shannon entropy to eliminate the effect of random noises on the coarse alignment firstly. Furthermore, a robust estimation method is proposed based on CKF and Huber filter to solve the influence of the nonlinearity and uncertainty on the accuracy of the fine alignment. Real-data is used to test the effectiveness of this novel initial alignment algorithm.
The rest is organized as follows. Section 2 indicates the relative background knowledge about initial alignment for unmanned vehicles. In Section 3, a novel initial alignment algorithm based on the improved EMD method and robust Huber filter is presented. The verification results are presented in Section 4 and the conclusions are drawn in Section 5.

Analytical Coarse Alignment Algorithm Based on the Solidification Coordinate Frame
The initial alignment is a key technology of the unmanned vehicle's SINS. The speed and accuracy of the initial alignment affect the system's starting speed and navigation accuracy, respectively. The initial alignment is divided into coarse alignment and fine alignment. Since there is no priori knowledge, only the measurement information from the accelerometer and the gyroscope can be used in the coarse alignment, the most common used coarse alignment method is the analytical method. In the analytical coarse alignment algorithm, a non-collinear vector is constructed firstly by utilizing the gravity vector and the earth's rotation vector measured by accelerometers and gyroscopes, respectively. Then, the strapdown attitude matrix is calculated analytically.
The traditional analytical coarse alignment requires the carrier to be static. However, SINS is inevitably subjected to various disturbance in the alignment process such as gusts, engine vibration, and maneuvering in dynamic conditions when the unmanned vehicle is moving. Thus, to extend the application and improve the accuracy of the coarse alignment, an analytical coarse alignment algorithm based on solidification coordinate frame was proposed. Its basic principle is illustrated as follows: the output of gyroscopes tracks the change of the inertial coordinate frame (i-frame), the output of accelerometers is projected on the i-frame, and, after isolating the vehicle's acceleration relative to the Earth, it can be observed that the gravity acceleration g due to the earth rotation is slowly drifting in the i-frame; the drift of g is in a cone in which the earth rotation axis is the main axis, and the geographical north can be determined from the g drift.
The solidification coordinate frame (i b0 -frame) is defined as an orthogonal reference frame nonrotating relative to the i-frame, which is formed by fixing the body coordinate frame (b-frame) at startup in the inertial space. The attitude matrix C n b is time-changing when the unmanned vehicle is moving. The attitude matrix C n b (t) at time t can be decomposed as follows: where n-frame means the navigation coordinate frame and e-frame denotes the earth coordinate frame; C n e denotes the coordinate transformation matrix between the e-frame and the n-frame; C e i is the coordinate transformation matrix between the i-frame and the e-frame; C i i b0 indicates the coordinate transformation matrix between i b0 -frame and i-frame; C n e , C e i and C i b0 b are known previously and expressed as follows: wherein L denotes the local latitude; ω ie is the Earth's rotation angular rate; t 0 denotes the starting time of initial alignment and t indicates the present time; ω b ib is the angular rate that is measured by the gyroscope directly.
Since only C i i b0 is not determined in Equation (1), the calculation of C n b can be translated into the calculation of C i i b0 . According to the principle of the double-vector attitude, C i i b0 can be determined by two non-collinear vectors whose projections in the i-frame and the i b0 -frame are already known.
The moving trajectory of the gravity vector in the i-frame is shown as Figure 1. It is known from this figure that the projection in the i-frame of the gravity vector trajectory, which, changing with time, is a cone. Since these projections are non-collinear at different time epochs, two non-collinear vectors a and b can be constructed by utilizing the gravity vector and their projections in the i-frame are expressed as a i and b i . Thus, it is obvious that Define a vector c = a × b satisfying is calculated from Equations (5)-(7) as Therefore, the matrix C n b can be calculated by Equations (1) and (8), accomplishing the coarse alignment based on the solidification coordinate frame.
There are two methods to construct the non-collinear vectors in the coarse alignment algorithm based on the solidification coordinate frame, which are velocity-based method and position-based method. In addition, in this manuscript, the position-based method is utilized: b is the actual value of C i b0 b and obtained from the following equation: whereinω b ib is the actual measured value of gyroscopes and expressed as wherein δω b ib is the measurement uncertainty of gyroscopes. In addition,Ṽ i b0 (t kj ) in Equation (9) is expressed asṼ wheref b is the actual measured value of accelerometers and expressed as wherev i b0 is the unmanned vehicle's acceleration projected to the i b0 -frame in the moving base and v i b0 is the unmanned vehicle's velocity projected to the i b0 -frame; δ f b is the measurement uncertainty of accelerometers. It is known from the previous analysis that the measurement uncertainty of gyroscopes δω b ib affects the calculation result of C i b0 b ; and the measurement uncertainty of accelerometers δ f b affects the calculation result of the gravity vector. Meanwhile, δω b ib and δ f b will lead the bias between the calculated value and the expected value of a i b0 and b i b0 , affecting the accuracy of coarse alignment.

Nonlinear Model for Integrated Fine Alignment of SINS/ GNSS Integrated Navigation Systems
In order to improve the initial alignment accuracy and speed of the unmanned vehicle's SINS, the localization information from the external GNSS is provided and the state vector of SINS/GNSS integrated navigation systems is estimated in real time by utilizing estimation methods to complete the fine alignment of the unmanned vehicle. Considering the nonlinear problem in the actual system, the nonlinear model of SINS/GNSS integrated navigation systems is established firstly.
In the literature [27,28], the nonlinear model of SINS/GNSS integrated navigation systems is derived in detail. No longer exhaustive, only the nonlinear equations are given here. The attitude error of SINS is expressed as follows based on the Euler error angle method: wherein φ = [ φ x φ y φ z ] is the Euler error angle vector, n denotes the calculation navigation coordinate system of SINS, and the direction cosine matrix from n to n is C n n ; C n b denotes the direction cosine matrix from b to n ; ε b and w b g are the gyro constant drift vector and the zero-mean Gaussian white noise vector, respectively;ω n in is the gyro measurement vector; ω n in is the rotating angular rate vector of n relative to i; δω n in is the calculated error vector of ω n in . The gyro measurement vector is equal toω n in = ω n in + δω n in . C ω is an intermediate matrix as follows: The velocity error equation is given by: whereinf b denotes the specific force vector;ω n ie andω n en are the calculated Earth's rotating angular rate and calculated angular rate, δω n ie and δω n en indicate the error vectors ofω n ie andω n en , respectively; v n and δv n denote velocity measurement vector and its corresponding error vector.
The longitude error δλ and the latitude error δϕ: wherein R M and R N are the Earth's radii of the meridian circle and the prime vertical circle, respectively; λ and ϕ are the longitude and latitude of a point of interest; v x and v y are the east and north velocities with their velocity errors δv x and δv y , respectively.
The differential equation of inertial sensors are: The state function is expressed as follows: where the state vector and the process noise are expressed as follows: where Q(t) denotes the covariance matrix of the process noise.
The difference between the localization information provided by GNSS and the localization information calculated by SINS is used as the observation and the observation function is expressed as wherein observation matrix and observation noise are expressed as follows: where R(t) is the covariance matrix of the observation noise.

Nonlinear Filter Algorithm Based on CKF
As we all know that CKF based on the spherical-radial cubature criterion is one of the most mature practical nonlinear filters. The same as the EKF, it is also based on Bayesian estimation. However, unlike the Taylor expression, which is used for linear approximation in EKF, in CKF, a set of Cubature points and corresponding weights are utilized to approximate the mean and variance of the probability distribution [29]. In this nonlinear filter, the Cubature points and weights are set by [ξ i , ω i ] firstly as follows: wherein ξ i is the i-th cubature point and ω i is the corresponding weight; i = 1, 2, . . . , 2N, and N is the dimension of the nonlinear system.

end if
Since the CKF is under the Bayesian estimation framework, we suppose that the initial and the prior density are already known accurately. Similar to EKF, the filtering operation can be divided into two parts: the time update and the measurement update [29]. The implementation of the CKF is given in detail as Algorithm 1.

Improved Denoising Method Based on the ELM and EMD-Shannon Method
It is known from the analysis results of Section 2.1 that the accuracy of the coarse alignment is affected by the measurement uncertainty of SINS inertial sensors. However, due to a vehicle's maneuvering, environment disturbance or mechanical noises in actual applications, the output signal is large nonlinearity and contains non-stationary random noises, which seriously affects the measurement accuracy of inertial sensors. In order to improve the accuracy of the coarse alignment, a denoising method based on the ELMEMD-Shannon method is proposed in this paper to preprocess the output signal of inertial sensors.

A Brief Review of the EMD Method
The EMD-based denoising method is one of the most common used signal denoising methods. The EMD-based method relies on the time scale of the signal itself to decompose adaptively multiple times and does not need the basis function. Therefore, the EMD-based method has great advantages in dealing with nonlinear and non-stationary random signals. Considered the output signal characteristics of inertial sensors, the EMD-based method is adopted to denoise the low-precision inertial measurement unit (IMU) signal in this paper.
The EMD, an effective analysis method for nonlinear and non-stationary signals, decomposes the signal into several Intrinsic Mode Functions (IMFs) and a residue adaptively based on the intrinsic characteristics of the signal. The so-called IMF is a function or signal that satisfies the following two conditions:

•
In the entire data set, the difference between the number of extreme values and the number of zero crossings must not be greater than one; • At any point of the data set, the mean value of the envelope defined by the local extrema is all zero.
For a given signal x(t), the process of EMD decomposition is illustrated as follows: • First of all, all local extreme values of the signal should be found out and identified. The cubic spline line is used to connect all the local maxima and all the local minima, producing the upper envelope and the lower envelope, respectively. Thus, all of the signal data x(t) should be covered by the upper and lower envelopes. We suppose that m 1 is the mean of the covered data by the envelopes, so the difference between the signal x(t) and the mean m 1 can be taken as a new signal, indicated as h 1 , named the first component: • In general, we can not guarantee that h 1 is a stationary data sequence, so we should repeat the above operation. Now, h 1 is taken as a new signal h 10 and its envelope mean is m 11 . Thus, the data sequence after removing the low-frequency components represented by m 11 is h 11 : Repeating the above operation up to k times, we can obtain the signal and the first IMF as follows: • Finally, the residual is a new signal that removed the high frequency component from the original signal: • Then, we can deal with the residual iteratively to get the other IMFs. The stop iterating condition is that when the residue r n becomes a monotonic function or a function with only one extremum. It means that no more IMF can be extracted from the residual signal r n . Thus, after the decomposition, x(t) is decomposed into several IMFs and a residual:

Improved EMD Denoising Method Based on ELM and Shannon Entropy
Due to the uncertainty of the basis function, the EMD method has the endpoint effect when decomposing the limited signal sequence, which will seriously affect the reliability of the EMD method. Instead of training the parameters, the ELM method uses the minimum-norm least-squares solution as the output weight of the network by solving the linear equations [30]. Thus, the ELM method has faster speed and better performance. Therefore, in order to solve the endpoint effect of the traditional EMD denoising method, the ELM method is utilized to predict and extend the inertial sensor data in this paper.
According to the basic properties of the EMD decomposition, the active ingredient concentration of signal increases with the index of IMF. The first few IMFs, especially the first IMF component, consist almost entirely of high-frequency noise. As a result, there should be a sudden change in the ratio of the high-frequency noise and the effective signal in the IMF components, resulting in a sudden change in the probability distribution of the IMF component. The effective signal can be separated by determining this mutation point.
The Shannon entropy is a method to quantify the information. Suppose that the probability distribution of a discrete variable is (p 1 , p 2 , . . . , p n ) and its Shannon entropy is defined as Suppose that x(t) is the sampled signal and it is expressed as where y(t) is the original signal and z(t) is the noise signal. The original signal has been decomposed into n IMFs and a residual by utilizing the EMD method. In order to reduce the mixing of noise, only a few IMFs and residuals are superimposed in the process of signal restoring by using the IMFs: A noise separation method based on the Shannon entropy is proposed in this paper. Suppose that the entropy of each IMF is S i (i = 1, 2, . . . , n) and the entropy variation between adjacent IMFs is expressed as: Thus, the corresponding index j s can be expressed as Therefore, an improved EMD denoising method based on the ELM method and the Shannon entropy is proposed in this paper. In this novel algorithm, the IMU signal is decomposed by using the ELMEMD method to obtain a series of IMFs firstly; then, after calculating the Shannon entropy of each IMF, the index corresponding to the maximum adjacent Shannon entropy variation is determined; on this basis, the inertial sensor signal is reconstructed so as to effectively suppress the influence on the signal quality of the inertial sensor.
The process of this novel denoising method is shown as Figure 2. The specific steps of the denoising method are as follows: Step 1: Extend the time series to the right and seven adjacent samples are used as the input of the ELM method. Use the adjacent right (or left) samples as a training sample.
Step 2: Add the previous prediction value into each new learning before each step of learning.
Repeatedly training and learning, obtain all the required extension sequence according to the required extension of the extreme points.
Step 3: Decompose the inertial sensor signal into several IMFs and residuals by using the EMD method.
Step 4: Calculate the Shannon entropy of each IMF S i .
Step 5: Calculate the adjacent Shannon entropy variation ∆S i .
Step 6: Determine the value of j s based on Step 5.
Step 7: Reconstruct the signal based on the value of j s .

Improved Robust Filter based on the RHCKF Method for Fine Alignment
The filtering performance directly affects the estimation accuracy of the system state vector. Compared with the commonly used filters, such as KF, EKF, and UKF, CKF can not only be used in nonlinear systems but also obtain better filtering accuracy. However, as we all know, the CKF is based on Bayesian estimation. When the system model is known exactly in advance and the external noise signal is a Gaussian noise, the CKF can obtain the optimal estimations; otherwise, the State estimation may be suboptimal. Because of vehicle's maneuvering and external airflow impact, the condition of optimal estimation is difficult to ensure in practical applications, resulting in filtering accuracy decreasing or even filtering divergence.
From the perspective of the approximate Bayesian estimation, the essence of the Huber filter is adding a weight matrix before the innovation, to truncate the average of filter innovations, thereby suppressing the effect of interference noise or outliers in system observation information and enhancing its robustness [22,31,32].
Assumed in a nonlinear system, the transformation innovation probability density based on the Huber cost function can be used to calculate the transformation innovation η k wherein P zz is the autocorrelation covariance matrix of the observation, z k andẑ k|k−1 are the actual observation and the predicted observation.
Calculating the transformation innovation function ϕ(η k,i ) wherein γ is the adjustment factor of the cost function in the Huber filter, and the intermediate variable ζ i can be calculated by the following equation: In Equation (39), er f (x) is the error function and: Then, the innovation weight matrix Θ k and ρ(k) can be obtained: wherein 0 ≤ ε ≤ 1 and is a standard Gaussian distribution function. The state vector and the state covariance matrix can be estimated by Equation (42), suppressing the effect of outliers and accomplishing the high-precision robust estimation: Therefore, the Huber filter is introduced into the CKF method and an improved robust Huber-CKF (RHCKF for short) algorithm is constructed in this manuscript. The flow chart of this improved algorithm is shown in Figure 3.  In this improved algorithm, the Huber filter is used to preprocess observations based on CKF framework, and a CKF measurement update is used to complete the nonlinear system estimation, avoided the linearized approximation of nonlinear equations and implemented the robustness of the algorithm.

Improved Initial Alignment Algorithm Based on ELMEMD-Shannon and RHCKF Methods
The initial alignment is a prerequisite for the normal operation of the unmanned vehicle's SINS and the alignment accuracy and convergence speed are two important performance indicators. Due to the limitation of installation space and weight in unmanned vehicles, the inertial sensor is generally smaller and low accurate, which are easily interfered by external factors, especially under dynamic conditions. In order to improve the alignment performance of unmanned vehicles in dynamic base, this paper presents a novel initial alignment algorithm. Firstly, the signal of inertial sensors is denoised by the ELMEMD-Shannon method, and, on this basis, a coarse alignment based on the solidification coordinate frame is used, suppressing the alignment error caused by dynamic noises; secondly, nonlinear system equations of the SINS/GNSS integrated fine alignment are established; meanwhile, the RHCKF filter is used to complete the state estimation, inhibited the impact of system nonlinearity and uncertainty, and ultimately improved the alignment accuracy and robustness of unmanned vehicles under dynamic conditions. The relative pseudo-code is illustrated as Algorithm 2.

Algorithm 2 Improved Initial Alignment Algorithm
Require: k = 0; Coarse alignment time T coarse ; Total alignment time T 1: if k ≥ 1 then 2: if k < T coarse then 3: Denoise the inertial sensors' information with the ELMEMD-Shannon method; 4: Coarse alignment based on solidification coordinate frame; 5: else T coarse < k < T

6:
Fine alignment with the RHCKF filter; 7: end if 8: end if 9: return Output the alignment results

Test Environment Establishment
In order to verify the performance of the denoising method proposed in this paper, the static test in the lab has been done firstly. A fiber optic gyro (FOG) IMU, developed by the Navigation Instrument Institute, Harbin Institute of Technology, was used in this test, as shown in Figure 4a. In order to isolate the interference from external environment, the IMU was installed on a metal plate and placed on a vibration isolation platform. The source data was collected by a laptop, shown as Figure 4b. The sampling frequency was 100 Hz and the test lasted 90 min. The source data was processed offline to verify the effectiveness of the proposed algorithm.
(a) (b) Figure 4. Static test in the laboratory. In this test, a FOG IMU, developed by the Navigation Instrument Institute, Harbin Institute of Technology, is used as shown in (a). This FOG IMU is placed on a stable and leveled marble platform to isolate the influence from external disturbance; In (b), a rugged laptop is used to collected the source date.

Static Test Results and Analysis
In this part, three different denoising methods, including the wavelet denoising method, traditional EMD denoising method and the improved ELMEMD-Shannon denoising method, are used to process the static data, verifying the performance of the proposed denoising method. In addition, the denoising results are shown in Figures 5 and 6. In these figures, the blue solid line indicates the original signal, the red-brown dashed line represents the signal with the wavelet denoising method, and the green dot-dashed line represents the signal with the traditional EMD denoising method while the red dotted line denotes the signal with the proposed ELMEMD-Shannon (improved EMD for short) denoising method. In order to analyze the noise reduction effect of the signal more intuitively, the result was locally amplified from 2018 s to 2019 s shown in the corresponding right subfigures of Figures 5 and 6.
It is obviously known from Figures 5 and 6 that amplitudes of inertial sensor signals are decreased with these three denoising methods while the ones of three axes are all much smaller with the improved ELMEMD-Shannon denoising method than the ones with the other two methods.
In order to quantitatively analyze the performance of the proposed denoising method, the Allan variance of gyro signals is calculated and corresponding errors, including the quantizing noise (QN), Random angle walk (RAW), Bias instability (BI), Angular rate walk (ARW) and Rate Ramp (RR), are shown in Tables 1-3.  From Tables 1-3, errors of original signals are the largest while ones with these three denoising methods are reduced to some extent and errors with ELMEMD-Shannon denoising method is the smallest. Thus, it is clear that gyroscopes' noise can be eliminated with these three denoising methods, and, moreover, the ELMEMD-Shannon denoising method has better performance.   In order to test the performance of the novel initial alignment algorithm, a dynamic test platform with a test vehicle has been established to simulate the unmanned vehicle dynamic environment, shown as Figure 7. Besides the Global Positioning System (GPS), two FOG IMUs developed by our own lab were used in this test. As shown in Figure 7, a high-precision FOG SINS was used as a reference to provide the attitude reference while a low-accuracy FOG IMU was used to be tested. In this experiment, the initial latitude is 45.7347 • N, the total time and the sampling frequency are 1 h and 100 Hz, respectively.

Test IMU
Reference IMU Figure 7. Dynamic test in a car. In this test, a high-precision FOG SINS was used as a reference to provide the attitude reference while a low-accuracy FOG IMU was used to be tested. The sampling frequency is 100 Hz and the testing time is about 1 h.

Result and Analysis
In this part, three denoising methods, the Wavelet denoising method, the EMD denoising method, and the ELMEMD-Shannon denoising method, were used to preprocess original signals firstly. In addition, the power spectrum of gyroscope and accelerometer signals is given, shown in Figures 8 and 9. From them, we can see that the power spectrum at high frequency with the ELMEMD-Shannon denoising method is smaller than the ones with other methods. This means that the high frequency noise can be eliminated with the proposed denoising method in this experiment.  In order to verify the robustness of the proposed initial alignment algorithm, different measurement noise is added to the system measurements, simulating the measurement uncertainty as follows: where R 0 is the initial measurement noise matrix.
To verify the performance of the proposed initial alignment algorithm, the initial alignment results were analyzed. Firstly, the analytical coarse alignment method based on the solidification coordinate frame is carried out by utilizing the denoised inertial information; and then two fine alignment algorithms are compared. One algorithm is the fine alignment based on the CKF method while the other is the fine alignment based on the RHCKF method, called 'Algorithm 1' and 'Algorithm 2', respectively. The time of the coarse alignment was set as two minutes and the fine alignment is carried out for the rest of the time.
The initial alignment errors are shown in Figure 10, and the error curves of the pitch, roll and yaw angle are shown in Figure 10a-c, respectively. In these subfigures, the red solid line denotes the angle error curve with the CKF fine alignment algorithm while the blue dashed line denotes the angle error curve with the RHCKF fine alignment algorithm. From Figure 10, it is obvious that the errors with Algorithm 1 is larger than the ones with Algorithm 2, and, especially, the yaw error with Algorithm 1 is more volatile. In order to further analyze the accuracy of the alignment algorithm, means and standard deviations of the errors from 1000 s to the end in each axis were calculated as shown in Figure 11. Considered the mean and standard deviation of the horizontal attitude angles, although Algorithm 1 and Algorithm 2 have little difference, the errors are smaller with Algorithm 2. Focused on the azimuth angle, the mean and standard deviation are −0.215 • and 0.042 • with Algorithm 1 while the mean and standard deviation are −0.196 • and 0.033 • , respectively. Thus, we can see that the initial alignment algorithm based on the robust RHCKF method can achieve better performance.

Conclusions
In order to enhance the accuracy of SINS initial alignment for the unmanned vehicle, a novel initial alignment algorithm is proposed in this manuscript. Considering the influence of high frequency noise from inertial sensors on the coarse alignment accuracy, a signal denoising algorithm based on the ELMEMD method and Shannon entropy for inertial sensors is proposed firstly to reduce the impact of high frequency noises. Considering the influence of the system model uncertainty caused by the maneuvering on the fine alignment accuracy, a novel fine alignment based on the robust RHCKF method to enhance the accuracy and robustness of the fine alignment. The test results showed that the proposed initial alignment algorithm is significantly superior to original algorithms.