Adaptive Dual Extended Kalman Filter Based on Variational Bayesian Approximation for Joint Estimation of Lithium-Ion Battery State of Charge and Model Parameters

An accurate state of charge (SOC) estimation is vital for the safe operation and efficient management of lithium-ion batteries. At present, the extended Kalman filter (EKF) can accurately estimate the SOC under the condition of a precise battery model and deterministic noise statistics. However, in practical applications, the battery characteristics change with different operating conditions and the measurement noise statistics may vary with time, resulting in nonoptimal and even unreliable estimation of SOC by EKF. To improve the SOC estimation accuracy under uncertain measurement noise statistics, a variational Bayesian approximation-based adaptive dual extended Kalman filter (VB-ADEKF) is proposed in this paper. The variational Bayesian inference is integrated with the dual EKF (DEKF) to jointly estimate the lithium-ion battery parameters and SOC. Meanwhile, the measurement noise variances are simultaneously estimated in the SOC estimation process to compensate for the model uncertainties, so that the adaptability of the proposed algorithm to dynamic changes in battery characteristics is greatly improved. A constant current discharge test, a pulse current discharge test, and an urban dynamometer driving schedule (UDDS) test are performed to verify the effectiveness and superiority of the proposed algorithm by comparison with the DEKF algorithm. The experimental results show that the proposed VB-ADEKF algorithm outperforms the traditional DEKF algorithm in terms of SOC estimation accuracy, convergence rate, and robustness.


