State of Charge Estimation of Lithium-Ion Batteries Based on an Improved Sage-Husa Extended Kalman Filter Algorithm

<jats:p>An improved Sage-Husa extended Kalman filter (SHEKF) algorithm is intended to improve the accuracy and stability of SOC prediction. In this paper, two different exponential weighting algorithms are used to adaptively select the forgetting factor for adaptive noise estimation. Moreover, the OCV-SOC curve is obtained using a 7-segment linear fitting method before the algorithms estimate the SOC. In addition, by combining this improved method with a third-order RC equivalent circuit model in the dynamic stress test (DST) case the convergence time is reduced by 0.15 s compared to the second-order RC equivalent circuit model. Following that, four different types of comparison experiments are carried out by comparing the improved algorithm to EKF and other SHEKF algorithms.The estimation accuracy under DST conditions of 0 °C, 25 °C and 45 °C is approximately 0.5%, 2.2% and 1.3% improvement compared to the EKF algorithm.</jats:p>

section_content: New energy vehicles are gradually replacing conventional internal combustion enginepowered vehicles as a result of the growing energy crisis and environmental pollution [1]. One of the most promising new energy vehicles is electric vehicles (EVs), which have energy storage systems made up of several lithium-ion batterie packs. These days, popular research topics include the prediction of lithium-ion battery internal failures [2], the prediction of cycle life [3], the prediction of technology for the various stages of the battery life cycle, including production, use, secondary use, and material recycling [4], and the prediction of battery charge state. The goal of the research presented in this paper is to precisely predict the charge state of lithium-ion batteries. One of them, to ensure the stable operation of the battery management system (BMS), an accurate calculation of the SOC value is necessary. One of the most important things the BMS does is check the state of charge (SOC), which shows what percentage of the nameplate capacity is actually being used [5]. A precise estimate of SOC makes it possible to determine the amount of power that is still available in real time, obtain an exact picture of the state in which the battery is working, and select the charging strategy in a reliable manner. Additionally, the SOC can only be calculated using the pertinent factors because it cannot be directly measured [6]. Simultaneously, several variables, including the algorithm, the model, and the working environment, have an impact on the accuracy of estimation. 

The open-circuit voltage (OCV) method [7], the ampere-hour (Ah) integration method, the Kalman filter (KF) and its improved algorithms [8], the fuzzy logic algorithm [9], and numerous additional potential SOC estimation methods. One of the most popular ways is the OCV method, but one of its main drawbacks is how long it takes to get the OCV of batteries. As a result, there is a disadvantage to this method when it comes to using OCV in practice. Additionally, although the Ah integral approach can directly determine SOC, it is challenging to apply alone since it requires the initial SOC value and high accuracy current measurement. The simplicity, clarity, and ease of understanding of these two approaches are advantages, but accuracy is a drawback that is hard to ensure. 

Methods that are based on the model [10] constitute an additional subfield of strategy. Those methods rely heavily on the battery model [11]. Therefore, it is essential to create a reliable and accurate battery model [12]. Another, those methods can deal with unpredictability and disturbances by self-correcting within a closed-loop system, which improves the precision of estimating the SOC. The noise is unknown in reality, contrary to what the typical EKF or the unscented Kalman filter (UKF) expects. Numerous adaptive filtering algorithms that enable the estimate of noise have been thoroughly researched for battery SOC estimation to address this weakness. A filtering technique called adaptive filtering has the effect of suppressing filter divergence [13]. The Sage-Husa algorithm is a fine example of an algorithm for the measurement noise R and uncertain system noise Q. Although this algorithm is useful, it also has flaws such as excessive computing effort caused by ongoing noise growth and simple filter scattering. 

