A New Method for State of Charge Estimation of Lithium-Ion Battery Based on Strong Tracking Cubature Kalman Filter

The estimation of state of charge (SOC) is a crucial evaluation index in a battery management system (BMS). The value of SOC indicates the remaining capacity of a battery, which provides a good guarantee of safety and reliability of battery operation. It is difficult to get an accurate value of the SOC, being one of the inner states. In this paper, a strong tracking cubature Kalman filter (STCKF) based on the cubature Kalman filter is presented to perform accurate and reliable SOC estimation. The STCKF algorithm can adjust gain matrix online by introducing fading factor to the state estimation covariance matrix. The typical second-order resistor-capacitor model is used as the battery’s equivalent circuit model to dynamically simulate characteristics of the battery. The exponential-function fitting method accomplishes the task of relevant parameters identification. Then, the developed STCKF algorithm has been introduced in detail and verified under different operation current profiles such as Dynamic Stress Test (DST) and New European Driving Cycle (NEDC). Making a comparison with extended Kalman filter (EKF) and CKF algorithm, the experimental results show the merits of the STCKF algorithm in SOC estimation accuracy and robustness.


Introduction
In recent years, energy-saving and emission reduction attract special attentions and environmental pollution have become more and more critical.Electric vehicles (EVs), characterized by environmentally benign property and zero-emissions, are with people's horizons.Research and development on electric vehicle batteries have been greatly active, especially because of the increasing cost of energy consumption.Lithium-ion batteries conform to the demands of EVs and hybrid electric vehicles (HEVs) for their high energy density, safety, low self-discharge and long cycle life.The system that regulates the energy flow in a battery pack with respect to voltages of individual cells, temperature, state of charge and health is defined as the Battery Management System (BMS) [1,2].The main task of BMS is to maintain a safe operating environment for the battery system, and to protect it from damage.The state of charge (SOC), defined as the ratio of the remaining capacity to the rated capacity, is the core part of BMS.An accurate SOC estimation prevents over-charging or over-discharging time and prolongs the life of batteries on the one hand, and, on the other hand, it helps drivers effectively grasp running condition of their electric vehicle.Because the SOC is Energies 2015, 8,[13458][13459][13460][13461][13462][13463][13464][13465][13466][13467][13468][13469][13470][13471][13472] In 2009, a new nonlinear filter based on KF called as cubature Kalman filter (CKF) was proposed by Arasaratnam and Haykin [21].The key idea of this method is using the radial-spherical cubature rule to capture the mean and covariance of the states of a nonlinear system with additive Gaussian noise.The CKF uses a set of 2n points to predict the SOC of the battery, while the UKF uses a set of 2n + 1 points, where n represents the state-vector dimension.It has already been proven more accurate and stable than the UKF [22,23].The CKF has been successfully applied in many fields, such as moving-target tracking [24] and spacecraft attitude estimation [25].
In this paper, a new strong tracking cubature filter (STCKF) method for SOC estimation is presented to improve the robustness and accuracy of the CKF method.Two typical driving cycles, including Dynamic Stress Test (DST) and NEDC, are applied to assess the performance of the proposed method.The experimental results show that this method can obtain super performance in terms of accuracy and robustness, better than EKF and CKF.
The remainder of this paper is organized as follows.In Section 2, the second-order resistor-capacitor equivalent circuit model is used and the model parameters are determined by exponential-function fitting method.The CKF and STCKF methods are illustrated in detail in Section 3. Section 4 describes the experimental configurations.The verification results and comparisons of different approaches are discussed in Section 5, and conclusions of the paper are made in Section 6.

Battery Model
Equivalent circuit models, which consist of basic circuit elements (e.g., resistors, capacitors and inductors), are commonly used to simulate the voltage-current behaviors of a battery [26][27][28].Building an effective equivalent circuit model (ECM) makes a great difference for researching battery management system.Estimation accuracy is certainly based on the precise ECM for SOC estimation algorithms.Then, the application of ECM prevents batteries from unnecessary and permanent damage in actual operation.
When choosing ECM, the following several aspects must be considered.First of all, the model can clearly imitate the dynamic behaviors of the battery.Secondly, it must not be too complicated to establish the state-space equations so complex computations do not increase the burden on the CPU.The complicated ECM can well simulate dynamic behaviors while it increases the computation cost.Therefore, a tradeoff must be made to meet the above requirements.
Hence, this paper is based on a second-order equivalent circuit battery model.As shown in Figure 1, the second-order RC battery model consists of an open-circuit voltage U oc (SOC) that has a certain relationship with SOC at the same temperature, a resistor R o , and two parallel resistor-capacitor (RC) networks connected in series (i.e., R 1 -C 1 and R 2 -C 2 ).The resistor R o is the ohmic resistance caused by the accumulation and dissipation of charge in the electrical double-layer, R 1 and C 1 are the activation polarization resistance and capacitance, respectively, while R 2 and C 2 separately are the concentration polarization resistance and capacitance, respectively.and stable than the UKF [22,23].The CKF has been successfully applied in many fields, such as moving-target tracking [24] and spacecraft attitude estimation [25].
In this paper, a new strong tracking cubature filter (STCKF) method for SOC estimation is presented to improve the robustness and accuracy of the CKF method.Two typical driving cycles, including Dynamic Stress Test (DST) and NEDC, are applied to assess the performance of the proposed method.The experimental results show that this method can obtain super performance in terms of accuracy and robustness, better than EKF and CKF.
The remainder of this paper is organized as follows.In Section 2, the second-order resistor-capacitor equivalent circuit model is used and the model parameters are determined by exponential-function fitting method.The CKF and STCKF methods are illustrated in detail in Section 3. Section 4 describes the experimental configurations.The verification results and comparisons of different approaches are discussed in Section 5, and conclusions of the paper are made in Section 6.

Battery Model
Equivalent circuit models, which consist of basic circuit elements (e.g., resistors, capacitors and inductors), are commonly used to simulate the voltage-current behaviors of a battery [26][27][28].Building an effective equivalent circuit model (ECM) makes a great difference for researching battery management system.Estimation accuracy is certainly based on the precise ECM for SOC estimation algorithms.Then, the application of ECM prevents batteries from unnecessary and permanent damage in actual operation.
When choosing ECM, the following several aspects must be considered.First of all, the model can clearly imitate the dynamic behaviors of the battery.Secondly, it must not be too complicated to establish the state-space equations so complex computations do not increase the burden on the CPU.The complicated ECM can well simulate dynamic behaviors while it increases the computation cost.Therefore, a tradeoff must be made to meet the above requirements.
Hence, this paper is based on a second-order equivalent circuit battery model.As shown in Figure 1, the second-order RC battery model consists of an open-circuit voltage Uoc(SOC) that has a certain relationship with SOC at the same temperature, a resistor Ro, and two parallel resistor-capacitor (RC) networks connected in series (i.e., R1-C1 and R2-C2).The resistor Ro is the ohmic resistance caused by the accumulation and dissipation of charge in the electrical double-layer, R1 and C1 are the activation polarization resistance and capacitance, respectively, while R2 and C2 separately are the concentration polarization resistance and capacitance, respectively.The cell SOC is defined as the coefficient of the remaining capacity to the nominal capacity.The SOC can be expressed as the following equation: The cell SOC is defined as the coefficient of the remaining capacity to the nominal capacity.The SOC can be expressed as the following equation: where i(t) represents the value of current (defined to be negative for discharging and positive for charging).Q N is the nominal capacity.SOC(t 0 ) represents the SOC value at time t 0 .
The derivation of SOC based on Equation (1) can be formulated as: According to the circuit theory, the electrical behavior of the second-order RC battery model can be expressed as: where U 1 and U 2 denote the terminal voltage of C 1 and C 2 , respectively; and U t and i(t) denote individually the value of terminal voltage and current, respectively.U oc denotes the open-circuit voltage, which is usually a nonlinear function of SOC at the same temperature.

Parameters Identification
For the second-order RC equivalent circuit model indicated in Figure 1, the value of parameters {U oc (SOC), R o , R 1 , R 2 , C 1 , C 2 } and the relationship between OCV and SOC have to be identified.These parameters can be divided into two parts.The first part is U oc , which has a certain relationship with SOC.The second is comprised of {R o , R 2 , R 2 , C 1 , C 2 }, which describe the dynamic behaviors of battery.To identify these parameters, some experiments have been carried out on the ICR18650-22F lithium-ion battery (Samsung, Seoul, Korea).
To collect data for determining the relationship of OCV versus SOC, a consecutive test was performed on the battery.First, the battery is charged to the fully charged state with the standard charging pattern and fully discharged at the same temperature twice.Then, after 2 h rest, the fully charged battery that has reached steady state is discharged with a constant current of 0.1 C by 10% of the nominal capacity and it is left in the open-circuit for condition for 5 h.Repetitive operation of the above discharging with a constant current and rest is executed until the batter achieves the fully discharged state.The entire process of experiment is presented in Figure 2.
According to the measured data from the above test, sixth-order polynomial fitted curve are shown in Figure 3, which can well describe the nonlinear relationship between the OCV and the SOC: The exponential-function fitting method [28][29][30] The details of parameters identification can be found in [29], which has been published by our research group.
The identified parameters of second-order RC ECM are listed in Table 1.

Cubature Kalman Filter Algorithm (CKF)
In 2009, a new nonlinear filter called the cubature Kalman filter (CKF) [21,22] was proposed by Arasaratnam and Haykin.On the basis of the Bayesian filter, the CKF uses a set of 2n points to capture the mean and covariance of the states with additive Gaussian noise by the three-degree radial-spherical cubature rule.
The basic idea of the Bayesian filter is calculating the multi-dimensional weighted integral, which can be expressed as:

Cubature Kalman Filter Algorithm (CKF)
In 2009, a new nonlinear filter called the cubature Kalman filter (CKF) [21,22] was proposed by Arasaratnam and Haykin.On the basis of the Bayesian filter, the CKF uses a set of 2n points to capture the mean and covariance of the states with additive Gaussian noise by the three-degree radial-spherical cubature rule.
The basic idea of the Bayesian filter is calculating the multi-dimensional weighted integral, which can be expressed as: Table 1.Identified parameters of second-order RC ECM.

Cubature Kalman Filter Algorithm (CKF)
In 2009, a new nonlinear filter called the cubature Kalman filter (CKF) [21,22] was proposed by Arasaratnam and Haykin.On the basis of the Bayesian filter, the CKF uses a set of 2n points to capture the mean and covariance of the states with additive Gaussian noise by the three-degree radial-spherical cubature rule.
The basic idea of the Bayesian filter is calculating the multi-dimensional weighted integral, which can be expressed as: Energies 2015, 8, 13458-13472 where f (¨) is some nonlinear function, D R n is the region of integration, and the known weighting function w(x) ě 0 for all x D. In a Gaussian-weighted integral, w(x) is a Gaussian density and satisfies the non-negativity condition in the entire region.According to [21], it can be obtained that: where N(x; 0, I) is Standard Gaussian distribution, n represents the state-vector dimension, [1] i represents the i-th column vector.The third-degree cubature rule is exact for Gaussian-weighted integrals whose integrands are written in the form of a linear combination of monomials up to the third degree.
A discrete-time nonlinear dynamical system with additive process and measurement noises is as follows: where X k R n is the state vector at time k; y k R m is the measurement vector at time k; f (¨) and h(¨) are known nonlinear functions; Γ k is the discrete-time process noise distribution matrix; w k and v k are independent Gaussian white process noise and measurement noise with covariance Q k and R k separately.The process of the CKF algorithm for battery SOC estimation is summarized as follows.
(a) Initialization: x0 " Erx 0 s, P 0 " Erpx 0 ´x 0 qpx 0 ´x 0 q T s ( (b) Time update (1) Calculate the cubature points: where n represents the state-vector dimension and ξ is the set of standard cubature points, which is shown by: where [1] i represents the identity matrix and [1] (i) denotes its i-th column vector.(2) Calculate the propagated cubature points: (3) Calculate the predicted state and covariance: where Q k is the process noise covariance matrix at time step k.
(2) Calculate the propagated cubature points: (3) Calculate the predicted measurement and covariance: where R k is the measurement noise covariance matrix at time step k.

Strong Tracking Cubature Kalman Filter (STCKF)
Although the EKF has been widely used, it is limited for the following drawbacks: (1) Poor robustness against model uncertainties.
(2) Loss of tracking ability for sudden changes of the state when it has reached steady state.
(3) Cannot be used to estimate time-varying parameters.
In order to improve performance of EKF, strong tracking filter (STF) was proposed [31,32] and used for battery SOC estimation [33].In the actual design, STF still continues to use framework of Kalman filter because of conciseness and recursive property.The framework of Kalman filter can be expressed as: xk where γ k`1 is residual sequence of measurement.Hk is H k " when the system is nonlinear.Hence, STF is transformed into seeking the optimal filter gain matrix K k+1 that can be deduced by orthogonality principle [32][33][34].
When system model and real system completely match and modeling error does not exist, residual sequence of Kalman filter output is white Gaussian noise with zero mean, so these residual Energies 2015, 8, 13458-13472 sequences are mutually orthogonal.Equation (25), which is the performance index of Kalman filter, denotes the minimum error estimation covariance and Equation (26) represents that these residual sequences are mutually orthogonal.Since the system model is not very accurate in fact, the mean of residual sequence will show this problem when state estimation of filter deviates from the system state.Thus, the STF can adjust gain matrix online by introducing fading factor to the state estimation covariance matrix.
where λ k+1 is the fading factor; P x k`1|k is state estimation covariance matrix before introducing fading factor; P k`1|k is state estimation covariance matrix after introducing fading factor; Q k is covariance of process noise.The fading factor can be obtained by following equations: where H k and F k are the measure matrix and process matrix, respectively, and λ k is called the fading factor which adjusts gain matrix to realize orthogonality principle.Q k and R k are covariance of process noise and measure noise, respectively.The matrix N k+1 and M k+1 are defined by Equations ( 30) and (31), which are used for calculating the fading factor λ k ; tr[N k ] and tr[M k ] are the trace of matrix N k and M k , respectively.By combining the above cubature Kalman filter and strong tracking filter, we can deduce the strong tracking cubature Kalman filter (STCKF) algorithm in detail, as shown in Table 2.
In Table 2, after initialization, the estimated state vector xk+1|k and the estimation error covariance P k+1|k can be firstly achieved according to the time and measurement update processes.Then, the fading factor λ k+1 is calculated.The estimated state vector xk`1 and the estimation error covariance P k+1 after adding the fading factor are obtained though update processes.xk`1 , P k+1 and λ k`1 are used for the next prediction and update processes.The battery SOC can be recursively estimated by repeating the above procedure.

Initialization
x0 " Erx 0 s, P 0 " Erpx 0 ´x 0 qpx 0 ´x 0 q T s (1) Time update (a) The cubature points (3) The fading factor (4) Update after add fading factor (a) The cubature points P k`1|k " λ k`1 P x k`1|k `Qk (c) Measurement and error covariance x i k`1|k pγ i k`1 q T ´xk`1|k y k`1 T (d) Estimate the Kalman gain, updated state and error covariance

Experimental Configurations
Figure 4 shows the schematic of battery test bench.It consists of (1) tested Samsung lithium-ion batteries; (2) a control board for controlling battery charge/discharge, as well as battery voltage and current sampling with a period of 1 s; (3) a host computer with monitoring software for data sampling and MATLAB R2010a for data analysis; (4) a DC contactor for charge/discharge switching; (5) a programmable power supply for cell charging; and (6) a programmable electric load for cell discharging.The LiNi x Co y Mn z O 2 (x + y + z = 1) lithium-ion battery, which will be the next generation of the mainstream batteries for EVs, has a distinctive features in energy density, power density, cycle life and so on.The batteries used in this test are Samsung ICR18650-22P lithium-ion battery, whose nominal capacity and nominal voltage are 2.2 Ah and 3.62 V, respectively.and MATLAB R2010a for data analysis; (4) a DC contactor for charge/discharge switching; (5) a programmable power supply for cell charging; and (6) a programmable electric load for cell discharging.The LiNixCoyMnzO2 (x + y + z = 1) lithium-ion battery, which will be the next generation of the mainstream batteries for EVs, has a distinctive features in energy density, power density, cycle life and so on.The batteries used in this test are Samsung ICR18650-22P lithium-ion battery, whose nominal capacity and nominal voltage are 2.2 Ah and 3.62 V, respectively.

Estimation Results under Dynamic Stress Test (DST) Cycle
A good SOC estimation technique should be able to perform effectively under different loading profiles.Thus, we collected battery discharge data using dynamic stress testing (DST) to verify the applicability of the developed method during working conditions, which is a step charge/discharge profile.The current profile and zoom profile for DST are shown in Figure 5a,b, respectively.

Estimation Results under Dynamic Stress Test (DST) Cycle
A good SOC estimation technique should be able to perform effectively under different loading profiles.Thus, we collected battery discharge data using dynamic stress testing (DST) to verify the applicability of the developed method during working conditions, which is a step charge/discharge profile.The current profile and zoom profile for DST are shown in Figure 5a,b, respectively.
11 current sampling with a period of 1 s; (3) a host computer with monitoring software for data sampling and MATLAB R2010a for data analysis; (4) a DC contactor for charge/discharge switching; (5) a programmable power supply for cell charging; and (6) a programmable electric load for cell discharging.The LiNixCoyMnzO2 (x + y + z = 1) lithium-ion battery, which will be the next generation of the mainstream batteries for EVs, has a distinctive features in energy density, power density, cycle life and so on.The batteries used in this test are Samsung ICR18650-22P lithium-ion battery, whose nominal capacity and nominal voltage are 2.2 Ah and 3.62 V, respectively.

Estimation Results under Dynamic Stress Test (DST) Cycle
A good SOC estimation technique should be able to perform effectively under different loading profiles.Thus, we collected battery discharge data using dynamic stress testing (DST) to verify the applicability of the developed method during working conditions, which is a step charge/discharge profile.The current profile and zoom profile for DST are shown in Figure 5a,b, respectively.The estimation curve of SOC is presented in Figure 6a, where the black dotted-line is the reference SOC values computed using Coulomb counting method with accurate current values and initial SOC value, while the red solid-line and black solid-line is the estimated results with STCKF and EKF method, respectively.The initial SOC value for STCKF and EKF method was set at 95% for two reasons: (1) to simulate preferably the actual conditions, where the exact SOC at the beginning of usage is unknown or inaccurate, and (2) to measure the self-correction capability of the STCKF method.The SOC estimation error of STCKF and EKF method is indicated in Figure 6b with the red Energies 2015, 8, 13458-13472 solid-line and black solid-line.From Figure 6c, we can easily get that the value of SOC estimation error with STCKF is a little bigger than the result with EKF method at the beginning.Then, SOC estimation error deceases and approaches a very small value gradually; however, estimation error using EKF method changes little.The SOC estimation error from 4000 to 5000 s is shown in Figure 6c, which indicates that the maximum error of STCKF is less than 1% and the maximum error of EKF is about 3%.The battery terminal voltage versus time is shown in Figure 7a, where the blue solid-line is the measured terminal voltage with a high precision voltage sensor and the red dotted-line is the estimated voltage by the developed method.It can be seen that the terminal voltage shows serious fluctuations due to the sharply variable current.However, the voltage estimation error has small fluctuations and the average value of estimation error is about 0.02 V, as shown in Figure 7b.The estimation results of SOC and terminal voltage illustrate that the proposed method for SOC estimation has good performance in terms of estimation accuracy and self-correction capability.

12
The estimation curve of SOC is presented in Figure 6a, where the black dotted-line is the reference SOC values computed using Coulomb counting method with accurate current values and initial SOC value, while the red solid-line and black solid-line is the estimated results with STCKF and EKF method, respectively.The initial SOC value for STCKF and EKF method was set at 95% for two reasons: (1) to simulate preferably the actual conditions, where the exact SOC at the beginning of usage is unknown or inaccurate, and (2) to measure the self-correction capability of the STCKF method.The SOC estimation error of STCKF and EKF method is indicated in Figure 6b with the red solid-line and black solid-line.From Figure 6c, we can easily get that the value of SOC estimation error with STCKF is a little bigger than the result with EKF method at the beginning.Then, SOC estimation error deceases and approaches a very small value gradually; however, estimation error using EKF method changes little.The SOC estimation error from 4000 to 5000 s is shown in Figure 6c, which indicates that the maximum error of STCKF is less than 1% and the maximum error of EKF is about 3%.The battery terminal voltage versus time is shown in Figure 7a, where the blue solid-line is the measured terminal voltage with a high precision voltage sensor and the red dotted-line is the estimated voltage by the developed method.It can be seen that the terminal voltage shows serious fluctuations due to the sharply variable current.However, the voltage estimation error has small fluctuations and the average value of estimation error is about 0.02 V, as shown in Figure 7b.The estimation results of SOC and terminal voltage illustrate that the proposed method for SOC estimation has good performance in terms of estimation accuracy and self-correction capability.

Estimation Results under New European Driving Cycle (NEDC) Cycle with Voltage Noise
In practice, it is difficult to always obtain the precise measurement of current and voltage, because of the noise interference caused by incorrect current and voltage sensor error and electromagnetic interference.In order to simulate the practical condition, a sequence of voltage noise was attached to measured voltage to verify the robustness of the developed method against measurement noise.The comparison with the frequently-used extended Kalman filter (EKF) and cubature Kalman filter (CKF) can show the characteristics in terms of accuracy and robustness.As is shown in Figure 8a, the absolute value of voltage noise is less than 0.08 V and Figure 8b illustrates more details.
The comparison results of SOC estimation are shown in Figure 9.As indicated in the Figure 9a, the black solid-line is the reference SOC value calculated by Coulomb counting method with accurate current values and initial SOC value, the red solid-line is SOC estimation using EKF method, the blue solid-line shows SOC estimation result by CKF method and the green solid-line represents the SOC estimation computed by STCKF method.From Figure 9b, it is clear that the SOC estimation error calculated by STCKF method has a smaller fluctuation than estimation error using EKF or CKF method.In order to further evaluate merits of the STCKF, the comparison of root mean square error (RMSE), maximum estimation error and execution time with different method are shown in Table 3.It is shown that the RMES of EKF method and CKF method are 0.0233 and 0.0157, respectively, while the RMSE computed by STCKF method is 0.0133.The maximum errors of EKF and CKF are 6.13% and 5.28%, respectively, while that of the STCKF is only 4.17%.As shown in Table 3, execution time of EKF method and CKF method are 0.89 and 1.55 s, respectively.However, execution time of STCKF method is 2.58 s.Therefore, it can be concluded that although the STCKF takes more computational cost, it is more robust to measurement noise compared with the EKF and CKF algorithms.

Estimation Results under New European Driving Cycle (NEDC) Cycle with Voltage Noise
In practice, it is difficult to always obtain the precise measurement of current and voltage, because of the noise interference caused by incorrect current and voltage sensor error and electromagnetic interference.In order to simulate the practical condition, a sequence of voltage noise was attached to measured voltage to verify the robustness of the developed method against measurement noise.The comparison with the frequently-used extended Kalman filter (EKF) and cubature Kalman filter (CKF) can show the characteristics in terms of accuracy and robustness.As is shown in Figure 8a, the absolute value of voltage noise is less than 0.08 V and Figure 8b illustrates more details.
The comparison results of SOC estimation are shown in Figure 9.As indicated in the Figure 9a, the black solid-line is the reference SOC value calculated by Coulomb counting method with accurate current values and initial SOC value, the red solid-line is SOC estimation using EKF method, the blue solid-line shows SOC estimation result by CKF method and the green solid-line represents the SOC estimation computed by STCKF method.From Figure 9b, it is clear that the SOC estimation error calculated by STCKF method has a smaller fluctuation than estimation error using EKF or CKF method.In order to further evaluate merits of the STCKF, the comparison of root mean square error (RMSE), maximum estimation error and execution time with different method are shown in Table 3.It is shown that the RMES of EKF method and CKF method are 0.0233 and 0.0157, respectively, while the RMSE computed by STCKF method is 0.0133.The maximum errors of EKF and CKF are 6.13% and 5.28%, respectively, while that of the STCKF is only 4.17%.As shown in Table 3, execution time of EKF method and CKF method are 0.89 and 1.55 s, respectively.However, execution time of STCKF method is 2.58 s.Therefore, it can be concluded that although the STCKF takes more computational cost, it is more robust to measurement noise compared with the EKF and CKF algorithms.

Conclusions
In this paper, the STCKF algorithm is proposed for SOC estimation of lithium-ion batteries in EVs.The STCKF algorithm is based on the cubature Kalman filter and satisfies the orthogonality principle by introducing fading factor to the state estimation covariance matrix, improving the accuracy and robustness.The frequently-used second-order RC equivalent circuit model is selected to describe the nonlinear dynamic behaviors of lithium-ion batteries to realize the trade-off between

Conclusions
In this paper, the STCKF algorithm is proposed for SOC estimation of lithium-ion batteries in EVs.The STCKF algorithm is based on the cubature Kalman filter and satisfies the orthogonality principle by introducing fading factor to the state estimation covariance matrix, improving the accuracy and robustness.The frequently-used second-order RC equivalent circuit model is selected to describe the nonlinear dynamic behaviors of lithium-ion batteries to realize the trade-off between

Conclusions
In this paper, the STCKF algorithm is proposed for SOC estimation of lithium-ion batteries in EVs.The STCKF algorithm is based on the cubature Kalman filter and satisfies the orthogonality Energies 2015, 8, 13458-13472 principle by introducing fading factor to the state estimation covariance matrix, improving the accuracy and robustness.The frequently-used second-order RC equivalent circuit model is selected to describe the nonlinear dynamic behaviors of lithium-ion batteries to realize the trade-off between computation cost and accurate dynamic characteristics.The exponential-function fitting method and sixth-order polynomial fitting are used to implement parameters identification.The Dynamic Stress Testing with 5% initial SOC error is applied to validate performance of STCKF algorithm.The experiment results of SOC and terminal voltage illustrate that the proposed method for SOC estimation performs well in terms of estimation accuracy and self-correction capability to initial SOC error.Furthermore, comparison with EKF and CKF method under NEDC cycle test with a sequence of voltage noise reveals the good property in SOC estimation accuracy and robustness.However, there are also some issues in the proposed method that need further investigations.Constant model parameters are applied to obtain SOC information and the SOC dependency of the model parameters is not studied.

Figure 1 .
Figure 1.Schematic diagram of the second-order resistor-capacitor (RC) equivalent circuit model.

Figure 1 .
Figure 1.Schematic diagram of the second-order resistor-capacitor (RC) equivalent circuit model.

Energies 2015, 8 , 5 Figure 2 .
Figure 2. Flow chart of an open-circuit voltage (OCV) test for determining the relationship of OCV and state of charge (SOC).

Figure 3 .
Figure 3.The sixth-order polynomial fitted curve between OCV and SOC.

Figure 2 . 5 Figure 2 .
Figure 2. Flow chart of an open-circuit voltage (OCV) test for determining the relationship of OCV and state of charge (SOC).

Figure 3 .
Figure 3.The sixth-order polynomial fitted curve between OCV and SOC.

Figure 3 .
Figure 3.The sixth-order polynomial fitted curve between OCV and SOC.

Figure 5 .
Figure 5. Dynamic Stress Testing: (a) current profile and (b) zoom plot of current profile.

Figure 5 .
Figure 5. Dynamic Stress Testing: (a) current profile and (b) zoom plot of current profile.Figure 5. Dynamic Stress Testing: (a) current profile and (b) zoom plot of current profile.

Figure 5 .
Figure 5. Dynamic Stress Testing: (a) current profile and (b) zoom plot of current profile.Figure 5. Dynamic Stress Testing: (a) current profile and (b) zoom plot of current profile.

Figure 6 .
Figure 6.The estimation results under Dynamic Stress Test (DST): (a) SOC estimation profile; (b) SOC estimation error profile; and (c) Zoom plot of SOC estimation error profile.

Figure 6 . 13 Figure 7 .
Figure 6.The estimation results under Dynamic Stress Test (DST): (a) SOC estimation profile; (b) SOC estimation error profile; and (c) Zoom plot of SOC estimation error profile.

Figure 7 .
Figure 7. Terminal voltage estimation under DST: (a) Voltage profile; and (b) Voltage estimation error profile.

Figure 8 .
Figure 8. Random voltage noise: (a) Random noise profile; (b) Zoom plot of random noise profile.

Figure 9 .
Figure 9.The SOC estimation results under New European Driving Cycle (NEDC) test: (a) SOC estimation profile; (b) SOC estimation error profile.

Figure 9 .
Figure 9.The SOC estimation results under New European Driving Cycle (NEDC) test: (a) SOC estimation profile; (b) SOC estimation error profile.

Figure 9 .
Figure 9.The SOC estimation results under New European Driving Cycle (NEDC) test: (a) SOC estimation profile; (b) SOC estimation error profile.

Table 1 .
Identified parameters of second-order RC ECM.

Table 1 .
Identified parameters of second-order RC ECM.

Table 3 .
Comparison of SOC estimation method under NEDC test.

Table 3 .
Comparison of SOC estimation method under NEDC test.