A Novel Hybrid of a Fading Filter and an Extreme Learning Machine for GPS/INS during GPS Outages

In this paper, a novel algorithm based on the combination of a fading filter (FF) and an extreme learning machine (ELM) is presented for Global Positioning System/Inertial Navigation System (GPS/INS) integrated navigation systems. In order to increase the filtering accuracy of the model, a variable fading factor fading filter based on the fading factor is proposed. It adjusts the fading factor by the ratio of the estimated covariance before and after the moment which proves to have excellent performance in our experiment. An extreme learning machine based on a Fourier orthogonal basis function is introduced that considers the deterioration of the accuracy of the navigation system during GPS outages and has a higher positioning accuracy and faster learning speed than the typical neural network learning algorithm. In the end, a simulation and real road test are performed to verify the effectiveness of this algorithm. The results show that the accuracy of the fading filter based on a variable fading factor is clearly improved, and the proposed improved ELM algorithm can provide position corrections during GPS outages more effectively than the other algorithms (ELM and the traditional radial basis function neural network).


Introduction
With the rapid development of the intelligent transportation system (ITS), vehicle navigation and positioning has attracted more and more researchers in recent years. In most systems, a combination of the Global Positioning System (GPS) and the Inertial Navigation System (INS) is used as the main positioning system. GPS receivers can provide high-precision navigation and positioning information by tracking at least four satellites. However, the performance of standalone GPS receivers may deteriorate under conditions such as GPS signal outages due to multipath effects [1,2]. However, INS is a self-contained system with excellent concealment, which makes it free from external environment interference. It obtains the position and velocity of a vehicle using the inertial measurement unit (IMU), which consists of three-axis accelerometers and three-axis rate gyros [3]. Unfortunately, errors in INS may increase over time because of noise, bias instability errors, dependent random noise, and random-work errors. Therefore, GPS is usually integrated with INS to restrain the accumulated positioning errors [4]. Currently, GPS/INS is considered to be the best style of vehicle navigation system.
GPS and INS data are processed synchronously by the fusion algorithm, which plays an important role in the accuracy of GPS/INS integrated navigation systems. As a linear optimal estimation algorithm, the Kalman filter (KF) is known as the most popular algorithm for GPS/INS integrated navigation at present. Scholars have conducted many studies based on KF to enhance the performance of combined navigation systems. So, an extended Kalman filter (EKF) algorithm was introduced in

Error Modeling
Select the "East, North, Up (ENU)" geographical coordinate system (g) as the inertial navigation system navigation reference frame, remembered as the n department. The n department, as the reference frame of the differential equation of attitude, can be expressed as where b denotes the body frame; C n b denotes the attitude matrix or direction cosine matrix, which can be used to transform from frame b to frame n; ω b nb denotes the angular rate vector of frame b with respect to the navigation frame n; (·×) denotes the skew-symmetric matrix. Since the gyroscope outputs the angular rate vector of frame b to the inertial frame n, the velocity information ω b nb cannot be directly measured. So, Equation (1) needs to be transformed, as follows [23]: where ω b ib denotes the output of the gyroscope, and ω n in is the rotation of frame n to inertial frame i, which can be expressed as ω n in = ω n ie + ω n en (3) where ω n ie denotes the earth self-rotation rate in frame n; ω n en denotes the angle rate of frame n relative to frame e in frame n; e denotes the earth frame; L and h denote the local latitude and altitude, respectively; V N and V E denote the north and east speeds of the vehicle, respectively; R M denotes the radius of curvature in the meridian; and R N denotes the radius of the curvature the prime vertical direction.
The specific force equation of INS is another basic equation, which can be expressed as . V n en = C n b f b s f − (2ω n ie + ω n en ) × V n en + g n (5) where f b s f is the vehicle ground specific force measured by the accelerometer; V n en is the vehicle velocity in frame n; and g n is the gravitational acceleration.
Last, the position information of INS can be given follows: where: where λ is the local longitude, which is provided by GPS; d is the Earth's elliptical eccentricity; and f is the flattening of the earth. By ignoring the effects of some small quantities, the linear and simple error equations of INS can be expressed as . φ = φ × ω n in + δω n in − δω n ib (8) δ .