Under the circumstances of abrupt current shifts, the Sage-Husa algorithm is prone to lose the positive characterization of the covariance matrix, which results in filtering scattering. As a consequence of this, Lv et al. [14] demonstrated that the Sage-Husa method is unable to independently estimate the measurement noise and the system noise, and they also examined the reasons for the algorithm's filter scattering. Measurement noise R is thought to be known, whereas system noise Q is unknowable. Additionally, it is suggested that addition and subtraction operations be omitted from the corrected estimation of system noise Q in favor of just using the matrix's multiplication and inverse operations. Moreover, the system noise needs to be guaranteed to be positive. Liu et al. performed a similar reduction of the system noise equation for the Sage-Husa [15] method. The distinction is that only an adaptive prediction of measurement noise is carried out in this research, which presupposes that system noise is known. Generally speaking, it is possible to minimize the complexity of the covariance matrix, but doing so would decrease the estimation accuracy. According to Li et al. [16], the system and measurement noise covariance matrices were both diagonalized. Although accuracy and stability have improved somewhat, the algorithm still has an average error of about 3. 38 %. It demonstrates that there is room for the accuracy of this algorithm to be improved. To evaluate the gyroscope attitude angle, Jie et al. [17] used the singular value decomposition (SVD) Sage-Husa unscented Kalman filter (UKF), which deconstructs the derived filter covariance matrix P into its parts. The enhanced Sage-Husa UKF technique guarantees the accuracy standards and prevents mistakes from dispersing with time accumulation. Both the SVD and the square-root (SR) decomposition can increase the algorithm's robustness and ensure that the covariance matrix is positive and definite. However, their algorithm complexity has significantly grown, which has resulted in a longer convergence time. In addition, Wang et al. [18] used a technique called covariance matching, as well as the adaptive filtering method developed by Sage-Husa. With this approach, computation is minimized while real-time performance is increased. Nevertheless, the accuracy of the algorithm can still be enhanced. 

To tackle the above problems, the following enhancements are made in this paper. Section 2 describes the modified second-order RC model. Section 3 selects distinct forgetting factors for the two noises according to the condensed Sage-Husa algorithm. Both system and measurement noise are taken into account while performing adaptive estimate. To verify the correctness and stability of the improved algorithm, four types of methods are described in Section 4 of the paper. Finally, in Section 5, it is concluded that the instability problem of Sage-Husa algorithm is solved and the accuracy of the algorithm is increased.
section_name: Introduction

section_name: Equivalent Circuit Model and Parameter Identification

section_content: From the literature [19, 20], it can be seen that considering the influence of temperature on the battery, the electro-thermal model will become the focus of research. If temperature is not a consideration, the equivalent circuit model (ECM) is the best choice. Additionally, it stimulates the dynamic properties of the battery using a resistor-capacitor and voltage source, which has the important benefits of a straightforward principle and simple parameter identification. Because of its accuracy, complexity, and versatility, the second-order RC equivalent circuit model has become a widely used instrument in the field of academic research. However, the addition of the RC network to the third-order model serves to further highlight the effect of the cell's polarization impedance on its performance. In references [21, 22], the benefits and drawbacks of various equivalent models are enumerated. Figure 1 depicts the model of the equivalent circuit. The criteria for the onboard system are extremely strict and challenging to implement, while the higher-order RC network model is excessively complex. Here, the third-order model is adopted. 

Battery pack Equivalent Circuit Model In Figure 1, the U oc is used to indicate the open-circuit voltage. U stands for the terminal voltage of battery. R 0 is a depiction of the battery's ohmic resistance. The Ref. [23] provides a comprehensive explanation of the first two sessions of the RC. Additionally, R 3 and C 3 describe the concentration difference polarization resistance and capacitance, respectively. The absence of the RC network that would have represented the polarization resistor capacitance is the only distinction between the third-order model and the secondorder model. Due to a lack of available space, only the third-order model's depiction will be shown below. The third RC network is only eliminated in the second-order model. The RC model of the third-order produces the following mathematical relationship:

where U1, U2 and U3 denote the derivative of U 1, U 2 and U 3 respectively. The terminal voltage can be expressed as follows:

where the voltage reading at R 1 and C 1 's ends is the U 1. And the voltage reading at the endpoints of R 2 and C 2 is known as the U 2. Particularly, the voltage reading at R 3 and C 3 's ends is the U 3. Following are the components of the third-order model that must be determined. Where θ is the variable that stands in for the parameter that has to be identified.
section_name: Equivalent Circuit Model

section_content: When performing parameter identification on the model, two types of methods are available: online parameter identification [24] [25] [26] and offline parameter identification. In this paper, the latter will be used for parameter identification of the model with three RC links. By extracting features from the pulse discharge experimental voltage fluctuation curve, the model parameters may be determined (as in Figure 2 ). The specific steps are listed below. • Identification of R 0 : R 0 is mainly composed of electrolyte resistance. Once the discharge current is performed or stopped, the terminal voltage of the battery model presented in Figure 1 decreases or increases. According to the discharge-static experiment, in the first moment of battery resting and the last moment of applying discharge current, the instantaneous drop of terminal voltage (U a drops rapidly to U b, U c rises sharply to U d ) reflects the ohmic characteristics of the battery. Hence, the ohmic resistance R 0 is calculated by the Equation (4). 

