Next Article in Journal
On the Mathematical Modeling of Line-Start Permanent Magnet Synchronous Motors under Static Eccentricity
Next Article in Special Issue
Development of Automotive Permanent Magnet Alternator with Fully Controlled AC/DC Converter
Previous Article in Journal
Coordinated Engine-Start Control of Single-Motor P2 Hybrid Electric Vehicles with Respect to Different Driving Situations
Previous Article in Special Issue
Stochastic Model Predictive Fault Tolerant Control Based on Conditional Value at Risk for Wind Energy Conversion System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Method for State of Charge Estimation of Lithium-Ion Batteries Using Square Root Cubature Kalman Filter

1
State Key Laboratory of Advanced Design and Manufacture for Vehicle Body, Hunan University, Changsha 410082, China
2
Haima Automobile Group Co., Ltd., Haikou 570216, China
3
Hubei Key Laboratory of Advanced Technology for Automotive Components, Wuhan University of Technology, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
Submission received: 22 December 2017 / Revised: 2 January 2018 / Accepted: 8 January 2018 / Published: 15 January 2018

Abstract

:
State of charge (SOC) is a key parameter for lithium-ion battery management systems. The square root cubature Kalman filter (SRCKF) algorithm has been developed to estimate the SOC of batteries. SRCKF calculates 2n points that have the same weights according to cubature transform to approximate the mean of state variables. After these points are propagated by nonlinear functions, the mean and the variance of the capture can achieve third-order precision of the real values of the nonlinear functions. SRCKF directly propagates and updates the square root of the state covariance matrix in the form of Cholesky decomposition, guarantees the nonnegative quality of the covariance matrix, and avoids the divergence of the filter. Simulink models and the test bench of extended Kalman filter (EKF), Unscented Kalman filter (UKF), cubature Kalman filter (CKF) and SRCKF are built. Three experiments have been carried out to evaluate the performances of the proposed methods. The results of the comparison of accuracy, robustness, and convergence rate with EKF, UKF, CKF and SRCKF are presented. Compared with the traditional EKF, UKF and CKF algorithms, the SRCKF algorithm is found to yield better SOC estimation accuracy, higher robustness and better convergence rate.

Graphical Abstract

1. Introduction

Increasingly along with the fossil energy depletion, air pollution and more and more serious global climate changes, people have begun to realize the great importance of the utilization and development of non-fossil energy [1]. Governments all over the world have introduced a variety of incentives to reduce pollutant emissions and greenhouse gas emissions [2,3,4,5,6]. Since transportation consumes a large amount of energy, it is necessary to develop and utilize electric vehicles (EVs) to realize green mobility. Lithium-ion batteries (LIBs) have many advantages, such as high power density and durability, and therefore they have been widely used in EVs. However, when overcharging occurs, the LIBs are more likely to burn and explode than other batteries, necessitating higher requirements for battery management system (BMS) [7,8,9,10]. A few parameters, such as maximum discharge current, specific energy of weight, specific energy of volume, and power density, determine an EV’s performance [11,12,13,14,15]. The most important index for a BMS is state of charge (SOC). Because of the inherently time-varying and non-linearity characteristics of LIB under working conditions, accurate estimation of SOC remains a challenge [16,17,18,19].

1.1. Review of Methods

Some popular methods for SOC estimation are discussed below.
The open-circuit voltage method is simple and easy to use. SOC estimation in the late and early stages of charging is effective, but batteries must be tested online for a long time. This approach is used to estimate the charging efficiency of battery packs in laboratories and the accuracy of SOC estimation [10,19,20,21].
The neural network model method can realize basic nonlinear characteristics and can give corresponding outputs to external excitation. This method can also be used to simulate the characteristics of the nonlinear dynamic of batteries. The main problem with this method is that a large amount of reference data are required to be trained, and the estimation error of SOC is significantly impacted by the approach and data. Other flaws are that this method consumes considerable computer time in BMS applications, high demands for the BMS design, and high costs to adopt high-performance chips, which obviously cannot meet the requirements in the competitive automobile industry; consequently, the method cannot be used in mass production of electric vehicles [22].
Particle filter (PF) [23] and unscented PF [24,25] perform efficiently in nonlinear and non-Gaussian estimation. However, the major issue is that a complex BMS environment requires considerable sample data to describe posterior probability distribution. Other flaws are similar to those of neural network approaches [26,27,28,29,30].
The extended Kalman filter (EKF) is a typical representative of traditional nonlinear estimation approaches. The basic idea is to locally linearize a nonlinear state and the functions of measurement to perform a first-order Taylor polynomial expansion and then apply the Kalman filter formula of linear systems. EKF requires that the state and measurement functions of nonlinear systems be continuously differentiable, which limits its application scope. The first-order linearization approximation of nonlinear functions leads to low estimation accuracy, especially when systems are strongly nonlinear and the EKF estimation precision is severely reduced or even divergent. The Jacobin matrix of nonlinear functions should be calculated, which can easily cause differential stability of EKF and calculation divergence [31,32,33,34,35,36].
Unscented Kalman filter (UKF) [37] differs from EKF in that it replaces the first-order Taylor expansion in EKF with unscented transform (UT). UKF selects 2n + 1 sigma points with weight according to UT to approximate the mean of state variables. After the sigma points are propagated by nonlinear functions, the mean and the variance of the capture can achieve third-order precision of the real values of the nonlinear functions; thus, the precision of UKF is higher than that of EKF [38]. Nevertheless, UKF needs to compute and transform the sigma points, which leads to a high computational complexity and cost of BMS.
In cubature Kalman filter (CKF), the mean and the variance of the system state are propagated through 2n equal-weight cubature points, which enables it to obtain high estimation accuracy. CKF, on the basis of cubature criteria, gives new points after the transformation of nonlinear system equations from 2n cubature points that have equal weight values to give the estimation of the system state at the next time. The cubature point and its weight are only determined by the state dimension and can be calculated and stored in advance. CKF can reduce computation time while maintaining estimation accuracy because of its use of minimal sampling points, and it is gradually replacing UKF as the filtering method for embedded system applications [39]. The CKF algorithm in engineering applications is susceptible to abnormal perturbation, inexact initial values, and difficulty in decomposing non-semi-definite matrices with Cholesky, which leads to system divergence [40].
In addition to the above methods, data-driven approaches, such as sliding mode observer [41], H-infinity observer [42], recursive further squares [43], EKF-based methods with online parameter identification [44,45], the Lyapunov-based estimator [46], and an estimation of the SOC and the model parameters simultaneously with paralleled filters [47,48,49,50] have also been studied. Most of these methods require adequate incentives. If the incentives are not obvious (such as charging condition), these methods are not very good for SOC estimation. These methods can solve some particular problems with certain requirements, but these methods cannot be good for the robustness [13,51] and anti-measurement noise interference. In addition, there are other methods based on electrochemical models [52,53] that have been reported. However, because of the complexity (strict electrode dynamic theoretical framework and a large number of battery data accumulation) of these methods, they are not suitable for real BMS.
In real conditions, most controllers use embedded micro-controller units to compute the algorithm and control logic. Under the finite computing resources of embedded micro-controller units, the method for SOC estimation must not be complex and must have a rapid convergence speed and anti-measurement of noise interference to increase the robustness of BMS. Through what has been reviewed above, some requirements have not been fully met. We proposed a new method for SOC estimation using the square root cubature Kalman filter (SRCKF) algorithm in this paper. The SRCKF is based on CKF, directly propagates and updates the square root of the state covariance matrix in the form of Cholesky decomposition, obtains higher computational efficiency, guarantees the non-negative qualitative of the covariance matrix, avoids the divergence of the algorithm, and improves the convergence speed and stability of the filter.