Introduction
Electric vehicles (EVs) are believed worldwide to be one of the most important development directions in the vehicle industry because of their advantages in low pollution and energy saving.Lithium-ion batteries, by virtue of their high energy and power density, are the fundamental power source of EVs [1].Nevertheless, lithium-ion batteries have hidden dangers in safety.Once over-charged or over-discharged, the battery capacity will drop, leading to reduced lifetime or damage, or even explosion.Moreover, the battery characteristics change with different operating conditions.Therefore, to ensure the safety, reliability, and efficiency of the battery system, a battery management system (BMS) is developed and utilised.The BMS is used to control and monitor the battery operating conditions to guarantee the safe operation and longevity of the battery.
The state of charge (SOC) is a significant parameter which indicates the amount of remaining energy in the battery for further service.It needs to be estimated accurately online in order to prevent any over-charge and over-discharge, while providing the information and support for the effective and flexible operation of the vehicle.However, it is difficult to acquire precise SOC estimates since the battery itself is a highly nonlinear system and has a lot of uncertainties.
Extensive studies have been made into SOC estimation.Hannan et al. [2] gave a comprehensive review of these research fruits.Overall, these methods can be mainly classified into four categories: (1) the open-circuit voltage (OCV) method; (2) current integration or coulomb counting (CC) method; (3) data-driven methods; and (4) model-based methods.
The OCV method is based on the one-to-one corresponding relationship of open-circuit voltage and SOC.It is simple but susceptible to temperature, battery aging, and other factors, and the acquisition of accurate measurement of OCV needs lots of rest time, making it almost impossible for moving vehicles.The CC method is widely used in many applications, but it has two weak points.First, its estimation accuracy is strongly dependent on the correctness of the initial SOC value, and second, it can easily diverge due to error accumulation in the current measurement.
The data-driven methods include neural network (NN) algorithms [3][4][5], fuzzy logic (FL) algorithms [6][7][8], support vector machines (SVMs) [9][10][11], and so on.For instance, He et al. [3] employed a multilayer feedforward neural network to estimate the SOC of Lithium-ion batteries by use of the discharge current, terminal voltage, and temperature of the battery as input.Zhao et al. [4] combined a back propagation neural network and adaptive Kalman filter to estimate the SOC, and used a forgetting-factor recursive least-squares (FFRLS) algorithm to identify the time-varying battery model parameter.In [6], a fuzzy logic algorithm was applied to estimate SOC of lithium-ion batteries for the application of a portable defibrillator.AC impedance and voltage recovery measurements were used as the input parameters for the FL model.In [7], a more advanced algorithm named adaptive neuro-fuzzy inference system (ANFIS) was developed to estimate SOC.In addition, Hu et al. [9] proposed an SOC estimation based on an optimized support vector machine for regression with double search optimization process.In [10], an SOC estimation method based on fuzzy least-squares SVM was proposed.These data-driven algorithms do not need to know any battery characteristics and have a good ability of nonlinear mapping and self-learning.However, they require a large amount of experimental data to train the intelligent model beforehand.Different types of training data and training methods exert a great influence on the model error.When the training data cannot cover all the operating conditions, for example, the experiments are incomplete or the battery characteristics have changed, the SOC estimation accuracy will decrease.
The model-based methods are dominated by Kalman filter and its derivatives.Chen et al. [12] used an extended Kalman filter (EKF) along with a nonlinear battery model to estimate the SOC of the lithium-ion battery.In order to further improve the accuracy, adaptive extended Kalman filter (AEKF), unscented Kalman filter (UKF), and adaptive unscented Kalman filter (AUKF) are proposed.He et al. [13] identified the parameters of an improved Thevenin battery model using EKF, and then adopted an AEKF for obtaining correct and robust SOC of the lithium-ion batteries.In [14], a UKF-based method was used to self-adjust the model parameters and provide state of charge estimation of the battery.References [15,16] adopted AUKF to realize SOC estimation.In addition, an adaptive cubature Kalman filter (CKF) was proposed in [17] to improve the convergence rate and SOC estimation accuracy.To overcome the accuracy degradation caused by non-Gaussian noise, particle filters were utilized to estimate the SOC of batteries in [18,19].To compensate for the time-variability of battery parameters due to variational operating conditions and battery aging, a dual EKF and a dual UKF were employed in [20,21] for simultaneous SOC and parameters estimation, respectively.The model-based methods eliminate the need for correct initial SOC values and accurate measurement, which are requisite for the CC method, and have no demand for a large amount of train data, so it is widely applied and is studied extensively.However, the estimation accuracy is strongly affected by the battery model and parameters.Once there is a model mismatch, the estimation performance will rapidly decline.
In practical applications, the battery model parameters will change with SOC, temperature, and the degree of battery aging.Moreover, the statistic information of the process noise and sensor noise may be unknown or time-varying.Under these situations, the traditional KF-based methods will have low estimation accuracy and poor robustness.To address this issue, some researchers [22][23][24] resort to a H∞ filter, which takes the time-varying battery parameter into account and has no need to know the statistics of the process noise and the measurement noise.It has strong robustness under uncertain conditions.However, the H∞ filter is a tradeoff between robustness and accuracy, so the SOC estimation accuracy is sacrificed to some extent for the robustness.Liu et al. [25] combined the idea of square root filter and adaptive unscented Kalman filter (ASRUKF) algorithm based on improved Sage-Husa estimation to estimate the SOC of a lithium-ion battery.This method adaptively adjusted the values of the process and measurement covariances in the estimation process to improve the accuracy of SOC estimation.However, it does not consider the uncertainties brought by varying battery model parameters.EI Din et al. [26] proposed a multiple-model EKF (MM-EKF) and an autocovariance least squares (ALS) method for estimating the SOC of lithium-ion battery cells.MM-EKF reduced the dependence of the EKF algorithm on the correct assumptions of the measurement noise statistics by weighted summation of the estimates of multiple hypothesized EKFs.The ALS method extracted the possible correlation in the innovation sequence to estimate the measurement noise covariance.However, both methods leave out consideration of the time-varying battery parameters.Furthermore, the computation load is larger compared with the conventional EKF algorithm.
In fact, the Bayesian approach is the most general approach of solving the problem with uncertain parameters.But it is hard to get the analytical solution for most Bayesian approaches due to complex probability density function or high dimension of integration.Recently, the variational Bayesian (VB) inference method [27][28][29][30][31] has drawn extensive attentions, which utilizes a new simpler, analytically tractable distribution to approximate the true posterior distribution so as to avoid the direct complex calculation of multi-dimensional probability density function.Sarkka et al. [27] adopted the VB method for joint recursive estimation of the dynamic state and the time-varying measurement noise parameters in linear state space models.Li et al. [28] employed VB approximation for the unscented Kalman filter to estimate the time-varying measurement noise covariance so as to improve algorithm adaptability.Sun et al. [29] proposed a VB method to estimate the system states with unknown inputs.Hou et al. [30] combined the VB method with the shifted Rayleigh filter to jointly estimate the target position and the clutter probability so that improving the performance of bearings-only target tracking.In [31], the VB approach was applied to estimate the ARX model parameters along with time delays.
In this paper, we combine the idea of variational Bayesian inference with the dual EKF algorithm (VB-ADEKF) to jointly estimate the battery parameters and SOC of lithium-ion batteries of electric vehicles.Meanwhile, the measurement noise variances are simultaneously estimated with the state estimation to account for the battery model uncertainties and measurement uncertainties, minimizing the impact of model mismatch.The effectiveness of the proposed algorithm has been verified through experiments under different operating conditions.The results show that the proposed VB-ADEKF algorithm outperforms the dual EKF algorithm in terms of SOC estimation accuracy, convergence rate, and robustness.The mean SOC estimation error of VB-ADEKF is under 1% for most of the time.
The remainder of this paper is organized as follows: Section 2 describes the battery model, the definition of SOC and establishes the state space models for SOC estimation and battery parameter estimation.Section 3 illustrates the variational Bayesian approximation-based adaptive Kalman filter and Section 4 presents the proposed variational Bayesian approximation-based adaptive dual extended Kalman filter.The experimental verification and analysis are presented in Section 5. Finally, Section 6 provides a conclusion.