• Parameter identification of RC parallel circuits: To begin, it is necessary to establish the values for the time constants τ 1, τ 2, and τ 3. And then based on the determined time constants, the parameters of R 1, R 2, R 3, C 1, C 2 and C 3 are determined in detail. Furthermore, crucial for identification is the response of a first-order RC circuit classically consisting of a resistor R, a capacitor C, and a constant current I, which is given by Equation ( 5 ) below:

•

The previous equation states that τ = RC, and t 0 is the time that represents the beginning of the process. The following is the specific process to identify the parameters R 1, R 2, R 3, C 1, C 2, and C 3. Firstly, identify the time constants during the process U b -U c -U d -U e. Take note that there is no current flowing at all during the relaxation stage of the procedure. Then, the voltages U 1, U 2, and U 3 can then be computed using the following calculation based on Equation (5):

We can get Equation ( 7 ) from the output Equation (2). 

Additionally, rewrite Equation (7) as Equation (8). 

where t c signifies the next period. And α 1, α 2, α 3, α 4, β 1, β 2, and β 3 are unknown coefficients, α 1 = U b (∞) which is measured at the end of the polarization recovery process. The ideal coefficients α 2, α 3, α 4, β 1, β 2, and β 3 can be produced by using the MATLAB function "Custom Equation" in the Curve Fitting Toolbox. Therefore, the time constants τ 1, τ 2, τ 3, and the voltages U 1 (t), U 1 (t), and U 1 (t) can be identified. The parameters R 1, R 2, R 3, C 1, C 2 and C 3 can then be determined using Equation (9). 

• Parameter identification of OCV-SOC curve: The steady voltage measured at both ends of the positive and negative electrodes when the battery is left in one location for an extended period of time is the open-circuit voltage (OCV). However, the functional connection between OCV and SOC is not linear. The analytical complexity of SOC estimation is increased by this nonlinearity. In most publications, polynomial fitting is used to obtain OCV-SOC curves. The precision of the SOC estimation is directly effected by the fitting order chosen because it determines the quality of the fitting effect. To enhance the precision of SOC estimation, consider the functional relationship between OCV and SOC as a piecewise linear relationship [27], and refine the fitting curve. The test points are divided into seven segments, each of which is a function of SOC in relation to OCV. The following linear Equation ( 10 ) can be used to represent each component. Table 1 shows the values of k1 and k 2 that correspond to each curve. The 8th-order polynomial fitting curve was compared to the 7-segment linear fitting curve. As shown in Figure 3, the 7-segment linear fitting outperforms the 8th-order polynomial fitting at the start and finish of the discharge.
section_name: Parameter Identification

section_content: Kalman filter, extended Kalman filter and its improved unscented Kalman algorithms need to be used when their system noise and measurement noise are known and constant from beginning to end. Nevertheless, the SOC of the battery is impacted by a wide variety of unknown elements, including temperature or the conditions of the road, while actual driving of electric vehicles takes place. Consequently, noise is something that needs to be determined in an ongoing manner. Sage-Husa is one of the adaptive algorithms. It has a mediocre unbiased estimator of significant posterior noise. The algorithm's core consists of dynamically estimating the statistical features of measurement noise and system noise in actual time utilizing data from measurement variables during the filtering process. Subsequently, the state parameter calculation for EKF makes use of the estimated noise parameters. The system's estimated final state is calculated.
section_name: Sage-Husa Extended Kalman Filter Algorithm

section_content: Since Equations ( 1 ) and ( 2 ) in Section II should always be modified because the battery is a nonlinear system [28]. The rewritten equations of state and observation are presented in the following set of equations. 

where A stands for the state transfer matrix, B represents the state control matrix, H refers to the observation transfer matrix, and D describes the observation control matrix. The system noise matrix is denoted by w k, and the measurement noise matrix is denoted by v k. 

where the measurement noise covariance matrix is denoted by R k, and Q k embodies the system noise covariance matrix. q k represents the mean of w k. r k indicates the mean of v k. 

The SHEKF algorithm adds to the basic equations of the discrete extended Kalman filter [29] the equations for real-time updating of R k and Q k. Below are the formulas that have been added. 