1.2. Contribution of This Study

This paper aims to develop a new method based on the Thevenin model. The new method can improve the accuracy, robustness, and convergence rate of SOC estimation. The SRCKF algorithm is different from EKF and UKF algorithms, uses the cubature principle for nonlinear state particles and presents robustness in extreme nonlinear conditions. Compared with the EKF method, the present method can enhance SOC estimation accuracy and robustness and improve the convergence rate. When compared with the UKF method, the computation cost can be reduced with better estimation accuracy of SOC and robustness. Compared with the CKF approach, the present method can enhance the robustness and convergence rate with the same SOC estimation. Furthermore, the SRCKF algorithm has a better capability to overcome disturbances of unstable noises for SOC estimation. Thus, it can ensure the robustness and estimation accuracy of the algorithm in the case of limited computing resource requirements and can be applied to new BMS with online SOC estimation easily.

1.3. Organization of This Paper

The remaining sections of this paper are structured as follows. Section 2 discusses the parameters of the dynamic LIB model. Section 3 elucidates the principles of the SRCKF algorithm of SOC estimation. Section 4 shows the experimental configurations and the results of the parameters of the lithium-ion battery model identified offline. Experimental results and discussion are presented in Section 5, followed by conclusions in Section 6.

2. Model of the Battery

Lumped-Parameter Model of the Battery