Battery Model
For the accurate estimation of SOC, a reliable battery model is required.Considering the model accuracy, the structure complexity, and the computation time, the first order resistor-capacitor (RC) model as shown in Figure 1 is adopted to model the lithium-ion battery.The electrical behavior of the model can be written as follows: where U t denotes the terminal voltage of the battery, U oc is the open-circuit voltage, U 1 is the polarization voltage of the RC network, I L is the load current, R 0 represents the ohmic internal resistance, and R 1 and C 1 represent the polarization resistance and polarization capacitance, respectively.The nonlinear relationship between the open-circuit voltage and the SOC is described using the fifth-order polynomial model as: where k 0 , k 1 , k 2 , k 3 , k 4 , k 5 are the parameters to be identified.

Definition of State of Charge
The state of charge (SOC) is defined as the ratio of the remaining capacity in a battery over the rated battery capacity.Using the CC method, the battery SOC can be calculated as follows: where η c is the coulomb efficiency, I L,t is the load current at time t.Q rate is the rated capacity of the battery.

State Space Model for SOC Estimation
Taking X = [SOC, U 1 ] T as the state vector, the load current I L as the input and the terminal voltage U t as the output, we can obtain the discrete state-space model as where T represents the time-varying battery model parameter vector; The initial state has a Gaussian prior distribution X 0 ∼ N ( X0 , P 0 ).The state prior and process noise are assumed to be known, while the measurement noise variance Σ x k is assumed to be unknown.In addition, the process noise, measurement noise, and initial state value are independent of each other.
f (•) and h(•) represent the nonlinear functions of state vector X k , input I L,k , and battery model parameter vector θ k .Their mathematical expressions are where ∆t is the sampling interval of the current.
The state transition matrix and the input control matrix are respectively Linearizing the measurement function, we can get the Jacobian measurement matrix as

State Space Model for Battery Parameter Estimation
Because the battery model parameters vary with changes in the batteries' SOC, degree of aging, and environmental temperature, online recursive battery parameter estimation is needed.So here we establish the state-space equations of the battery parameters as: where is a small white noise with covariance Q θ k that reflects the time-varying parts of the parameters, d k is a measurement function of θ k , and e k ∼ N (0, Σ θ k ) is the measurement noise to account for the sensor noise and modeling uncertainties.Σ θ k is assumed unknown here.The Jacobian measurement matrix of θ k is calculated as follows [20]: where