where the forgetting factor is denoted by the letter b, and the weighting factor is characterized by the letter d. And the error covariance is represented by P k. V k represents the difference between the terminal voltage measurement and the estimated value. After that, a short assessment of the Sage-Husa algorithm is provided below. 

1. Write the recursive equation from Equation ( 16 ):

From basic mathematics, it is known that d is gradually decreasing and approaching 1b as k keeps increasing. Therefore, the predetermined initial value Q 0 is the quantity whose contribution to the estimated quantity Q decays fastest as k increases. In addition, the weight of the current time correction is the largest in the estimation of Q k. Namely, the magnitude of the effect that the present instant has on the value of Q is a scaled version of the value of 1b. 

2. With the assumption that the measurement noise's predicted value is, the observation equation can be written as:

rk is calculated as follows:

A benchmarking with Equation ( 13 ) exposes that the final estimate is replaced by a one-step prediction during the estimation of the measurement noise expectation, and thus this is a suboptimal algorithm. Expectation estimation errors may gradually build up as a result of the prediction of the system noise and measurement noise expectation being suboptimal. Additionally, it causes a significant bias in the particular iteration. Even the evaluation of measurement noise or the variance in system noise may be impacted. Therefore, this component of the calculation ought to be disregarded. This is consistent with the study that was conducted and documented in the literature [30]. 

To sum up, for statistical computing of noise, the Sage-Husa extended Kalman filter algorithm is improved over the standard Kalman filter approach. However, the conventional Sage-Husa algorithm is prone to filtering divergence and instability under conditions of severe current variations. The covariance matrix loses its positivity as a result, which is the cause. The SHEKF algorithm enhancement concept was put up to address this issue.
section_name: SHEKF Algorithm

section_name: Improvement of the SHEKF Algorithm

section_content: Equation (16) reveals that in the first half term (1d k ) Rk-1, as long as the initial value of Rk is positive definite, the multiplication result is also positive definite. This is because the range of the forgetting factor d is limited. To update the Rk to be calculated, the second half term

When the filtering process is convergent, the error covariance P k will decrease, and if H is a specific value, HP k,k-1 H T will gradually decline and tend to zero. The noise covariance estimate barely makes a difference. The first half of Equation ( 17 ) has the same basic architecture as Equation (16). Additional, since A has a specific value as well, the two equations can be processed in the same way. In summary, the stability of the filtering is chosen in this paper at the expense of the measurement noise and system noise covariance matrix estimation's objectivity. The adjusted noise equation is displayed by Equations ( 20 ) and (21).
section_name: Simplifying Noise Covariance Matrix

section_content: Because of the emphasis on the role of recent data in time-varying noise statistics, old data should be gradually forgotten, which can be accomplished by using an exponentially weighted decay method. In comparison to other adaptive extended Kalman filtering (AEKF) methods, the SHEKF algorithm includes a forgetting component. SHEKF can change the weights of the inputs at different times due to the forgetting factor, reducing their memory size. This prevents inaccurate and out-of-date data from influencing the filter and ensures that the most recent observations are fully utilized. The forgetting factor is usually represented by the symbol b, and its value varies depending on whether the system error becomes faster or slower, resulting in a different value. A value of b = 0 indicates complete forgetting, whereas a value of b = 1 indicates no forgetting at all. The lower the b, the greater the loss of knowledge. 

In this paper, depending on how quickly or slowly the statistical characteristics of the system noise and measurement noise change, we select various exponential weighting methods. This is due to the fact that model error is the primary source of system noise. Solving this issue is necessary if the system is to be reliable. The fading memory exponential weighting method is employed due to the slow rate of change in the statistical characteristics of the system noise. The fading memory still retains all of the historical data from the past, but the weighting factor and impact on the total weighted sum decrease with age. To prevent being replaced by new data information, it is advised to forget the old data that contributes to system noise for a longer period of time. As a result, the system noise forgetting factor's starting value can be decreased. 

Equation ( 23 ) is generated by substituting Equation (22) into Equation (20). 