An equivalent circuit model is developed on the basis of the dynamic features and operation mechanism of batteries [54]. The result indicates that 1RC and 2RC network-based lumped-parameter models of battery have good performance in terms of their estimation precision and model complexity [8]. In this paper, the 1RC network-based lumped-parameter model (Thevenin model) of battery, shown in Figure 1, is selected in consideration of calculation accuracy, computational and parameter identification complexities.
An equivalent circuit model has these components: the open circuit voltage U o c v represents the voltage source and describes the static characteristic of the battery, an RC network consists of polarization resistance R p and polarization capacitance C p describes the polarization phenomenon of the battery, resistance R 0 represents the internal resistance of the battery.
According to the analysis of the electric circuit [20], the electrical characteristics of the model of battery can be described as:
U ˙ p = U p / C p R p + i / C p
U t = U o c v U p U 0
{ U 0 = i R d c h g ( d i s c h a r g e ) U 0 = i R c h g ( c h a r g e )
where the load current is represented with i (charge describes negative, discharge describes positive); the open-circuit voltage (OCV) is represented with U o c v ; U t and U p respectively denote the terminal voltage and the polarization voltage; R p and C p represent the polarization resistance and capacitance of the RC network, respectively; U 0 represents the internal resister voltage; R c h g and R d c h g denote the internal resistance when the battery is charging or discharging, respectively.
The SOC of a battery denotes the ratio of the remaining capacity to the total capacity. After the current integral and discretization, SOC can be obtained using the following equation:
SOC ( t ) = SOC ( t 0 ) + ( η t 0 t i ( t ) d t ) / C N
where SOC ( t ) is the SOC of the battery at time ( t ), SOC ( t 0 ) is the SOC of the battery at the previous point ( t 0 ), η k is the Coulomb efficiency, and C N is the battery-rated capacity.
After the discretization of (2), we can obtain the following formula:
U p , k U p , k 1 / T s = U p , k 1 / C p R p + I k 1 / C p
Subsequently, we can acquire the following formula:
U p , k = ( 1 T s / C p R p ) U p , k 1 + T s / C p I k 1
where U p , k denotes the U p at time ( k ), U p , k 1 denotes the U p at time ( k 1 ), T s represents the sample time, and I k 1 represents the current at time ( k 1 ).
U p and SOC are selected as system state variables, i is selected as an input variable, and U t is selected as an output variable on the basis of Formulas (1)–(4) and (6). Then, the following formulas can be obtained:
( U p , k SOC k ) = ( ( 1 T s / C p R p ) 0 0 1 ) ( U p , k 1 SOC k 1 ) + ( T s / C p η k T s / C N ) ( I k 1 ) + ( ω 1 , k 1 ω 2 , k 2 )
{ ( U t , k ) = ( 1 0 ) T ( U P , k SOC k ) + ( R d c h g ) ( I k ) + U O C V + ( ν k ) ( d i s c h a r g e ) ( U t , k ) = ( 1 0 ) T ( U P , k SOC k ) ( R c h g ) ( I k ) + U O C V + ( ν k ) ( c h a r g e )
Formula (7) describes the state equation of the Thevenin model, and Formula (8) describes the observation equation of the Thevenin model. At time ( k ), I k , U t , k and ν k denote the current, terminal voltage and measurement noise, respectively, and ω 1 , k 1 , ω 2 , k 2 represent the process noise.
Define:
x k = ( U p , k SOC k ) , A k = ( ( 1 T s / C p R p ) 0 0 1 ) , B k = ( T s / C p η k T s / C N )
The equation of time state is:
f ( x k , u k ) = A k x k + B k I k
With battery voltage as an output variable, we can generate the following observation equation:
g ( x k , u k ) = U o c v + U p , k + R d c h g I k ( d i s c h a r g e ) g ( x k , u k ) = U o c v + U p , k R c h g I k ( c h a r g e )

3. Cubature Integral Approximation Method and SRCKF Algorithm

3.1. Cubature Using Numerical Integral Approximation Method

The numerical method of calculating cubature is introduced before elaborating the SRCKF algorithm, considering that SRCKF is based on the cubature integral approximation method. We take the n-dimensional integral of the following Gaussian density function [40]:
I ( f ) = R n f ( x ) exp ( x T x ) d x
where f ( x ) denotes an arbitrary function and R n denotes the integral region. The key of the Bayesian theory to solve nonlinear Gaussian filtering is the integral of Formula (11). The integral term f ( x ) exp ( x T x ) can be summarized as multidimensional integrals with the form of “nonlinear function × Gaussian probability density”.
We define x = r y , y T y = 1 , r [ 0 , ) , and Formula (11) can be expressed as:
I ( f ) = 0 U n f ( r y ) r n 1 exp ( r 2 ) d σ ( y ) d r
U n = { y R n | y T y = 1 } is a supersphere with a radius of 1, and σ ( · ) is a spherical metric unit. Formula (12) can be simplified as
I ( f ) = 0 S ( r ) r n 1 exp ( r 2 ) d r
S ( r ) = U n f ( r y ) d σ ( y )
Therefore, the n-dimensional integral Formula (11) can be transformed into a spherical radial integral form, as shown in Formulas (13) and (14). The spherical integral principle and the radial principle can be used to solve this equation.
U n f ( y ) d σ ( y ) ω i = 1 2 n f [ u ] i
a b f ( x ) ω ( x ) d x ω i i = 1 m f ( x i )
Formula (15) shows a third-order spherical integral structure under the spherical integral principle, and [ u ] i is the i element of the generated operator. Formula (16) presents the radial principle of the sphere, which means that m point Gaussian integral can be equivalent to the sum of ( 2 m 1 ) polynomials. ω ( x ) is a nonnegative weight function in the integral interval [ a , b ] .
SRCKF uses the spherical radial principle to select 2n (n is state dimension) with corresponding weight point set ( ω i , ζ i ) , which can be used to approximate the posterior mean value and covariance of the nonlinear state.
I ( f ) = i = 1 m ω i f ( ζ i )
{ ζ i = n [ 1 ] i , ω i = 1 2 n , i = 1 , 2 , , 2 n .
where [ 1 ] i represents the column i of set [ 1 ] . We define n = 2 and obtain
[ 1 ] = { [ 1 , 0 ] T , [ 0 , 1 ] T , [ 1 , 0 ] T , [ 0 , 1 ] T }

3.2. SRCKF Algorithm

We define n-dimensional nonlinear discrete-state space mode and acquire
x k = f ( x k 1 , u k 1 ) + ω k 1 , k = 0 , 1 ,
z k = h ( x k , u k ) + υ k , k = 0 , 1 ,
where u k is the known input, ω k 1 is the system noise, and υ k is the noise of measurement.
The following three steps describe the SRCKF algorithm:
  • Step 1. Initialization of state estimation. We set the initial value of state to x ^ 0 | 0 and obtain
    { x ^ 0 | 0 = E ( x 0 ) S 0 | 0 = c h o l ( E [ ( x 0 x ^ 0 | 0 ) ( x 0 x ^ 0 | 0 ) T ] )
    where U = c h o l ( A ) is the Cholesky decomposition of matrix A , A = U T U , U is the upper triangular matrix, and E ( ) is the expectation.
  • Step 2. State prediction ( k = 1 , 2 , 3 , ).
    • Step 2.1 n cubature points for battery states at time k 1 ( i = 1 ,   2 , ,   2 n ) are calculated.
      X i , k 1 | k 1 = S k 1 | k 1 ξ i + x ^ k 1 | k 1
    • Step 2.2: The battery state Formula (7) is used to propagate cubature points and generate new points.
      X i , k | k 1 = f ( X i , k 1 | k 1 , u k 1 )
    • Step 2.3: After the newly generated points are weighted, the sum of the weighted points is determined and then the prediction value of the battery state at time k can be estimated (SRCKF uses equal weights).
      x ^ k | k 1 = 1 / 2 n i = 1 2 n X i , k | k 1
    • Step 2.4: The square root of the covariance matrix of the prediction value is estimated.
      S k | k 1 = T r i a ( [ χ k | k 1 S Q , k 1 ] )
      where Q k 1 = S Q , k 1 S Q , k 1 T .
      χ k | k 1 = 1 / 2 n [ X 1 , k 1 x ^ k | k 1 , X 2 , k 1 x ^ k | k 1 , , X 2 n , k 1 x ^ k | k 1 ]
S = T r i a ( A ) represents the QR decomposition of matrix A , which obtains a normal orthogonal matrix B and an upper triangular matrix C . We define S = C T , and S is the upper triangular matrix.
Step 3. Measurement update.
After the state prediction value of time k is obtained, the prediction value can be updated with the observation value at this moment to obtain the optimal state estimation value.
  • Step 3.1: A set of equal-weight cubature points by the square root matrix of the state prediction value and its covariance matrix is generated using the spherical radial rule ( i = 1 ,   2 , ,   2 n ).
    X i , k | k 1 = S k | k 1 ξ i + x ^ k | k 1
  • Step 3.2: The battery observation Formula (8) is used to propagate the cubature points.
    Z i , k | k 1 = h ( X i , k | k 1 , u k )
  • Step 3.3: The prediction value of observation at time k is generated.
    z ^ k | k 1 = 1 / 2 n i = 1 2 n Z i , k | k 1
  • Step 3.4: The square root of the covariance matrix of the estimated prediction value of observation is estimated.
    S z z , k | k 1 = T r i a ( [ Z k | k 1    S R , k ] )
    where R k = S R . k S R , k T , and
    Z k | k 1 = 1 / 2 n [ Z 1 , k | k 1 z ^ k | k 1 , Z 2 , k | k 1 z ^ k | k 1 , , Z 2 n , k | k 1 z ^ k | k 1 ]
  • Step 3.5: The observation and prediction values of the square root of the covariance matrix of each other are estimated.
    P x z , k | k 1 = χ k | k 1 Z k | k 1 T
    where
    χ k | k 1 = 1 / 2 n [ X 1 , k | k 1 x ^ k | k 1 , X 2 , k | k 1 x ^ k | k 1 , , X 2 n , k | k 1 x ^ k | k 1 ]
  • Step 3.6: The SRCKF filter gain matrix is solved.
    W k = ( P x z , k | k 1 / S z z , k | k 1 T ) / S z z , k | k 1
  • Step 3.7: The best state value estimate of the battery at time k is calculated.
    x ^ k | k = x ^ k | k 1 + W k ( z k z ^ k | k 1 )
  • Step 3.8: The square root of the error covariance matrix of the optimal state estimation of the battery at time k is calculated.
    S k | k = T r i a ( [ χ k | k 1 W k Z k | k 1    W k S R , k ] )
After the value of the state estimation x ^ k | k and its covariance square root S k | k are obtained at time k , Step 2 is repeated at time k + 1 , and the next iteration is performed. Flowchart of the method using SRCKF is shown in Figure 2.

4. Experimental Configurations, HPPC Test and Off-Line Identification of Relevant Parameters

According to the algorithm model, a test bench for the cells was built, which included cells of the battery, charging and discharging equipment (NBT BTS5200C4, Newaresles Limited, Shenzhen, China), monitoring platform, environment oven, testing cable, the hardware of the battery management unit (BMU) and local electronic control unit (LECU), and communication cables. The test bench can not only conduct Hybrid pulse power characterization (HPPC) tests to identify the relevant parameters of the battery off-line but also it can be used for simulation and verification according to the condition of a real environment.
HPPC tests must be conducted at various temperatures and different state of health (SOH) to identify the parameters of U o c v , R d c h g ( R c h g ) , C p , R p , and SOC and input them into the table to be used by the SRCKF algorithm to meet the requirements of the vehicle at various temperatures (−20 °C to 55 °C every 5 °C) and SOH (1 to 0.8 every 0.05). The specifications of the battery cell are shown in Table 1, the test equipment specifications are shown in Table 2, and Figure 3 shows the configuration of the battery test bench.
It is noted that, because of the high accuracy of NBT BTS5200C4, the reference values of SOC can be acquired by the current integral method with the measured value. To get a precise initial SOC, the battery cell is fully charged before test with a constant current, firstly. And to obtain an accurate terminal SOC, the battery cell is fully discharged after test, finally.
The HPPC test cycle of the battery cell is shown in Figure 4. There are various methods to fit functions in the MATLAB (R2016b, The MathWorks, Inc., Natick, MA, USA) simulation environment, such as polyfit command, is curvefit command, and cftool [55]. Among these methods, cftool is chosen to fit functions in this study, because it provides a visual graphical interface and powerful curve-fitting ability. The least square method [56] is used on the basis of the Formulas (1) and (2) to identify unknown parameters. The fitting curves of U o c v , R d c h g , C p , R p , and SOC equations are provided in Figure 5 (Temperature = 25 °C, SOH = 1) and Figure 6 (Temperature = 15 °C, SOH = 0.9). The coefficients of the corresponding equation ( y = f ( x ) = a x 6 + b x 5 + c x 4 + d x 3 + e x 2 + f x 1 + g ) are indicated in Table 3 and Table 4.

5. Results and Discussion

Under the same conditions, three different experiments were performed to test the performance of the proposed algorithm: constant current discharge test, dynamic stress test (DST), and DST test with noise measurement. At a temperature of 25 °C and SOH = 1, by using Figure 3 as the test bench, and the results of the battery parameters shown in Table 2, the EKF, UKF, CKF, and SRCKF algorithms were compared in the experiment. The comparison results were given in the performances of estimation accuracy, convergence speed, robustness, and calculation cost. Moreover, under real conditions, the actual current and voltage showed deviation from measured results due to the measurement errors of the voltage and current sensors and interference of corresponding hardware circuits (such as EMI). To verify the proposed algorithm robustness against the noise of measurement, a random interference to the system was conducted. Simulations of 100 cycles of the four algorithms were run in the Intel Core i7-4600U (2.1 GHz, 2.7 GHz, Intel Corporation, Santa Clara, CA, USA) to obtain the calculation cost of the algorithms. The test results are shown in Table 5. The comparison results showed that the proposed method had lesser computation cost than the other methods. Hence, because of the difference of calculation cost of the four algorithms, a cycle of the algorithms was run every 100 ms with BMU to ensure that the cycle was fully completed.

5.1. Experiment A: Test with Constant Discharge Current

In this experiment, with a constant current of 34.95 A (approximately 1 C), the cell of the battery was discharged from the initial state (SOC = 1) to the final state (SOC = 0.00027). The voltage of terminal decreased from 4.14 V to 2.5105 V. The curve of the test is shown in Figure 7. The current of the battery is constant with the car running in a constant speed. Thus, this experiment can be used to simulate this condition. The experimental results of convergence speed are shown in Figure 8. The cyan line is the reference value of SOC calculated by the BTS5200C4 with a precise initial value of SOC, and the estimated SOC is described with the red line. The initial value of real SOC is 1, and the estimated one is 0.5. A comparison of results from the EKF, UKF, CKF, and SRCKF estimation algorithms are shown in Figure 8a–d, respectively. Figure 9 shows the estimation SOC errors of the EKF, UKF, CKF, and SRCKF estimation algorithms, and the results of probability of estimation SOC errors are shown in Figure 10. Table 6 shows the data of the root mean square errors (RMSE), mean errors, maximum errors, and convergence time. The estimation SOC errors of the EKF, UKF, CKF, and SRCKF estimation algorithms are shown in Figure 9a–d, respectively. Similarly, the results of the probability of estimation SOC errors of the EKF, UKF, CKF, and SRCKF estimation algorithms are shown in Figure 10a–d, respectively. From the experiment results of Experiment A, the following conclusions are obtained. Under constant current discharge, (1) the proposed method had a faster convergence rate in the case of a higher accuracy of SOC estimation than the other three methods; (2) the estimation SOC error of the proposed method was distributed in the small region near zero, which showed good stability of the SRCKF algorithm; and (3) the proposed method met the requirement of the constant current discharge condition.

5.2. Experiment B: DST

The DST is simplified from actual urban driving cycles (AUDC). The current of battery with the DST in Experiment B is shown in Figure 11, with battery charge denoted with the positive current and battery discharge denoted with the negative current. Figure 11a shows one test cycle of DST and Figure 11b shows 20 test cycles of DST that were used in Experiment B. Experimental results of the convergence speed are shown in Figure 12. The cyan line is the reference value of SOC calculated by the BTS5200C4 with a precise initial value of SOC, and the estimated SOC is plotted by the red line. The initial value of real SOC was 1, and the estimated one was 0.5. The results of the EKF, UKF, CKF, and SRCKF estimation algorithms are shown in Figure 12a–d, respectively. To maintain proper estimation accuracy and convergence speed, the parameters of the EKF were adjusted to avoid the divergence of the EKF algorithm in Experiment B. However, the EKF algorithm still had the problem of divergence and filtering interruption in Experiment B. The curve in the purple circles represents this problem in Figure 12a. Figure 13 shows the estimation SOC errors of the EKF, UKF, CKF, and SRCKF estimation algorithms, and the results of probability of estimation SOC errors are shown in Figure 14. Table 7 shows the data of the root mean square errors (RMSE), mean errors, maximum errors, and convergence time. In comparison with Experiment A, the accuracy of the SOC estimation of the EKF, UKF, CKF, and SRCKF estimation algorithms decreases because the voltage and current have strong nonlinear changes. At the same time, because of the adjustment of the parameters of the EKF algorithm, the convergence speed time is long. Hence, the following conclusions are obtained from Experiment B. Under the DST condition, (1) the SRCKF algorithm had a faster convergence rate and better accuracy of SOC estimation than the other three methods; (2) the estimation error of the SRCKF algorithm was distributed in the small region near zero, which showed the good stability of the SRCKF algorithm; (3) the SRCKF algorithm met the requirement of the DST condition.

5.3. Experiment C: DST with Noise Measurement

In Experiment C, the conditions were the same as Experiment B, as shown in Figure 11. To verify the robustness of the proposed algorithm against the voltage and current errors, stochastic normal distributed noises were added to the measurement voltage and current. The standard deviation of the noise was calculated as follows:
θ = α N max / 3
where θ denotes the arithmetic mean, α denotes a scaling factor, and N max denotes the maximum voltage or the maximum current. Conditions of α = 1%, 2.5%, 5% are studied in this paper. The results of the current and voltage in Experiment C, with both voltage and current noises, are shown in Figure 15. Similar to Experiment B, the initial value of real SOC was 1, and the estimated one was 0.5. Because of the noise, the methods of EKF, UKF, and CKF diverged. The filters were interrupted in approximately 120 s. The curve in blue circle represents this problem in Figure 16. Evidently, the proposed method worked well. Hence, the RMSEs, mean errors, maximum errors, and convergence time of EKF, UKF, and CKF were not obtained in Experiment C. They were obtained in the SRCKF method with the different noises. The results are shown in Table 8, and the results of the probability of the corresponding estimation SOC errors are shown in Figure 17. In the real EVS, the error of current and voltage was not more than 2.5% after the signals were filtered through the hardware and software. As shown in Table 8, in the case of error ≤2.5%, the proposed method still maintained a high estimation accuracy, and its RMSEs and mean error values were still small enough to be applied to the actual EVS. Although the maximum error of the estimation SOC was bound within 8%, the RMSEs and mean errors of the proposed method were still small. From the point of error distribution shown in Figure 17d, the majority of the error distribution is within 5%. The conclusions of Experiment C, in the DST condition with noise (≥1%), are the following: (1) The proposed method had a faster convergence rate (approximately 25 s) and better accuracy of estimation SOC with noise than the other three methods; (2) the estimation error of the SRCKF algorithm was distributed in the small region near zero, which showed the good stability of the proposed method; (3) the proposed method had a weaker robustness against the disturbance than the methods of EK, UKF, and CKF and met the requirement of the DST condition with noises.

6. Conclusions

In this paper, a new SOC estimation method using the SRCKF algorithm is proposed to simultaneously improve the SOC estimation accuracy and reduce the computation cost. Under strong nonlinear working conditions, filtering divergence can be avoided. Some findings from the present study are summarized as follows.
(1)
The shortcomings of several common algorithms for battery SOC are investigated, and the SRCKF algorithm is recommended for its accuracy and speed. To balance between model accuracy and computation cost, the Thevenin model was applied to simulate the dynamic characteristics of the LIB, based on which Equations (9) and (10) were derived. With this model, experimental and parameter identification methods for off-line identification of battery parameters were proposed, and four Simulink algorithm models were established according to the EKF, UKF, CKF, and SRCKF algorithms, respectively.
(2)
Constant current discharge tests and DST cycles were conducted to verify the performance of the proposed method against other three methods (i.e., EKF, UKF, and CKF). The results of experiments showed that the estimation SOC error with the proposed method quickly converged to 2% within approximately 25 s, whereas the initial SOC error reached 50%. The RMSEs (0.00892), maximum error (0.02469), and mean error (−0.00704) of the proposed method were the best among all the methods in DST. Moreover, the maximum error was less than 5.5% even when the measurement noise of the voltage and the current was up to 2.5% in DST with noise.
(3)
Given the insufficient sensor accuracy and inaccurate hardware circuit, the signal of voltage and current are often interfered. Hence, a real EV running in a complex working condition (such as DST) with noises is very important. The method should be able to handle the nonlinear changes and random disturbance. In comparison with other methods (EKF, UKF, and CKF), only the proposed method avoided the filtering divergence in DST with noise (up to 5%).
Therefore, the proposed method ensured the robustness and estimation accuracy of the algorithm in the case of limited computing resource requirements and can be applied to new BMS with online SOC estimation easily. In the future, a new BMS with the proposed method to achieve good performance can be developed.

Acknowledgments

This work was sponsored by China’s National “863” Project of the Twelfth Five-year Plan “Research and industrialization of pure electric vehicles” (No. 2011AA11A215). The authors would also like to thank the Wuhan University of Technology for the support.

Author Contributions

Xiangyu Cui planned the initial aide. Xiangyu Cui and Yazhou Guo designed the algorithm. Xiangyu Cui analyzed the experimental date and wrote the original manuscript. Maji Luo and Zhu Jing and Huimin Qiao established the experiment bench. Xiangyu Cui and Zhu Jing and Maji Luo revised the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

iload current
Uocvopen-circuit voltage
Utterminal voltage
Uppolarization voltage
Rppolarization resistance
Cppolarization capacitance
U0internal resister voltage
Rchginternal resistance when discharging
Rdchginternal resistance when discharging
Rnintegral region
ukknown input
E(·)expectation
Nmaxmax voltage or the max current

Greek Symbol

σ(·)spherical metric unit
ω(x)nonnegative weight function
(ωii)weight point set
ωk−1system noise
νknoise of measurement
θarithmetic mean
αscaling factor

Acronyms and Abbreviations

OCVopen-circuit voltage
SOCstate of charge
SRCKFsquare root cubature Kalman filter
EKFExtended Kalman filter
UKFunscented Kalman filter
CKFcubature Kalman filter
EVelectric vehicle
LIBlithium-ion battery
BMSbattery management system
PFparticle filter
UTunscented transform
LECUlocal electronic control unit
BMUbattery management unit
HPPChybrid pulse power characteristic
SOHstate of health
DSTdynamic stress test
EMIelectromagnetic interference
RMSEroot mean square error
Ahampere-hour
Cdischarge rate
MATLABMATrix LABoratory

References

  1. Chong, U.; Yim, S.H.L.; Barrett, S.R.H.; Boies, A.M. Air quality and climate impacts of alternative bus technologies in greater London. Environ. Sci. Technol. 2014, 48, 4613–4622. [Google Scholar] [CrossRef] [PubMed]
  2. Remmlinger, J.; Buchholz, M.; Meiler, M.; Bernreuter, P.; Dietmayer, K. State-of-health monitoring of lithium-ion batteries in electric vehicles by on-board internal resistance estimation. J. Power Source 2011, 196, 5357–5363. [Google Scholar] [CrossRef]
  3. Gerssen-Gondelach, S.J.; Faaij, A.P.C. Performance of batteries for electric vehicles on short and longer term. J. Power Source 2012, 212, 111–129. [Google Scholar] [CrossRef]
  4. Lu, L.; Han, X.; Li, J.; Hua, J.; Ouyang, M. A review on the key issues for lithium-ion battery management in electric vehicles. J. Power Source 2013, 226, 272–288. [Google Scholar] [CrossRef]
  5. Xiong, R.; Sun, F.; Gong, X.; He, H. Adaptive state of charge estimator for lithium-ion cells series battery pack in electric vehicles. J. Power Source 2013, 242, 699–713. [Google Scholar] [CrossRef]
  6. Adany, R.; Aurbach, D.; Kraus, S. Switching algorithms for extending battery life in Electric Vehicles. J. Power Source 2013, 231, 50–59. [Google Scholar] [CrossRef]
  7. Li, J.; Klee Barillas, J.; Guenther, C.; Danzer, M.A. Multicell state estimation using variation based sequential Monte Carlo filter for automotive battery packs. J. Power Source 2015, 277, 95–103. [Google Scholar] [CrossRef]
  8. Sun, F.; Xiong, R.; He, H. A systematic state-of-charge estimation framework for multi-cell battery pack in electric vehicles using bias correction technique. Appl. Energy 2016, 162, 1399–1409. [Google Scholar] [CrossRef]
  9. Zhong, L.; Zhang, C.; He, Y.; Chen, Z. A method for the estimation of the battery pack state of charge based on in-pack cells uniformity analysis. Appl. Energy 2014, 113, 558–564. [Google Scholar] [CrossRef]
  10. Xiong, R.; Sun, F.; He, H.; Nguyen, T.D. A data-driven adaptive state of charge and power capability joint estimator of lithium-ion polymer battery used in electric vehicles. Energy 2013, 63, 295–308. [Google Scholar] [CrossRef]
  11. Fathabadi, H. A novel design including cooling media for Lithium-ion batteries pack used in hybrid and electric vehicles. J. Power Source 2014, 245, 495–500. [Google Scholar] [CrossRef]
  12. Hooper, J.M.; Marco, J. Characterising the in-vehicle vibration inputs to the high voltage battery of an electric vehicle. J. Power Source 2014, 245, 510–519. [Google Scholar] [CrossRef]
  13. Waag, W.; Fleischer, C.; Sauer, D.U. Critical review of the methods for monitoring of lithium-ion batteries in electric and hybrid vehicles. J. Power Source 2014, 258, 321–339. [Google Scholar] [CrossRef]
  14. Panchal, S.; Dincer, I.; Agelin-Chaab, M.; Fraser, R.; Fowler, M. Experimental and simulated temperature variations in a LiFePO4 -20 Ah battery during discharge process. Appl. Energy 2016, 180, 504–515. [Google Scholar] [CrossRef]
  15. Panchal, S.; McGrory, J.; Kong, J.; Dincer, I.; Agelin-Chaab, M.; Fraser, R.; Fowler, M. Cycling degradation testing and analysis of a LiFePO4 battery at actual conditions. Willey Int. J. Energy Res. 2017, 41, 2565–2575. [Google Scholar] [CrossRef]
  16. Unterrieder, C.; Zhang, C.; Lunglmayr, M.; Priewasser, R.; Marsili, S.; Huemer, M. Battery state-of-charge estimation using approximate least squares. J. Power Source 2015, 278, 274–286. [Google Scholar] [CrossRef]
  17. Hu, C.; Youn, B.D.; Chung, J. A multiscale framework with extended Kalman filter for lithium-ion battery SOC and capacity estimation. Appl. Energy 2012, 92, 694–704. [Google Scholar] [CrossRef]
  18. Hung, M.-H.; Lin, C.-H.; Lee, L.-C.; Wang, C.-M. State-of-charge and state-of-health estimation for lithium-ion batteries based on dynamic impedance technique. J. Power Source 2014, 268, 861–873. [Google Scholar] [CrossRef]
  19. Xiong, R.; Gong, X.; Mi, C.C.; Sun, F. A robust state-of-charge estimator for multiple types of lithium-ion batteries using adaptive extended Kalman filter. J. Power Source 2013, 243, 805–816. [Google Scholar] [CrossRef]
  20. Ng, K.S.; Moo, C.-S.; Chen, Y.-P.; Hsieh, Y.-C. Enhanced coulomb counting method for estimating state-of-charge and state-of-health of lithium-ion batteries. Appl. Energy 2009, 86, 1506–1511. [Google Scholar] [CrossRef]
  21. Xing, Y.; He, W.; Pecht, M.; Tsui, K.L. State of charge estimation of lithium-ion batteries using the open-circuit voltage at various ambient temperatures. Appl. Energy 2014, 113, 106–115. [Google Scholar] [CrossRef]
  22. Kim, J.; Lee, S.; Cho, B.H. Discrimination of Li-ion batteries based on Hamming network using discharging–charging voltage pattern recognition for improved state-of-charge estimation. J. Power Source 2011, 196, 2227–2240. [Google Scholar] [CrossRef]
  23. Wang, Y.; Zhang, C.; Chen, Z. A method for joint estimation of state-of-charge and available energy of LiFePO4 batteries. Appl. Energy 2014, 135, 81–87. [Google Scholar] [CrossRef]
  24. He, Y.; Liu, X.; Zhang, C.; Chen, Z. A new model for State-of-Charge (SOC) estimation for high-power Li-ion batteries. Appl. Energy 2013, 101, 808–814. [Google Scholar] [CrossRef]
  25. Liu, X.; Chen, Z.; Zhang, C.; Wu, J. A novel temperature-compensated model for power Li-ion batteries with dual-particle-filter state of charge estimation. Appl. Energy 2014, 123, 263–272. [Google Scholar] [CrossRef]
  26. Arulampalam, M.S.; Maskell, S.; Gordon, N.; Clapp, T. A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking. IEEE Trans. Signal Process. 2001, 50, 174–188. [Google Scholar] [CrossRef]
  27. Djuric, P.M.; Kotecha, J.H.; Zhang, J.; Huang, Y.; Ghirmai, T.; Bugallo, M.F.; Miguez, J. Particle filtering. IEEE Signal Process. Mag. 2003, 20, 19–38. [Google Scholar] [CrossRef]
  28. Kotecha, J.H.; Djuric, P.M. Gaussian particle filtering. IEEE Trans. Signal Process. 2003, 51, 2592–2601. [Google Scholar] [CrossRef]
  29. Kotecha, J.H.; Djuric, P.M. Gaussian sum particle filtering. IEEE Trans. Signal Process. 2003, 51, 2602–2612. [Google Scholar] [CrossRef]
  30. Schwunk, S.; Armbruster, N.; Straub, S.; Kehl, J.; Vetter, M. Particle filter for state of charge and state of health estimation for lithium–iron phosphate batteries. J. Power Source 2013, 239, 705–710. [Google Scholar] [CrossRef]
  31. Plett, G.L. Extended Kalman filtering for battery management systems of LiPB-based HEV battery packs: Part 1. J. Power Source 2004, 134, 252–261. [Google Scholar] [CrossRef]
  32. Plett, G.L. Extended Kalman filtering for battery management systems of LiPB-based HEV battery packs: Part 2. J. Power Source 2004, 134, 62–276. [Google Scholar] [CrossRef]
  33. Plett, G.L. Extended Kalman filtering for battery management systems of LiPB-based HEV battery packs: Part 3. J. Power Source 2004, 134, 277–292. [Google Scholar] [CrossRef]
  34. Wang, Y.; Zhang, C.; Chen, Z. A method for state-of-charge estimation of Li-ion batteries based on multi-model switching strategy. Appl. Energy 2015, 137, 427–434. [Google Scholar] [CrossRef]
  35. Xiong, R.; Sun, F.; Chen, Z.; He, H. A data-driven multi-scale extended Kalman filtering based parameter and state estimation approach of lithium-ion olymer battery in electric vehicles. Appl. Energy 2014, 113, 463–476. [Google Scholar] [CrossRef]
  36. Xiong, R.; Sun, F.C.; He, H.W. Data-driven state-of-charge estimator for electric vehicles battery using robust extended Kalman filter. Int. J. Automot. Technol. 2014, 15, 89–96. [Google Scholar] [CrossRef]
  37. Tian, Y.; Xia, B.; Sun, W.; Xu, Z.; Zheng, W. A modified model based state of charge estimation of power lithium-ion batteries using unscented Kalman filter. J. Power Source 2014, 270, 619–626. [Google Scholar] [CrossRef]
  38. Alkaya, A. Unscented Kalman filter performance for closed-loop nonlinear state estimation: A simulation case study. Electr. Eng. 2014, 96, 299–308. [Google Scholar] [CrossRef]
  39. Arasaratnam, I.; Haykin, S. Cubature Kalman filters. IEEE Trans Autom. Control 2009, 54, 1254–1269. [Google Scholar] [CrossRef]
  40. Crisan, D.; Doucet, A. A survey of convergence results on particle filtering methods for practitioners. IEEE Trans. Signal Process. 2002, 50, 736–746. [Google Scholar] [CrossRef]
  41. Chen, X.; Shen, W.; Dai, M.; Cao, Z.; Jin, J.; Kapoor, A. Robust adaptive sliding-mode observer using RBF neural network for Lithium-ion battery State of Charge estimation in electric vehicles. IEEE Trans. Veh. Technol. 2016, 65, 1936–1947. [Google Scholar] [CrossRef]
  42. Lin, C.; Mu, H.; Xiong, R.; Shen, W. A novel multi-model probability battery state of charge estimation approach for electric vehicles using H-infinity algorithm. Appl. Energy 2016, 166, 76–83. [Google Scholar] [CrossRef]
  43. Duong, V.-H.; Bastawrous, H.A.; Lim, K.; See, K.W.; Zhang, P.; Dou, S.X. Online state of charge and model parameters estimation of the LiFePO4 battery in electric vehicles using multiple adaptive forgetting factors recursive least-squares. J. Power Source 2015, 296, 215–224. [Google Scholar] [CrossRef]
  44. He, H.; Xiong, R.; Guo, H. Online estimation of model parameters and state-of-charge of LiFePO4 batteries in electric vehicles. Appl. Energy 2012, 89, 413–420. [Google Scholar] [CrossRef]
  45. Xiong, R.; He, H.; Sun, F.; Zhao, K. Evaluation on State of Charge Estimation of batteries with adaptive extended Kalman filter by experiment approach. IEEE Trans. Veh. Technol. 2013, 62, 108–117. [Google Scholar] [CrossRef]
  46. Chaoui, H.; Golbon, N.; Hmouz, I.; Souissi, R.; Tahar, S. Lyapunov-based adaptive state of charge and state of health estimation for lithium-ion batteries. IEEE Trans. Ind. Electron. 2015, 62, 1610–1618. [Google Scholar] [CrossRef]
  47. Xiong, R.; Sun, F.; Gong, X.; Gao, C. A data-driven based adaptive state of charge estimator of lithium-ion polymer battery used in electric vehicles. Appl. Energy 2014, 113, 1421–1433. [Google Scholar] [CrossRef]
  48. Sun, F.; Xiong, R. A novel dual-scale cell state-of-charge estimation approach for series-connected battery pack used in electric vehicles. J. Power Source 2015, 274, 582–594. [Google Scholar] [CrossRef]
  49. Wei, Z.; Meng, S.; Xiong, B.; Ji, D.; Tseng, K.J. Enhanced online model identification and state of charge estimation for lithium-ion battery with a FBCRLS based observer. Appl. Energy 2016, 181, 332–341. [Google Scholar] [CrossRef]
  50. Dai, H.; Xu, T.; Zhu, L.; Wei, X.; Sun, Z. Adaptive model parameter identification for large capacity Li-ion batteries on separated time scales. Appl. Energy 2016, 184, 119–131. [Google Scholar] [CrossRef]
  51. Zou, C.; Manzie, C.; Nešić, D. A framework for simplification of PDE-based lithium-ion battery models. IEEE Trans. Control Syst. Technol. 2016, 5, 1594–1609. [Google Scholar] [CrossRef]
  52. Doyle, M.; Newman, J.; Gozdz, A.S.; Schmutz, C.N.; Tarascon, J. Comparison of modeling predictions with experimental data from plastic lithium ion cells. J. Electrochem. Soc. 1996, 143, 1890–1903. [Google Scholar] [CrossRef]
  53. Kumaresan, K.; Sikha, G.; White, R.E. Thermal model for a Li-ion cell. J. Electrochem. Soc. 2008, 155, A164–A171. [Google Scholar] [CrossRef]
  54. He, H.; Xiong, R.; Guo, H.; Li, S. Comparison study on the battery models used for the energy management of batteries in electric vehicles. Energy Convers. Manag. 2012, 64, 113–121. [Google Scholar] [CrossRef]
  55. Ljung, L. Convergence Analysis of Parametric Identification Methods. IEEE Trans Autom. Control 1978, 23, 770–783. [Google Scholar] [CrossRef]
  56. Li, D.; Ouyang, J.; Li, H.Q.; Wan, J.F. State of charge estimation for LiMn2O4 power battery based on strong tracking sigma point Kalman filter. J. Power Source 2015, 279, 439–449. [Google Scholar] [CrossRef]
Figure 1. Model for a lithium-ion battery (LIB).
Figure 1. Model for a lithium-ion battery (LIB).
Energies 11 00209 g001
Figure 2. Flowchart of the method using square root cubature Kalman filter (SRCKF).
Figure 2. Flowchart of the method using square root cubature Kalman filter (SRCKF).
Energies 11 00209 g002
Figure 3. Configuration of the battery cell test bench.
Figure 3. Configuration of the battery cell test bench.
Energies 11 00209 g003
Figure 4. Hybrid pulse power characterization (HPPC) test cycle of the battery (at 1 C).
Figure 4. Hybrid pulse power characterization (HPPC) test cycle of the battery (at 1 C).
Energies 11 00209 g004
Figure 5. Relationship between state of change (SOC) and U o c v , R d c h g , C p , R p (temperature = 25 °C, state of health (SOH) = 1).
Figure 5. Relationship between state of change (SOC) and U o c v , R d c h g , C p , R p (temperature = 25 °C, state of health (SOH) = 1).
Energies 11 00209 g005
Figure 6. Relationship between SOC and U o c v , R d c h g , C p , R p (temperature = 15 °C, SOH = 0.95).
Figure 6. Relationship between SOC and U o c v , R d c h g , C p , R p (temperature = 15 °C, SOH = 0.95).
Energies 11 00209 g006
Figure 7. Test with constant discharge current using battery cell.
Figure 7. Test with constant discharge current using battery cell.
Energies 11 00209 g007
Figure 8. Results of the convergence speed in experiment A: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Figure 8. Results of the convergence speed in experiment A: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Energies 11 00209 g008
Figure 9. Results of SOC errors in experiment A: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Figure 9. Results of SOC errors in experiment A: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Energies 11 00209 g009
Figure 10. Results of the probability of SOC errors in experiment A: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Figure 10. Results of the probability of SOC errors in experiment A: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Energies 11 00209 g010
Figure 11. Dynamic stress test (DST) of the battery: (a) One test cycle; (b) 20 test cycles.
Figure 11. Dynamic stress test (DST) of the battery: (a) One test cycle; (b) 20 test cycles.
Energies 11 00209 g011
Figure 12. Results of the convergence speed in Experiment B: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Figure 12. Results of the convergence speed in Experiment B: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Energies 11 00209 g012
Figure 13. Results of SOC errors in Experiment B: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Figure 13. Results of SOC errors in Experiment B: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Energies 11 00209 g013
Figure 14. Results of the probability of SOC errors in Experiment B: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Figure 14. Results of the probability of SOC errors in Experiment B: (a) EKF, (b) UKF, (c) CKF, (d) SRCKF.
Energies 11 00209 g014
Figure 15. Results of the current and voltage with noise in Experiment C: (a) current, (b) voltage.
Figure 15. Results of the current and voltage with noise in Experiment C: (a) current, (b) voltage.
Energies 11 00209 g015
Figure 16. Results of SOC estimation in Experiment C.
Figure 16. Results of SOC estimation in Experiment C.
Energies 11 00209 g016
Figure 17. Results of the probability of SOC errors with the proposed method in Experiment C: (a) no noise, (b) 1% noise, (c) 2.5% noise, (d) 5% noise.
Figure 17. Results of the probability of SOC errors with the proposed method in Experiment C: (a) no noise, (b) 1% noise, (c) 2.5% noise, (d) 5% noise.
Energies 11 00209 g017
Table 1. Specifications of the battery cell.
Table 1. Specifications of the battery cell.
ItemValue
Cell Voltage3.7 V
Cell Capacity35 Ah
Maximum charging rates2 C
Maximum discharging rates4 C
Table 2. Specifications of the NBT BTS5200C4.
Table 2. Specifications of the NBT BTS5200C4.
ItemValue
Max value of discharging current200 A
Max value of charging current200 A
Range of voltage measurement0~5 V
Current measurement error0.1%F.S. but not better than ±20 mA
Voltage measurement error0.1%F.S. but not better than ±20 mV
Temperature measurement error±1 °C
Table 3. Equation coefficients at temperatures = 25 °C and SOH = 1.
Table 3. Equation coefficients at temperatures = 25 °C and SOH = 1.
Coefficientsgfedcba
Rp0.004107−0.0528420.316150−0.9176901.395537−1.0624470.318105
Rdchg0.001812−0.0066020.041589−0.1440260.258410−0.2266390.076809
Uocv2.9711247.396502−37.015633101.27299−143.46866100.40782−27.381834
Cp383986,475−502,0832,080,696−4,107,2963,696,315−1,244,161
Table 4. Equation coefficients at temperatures = 15 °C and SOH = 0.95.
Table 4. Equation coefficients at temperatures = 15 °C and SOH = 0.95.
Coefficientsgfedcba
Rp0.008551−0.1089540.621934−1.7580402.593975−1.9107630.554538
Rdchg0.005575−0.0272670.171668−0.5384370.840745−0.6366560.187370
Uocv2.71202310.728256−55.006330149.858773−213.2116151.254438−42.189757
Cp2677−51,640397,937−773,697312,761400,225−282,581
Table 5. Simulation results of the calculation cost.
Table 5. Simulation results of the calculation cost.
MethodsEKFUKFCKFSRCKF
Computation cost (s)0.00230.017880.0067750.00774
Table 6. Experimental results of Experiment A.
Table 6. Experimental results of Experiment A.
MethodsConvergence Time (s)RMSEMax ErrorMean Error
EKF250.012110.027270.00270
UKF1100.009700.01701−0.00268
CKF600.006020.012330.00008
SRCKF250.005990.013080.00015
Table 7. Experimental results of Experiment B.
Table 7. Experimental results of Experiment B.
MethodsConvergence Time (s)RMSEMax ErrorMean Error
EKF650.015600.05064−0.01443
UKF850.008640.03058−0.01047
CKF600.009120.02665−0.00705
SRCKF250.008920.02469−0.00704
Table 8. Experimental results of Experiment C with the proposed method.
Table 8. Experimental results of Experiment C with the proposed method.
SRCKFRMSEMax ErrorMean Error
1% noise0.010810.03489−0.00273
2.5% noise0.016930.053470.00055
5% noise0.020040.079710.00279

Share and Cite

MDPI and ACS Style

Cui, X.; Jing, Z.; Luo, M.; Guo, Y.; Qiao, H. A New Method for State of Charge Estimation of Lithium-Ion Batteries Using Square Root Cubature Kalman Filter. Energies 2018, 11, 209. https://0-doi-org.brum.beds.ac.uk/10.3390/en11010209

AMA Style

Cui X, Jing Z, Luo M, Guo Y, Qiao H. A New Method for State of Charge Estimation of Lithium-Ion Batteries Using Square Root Cubature Kalman Filter. Energies. 2018; 11(1):209. https://0-doi-org.brum.beds.ac.uk/10.3390/en11010209

Chicago/Turabian Style

Cui, Xiangyu, Zhu Jing, Maji Luo, Yazhou Guo, and Huimin Qiao. 2018. "A New Method for State of Charge Estimation of Lithium-Ion Batteries Using Square Root Cubature Kalman Filter" Energies 11, no. 1: 209. https://0-doi-org.brum.beds.ac.uk/10.3390/en11010209

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