Variational Bayesian Approximation-Based Adaptive Kalman Filter Algorithm
As we know, when the measurement noise covariance is unknown or time varying, the classical approach of solving the problem is to use adaptive filters.The Bayesian approach is a typical choice.But it is usually hard to get the analytical form due to complex probability density function or high dimension of integration.In [27], variational Bayesian (VB) approximation was firstly used for the Kalman filter to estimate the joint posterior distribution of the states and noise covariances.
Given a discrete-time linear state space model as: where Assuming the state vector and the measurement noise covariance are independent, the joint posterior distribution of the state and covariance p(X k , Σ k |y 1:k ) can be solved by VB approximation as follows: The VB approximation can now be formed by minimizing the Kullback-Leibler (KL) divergence between the separable approximation and the true posterior: Minimizing the KL divergence with respect to the probability densities Q x (X k ) and Q Σ (Σ k ) in turn, while keeping the other fixed, we can get the following equations: Computing the above equations, we can get the following densities [27]: where the parameters m k , P k , α k,i , β k,i can be calculated as: where i = 1, ..., d and d denote the dimensionality of the measurement.m − k and P − k are the predicted state estimate and its covariance, respectively.Here we assume each component of the measurement noise variance is mutually independent, and then the covariance matrix is diagonal, estimated as: Furthermore, in order to describe the dynamics of the measurement noise variance, the inverse Gamma distribution parameters are assumed to change by a scale factor ρ ∈ (0, 1].The formulas are given as: Note that the value ρ = 1 corresponds to stationary variances and lower values represent larger time-fluctuations.
The above are the basic equations of adaptive Kalman filter based on VB approximation.Moreover, when the state equation and the measurement equation are nonlinear, the VB method can be rewritten in the EKF framework.

Variational Bayesian Approximation-Based Adaptive Dual Extended Kalman Filter
In order to handle the joint estimation of the SOC and the battery model parameters as well as the unknown measurement noise covariances, we propose a variational Bayesian approximation-based dual extended Kalman filter (VB-ADEKF) in this paper.
First, let us rewrite the state-space equations for SOC estimation as and the state-space equations for battery model parameters as where the process noise covariances Q x k , Q θ k are assumed known, and the measurement noise variances Σ x k , Σ θ k are unknown, being assumed as stochastic variables.
Then, VB-ADEKF is to alternatively solve the joint posterior distribution p(X k , Σ x k |y 1:k ) of the SOC and its measurement noise variance and the joint posterior distribution p(θ k , Σ θ k |y 1:k ) of the battery parameter and its measurement noise variance by VB approximation as follows: where ) can be solved by Equations ( 23) and (24).It should be noted that the parameters in these approximating densities are actually obtained using the extended Kalman filters since the measurement function h(•) is nonlinear.The filtering procedure of the VB-ADEKF algorithm is summarized in Algorithm 1.
(3) Update: the update of VB-ADEKF utilizes iterate filtering framework. x,k For n = 1 : N, iterate the following N (N denotes iterated times) steps: • Measurement variances: • State estimate and its covariance: • Battery parameters estimate and its covariance: • Parameters for the measurement noise variances estimation: β(n+1) x,k By alternatively using two VB-based extended Kalman filters for online estimation of the battery SOC and model parameters, while compensating for the uncertainties in the model parameters by simultaneous estimation of the measurement noise variances, the adaptability of the proposed algorithm to dynamic changes in battery characteristics is greatly improved.Hence, it is very promising to further increase the SOC estimation accuracy and robustness.