The Model of the GPS/INS Loosely Coupled Integrated System
The navigation system based on GPS/INS can overcome the shortcomings of each navigation device and enhance the performance of the overall system. At present, there are two commonly used navigation methods: loosely coupled and tightly coupled. Considering that the loosely coupled method is easy to implement and the calculation process is simple [24], this method was chosen as the research object in this paper. The architecture of the loosely coupled GPS/INS integrated system is shown as Figure 1. The velocity difference and position difference, which are calculated by the information between GPS and INS, are the KF inputs. The INS navigation information is corrected by the output information of KF, which includes the velocity difference, position difference, gyro biases, and accelerometer biases.
where λ is the local longitude, which is provided by GPS; d is the Earth's elliptical eccentricity; and f is the flattening of the earth.
By ignoring the effects of some small quantities, the linear and simple error equations of INS can be expressed as where Equation (8)

The Model of the GPS/INS Loosely Coupled Integrated System
The navigation system based on GPS/INS can overcome the shortcomings of each navigation device and enhance the performance of the overall system. At present, there are two commonly used navigation methods: loosely coupled and tightly coupled. Considering that the loosely coupled method is easy to implement and the calculation process is simple [24], this method was chosen as the research object in this paper. The architecture of the loosely coupled GPS/INS integrated system is shown as Figure 1. The velocity difference and position difference, which are calculated by the information between GPS and INS, are the KF inputs. The INS navigation information is corrected by the output information of KF, which includes the velocity difference, position difference, gyro biases, and accelerometer biases.  According to the loosely coupled integrated system based on the Kalman filter, the state equation and measurement equation can be expressed as where F represents the state transition matrix; G represents the state noise matrix; W represents the process noise vector; Z represents the measurement vector; H represents the measurement matrix; and V represents the measurement noise matrix. A 15-state vector was proposed for the experiment in this paper. The state vector X is given by The details of the state vectors are listed in Table 1.
Velocity errors Geographical frame g δL, δλ, δh Position errors Geographical frame g ∇ x , ∇ y , ∇ z Accelerometer biases Body Frame b ε x , ε y , ε z Gyro biases Body Frame b A 6-dimension measurement vector was designed in the loosely coupled GPS/INS integrated system in this study, which includes a 3-dimension speed error and 3-dimension position error between INS and GPS. The measurement vector Z can be expressed as where V n I NS and P I NS denote the speed and location of the inertial navigation system, and V n GPS and P GPS denote the speed and location information provided by GPS.

Proposed Algorithms
In this section, the principle of the fading filter is described, and a new, improved fading factor algorithm is proposed. Then, an ELM algorithm based on a single hidden layer feed-forward neural network is discussed, while an improved ELM algorithm with a Fourier orthogonal basis function is introduced.

Principles of the Fading Filter
All of the historical measurements are utilized comprehensively in the standard Kalman filter, and the optimal estimation of the state can be obtained theoretically when the filtering model is accurate. After long-time filtering, the filtering gain calculation loop generally converges gradually and the filtering gain decreases, which makes the inertia of the filter increase, and the correction effect of the new measurement value on the state estimation decreases gradually. In order to overcome this problem, researchers proposed a fading filter algorithm to modify the system noise and the weight of measurement noise in the filtering process to gradually reduce the weight of the historical information and achieve the purpose of reducing the filter inertia. A fading filter improves the filtering accuracy under the condition of inaccurate system modeling; it is a sub-optimal filtering algorithm.
The loosely coupled integrated system model in Equation (13) is transformed into the discrete time formula: The time update isX whereX k,k−1 denotes the predicted state estimate; Φ k,k−1 denotes the state transition matrix; P k,k−1 denotes the predicted estimate covariance; S k denotes the fading factor; and S k ≥ 1. The measurement update is where K k represents the Kalman filter gain, and Q k and R k represent the covariance matrices of the state noise and measurement noise, respectively, which can be calculated as All of the above processes make up the fading filter. It can be regarded as a traditional Kalman filter when S k = 1 is satisfied in the fading filter. If S k is greater than 1, the historical information is forgotten faster in the fading filter [25].

