A Hybrid Intelligent Multisensor Positioning Methodology for Reliable Vehicle Navigation

With the rapid development of intelligent transportation systems worldwide, it becomes more important to realize accurate and reliable vehicle positioning in various environments whether GPS is available or not. This paper proposes a hybrid intelligent multisensor positioning methodology fusing the information from low-cost sensors including GPS, MEMS-based strapdown inertial navigation system (SINS) and electronic compass, and velocity constraint, which can achieve a significant performance improvement over the integration scheme only including GPS and MEMS-based SINS. First, the filter model of SINS aided by multiple sensors is presented in detail and then an improved Kalman filter with sequential measurement-update processing is developed to realize the filtering fusion. Further, a least square support vector machine(LS SVM-) based intelligent module is designed and augmented with the improved KF to constitute the hybrid positioning system. In case of GPS outages, the LS SVMbased intelligent module trained recently is used to predict the position error to achieve more accurate positioning performance. Finally, the proposed hybrid positioning method is evaluated and compared with traditional methods through real field test data. The experimental results validate the feasibility and effectiveness of the proposed method.


Introduction
With the rapid development of intelligent transportation systems, global positioning system (GPS) has been widely utilized in many positioning and navigation applications for land vehicles [1,2].The major drawback related to GPS is satellite signal blockage in urban areas, high buildings, and other GPS-denied environments, where its positioning accuracy, reliability, and integrity will be deteriorated [2,3].
In this situation, GPS is usually incorporated with strapdown inertial navigation system (SINS) to achieve accurate and reliable positioning.SINS is a self-contained system that is not affected by external disturbances.However, the standalone SINS solution drifts with time.Through integrating GPS and SINS, one can achieve superior performance over a standalone GPS or SINS [4][5][6].SINS consists of a group of inertial measurement units (IMU), that is, three accelerometers and three gyroscopes, which are mounted on a moving platform to measure linear accelerations and angular velocities.Given the initial values of navigation parameters, the measurements from IMU can be processed to determine current position, velocity, and attitude of the moving platform with respect to a certain frame of reference [7,8].Note that the SINS performance depends on its grade.The high-end navigation or tactical grade SINS can maintain proper positioning performance even after a relatively long period of GPS outage.Unfortunately, the high-end SINS is very expensive and thus not suitable for low-cost applications such as those for typical passenger cars.The SINS utilizing microelectromechanical system (MEMS) technology has been developed rapidly in recent years, which provides a viable solution to achieve mass-market applications in land vehicles for its low cost, low power consumption, and small size [2,7,9].However, the MEMS-based SINS sensors may suffer from large noise, bias, and drift errors, which are much more serious than the high-grade sensors [7,9,10].Even when GPS outage occurs for a short time, the MEMS-based SINS performance will degrade rapidly.
For the SINS/GPS integration, the Kalman filter (KF) has been widely adopted as the standard optimal estimation tool.However, the limitations of the KF-based SINS/GPS integration approach are the necessity of stochastic modeling of sensor errors, which are difficult to determine especially for MEMS-based gyroscopes and accelerometers, and the need for accurate a priori information of the covariance matrices of the noises associated with both SINS and GPS [8,10,11].
Generally, although the AI-based methods above perform well for the navigation or tactical grade SINS, they have shown a very limited success when applied to a MEMS-based SINS/GPS navigation system [12].The main reason is due to the high noise level and bias instability of MEMS-based inertial sensors, which increases the nonlinear complexity of the input/output functional relationship being modeled and limits the AI generalization ability that in turn affects its prediction accuracy.To overcome the standalone MEMSbased SINS limitation when GPS outage occurs, an effective and viable solution is to fuse other sensors, especially lowcost ones in land vehicle applications [2,9].Through fusing more information from other sensors or sources, one can take advantage of the complementary properties of different sensors and thus enhance the accuracy and reliability of the positioning system.For instance, if the yaw angle information from low-cost electronic compass is fused into the system during GPS outages, the accumulative error of yaw angle caused by the integration of low-performance MEMS gyroscope will be eliminated, and then the horizontal positioning performance will be improved significantly because it is affected largely by the azimuth accuracy.
In this paper, a hybrid intelligent multisensor positioning methodology is proposed, which integrates the information from low-cost sensors including GPS, MEMS-based IMU and electronic compass, and velocity constraint.An improved Kalman filter (KF) with sequential measurement-update processing is first developed to fuse different information.Further, a least square support vector machine-(LS SVM-) based intelligent module is designed and augmented with the improved KF to form the hybrid positioning system.The hybrid positioning system has two modes of operation, that is, the update and prediction modes, which are dependent on the situation of GPS outage.The novel aspects of this paper can be summarized as follows: (1) Instead of only integrating SINS and GPS as presented in much AI-based literature, the proposed hybrid positioning methodology has the capability of fusing more information from other low-cost sensors or sources through combining the improved KF with the LS SVM-based intelligent algorithm.The distinct advantage is that, without adding much cost, the proposed methodology can significantly improve the positioning performance.
(2) The proposed methodology exhibits strong generalization ability even when the GPS outage duration exceeds the length of the sliding window designed for training the LS VSM, which has been verified by real field tests.To the authors' knowledge, the performance of the prediction stability is seldom evaluated and discussed in such cases in the existing AI-based literature.
This paper is organized as follows.Section 2 gives an overview of the proposed hybrid positioning methodology for reliable land vehicle navigation.The improved Kalman filter is presented in detail in Section 3. Section 4 describes the design of LS SVM-based intelligent module.Field tests and experimental results are provided in Section 5. Section 6 is devoted to the conclusion.