Experimental Settings
The proposed method is experimentally evaluated in this section.The experimental setup for the tests is shown in Figure 2. The setup consists of three lithium-ion battery cells connected in series, an electronic load, and a host computer.The tested lithium-ion battery cells are type 18,650, whose nominal capacity is 2200 mAh, nominal voltage is 3.7 V, charging and discharging cutoff voltages are 4.2 V and 3 V, respectively.The type of the electronic load is IT8516S produced by ITECH, whose current measurement accuracy is ±(0.1% + 0.1% full scale), and voltage measurement accuracy is ±(0.02% + 0.02% full scale).The load current, terminal voltage, and SOC can be recorded via the host computer during the discharge test.
First, a sequence of pulsed discharging experiments were implemented to determine the relationship of open circuit voltage (OCV) and SOC.By using the curve-fitting toolbox in MATLAB, k 0 , k 1 , k 2 , k 3 , k 4 , k 5 were identified, which are shown in Table 1.The measured data and fitted curve are presented in Figure 3.The R-square is used to represent the goodness of fit.The normal value range of the R-square is 0-1 and a value closer to 1 indicates a better fitting curve [32].It can be seen that the curve fits well with the measurement data, indicating that the selected fifth-order polynomial model can describe the relationship of OCV and SOC very well.Then, a constant current discharge test, a pulse current discharge test, and an urban dynamometer driving schedule (UDDS) test were performed.They are commonly used to verify the performances of SOC estimation methods in EVs.In the constant current discharge test, the current keeps invariant but the terminal voltage declines continuously.In the pulse current discharge test, the current stays at 1 A for 10 min and then decreases to 0 and lasts for 30 min.The process is repeated until the battery reaches the lower cut-off voltage.UDDS, also known as FTP72, was used as a test procedure to certify vehicle emissions by the US Environmental Protection Agency.It can simulate the actual driving conditions of vehicles on the road.The battery current and voltage are both sampled at 1 s.In each test, the true SOC was obtained using CC method.The estimation accuracy, convergence rate, and robustness of the proposed VB-ADEKF are evaluated by comparison with DEKF under different tests.

Constant Current Discharge Test
The experiment was performed with a constant discharge current of 1 A. The initial SOC value is set to be 0.8, rather than the real SOC of 1.0.The process noise covariances are set as The measurement noise variances used for DEKF are Σ x k = 0.001 and Σ θ k = 0.0005.The scale factors for VB-ADEKF are set to 1 × 10 −4 .The initial parameters α 0 and β 0 for battery parameters and SOC are both set as 10 and 0.001, respectively.
The estimated values of the battery model parameters by VB-ADEKF are presented in Figure 4.The ohmic resistance is stable at the beginning of the discharge and increases at the end of the discharge.The polarization resistance decreases at first and increases later with the depth of the discharge and again decreases in the last 1000 s.The polarization capacitance shows a declining trend overall in the process of discharge.Furthermore, to verify the accuracy of online identification of the battery parameters by the VB-ADEKF algorithm, the measured terminal voltage and the estimated terminal voltage by VB-ADEKF were compared, as shown in Figure 5.The maximum absolute error was 0.023 V, except for the first big error caused by an incorrect initial SOC value.The mean absolute error was 0.0050 V and the relative mean absolute error was 0.052%.It is clear that the estimated terminal voltage agrees well with the measured voltage.This illustrates the effectiveness of the battery model whose parameters are identified by the proposed VB-ADEKF method.The SOC estimation results using VB-ADEKF and DEKF are shown in Figure 6.Clearly, VB-ADEKF has a much more accurate SOC estimation than DEKF.The SOC estimation error of VB-ADEKF is bounded within ±1% for most of the time, but DEKF goes outside of this interval.The detailed error values are shown in Table 2.The maximum absolute estimation error of VB-ADEKF is 1.28% and the mean absolute error is 0.64% after discharge for 12 min.Both are significantly smaller than the errors of DEKF, which are correspondingly 2.76% and 1.39%.Meanwhile, from the figure it can be seen that VB-ADEKF converges much faster than DEKF under initial SOC errors.This is also verified in Table 2.The convergence time, which is defined as the first time instant at which the SOC estimation error decreases to ±5%, are 10 s and 335 s for VB-ADEKF and DEKF, respectively.

