Online State of Charge Estimation for Lithium-Ion Battery by Combining Incremental Autoregressive and Moving Average Modeling with Adaptive H-Infinity Filter

The state of charge (SOC) estimation is one of themost important features in batterymanagement system (BMS) for electric vehicles (EVs). In this article, a novel equivalent-circuit model (ECM) with an extra noise sequence is proposed to reduce the adverse effect of model error. Model parameters identification method with variable forgetting factor recursive extended least squares (VFFRELS), which combines a constructed incremental autoregressive and moving average (IARMA) model with differential measurement variables, is presented to obtain the ECM parameters. The independent open circuit voltage (OCV) estimator with error compensation factors is designed to reduce the OCV error of OCV fittingmodel. Based on the IARMAbatterymodel analysis and the parameters identification, an SOC estimator by adaptive H-infinity filter (AHIF) is formulated.The adaptive strategy of the AHIF improves the numerical stability and robust performance by synchronous adjusting noise covariance and restricted factor.The results of experiment and simulation have verified that the proposed approach has superior advantage of parameters identification and SOC estimation to other estimation methods.


Introduction
In the field of electric vehicles (EVs), the Lithium-ion batteries have been widely researched and quickly developed in battery management system (BMS) at present [1,2].Onboard BMS is used to guarantee the safety of the work and reliable battery application in high energy EVs.An accurate SOC is one of the core functions in the BMS, which provides the fundamental principle of the battery available capacity for maintaining a longer lasting time.In many cases, other additional functions in the BMS also rely on precise SOC estimation [3,4].Beyond that, the accurate SOC estimation is dramatically impacted on the time scale for charge and discharge; otherwise overcharge and overdischarge will have serious and negative impact on battery life time and EVs driving safety.Therefore, a reliable and high-precision SOC prediction method is significant for the optimization of BMS in EVs control strategy [5,6].
Various efforts focusing on the estimation of SOC have been made in the last few years.In general, the known estimation methods could be classified into two categories: model-free methods and model-based methods.The classical model-free approach is the Ampere-Hour integral (AH) [7,8] that is the remainder charge is equal to the integral of charge or discharge current.Nevertheless, two preconditions should be satisfied simultaneously, which are exact initialization SOC data and the high-precision sensor.In the process of current integration, AH suffers from the accumulative error and rounding error in disturbances condition, such as current or temperature fluctuations, etc.The OCV measurement [9][10][11] also is commonly used technique for SOC estimation based on the nonlinear OCV-SOC fitting curve under incremental OCV test or low-current OCV test.However, it is a poor choice for real-time SOC estimation in EVs ongoing application because the battery's steady state is not achieved until the long rest time.In addition, the machine learning approach (MLA), such as the artificial neural network (ANN) [12,13], the fuzzy logic [14,15], and the support vector machine (SVM) [16,17], takes a Lithium-ion battery as a black-box model.Although MLA does not require the exact current-voltage information of battery dynamics, it is almost impossible to cover all the battery loading data through studying the training dataset.And because of large amount of time consumption, it is not a better design for BMS.
In the model-based estimation methods, the choice of model should be firstly determined.Equivalent-circuit model (ECM) is the most widely used model for battery applications.The ECM can be classed as one reduced form of complex electrochemical model, which is applied to approximate the battery dynamic behavior with identified parameters.In [18], the dynamic response of battery can be simulated by oneorder ECM, and a recursive least squares (RLS) method is used for parameters identification.Other types of ECM which introduce two or three order RC network are presented in [19,20].Although multiorder ECM can lead to better performance in representing battery dynamic behavior than the one-order ECM, the more ECM parameters need to be identified with deepening of algorithm complexity.In [21,22], an additional noise is added to one-order ECM, but the RLS with constant forgetting factor may encounter the difficulties of compromise between stability and convergence.Therefore, the effective ECM is largely dependent on the balance between complexity and precision of corresponding algorithm.
In the ECM-based state observers, the flexible usage of Kalman filter has been generally accepted for its better performance of optimal estimation, which assumes some special properties of the problem formulation, such as model error-free and known noise statistics.For example, the extended Kalman filter (EKF) algorithm is applied for Lithium-ion battery state estimation to obtain a desired resistance estimation result [23,24].But with the strong fluctuating current problems and the undeniable error of the first-order Taylor series approximation, the convergence speed of EKF is slow or even cannot get a convergence solution in time.The unscented Kalman filter (UKF) with unscented transformation (UT) method is utilized to handle the SOC estimation problem in strong nonlinear battery systems [25,26], and the estimation results show that the UKF has better robustness and higher precision than the EKF.The adaptive EKF (AEKF) and adaptive UKF (AUKF) based on the dynamic adaptive strategy of adjusting noise were adopted in state estimator to achieve the goal of highprecision and better stability than EKF and UKF [27,28].In [29][30][31][32][33], an H-Infinity Filter (HIF) based on the minimax residual error criterion observer is designed for model parameters identification and state estimation, which has more robust to model error uncertainty without assumptions of accurate battery model.Generally, the HIFs are more robust than KFs [34,35], but they do not sufficiently consider the effect of the error of process noise and measurement noise.In other words, the uncertain covariance of noise can result in unavoidable estimation error of the HIF.One common feather of the aforementioned methods is that OCV can be calculated from OCV-SOC fitting model [23][24][25][26][27][28][29][30][31][32][33][34][35].However, due to fitting error in OCV-SOC curve [9], any small perturbation in OCV calculation by OCV-SOC fitting model may cause a large deviation in SOC estimation, especially a wide flat area existing in the OCV-SOC curve with a strong nonlinear relationship.
The main contributions of this paper include the following: (1) a novel IARMA battery model based on VFFRELS method is designed to identify the battery model parameters; (2) an error compensation based independent OCV estimator is adopted to accurately capture the OCV; (3) with design and application of adaptive HIF (AHIF) including noise covariance and restricted factor adjustment to the proposed model, the SOC is directly estimated with reasonable performance.Compared with the wide usage of AEKF, the proposed method has the features of accurate model, strong robustness, and high precision.
The remainder of this paper is arranged as follows.Section 2 presents details of the IARMA-based battery ECM and the model parameters identification based on VFFRELS algorithm.The error feed forward compensation based independent OCV estimator is analyzed in Section 3. Then based on the IARMA battery model and identified parameters, an AHIF-based SOC estimator is proposed in Section 4. Section 5 shows the experimental and simulation results under the different conditions.Finally, some conclusions are drawn in Section 6.

