Next Article in Journal
JMAC Protocol: A Cross-Layer Multi-Hop Protocol for LoRa
Previous Article in Journal
Application of Correlation Analysis for Assessment of Infrasound Signals Emission by Wind Turbines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive Multi-Dimensional Vehicle Driving State Observer Based on Modified Sage–Husa UKF Algorithm

1
Chongqing University-University of Cincinnati Joint Co-op Institute, Chongqing University, Chongqing 400030, China
2
State Key Laboratory of Power Transmission Equipment & System Security and New Technology, Chongqing University, Chongqing 400030, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(23), 6889; https://doi.org/10.3390/s20236889
Submission received: 30 October 2020 / Revised: 21 November 2020 / Accepted: 30 November 2020 / Published: 2 December 2020
(This article belongs to the Section Intelligent Sensors)

Abstract

:
An accurate vehicle driving state observer is a necessary condition for a safe automotive electronic control system. Vehicle driving state observer is challenged by unknown measurement noise and transient disturbances caused by complex working conditions and sensor failure. For the classical adaptive unscented Kalman filter (AUKF) algorithm, transient disturbances will cause the failure of state estimation and affect the subsequent process. This paper proposes an AUKF based on a modified Sage–Husa filter and divergence calculation technique for multi-dimensional vehicle driving state observation. Based on the seven-degrees-of-freedom vehicle model and the Dugoff tire model, the proposed algorithm corrects the measurement noise by using modified Sage–Husa maximum posteriori. To reduce the influence of transient disturbance on the subsequent process, covariance matrix is updated after divergence is detected. The effectiveness of the algorithm is tested on the double lane change and Sine Wave road conditions. The robustness of the algorithm is tested under severe transient disturbance. The results demonstrate that the modified Sage–Husa UKF algorithm can accurately detect transient disturbance and effectively reduce the resulted accumulated error. Compared to classical AUKF, our algorithm significantly improves the accuracy and robustness of vehicle driving state estimation. The research in this paper provides a reference for multi-dimensional data processing under changeable vehicle driving states.

1. Introduction

According to the released Global Status Report on Road Safety 2018 [1] by WHO, the annual road traffic deaths reached 1.35 million in 2016. An effective solution towards safe driving is a worldwide problem. The stability-control technique is one of the most important techniques in vehicle active security [2].
Since the 1990s, automotive electronics have developed rapidly. Many methods have been proposed to improve automotive safety, including anti-lock brake system (ABS) based on tire dynamic model [3], traction control system (TCS) based on slip ratio [4], electric stability program (ESP) based on measured rolling angle speed, acceleration, and yaw moment of a vehicle [5]. An accurate vehicle driving state is a precondition for these systems.
However, vehicle driving state cannot be obtained by vehicular sensors and requires additional sensors because of complex vehicle driving conditions, low accuracy, and high cost of vehicular sensors. Even if enough high-accuracy sensors are installed regardless of the cost, the changing accuracy and transient disturbance cannot be avoided in vehicle driving. Therefore, the accuracy and robustness of the vehicle driving state observer is a key problem in the study of the automotive electronic control system. It is also the prerequisite and essential condition of closed-loop control.
The observation of vehicle driving states refers to the process of using the data measured by sensors to estimate other parameters of the vehicle system. The estimated parameters are usually difficult or inaccurate to measure. Most studies in this field focus on obtaining more accurate estimation based on collecting multi-dimensional information by different information fusion technologies.
There are two kinds of commonly used observation methods. The first one is based on dynamics, including Kalman filter [6,7,8] and neural network [9,10]. As for the neural network, the internal mapping relation is closely related to training data. It is unrealistic to apply all working conditions training on neural networks. The second estimation method is based on kinematics [11]. Sensors are used to measure vehicle data and estimate target parameters by following kinematic relationships. This method requires sensors with high accuracy and the error will accumulate in the whole system, which is with low robustness.
When the target state is difficult to be directly detected by a sensor, there are two common options. Firstly, we can design an optimized sensing method and indirectly detect the state. Secondly, we can utilize some measurement-friendly states to estimate the target state. Automobile scientists keep working on using the second method to estimate vehicle driving states. Instead of hardware measurement, a proper and accurate estimation can not only save costs but also adaptive to different vehicles. Scientists are optimistic about the application of the Kalman filter because it is a process of estimation and correction and does not need to store or analyze raw observation data. The development of the Kalman filter and verification of real vehicles proves this idea.
In 1995, L.R.Ray et al. [12] developed an extended Kalman filter (EKF) to estimate non-linear vehicle driving state based on wheel angular velocity, longitudinal and lateral acceleration, and yaw rate. Although this method was not widely spread because of the immaturity, it promoted the usage of Kalman filter in vehicle driving state estimation.
Vehicle states change violently and non-linearly in normal running. Based on the unscented Kalman filter (UKF), Doumiati M et al. [13] designed a state estimator to solve problems when applying KF and EKF in nonlinear systems. They accurately estimated tire forces and sideslip angles. Their method proved the feasibility of using UKF soft sensing to replace previous high-cost measurements. However, they did not consider the driving states under realistic road conditions. Shengbo Eben Li et al. [14] uses linearly arrayed ultrasonic sensors to track moving objects around the vehicle, which focuses on the perception of the surrounding environment.
The application of GPS/IMU Fusion in vehicle state estimation improves the accuracy but relies on GPS signals and with high cost. When GPS signals face severe transient disturbances or loss, the filter will diverge [15], which is with low robustness.
Vincent Havyarimana et al. [16] proposed a hybrid approach combining square root unscented Kalman filter and particle filter to determine vehicle states. The approach considered noise as a limited gaussian kernel mixture and estimated states by sparse gaussian kernel model. It solves the problem of filter diverging when a transient disturbance or weak GPS signal occurs. However, it is still at a high cost due to the usage of the GPS/INS component.
Sage–Husa adaptive filtering is applied in the unscented transformation process. It can optimize the mean and the error covariance of the measurement and process noise. Sage–Husa adaptive filter is widely applied in target tracking. Based on SINS (Strapdown Inertial Navigation System), CNS (Celestial Navigation System), and GNSS (Global Navigation Satellite System), Shuqing Xu et al. [17] proposed a method to track the target in a multi-system and multi-model way. It focused on filter fusion and significantly improved tracking accuracy. To solve the undetermined measurement noise and system noise, a classical Sage–Husa adaptive filter was involved. Combining the traditional UKF and traditional Sage–Husa filter, Bin Huang et al. [18] proposed a method based on a four-degrees-of-freedom nonlinear vehicle dynamic model and a simplified Magic Formula tire model. Their method could effectively detect the change in vehicle state errors. To estimate process noise, Qingya Wu et al. [19] applied Sage–Husa sub-optimal maximum posteriori noise estimator to adaptive simplified spherical simplex unscented Kalman filter. Their approach was more efficient than AUKF and effectively improved computational accuracy. To further improve the Sage–Husa algorithm, Gang Li et al. [20] eliminated the estimation of process noise and the error correction term of measurement noise estimation. However, their algorithm was applied in an optimized system and did not consider transient disturbances.
In practical application, vehicle driving states are affected by many factors including working conditions and weather. The prior variables of the measurement noise are unknown. Although UKF and AUKF can track the changing state, the estimation results will lag true values because of the system’s internal unchangeable process noise [21]. When working conditions change drastically, the measurement information transmitted to the observer is inaccurate. The disturbance caused by the changing working condition will cause the failure of state estimation and affect the subsequent process. When prior noise statistics are unknown or inaccurate, the Sage–Husa estimator and Sage–Husa maximum posterior have shown their potentialities in estimating the measurement and process noise in various fields [22,23,24,25]. However, the Sage–Husa estimator and Sage–Husa maximum posterior will diverge after being disturbed.
Therefore, we propose an AUKF based on a modified Sage–Husa filter and divergence calculation technique for multi-dimensional vehicle driving state observation. The noise estimator is based on a modified Sage–Husa maximum posteriori estimation. Its purpose is to perform an online estimation on the mean and error covariance of measurement noise. Besides, we use the divergence calculation technique to modify the covariance matrix when a disturbance occurs. To evaluate the practical value of our method, we add a severe disturbance of 0.1 s in the robustness test. The result shows that the modified Sage–Husa UKF algorithm can effectively reduce the impact of disturbances.
The main contributions of our work are shown as follows:
  • We propose a modified Sage–Husa UKF algorithm. It effectively increases the filtering convergence speed and reduces the influence of transient disturbance on subsequent estimation.
  • We introduce divergence calculation to detect transient disturbances and solve the divergence problem of Sage–Husa maximum posteriori. Meanwhile, the introduction of divergence calculation technique reduces the amount of calculation.
  • We apply the modified Sage–Husa UKF algorithm on vehicle driving state estimation. The simulation result proves the robustness of our approach, which provides a reference for multi-dimensional data processing under changeable vehicle driving states.
