An Optimal Radial Basis Function Neural Network Enhanced Adaptive Robust Kalman Filter for GNSS/INS Integrated Systems in Complex Urban Areas

Inertial Navigation System (INS) is often combined with Global Navigation Satellite System (GNSS) to increase the positioning accuracy and continuity. In complex urban environments, GNSS/INS integrated systems suffer not only from dynamical model errors but also GNSS observation gross errors. However, it is hard to distinguish dynamical model errors from observation gross errors because the observation residuals are affected by both of them in a loosely-coupled integrated navigation system. In this research, an optimal Radial Basis Function (RBF) neural network-enhanced adaptive robust Kalman filter (KF) method is proposed to isolate and mitigate the influence of the two types of errors. In the proposed method, firstly a test statistic based on Mahalanobis distance is treated as judging index to achieve fault detection. Then, an optimal RBF neural network strategy is trained on-line by the optimality principle. The network’s output will bring benefits in recognizing the above two kinds of filtering fault and the system is able to choose a robust or adaptive Kalman filtering method autonomously. A field vehicle test in urban areas with a low-cost GNSS/INS integrated system indicates that two types of errors simulated in complex urban areas have been detected, distinguished and eliminated with the proposed scheme, success rate reached up to 92%. In particular, we also find that the novel neural network strategy can improve the overall position accuracy during GNSS signal short-term outages.


Introduction
GNSS is a commonly used technology in location-based services (LBS). A typical land vehicle navigation system (LVNS) based on GNSS has to operate in dense urban areas where GNSS signals are either blocked or severely degraded by phenomena such as cycle slips or multipath effects, which limit its capability to achieve satisfactory accuracy and positioning reliability [1]. An Inertial Navigation System (INS) can provide continuous position information and accurate attitude information during GNSS signal failure [2]. It is clear that integrating GNSS and INS can deliver an enhanced performance adaptive Kalman filtering autonomously for the purpose of adjusting the contribution of observations and dynamical models to the navigation result. Lastly, in the absence of GNSS signals, this model operates in the prediction mode to generate and supply the estimated GNSS position difference data to prevent the vehicle from leaving its path. In order to verify the effectiveness of the proposed method, a GPS/BDS real-time kinematic (RTK)/MEMS IMU integrated hardware and software have been developed, then a land vehicle test has been performed.
The structure of this paper is organized as follows: the GPS/BDS/INS integrated navigation model is clearly described in Section 2. The optimal RBF neural network aided navigation fault detection, recognition and reduction technology are discussed in Section 3. Experiments and results of the proposed method are illustrated in Section 4. Finally, Section 5 draws the conclusions and ends the paper.

GPS/BDS/INS Integrated Navigation Model
In this research, a Loosely Coupled (LC) EKF is applied for GPS/BDS/INS integration. A 24-states EKF [23,24] is used for describe the system state, which can be can be described by Equations (1): where X nav are nine navigation solution errors of three dimensional position, velocity and attitude in the north-east-down navigation frame (n-); X acc are six accelerometer error modeling parameters (bias and scale factors for each axis) in the body frame (b-), X gyro are three gyro drifts in b-, X ant are three lever arm errors in b-, respectively; X grav are three gravity uncertainty errors in the n-frame.
The discrete-time form of the dynamic model of the system has the following form: where X k and X k−1 are the state vector at epoch k and k − 1, respectively; Φ k,k−1 is the state transition matrix (see [25]) from epoch k − 1 to k; w k−1 is uncorrelated white Gaussian noise sequences. w k should satisfy the Equations (3): where E{·} denotes the expectation function. Q k is the covariance matrix of process noise. The antenna phase center in Earth Centered Earth Fixed (ECEF) (e-) frame in consideration of the deviation between the antenna phase center and the INS reference center can be written as: where C n b is the rotation matrix from b-frame to e-frame. In this paper, the difference in position between GNSS measurements and INS measurements in the n-frame is regarded as measurements, so the integrated navigation observation model can be written as: where Z k is the observation vector at epoch k, H is the measurement matrix with the following form, V k represents the measurement noise vector: w k should satisfy the following conditions: where R k is the covariance matrix of measurement noise. The state prediction and state prediction covariance [26] are given by: The measurement update step [26] is provided as: where K k is the Kalman gain matrix; the symbols "ˆ" and "~" above a variable represent an estimate and a measurement, while the superscripts "−" and "+" represent the a priori and a posteriori estimates, respectively.