Battery Model and Parameters Identification
2.1.Battery Model Derivation.In a simplified model sense, the dynamic electrical characteristic of the Lithium-ion battery can be commonly represented by the one-order ECM [18], as shown in Figure 1.
One-order ECM is expressed as follows: where   () is used to describe the SOC dependent OCV which is associated with OCV-SOC fitting curve,  0 indicates the ohmic resistance, the parallel     composed of a polarization resistance and a polarization capacitance,   () represents polarization voltage across parallel     , () is the current in the charging and discharging status, and   () represents the terminal voltage.The one-order ECM cannot completely describe battery dynamic nature, since the model random noise characteristic is not included in (1).As illustrated in Figure 2, a colored noise   has been added as a compensation module in oneorder ECM, and   is to cover the random error which is not considered in the one-order ECM [21,22].
Based on Figure 2, the mathematical model of one-order ECM with additive colored noise can be derived as In ( 2), the colored noise   () is a noise sequence driven by a series of white noise sequence, and (), (−1), ⋅ ⋅ ⋅ , (−   ) are its regression of data; thus, The nonlinear time series model is adopted to analyze battery ECM; let u and y be the charge/discharge current and the terminal voltage of battery ECM, respectively.The ECM in (2) to be identified can be redefined as a multivariable regression expression: where  and  are the orders of regression expression;  1 ,  2 , ...  and  1 ,  2 , ...  are the undetermined coefficient.