As for the rest of this paper, Section 2 introduces the seven-degrees-of-freedom vehicle model and the Dugoff tire model. Section 3 introduces the proposed approach of vehicle driving state estimation based on the Modified Sage–Husa UKF algorithm. Section 4 introduces the simulation and robustness test results of our proposed approach. Section 5 concludes the paper.

2. Models

Our proposed method is based on the nonlinear seven-degrees-of-freedom dynamic vehicle model and the Dugoff tire model. The seven-degrees-of-freedom vehicle model contains all the variables we are concerned with. It was established to connect the measured variables from sensors with the obtained state variables. However, some new state variables were involved in the vehicle model, for example, the longitudinal and lateral tire forces. Since these variables could not be measured directly, the Dugoff tire model was involved to calculate and obtain these variables.

2.1. Seven-Degrees-of-Freedom Vehicle Model

In this paper, we apply the estimation based on the seven-degrees-of-freedom vehicle model. The established seven-degrees-of-freedom dynamic vehicle model contains three longitudinal, lateral, and yaw degrees of freedom (DOFs), and four rotational DOFs of the wheels [26]. We can obtain longitudinal and lateral motion including speed and acceleration, yaw rate, side slip angle, and torque from this model. In Figure 1, the ISO coordinate system is set with the mass center of a vehicle as the origin, where the right part of the x-axis and the upper part of the y-axis is positive. The vehicle is longitudinal symmetry about the x-axis. The set of dynamic equations obtained from the model is shown as follows:
{ u ˙ = a x + v ω v ˙ = a y + v ω ω ˙ = Γ I z Γ = a ( F x 1 sin δ 1 + F y 1 cos δ 1 ) t f 2 ( F x 1 cos δ 1 F y 1 sin δ 1 ) + a ( F x 2 sin δ 2 + F y 2 cos δ 2 ) + t f 2 ( F x 2 cos δ 2 F y 2 sin δ 2 ) b ( F x 3 sin δ 3 + F y 3 cos δ 3 ) t f 2 ( F x 3 cos δ 3 F y 3 sin δ 3 ) b ( F x 4 sin δ 4 + F y 4 cos δ 4 ) + t f 2 ( F x 4 cos δ 4 F y 4 sin δ 4 ) + β = v u a x = 1 m ( F x 1 cos δ 1 F y 1 sin δ 1 + F x 2 cos δ 2 F y 2 sin δ 2 + F x 3 cos δ 3 F y 3 sin δ 3 + F x 4 cos δ 4 F y 4 sin δ 4 F w ) a y = 1 m ( F x 1 sin δ 1 + F y 1 cos δ 1 + F x 2 sin δ 2 + F y 2 cos δ 2 + F x 3 sin δ 3 + F y 3 cos δ 3 + F x 4 sin δ 4 + F y 4 cos δ 4 ) ,
where m is the total mass of the vehicle; β is side slip angle; u is longitudinal speed; v is lateral speed; ω is yaw rate; a x is longitudinal acceleration; a y is lateral acceleration; Γ is torque around the z-axis; I z is car inertia around z-axis; δ is wheel steering angle; F x i is tire longitudinal force; F y i is tire lateral force; i = 1 ,   2 ,   3 ,   4 represents the wheel at the front left, front right, rear left and rear right; δ i represents the steering angle of the specific wheel responding to the wheel steer angle   δ ; F w is wind resistance; a , b represent the distance from the mass center to the front axis and the back axis; t f , t r represent the distance of wheels at the front axis and back axis.

2.2. Dugoff Tire Model

Pacejka Magic Formula [27], Dugoff model [28], and Burckhardt model [29] is mostly used in these years to predict the behavior of tires. Comparing with the other two models, the Dugoff model is more suitable for UKF, and UKF can perfectly track the road adhesion coefficient. In this paper, we use the Dugoff tire model to simulate tire longitudinal ( F x ) and lateral force ( F y ). The schematic is shown in Figure 2. In this tire model, the tire longitudinal and lateral forces are closely related to parameters including longitudinal slip ratio λ , tire longitudinal ( C x ) and lateral stiffness ( C y ).
The tire longitudinal and lateral force can be represented by Equations (2) and (3).
F x = C x λ 1 λ f ( L )
F y = C y tan α 1 λ f ( L ) ,
where
f ( L ) = { L ( 2 L ) , L < 1 1                             , L 1
L = μ F z C x 2 λ 2 + C y 2 tan 2 α   ( 1 λ ) ( 1 ε u C x 2 λ 2 + C y 2 tan 2 α ) ,
and μ is the current road adhesion coefficient; F z is the vertical load of each tire; L  is the introduced boundary condition; ε is the speed influence factor. We obtain a road adhesion coefficient based on [30].
State variable x ( t ) is represented as x ( t ) = [ v x v y r Γ a x a y ] T . Observation variable y ( t ) is represented as y ( t ) = [ a x a y r ] T . Control input u ( t ) = [ δ ] is represented as. Based on the seven-degrees-of-freedom vehicle model and Dugoff tire model, state equation can be written as:
X k + 1 = f ( X k , u k ) + Λ W k = [ v x v y r 0 0 0 ] k 1 + [ a x + v y r a y v x r Γ / I z 0 0 0 ] k 1 · t + [ 0 0 0 Γ a x a y ] k + Λ W k .
System observation equation can be written as:
Y k + 1 = h ( X k , u k ) + V k = [ 0 0 0 0 1 0 0 0 0 0 0 1 0 0 1 0 0 0 ] · [ v x v y r Γ a x a y ] k + V k ,
where k is the discrete time; Λ is the noise driving matrix; X ( k ) is the state at time k ; Y ( k )   is the observation matrix at time k ; W ( k )   is input white noise; and V ( k ) is observation noise.

3. Methods

3.1. Adaptive Unscented Kalman Filter (AUKF)

Compared to the unscented Kalman filter (UKF), AUKF is adaptive to process noise. UKF algorithm obtains the probability distribution of state variables by unscented transformation (UT). The nonlinear function is approximated based on probability density. The algorithm avoids the linearization of nonlinear systems [31].
The state and observation error equations can be represented as:
x k = f ( x k 1 ) + Λ W k
y k = h ( x k ) + V k ,
where x k is the L-dimension state parameter of t k epoch;   Λ is the noise driving matrix; W k is system noise matrix; Q k , R k is the covariance matrices of W k , V k .
The description of AUKF is shown as follows [32,33]:
1.
Sampling points construction. For L-dimension column-vector X , sampling points χ k 1 can be constructed by its estimation value χ ^ k 1 and variance P k 1 .
{ χ i , k 1 = x ^ k 1 i = 0 χ i , k 1 = x ^ k 1 + [ ( L + λ ) P k 1 ] i i = 1 , , L χ i , k 1 = x ^ k 1 [ ( L + λ ) P k 1 ] i i = L + 1 , , 2 L ,
where λ = α 2 ( L + κ ) L , the selection of κ should ensure ( L + κ ) P k 1 to be a positive semidefinite matrix, α affects the distribution of sampling points and should be within 0 to 1. [ ( L + λ ) P k 1 ] i is the ith column of the matrix [ ( L + λ ) P k 1 ] . In this paper, we assign α = 0.01 , κ = 0 .
2.
Time update.
χ i , k | k 1 = f ( χ i , k 1 ) , i = 0 , 1 , , 2 L
x ^ k | k 1 = i = 0 2 L ω i χ i , k | k 1
P k | k 1 = i = 0 2 L ω i c ( χ i , k | k 1 x ^ k | k 1 ) ( χ i , k | k 1 x ^ k | k 1 ) T ,
where ω i , ω i c are the weights of sampling points. Equations (11)–(13) calculates the one-step prediction of Sigma points, the one-step prediction of system state variables, and the covariance matrix of system state variables. Another set of Sigma points needed to be constructed based on Equation (10).
3.
Measurement update.
y ^ i , k | k 1 = h ( χ i , k | k 1 )
y ^ k = i = 0 2 L ω i y ^ i , k | k 1 + σ ^ k 1
P y y = i = 0 2 L [ ω i ( y ^ i , k k 1 y ^ k ) ( y ^ i , k k 1 y ^ k ) T ] + R k
P x y = i = 0 2 L [ ω i ( χ i , k | k 1 x ^ k | k 1 ) ( y ^ i , k | k 1 y ^ k ) T ]
K k + 1 = P x y P y y 1
x ^ k = x ^ k | k 1 + K k ( y k y ^ k )
P k = P k | k 1 K k P y y K k T .
Equations (14)–(20) calculate the estimated observation variable, the mean of estimated observation variable, the covariance matrix of estimated observation variable, the covariance matrix of true observation variable, Kalman gain, state update, and covariance matrix update. σ k is the average value of noise V k , and is usually assigned to be zero in AUKF.
4.
Noise characteristics update.
The innovation is defined as:
ε k = y k y ^ k | k 1 .
The real-time estimated covariance matrix of the innovation sequence is calculated as:
P ε ( k ) = 1 M m = 0 M 1 [ ε ( k m ) ε T ( k m ) ]   ,
where P ε ( k ) is the covariance and M is the scale coefficient.
The update of the R k is
R k = P ε ( k ) i = 0 2 L [ ω i ( χ i , k | k 1 y ^ k | k 1 ) ( χ i , k | k 1 y ^ k | k 1 ) T ]   .
This process estimates measurement noise and shows adaptiveness.

3.2. Divergence Calculation

Suppose there is a system described by the following state-space model
X k + 1 = ϕ X k + Λ W k
Y k = H X k + V k ,
where k is the discrete time; ϕ is the state transition matrix; H is the observation matrix; Λ is the noise driving matrix; X k is the state at time k ; Y k   is the observation matrix at time k ; W k   is input white noise; V k is observation noise.
We assume that W k and V k are zero mean; the covariance matrices of W k and V k is Q k and R k ; E [ W k , V j ] = 0 , where j represents random time.
Innovation is defined as
ε k + 1 = Y k + 1 Y ^ k + 1 | k .
We have
Y ^ k + 1 = H X ^ k + 1 | k .
We can obtain Equation (28) if we substitute Equations (25) and (27) into Equation (26):
ε k + 1 = H ( X k + 1 X ^ k + 1 | k ) + V k + 1 .
The estimated variance of the measured value is
X ˜ k + 1 | k = X k + 1 X ^ k + 1 | k .
Then
ε k + 1 = H X ˜ k + 1 | k + V k + 1 .
The innovation covariance matrix can be expressed as
C o v ( ε k + 1 · ε k + 1 T ) = H P k + 1 H T + R k + 1
where P k + 1 is the covariance matrix of the predicted error.
When the estimation converges, the innovation ε k is white noise that obeys the normal distribution. Hence, we can obtain Equation (32), where ζ x y is the correlation coefficient, and the value is between zero and one. The closer ζ x y approaches to zero, the weaker the correlation.
ζ x y = C o v ( X , Y ) V a r ( X ) · V a r ( Y ) .
The correlation coefficient of innovation can be used to judge the divergence. The critical value for identifying the filter divergence in this paper is ζ 0 = 0.4 . If Equation (33) is satisfied, the filter diverges.
ζ 0 · V a r ( X ) · V a r ( Y ) > C o v ( X , Y ) .
Therefore, the discriminant formula can be obtained by substituting Equation (31) into Equation (33):
ζ 0 · V a r ( ε k + 1 ) < H P k + 1 H T + R k + 1 .
Since V k + 1 and H X ˜ k + 1 | k are not related, then
P { | H X ˜ k + 1 | k | 2.5 ( H P k + 1 H T ) } = 98.8 %
P { | V k + 1 | 2.5 ( R k + 1 ) } = 98.8 % .
When ζ 0 = 0.4 , there is at least a 97.61% probability for successful recognition, and the error rate is 0.04–2.39%.

3.3. Modified Sage–Husa Unscented Kalman Filter

Based on the traditional UKF, Sage–Husa algorithm, and divergence calculation technique, we propose a modified Sage–Husa UKF algorithm. The algorithm is designed to be more applicable in the field of vehicle driving state estimation than classical AUKF. The description of the proposed algorithm is shown as follows:
1.
Sampling points construction and time update. These two stages are the same as the classical AUKF mentioned above.
2.
Divergence calculation. Based on covariance matching, the filter is convergent if and only if
ζ 0 · V a r ( ε k ) < H P k H T + R k ,
where ζ 0 = 0.4 is the predetermined correlation coefficient; R k is the covariance matrix of V k .
If the filter diverges, we should correct the covariance matrix P k | k 1 [34]:
P k | k 1 = λ k i = 0 2 n [ W i c ( x k | k 1 i x ^ k | k 1 ) ( x k | k 1 i x ^ k | k 1 ) T ] + Q k 1
λ k = { λ 0 λ 0 > 1 1 λ 0 1
λ 0 = t r [ C 0 , k R k ] t r ( i = 0 2 L ω i ( y ^ i , k k 1 y ^ k ) ( y ^ i , k k 1 y ^ k ) T )
C 0 , k = { ε k ε k T k = 0 ρ C 0 , k + ε k ε k T 1 + ρ k 1
where ρ is within zero to one and assigned to be 0.95. A large ρ will reduce the residual information effect before time k , and strengthen the effect of the current residual vector.
If the filter is convergent, then
P k | k 1 = i = 0 2 L [ ω i c ( χ i , k | k 1 x ^ k | k 1 ) ( χ i , k | k 1 x ^ k | k 1 ) T ] + Q k ,
which is the same as classical AUKF.
3.
Measurement update is the same as classical AUKF.
4.
Noise characteristics update [35].
R ^ k = ( 1 d k ) R ^ k 1 + d k ( v k v k T H k P k | k 1 H k T )
σ ^ k = ( 1 d k ) r ^ k 1 + d k [ y k i = 0 2 n [ ω i m h ( x i , k + 1 | k , u k ) ] ]
d k = ( 1 b ) ( 1 b k + 1 ) ,
where H k is the transition matrix of the nonlinear equation measuring equation.   d k ( x ^ k | k i = 0 2 L ω i m f ( x k | k 1 i ) ) only has a small effect on the system and can be ignored.

4. Experiments and Analysis

The experiments are performed on CarSim and Simulink. We establish a simulation model based on the models and algorithms mentioned in Section 2 and Section 3. The Modified Sage–Husa UKF algorithm is tested on both double lane change and Sine Wave road conditions. We focus on five parameters (longitudinal speed, lateral speed, yaw rate, longitudinal acceleration, and lateral acceleration) because these parameters are sensitive to changes and more state variables can be further calculated based on these five parameters. We first present normal simulation tests to verify that the proposed approach can reach or even superior to the performance of classical AUKF. Since our approach mainly aims at the weak robustness and the slow adaptive speed of AUKF, we present a robustness test to compare the performances of two methods under severe disturbances.

4.1. Simulation Platform

CarSim released by MSC in 1996 is a vehicle dynamics simulation software. It can simulate the vehicle response to the road, air, and driver behavior in the real driving process. Carsim can simulate three to six times the speed of actual motion and allows users to manually set parameters. It can directly establish a Simulink model, which is friendly with MATLAB users. This paper verifies our proposed approach based on CarSim and Simulink co-simulation. Driving state data is provided by CarSim and the on-line estimation on the modified Sage–Husa UKF algorithm is performed on Simulink.
We first set the road condition as double lane change in Carsim, with a maximum speed of 60 km/h. We then test the performance of the algorithm under Sine Wave road condition, with a maximum speed of 90 km/h. The parameters of the vehicle model are shown in Table 1.
We import the simulation data from CarSim to Simulink and co-simulate with the established Dugoff tire model and seven-degrees-of-freedom vehicle model.

4.2. Simulation Model

The established simulation model is shown in Figure 3. The output state variables of CarSim is time linear, including wheel steering angle δ , longitudinal acceleration a x , lateral acceleration a y , and yaw rate r. We can calculate the current tire longitudinal and lateral force F x i and F y i by inputting longitudinal and lateral acceleration, yaw rate, wheel steering angle, and estimated longitudinal and lateral speed into the Dugoff tire model. We can obtain current estimated longitudinal and lateral acceleration, yaw rate, and longitudinal and lateral speed by inputting F x i , F y i and current measured accelerate and yaw rate into Vehicle Driving State Observer based on Modified Sage–Husa UKF algorithm. Longitudinal and lateral speed are regarded as output values and are involved with tire force calculation of the next moment. Here we suppose the acceleration only changes a little and the simulation result proves our hypothesis.
The main difference between the classical AUKF vehicle driving state estimation model and our model is the adaptive part. Instead of the classical AUKF algorithm, the proposed Vehicle Driving State Observer is based on the Modified Sage–Husa UKF algorithm.
We compare the proposed Modified Sage–Husa UKF algorithm and classical AUKF based on double lane change and Sine Wave road conditions. The trajectories of the two road conditions are shown in Figure 4.

4.3. Simulation Results and Analysis

In the simulation, the initial measurement noise covariance matrix R 0 is
R 0 = [ 0.1 0.1 0.1 ]
The initial process noise covariance matrix Q 0 is
Q 0 = [ 0.1 0.1 0.1 0.1 1 0.1 ]
The modified Sage–Husa UKF algorithm will adaptively iterate covariance matrix R. Hence the initial condition can be randomly selected within the proper range and the matrix Q must be positive definite.

4.3.1. Double Lane Change Road Condition

Figure 5a, Figure 6a, Figure 7a, Figure 8a, Figure 9a demonstrate the AUKF and modified Sage–Husa UKF algorithm estimation results. Figure 5b, Figure 6b, Figure 7b, Figure 8b, Figure 9b indicate the errors of two methods compared with true values. We apply RMSE (root mean square error) and MAE (mean absolute error) error evaluation and the results are shown in Table 2.
According to Table 2, the proposed Modified Sage–Husa UKF algorithm is better than AUKF in all five state variables and both two aspects. The RMSE values of longitudinal speed, lateral speed, yaw rate, longitudinal and lateral acceleration are improved by 0.71%, 2.4%, 30.7%, 4.9%, and 17.8% respectively compared to AUKF. The MAE values of these five state variables are improved by 0.68%, 2.4%, 30.0%, 3.4% and 19.7% respectively.

4.3.2. Sine Wave Road Condition

Figure 10a, Figure 11a, Figure 12a, Figure 13a, Figure 14a demonstrate the AUKF and modified Sage–Husa UKF algorithm estimation results. Figure 10b, Figure 11b, Figure 12b, Figure 13b, Figure 14b indicate the errors of two methods compared with true values. We apply RMSE (root mean square error) and MAE (mean absolute error) error evaluation and the results are shown in Table 3.
According to Table 3, the proposed Modified Sage–Husa UKF algorithm is better than AUKF in all five state variables and both two aspects. The RMSE values of longitudinal speed, lateral speed, yaw rate, longitudinal and lateral acceleration are improved by 1.19%, 0.63%, 1.21%, 53.45%, and 39.11% respectively compared to AUKF. The MAE values of these five state variables are improved by 1.97%, 2.55%, 2.38%, 40.06%, and 35.37% respectively.

4.3.3. Discussion

Based on the results in Section 4.3.1 and Section 4.3.2, it can be seen that:
  • Our approach is with high accuracy and easy to converge. The MAE values near zero indicate small average deviations. Besides the involvement of the divergence calculation technique, we improve the filter as well. It should be noted that AUKF is so compact that even a little modification on the filter may result in lower accuracy. According to Table 2 and Table 3, in the normal road condition, the proposed Modified Sage–Husa UKF Algorithm can rival and even better than the performance of classical AUKF.
  • Our approach is adaptive to the changing system and remains excellent performance when the road condition changes. The initial conditions in Section 4.3.1 and Section 4.3.2 are the same. Since the initial matrices R 0 and Q 0 are set under double lane change road condition, AUKF performs undesirable in the Sine Wave road condition. Especially in Figure 13b and Figure 14b, the proposed approach demonstrates its superiority.
  • Our approach can update the covariance matrix. The introduction of the divergence calculation in our approach effectively improves the adaptivity and anti-interference ability. According to Figure 13b, the error curve of AUKF estimation is rough, which is with changeable extremum. This is caused by the comparatively more extreme road condition, which challenges the sensors and observers. AUKF performs undesirable because of its slow updating of the noise covariance matrix. The involvement of divergence calculation in the proposed method accelerates the process of convergence and performs desirable when the road condition changes.
  • It can be noticed that the estimation is undesirable in Figure 11a. This is because the Sine Wave road condition under 90 km/h is so extreme that most vehicle states reach their extreme values. Under this road condition, the direction of lateral speed changes with high frequency and small amplitude. The undesirable estimation of tire forces negatively affects the estimation of lateral speed. Nevertheless, the proposed approach tracks the changing trend and keeps convergent, which demonstrates its stability. Besides, it is comparatively unsafe to drive under 90 km/h in a real scenario. The test in Section 4.3.1 indicates that our approach can accurately estimate vehicle driving states under safe driving conditions.

4.4. Robustness Test

We simulate the change of measurement error or system error in the driving process. The purpose of this test is to simulate transient disturbances caused by an onboard processor failure in extreme conditions.

4.4.1. Double Lane Change Road Condition

We add a drastic disturbance with a length of 0.1 s at the 6.58th second. In this 0.1 s, longitudinal, lateral speed, yaw rate, longitudinal and lateral acceleration are set to be 15 m/s, 0.2 m/s, 0.2 rad/s, 15 m/s2, 15 m/s2 respectively when the real value is 16.67 m/s, 0.01 m/s, 0.01 rad/s, 0.21 m/s2, −0.12 m/s2 respectively.
Figure 15a, Figure 16a, Figure 17a, Figure 18a, Figure 19a demonstrate the AUKF and Modified Sage–Husa UKF algorithm estimation results. Figure 15b, Figure 16b, Figure 17b, Figure 18b, Figure 19b indicate the errors of the two methods compared to the true values. We apply RMSE and MAE error evaluation and the results are shown in Table 4.
The proposed Modified Sage–Husa UKF algorithm is better than AUKF in most state variables. The RMSE values of longitudinal speed, yaw rate, longitudinal and lateral acceleration are improved by 37.1%, 20.8%, 44.5%, and 76.1% respectively compared to AUKF. The MAE values of longitudinal speed, yaw rate, longitudinal and lateral acceleration are improved by 3.9%, 10.3%, 22.6%, and 19.4% respectively. As shown in Figure 15b, Figure 16b, Figure 17b, Figure 18b, Figure 19b, the Modified Sage–Husa UKF algorithm and AUKF both converge quickly to the true value after being disturbed. For the states shown in Figure 17a, Figure 18a, Figure 19a which exist measured values, the convergence speed of our approach is much faster than AUKF. In these three situations, our approach has improved convergence speed by 75%, 33%, and 67.2% respectively.

4.4.2. Sine Wave Road Condition

We add a drastic disturbance with a length of 0.015 s at the 10th second. In this 0.015 s, longitudinal, lateral speed, yaw rate, longitudinal and lateral acceleration are set to be 22 m/s, 1 m/s, 0.2 rad/s, −17 m/s2, −8 m/s2 respectively when the real value is 24.58 m/s, 0.04 m/s, 0.02 rad/s, −4.03 m/s2, −0.27 m/s2 respectively.
Figure 20a, Figure 21a, Figure 22a, Figure 23a, Figure 24a demonstrate the AUKF and Modified Sage–Husa UKF algorithm estimation results. Figure 20b, Figure 21b, Figure 22b, Figure 23b, Figure 24b indicate the errors of the two methods compared to the true values. We apply RMSE and MAE error evaluation and the results are shown in Table 5.
According to Table 5, the proposed Modified Sage–Husa UKF algorithm is better than AUKF in all five state variables and both two aspects. The RMSE values of longitudinal speed, lateral speed, yaw rate, longitudinal and lateral acceleration are improved by 38.05%, 6.35%, 14.40%, 41.31%, and 35.90% respectively compared to AUKF. The MAE values of these five state variables are improved by 12.84%, 4.22%, 16.67%, 3.61%, and 21.56% respectively.

4.4.3. Discussion

Based on the results in Section 4.4.1 and Section 4.4.2, it can be seen that:
  • The 57.27% MAE value improvement of lateral speed under double lane change road condition demonstrates that the transient disturbance impacts our approach less and does not affect the subsequent estimation.
  • Our approach has a better anti-interference ability than AUKF. According to Figure 15b, Figure 16b, Figure 17b, Figure 18b, Figure 19b, Figure 20b, Figure 21b, Figure 22b, Figure 23b, Figure 24b, the AUKF is sensitive to disturbances and takes a longer time at adapting.
  • As shown in Figure 9b, Figure 10b, Figure 11b, Figure 12b, Figure 13b, the accuracies of our approach are more stable when facing severe disturbances. This indicates that our approach can adapt to the most optimistic measurement noise covariance matrix.
  • As shown in the test under double lane change road condition, our approach can reflect the small state changes during the driving process, which is meaningful and important to the vehicle safety control system. This is because AUKF uses a windowing mask to eliminate the influence of state change on current estimation. The proportional coefficient M is set according to experiences. M cannot be determined or will be inaccurate if unknown changes occur.
  • According to Table 4, it can be noticed that in the lateral speed comparison, the RMSE value of the proposed approach is larger than that of the classical AUKF. There are two reasons. Firstly, the larger RMSE value here indicates that the proposed approach is more sensitive to the state change. It needs to be acknowledged that since we are performing on-line estimation when the disturbance is added, the system cannot distinguish whether the changing observation value is caused by a sudden change of vehicle state or just an exerted disturbance. The standard error is sensitive to the extremely large value or the extremely small value during the measurement. The proposed approach is sensitive to the state change, and the exerted disturbance is twenty times of the true value at that moment. Secondly, our approach is with a smaller MAE value under the circumstance of a larger RMSE value. This indicates our approach is with stronger robustness than classical AUKF. MAE value is the absolute value summation of the difference between the true value and the estimated value, which reflects the stability of the system.
  • According to Table 4 and Table 5, the proposed approach remains high accuracy and stability when the transient disturbance occurs, and it is with stronger robustness comparatively.

5. Conclusions

When the road condition changes frequently and drastically, traditional vehicle driving state estimation methods perform undesirable. The unknown measurement noise and transient disturbances will cause the failure of state estimation and affect the subsequent process. We propose an AUKF based on a modified Sage–Husa filter and divergence calculation technique to solve these problems. We mainly improve the classical AUKF algorithm in two aspects. First, we introduce the divergence calculation technique to detect transient disturbances and solve the divergence problem of Sage–Husa maximum posteriori. Second, we modify the covariance matrix when the transient disturbance occurs.
Our approach can excellently face transient disturbances and has strong robustness through standard simulation tests and robustness tests under double lane change and Sine Wave road conditions. The proposed Modified Sage–Husa UKF algorithm fully satisfies the actual engineering precision requirements under simulation. The results of our approach rival the measurements obtained from high-cost and high-accuracy sensors. Besides, the simulation results show the feasibility of applying our approach to vehicle driving state estimation.
In our future work, we will first improve our algorithm, including the verification of more complicated working conditions and applications in a real car. Then we will add application scenarios for vehicle driving state estimation. The estimated data can combine with IoT (Internet of Things) and big data. For example, the vehicle driving states of a car getting out of control can share with surrounding cars. This will avoid secondary accidents. Data sharing could also share those repeated measuring state variables. It will save computational power and improve accuracy.

Author Contributions

Conceptualization, methodology, software, validation, formal analysis, investigation, data curation, Z.L.; software, validation, formal analysis, investigation, data curation, visualization, writing—original draft preparation, Z.F.; resources, writing—review and editing, supervision, project administration, funding acquisition, Q.X. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Chongqing Science and Technology Commission of China, grant number cstc2018jcyjA3148.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Global Status Report on Road Safety 2018. Available online: https://www.who.int/violence_injury_prevention/road_safety_status/2018/en/ (accessed on 26 October 2020).
  2. Chakraborty, S.; Sen, S.; Sutradhar, A.; Sengupta, A. Estimation of tire-road friction coefficient and frictional force for Active Vehicle safety system. In Proceedings of the International Conference on Industrial Instrumentation and Control, Pune, India, 28–30 May 2015; pp. 674–679. [Google Scholar]
  3. Sun, W.; Zhang, J.; Liu, Z. Two-Time-Scale Redesign for Antilock Braking Systems of Ground Vehicles. IEEE Trans. Ind. Electron. 2018, 66, 4577–4586. [Google Scholar] [CrossRef]
  4. Kang, M.; Li, L.; Li, H.; Song, J.; Han, Z. Coordinated vehicle traction control based on engine torque and brake pressure under complicated road conditions. Veh. Syst. Dyn. 2012, 50, 1473–1494. [Google Scholar] [CrossRef]
  5. Wang, P.; Wang, P.; Hu, Y.; Chen, H. A stability-guaranteed and energy-conserving torque distribution strategy for electric vehicles under extreme conditions. Appl. Energy 2020, 259, 114162. [Google Scholar] [CrossRef]
  6. Khalkhali, M.B.; Vahedian, A.; Sadoghi-Yazdi, H. Multi-Target State Estimation Using Interactive Kalman Filter for Multi-Vehicle Tracking. IEEE Trans. Intell. Transp. Syst. 2020, 21, 1131–1144. [Google Scholar] [CrossRef]
  7. Kim, B.; Yi, K.; Yoo, H.-J.; Chong, H.-J.; Ko, B. An IMM/EKF Approach for Enhanced Multitarget State Estimation for Application to Integrated Risk Management System. IEEE Trans. Veh. Technol. 2015, 64, 876–889. [Google Scholar] [CrossRef]
  8. Woo, R.; Yang, E.-J.; Seo, D.-W. A Fuzzy-Innovation-Based Adaptive Kalman Filter for Enhanced Vehicle Positioning in Dense Urban Environments. Sensors 2019, 19, 1142. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, Y.; Li, J.; Guo, Y.; Xu, C.; Bao, J.; Song, Y. Vehicle Driving Behavior Recognition Based on Multi-View Convolutional Neural Network With Joint Data Augmentation. IEEE Trans. Veh. Technol. 2019, 68, 4223–4234. [Google Scholar] [CrossRef]
  10. Guzmán, J.G.; González, L.P.; Redondo, J.P.; Martínez, M.M.M.; Boada, M.J.L. Real-Time Vehicle Roll Angle Estimation Based on Neural Networks in IoT Low-Cost Devices. Sensors 2018, 18, 2188. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Wu, L.-J. Experimental study on vehicle speed estimation using accelerometer and wheel speed measurements. In Proceedings of the Second International Conference on Mechanic Automation and Control Engineering, Inner Mongolia, China, 15–17 July 2011; pp. 294–297. [Google Scholar]
  12. Ray, L. Nonlinear state and tire force estimation for advanced vehicle control. IEEE Trans. Control Syst. Technol. 1995, 3, 117–124. [Google Scholar] [CrossRef]
  13. Doumiati, M.; Victorino, A.C.; Charara, A.; Lechner, D. Onboard Real-Time Estimation of Vehicle Lateral Tire–Road Forces and Sideslip Angle. IEEE/ASME Trans. Mechatron. 2011, 16, 601–614. [Google Scholar] [CrossRef]
  14. Li, S.E.; Li, G.; Yu, J.; Liu, C.; Cheng, B.; Wang, J.; Li, K. Kalman filter-based tracking of moving objects using linear ultrasonic sensor array for road vehicles. Mech. Syst. Signal Process. 2018, 98, 173–189. [Google Scholar] [CrossRef]
  15. Saadeddin, K.; Abdel-Hafez, M.F.; Al-Jarrah, M.A. Estimating Vehicle State by GPS/IMU Fusion with Vehicle Dynamics. J. Intell. Robot. Syst. 2013, 74, 147–172. [Google Scholar] [CrossRef]
  16. Havyarimana, V.; Xiao, Z.; Wang, D. A Hybrid Approach-based Sparse Gaussian Kernel Model for Vehicle State Estimation during the Free and Complete GPS Outages. ETRI J. 2016, 38. [Google Scholar] [CrossRef]
  17. Xu, S.; Zhou, H.; Wang, J.; He, Z.; Wang, D. SINS/CNS/GNSS Integrated Navigation Based on an Improved Federated Sage–Husa Adaptive Filter. Sensors 2019, 19, 3812. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Huang, B.; Wu, S.; Fu, X.; Luo, J. State estimation of four-wheel independent drive electric vehicle based on adaptive unscented Kalman filter. Int. J. Electr. Hybrid Veh. 2017, 9, 151–168. [Google Scholar] [CrossRef]
  19. Wu, Q.; Jia, Q.; Shan, J.; Meng, X. Angular velocity estimation based on adaptive simplified spherical simplex unscented Kalman filter in GFSINS. Proc. Inst. Mech. Eng. Part G J. Aerosp. Eng. 2014, 228, 1375–1388. [Google Scholar] [CrossRef]
  20. Li, G.; Zhao, D.; Xie, R.; Han, H.; Zong, C. Vehicle State Estimation Based on Improved Sage–Husa Adaptive Extended Kalman Filtering. Automot. Eng. 2015, 37, 1426–1432. [Google Scholar]
  21. Lin, M.; Yoon, J.; Kim, B. Self-Driving Car Location Estimation Based on a Particle-Aided Unscented Kalman Filter. Sensors 2020, 20, 2544. [Google Scholar] [CrossRef]
  22. Guo, H.; Guo, J.; Yu, M.; Hong, H.; Xiong, J.; Tian, B. A weighted combination filter with nonholonomic constrains for integrated navigation systems. Adv. Space Res. 2015, 55, 1470–1476. [Google Scholar] [CrossRef]
  23. Zhang, J.; Sun, T.; Liu, Z. Robust model predictive control for path-following of underactuated surface vessels with roll constraints. Ocean Eng. 2017, 143, 125–132. [Google Scholar] [CrossRef]
  24. Gao, X.; You, D.; Katayama, S. Seam Tracking Monitoring Based on Adaptive Kalman Filter Embedded Elman Neural Network During High-Power Fiber Laser Welding. IEEE Trans. Ind. Electron. 2012, 59, 4315–4325. [Google Scholar] [CrossRef]
  25. Wang, Y.; Sun, Y.; Dinavahi, V. Robust Forecasting-Aided State Estimation for Power System against Uncertainties. IEEE Trans. Power Syst. 2020, 35, 691–702. [Google Scholar] [CrossRef]
  26. Zhao, Z.; Chen, H.; Yang, J.; Wu, X.; Yu, Z. Estimation of the vehicle speed in the driving mode for a hybrid electric car based on an unscented Kalman filter. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2014, 229, 437–456. [Google Scholar] [CrossRef]
  27. Pacejka, H.B.; Besselink, I.J.M. Magic Formula Tyre Model with Transient Properties. Veh. Syst. Dyn. 1996, 27, 243–249. [Google Scholar] [CrossRef]
  28. Dugoff, H.; Fancher, P.S.; Segel, L. An Analysis of Tire Traction Properties and Their Influence on Vehicle Dynamic Performance. SAE Tech. Pap. Ser. 1970, 1219–1243. [Google Scholar] [CrossRef]
  29. Li, J.; Zhang, Y.; Yi, J. A Hybrid Physical-Dynamic Tire/Road Friction Model. J. Dyn. Syst. Meas. Control 2012, 135, 011007. [Google Scholar] [CrossRef]
  30. Han, K.S.; Lee, E.; Choi, S. Estimation of the maximum lateral tire-road friction coefficient using the 6-DoF sensor. In Proceedings of the 15th International Conference on Control, Automation and Systems, Busan, Korea, 13–16 October 2015; pp. 1734–1738. [Google Scholar]
  31. Yang, Y.; Gao, W. An Optimal Adaptive Kalman Filter. J. Geod. 2006, 80, 177–183. [Google Scholar] [CrossRef]
  32. Costanzi, R.; Fanelli, F.; Meli, E.; Ridolfi, A.; Caiti, A.; Allotta, B. UKF-Based Navigation System for AUVs: Online Experimental Validation. IEEE J. Ocean. Eng. 2018, 44, 633–641. [Google Scholar] [CrossRef]
  33. Zheng, B.; Fu, P.; Li, B.; Yuan, X. A Robust Adaptive Unscented Kalman Filter for Nonlinear Estimation with Uncertain Noise Covariance. Sensors 2018, 18, 808. [Google Scholar] [CrossRef] [Green Version]
  34. Gao, S.; Zhu, H.; Guo, Z. Adaptive Compensation Method Based on STCKF for AHRS. In Proceedings of the International Conference on Mechatronics Engineering and Information Technology, Xi’an, China, 27–28 August 2016. [Google Scholar]
  35. Sage, A.P.; Husa, G.W. Adaptive filtering with unknown prior statistics. In Joint Automatic Control Conference; IEEE: Piscataway, NJ, USA, 1969; Volume 7, pp. 760–769. [Google Scholar]
Figure 1. Schematic of seven-degrees-of-freedom vehicle model.
Figure 1. Schematic of seven-degrees-of-freedom vehicle model.
Sensors 20 06889 g001
Figure 2. Schematic of Dugoff tire model.
Figure 2. Schematic of Dugoff tire model.
Sensors 20 06889 g002
Figure 3. Schematic of the estimation process.
Figure 3. Schematic of the estimation process.
Sensors 20 06889 g003
Figure 4. Schematic of testing road condition trajectory. (a) The trajectory of double lane change road condition. (b) The trajectory of Sine Wave road condition.
Figure 4. Schematic of testing road condition trajectory. (a) The trajectory of double lane change road condition. (b) The trajectory of Sine Wave road condition.
Sensors 20 06889 g004
Figure 5. The simulation results of longitudinal speed under double lane change road condition. (a) The estimation comparison of longitudinal speed. (b) The error comparison of longitudinal speed estimation.
Figure 5. The simulation results of longitudinal speed under double lane change road condition. (a) The estimation comparison of longitudinal speed. (b) The error comparison of longitudinal speed estimation.
Sensors 20 06889 g005
Figure 6. The simulation results of lateral speed under double lane change road condition. (a) The estimation comparison of lateral speed. (b) The error comparison of lateral speed estimation.
Figure 6. The simulation results of lateral speed under double lane change road condition. (a) The estimation comparison of lateral speed. (b) The error comparison of lateral speed estimation.
Sensors 20 06889 g006
Figure 7. The simulation results of yaw rate under double lane change road condition. (a) The estimation comparison of yaw rate. (b) The error comparison of yaw rate estimation.
Figure 7. The simulation results of yaw rate under double lane change road condition. (a) The estimation comparison of yaw rate. (b) The error comparison of yaw rate estimation.
Sensors 20 06889 g007
Figure 8. The simulation results of longitudinal acceleration under double lane change road condition. (a) The estimation comparison of longitudinal acceleration. (b) The error comparison of longitudinal acceleration estimation.
Figure 8. The simulation results of longitudinal acceleration under double lane change road condition. (a) The estimation comparison of longitudinal acceleration. (b) The error comparison of longitudinal acceleration estimation.
Sensors 20 06889 g008
Figure 9. The simulation results of lateral acceleration under double lane change road condition. (a) The estimation comparison of lateral acceleration. (b) The error comparison of lateral acceleration estimation.
Figure 9. The simulation results of lateral acceleration under double lane change road condition. (a) The estimation comparison of lateral acceleration. (b) The error comparison of lateral acceleration estimation.
Sensors 20 06889 g009
Figure 10. The simulation results of longitudinal speed under Sine Wave road condition. (a) The estimation comparison of longitudinal speed. (b) The error comparison of longitudinal speed estimation.
Figure 10. The simulation results of longitudinal speed under Sine Wave road condition. (a) The estimation comparison of longitudinal speed. (b) The error comparison of longitudinal speed estimation.
Sensors 20 06889 g010
Figure 11. The simulation results of lateral speed under Sine Wave road condition. (a) The estimation comparison of lateral speed. (b) The error comparison of lateral speed estimation.
Figure 11. The simulation results of lateral speed under Sine Wave road condition. (a) The estimation comparison of lateral speed. (b) The error comparison of lateral speed estimation.
Sensors 20 06889 g011
Figure 12. The simulation results of yaw rate under Sine Wave road condition. (a) The estimation comparison of yaw rate. (b) The error comparison of yaw rate estimation.
Figure 12. The simulation results of yaw rate under Sine Wave road condition. (a) The estimation comparison of yaw rate. (b) The error comparison of yaw rate estimation.
Sensors 20 06889 g012
Figure 13. The simulation results of longitudinal acceleration under Sine Wave road condition. (a) The estimation comparison of longitudinal acceleration. (b) The error comparison of longitudinal acceleration estimation.
Figure 13. The simulation results of longitudinal acceleration under Sine Wave road condition. (a) The estimation comparison of longitudinal acceleration. (b) The error comparison of longitudinal acceleration estimation.
Sensors 20 06889 g013
Figure 14. The simulation results of lateral acceleration under Sine Wave road condition. (a) The estimation comparison of lateral acceleration. (b) The error comparison of lateral acceleration estimation.
Figure 14. The simulation results of lateral acceleration under Sine Wave road condition. (a) The estimation comparison of lateral acceleration. (b) The error comparison of lateral acceleration estimation.
Sensors 20 06889 g014
Figure 15. The simulation results of longitudinal speed in robustness test under double lane change road condition. (a) The estimation comparison of longitudinal speed in robustness test. (b) The error comparison of longitudinal speed estimation in robustness test.
Figure 15. The simulation results of longitudinal speed in robustness test under double lane change road condition. (a) The estimation comparison of longitudinal speed in robustness test. (b) The error comparison of longitudinal speed estimation in robustness test.
Sensors 20 06889 g015
Figure 16. The simulation results of lateral speed in robustness test under double lane change road condition. (a) The estimation comparison of lateral speed in robustness test. (b) The error comparison of lateral speed estimation in robustness test.
Figure 16. The simulation results of lateral speed in robustness test under double lane change road condition. (a) The estimation comparison of lateral speed in robustness test. (b) The error comparison of lateral speed estimation in robustness test.
Sensors 20 06889 g016
Figure 17. The simulation results of yaw rate in robustness test under double lane change road condition. (a) The estimation comparison of yaw rate in robustness test. (b) The error comparison of yaw rate estimation in robustness test.
Figure 17. The simulation results of yaw rate in robustness test under double lane change road condition. (a) The estimation comparison of yaw rate in robustness test. (b) The error comparison of yaw rate estimation in robustness test.
Sensors 20 06889 g017
Figure 18. The simulation results of longitudinal acceleration in robustness test under double lane change road condition. (a) The estimation comparison of longitudinal acceleration in robustness test. (b) The error comparison of longitudinal acceleration estimation in robustness test.
Figure 18. The simulation results of longitudinal acceleration in robustness test under double lane change road condition. (a) The estimation comparison of longitudinal acceleration in robustness test. (b) The error comparison of longitudinal acceleration estimation in robustness test.
Sensors 20 06889 g018
Figure 19. The simulation results of lateral acceleration in robustness test under double lane change road condition. (a) The estimation comparison of lateral acceleration in robustness test. (b) The error comparison of lateral acceleration estimation in robustness test.
Figure 19. The simulation results of lateral acceleration in robustness test under double lane change road condition. (a) The estimation comparison of lateral acceleration in robustness test. (b) The error comparison of lateral acceleration estimation in robustness test.
Sensors 20 06889 g019
Figure 20. The simulation results of longitudinal speed in robustness test under sine wave road condition. (a) The estimation comparison of longitudinal speed in robustness test. (b) The error comparison of longitudinal speed estimation in robustness test.
Figure 20. The simulation results of longitudinal speed in robustness test under sine wave road condition. (a) The estimation comparison of longitudinal speed in robustness test. (b) The error comparison of longitudinal speed estimation in robustness test.
Sensors 20 06889 g020
Figure 21. The simulation results of lateral speed in robustness test under sine wave road condition. (a) The estimation comparison of lateral speed in robustness test. (b) The error comparison of lateral speed estimation in robustness test.
Figure 21. The simulation results of lateral speed in robustness test under sine wave road condition. (a) The estimation comparison of lateral speed in robustness test. (b) The error comparison of lateral speed estimation in robustness test.
Sensors 20 06889 g021
Figure 22. The simulation results of yaw rate in robustness test under sine wave road condition. (a) The estimation comparison of yaw rate in robustness test. (b) The error comparison of yaw rate estimation in robustness test.
Figure 22. The simulation results of yaw rate in robustness test under sine wave road condition. (a) The estimation comparison of yaw rate in robustness test. (b) The error comparison of yaw rate estimation in robustness test.
Sensors 20 06889 g022
Figure 23. The simulation results of longitudinal acceleration in robustness test under sine wave road condition. (a) The estimation comparison of longitudinal acceleration in robustness test. (b) The error comparison of longitudinal acceleration estimation in robustness test.
Figure 23. The simulation results of longitudinal acceleration in robustness test under sine wave road condition. (a) The estimation comparison of longitudinal acceleration in robustness test. (b) The error comparison of longitudinal acceleration estimation in robustness test.
Sensors 20 06889 g023
Figure 24. The simulation results of lateral acceleration in robustness test under sine wave road condition. (a) The estimation comparison of lateral acceleration in robustness test. (b) The error comparison of lateral acceleration estimation in robustness test.
Figure 24. The simulation results of lateral acceleration in robustness test under sine wave road condition. (a) The estimation comparison of lateral acceleration in robustness test. (b) The error comparison of lateral acceleration estimation in robustness test.
Sensors 20 06889 g024
Table 1. Parameters of the simulation vehicle model.
Table 1. Parameters of the simulation vehicle model.
ParameterSymbolValue
Mass of the vehicle (kg) m 1200
Height of mass center (m) h c 0.375
Distance from the front axis to the mass center (m) a 1.455
Distance from the rear axis to the mass center (m) b 1.195
Length of the front axis (m) T f 1.675
Length of the rear axis (m) T r 1.675
Inertia around the z-axis ( k g · m 2 ) I z 1.652
Table 2. Error evaluation under double lane change road condition.
Table 2. Error evaluation under double lane change road condition.
State VariablesAUKFModified Sage–Husa UKF
RMSEMAERMSEMAE
Longitudinal speed0.11240.08810.11160.0875
Lateral speed0.02070.01240.02020.0121
Yaw rate0.01400.01100.00970.0077
Longitudinal acceleration0.09640.07550.09170.0729
Lateral acceleration0.09260.07450.07610.0598
Table 3. Error evaluation under Sine Wave road condition.
Table 3. Error evaluation under Sine Wave road condition.
State VariablesAUKFModified Sage–Husa UKF
RMSEMAERMSEMAE
Longitudinal speed0.06710.05590.06630.0548
Lateral speed0.30260.21610.30070.2106
Yaw rate0.02480.02100.02450.0205
Longitudinal acceleration0.21870.13330.10180.0799
Lateral acceleration0.03810.03110.02320.0201
Table 4. Error evaluation in robustness test under double lane change road condition.
Table 4. Error evaluation in robustness test under double lane change road condition.
State VariablesAUKFModified Sage–Husa UKF
RMSEMAERMSEMAE
Longitudinal speed0.16620.10300.10450.0990
Lateral speed0.08630.04190.10180.0179
Yaw rate0.03270.02330.02590.0209
Longitudinal acceleration0.18830.10230.10450.0792
Lateral acceleration0.42680.09530.10180.0768
Table 5. Error evaluation in robustness test under sine wave road condition.
Table 5. Error evaluation in robustness test under sine wave road condition.
State VariablesAUKFModified Sage–Husa UKF
RMSEMAERMSEMAE
Longitudinal speed0.14900.06930.09230.0604
Lateral speed0.36240.25610.33940.2453
Yaw rate0.02570.02160.02200.0180
Longitudinal acceleration0.16850.08030.09890.0774
Lateral acceleration0.03120.02180.02000.0171
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Luo, Z.; Fu, Z.; Xu, Q. An Adaptive Multi-Dimensional Vehicle Driving State Observer Based on Modified Sage–Husa UKF Algorithm. Sensors 2020, 20, 6889. https://doi.org/10.3390/s20236889

AMA Style

Luo Z, Fu Z, Xu Q. An Adaptive Multi-Dimensional Vehicle Driving State Observer Based on Modified Sage–Husa UKF Algorithm. Sensors. 2020; 20(23):6889. https://doi.org/10.3390/s20236889

Chicago/Turabian Style

Luo, Zeyuan, Zanhao Fu, and Qiwei Xu. 2020. "An Adaptive Multi-Dimensional Vehicle Driving State Observer Based on Modified Sage–Husa UKF Algorithm" Sensors 20, no. 23: 6889. https://doi.org/10.3390/s20236889

APA Style

Luo, Z., Fu, Z., & Xu, Q. (2020). An Adaptive Multi-Dimensional Vehicle Driving State Observer Based on Modified Sage–Husa UKF Algorithm. Sensors, 20(23), 6889. https://doi.org/10.3390/s20236889

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop