Next Article in Journal
Design and Analysis of Electrical Distribution Networks and Balancing Markets in the UK: A New Framework with Applications
Next Article in Special Issue
Modeling and Simulation of the Thermal Runaway Behavior of Cylindrical Li-Ion Cells—Computing of Critical Parameters
Previous Article in Journal
One-Dimensional Modeling of an Entrained Coal Gasification Process Using Kinetic Parameters
Previous Article in Special Issue
A Fuzzy-Logic Power Management Strategy Based on Markov Random Prediction for Hybrid Energy Storage Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Estimation of the Electric Vehicle Power Battery State of Charge Based on the Least Squares Method and the Kalman Filter Algorithm

1
New Energy Research Center of Electric Power College, South China University of Technology, Guangzhou 510640, China
2
Guangdong Key Laboratory of Clean Energy Technology, South China University of Technology, Guangzhou 510640, China
3
College of Electrical Engineering and Automation, Henan Polytechnic University, Jiaozuo 454000, China
*
Author to whom correspondence should be addressed.
Energies 2016, 9(2), 100; https://doi.org/10.3390/en9020100
Submission received: 6 October 2015 / Revised: 27 December 2015 / Accepted: 22 January 2016 / Published: 8 February 2016
(This article belongs to the Special Issue Electrochemical Energy Storage - 2015)

Abstract

:
An estimation of the power battery state of charge (SOC) is related to the energy management, the battery cycle life and the use cost of electric vehicles. When a lithium-ion power battery is used in an electric vehicle, the SOC displays a very strong time-dependent nonlinearity under the influence of random factors, such as the working conditions and the environment. Hence, research on estimating the SOC of a power battery for an electric vehicle is of great theoretical significance and application value. In this paper, according to the dynamic response of the power battery terminal voltage during a discharging process, the second-order RC circuit is first used as the equivalent model of the power battery. Subsequently, on the basis of this model, the least squares method (LS) with a forgetting factor and the adaptive unscented Kalman filter (AUKF) algorithm are used jointly in the estimation of the power battery SOC. Simulation experiments show that the joint estimation algorithm proposed in this paper has higher precision and convergence of the initial value error than a single AUKF algorithm.

1. Introduction

In an electric vehicle, the power battery State of Charge (SOC), an important parameter of the battery state, is used to directly reflect the remaining capacity of the battery and provide a basis for the formulation of an optimal energy management strategy for the vehicle control system. An inaccurate SOC will result in a reduced performance of the vehicle and lead to potential damage to the battery system; therefore, it is critical to develop algorithms that can accurately estimate the battery SOC in real time.
An accurate estimation of the SOC is important to prolong the battery life and improve the performance of the electric vehicle [1,2]. However, because the battery is a strongly nonlinear and time-variable system, in practical applications it is hard to measure the SOC directly due to its complicated electrochemical processes and the influence of various factors [3]. At present, the most common methods of estimation [4,5,6,7,8] can be roughly divided into two main categories.
One main category is based on the relationship between energy conservation and the physical properties of the battery. For example, the most commonly used methods in this category include the open circuit voltage method and the ampere-hour integral method, among others, in which the battery charge and discharge current or the open circuit voltage are used to calculate the residual capacity of the battery.
The open circuit voltage method, when used alone, can only be applied to an electric vehicle in a non-moving state. It cannot provide a real time dynamic estimation and is therefore usually used to provide a rough SOC initial value for other methods.
The ampere-hour integral method calculates the accumulated charge of the battery during charging or discharging and among other advantages, it is economical and easy to conduct. However, when it is applied in electric vehicles, the following main problems result: (1) the SOC initial value must be obtained by other methods; (2) a higher current measurement accuracy is needed, because the accuracy of the SOC estimation is largely determined by the current measurement accuracy; and (3) the accumulative errors cannot be eliminated, and as the charging or discharging time increase, the accumulative errors may get out of control.
Another main category of methods for SOC estimation is by first establishing a mathematical model of the battery, and then the battery SOC can be estimated indirectly based on the established model, the measured charge or the discharge current, and the terminal voltage. Common methods in this category include the neural network and the Kalman filter (KF) methods, among others.
The neural network method utilizes a complex nonlinear system, i.e., a neural network, which is composed of a large number of simple neurons with extensive connections. The neural network can automatically induce, organize and study the collected data to obtain the inner rules of these data. The neural network also has the ability to map a nonlinear system and thus can better reflect the dynamic characteristics of a battery. The disadvantage of the neural network method is that a large amount of data is needed for training, and the SOC estimation accuracy is greatly influenced by the training methods and the training data.
The main idea behind the Kalman filter method is to make an optimal estimation of the minimum mean square sense of the dynamic system states. This method has strong error correction ability and requires a highly accurate battery model. When the KF method is used to estimate the SOC, the general mathematical form of the battery model can be expressed as:
State equation : x k + 1 = A x x k + B k u k + w k
Observation equation : y k = C k x k + D k u k + v k
where uk is the input of the system, including the battery current, residual capacity, battery temperature, among other variables, and yk is the output of the system, which usually indicates the terminal voltage of the battery. The most difficult challenge of the KF method is determining the state equation and the observation equation.
In this paper, using the second–order RC circuit as the equivalent model of the power battery, the online parameters of the circuit model are identified and the SOC is estimated based on the least squarea (LS) method with a forgetting factor and adaptive unscented Kalman filtering (AUKF). The comparison of these two algorithms is given, and a novel joint estimation algorithm of the power battery SOC based on the LS and the KF is proposed. The joint algorithm has the characteristics of a high estimation precision and good convergence to the initial value error. Furthermore, the advantages of the proposed algorithm are demonstrated by simulation experiments.
The structure of this paper is arranged as follows: in Section 1, the most commonly used methods for SOC estimation are introduced, and the proposed method of this paper is briefly described. In Section 2, the proposed equivalent circuit structure is determined. Dynamic parameter identification and model verification of battery model are described in Section 3. In Section 4, the AUKF algorithm is presented. In Section 5, using the LS with a forgetting factor and the AUKF algorithm to jointly estimate the power battery SOC, the advantages of the proposed algorithm is demonstrated by simulation experiments. Finally, in Section 6, the research results of this paper are summarized, and future research directions are provided.