And because measurement noise is primarily a product of its surroundings, its characteristics are poorly defined. If the value of the next moment is always very close to the value of the preceding moment, the filter instability caused by significant noise changes can be reduced. When the measurement noise changes quickly or abruptly, the filtered difference data that is too old has a small role in estimating the system error at the current time. When the sampling time is long, it is too old due to the relationship between the forgetting factor and the weighting coefficient. The weight assigned to the filtered difference data of s is very small, which not only has little effect on the estimation of the current moment, but also increases the calculation time due to the frequent and repeated weight factor assignment at each moment. To improve the estimation accuracy of time-variant measurement noise, on the one hand, and to reduce calculation time, only use data from a fixed period of time before the current moment to estimate the error, that is, apply the limited exponential weighting method to SOC estimation. The limited memory index weighting method must take into account the memory length limitation. The memory length in this paper is set to m, where m is a fixed natural number. That is, starting with the most recent data and moving back m generations, the newer data set is multiplied by a large forgetting factor, the older data set is multiplied by a relatively small forgetting factor, and the data before the m generation data is multiplied by a large forgetting factor in this m generation data. The data is completely erased, and a new matrix is created. 

Equation ( 25 ) is obtained by substituting Equation (24) into Equation (21). 

In order to obtain better results, this paper also makes certain differences in the initial values of the two forgetting factors. It is still based on the speed of change in the statistical characteristics of noise. According to the characteristics of the exponential function, the initial value of b 1 greater than the initial value of b 2. 

As a result, Figure 4 depicts the summary flowchart of the modified SHEKF algorithm. 

[ ]

(1 ) / ( 1 )
section_name: Setting Two Improved Forgetting Factors

section_content: The battery research team at the Center for Advanced Life Cycle Engineering Research (CALCE) at the University of Maryland provided the test results for the batteries used in this work [31]. The literature [32] provides a detailed description of the battery test system utilized by CALCE to get these statistics. The battery model was INR18650 20R, and it had a 2000 mAh capacity rating. Additionally, the pulse current multiplier was 0. 5C and the test ambient temperature was 25 °C. After receiving the data on the discharge of the battery, five different algorithms were compared. 

Among these, SHEKF1 is an example of the revised version of the Sage-Husa algorithm that was presented in this paper. The fact that there is just one forgetting factor distinguishes SHEKF2 [33] from its predecessor, SHEKF1. In addition, the Sage-Husa algorithm is represented by SHEKF3 [34], which includes an adaptive estimation of the system noise Q in situations where the measurement noise R is already known. The Sage-Husa algorithm with adaptive estimation of R when Q is known is outlined in SHEKF4 [35]. The final algorithm is known as the EKF algorithm [36], therefore it does not include the adaptive estimation of noise. The end voltage value from the pulse discharge experiment serves as the voltage reference value. Additionally, the reference value for algorithm verification is the SOC value computed using the ampere-hour (Ah) integration method during a lengthy resting period. As evaluation metrics, the root means square error (RMSE) and the mean absolute error (MAE) are utilized to quantify the estimation mistakes produced by each of the five algorithms. The following definitions of MAE and RMSE are presented.
section_name: Experiments and Analysis

section_content: The literature frequently makes use of the second-order RC model. However, the thirdorder model is not widely used since, despite increasing accuracy, the convergence time is proportionately longer. In this paper, the enhanced methods-a streamlined algorithm coupled with a third-order model are employed to resolve the issue. 

The model validation approach of this paper is separated into two parts. In the preliminary stage, the model's correctness is examined. By incorporating the specified model parameters into the end voltage solution formula, the simulated value of the end voltage is achieved. The size of the discrepancy between the measured and simulated values can be used to gauge the precision of model parameter identification. These two models are contrasted in this paper for the pulsed discharge and DST conditions, respectively. 

To start, the voltage measurement is obtained using pulse discharge experiments. Figure 5 displays the outcomes of the simulated values and the measured values of the cells for both models after substituting the parameter identification findings of the models. The highest inaccuracy of the third-order model and the second-order model can both be observed in the error plot to be less than 3. 5%. It is shown that both models accurately represent the real battery situation. After that, an experiment using DST working conditions was used to get the voltage data. The simulated values and measured values of the battery are given in Figure 6 after the parameter identification findings of the two models were again swapped. The greatest error of the models is less than 3%, as can be seen from the error plot. Additionally, we can deduce from Figure 7 that the third-order model's mean error is significantly lower than the second-order model's. Furthermore, we can infer from Figure 7 that the third-order model's mean error is roughly 0. 33%. The second-order model's mean inaccuracy is roughly 1%. 

The convergence time of model is then verified as the second step. To guarantee the validity of the measurements, the experiment was conducted five times, and the findings were averaged. The results are displayed in Table 2 below. According to Table 2, the third-order model's convergence time is 0. 0482 s longer than that of the second-order under the pulse discharge condition, whereas it is 0. 154 s shorter under the DST condition. Without taking into account the impact of computational mistakes, it can be said that the simplification approach can be used to get over the higher-order model's issue with a slow convergence rate. In conclusion, the third-order model has a higher level of model accuracy than the second-order model. When the battery output current value fluctuates frequently, the thirdorder RC model performs better and is more stable. Figure 7 provides a clearer visual example of this. In other words, when the operating conditions are more complicated, the third-order model is preferable to the second-order model.
section_name: Model Verification

section_content: To verify the accuracy and consistency of the proposed algorithms that was presented in this paper, which is the primary objective of this section. The values corresponding to the true values are assigned to the initial values of SOC. Under the same initial circumstances, the five algorithms are contrasted. Figure 8 displays the outcomes of SOC estimate utilizing the different algorithms as well as any errors that may have occurred. Figure 8 demonstrates that the improved algorithm estimates SOC more accurately than the other algorithms. There is a greater estimation error in the first stage because SHEKF1 and SHEKF2 integrate adaptive estimation of the measurement noise. Adaptive estimate on the measurement noise R is exceedingly difficult because of its extreme instability. Because of this, many academics contend that adaptive estimation of Q and R cannot be done at the same time. Figure 8 shows that when the initial SOC value is accurate, the revised algorithm is noticeably more stable. And the algorithm proposed in this paper can converge to the true value when the sampling time is about 870 s. When the SOC value at startup is accurate, the MAE and RMSE for each algorithm are listed in Table 3. According to the table, the revised SHEKF1 algorithm reduces the mean absolute error by 0. 23% when compared to EKF. The SHEKF3 algorithm, which has been used extensively in publications, has a forecast error that is 0. 05% less than that of EKF. The comparison of SHEKF1 and SHEKF3 shows that the improved algorithm proposed outperforms the commonly used adaptive algorithm for estimating EV battery state of charge SOC. SOC estimation accuracy is lower for the SHEKF2 algorithm, which estimates both measurement noise and system noise, and the SHEKF4 algorithm, which only estimates measurement noise. This also demonstrates that the method of simultaneously simplifying the noise covariance matrix and estimating only R will sacrifice some accuracy, which is not applicable in estimating the SOC of lithium batteries.
section_name: Evaluation of Correct Initial SOC Value

section_content: Incorrect initial values will directly impact the estimation precision of the algorithm. During the operation of an electric vehicle, the preliminary SOC worth is not constant. Consequently, the estimation algorithm must be robust to the initial SOC value. Given that the electric vehicle charge is between 50% and 20%, it is even more important to ensure that the battery charge is correct. As a result, this paper uses the incorrect initial SOC values of 0. 5 and 0. 2 to test the algorithm's effectiveness. Figure 9a, b show the simulation results when SOC is set to 0. 5 as the starting value. Furthermore, the simulation results are shown in Figure 10a,b, and SOC is set to 0. 2 at the start. Figures 9a and 10a illustrate that the presented algorithm that this paper proposes is more resistant to the initial setting of the SOC. The improved algorithm can converge to the true value after approximately 400 s of sampling. Even in low SOC time periods, the improved algorithm can improve estimation. Furthermore, as shown in Figures 9b and 10b, the improved algorithm's maximum error does not exceed 0. 4%. The maximum error between the two algorithms SHEKF2 and SHEKF4 is greater than 3%. 

To put it another way, the numerical results obtained by using each of the five different algorithms are depicted in Figure 11. Figure 11 clearly shows that regardless of the initial SOC value, the improved algorithm proposed in this paper outperforms other algorithms. It can be deduced that a larger difference between the initial and accurate SOC values has a greater impact on SOC estimation. Although the EKF algorithm, which is directly based on the SHEKF algorithm, removes adaptive noise estimation, the effect is superior to the SHEKF2 and SHEKF4 algorithms. In the case of different initial SOC values, the EKF algorithm's estimation error varies by a maximum of 0. 13%. This demonstrates that the inaccurate SOC initial value is equivalent to a portion of the noise, and the EKF algorithm cannot perform adaptive noise estimation. To improve the accurateness of the SOC, an algorithm that is not robust to the initial value of the SOC should permit the initial value to be set as close to the actual value as possible.
section_name: Evaluation of Incorrect Initial SOC Value