The Improved Fading Factors
The choice of fading factor plays a key role in the performance of the fading filter. The fading factor of conventional fading filtering algorithms is usually chosen empirically. In reference [26], a recursive least squares (RLS) based variable fading factor algorithm was proposed. Reference [27] introduced a multiple fading factor calculation method. It is very difficult to apply these methods to the GPS/INS system because of the complex computation required. So, a simplified and well-behaved algorithm is proposed in this section.
The covariance matrix of the measurement prediction error σ 2 k is defined as When σ 2 k increases, it satisfies σ 2 k − σ 2 k−1 > 0; at the same time, the filter is in a divergent state. We should increase the Kalman filter gain and the system noise variance matrix to emphasize the effect of the new data on the filtering-the larger the σ 2 k value is, the greater the correction must be. According to the above description, an algorithm for calculating variable fading factor is proposed: where u represents the filter step, and β represents the ratio of estimated covariance before and after the moment. In the filtering process, β = −1 is taken if the covariance of the estimated error σ 2 k tends to increase and the ratio of the previous two times exceeds the limits, otherwise, β = 1.
On the other hand, the value of S k needs to be corrected in order to prevent the filter from diverging out of the limited range.
In summary, in the filtering process, if σ 2 k is not high, then S k will eventually approach 1, and the filter will be in a steady state. If there is a large deviation σ 2 k , S k is reduced to achieve the purpose of highlighting the new data correction. Figure 2 shows the proposed improved fading filtering algorithm (IFF).

Extreme Learning Machine
An extreme learning machine is a single hidden layer feed-forward neural network (SLFNN) algorithm, the most prominent feature of which is that it can be faster than the traditional neural network algorithm under the premise of ensuring learning accuracy. Moreover, it can randomly generate the link weights between the input layer and the hidden layer as well as the thresholds of the hidden layer neurons which need not be adjusted in the training process so that the optimal solution can be obtained by only setting the number of neurons in the hidden layer [28,29]. The architecture of the ELM is illustrated in Figure 3.  There are three layers in ELM-the input layer, the hidden layer, and the output layer from , , , , , , , an ELM with I hidden neurons can be expressed as There are three layers in ELM-the input layer, the hidden layer, and the output layer from hidden neurons can be expressed as where g(x) is the activation function, which is often used as a sigmoid function in traditional ELM and is also selected in this paper; w i = [w i1 , w i2 , · · · , w iM ] T denotes the weight between the input neurons and the hidden neurons; β i = [β i1 , β i2 , · · · , β iI ] T represents the weight connecting the hidden neurons and the output neurons; and b i is the bias of the hidden neuron. The stand SLFNN tries to minimize the difference between o j and t j , which can be expressed as At the same time, Equation (30) can also be expressed as a matrix: where H is the output of the hidden layer neuron; β is the output weight; and T is the target matrix of the N training samples.
The training target is to meet the following requirements: Thus, the output weight vector can be calculated by the smallest norm least square solution as follows: where H † is the Moore-Penrose generalized inverse of matrix H. The specific steps of the extreme learning machine are shown in Table 2. Table 2. Extreme learning machine algorithm.

Number Content
Step 1 Giving the training samples set (x i , t i ), activation function g(x), and hidden neurons I Step 2 Randomly generate hidden layer node parameters (w i , b i ) Step 3 Calculate the hidden layer output matrix H Step 4 Calculate the output weight vector β = H † T

The Improved Extreme Learning Machine
One of the biggest problems for the traditional extreme learning machine is that the activation function is fixed which leads to poor convergence of network training [30]. Therefore, an improved extreme learning machine is proposed in this section, which uses a Fourier orthogonal basis function instead of a sigmoid function for network activation function. Any nonlinear function y = f (x) can be represented linearly by a set of orthogonal basis functions: where G(x) is the orthogonal basis function; W is the correlation coefficient; and R(x) is the approximation accuracy error. According to Equation (35), the ELM mathematical model based on a Fourier orthogonal basis function can be expressed aŝ . Substituting a Fourier orthogonal basis function for a sigmoid activation function in Formula (29) gives where i denotes ith hidden neuron. The IELM has a different activation function for each neuron, which improves the training convergence rate while ensuring the training accuracy.

System Structure for GPS and INS
In the proposed system structure for GPS and INS, the IELM works in two modes. The IELM works in training mode when GPS signals are available, as shown in Figure 4. Through the specific force f b ib and angular velocity ω b ib measured from IMU, the attitude A I NS , speed V I NS , and position P I NS information of the vehicle's motion are calculated by the INS. Meanwhile, the heading H, speed V x , and speed V y are selected as the inputs of the IELM. In this system, the GPS provides the pseudo-range position and speed, which are used to loosely couple the navigation with the INS information, while the pseudo-GPS position and pseudo-GPS velocity information are collected as the inputs of the IELM. The proposed IFF algorithm is used to process the speed errors and position errors between INS and GPS. Finally, V x , V y , P E and P N are used in the final output of the system as the inputs of the IELM. The velocity error and the position error are introduced into the IELM as the target vectors of the training. Another mode is the prediction mode, which is based on the IELM in the presence of GPS outage, as shown in Figure 5. δV x , δV y , δP E , and δP N are predicted by IELM as input information of the IFF algorithm. We can obtain the final output through the INS and the IFF-processed δV, δP during GPS signals outages. , , x y E V V P , and δ N P are predicted by IELM as input information of the IFF algorithm. We can obtain the final output through the INS and the IFF-processed δ δ , V P during GPS signals outages. Figure 4. Training mode based on the improved extreme learning machine (IELM) when GPS data is available. IFF: improved fading filter. , , x y Figure 5. Prediction mode based on the IELM during GPS outages.

Simulation Test
The proposed algorithm was simulated under the GPS/INS loosely coupled mode. The biases caused by drift and the random walk noise of the accelerometer were set as μ 100 g and  Table 3.

Simulation Test
The proposed algorithm was simulated under the GPS/INS loosely coupled mode. The biases caused by drift and the random walk noise of the accelerometer were set as 100 µg and 100 µg/ √ Hz, respectively. The biases and random walk noise of the gyroscope were set as 0.02 • /h and 0.02 • / √ h, respectively. The initial misalignment angle was set as 0.01 • for heading, pitching, and roll. The GPS speed and position errors were set as 1 m and 0.1 m/s, respectively. The out frequency of the inertial sensors and GPS receiver were set as 100 Hz and 1 Hz, respectively. The vehicle movement start position was set to latitude 32.05 • N and longitude 118.79 • E. The process of the vehicle's motion is listed in Table 3. Table 3. The process of the vehicle's motion. The moving speeds of the vehicle and trajectories are shown in Figure 6a,b; the velocities of the vehicle in the east and north directions were less than 20 m/s, respectively. The whole trajectory had two turns representing rotating speeds of 9 • /s and 10 • /s, respectively. In order to verify the performance of GPS/INS integrated positioning when GPS was out of lock, two simulated GPS outages are marked by purple lines in Figure 6b, which represent GPS outages times of 50 s (350~400 s) and 100 s (450~550 s). The moving speeds of the vehicle and trajectories are shown in Figure 6a,b; the velocities of the vehicle in the east and north directions were less than 20 / m s , respectively. The whole trajectory had two turns representing rotating speeds of 9°/s and 10°/s, respectively. In order to verify the performance of GPS/INS integrated positioning when GPS was out of lock, two simulated GPS outages are marked by purple lines in Figure 6b, which represent GPS outages times of 50 s (350~400 s) and 100 s (450~550 s). First, the IFF algorithm proposed in this paper was verified. In the IFF algorithm, the select step size was set as = 0.01 u , error variance ratio threshold was = 1.3 m , and the minimum value of fading factor was = min 0.6 S . In the conventional fading filter (FF), the fading factor was selected to be = 0.9 S . The curves of the east and north position errors for the fading filter (FF) and IFF algorithms are shown in Figure 7, in which the FF and IFF are marked by black and red lines, respectively, and the simulation time is 0s to 350s. In Figure 7, we can intuitively see that the red line is less volatile than the black line, so the IFF performed better than the FF. To compare the performance of each algorithm in a clearer way, the root-mean-square errors (RMSEs) of the east and north positions for each algorithm were calculated, and they are listed in Table 4, showing that the RMSE position of the IFF algorithm was about half the RMSE position of the FF algorithm.  First, the IFF algorithm proposed in this paper was verified. In the IFF algorithm, the select step size was set as u = 0.01, error variance ratio threshold was m = 1.3, and the minimum value of fading factor was S min = 0.6. In the conventional fading filter (FF), the fading factor was selected to be S = 0.9. The curves of the east and north position errors for the fading filter (FF) and IFF algorithms are shown in Figure 7, in which the FF and IFF are marked by black and red lines, respectively, and the simulation time is 0s to 350s. In Figure 7, we can intuitively see that the red line is less volatile than the black line, so the IFF performed better than the FF. To compare the performance of each algorithm in a clearer way, the root-mean-square errors (RMSEs) of the east and north positions for each algorithm were calculated, and they are listed in Table 4, showing that the RMSE position of the IFF algorithm was about half the RMSE position of the FF algorithm.