2. Building the Battery Model

Generally, a good battery model provides an accurate description of the dynamic and static characteristics of the battery, has a relatively simple model structure making analysis and calculations easy and is not difficult to implement for a project. Currently, four main equivalent circuit models, i.e., the Rint model, Thevenin model, PNGV model, and a multi-order RC circuit model, are widely used in electric vehicle simulation [9,10,11,12,13].
The first three models have simple structures but low precision performances. In the multi-order RC circuit model, as the model order increases, the model precision will increase. However, with increasing model order, the model will not be pragmatic due to its computational complexity.
Thus, this article uses the second-order RC equivalent circuit as the battery model, as shown in Figure 1. Figure 2 shows the terminal voltage response of a 2.6 Ah Sanyo ternary lithium battery after completing a discharge cycle.
Figure 1. The second-order RC equivalent circuit.
Figure 1. The second-order RC equivalent circuit.
Energies 09 00100 g001
Figure 2. The terminal voltage response after the completing a discharge cycle.
Figure 2. The terminal voltage response after the completing a discharge cycle.
Energies 09 00100 g002
This model includes three parts:
  • Voltage source: Using the open circuit voltage Voc as the power battery electromotive force, this work neglects the temperature and SOH influences on the open circuit voltage (OCV) and only studies the relationship between Voc and the battery SOC at a constant temperature (25 °C) and constant SOH (i.e., a new battery).
  • Ohmic resistance: The battery’s resistance consists of electrode materials, electrolytes and other resistors. The change in voltage in region ① in Figure 2 is due to the influence of the ohmic resistance R.
  • RC loop circuit: Two links of a resistor and a capacitor superpose to simulate battery polarization [14], which is used to simulate the process voltage stabilization after discharge. The region ② of Figure 2 shows the change in voltage influenced by the RC loop circuit.
Equation (3) shows the function relation of the equivalent circuit model in Figure 1:
{ E ( t ) = i R + u s + u p + U ( t ) = F ( S O C ( t ) ) i = u s R s + C s d u s d t i = u p R p + C p d u p d t
We can then discretize Equation (3) and solve the state equation as follows:
[ u s , k u p , k ] = [ a s 0 0 a p ] [ u s , k 1 u p , k 1 ] + [ b s b p ] I k 1 + [ w 3 ( k ) w 5 ( k ) ]
U k = E k I k R U s , k U p , k + v ( k ) = F ( S O C k ) I k R U s , k U p , k + v ( k )
where:
{ a s = e T R s C s ,   b s = R s R s e T R s C s   a p = e T R p C p ,   b p = R p R p e T R p C p

3. Identification and Verification of Dynamic Parameters of the Battery Model

Figure 3 shows the flow chart for identifying the dynamic parameters and verifying the model. According to Figure 3, based on the LS with a forgetting factor, the dynamic parameters of an actual battery is identified and the established model is verified, in combination with the corresponding relation of the battery OCV and the SOC.
Figure 3. Parameter identification flow chart.
Figure 3. Parameter identification flow chart.
Energies 09 00100 g003

3.1. OCV-SOC Calibration Experiment

In this paper, the discharge experiment, conducted at a constant temperature (25 °C) under intermittent discharge conditions with constant current and capacity, calibrates the OCV-SOC curve with 0.2 C, 0.3 C, 0.4 C, 0.5 C, 0.6 C, 0.75 C, and 1 C. Figure 4 shows the calibration steps of the battery.
Figure 4. Calibration steps of the OCV-SOC.
Figure 4. Calibration steps of the OCV-SOC.
Energies 09 00100 g004
The corresponding relationships between OCV and SOC were recorded for x = 2, 3, 4, 5, 6, 7.5, 10, respectively. Figure 5 shows the different OCV-SOC relationships with sixth order polynomial fittings.
Figure 5. Different OCV-SOC relationships with sixth order polynomial fittings.
Figure 5. Different OCV-SOC relationships with sixth order polynomial fittings.
Energies 09 00100 g005
From Figure 5, when the SOC is above 10%, all the relationships seem to be superposed. This indicates that at the same temperature and the same SOH, any of the curves can be chosen to represent the OCV-SOC relationship. However, a smaller current leads to a smaller change of the battery characteristics. The 0.2 C constant current intermittent discharging OCV-SOC relationship is selected as the reference curve, and the open circuit voltage of the battery as a function of SOC can be represented by Equation (7):
Voc = b1 × SOC6 + b2 × SOC5 + b3 × SOC4 + b4 × SOC3 + b5 × SOC2 + b6 × SOC + b7
where a1 to a7 are coefficients obtained by the sixth order polynomial fitting giving b1 = −34.72, b2 = 120.7, b3 = −165.9, b4 = 114.5, b5 = −40.9, b6 = 7.31, and b7 = 3.231.

3.2. Application of LS with a Forgetting Factor

From Formulas (4) and (5), the Laplace equation for the battery model can be deduced:
E ( s ) U ( s ) = I ( s ) ( R + R s 1 + R s C s s + R p 1 + R p C p s )
Therefore:
G ( s ) = ( R + R s 1 + R s C s s + R p 1 + R p C p s ) = R τ s τ p s 2 + ( R τ s + R τ p + R p τ s + R s τ p ) s + R + R p + R s τ s τ p s 2 + ( τ s + τ p ) s + 1 = R s 2 + 1 τ p τ s ( R τ s + R τ p + R s τ p + R p τ s ) s + R + R s + R p τ p τ s s 2 + ( τ p + τ s ) τ p τ s s + 1 τ p τ s
where τ s is the time constant of Rs, Cs, and τ p is the time constant of Rp, Cp.
Using bilinear transform to discretize Equation (9), s = 2 T 1 z 1 1 + z 1 is obtained, and the discrete transfer function is as follows:
G ( z 1 ) = a 3 + a 4 z 1 + a 5 z 2 1 a 1 z 1 a 2 z 2
where a1, a2, a3, a4, a5 are the corresponding constant coefficients. Formula (9) can then be converted to a differential equation:
y ( k ) = E ( k ) U ( k ) = a 1 y ( k 1 ) + a 2 y ( k 2 ) + a 3 I ( k ) + a 4 I ( k 1 ) + a 5 I ( k 2 )
where I(k) and y(k) indicate the system input and output, respectively, and subsequently gives:
φ ( k ) = [ y ( k 1 )   y ( k 2 )   I ( k )   I ( k 1 )   I ( k 2 ) ] T ,   θ = [ a 1   a 2     a 3     a 4     a 5 ] ;
If we assume the k moment sensor sampling error is e(k), then:
y ( k ) = φ T ( k ) θ + e ( k )
and expanding ϕ   ( k ) to N-dimensional, where k = 1, 2, 3…N+ n, n = 2, the following equation is deduced:
{ Φ ( k ) = [ y ( 2 ) y ( 1 ) I ( 3 ) I ( 2 ) I ( 1 ) y ( 3 ) y ( 2 ) I ( 4 ) I ( 3 ) I ( 2 ) y ( N + 1 ) y ( N ) I ( N + 2 ) I ( N + 1 ) I ( N ) ] Y = [ y ( 3 ) y ( 4 ) y ( 5 ) y ( N + 2 ) ] T e = [ e ( 3 ) e ( 4 ) e ( 5 ) e ( N + 2 ) ] T k 3
Next, taking the criterion Equation J ( θ ) :
J ( θ ) = i = 1 N ( Y Φ θ ) 2 = i = 1 N ( e ( i + 2 ) ) 2
and considering the least squares method is to take the J ( θ ) minimum:
J ( θ ) ( θ ) = ( θ ) [ ( Y Φ θ ) T ( Y Φ θ ) ] = 0
which gives:
θ ^ = [ Φ T Φ ] 1 Φ T Y
The above equations constitute a one-time least square calculation. However, for the actual system, a one-time calculation does not give an estimated value close to the true value. Thus, a recursive least square method is introduced:
{ θ ^ ( k + 1 ) = θ ^ ( k ) + K ( k + 1 ) [ y ( k + 1 ) Φ T ( k + 1 ) θ ^ ( k ) ] K ( k + 1 ) = P ( k + 1 ) Φ ( k + 1 ) P ( k + 1 ) = P ( k ) P ( k ) Φ ( k + 1 ) Φ T ( k + 1 ) P ( k ) I + Φ T ( k + 1 ) P ( k ) Φ ( k + 1 )
where θ ^ ( k ) is the system-estimated reference value of the previous cycle, Φ T ( k + 1 ) θ ^ ( k ) is the observed value of this cycle, y ( k + 1 ) is the actual observed value of the system, Φ T ( k + 1 ) θ ^ ( k ) is the prediction error, and K ( k + 1 ) is the corrected prediction value. To obtain the optimal estimation of the present cycle, θ ( 0 ) and P ( 0 ) must be first provided to meet the requirements, K ( k + 1 ) is then obtained, and the least square method is executed. Generally, θ ^ ( 0 ) can be any value, and P ( 0 ) = α I , where α is a positive real number, and I is a unit matrix.
The recursive least squares method has an unlimited memory, i.e., as the length K increases, the older data accumulates, making new data difficult to substitute into the least-square steps. This will subsequently affect parameter estimation, especially in time-varying systems. Because a large amount of accumulated old data creates an imbalance with the new data, the newly estimated parameters cannot accurately reflect the characteristics of the system at a current moment. Thus, to avoid the above situation [15], a forgetting factor λ , where 0 < λ < 1 , is introduced:
P 1 ( k + 1 ) = λ P 1 ( k ) + Φ ( k + 1 ) Φ T ( k + 1 )
Thus, even when ( N + 1 ) is large, P ( N + 1 ) does not go to 0, and “data saturation” can be eliminated.
The steps of the least square algorithm with the forgetting factor are as follows:
{ θ ^ ( k + 1 ) = θ ^ ( k ) + K ( k + 1 ) [ y ( k + 1 ) Φ T ( k + 1 ) θ ^ ( k ) ] K ( k + 1 ) = P ( k ) Φ ( k + 1 ) [ λ + Φ T ( k + 1 ) P ( k ) Φ ( k + 1 ) ] 1 P ( k + 1 ) = 1 λ [ I K ( k + 1 ) Φ T ( k + 1 ) ] P ( k )
In Equation (19), λ = 1 the most common least squares. When λ is smaller, the tracking ability is stronger, but the volatility is greater; hence, generally 0.95 < λ < 1.

3.3. Dynamic Parameter Identification

Parameter identification is based on the information of the measurement system and provides guidelines to estimate the model structure and unknown parameters. According to the value of θ derived from the previous algorithm:
z 1 = 1 T 2 s 1 + T 2 s
Substituting this into Equation (10):
G ( s ) = a 3 a 4 + a 5 1 + a 1 a 2 s 2 + 4 ( a 3 a 5 ) T ( 1 + a 1 a 2 ) s + 4 ( a 3 a 4 + a 5 ) T 2 ( 1 + a 1 a 2 ) s 2 + 4 ( 1 + a 2 ) T ( 1 + a 1 a 2 ) s + 4 ( 1 a 1 a 2 ) T 2 ( 1 + a 1 a 2 )
Because the corresponding coefficients of Equations (9) and (21) are equal, we can obtain:
{ R = a 3 a 4 + a 5 1 + a 1 a 2 τ s τ p = T 2 ( 1 + a 1 a 2 ) 4 ( 1 a 1 a 2 ) τ s + τ p = T ( 1 + a 2 ) 1 a 1 a 2 R + R s + R p = a 3 a 4 + a 5 1 a 1 a 2 R τ s + R τ p + R p τ s + R s τ p = 4 ( a 3 a 5 ) T ( 1 + a 1 a 2 )
The coefficients on the right-hand side of Equation (22) can be obtained by a recursive algorithm, and the variables on the left-hand side are the unknown parameters of the battery model. This completes the process of parameter identification.
In the process of identifying model parameters, the known variables are V(k), I(k), V(k−1), I(k−1), SOC(k−1) and V(k−2), I(k−2), and the unknown variable is θ = [ a 1   a 2   a 3   a 4   a 5 ] T . The steps of using the LS method with a forgetting factor for the identification of dynamic parameters of the battery are as follows:
  • Identification initialization using sampling time T = 1s, and SOC (0) = 90%.
  • Calculate V o c ( k ) V ( k ) each time and obtain the input Φ ( k ) and output y ( k ) of the identification process accurately.
  • Initialize θ ( 0 ) , P ( 0 ) and the forgetting factor λ , and start the forgetting factor least square parameter identification; in this paper, α = 5000 , λ = 0.96 .
Using this process, the value of θ can be obtained, and then according to Equation (22), R, Rs, Rp, Cs, and Cp can consequently be obtained; hence, the dynamic real-time update of the battery model parameters is realized, along with an accurate description of the dynamic response of the battery. Furthermore, the accuracy of the battery model is improved, and the basis for estimating the battery SOC accurately is provided in latter sections.

3.4. Model Verification

After the dynamic parameters of the battery model are determined, the next step is to verify the accuracy of the model using Hybrid Pulse Power Characterization (HPPC) [16]. Here, the initial SOC is set to 0.5, and the input current waveform is shown in Figure 6 with a pulse current size of 1 C (2.6 A).
Figure 6. HPPC pulse current.
Figure 6. HPPC pulse current.
Energies 09 00100 g006
A comparison of voltage responses is shown in Figure 7, and Figure 8 depicts the voltage error, i.e., the differences between the measured and estimated voltages.
Figure 7. Comparison of voltage responses.
Figure 7. Comparison of voltage responses.
Energies 09 00100 g007
Figure 8. Differences between measured and estimated voltages.
Figure 8. Differences between measured and estimated voltages.
Energies 09 00100 g008
Figure 7 and Figure 8 show that when the current suddenly changes, the model estimated voltage can track the actual voltage well, and the error remains at approximately 0.01 V. Hence, this model can be used to verify the algorithms of the SOC estimation in this paper.

4. Establishment of the AUKF Algorithm

The main idea of the Kalman filter is to make an optimal estimation of the minimum mean square value, which includes the following two stages: prediction and updating. In the prediction stage, the filter makes an estimate of the current state according to the value of the last state. In the updating stage, according to the observed value of the current state, the filter optimizes the predicted value from the prediction stage to obtain a more accurate estimation of the current state.
It is important to note that the Kalman filter is mainly used in linear systems, while a battery system reflects complex nonlinear characteristics. Some people [4,7,8] have used the extended Kalman filter (EKF) for SOC estimation, and while some good results have been achieved, linearization errors are inevitable, and the Jacobian matrix is also difficult to estimate. In recent years, a new nonlinear filtering method has emerged, collectively referred to as the sigma point Kalman filter, including the unscented Kalman filter (UKF). UKF does not require Taylor approximations of nonlinear equations; instead, the nonlinear unscented transform (UT) technique is used directly, and thus the mean and the variance of the nonlinear system states can be mapped directly to achieve a higher estimation accuracy.
In a normal UKF algorithm [17,18], the covariance is a constant and cannot satisfy the real-time dynamic characteristics of the noise, which has a certain impact on the accuracy. In this paper, to eliminate this effect, the normal UKF algorithm is improved by updating the covariance in real-time and thus improving the accuracy of the UKF. This type of algorithm is called the adaptive unscented Kalman filter (AUKF) algorithm. The establishment process of the algorithm is described as follows.
A discrete-time controlled system is governed by the equation of state and the observation equations are shown in Equation (23):
{ x k + 1 = A x x k + B k u k + w k y k = C k x k + D k u k + v k
where the random variables wk and vk represent the process and measurement noise, respectively. As for the UKF, the iteration equation is based on a certain set of sample points, which is chosen to make their mean value and variance consistent with the mean value and variance of the state variables. Then, these points will recycle the equation of the discrete-time process model to produce a set of predicted points. After that, the mean value and the variance of the predicted points will be calculated to modify the results, and the mean value and the variance will be estimated. Before the UKF recursion, the state variables must be modified in a superposition of the process noise and the measurement noise of the original states. The SOC of the Li-ion battery pack can be calculated using the ampere-hour integral method:
S O C ( t ) = SOC ( t ) 1 C N t t η i d t
In Equation (24), η = k t k i k c , where ki is the compensation coefficient of charge and discharge rate, kt is the temperature compensation coefficient, and kc is the compensation coefficient of cycles. CN is the actual available battery capacity.
The Li-ion battery equation of state can be obtained from Equations (4) and (24):
[ S O C k u s , k u p , k ] = [ 1 0 0 0 a s 0 0 0 a p ] [ S O C k 1 u s , k 1 u p , k 1 ] + [ η T / C N b s b p ] I k 1 + [ w 1 ( k ) w 3 ( k ) w 5 ( k ) ]
U k = E k I k R U s , k U p , k + v ( k ) = F ( S O C k ) I k R U s , k U p , k + v ( k )
For the circuit model shown in the Equations (25) and (26):
X k = [ x k T ,   w k T ,   v k T ] T = [ S O C k ,   u s , k ,   u p , k ,   w s o c , k ,   w s , k ,   w p , k ,   v k ] T
To facilitate the distinction, we can take x k = [ S O C k ,     U s , k ,   U p , k ] as the initial state of the system, yk as the raw output (its corresponding symbol is Uk in the circuit model of the Li-ion battery), uk as the control variable (its corresponding symbol is Ik), and make Ψ = [ y 1 , y 2 , y 3 , , y k   ] . The operations of an ordinary UKF are as follows.
First, select (2L + 1) sampling points, and make Sample= { z i , X k 1 , i } , i =0, 1, 2, , 2 L + 1 , where X k 1 , i is the selected points and z i is the corresponding weighting value. Then, select the points in the following manner:
{ X k 1 , 0 = X ^ k 1   X k 1 , i = X ^ k 1 + ( ( L + λ ) P X , k 1 ) i   i = 1 ~ L X k 1 , i = X ^ k 1 ( ( L + λ ) P X , k 1 ) i   i = ( L + 1 ) ~ 2 L  
The corresponding weighting values are:
{ z 0 ( m ) = λ L + λ z 0 ( c ) = λ L + λ + ( 1 + α 2 + β ) z i ( m ) = z i ( c ) = 1 2 ( L + λ )   i = 1 ~ 2 L
where λ = α 2 ( L + t ) L , z ( m ) is the corresponding weighting value of the mean, z ( c ) is the corresponding weighting value of the variance, and ( ( L + λ ) P X , k 1 ) i denotes the column i of the square-rooting matrix ( L + λ ) P X , k 1 . To ensure that the covariance matrix is definitely positive, we must take t 0 ; α controls the distance of the selected points, with 10 2 α 1 , and β is used to reduce the error of the higher-order terms. For a Gaussian, the optimal choice is β = 2 in this paper, along with t = 0 and α = 1 . X ^ k 1 , i consists of X ^ k 1 , i x ,   X ^ k 1 , i w and X ^ k 1 , i v .
The time updates of the iterative process are as follows:
{ X ^ k | k 1 = E { [ f ( x k 1 , u k 1 ) + ω k 1 ] | Ψ k 1 } = i = 0 2 L z i ( m ) [ A k 1 x k 1 , i x + B k 1 u k 1 + x k 1 , i w ] = i = 0 2 L z i ( m ) X k | k 1 , i x P x , k | k 1 = E [ ( x k x ^ k | k 1 ) ( x k x ^ k | k 1 ) T ] = i = 0 2 L z i ( c ) ( X k | k 1 , i x x ^ k | k 1 ) ( X k | k 1 , i x x ^ k | k 1 ) T y ^ k = E { [ h ( x k , u k ) + ν k ] | Ψ k 1 } = i = 0 2 L z i ( m ) [ h ( X k | k 1 , i x , u k ) + x k 1 , i ν ] = i = 0 2 L z i ( m ) y k | k 1 , i
The gain matrix of the filter is:
L k = P x y , k P y , k 1 = i = 0 2 L z i ( c ) ( X k | k 1 , i x x ^ k | k 1 ) ( y k | k 1 , i y ^ k ) T [ i = 0 2 L z i ( c ) ( y k | k 1 , i y ^ k ) ( y k | k 1 , i y ^ k ) T ] 1
The estimates of the states and the mean square error are as follows:
{ x k = x ^ k | k 1 + L k ( y k y ^ k ) P x , k = P x , k | k 1 L k P y , k L k T
where yk is the actual measured value of the system output. As the process noise and measurement noise are real-time and to make the covariance of the process noise and measurement noise update in real time, the following is needed:
{ μ k = y k H [ x ^ k , u k ] F k = μ k μ k T R k v = ( F k + i = 0 2 L z i c ( y k | k 1 , i y k ) ( y k | k 1 , i y k ) T ) / 2 R k w = L k F k L k T
where μ k is the residual error of the system measured output and y k | k 1 , i is the residual error of the system measured output estimated by the sigma points.

5. Experiments with the SOC Estimation Algorithm

5.1. Reasons and Steps of the Joint Estimation Algorithm

The LS method with a forgetting factor undertakes the work of parameter identification, and the AUKF functions in the progress of the SOC estimation. The characteristics of these two algorithms [17,18,19,20] are shown in Table 1.
Table 1. The characteristics of the two algorithms.
Table 1. The characteristics of the two algorithms.
AlgorithmAdvantagesDisadvantages
LS with a forgetting factorThis algorithm does not require the observation data to provide the probability and statistics of the noise under random conditions; the statistical properties are quite good.This algorithm cannot identify an unbiased, coherent parameter with colored noise.
AUKFThis algorithm has a strong immunity to the disturbance of the initial value; the iterative calculations ensure the acquirement of the desired value.In theory, the minimum variance estimation can be obtained only when the statistical properties are known.
Based on the merits and the drawbacks mentioned above, a joint algorithm of these two algorithms is proposed. The implementation of the joint algorithm can be divided briefly into two steps. First, the Kalman filter model updates parameters using the data provided by the LS method with a forgetting factor. Then, the filter generates the SOC, which will be used to deduce the OCV. Second, the OCV combines the measured voltage and the current value to update the LS estimation result for the next reiteration. Figure 9 illustrates the steps of the joint algorithm.
Figure 9. The combination of the LS method with a forgetting factor and the AUKF.
Figure 9. The combination of the LS method with a forgetting factor and the AUKF.
Energies 09 00100 g009
The details of the implemented algorithm are as follows:
  • The BMS measures the voltage of the Li-ion battery in the static state, according to the function of OCV-SOC, and the initial value of SOC (0) is calculated.
  • The initial values of the model, i.e., R(0), Rs(0), Rp(0), Cs(0), and Cp(0), are estimated, according to the current and voltage responses in the early stage of battery operation.
  • The initial values of the model are used to calculate the initial coefficients, and then the adaptive unscented Kalman Filter will be used to obtain the SOC value at the current moment.
  • The open circuit voltage Voc at that time is calculated according to the function relationship of OCV-SOC. Then, the parameters of the model at the current instance are obtained using the LS method with a forgetting factor.
  • The model parameters are utilized to update the corresponding coefficients, and the AUKF is used again to calculate the estimated value of SOC in the next instance, and step 4 is repeated.
Step 4 is applied to calculate the model parameters, and step 5 is utilized to estimate the SOC values. These two steps are repeated, and the Li-ion battery parameters and the estimated SOC values at every instance will be obtained.

5.2. Experimental Analysis of the Joint Estimation Algorithm

In this section, simulations of the joint algorithm mentioned above and a single AUKF are presented, and analysis will be made on the accuracy of the algorithms and the convergence to the initial value error.
In the process of setting the battery discharging and charging states, to correspond to the experimental object in the OCV-SOC calibration, the current signal in Figure 10 is adopted to describe the increase or decrease of the current in the discharging or charging process of the power battery. In one period, the average output current is 1.77 A, and the maximum discharging current is 5.28 A, the maximum charging current is 2.42 A. Each period is 1367 s, and the condition lasts two periods.
Figure 10. The input current waveform.
Figure 10. The input current waveform.
Energies 09 00100 g010

5.2.1. Accuracy Verification

In this simulation model, the input current I(k) is integrated using the Ampere-hour (Ah) integral method. As there is no error in the current measurement due to outside disturbances, no accumulative error exists; thus, the integration of the current in the simulation model could be regarded as the theoretical value of the SOC. Figure 11 displays a comparison of the results generated by the two algorithms and the theoretical value. Figure 12 compares the error of these two algorithms.
Figure 11 and Figure 12 illustrate that both algorithms are able to follow the theoretical values of SOC. The maximum error of the joint estimation algorithm is 1.3% and that of AUKF is 2.5%, which indicates a greater accuracy of the joint estimation algorithm.
Figure 11. Comparison of the SOC estimates.
Figure 11. Comparison of the SOC estimates.
Energies 09 00100 g011
Figure 12. Comparison of the SOC estimation error.
Figure 12. Comparison of the SOC estimation error.
Energies 09 00100 g012

5.2.2. Convergence to the SOC Initial Value Error

In the early operational stage of the Li-ion battery, errors exist between the measured and real voltage and current values, which will subsequently lead to an erroneous SOC value. Taking this condition into consideration, the convergence to the initial value error is a vital performance index of an algorithm. This section will discuss this issue. We assume that the real initial value of the SOC is 1, while a SOC value of 0.96 is deduced from the measurements. Figure 13 shows the theoretical value and results of the two algorithms, and Figure 14 is the error comparison of these two algorithms in a time period of 100 s.
Figure 13 and Figure 14 indicate that both algorithms have the ability to converge to the SOC initial value error. At 100 s, the error percentage of the joint estimation algorithm is 1.628%, while that of the AUKF is higher at 3.054%; hence, the joint estimation algorithm is better at converging to the SOC initial value error.
Figure 13. Results of the theoretical value and the two algorithms.
Figure 13. Results of the theoretical value and the two algorithms.
Energies 09 00100 g013
Figure 14. Error comparison in a time period of 100 s.
Figure 14. Error comparison in a time period of 100 s.
Energies 09 00100 g014

6. Conclusions

Power battery SOC is a vital state information of an electrical vehicle and is strongly nonlinear and time-varying. The proposed joint estimation algorithm combines the LS method with a forgetting factor and the AUKF method. This work can be summarized as follows:
  • A battery model was built according to its external characteristics, and the parameters were identified and verified at the same time.
  • The advantages and disadvantages of the LS method with a forgetting factor and the AUKF method were analyzed, and a joint estimation algorithm was proposed.
  • Comparison of the joint estimation algorithm and the AUKF, including the accuracy and the ability to converge to initial value errors, were conducted, and it was concluded that the joint estimation was better than the AUKF.
Suggested future research directions building off of the current data are as follows:
  • All of the current data are obtained in a constant temperature environment. Thus, data of the battery operating in a variable temperature environment should be obtained in the future, to evaluate the exact relationship of SOC-OCV.
  • The calculations in this paper do not consider the battery health; hence, the joint estimation algorithm should be applied on batteries with different battery health conditions.

Acknowledgments

This work is supported by DongGuan Innovative Research Team Program (No. 201460711900131) and the National Natural Science Foundation of China (Grant No. 51377058). I would like to express my deepest gratitude to my supervisor, Longyun Kang, who has provided me with valuable guidance at every stage of writing this paper. I would also like to thank the anonymous reviewers for dedicating the time to review my paper despite their busy schedules.

Author Contributions

This research article has five authors. The circuit structure was designed by Xiangwei Guo and Longyun Kang. Xiangwei Guo and Zhizhen Huang designed the research methods and control strategies. Xiangwei Guo, Yuan Yao and Wenbiao Li designed and performed the experiments. Longyun Kang provided the experimental environment. Xiangwei Guo wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Piller, S.; Perrin, M.; Jossen, A. Methods for state-of-charge determination and their applications. J. Power Sources 2001, 96, 113–120. [Google Scholar] [CrossRef]
  2. Álvarez Antón, J.C.; García Nieto, P.J.; de Cos Juez, F.J. Battery State-of-Charge Estimator Using the MARS Technique. IEEE Trans. Power Electron. 2013, 28, 3798–3805. [Google Scholar] [CrossRef]
  3. He, H.; Xiong, R.; Fan, J. Evaluation of Lithium-Ion Battery Equivalent Circuit Models for State of Charge Estimation by an Experimental Approach. Energies 2011, 4, 582–598. [Google Scholar] [CrossRef]
  4. Chen, Z.; Fu, Y.; Mi, C.C. State of Charge Estimation of Lithium-Ion Batteries in Electric Drive Vehicles Using Extended Kalman Filtering. IEEE Trans. Veh. Technol. 2013, 62, 1020–1030. [Google Scholar] [CrossRef]
  5. 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]
  6. Alvarez Anton, J.C.; Garcia Nieto, P.J.; Blanco Viejo, C.; Vilan Vilan, J.A. Support Vector Machines Used to Estimate the Battery State of Charge. IEEE Trans. Power Electron. 2013, 28, 5919–5926. [Google Scholar] [CrossRef]
  7. Sepasi, S.; Roose, L.; Matsuura, M. Extended Kalman Filter with a Fuzzy Method for Accurate Battery Pack State of Charge Estimation. Energies 2015, 8, 5217–5233. [Google Scholar] [CrossRef]
  8. Yu, Z.; Huai, R.; Xiao, L. State-of-Charge Estimation for Lithium-Ion Batteries Using a Kalman Filter Based on Local Linearization. Energies 2015, 8, 7854–7873. [Google Scholar] [CrossRef]
  9. Rodrigues, S.; Munichandraiahb, N.; Shuklaa, A.K. A review of state-of-charge indication of batteries by means of a.c. impedance measurements. J. Power Sources 2000, 87, 12–20. [Google Scholar] [CrossRef]
  10. Watrin, N.; Ostermann, H.; Blunier, B.; Miraoui, A. Multiphysical Lithium-Based Battery Model for Use in State-of-Charge Determination. IEEE Trans. Veh. Technol. 2012, 61, 3420–3429. [Google Scholar] [CrossRef]
  11. Corno, M.; Bhatt, N.; Savaresi, S.M.; Verhaegen, M. Electrochemical Model-Based State of Charge Estimation for Li-Ion Cells. IEEE Trans. Control. Syst. 2015, 23, 117–127. [Google Scholar] [CrossRef]
  12. Hu, Y.; Wang, Y. Two Time-Scaled Battery Model Identification with Application to Battery State Estimation. IEEE Trans. Control. Syst. 2015, 23, 1180–1188. [Google Scholar] [CrossRef]
  13. Tong, S.; Klein, M.P.; Park, J.W. On-line optimization of battery open circuit voltage for improved state-of-charge and state-of-health estimation. J. Power Sources 2015, 293, 416–428. [Google Scholar] [CrossRef]
  14. Hu, Y.; Yurkovich, S. Linear parameter varying battery model identification using subspace methods. J. Power Sources 2011, 196, 2913–2923. [Google Scholar] [CrossRef]
  15. Duong, V.; 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 Sources 2015, 296, 215–224. [Google Scholar] [CrossRef]
  16. Ranjbar, A.H.; Banaei, A.; Khoobroo, A.; Fahimi, B. Online Estimation of State of Charge in Li-Ion Batteries Using Impulse Response Concept. IEEE Trans. Smart Grid 2012, 3, 360–367. [Google Scholar] [CrossRef]
  17. Partovibakhsh, M.; Liu, G. An Adaptive Unscented Kalman Filtering Approach for Online Estimation of Model Parameters and State-of-Charge of Lithium-Ion Batteries for Autonomous Mobile Robots. IEEE Trans. Control Syst. 2015, 23, 357–363. [Google Scholar] [CrossRef]
  18. Aung, H.; Soon Low, K.; Ting Goh, S. State-of-Charge Estimation of Lithium-Ion Battery Using Square Root Spherical Unscented Kalman Filter (Sqrt-UKFST) in Nanosatellite. IEEE Trans. Power Electron. 2015, 30, 4774–4783. [Google Scholar] [CrossRef]
  19. Zhang, C.; Wang, L.Y.; Li, X.; Chen, W.; Yin, G.G.; Jiang, J. Robust and Adaptive Estimation of State of Charge for Lithium-Ion Batteries. IEEE Trans. Ind. Electron. 2015, 62, 4948–4957. [Google Scholar] [CrossRef]
  20. Rahimi-Eichi, H.; Baronti, F.; Chow, M. Online Adaptive Parameter Identification and State-of-Charge Coestimation for Lithium-Polymer Battery Cells. IEEE Trans. Ind. Electron. 2014, 61, 2053–2061. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Guo, X.; Kang, L.; Yao, Y.; Huang, Z.; Li, W. Joint Estimation of the Electric Vehicle Power Battery State of Charge Based on the Least Squares Method and the Kalman Filter Algorithm. Energies 2016, 9, 100. https://doi.org/10.3390/en9020100

AMA Style

Guo X, Kang L, Yao Y, Huang Z, Li W. Joint Estimation of the Electric Vehicle Power Battery State of Charge Based on the Least Squares Method and the Kalman Filter Algorithm. Energies. 2016; 9(2):100. https://doi.org/10.3390/en9020100

Chicago/Turabian Style

Guo, Xiangwei, Longyun Kang, Yuan Yao, Zhizhen Huang, and Wenbiao Li. 2016. "Joint Estimation of the Electric Vehicle Power Battery State of Charge Based on the Least Squares Method and the Kalman Filter Algorithm" Energies 9, no. 2: 100. https://doi.org/10.3390/en9020100

APA Style

Guo, X., Kang, L., Yao, Y., Huang, Z., & Li, W. (2016). Joint Estimation of the Electric Vehicle Power Battery State of Charge Based on the Least Squares Method and the Kalman Filter Algorithm. Energies, 9(2), 100. https://doi.org/10.3390/en9020100

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