section_content: The inclusion of a forgetting factor distinguishes the Sage-Husa algorithm from other adaptive algorithms. This adaptive algorithm's basic goal is to emphasize the importance of recent data while gradually forgetting the importance of earlier data. However, the forgetting factor must be chosen in accordance with the system's actual working environment. Relying on the fast and slow change degrees of these two types of noise, this paper's forgetting factor for system noise and measurement noise is determined. The reasons for this are discussed further in the section III. Some simulation experiments show that the algorithm's accuracy and stability increase as the difference between the forgetting factor of the system noise and the initial value of the measured noise increases. Furthermore, in this paper, the optimal forgetting factor initial value difference was 0. 045 (i. e., b1 = 0. 995, b2 = 0. 95). 

Figure 12 shows that as the difference between the initial values of the two forgetting factors decreases, so do the algorithm's maximum error and stability. The RMSE and MAE of the difference between the various forgetting factors of the modified algorithm are summarized in Table 4. The data in Table 4 's first and second columns can be used to demonstrate that, when it exceeds 0. 045, the algorithm is not significantly affected by the difference between the two forgetting factors' initial values. Furthermore, the simulation results show that when the difference between the initial values of the two forgetting factors is 0. 03, the estimation error is 0. 5% higher than when the difference is 0. 045. It is advantageous to improve the algorithm's accuracy by treating the forgetting factors of the two types of noise differently.
section_name: Impact of the Forgetting Factor on SOC Estimation

section_content: The capacity of Li-ion batteries can be impacted by temperature. As a result, one of the determining factors in estimating SOC is temperature. This section simulates the operation of a Li-ion battery power system using the US dynamic stress test (DST). To test the effectiveness of the enhanced algorithm, three different kinds of temperature data representing low temperature (0 °C), normal temperature (25 °C) and high temperature (45 °C) were also selected. The precise current and amplification curves for the DST operating conditions are shown in Figure 13. Each current is made up of many tiny cycles that operate over a very long dynamic period. As a result, the DST operating environment can represent the battery usage more accurately. In this discharge case, five algorithms are used to calculate the SOC, with an initial SOC value of 0. 8 for each algorithm. The estimated SOC curve and the estimation error at 25 °C are shown in Figure 14. The SHEKF1 algorithm has the smallest MAE, which is about 2% smaller than the error of the EKF algorithm, according to a comparison of the data in each column of Table 5 in the first place. 14b, 15b and 16b show a clear correlation between temperature change and an increase in the SOC estimation algorithm's error. To make it easier to understand the benefits and drawbacks of each algorithm, the error data for the three temperatures are condensed in Tables 5 6 7. Similar to the 25 °C pulse discharge condition, the MAE of the SHEKF3 algorithm is almost 0. 8% less than that of the SHEKF1 and SHEKF4 algorithms. Table 6 's comparison of the data in each column still demonstrates the SHEKF1 algorithm's superiority. The SHEKF1 algorithm has a 1. 4% smaller error than the SHEKF1 algorithm, which has the largest error. five algorithms for SOC estimation at the aforementioned three temperatures. In addition, the SHEKF1 algorithm has roughly 1. 3% less error than the EKF algorithm, as shown by a comparison of the data in each column of Table 7.
section_name: The Algorithm Performance under DST Working Condition

section_content: Analyzing the aforementioned simulation experimental data can lead to the following conclusions:

i. SOC estimation algorithms based on models, such as the Sage-Husa and EKF algorithms, are highly model-dependent. If the model is not accurate enough, none of these algorithms can calculate SOC accurately enough. The accuracy of SOC estimation may be jeopardized if the model fails. As a result, it is critical to validate the developed models before performing SOC estimation. Furthermore, the superiority of the third-order model over the second-order model is not readily apparent under operating conditions with reasonably consistent current variation. However, when dealing with complex current variations, the third-order model outperforms the second-order model significantly. The reason for this is that sudden changes in current can exacerbate cell polarization. The former increases the number of RC links to more accurately describe the actual state of the battery polarization phenomenon. In terms of convergence time, this paper believes that combining the simplified algorithm with the third-order RC model can reduce the algorithm's convergence time. The convergence time of the simplified and improved algorithm SHEKF1 combined with the third-order RC model is 0. 15 s faster under DST conditions than the second-order RC model. This suggests that combining the higher-order model with the simplified algorithm can result in high accuracy without increasing the algorithm's convergence time. 