Time (s)
The moving speeds of the vehicle and trajectories are shown in Figure 6a,b; the velocities of the vehicle in the east and north directions were less than 20 / m s , respectively. The whole trajectory had two turns representing rotating speeds of 9°/s and 10°/s, respectively. In order to verify the performance of GPS/INS integrated positioning when GPS was out of lock, two simulated GPS outages are marked by purple lines in Figure 6b, which represent GPS outages times of 50 s (350~400 s) and 100 s (450~550 s). First, the IFF algorithm proposed in this paper was verified. In the IFF algorithm, the select step size was set as = 0.01 u , error variance ratio threshold was = 1.3 m , and the minimum value of fading factor was = min 0.6 S . In the conventional fading filter (FF), the fading factor was selected to be = 0.9 S . The curves of the east and north position errors for the fading filter (FF) and IFF algorithms are shown in Figure 7, in which the FF and IFF are marked by black and red lines, respectively, and the simulation time is 0s to 350s. In Figure 7, we can intuitively see that the red line is less volatile than the black line, so the IFF performed better than the FF. To compare the performance of each algorithm in a clearer way, the root-mean-square errors (RMSEs) of the east and north positions for each algorithm were calculated, and they are listed in Table 4, showing that the RMSE position of the IFF algorithm was about half the RMSE position of the FF algorithm.   In order to verify the validity of the IELM+IFF algorithm proposed in this paper during GPS outages, two GPS outages were simulated for (#1) 350s~400 s and (#2) 450s~550 s. Figure 8a shows the east and north velocity errors for pure INS, ELM-IFF, and IELM-IFF. Meanwhile, Figure 8b displays the east and north position errors for pure INS, ELM+IFF, and IELM-IFF. Figure 8a,b clearly demonstrate that when the IFF algorithm is used at the same time, the improved ELM algorithm for the suppression of speed errors and position errors is obviously superior to the traditional ELM algorithm during GPS outages. In order to verify the validity of the IELM+IFF algorithm proposed in this paper during GPS outages, two GPS outages were simulated for (#1) 350s~400 s and (#2) 450s~550 s. Figure 8a shows the east and north velocity errors for pure INS, ELM-IFF, and IELM-IFF. Meanwhile, Figure 8b displays the east and north position errors for pure INS, ELM+IFF, and IELM-IFF. Figures 8a,b clearly demonstrate that when the IFF algorithm is used at the same time, the improved ELM algorithm for the suppression of speed errors and position errors is obviously superior to the traditional ELM algorithm during GPS outages.    Table 5. The result show that the proposed IELM-IFF algorithm was more effective than ELM-IFF as it decreased the RMSEs of the east and north velocities by about 38% and 60%, while the RMSEs of the east and north positions decreased by about 45% and 43%, respectively.

Real Road Test
To evaluate the performance of the proposed algorithm compared to the conventional counterparts in practical applications, a real road test was designed and is detailed in this subsection. The vehicle test equipment included an inertial measurement unit (IMU), a GPS receiver, PHINS, and a computer. The IMU consisted of three fiber optic gyroscopes and three accelerometers; the GPS receiver used the FlexPark6 GPS receiver, PHINS; from the French IXBLUE Inertial Navigation system; and the computer used was PC104. PHINS was used to provide accurate navigation reference information. The detailed performance parameters of IMU and GPS receiver are shown in Table 6, and the PHINS specifications are listed in Table 7. As Figure 9a shows, the PHINS and IMU were installed together, Figure 9b shows the experimental vehicle, and Figure 9c shows the structure of the vehicle test equipment. We can see from Figure 9b that the outputs of GPS provided the time-synchronization signal for PHINS and IMU. The raw data of the outputs of the IMU were transferred via an RS422 port, and the PHINS data were collected via Ethernet from the computer. The GPS data was acquired through an RS232 communication interface. At the same time, a real-time operation system, VxWorks, was embedded in the computer.  Figure 10 shows the vehicle trajectory, which was tested at the Jiulonghu campus of Southeast University in Nanjing. Figure 10a shows the Google map of the reference trajectory. Meanwhile, Figure 10b shows the coordinates of the reference trajectory.  The initial alignment time of the system was 0~300 s; after 300s, the whole system worked under the GPS/INS loosely coupled mode. The entire testing process took 1850 s and the GPS signal was good under the test environment. The yaw angle, east velocity, and north velocity information for the  Figure 10 shows the vehicle trajectory, which was tested at the Jiulonghu campus of Southeast University in Nanjing. Figure 10a shows the Google map of the reference trajectory. Meanwhile, Figure 10b shows the coordinates of the reference trajectory.  Figure 10 shows the vehicle trajectory, which was tested at the Jiulonghu campus of Southeast University in Nanjing. Figure 10a shows the Google map of the reference trajectory. Meanwhile, Figure 10b shows the coordinates of the reference trajectory.  The initial alignment time of the system was 0~300 s; after 300s, the whole system worked under the GPS/INS loosely coupled mode. The entire testing process took 1850 s and the GPS signal was good under the test environment. The yaw angle, east velocity, and north velocity information for the  The initial alignment time of the system was 0~300 s; after 300s, the whole system worked under the GPS/INS loosely coupled mode. The entire testing process took 1850 s and the GPS signal was good under the test environment. The yaw angle, east velocity, and north velocity information for the entire exercise are shown in Figure 11, where the yaw angle is depicted by the blue line, and the east velocity and the north velocity are described by the black and red lines, respectively.  Figure 11, where the yaw angle is depicted by the blue line, and the east velocity and the north velocity are described by the black and red lines, respectively. The proposed algorithm was verified by the data collected in the above real road test. First, the IFF algorithm and FF algorithm were verified by the east position error and the north position error of GPS/INS integrated navigation from 350 s to 650 s, which is shown in Figure 12. In the IFF algorithm, the step size was = 0.005 u , the error variance ratio threshold was 1.3 m = , and the minimum value of the fading factor was = min 0.6 S , while in FF, the fading factor was = 0.9 S . In order to more accurately illustrate the superiority of the IFF algorithm, the RMSEs of the east and north positions were obtained by calculating the root-mean-square error of the error data in Figure 12, which are listed in Table 8. Compared with the FF algorithm, the RMSEs of the east position and north position of the IFF algorithm (0.6319 and 0.9639) reduced by about 38% and 15%, respectively.   The proposed algorithm was verified by the data collected in the above real road test. First, the IFF algorithm and FF algorithm were verified by the east position error and the north position error of GPS/INS integrated navigation from 350 s to 650 s, which is shown in Figure 12. In the IFF algorithm, the step size was u = 0.005, the error variance ratio threshold was m = 1.3, and the minimum value of the fading factor was S min = 0.6, while in FF, the fading factor was S = 0.9. In order to more accurately illustrate the superiority of the IFF algorithm, the RMSEs of the east and north positions were obtained by calculating the root-mean-square error of the error data in Figure 12, which are listed in Table 8. Compared with the FF algorithm, the RMSEs of the east position and north position of the IFF algorithm (0.6319 and 0.9639) reduced by about 38% and 15%, respectively.  Figure 11, where the yaw angle is depicted by the blue line, and the east velocity and the north velocity are described by the black and red lines, respectively. The proposed algorithm was verified by the data collected in the above real road test. First, the IFF algorithm and FF algorithm were verified by the east position error and the north position error of GPS/INS integrated navigation from 350 s to 650 s, which is shown in Figure 12. In the IFF algorithm, the step size was = 0.005 u , the error variance ratio threshold was 1.3 m = , and the minimum value of the fading factor was = min 0.6 S , while in FF, the fading factor was = 0.9 S . In order to more accurately illustrate the superiority of the IFF algorithm, the RMSEs of the east and north positions were obtained by calculating the root-mean-square error of the error data in Figure 12, which are listed in Table 8. Compared with the FF algorithm, the RMSEs of the east position and north position of the IFF algorithm (0.6319 and 0.9639) reduced by about 38% and 15%, respectively.  Second, to validate the performance of the IELM algorithm, the results were compared with the ELM and RBF neural network algorithms. The number of hidden layer nodes in the IELM, ELM, and RBF neural network was all 20. The training time was 300 s to 650 s, and the simulation set the  Second, to validate the performance of the IELM algorithm, the results were compared with the ELM and RBF neural network algorithms. The number of hidden layer nodes in the IELM, ELM, and RBF neural network was all 20. The training time was 300 s to 650 s, and the simulation set the GPS outage time from 650 to 750 s. The east velocity error and north velocity error for the pure INS, RBF-IFF, ELM-IFF, and IELM-IFF of GPS outages are shown in Figure 13a Figure 14a,b. Compared with other algorithms, the IELM-IFF algorithm had the smallest position error when the GPS signal lost lock. To compare the performance of each algorithm in a clearer way, the root-mean-square error and maximum error of the velocity and position information for each algorithm during the GPS outages are listed in Table 9. The maximum errors and RMSEs of the velocity and position for IELM-IFF were the smallest compared with the other algorithms. To compare the performance of each algorithm in a clearer way, the root-mean-square error and maximum error of the velocity and position information for each algorithm during the GPS outages are listed in Table 9. The maximum errors and RMSEs of the velocity and position for IELM-IFF were the smallest compared with the other algorithms.   Another existing algorithm was compared to demonstrate the effectiveness of the proposed approach. The artificial-intelligence-based segmented forward predictor (ASFP) was developed in reference [31] and uses two RBFNNs and a forward prediction algorithm. The ASFP algorithm was compared with the proposed IELM-IFF algorithm. Figure 15   Another existing algorithm was compared to demonstrate the effectiveness of the proposed approach. The artificial-intelligence-based segmented forward predictor (ASFP) was developed in reference [31] and uses two RBFNNs and a forward prediction algorithm. The ASFP algorithm was compared with the proposed IELM-IFF algorithm. Figure 15 presents the results during GPS outages from 650 to 750 s, where Figure 15a Third, Figure 16 shows the position prediction errors for ELM-IFF and IELM-IFF in the case of good GPS signal with training times from 300 to 525 s and a forecast period of 525 to 570 s. In addition, the RMSEs of the position for these two algorithms under good GPS signal are listed in Table 10. The RMSEs of the east position and north position predictions were 0.7346 and 1.8919, which is better than the results for the ELM-IFF algorithm.   [32]. Figure 17 displays the horizontal velocity error and the horizontal position Third, Figure 16 shows the position prediction errors for ELM-IFF and IELM-IFF in the case of good GPS signal with training times from 300 to 525 s and a forecast period of 525 to 570 s. In addition, the RMSEs of the position for these two algorithms under good GPS signal are listed in Table 10. The RMSEs of the east position and north position predictions were 0.7346 and 1.8919, which is better than the results for the ELM-IFF algorithm. Third, Figure 16 shows the position prediction errors for ELM-IFF and IELM-IFF in the case of good GPS signal with training times from 300 to 525 s and a forecast period of 525 to 570 s. In addition, the RMSEs of the position for these two algorithms under good GPS signal are listed in Table 10. The RMSEs of the east position and north position predictions were 0.7346 and 1.8919, which is better than the results for the ELM-IFF algorithm.   [32]. Figure 17 displays the horizontal velocity error and the horizontal position  Fourth, to evaluate the performance of the proposed model, the O I NS − X k model was utilized as a comparison [32]. Figure 17 displays the horizontal velocity error and the horizontal position error, where Figure 17a,b represent good GPS signals, and Figure 17c,d represent GPS outages from 650 to 750 s. When the GPS signal was good, the proposed model showed better accuracy than the O I NS − X k model in most cases. When the GPS signal was unavailable for 650 to 750 s, the two models showed different results. It is obvious that the proposed model outperformed the O I NS − X k model. During the GPS outages for a period of 650 to 680 s, the velocity errors of the O I NS − X k model and proposed model were similar. This means that during short GPS outages, both models could be employed to reduce the velocity error. However, when the GPS outage becomes longer, the proposed model achieves higher accuracy than the O I NS − X k model.
The reason why the proposed model performs better than the O I NS − X k model is that the predicted pseudo GPS position only relates to the output of INS, while the predicted state vector X k is influenced by both the INS information and the accuracy of the loosely-coupled KF. When the last estimation of the KF is correct, the O I NS − X k can be well utilized. However, the GPS/INS integrated system cannot ensure absolute accuracy, which means there are always small errors in the estimation of the KF.  Fifth, in order to compare the computational costs between different algorithms, RBF, ELM-IFF, ASFP, and IELM-IFF with the same number of nodes were investigated. To facilitate the observation of the computational cost, the number of pieces of training sample data was set as 100. The average time consumption of training procedures for each algorithm in the simulation is listed in Table 11, showing that the ELM-IFF, ASFP, and IELM-IFF performed faster than RBF. The computational cost of the proposed IELM-IFF algorithm was 4.71ms, which is almost similar to the Fifth, in order to compare the computational costs between different algorithms, RBF, ELM-IFF, ASFP, and IELM-IFF with the same number of nodes were investigated. To facilitate the observation of the computational cost, the number of pieces of training sample data was set as 100. The average time consumption of training procedures for each algorithm in the simulation is listed in Table 11, showing that the ELM-IFF, ASFP, and IELM-IFF performed faster than RBF. The computational cost of the proposed IELM-IFF algorithm was 4.71ms, which is almost similar to the ELM-IFF algorithm. So, the IELM-IFF algorithm performed better than ELM-IFF and, at the same time, it did not increase the computational cost. In addition, the IELM-IFF method was shown to have a lower computational coast than the ASFP algorithm proposed in reference [31]. Finally, Figure 18 shows the curves of convergence performance for the ELM and IELM algorithms, in which the ordinate represents the RMSE values during training, and the abscissa indicates the training times. It can be seen from Figure 16 that the IELM algorithm achieved higher convergence accuracy than the ELM algorithm for the same training time. Finally, Figure 18 shows the curves of convergence performance for the ELM and IELM algorithms, in which the ordinate represents the RMSE values during training, and the abscissa indicates the training times. It can be seen from Figure 16 that the IELM algorithm achieved higher convergence accuracy than the ELM algorithm for the same training time.

Conclusions
In this paper, a new algorithm was proposed for GPS/INS integrated navigation during GPS outages based on a fading filter and an extreme learning machine, and a new training model strategy was developed. An improved fading filter algorithm was also proposed with the aim of adjusting the fading factor in the traditional forgetting filter. This algorithm can dynamically adjust the fading factor so that the fading filter can achieve a better filtering effect in real-time. In order to solve the problem of the rapid divergence of a GPS/INS loosely coupled navigation system during GPS outages, this paper introduced the ELM algorithm into the integrated navigation system, which greatly improved the speed of the network training compared with the traditional radical basis function (RBF) neural network. In order to solve the problems of the fixed activation function in the ELM algorithm and the slow convergence speed, this paper presented an improved ELM algorithm based on a Fourier orthogonal basis function.
In order to verify the performance of the proposed algorithm, this paper proposed simulation experiments and a real road vehicle test. Regarding the IFF algorithm verification, by comparing and analyzing the FF algorithm under good GPS signal conditions, the IFF algorithm was shown to have a better filtering effect. The performance of the IELM algorithm was verified by the training time, prediction accuracy, and convergence speed, and compared with a traditional RBF neural network. The results showed that compared with the RBF and ELM algorithms, IELM can reduce the divergence of inertial navigation errors and achieve higher positioning accuracy.

Conclusions
In this paper, a new algorithm was proposed for GPS/INS integrated navigation during GPS outages based on a fading filter and an extreme learning machine, and a new training model strategy was developed. An improved fading filter algorithm was also proposed with the aim of adjusting the fading factor in the traditional forgetting filter. This algorithm can dynamically adjust the fading factor so that the fading filter can achieve a better filtering effect in real-time. In order to solve the problem of the rapid divergence of a GPS/INS loosely coupled navigation system during GPS outages, this paper introduced the ELM algorithm into the integrated navigation system, which greatly improved the speed of the network training compared with the traditional radical basis function (RBF) neural network. In order to solve the problems of the fixed activation function in the ELM algorithm and the slow convergence speed, this paper presented an improved ELM algorithm based on a Fourier orthogonal basis function.
In order to verify the performance of the proposed algorithm, this paper proposed simulation experiments and a real road vehicle test. Regarding the IFF algorithm verification, by comparing and analyzing the FF algorithm under good GPS signal conditions, the IFF algorithm was shown to have a better filtering effect. The performance of the IELM algorithm was verified by the training time, prediction accuracy, and convergence speed, and compared with a traditional RBF neural network. The results showed that compared with the RBF and ELM algorithms, IELM can reduce the divergence of inertial navigation errors and achieve higher positioning accuracy.