Parameters Identification.
By comparing (2) with (4), the polarization voltage is firstly needed to eliminate from (2).The discretization form of (2) by bilinear transformation method:  = (2/  )((1 −  −1 )/(1 +  −1 )) ( is the discretization operator) is given as follows: where   (),   () and () is the OCV, terminal voltage, and load current at the kth sampling time, respectively, and the coefficient can be obtained by Then  0 ,   , and   can be obtained according to the inverse equations of  1 ,  1 , and  2 ; thus, By comparing ( 5) with (4), the immeasurable part [  () −  1   ( − 1)] is assumed as the residual model error (RME) of autoregression (AR) model.Since the above colored noise sequence   () can be used as a MA process, (5) can be considered as an autoregressive moving average model (ARMA) with external input system.
If the parameters identification algorithm is based on (5) or (7), whether the RME can be ignored compared with terminal voltage in ( 5) or (7) is the first consideration.The detailed property of two RME of ARMA and IARMA will be analyzed as a quantifiable form in Section 5.
From (7), it is seen that  and  of the IARMA model are Δ  and Δ, respectively.In accordance with nonlinear regression mode principle, (7) can be rewritten as leastsquares form: where In (8), () is the unknown parameters set, and () is the coefficient determined by known variable.Since the  time sequence of Δ() is immeasurable, () cannot be got directly; therefore some traditional least squares (LS) series online identification algorithms are unable to identify ().
Accordingly, the RELS algorithm [21,22] is introduced to identify the immeasurable noise terms first.Based on ( 5) and (8), () can be approximated to the difference between measured terminal voltage and identified terminal voltage; meanwhile, Δ() can be approximated to the difference between measured differential terminal voltage and identified differential terminal voltage, as follows.
With the updated φ() by (11), the regression data of colored noise   () can be obtained by ( 5)- (10).The procedures of IARMA and RELS based identification method are illustrated in Figure 3.
The RELS is used to solve the regression model described in Figure 3.However, the RELS with constant forgetting factor  may encounter the difficulties of balancing between stability and convergence if the model parameters change with different rates.The battery dynamics nature of one-order ECM can be distinguished into two slow parameters and one fast parameter [38][39][40].The fast parameter represents internal ohmic resistance  0 , while the slow parameters include   and   .In this regard, a large forgetting factor should be assigned to the parameters changing slowly to guarantee the stability of the algorithm, while a small forgetting factor is more appropriate for the tracking of fast varying parameter.In seeking to address this problem, the VFFRELS with multiple forgetting factors [41][42][43] for identification is applied in this paper.With the VFFRELS, the forgetting factors can be decoupled and tuned separately to improve the parameters stability and tracking accuracy of SOC estimation.

Independent OCV Estimator with Error Compensation
3.1.OCV Estimator.In this section, an independent OCV estimator is introduced.The OCV is observed with the aid of battery ECM characterization.  () in (1) can be rewritten in the discrete-time form as follows: From ( 1), obviously   () can be expressed in the following discrete-time relationship as follows: Substituting ( 13) into (12) yields the following: Since the OCV can be described as a slow time-varying variable, then the OCV estimation can be solved by the expression of  0 ,   , and   , as follows:

OCV Error Compensation.
The OCV-SOC curve can be easily obtained under incremental OCV test and low-current OCV test [35].The nonlinear characteristic relationship of OCV-SOC fitting model is built in processing the experimental data by the polynomial curve fitting methods.
where   is the polynomial fitting coefficient, N is the order of fitting model, and N is set to 8. The direct calculation of OCV by OCV-SOC fitting model contains perturbations due to the existence of fitting error.Therefore, the OCV error will largely decline the accuracy of SOC estimation.Besides, there exists a flat area in the OCV-SOC curve with a strong nonlinear relationship when the value of SOC varies from 0.2 to 0.8 [18].Even if the measurement error of voltage or current is relatively slight, the OCV error can lead to a relatively large OCV-based SOC estimation error.To solve this problem, OCV is corrected with two compensations by considering the intrinsic OCV variation.
First, the residual error of terminal voltage is introduced to approximately describe the behaviors of deviation of OCV.A microadjustment factor is adopted for compensating error caused by nonequivalence behavior between terminal voltage and OCV.The simple correction formula is shown as follows: where Δ V1 is used as first compensation, microadjustment factor  1 is below a certain range (0 <  1 < 1 in this paper), and the proportion coefficient  1 is 10 2 .The microadjustment factor  1 is determined by the slope of OCV-SOC curve.
Next, due to the fitting error in OCV-SOC fitting model, the second appropriate compensation is necessary to further eliminate OCV error.The second compensation Δ V2 is adopted to reduce the fitting error in OCV-SOC fitting model.Δ V2 for OCV is equal to the behaviors of the deviation between V () by OCV estimator and   () by OCV-SOC fitting model.
The block diagram of two OCV compensation methods is shown in Figure 4.The terminal voltage estimation value ( V ) and OCV estimation value ( V1 ) are utilized to the get corresponding compensation factor (Δ V1 and Δ V2 ) which are employed to correct the OCV by the feed forward method.