Fault Detection Based on Mahalanobis Distance
Reliable GNSS observation is a prerequisite to achieve high-accuracy with a GNSS/INS integrated navigation system. In ideal conditions, the filtering observation should be Gaussian distributed, so the standard EKF is carried out according to Equations (8)- (9). The square of the Mahalanobis distance [27] from observation Z k to its meanŶ − k is treated as the relevant test statistic witch can be described as: where M k is the Mahalanobis distance,Ŷ − k is the observation prediction vector stated above withP − Y k as its associate covariance matrix which can be written as follows: If dynamical model noise w k−1 and observation error are both Gaussian distribution, this test statistic γ k should meet Chi-square distributed with degree of freedom m: where m is the degree of freedom, i.e., the dimension of the observation, χ 2 a (m) is the threshold at significance level α. α is a small value, in this contribution 1% is adopted. If the actual judging index γ k is greater than α-quantile, some kinds of observation errors can be thought to exist.

Robust Kalman Filter
In this paper, a robust EKF with a scaling factor λ k is introduced to address the observation gross errors. The observation noise covariance R k can be adapted as: Then we have: so the following equation can be satisfied: Equation (16) is nonlinear in λ k , so it can be solved iteratively using Newton's method (we omitted the derivation process): In this paper, the α-quantile χ α of the Chi-square distribution is predetermined, e.g., that for the 6-degree-of-freedom Chi-square distribution with the significance level being 1%, it is 16.812.
As we can see from Equation (8), an inflated covariance of the observation noise will result in an inflated covariance of the observation prediction. So another scalar factor, used to adjust the covariance of the observation prediction, can also be introduced to ensure the robustness as follows: (19) and κ k can be calculated analytically with: The robust Kalman gain matrix changes to: Through this method, the actual observation is less weighted and the information from the dynamic model is more weighted. The effect of the actual observation gross errors is effectively resisted.

Adaptive Kalman Filter
Even if the observation is reliable, it still requires more accurate dynamic model for both INS and GPS errors, since it is usually difficult to set a certain stochastic model for each inertial sensor that works efficiently in all environments [28]. In a low-cost GNSS/MINS integrated system, inertial sensor also includes gross error due to unmeasurable external disturbances and high dynamics which against stochastic model, and may be harmful for state prediction vector and its covariance [6]. So an adaptive Kalman filter based on state prediction covariance should be constructed to adjust the contribution of the predicted states from the IMU sensors.
A conventional adaptive Kalman filter [29] is given by: where α k is adaptive factor, α k should less than 1 when there's large dynamical model error in state prediction. The fading factor filtering established by Xia [29] has the distinct advantage that it remains convergent and tends to be optimal in the presence of model errors. Therefore we bring it into Equation (22) as an adaptive factor: where tr is the trace of a matrix: where v k = Z k − H kX − k is residual sequence. By this way, the availability factor of historical states information has been reduced, in other words, the availability factor measurement information at present time has been increased. The effect of the dynamic model errors is therefore effectively resisted.

Radial Basis Function Neural Network Algorithm
RBF neural network is a popular kernel function used in various kernelized learning algorithms. It typically has three layers: an input layer, a hidden layer with a non-linear RBF activation function and a linear output layer [30,31]. The architecture of RBF neural network is shown in Figure 1, which can be considered as a mapping R r → R s . All inputs are connected to each hidden neuron (neural unit). Then a radial basis function is used as the activation function in the hidden neuron. Usually, the Gaussian function is preferred among all possible radial basis functions due to the fact that it is factorizable. Consequently, the norm is typically taken to be the Euclidean distance and the radial basis function is commonly taken to be Gaussian. The further distance neuron's input gets from the center of radial basis function, the lower level of the activation for the neurons become. Let r p x   be the input vector, the radial basis function is provided as [32]: All inputs are connected to each hidden neuron (neural unit). Then a radial basis function is used as the activation function in the hidden neuron. Usually, the Gaussian function is preferred among all possible radial basis functions due to the fact that it is factorizable. Consequently, the norm is typically taken to be the Euclidean distance and the radial basis function is commonly taken to be Gaussian. The further distance neuron's input gets from the center of radial basis function, the lower level of the activation for the neurons become. Let x p ∈ R r be the input vector, the radial basis function is provided as [32]: where P is the total number of samples, Euclidean norm; c i is the center of the Gaussian function, σ is the variance of the Gaussian function, respectively. Then the output of the network is denoted as: where ω ij is the link weight from the hidden layer to the output layer; i = 1, 2, 3, . . . , h is the number of neurons in the hidden layer; y j is the network's actual output of the jth node corresponding to the input sample. It is assumed that d is the desired output of the sample, so the variance of the odd function is given by: The RBF learning algorithm can be expressed as follows: Step 1: The center of the basis function c i is calculated with K-means clustering.
(1) Network initialization. h training samples are randomly selected as the center of clustering c i , in other words, the center of the Gaussian function. Step 2: Computing the variance of the Gaussian function σ i . The variance of the Gaussian function σ i can be satisfied by the expression: where c max represents furthest Euclidean distance in every c i .
Step 3: Computing the link weight of the neural unit between the hidden layer and the output layer, which is given by:

Optimal RBF Neural Network Aided Navigation Fault Recognition
We take the integral of raw acceleration f b ib and rotation rates ω b ib as training inputs of RBF algorithm, which considered as the summation of angle and velocity increments [17]. Thereafter, we take GPS/BDS position difference as training outputs, which is a three dimensional vector with north, east and down position as components. The expression can be written as: where str is the RBF network structure; f b ib is the raw acceleration of the accelerometer; ω b ib is the rotation rate of the gyro; ∆P is the position increments in the n-frame; ∆B, ∆L and ∆H are the differences of geodetic coordinates, respectively; R M is the meridian radius of curvature; R N is the transverse radius of curvature [33].
As shown in Figure 2, it's a nonlinear problem from IMU outputs to position changes. The output of the network may predict the optimal estimation of state parameters of dynamic model, provided that differential navigation solution is reliable and network learning is reasonable before k epoch, so a reliable observation is a key precondition to obtain optimal estimation. For this purpose, we propose the following strategies to improve the quality of observation information to obtain optimal RBF (ORBF) training results: Step 1: GPS/BDS abnormal observations are eliminated firstly. The standard deviation of residual sequence is given by: , both the GPS/BDS observation and corresponding residual sequence should be removed, where c is a constant which may be determined by 2.0 [11].
Step 2: Optimal spread factor of RBF network is adjusted online during integrated system dynamic initialization phase.
The root mean square (RMS) of RBF predicted result can be written as: where Z GNSS is GNSS observation during training phase, Z pre is predicted value of the network, N is the number of epochs. Then the network structure is trained by the following function, which can be described as: where newrbe represents RBFNN function, s is the spread factor of the network. After the network structure is determined, according to (32), the RBFNN output can be satisfied: Optimal spread factor will make sure of getting the optimal output. From our experience, the initial value of spread factor s is set at 0.5 and the maximum value is limited to 10. It can be solved iteratively by Equations (34)-(36) with following criterion: Step 3: The RBF network training is carried out with the rest of observation and optimal spread factor by sliding window method. Step 3: The RBF network training is carried out with the rest of observation and optimal spread factor by sliding window method. In consideration of the bias instability of low-cost IMU will change over time, in this research, the width of the sliding window is set at 50. By this way, the RBF neural network can be more efficient, reliable and stable after the preprocessing above is applied on the network inputs. The predicted position with optimal RBF can avoid effect of probably fault or abnormal information at k epoch and be used to detect and locate the current observation information. That is to say, the source of system failure can be distinguished by following decision threshold and robust or adaptive Kalman filtering is automatically chosen for the integrated system: Furthermore, during satellite signals outage periods, the ORBF aided intelligent navigation system architecture provides prediction of position differences that can also replace the Kalman filter for an intelligent navigation information support. The technical scheme of this paper is shown in Figure 3. In consideration of the bias instability of low-cost IMU will change over time, in this research, the width of the sliding window is set at 50. By this way, the RBF neural network can be more efficient, reliable and stable after the preprocessing above is applied on the network inputs. The predicted position with optimal RBF can avoid effect of probably fault or abnormal information at k epoch and be used to detect and locate the current observation information. That is to say, the source of system failure can be distinguished by following decision threshold and robust or adaptive Kalman filtering is automatically chosen for the integrated system: Furthermore, during satellite signals outage periods, the ORBF aided intelligent navigation system architecture provides prediction of position differences that can also replace the Kalman filter for an intelligent navigation information support. The technical scheme of this paper is shown in Figure 3.

Field Test Details
In order to evaluate the performance of the proposed method, a field test based on a GPS/BDS RTK/MEMS-IMU integrated system was conducted around the campus of China University of Mining and Technology (CUMT). A Trimble receiver was fixed as a reference station on the campus.