Proposed Hybrid Positioning Methodology
To achieve accurate and reliable positioning performance for land vehicles, a hybrid positioning methodology of augmenting improved KF with LS SVM is proposed, which fuses the information from low-cost sensors including GPS, MEMSbased IMU and electronic compass, and lateral velocity constraint.The augmented system comprises two modes of operation, that is, the update mode and the prediction mode.
When GPS is available, the system operates in the update (i.e., training) mode, as shown in Figure 1(a).It is noted that the Kalman filter is improved here to work with sequential measurement-update processing, which is different from the centralized processing as utilized in the common KF.The measurement update 1 is associated with the measurements of electronic compass and "virtual" lateral velocity (i.e., Ṽ  = 0), while the measurement update 2 is associated with GPS measurements.In actual implementation, two KF processing modules are executed.For the first module, simulated successive -s GPS signal outages are applied ( denotes the duration of simulated GPS outage) and the position error estimated after the measurement update 1 can be attainable.For the second module, the GPS measurements are available at all times and thus the position error estimated after two sequential measurement updates 1 and 2 is always obtained.The difference between two position errors above at each epoch is then transferred to the LS SVM-based intelligent module as its desired output during the training process.Meanwhile, the time elapsed since GPS is unavailable and the IMU output act as the inputs of LS VSM at the same epoch.Besides, considering the balance between desired accuracy and computation efficiency, a 40 s length sliding window of input/output data pattern (i.e., the window size  = 40) is chosen for the adaptive optimization of LS SVM training algorithm.

SINS output
Corrected output

Measurement
Improved KF with sequential measurement update Whenever the GPS signal is lost due to high buildings, tunnels, or other factors, the measurement update 2 cannot be executed but the measurement update 1 can still be carried out.In this case, the system operation automatically switches to the prediction mode shown in Figure 1(b), where the LS SVM-based intelligent module is utilized to predict the difference between the position error estimated after the measurement update 1 and that estimated after sequential measurement updates 1 and 2. Based on the time elapsed from the moment of GPS outage and the IMU output, the LS SVM uses the latest LS VSM parameters updated just before losing the satellite signals to predict the difference between two position errors mentioned above.Then, utilizing the available position error estimated after the measurement update 1, the position error contained in SINS can be acquired.Finally, such enhanced position error can be removed from the SINS position output to obtain more precise and reliable positioning solution for land vehicles.