Adaptive HIF Based SOC Estimator
4.1.HIF with Tunable Restricted Factor.In the strict sense, the accurate ECM and noise statistics quantity is the precondition to bring into playing KFs advantages.However, the KFs cannot be guaranteed to work stably all the time when these conditions cannot be satisfied.The H-Infinity Filter (HIF) based on the game theory which is designed to minimize the maximum estimation error can make up for the lack of the precondition; thus it is common in resolving the problem of system identification in [29][30][31].The SOC estimation in [32][33][34][35] shows that the HIF can provide the impressive results with inexact or unknown initialized state.The cost function of HIF is defined as follows: where   is the state value and x is state estimation;  0 and x0 are the initial states value.  and V  are the process noise vector and measurement noise vector, respectively. 0 ,   ,   , and   are the symmetric positive definite matrixes.Further, since it is difficult to directly minimize cost function  1 , a user-specified restricted factor  2 which has important influence on the precision and robustness is presetting to guarantee an optimized boundary constrained condition, making sure that  1 satisfies  1 <  2 .In [32,33], the restricted factor  2 is set to a fixed value based on the engineering practice experience, and three restricted factors of 1.531, 1.048, and 1.003 are adopted in state estimation.For ensuring both robustness and precision, the self-adapting restricted factor is incorporated into HIF by innovation mechanism.
By applying the matrix inversion lemma to  | , we obtain As  | is a covariance matrix, it therefore should be positive definiteness.This means that In other words, the restricted factor  should satisfy where  max () −1 denotes the greatest eigenvalue of matrix () −1 ; therefore, self-adapting method for restricted factor  is where  2 is used to correct  and  2 > 1; the quadratic sum of innovation error   from ( 27) is used as estimation error      .As the restricted factor  is inversely proportional to the estimation error, the correction coefficient expression is shown as follows: where  2 is an unknown coefficient associated with experiments and  2 > 0, once  2 is determined, the correction coefficient  2 is only dependent on innovation error   .
That is, the relationship between restricted factor  and estimation error      can be quantitatively depicted by ( 22)- (23).

HIF with Adaptive Strategy. Multiple adaptive ways have integrated into original
KFs such as the Bayes method, maximum likelihood criterion, correlation method, and covariance matching method.As mentioned before, the HIF approaches to KF when  → ∞.Hence it can be concluded that the noise covariance matrix   and   in HIF play the similar role in the KFs which should be adjusted properly.Therefore, an adaptive HIF with the maximum likelihood criterion [34] is designed to update the noise covariance at each stage of measurement correction to strengthen its jamming immunity of SOC estimation.
According to (2), the innovation error can be written as an expression of difference between measurement estimation and true measurement value: The measurement innovation covariance is deduced through the use of probability theory.
where    denotes the covariance of   , then the definition of innovation estimation variance P  is introduced to describe    through using moving window method of innovation.
By combining ( 28) and ( 29), the following observation noise covariance equation can be drawn: Similarly, the state noise   can be expressed as follows: Substituting the prior state estimation x|−1 into (31) yields Based on the principle of orthogonality between the innovation and the residual, the state noise covariance Q is taken on both sides of (32).
Based on the adaptive strategy analysis, the adaptive iteration approach which is detailed in (30) and ( 33) is utilized to compensate for process and measurement noises uncertainties in HIF algorithm.

SOC Estimator
Based on the AHIF.First, the ECM parameters and OCV are identified as two independent parameters set by using VFFRELS and compensation method in improved one-order ECM.Next, the HIF accompanied with tunable restricted factor and adaptive strategy is applied to SOC estimation.The formulas of state and observation are as follows: (34) where  = [,   ],   is the terminal voltage,   represents a linear combination of the state, and the user-defined matrix   = [1, 0].The system matrixes   ,   ,   ,   which describe the dynamics of the ECM are given as follows: With the equations of state and observation above, the battery SOC is estimated by using the improved AHIF.The process of the HIF with tunable restricted factor for the SOC estimation can be seen in ( 24)-( 26), the noise covariance adaptive adjusted strategy is given in (30) and (33), and the general flowchart of the IARMA-AHIF based parameters and SOC joint estimation method is illustrated in Figure 5.The SOC estimation is obtained by recursion procedure based on identified ECM parameters and OCV.

Results and Discussion
To evaluate the performance of parameters identification and SOC estimation, the 20Ah/24V Lithium-ion phosphate battery is selected as the test objects, in which the cathode and anode are formed from Lithium iron phosphate and graphite, respectively.In order to acquire high-precision experimental data such as charge/discharge current, voltage, available capacity, etc., the battery experimental platform is set up for the validation experiments.The reliable test platform is composed of a host computer which is used for data collection and data processing, a high precise power battery test system (Arbin Corp, EVTS) with the maximum charging/discharging current of 300A and maximum voltage  of 400V, and a programmable temperature chamber to provide controllable temperature.