Field Test Details
In order to evaluate the performance of the proposed method, a field test based on a GPS/BDS RTK/MEMS-IMU integrated system was conducted around the campus of China University of Mining and Technology (CUMT). A Trimble receiver was fixed as a reference station on the campus. During the test, the data from GPS and BDS dual constellation were used for analysis with a sampling rate of 1 Hz.
In the meantime, a GNSS/MEMS-IMU integrated navigation system with a SCC2230-E02 IMU and Unicorecomm-UB380 board as rover station with its Novatai antenna was used to perform the field test above the roof of the test vehicle. 1 Hz GPS/BDS carrier phase and pseudo-range data and 100 Hz INS raw data were received and stored in a laptop for post-processing. The hardware system configuration of the rover station is shown in Figure 4.

Field Test Details
In order to evaluate the performance of the proposed method, a field test based on a GPS/BDS RTK/MEMS-IMU integrated system was conducted around the campus of China University of Mining and Technology (CUMT). A Trimble receiver was fixed as a reference station on the campus. During the test, the data from GPS and BDS dual constellation were used for analysis with a sampling rate of 1 Hz.
In the meantime, a GNSS/MEMS-IMU integrated navigation system with a SCC2230-E02 IMU and Unicorecomm-UB380 board as rover station with its Novatai antenna was used to perform the field test above the roof of the test vehicle. 1 Hz GPS/BDS carrier phase and pseudo-range data and 100 Hz INS raw data were received and stored in a laptop for post-processing. The hardware system configuration of the rover station is shown in Figure 4.  The specifications for the low-cost MEMS IMU are given in Table 1. In addition, a magnetometer also be integrated into this system to help initialize the heading angle. The initial attitude for the IMU is given in Table 2.  The trajectory obtained with the standard Real Time Kinematic (RTK) algorithm by the modified GPStk software is given in Figure 5. The trajectory obtained with the standard Real Time Kinematic (RTK) algorithm by the modified GPStk software is given in Figure 5. A total of 3005 s GPS/BDS data were collected in this test in 10 November 2017, which started from GPS time 290,633 s and ended at 293,638 s. The observed satellites at the rover station are shown in Figure 6, which shows the visibility of GPS and BDS satellites. Figure 7 plots the number of visible satellites and the position dilution of precision (PDOP) variations for the combined GPS/BDS system in the test. The average PDOPs for the combined system can reach at the level of less than 1.5, which was evidently better than the individual system. There are no observation outliers in the GNSS observations during the period of the vehicle test. A total of 3005 s GPS/BDS data were collected in this test in 10 November 2017, which started from GPS time 290,633 s and ended at 293,638 s. The observed satellites at the rover station are shown in Figure 6, which shows the visibility of GPS and BDS satellites. Figure 7 plots the number of visible satellites and the position dilution of precision (PDOP) variations for the combined GPS/BDS system in the test. The average PDOPs for the combined system can reach at the level of less than 1.5, which was evidently better than the individual system. There are no observation outliers in the GNSS observations during the period of the vehicle test.    A loosely-coupled strategy is adopted to calculate the navigation solutions of the GPS/BDS/INS integrated system based on conventional EKF. Figure 8 illustrates the estimation of the accelerometer bias and the gyro bias for IMU sensors in integrated system, as expected, the sensor bias quickly converges to a stable value after the initial disturbances and the estimated sensor bias is consistent with the sensor specifications provided by the manufacturer.

Optimal RBF Neural Network Training
If GPS/BDS signals are available, the standard EKF strategy is adopted to calculate the navigation solutions of the GPS/BDS/INS integrated system. Simultaneously, the specific increments A loosely-coupled strategy is adopted to calculate the navigation solutions of the GPS/BDS/INS integrated system based on conventional EKF. Figure 8 illustrates the estimation of the accelerometer bias and the gyro bias for IMU sensors in integrated system, as expected, the sensor bias quickly converges to a stable value after the initial disturbances and the estimated sensor bias is consistent with the sensor specifications provided by the manufacturer.  A loosely-coupled strategy is adopted to calculate the navigation solutions of the GPS/BDS/INS integrated system based on conventional EKF. Figure 8 illustrates the estimation of the accelerometer bias and the gyro bias for IMU sensors in integrated system, as expected, the sensor bias quickly converges to a stable value after the initial disturbances and the estimated sensor bias is consistent with the sensor specifications provided by the manufacturer.

Optimal RBF Neural Network Training
If GPS/BDS signals are available, the standard EKF strategy is adopted to calculate the navigation solutions of the GPS/BDS/INS integrated system. Simultaneously, the specific increments

Optimal RBF Neural Network Training
If GPS/BDS signals are available, the standard EKF strategy is adopted to calculate the navigation solutions of the GPS/BDS/INS integrated system. Simultaneously, the specific increments of the GPS/BDS position are trained based on the ORBF neural network with corresponding acceleration and angular rate increments of the INS measurements as the input.
The feasibility of the ORBF neural network is verified using two trajectories which can be seen in Figure 9. For the authenticity of the test, one chosen route is a straight line and another one is a curve. The GPS/BDS position increments of the two routes produced by highly precise double-differenced (DD) carrier measurements are shown in Figure 10.
of the GPS/BDS position are trained based on the ORBF neural network with corresponding acceleration and angular rate increments of the INS measurements as the input.
The feasibility of the ORBF neural network is verified using two trajectories which can be seen in Figure 9. For the authenticity of the test, one chosen route is a straight line and another one is a curve. The GPS/BDS position increments of the two routes produced by highly precise double-differenced (DD) carrier measurements are shown in Figure 10. Test 1: Navigation solutions between the 1481th and 1580th seconds are provided when the vehicle moved along a straight line, 50 groups of data from the 1430th to 1480th seconds were chosen as the RBF training samples.
Test 2: Navigation solutions between the 1701th and 1780th seconds are provided when the vehicle moved along a curve, 50 groups of data from the 1650th to 1700th seconds were chosen as the RBF training samples as well.
What needs to be explained in advance is that the GPS/BDS observation has been preprocessed with the strategy mentioned in Section 3.3. In the interest of saving space, we omit discussing this step and the optimal spread factor is given directly i.e., 0.75 in our actual test. The prediction performance with ORBF network algorithm of the two tests is illustrated in Figure 11. To clarify the performance further, the prediction error of the two tests in n-frame is shown in Figure 12. The RMS error and prediction failure rate ( 3RMS  ) of the prediction errors of the two tests are summarized in Table 3. Test 1: Navigation solutions between the 1481th and 1580th seconds are provided when the vehicle moved along a straight line, 50 groups of data from the 1430th to 1480th seconds were chosen as the RBF training samples.
Test 2: Navigation solutions between the 1701th and 1780th seconds are provided when the vehicle moved along a curve, 50 groups of data from the 1650th to 1700th seconds were chosen as the RBF training samples as well.
What needs to be explained in advance is that the GPS/BDS observation has been preprocessed with the strategy mentioned in Section 3.3. In the interest of saving space, we omit discussing this step and the optimal spread factor is given directly i.e., 0.75 in our actual test. The prediction performance with ORBF network algorithm of the two tests is illustrated in Figure 11. To clarify the performance further, the prediction error of the two tests in n-frame is shown in Figure 12. The RMS error and prediction failure rate (≥ 3RMS) of the prediction errors of the two tests are summarized in Table 3.
The feasibility of the ORBF neural network is verified using two trajectories which can be seen in Figure 9. For the authenticity of the test, one chosen route is a straight line and another one is a curve. The GPS/BDS position increments of the two routes produced by highly precise double-differenced (DD) carrier measurements are shown in Figure 10. Test 1: Navigation solutions between the 1481th and 1580th seconds are provided when the vehicle moved along a straight line, 50 groups of data from the 1430th to 1480th seconds were chosen as the RBF training samples.
Test 2: Navigation solutions between the 1701th and 1780th seconds are provided when the vehicle moved along a curve, 50 groups of data from the 1650th to 1700th seconds were chosen as the RBF training samples as well.
What needs to be explained in advance is that the GPS/BDS observation has been preprocessed with the strategy mentioned in Section 3.3. In the interest of saving space, we omit discussing this step and the optimal spread factor is given directly i.e., 0.75 in our actual test. The prediction performance with ORBF network algorithm of the two tests is illustrated in Figure 11. To clarify the performance further, the prediction error of the two tests in n-frame is shown in Figure 12. The RMS error and prediction failure rate ( 3RMS  ) of the prediction errors of the two tests are summarized in Table 3.  As can be seen from the Table 3, the prediction error is significantly less than 0.5 m in test 1, and less than 1 m in the other test 2 if we take no account of several prediction gross errors. Hence, the accuracy of the predicted position in test 1 is obviously much higher than that in test 2 due to route type and smooth operation. The average RMS values are 0.099 m, 0.257 m and 0.012 m for north, east and down components respectively in test 1, 0.773 m, 0.590 m and 0.013 m for north, east and down components in test 2. The predicted anomaly is within a controllable range, which less than 7.5% in both of the test.
It is concluded that ORBF algorithm proposed in this study is able to predict GNSS position increments with sub-meter level precision, but due to biases, random walk error and stochastic noise of the low-cost IMU sensors, it will affect the high-dynamic position increments prediction. A high precision IMU may perform better owing to its low-noise characteristics.  As can be seen from the Table 3, the prediction error is significantly less than 0.5 m in test 1, and less than 1 m in the other test 2 if we take no account of several prediction gross errors. Hence, the accuracy of the predicted position in test 1 is obviously much higher than that in test 2 due to route type and smooth operation. The average RMS values are 0.099 m, 0.257 m and 0.012 m for north, east and down components respectively in test 1, 0.773 m, 0.590 m and 0.013 m for north, east and down components in test 2. The predicted anomaly is within a controllable range, which less than 7.5% in both of the test.
It is concluded that ORBF algorithm proposed in this study is able to predict GNSS position increments with sub-meter level precision, but due to biases, random walk error and stochastic noise of the low-cost IMU sensors, it will affect the high-dynamic position increments prediction. A high precision IMU may perform better owing to its low-noise characteristics.  As can be seen from the Table 3, the prediction error is significantly less than 0.5 m in test 1, and less than 1 m in the other test 2 if we take no account of several prediction gross errors. Hence, the accuracy of the predicted position in test 1 is obviously much higher than that in test 2 due to route type and smooth operation. The average RMS values are 0.099 m, 0.257 m and 0.012 m for north, east and down components respectively in test 1, 0.773 m, 0.590 m and 0.013 m for north, east and down components in test 2. The predicted anomaly is within a controllable range, which less than 7.5% in both of the test.
It is concluded that ORBF algorithm proposed in this study is able to predict GNSS position increments with sub-meter level precision, but due to biases, random walk error and stochastic noise of the low-cost IMU sensors, it will affect the high-dynamic position increments prediction. A high precision IMU may perform better owing to its low-noise characteristics.

Performance of the Proposed Method
By using the conventional EKF model, we obtained the estimated position error. The highly precise results from double-differenced (DD) carrier measurements are used only as "true values" for comparing with the results from the integrated measurements. Figure 13 plot the position error of the integrated system respect to the GPS/BDS RTK reference solution when using standard EKF model.

Performance of the Proposed Method
By using the conventional EKF model, we obtained the estimated position error. The highly precise results from double-differenced (DD) carrier measurements are used only as "true values" for comparing with the results from the integrated measurements. Figure 13 plot the position error of the integrated system respect to the GPS/BDS RTK reference solution when using standard EKF model. As can be seen from Figure 13, centimeter-level positioning accuracy is achievable through the GPS/BDS/INS integrated system after the dynamic initial alignment finished. In most cases, the accuracy is dramatically better than 3 cm. Then we simulate a specific complex environment to demonstrate the priorities of the proposed integration scheme. The complex urban areas environment can be described as: (1) 1 m, 2 m, 3 m, 4 m, 5 m, 10 m outliers, 6 in total, were intentionally given to the carrier phase measurement every other 200 epochs from 600th to 1600th epoch. (2) 2.8 g, 2.4 g, 2 g, 1.6 g, 1.2 g 0.8 g, 0.4 g outliers, 7 in total, were intentionally given to the Z-axis of the MEMS IMU measurement every other 200 epochs from 1800th to 3000th epoch.
The following five schemes are examined: The height position errors are plotted only to verify the algorithm performance. In this test, the size of the training sliding window is set to 50 epochs. From our theoretical derivation, actual test, analysis and comparison, the following conclusions can be drawn: (1) The two types of errors are obviously reflected in the results of the conventional EKF (Scheme 1, Figure 14). The filtering has no ability to resist the two types of outlier. The maximum integrated error caused by observation gross error reaches to 7.681 m, and the maximum integrated error caused by dynamic model error is −1.247 m.
(2) We recognize that the RKF (Scheme 2, Figure 15) does resist the influence of the GNSS observation gross error, but it cannot resist the influence of the dynamic model error. The filtering performs even worse than EKF when dynamic model errors occur. On the contrary, the AKF (Scheme 3, Figure 16 As can be seen from Figure 13, centimeter-level positioning accuracy is achievable through the GPS/BDS/INS integrated system after the dynamic initial alignment finished. In most cases, the accuracy is dramatically better than 3 cm. Then we simulate a specific complex environment to demonstrate the priorities of the proposed integration scheme. The complex urban areas environment can be described as: (1) 1 m, 2 m, 3 m, 4 m, 5 m, 10 m outliers, 6 in total, were intentionally given to the carrier phase measurement every other 200 epochs from 600th to 1600th epoch. (2) 2.8 g, 2.4 g, 2 g, 1.6 g, 1.2 g 0.8 g, 0.4 g outliers, 7 in total, were intentionally given to the Z-axis of the MEMS IMU measurement every other 200 epochs from 1800th to 3000th epoch.
The following five schemes are examined: • The height position errors are plotted only to verify the algorithm performance. In this test, the size of the training sliding window is set to 50 epochs. From our theoretical derivation, actual test, analysis and comparison, the following conclusions can be drawn: (1) The two types of errors are obviously reflected in the results of the conventional EKF (Scheme 1, Figure 14). The filtering has no ability to resist the two types of outlier. The maximum integrated error caused by observation gross error reaches to 7.681 m, and the maximum integrated error caused by dynamic model error is −1.247 m.
(2) We recognize that the RKF (Scheme 2, Figure 15) does resist the influence of the GNSS observation gross error, but it cannot resist the influence of the dynamic model error. The filtering performs even worse than EKF when dynamic model errors occur. On the contrary, the AKF (Scheme 3, Figure 16) can balance the contribution of updated parameters and the new measurements, but it needs the support of correctly observation. It should also be noted that AKF cannot resisted the influence of the observation gross errors.
(3) Comparing Schemes 2, 3, to the ARKF algorithm (Scheme 4, Figure 17), we find that ARKF cannot get reasonable results in the case where outliers exist. The ARKF strategy performs unsteadily, sometimes even anomalously due to the fact the two types of errors are treated as the same situation by ignoring characteristics of the errors.
(4) Among the above algorithms, the results from RRKF (Scheme 5, Figure 18) are the best. The two types of errors are effectively detected and identified by the optimal RBF neural network strategy. This algorithm can not only resist the impact of observation gross errors, but also reduce the dynamic model errors in time. This indicates that the approach of integrated navigation information fusion based on the RBF neural network aided Kalman Filter is more effective than the traditional method.
(5) However, as we can see in Figure 18, there's still an abnormal point marked in the figure even though RRKF performs much better than the other methods. It is caused by the ORBF neural network prediction gross error mentioned in Section 4.2. Speaking frankly, this error-detecting strategy has one shortcoming that it is not sensitive enough to quite small observation gross errors. So the 'success rate' of quite small observation outliers identification should be discussed in the actual test latter. measurements, but it needs the support of correctly observation. It should also be noted that AKF cannot resisted the influence of the observation gross errors.
(3) Comparing Schemes 2, 3, to the ARKF algorithm (Scheme 4, Figure 17), we find that ARKF cannot get reasonable results in the case where outliers exist. The ARKF strategy performs unsteadily, sometimes even anomalously due to the fact the two types of errors are treated as the same situation by ignoring characteristics of the errors.
(4) Among the above algorithms, the results from RRKF (Schemes 5, Figure 18) are the best. The two types of errors are effectively detected and identified by the optimal RBF neural network strategy. This algorithm can not only resist the impact of observation gross errors, but also reduce the dynamic model errors in time. This indicates that the approach of integrated navigation information fusion based on the RBF neural network aided Kalman Filter is more effective than the traditional method.
(5) However, as we can see in Figure 18, there's still an abnormal point marked in the figure even though RRKF performs much better than the other methods. It is caused by the ORBF neural network prediction gross error mentioned in Section 4.2. Speaking frankly, this error-detecting strategy has one shortcoming that it is not sensitive enough to quite small observation gross errors. So the 'success rate' of quite small observation outliers identification should be discussed in the actual test latter.   measurements, but it needs the support of correctly observation. It should also be noted that AKF cannot resisted the influence of the observation gross errors.
(3) Comparing Schemes 2, 3, to the ARKF algorithm (Scheme 4, Figure 17), we find that ARKF cannot get reasonable results in the case where outliers exist. The ARKF strategy performs unsteadily, sometimes even anomalously due to the fact the two types of errors are treated as the same situation by ignoring characteristics of the errors.
(4) Among the above algorithms, the results from RRKF (Schemes 5, Figure 18) are the best. The two types of errors are effectively detected and identified by the optimal RBF neural network strategy. This algorithm can not only resist the impact of observation gross errors, but also reduce the dynamic model errors in time. This indicates that the approach of integrated navigation information fusion based on the RBF neural network aided Kalman Filter is more effective than the traditional method.
(5) However, as we can see in Figure 18, there's still an abnormal point marked in the figure even though RRKF performs much better than the other methods. It is caused by the ORBF neural network prediction gross error mentioned in Section 4.2. Speaking frankly, this error-detecting strategy has one shortcoming that it is not sensitive enough to quite small observation gross errors. So the 'success rate' of quite small observation outliers identification should be discussed in the actual test latter.     For the purpose of evaluating the 'success rate' of the ORBF neural network strategy, 1 m dense observation gross errors, 25 in total, were intentionally added to the carrier phase measurement every other 50 epochs throughout the whole test process. As the result shown in Figure 19, the   For the purpose of evaluating the 'success rate' of the ORBF neural network strategy, 1 m dense observation gross errors, 25 in total, were intentionally added to the carrier phase measurement every other 50 epochs throughout the whole test process. As the result shown in Figure 19, the   For the purpose of evaluating the 'success rate' of the ORBF neural network strategy, 1 m dense observation gross errors, 25 in total, were intentionally added to the carrier phase measurement every other 50 epochs throughout the whole test process. As the result shown in Figure 19, the For the purpose of evaluating the 'success rate' of the ORBF neural network strategy, 1 m dense observation gross errors, 25 in total, were intentionally added to the carrier phase measurement every other 50 epochs throughout the whole test process. As the result shown in Figure 19, the reliability of the neural network proposed in this paper achieved was 92% for small observation outliers identification. reliability of the neural network proposed in this paper achieved was 92% for small observation outliers identification. Moreover, a GPS/BDS signal outages of 50 s was intentionally introduced to the raw navigation solution to demonstrate the superiorities of the ORBF-aided integrated navigation solution. Three kinds of ground track, which were generated by a standard integrated navigation solution and RBFaided integrated navigation solution are shown in Figure 20, respectively. The horizontal distance error of the INS-only estimation and ORBF-aided INS estimation are plotted in Figure 21, with the RTK fixed solution as the reference value.
It is clear that the proposed algorithm provides more accurate results in the absence of GPS/BDS signal. Without neural network aid, the INS position drift error will become larger rapidly with the increase of the outage duration and reaches 19.32 m after 50 s. With the help of the ORBF neural network algorithm, the horizontal distance error can be decreased to about 1 m during the complete GPS/BDS outage.  Moreover, a GPS/BDS signal outages of 50 s was intentionally introduced to the raw navigation solution to demonstrate the superiorities of the ORBF-aided integrated navigation solution. Three kinds of ground track, which were generated by a standard integrated navigation solution and RBF-aided integrated navigation solution are shown in Figure 20, respectively. The horizontal distance error of the INS-only estimation and ORBF-aided INS estimation are plotted in Figure 21, with the RTK fixed solution as the reference value.
It is clear that the proposed algorithm provides more accurate results in the absence of GPS/BDS signal. Without neural network aid, the INS position drift error will become larger rapidly with the increase of the outage duration and reaches 19.32 m after 50 s. With the help of the ORBF neural network algorithm, the horizontal distance error can be decreased to about 1 m during the complete GPS/BDS outage. reliability of the neural network proposed in this paper achieved was 92% for small observation outliers identification. Moreover, a GPS/BDS signal outages of 50 s was intentionally introduced to the raw navigation solution to demonstrate the superiorities of the ORBF-aided integrated navigation solution. Three kinds of ground track, which were generated by a standard integrated navigation solution and RBFaided integrated navigation solution are shown in Figure 20, respectively. The horizontal distance error of the INS-only estimation and ORBF-aided INS estimation are plotted in Figure 21, with the RTK fixed solution as the reference value.
It is clear that the proposed algorithm provides more accurate results in the absence of GPS/BDS signal. Without neural network aid, the INS position drift error will become larger rapidly with the increase of the outage duration and reaches 19.32 m after 50 s. With the help of the ORBF neural network algorithm, the horizontal distance error can be decreased to about 1 m during the complete GPS/BDS outage.