Improved Kalman Filter
In this section, we describe the construction of the improved Kalman filter that is adopted for this work.

System State Equation.
According to SINS mechanization, the dynamic error model of the SINS navigation parameters (i.e., position, velocity, and attitude) can be described as [19,20] where P = [  ℎ]  is the position error vector (latitude, longitude, and height), V  = [      ]  is the velocity error vector (east, north, and up),   = [      ]  is the attitude error vector, C   is the strapdown matrix,    is the angular rate vector of the rotation of the Earth relative to the inertial frame,     is the angular rate vector of the rotation of the navigation frame relative to the Earth, f  and V  are specific force and velocity vectors in the navigation frame, respectively, f  and    are accelerometer bias and gyro drift vectors in the body frame, respectively, and D 1 and D 2 are two 3 × 3 matrices whose nonzero elements are functions of the vehicle's latitude () and height (ℎ).
In addition, considering inertial sensor residuals, the state vector of Kalman filter (KF) is composed of fifteen states, which include nine navigation parameter errors above (i.e., position, velocity, and attitude), three accelerometer biases, and three gyro drifts: where ∇  ,∇  , and ∇  denote three accelerometer constant biases in the body frame, while   ,   , and   are three gyro constant drifts in the body frame.
Based on (1) and inertial sensor residual model [19,20], the continuous-time system state equation of KF can be expressed as where G is the system noise input matrix, W is the system noise vector, and F is the system state propagation matrix.The definitions of F, G, and W are given in Appendix.
According to the method described in [20], ( 3) is discretized.The discrete-time system state equation is where  refers to the discrete-time step, X() is the discretetime state vector, Φ( + 1, ) is the discrete-time state transition matrix, Γ() is the discrete system noise input matrix, and W() is the discrete system noise with covariance matrix Q().How to determine Φ( + 1, ), Γ(), and W() is also given in Appendix.