Mathematical Problems in Engineering
The data of two operating conditions including Dynamic Stress Test (DST) and Federal Urban Driving Schedule (FUDS) at a constant temperature 25 ∘ C for the battery are collected to evaluate the performance of model parameters identification and SOC estimation.In addition, the statistical index such as maximum absolute error (MAE) and average absolute error (AAE) are used to qualify performance of the identification and estimation algorithms.

Model Parameters Identification Results and Comparative
Analysis.First, the Lithium-ion battery is fully charged by standard charging method, and then it is left to rest for one hour before being discharged in working condition.The discharge current-time distribution of the DST and FUDS cycle in which the duration of uninterrupted working condition is 35000 seconds is shown in Figure 6.With the data flows of collected discharge current and terminal voltage, the VFFRELS algorithm and IARMA model are combined to identify the ECM parameters and terminal voltage.The identification results of  0 ,   , and   are shown in Figure 7.
As shown in Figures 7(a) and 7(b), the identification values of ECM parameters is able to converge to stable values rapidly in the DST and FUDS cycle, respectively.Among them, the ohmic resistance  0 and the polarization resistance   exhibit the similar tendencies except for a little fluctuation.Specifically,  0 keeps a higher stability and for the ohmic resistance is equal to the ratio of terminal voltage variation to transient current when the variable current is turned off.By contrast to  0 and   , the polarization capacitance   varies significantly because   and   are direct correlation to intricate electrochemical activity.Although the introduction of multiorder ECM can lead to better stable   and   in representing battery dynamic characteristics than the oneorder ECM in theory, the more ECM parameters need to be identified with increasing of algorithm complexity which is not suitable for real-time application.To summarize,  0 and one-order     of ECM have to real-time change with limited amplitude during operating condition to keep closer track of electrochemical activity.
With the identified  0 ,   , and   , the terminal voltage as ECM observation is recursive obtained in each sampling period.Figures 8(a IARMA ECM and the ARMA ECM can both stably track the measurement value of the terminal voltage.The estimation error of the IARMA ECM is obviously decreased compared with the ARMA ECM.As summarized in Tables 1 and 2, the maximum terminal voltage error (MTVE) of the ARMA ECM is up to 109mV and 138mV under the DST and FUDS cycle, respectively; correspondingly the MTVE of the IARMA ECM is only 18 mV and 25mV under the DST and FUDS cycle, respectively.The average terminal voltage error of the IARMA ECM is only 0.35 mV and 1.3mV under the DST and FUDS cycle, respectively, which can be neglected regarding the sensors precision of the EVTS.Besides, from ( 5) and ( 7), it can be inferred that the assumed RME indicates the effect of the OCV in the ARMA ECM is (1 −  1 )  ().Similarly, the RME which defines the distribution of OCV with incremental analysis method in the IARMA ECM can be represented as Δ[  () −  1   ( − 1)].As the simulation results in Figure 10, the RME value  in IARMA ECM is far smaller than that in ARMA ECM.Therefore, the proposed IARMA ECM can provide much better performance in parameters identification.
According to the above analysis, these comparing estimated results of the RME and terminal voltage error fully verify the capacity of the IARMA ECM to describe battery electric property.And the good performance further verifies the authenticity of the parameters identification strategy in practice.

SOC Estimation Results and Comparative Analysis.
Two operation conditions including DST and FUDS are applied to evaluate the performance of SOC estimation method at 25 ∘ C. The reference SOC is calculated by coulomb counting (CC) method in EVTS and the reference SOC is used to approximate true value.As the testing battery in non-lab operation condition, it is difficult to obtain an accurate initial SOC value; thus, it is assumed that the initial SOC is set to 90% (with 10% initial error).
The two different models are compared in regard to the estimation result and error in Figure 11 Figures 11(a)-11(d), where the black dashed-line presents the reference SOC with CC method and the red dotted-line represents the SOC estimation with ARMA-AHIF algorithm, while the blue solid-line indicates the SOC estimation with IARMA-AHIF algorithm.Among them, Figures 11(a To further explore the promoting effects by error compensation strategy, the IARMA-based AHIF and corresponding error compensation (ΔOCV) are utilized to estimate the battery SOC and their tracking results are shown in Figures 12(a)-12(d), respectively, where the black dashed-line shows the reference SOC and the red dashed-line represents the SOC estimation with single AHIF algorithm, while the blue solid-line means the SOC estimation using AHIF-ΔOCV.Among them, Figures 12(a error significantly in comparison with single AHIF, in which the AAE of SOC estimation is decreased to 0.11% under DST cycle and to 0.18% under FUDS cycle, respectively.It can be concluded that the error compensation strategy is conducive to improving SOC estimation precision in AHIF algorithm. Combining IARMA model and error compensation method in AHIF estimator, strengthening the robustness of SOC estimation is even more obvious.The AHIF-ΔOCV is compared with the AEKF-ΔOCV in terms of estimation precision and convergence.The comparison results, including tracking trajectory and estimation error under DST and FUDS cycles, are shown in Figure 13, where the black dasheddoted-line shows the reference SOC, the red dashed-line represents the SOC estimation using AEKF-ΔOCV, and the green doted-line denotes the SOC estimation using single AEKF, while the blue solid-line means the SOC estimation using AHIF-ΔOCV.As indicated in Figures 13(a) and 13(c), the SOC estimation results and reference SOC are subjected to DST and FUDS profile, respectively.The SOC estimation errors by those three algorithms are drawn in Figures 13(b) and 13(d).The estimated SOC by AHIF-ΔOCV tracks the reference SOC with a smoother path than by AEKF, which reflects the reliable correction capacity of the error compensation.The simulation results show that the three algorithms can always track the SOC reference trajectory after the convergence.Moreover, the comparison estimation results of the three algorithms mentioned above are summarized in Tables 3 and 4. From the new perspective on combining IARMA model and error compensation in AHIF and AEKF, an improvement of SOC estimation accuracy plays a certain significant role in the SOC estimation by Kalman filter family state observers.

Conclusion
As we all know that battery model uncertainty and OCV error have negative impacts on the battery SOC estimation.To overcome these shortcomings, the improved one-order ECM   OCV-based SOC estimator is sensitive to the change of OCV-SOC curve.Therefore, the OCV-SOC will be considered as a normalized function including temperature and SOH in the future work.The OCV-SOC normalized function effectively updating is the key factor in the improvement of the precision and robustness of SOC estimator.and Fundamental Ability Enhancement Project for Young and Middle-Aged University Teachers in Guangxi Province (2017KY0865).

Figure 2 :
Figure 2: Schematic diagram of improved one-order ECM.

Figure 3 :
Figure 3: Parameters identification algorithm for the improved one-order ECM.

Figure 4 :
Figure 4: The block diagram of OCV error compensation algorithm.
Figures11(a)-11(d), where the black dashed-line presents the reference SOC with CC method and the red dotted-line represents the SOC estimation with ARMA-AHIF algorithm, while the blue solid-line indicates the SOC estimation with IARMA-AHIF algorithm.Among them, Figures11(a) and 11(c) compare the observed SOC with the reference SOC and Figures 11(b) and 11(d) show the SOC estimation error in which the MAE is less than 1% except at the beginning of discharge.It shows that the two models based estimated SOC can both fast track the reference SOC at the inaccurate initial SOC.In particular, it is clear that the estimated SOC by proposed IARMA-based estimator has lower error in comparison with ARMA-based estimator since the former model can more accurately estimate the terminal voltage, and the difference between the measured and estimated terminal voltage causes a correction term to eliminate the SOC estimation error in a feedback control loop.In short, it is accordingly demonstrated that the IARMA-based AHIF

Table 1 :
Terminal voltage estimation under DST cycle.

Table 2 :
Terminal voltage estimation under FUDS cycle.

Table 3 :
Comparison for the results under DST cycle.

Table 4 :
Comparison for the results under FUDS cycle.These methods include the following: (1) the RELS with variable forgetting factor for parameters identification is applied in IARMA model; (2) the independent OCV estimator with error compensation is adopted to accurately capture the OCV; (3) the high-precision SOC estimation is recursively obtained with improved AHIF including noise covariance and restricted factor adjustment.A high-precision experimental platform has been established to gather reliable data of charge/discharge current and terminal voltage.Two typical operation conditions (DST and FUDS) are adopted to evaluate the performance of parameters identification and state estimation.The simulation results indicate the favorable performances whose average SOC estimation error is only 0.11% (DST) and 0.18% (FUDS) by the proposed integrated method.As a model-based approach, the proposed based on IARMA analysis, independent OCV estimator, and adaptive H-Infinity Filter are proposed for SOC estimation.