Pulse Current Discharge Test
In the pulse current discharge test, the initial SOC value is set to be 0.8 while the real SOC is 1.0.The process noise covariances are set as The measurement noise variances used for DEKF are Σ x k = 0.01 and Σ θ k = 0.0005.The initial parameter values of VB-ADEKF are set as α θ,0 = 100, β θ,0 = 0.0005, α x,0 = 10, β x,0 = 0.001.Figure 7 presents the estimated values of the battery model parameters.Clearly, there is a stepped increase in the ohmic resistance from 34 mΩ to 49 mΩ.The polarization resistance retains stable during the entire discharge process.The polarization capacitance first decreases rapidly then stays stable until the end of the discharge.Figure 8 shows the the measured terminal voltage and the estimated terminal voltage by VB-ADEKF to verify the effectiveness of the battery parameters identification.The maximum and mean absolute estimation errors are 0.086 V and 0.0027 V, respectively.It implies that the estimated terminal voltage has a good agreement with the measured voltage.It further shows the battery model parameters are well identified.The SOC estimation results of VB-ADEKF and DEKF are plotted in Figure 9.It shows that VB-ADEKF and DEKF have comparable performance in convergence rate.However VB-ADEKF has more accurate SOC estimation than that of DEKF.This is also verified in Table 2, in which the mean absolute error of VB-ADEKF is 0.68%, while it is 1.01% for DEKF after discharge for 12 min.

UDDS Test
To evaluate the SOC estimation performance under dynamic loading profiles, an UDDS cycle was performed on the battery cells.According to the actual tolerable currents of the lithium-ion battery cells, the loading currents are scaled down, as shown in Figure 10.The initial SOC value is set to 0.8.
The estimated values of the battery model parameters are presented in Figure 11.It can be seen that the range of parameter values of R 0 , R 1 and C 1 are consistent with these in constant current discharge test and pulse current discharge test, but the changing rules are slightly different.Figure 12 presents the the measured terminal voltage and the estimated terminal voltage by VB-ADEKF to verify the effectiveness of the battery parameters identification.The maximum and mean absolute estimation errors are 0.062 V and 0.0011 V, respectively.It is clear that the estimated terminal voltage has a good consistency with the measured voltage.It further shows the battery model parameters are well identified.Figure 13 presents the SOC estimation results of VB-ADEKF and DEKF.It shows that both VB-ADEKF and DEKF have good estimation accuracy when SOC is between 20% and 90%.But VB-ADEKF still outperforms the traditional DEKF in SOC estimation accuracy and convergence rate.A comparative summary of the two methods is given in Table 2.When SOC is decreased to 20%, the estimation error of both methods begins to increase, but is no larger than 5%.This may be because the polarization effect of the battery is further aggravated at lower SOC level.As a result, the terminal voltage measurement error has increased, thereby reducing the SOC estimation accuracy.

Convergence Ability with Initial SOC Error
Because it is difficult to determine the initial SOC value precisely in practical applications, it is important and indispensable for the algorithms to have the ability to correct the uncertainty brought by the initial SOC error.Therefore, the convergence rate with initial SOC error is adopted as another indicator for evaluating the SOC estimation algorithms.The true initial SOC value is 1.0.Figures 14-16 present the estimation results of VB-ADEKF and DEKF algorithms for different initial SOC values from 30% to 90% under the above three tests, respectively.Overall, the convergence time increases as the initial SOC error rises for both VB-ADEKF and DEKF.But the growth rates are very different for the two filters under different tests.Specifically, from the results of the convergence time shown in Table 3, it can be seen that the convergence time of VB-ADEKF is stabilized about 10 s into the constant current discharge test when the initial SOC value is larger than 50%, and increased to more than 100 s when the initial SOC values are 40% and 30%.But the convergence time of DEKF, which is between 300 s and 500 s, is much longer than VB-ADEKF.It shows that VB-ADEKF can quickly converge to the true SOC values without resulting in the accumulation of errors caused by the initial error of the SOC.In the pulse current discharge test, the convergence rates are comparable for VB-ADEKF and DEKF in the case of small initial SOC errors, for example, 10% or 20%.However, when the initial SOC error is relatively large, the convergence time of VB-ADEKF becomes much smaller than that of DEKF.In the UDDS test, VB-ADEKF only exhibits a slightly increasing trend with the increase of the initial SOC error.However, DEKF converges slower than VB-ADEKF, in the meantime, with a larger SOC estimation error.Its convergence time goes up quickly as the initial SOC error increases.This implies that the initial SOC error has a noticeable impact on the performance of DEKF.But, from the overall perspective, VB-ADEKF is not very sensitive to the initial SOC error.It shows that the proposed VB-ADEKF has better robustness for initial SOC errors than the traditional DEKF.