ii. The streamlined Sage-Husa algorithm continues to be utilized comparatively more frequently than the original Sage-Husa algorithm. Eliminating the factors that cause the noise covariance matrix to become non-negative definite can improve the filtering's stability, albeit at the expense of some precision. Additionally, the technique of directly removing the subtractive part that could result in filter divergence is incredibly effective despite being very basic and crude. As a result, it is essential to reduce the complexity of the covariance matrix created by the two noises and do away with the estimation of the expectation created by the two noises. The experiments in this paper show that the simplified SHEKF algorithm is very effective. 

iii. SHEKF1 and SHEKF2 comparison shows that the adjustment of the two forgetting factors can be more flexible. When compared to the Sage-Husa algorithm with one forgetting factor, the Sage-Husa algorithm with two forgetting factors significantly improves SOC estimation accuracy. Furthermore, the EKF algorithm is transformed into the EKF algorithm by removing the adaptive noise estimation based on the simplified SHEKF algorithm. This EKF algorithm outperforms the SHEKF2 and SHEKF4 algorithms in terms of stability. Many papers propose the SHEKF3 algorithm, which simplifies the SHEKF algorithm and adaptively estimates only the system noise and has an error of no more than 1. 5% less than the EKF. The improved SHEKF1 algorithm has a maximum error of 2. 2% less than the EKF and is resistant to the initial value of the SOC. To summarize, future improvements to the SHEKF algorithm should focus on different treatment of system noise and measurement noise in order to reduce computational degree and achieve higher precision. 

iv. As demonstrated in Section 4 of this paper, the magnitude of the selected forgetting factor value correlates directly with the precision of SOC estimation. However, there is no accepted formula to measure whether the set forgetting factor is reasonable or not. Therefore, the best results were obtained by setting the forgetting factor to 0. 95 for the Q and 0. 995 for the R through several simulations. In addition, the SOC estimation is thought to be more accurate when there is a greater disparity between the two. In the actual SOC estimation, it is necessary to consider both stability and real-time performance. Since the causes of system noise and measurement noise are distinct, two optimization algorithms, asymptotic memory index weighting and limited memory index weighting, are selected to make adaptive adjustments to the noise. Numerous experimental data demonstrate that this processing yields superior SOC estimation results.
section_name: Discussion

section_content: An estimate of SOC is obtained in this paper using a model-based methodology. The model part uses a third-order RC equivalent circuit instead of the second-order model that is typically used in literature. For the algorithm, two forgetting factors are used in a simplified Sage-Husa extended Kalman filter (SHEKF) algorithm. Using the limited memory index weighting method and the asymptotic memory index weighting method, the forgetting factors are chosen so that they adapt to the different levels of fast and slow changes in the measurement noise and the statistical properties of the system noise. The SOC estimation accuracy has increased by 0. 2% when contrasted with the EKF algorithm. The enhanced algorithm for estimating SOC outperforms other SHEKF algorithms in terms of robustness and accuracy. The improved algorithm for estimating SOC outperforms other SHEKF algorithms in terms of accuracy and persistence. Furthermore, by combining the third-order model with this simplified method, the problem of the higher-order model's long convergence time can be solved. The improved algorithm converges 0. 154 s quicker than the second-order model for the DST condition.
section_name: Conclusions

section_content: Funding: This research was funded by the Natural Science Foundation of Chongqing, China (No. cstc2021jcyj-msxmX0301, No. CSTB2022NSCQ-MSX1675, and No. 2022NSCQ-MSX3846 ), Science and Technology Project of Chongqing Municipal Education Commission ( KJZD-K202101202, KJZD-K201901203 ), Science and Technology Project of Wanzhou District, Chongqing, China and Postgraduate Research and Innovation Project of Chongqing Three Gorges University ( 072206 ). 

Institutional Review Board Statement: Not applicable for studies not involving humans or animals. 

Informed Consent Statement: Not applicable for studies not involving humans.

section_name: Conflicts of Interest:

section_content: Data Availability Statement: Not applicable.

section_content: The authors declare no conflict of interest.
section_name: Conflicts of Interest:

section_name: Abbreviations

section_content: The information in Tables 5 and 6 is compared for analysis. The temperature drops from 25 °C to 0 °C, while the RMSE and MAE of the SHEKF1 algorithm both rise by about 1. 5% and 0. 6%, respectively. The MAE of the SHEKF1 algorithm increased by about 0. 4% when comparing the data in Table 6 and Table 7. The superiority of the improved algorithm, which is the SHEKF1 algorithm in the calculation example, is identified after comparing the

section_content: The following abbreviations are used in this manuscript: BMS
section_name: Abbreviations