Measurement Model.
According to Figure 1, the KF measurements consist of three parts, that is, the measurement of electronic compass, the lateral velocity constraint, and the GPS measurements.
First, the measurement model about electronic compassaided SINS can be built up according to the coordinate transformation theory [21].Here, let   represent the vehicle's yaw angle measured by electronic compass and it can be described as where   denotes the true yaw angle and    is the compass measurement Gauss white noise.On the other hand, we suppose that   represents the yaw angle calculated by SINS and   denotes the subtraction between   and   : According to the inertial navigation principle, the relationship between the ideal strapdown matrix C   and the true attitude angles is where   and   are true pitch and roll angles of the vehicle, respectively.
Correspondingly, the calculated strapdown matrix C (the symbol  denotes the SINS's calculated platform frame) can be described as where   × is the skew-symmetric form of the vector   .Expanding (9) in terms of each corresponding element, then we obtain two useful equations; that is, where C   (, ) denotes the -row and -column element of C ( = 1, 2, 3 and  = 1, 2, 3).Further, based on (8), the following can be gotten: where C   (, ) denotes the -row and -column element of C   ( = 1, 2, 3 and  = 1, 2, 3).Perturbing (11), the yaw-angle error equation can be obtained as where , and   /  = −1.
Subtracting ( 5) from ( 6) and utilizing (12), we can get Then, we will establish the measurement equation about lateral velocity constraint.It is well-known that, in general driving conditions, the vehicle does not slip or does slip very little.Thus, we can make the reasonable assumption that, for land vehicle positioning applications, the virtual measurement of lateral velocity is zero (i.e., lateral velocity constraint) [22]: where Ṽ  denotes the virtually measured lateral velocity of the vehicle in the body frame.Besides, according to the coordinate transformation theory, we have where [C   ] 2,1:3 denotes a 1 × 3 new matrix composed of the 3 elements lying in the second row of C   .V  = [      ]  is the true velocity vector of the vehicle in the navigation frame.    is the virtual measurement noise that represents the error in the assumption above.
Denoting the vehicle's velocity vector calculated by SINS in the navigation frame as V = [ V V V ]  , the following equation can be yielded: where Note that the last item of ( 17) is a second-order small quantity and can be neglected.Based on the skew-symmetric matrix features, we subtract ( 15) from ( 17) and then obtain where V  × denotes the skew-symmetric matrix formed by the vector V  .For convenient discussion later, combining ( 13) and ( 18), the measurement model about SINS aided by both electronic compass and lateral velocity constraint can be obtained: where H 1 is the corresponding observation matrix, which can be determined according to (13) and (18).V 1 is the measurement noise vector with the covariance matrix R 1 .
Finally, the measurement model about GPS-aided SINS can be easily built by subtracting the GPS measurement information from the SINS's corresponding output information.When GPS is available, it can provide accurate position and velocity information for the vehicle, that is, P  and V  shown in Figure 1(a).Thus, the measurement model of GPS-aided SINS can be expressed as where Z 2 and H 2 are GPS measurement information and observation matrix, respectively, and they are discussed in detail in [19].V 2 is the corresponding measurement noise vector with the covariance matrix R 2 .
Combining (19) and (20), the overall measurement equation of the KF can be written as where is the observation noise vector with the covariance matrix R = diag(R 1 , R 2 ).

Improved Kalman Filtering
Algorithm.For the model described by the system state equation ( 4) and the measurement equation ( 21), we can execute the filtering fusion based on the standard KF algorithm, which comprises two phases, that is, time propagation and measurement update.
(1) Time Propagation.Utilizing (4), the equation for state prediction is The covariance of the state prediction P( + 1, ) is (2) Measurement Update.The Kalman gain can be calculated as ⋅ [H ( + 1) P ( + 1, ) H  ( + 1) + R ( + 1)] −1 . ( The estimated state vector after measurement update is and the updated covariance is From ( 24)-( 26), it can be seen that traditional measurement update processes all the measurement information in a centralized mode.The measurement-update results associated with different observation vectors (e.g., Z 1 and Z here) cannot be acquired in such centralized mode, and thus it is difficult to realize the mechanization proposed in Section 2. Besides, the centralized processing mode has another limitation that it is only suitable for the situation when GPS is available.When the GPS signal is missing, the GPS measurements cannot be obtained and its corresponding measurement update should be excluded, whereas the measurement update associated with electronic compass and lateral velocity constraint should still be executed.To overcome the limitations of traditional measurement-update processing, a sequential measurement-update processing algorithm is developed here.
To realize sequential measurement-update processing for different observation sources, that is, Z 1 and Z 2 , we divide the matrices involved in the overall measurement update into two corresponding subblocks before the measurement update is executed, which is actually equivalent to scalar measurement processing discussed in [23,24].Namely, let Based on the division above, the measurement update is executed in a sequential order.That is, the measurement update about Z 1 is performed first, and then the measurement update related to GPS (i.e., Z 2 ) is performed.Specifically, although the time propagation is still executed according to ( 22)-( 23), the measurement update is sequentially executed according to the following two steps.
Note that the improved measurement-update processing described above is corresponding to the situation when GPS is available.In such situation, the measurement updates about Z 1 and Z 2 are executed sequentially according to Steps 1 and 2 above.Thus, X 1 and X 2 can easily be obtained, which are utilized for the update mode as shown in Figure 1(a).When the GPS signal is lost, only the observation information about Z 1 is executed and so only the measurement-update algorithm described by Step 1 is executed.In this case, X 1 can still be acquired, which is utilized for the prediction mode as illustrated in Figure 1(b).Therefore, the improved measurement-update processing provides an effective and feasible solution for the proposed methodology.

Design of LS SVM-Based Intelligent Module
The least square support vector machine (LS SVM) is a statistical learning theory which adopts a least squares linear system as the loss function instead of the quadratic program [25,26].LS SVM is based on the principle of structure risk minimization, which minimizes an upper bound on the generalization error, and thus its generalization performance is superior to that of the common AI methods adopting empirical risk minimization [27,28].The excellent generalization ability of LS SVM has been verified to be very effective to solve both regression and classification problems, and it has also been successfully applied for solving many engineering problems.

LS SVM Algorithm. Consider a given training set of 𝑁
data points {x  ,   }  =1 with input data x  ∈ R  and output   ∈  where R  is the -dimensional vector space and  is the one-dimensional vector space.In the feature space, the LS SVM model takes the form where  ∈ R  ℎ is an adjustable weight vector, (⋅) : R  → R  ℎ is the nonlinear function that maps the input data into a higher dimensional feature space, and  ∈  is the scalar threshold.
In the LS SVM for function estimation, the following optimization problem is formulated: where   is the error variable and  ≥ 0 is a regularization constant.Smaller  can avoid the overfitting problem in cases where data are noisy.For (31), the Lagrangian is given by where   ∈  ( = 1, . . ., ) are the Lagrange multipliers.
Then, the conditions for optimality are Eliminating  and   will yield a linear system instead of a quadratic programming problem where Y = [ 1 , . . .,   ]  , 1  = [1, . . ., 1]  , and  = [ 1 , . . .,   ]  .I  is  ×  identity matrix.Ω ∈ R × is the kernel matrix defined by The radial basis function (RBF) has been used here as the kernel function (⋅, ⋅), which is given by where ‖x−x  ‖ 2 2 is the squared Euclidean distance between the two feature vectors and  is the width of RBF.
Further, the result of the LS SVM model for function estimation becomes where   and  are the solutions to (34).The design values of  and  can be determined during the training of LS SVM.

Input/Output Space Design of LS SVM.
For the highly nonlinear input/output modeling problem discussed here, the LS SVM algorithm can provide an effective and viable solution.Besides, for land vehicle applications, the horizontal positioning performance is generally the main concern, as presented in a large body of literature [4,9,10,12,13].Therefore, to simplify the LS VSM network structure and improve computational efficiency, the latitudinal and longitudinal components of the difference between the position error estimated after the measurement update 1 and that estimated after sequential measurement updates 1 and 2 will be selected as the outputs of LS SVM; that is, where and X 1 [2] denote the first and second elements of the state vector X 1 , respectively, while X 2 [1] and X 2 [2] are the first and second elements of the state vector X 2 , respectively.As discussed above, X  is the estimated state vector after the th measurement update ( = 1, 2).Also for simplifying the LS VSM structure, only the parameters that are the main factors affecting the outputs of LS SVM will be selected as its inputs.According to the mechanism of the inertial navigation, the fundamental error sources affecting the outputs (38) mainly include the measurement errors of IMU sensors, especially the longitudinal and lateral accelerometers and the yaw gyro and the time since the GPS signal is lost.Therefore, the inputs of LS SVM here are chosen as where   ,   , and   are the longitudinal and lateral accelerations and the yaw rate, respectively, which are output by IMU. denotes the time elapsed since the GPS signal is lost.From (39), it can be found that the adopted inputs have the advantage that they are all easy to determine.

Field Tests and Results
To evaluate the performance of the hybrid intelligent multisensor positioning methodology proposed in this paper, field tests were conducted on a Buick Sail SRV vehicle.It was equipped with a low-cost NovAtel Superstar II GPS receiver with 1 Hz update rate, Crossbow MEMS-based IMU VG440CA-200 inertial sensors sampled at 100 Hz, and SDI electronic compass that provides the yaw angle measurement at 1 Hz.The low-cost IMU VG440CA-200 adopted here includes three gyroscopes and three accelerometers, which are all made based on MEMS technology.The sensor accuracies (1) are 0.05 m/s and 3 m for the GPS velocity and position when its signal is available and 0.15 ∘ for the yaw angle of electronic compass.The field tests were carried out along a representative route in an open environment, as shown in Figure 2. Since GPS signals were available all the time, the SINS/GPS/Electronic Compass/Velocity Constraint integrated solution based on KF can be obtained and utilized  as a reference for performance evaluation.During the tests, all sensor data were recorded and then the positioning methods were evaluated using the logged data.
Figure 2 shows the trajectory of this selected test route, which includes typical driving conditions such as straight motion, curved driving, and 90 ∘ turns.Throughout the whole test, a sufficient number of satellites were available.To cover different driving conditions and dynamics and assess the positioning performance comprehensively, eight simulated 40 s GPS signal outages are intentionally introduced at eight typical locations along the trajectory, which are also marked in numeric labels and shown in Figure 2. The position error discussed in the following denotes the horizontal Euclidean distance error between the estimated position and the corresponding reference, which is the main concern for land vehicle positioning.

Performance Comparison between KF and Proposed
Method.Table 1 gives a quantitative comparison of the RMS and maximum position errors between the traditional KF method and the proposed method for the simulated eight 40 s GPS outages.From Table 1, it is clear that the proposed hybrid positioning method outperforms the KF and offers much better accuracy for all GPS outages, except for outage 4. For outage 4, the maximum position error of KF is 17.97 m, which is better than that of the proposed method, that is, 21.18 m, but the difference between them is small and only about 3 m.Besides, as far as the RMS value of the position error is concerned, the proposed method is superior to the KF, which achieves 38.1% improvement in accuracy during this GPS outage.For the remaining seven outages, the proposed method achieves the significant improvement in the positioning accuracy over the KF method.The proposed method can maintain the maximum position errors within the range of 5∼22 m (the RMS range is 3∼11 m), as opposed to the range of 17∼59 m (the RMS range is 12∼26 m) achieved by the traditional KF.For the remaining seven outages, the average improvement is more than 56% whether for the RMS error or for the maximum error.
To further illustrate the improvements of the proposed method over the traditional KF method, three representative GPS outages, that is, 1, 5, and 8, are chosen to show their  trajectories, as illustrated in Figures 3-5. Figure 3 shows their positioning results for GPS outage 1, which is introduced in a 90 ∘ turn where the driving maneuvers change greatly.During this GPS outage, the maximum value of the position error of traditional KF is 28.54 m, whereas that of the proposed method is reduced to 14.77 m, which means a 48.3% performance improvement.Figure 4 gives the positioning results for GPS outage 5, which is in a typical driving scenario where the vehicle goes through a straight-curved-straight road.The solution of the proposed method follows the true trajectory more closely and provides an approximately 80% better accuracy than the traditional KF for both the RMS and maximum values of the position error.Figure 5 shows the positioning results for GPS outage 8, during which the vehicle goes through a continuously curved road.Similarly, the proposed method provides more accurate positioning outputs than the traditional KF, with about 60% performance improvement.the positioning performance for the 40 s GPS outage duration, which is equal to the designed sliding window size (i.e., 40 s).However, in practice, the period of GPS outage cannot be determined in advance.Therefore, the positioning performance should also be evaluated in conditions where the GPS outage time exceeds the size of the window for training to verify the generalization ability of the positioning system.Such evaluation is seldom discussed in the existing literature.

Performance Evaluation in
To assess the performance in longer GPS outage conditions, the size of the sliding window for training is kept unchanged but the GPS outage duration is extended from 40 s to 50 s.Besides, to compare with traditional neural network methods, the radial basis function (RBF) network as utilized in much literature is also designed.Specifically, the RBF method here denotes that the LS SVM module in Figure 1 is replaced by the RBF network with the same inputs/outputs and other modules in Figure 1 are kept unchanged.The model parameters of RBF are determined mainly via trial to achieve high training accuracy and good computational efficiency, which include 0.01 for the MSE (mean squared error) goal, 3 for the spread of radial basis functions, 40 for the maximum number of neurons, and 25 for the number of neurons to add between displays.Three typical GPS outages above, that is, 1, 5, and 8, are still chosen to be analyzed.Figure 6 shows the position errors for three extended GPS outages.Figures 7-9 illustrate the positioning results for these three extended GPS outages (i.e., 1, 5, and 8), respectively.In Figures 8-9, the locally enlarged plots are also shown.
From Figure 6, both the proposed and RBF methods can significantly enhance the positioning accuracy compared with the traditional KF when the GPS outage is less than 40 s.Moreover, they exhibit similar positioning performance within 40 s GPS outage.However, once the GPS outage is longer than 40 s, the prediction stability of RBF deteriorates, while that of the proposed method is still good and reliable, which is especially remarkable for outages 5 and 8, as shown in Figures 6(b)-6(c), 8, and 9.That is, the generalization performance of the proposed method is obviously superior to that of traditional RBF method.The main reason is that the overfitting problem easily occurs for the RBF, which cannot ensure the generalization ability of RBF once the prediction  data is beyond the space covered by the training sample.For example, for the extended GPS outage 5, the maximum value of the position error of RBF reaches 51.42 m as indicated in Figure 6(b), which is very close to that of KF, that is, 55.40 m.In contrast, the maximum value of the position error of the proposed method is only 12.36 m, which achieves an accuracy improvement of 76% and 78% over RBF and KF, respectively.Figure 8 clearly shows that the proposed method can provide more accurate positioning solution, which is much closer to the true trajectory.

Conclusions
To achieve reliable positioning performance even in adverse GPS environments, a hybrid intelligent multisensor positioning methodology is proposed for land vehicles.The proposed method can fuse the information from low-cost sensors including GPS, MEMS-based SINS, electronic compass, and lateral velocity constraint.For such information, the improved KF with sequential measurement-update processing is first developed to realize the filtering fusion.Further, the LS SVM-based intelligent module is designed and augmented with the improved KF to form the hybrid positioning system.When GPS works well, the LS SVM module is trained and updated using the recent input/output data.Once GPS outage occurs, it switches to predict the position error to achieve more accurate positioning performance.
The field tests on typical roads have been performed for evaluation of the proposed methodology.The experimental results verify that the proposed method can significantly improve the positioning accuracy over the traditional KF method for most conditions under GPS outages.In addition, it has shown much stronger prediction stability than other AI-based methods such as RBF, which indicates that it can achieve more satisfactory generalization ability.

Appendix
where all the subcomponents in F are discussed in detail in [19].For instance,

Figure 1 :
Figure 1: Proposed hybrid positioning methodology.(a) The update mode; (b) the prediction mode.
] , (A.3)where   ,   , and   denote three accelerometer white noises, while   ,   , and   are three gyro white noises.One hasΦ ( + 1, ) ≈ I + F (  ) ( +1 −   ) , (A.4)where  +1 and   are the discrete time.ConsiderΓ () ≈ G (  ) ( +1 −   ) W () = 1  +1 −   ∫  +1  W () .(A.5) , and   are the yaw, pitch, and roll angles of the vehicle calculated by SINS, respectively.The relationship between C cos   sin   sin   − sin   cos   cos   sin   cos   + sin   sin   sin   cos   sin   sin   sin   + cos   cos   sin   sin   cos   − cos   sin   − sin   cos   sin   cos   cos represents the lateral velocity calculated by SINS in the body frame and V  = [      ]  is the difference between V and V  .[C   ] 2,1:3 is a 1 × 3 new matrix composed of the 3 elements lying in the second row of C   .

Table 1 :
Statistics of position errors for KF and proposed method.
Longer GPS Outage Conditions.Note that, in Section 5.1, we only discuss and analyze = (1 − 2 + 3sin 2 ),   = (1 + sin 2 ). is the Earth semimajor axis length, and  is the Earth ellipsoid oblateness.  is the Earth rotation rate. is the latitude of the vehicle location, and ℎ is the altitude of the vehicle location.  and   are the east and north velocity components, respectively.One has