Effect of Mistuning
If the working condition of the battery changes abruptly, the SOC measurement error would probably be varied largely with before, so the prior tuning of the measurement variance of the DEKF will not give an optimal estimate of the SOC.But there is no such issue in the proposed VB-ADEKF since the measurement variances are estimated online.This effects of mistuning brought about by inappropriate measurement variance of the traditional DEKF and VB-ADEKF are mainly reflected in the SOC estimation error, as shown in Figures 17-19.Here, the measurement noise variances of DEKF are mistuned to Σ x k = 0.01 and Σ θ k = 0.001 in constant current discharge test, and Σ x k = 0.1 and Σ θ k = 0.005 in the pulse current discharge test and UDDS test.The initial estimates of the measurement variances of VB-ADEKF are also correspondingly mistuned.From the results, we can see that the convergence rate of DEKF slows down and the SOC estimation accuracy also declines in the case of mistuning, while the SOC estimation performance of the proposed VB-ADEKF remains almost unchanged.This suggests that the proposed VB-ADEKF is more robust than the traditional DEKF.

Conclusions
To deal with the measurement statistical uncertainties and inaccurate battery model, a variational Bayesian approximation-based adaptive dual extended Kalman filter (VB-ADEKF) is proposed in this paper for SOC estimation of lithium-ion batteries.First, the variational Bayesian inference is integrated with the extended Kalman filter to jointly estimate the states and the measurement noise covariances.Then, two VB-based extended Kalman filters are alternatively used for online estimation of the battery SOC and model parameters, while simultaneously estimating the measurement noise variances to compensate for the uncertainties in the measurement and battery parameters.Therefore,

Figure 3 .
Figure 3.The relationship curve of open circuit voltage (OCV) versus state of charge (SOC).

Figure 4 .
Figure 4. Results of parameter identification using variational Bayesian approximation-based adaptive dual extended Kalman filter (VB-ADEKF) in the constant current discharge test.

Figure 5 .
Figure 5. Experimental terminal voltage results using VB-ADEKF in the constant current discharge test.

Figure 6 .
Figure 6.SOC estimation results using VB-ADEKF and dual extended Kalman filter (DEKF) in the constant current discharge test.

Figure 7 .
Figure 7. Results of parameters identification using VB-ADEKF in the pulse current discharge test.

Figure 8 .
Figure 8. Experimental terminal voltage results using VB-ADEKF in the pulse current discharge test.

Figure 9 .
Figure 9. SOC estimation results using VB-ADEKF and DEKF in the pulse current discharge test.

Figure 12 .
Figure 12.Experimental terminal voltage results using VB-ADEKF in the UDDS test.

Figure 13 .
Figure 13.SOC estimation results using VB-ADEKF and DEKF in the UDDS test.

Figure 16 .
Figure 16.SOC estimation results of the two filters with different initial SOC values in the UDDS test.

Figure 17 .
Figure 17.SOC estimation results using VB-ADEKF and DEKF in the case of mistuning in the constant current discharge test.

Figure 18 .Figure 19 .
Figure 18.SOC estimation results using VB-ADEKF and DEKF in the case of mistuning in the pulse current discharge test.

Table 2 .
Statistical analysis of state of charge (SOC) estimation error (after 12 min) in three tests.

Table 3 .
Convergence time (s) of variational Bayesian approximation-based adaptive dual extended Kalman filter (VB-ADEKF) and dual extended Kalman filter (DEKF) with different initial SOC values.