State of Charge Estimation of Lithium Battery Based on Improved Correntropy Extended Kalman Filter

: State of charge (SOC) estimation plays a crucial role in battery management systems. Among all the existing SOC estimation approaches, the model-driven extended Kalman ﬁlter (EKF) has been widely utilized to estimate SOC due to its simple implementation and nonlinear property. However, the traditional EKF derived from the mean square error (MSE) loss is sensitive to non-Gaussian noise which especially exists in practice, thus the SOC estimation based on the traditional EKF may result in undesirable performance. Hence, a novel robust EKF method with correntropy loss is employed to perform SOC estimation to improve the accuracy under non-Gaussian environments ﬁrstly. Secondly, a novel robust EKF, called C-WLS-EKF, is developed by combining the advantages of correntropy and weighted least squares (WLS) to improve the digital stability of the correntropy EKF (C-EKF). In addition, the convergence of the proposed algorithm is veriﬁed by the Cram é r–Rao low bound. Finally, a C-WLS-EKF method based on an equivalent circuit model is designed to perform SOC estimation. The experiment results clarify that the SOC estimation error in terms of the MSE via the proposed C-WLS-EKF method can e ﬃ ciently be reduced from 1.361% to 0.512% under non-Gaussian noise conditions.


Introduction
In recent years, electric vehicles (EV) have become a trend in the automotive industry due to their advantages of no emissions, low energy consumption and low noise [1]. The battery management system (BMS) is the core device of EV, in which the state of charge (SOC) is an important parameter reflecting the state of the battery residual capacity. Hence, to ensure the scientific and effective operation of the battery, designing an effective SOC estimation method is an essential part of the BMS [2][3][4]. Most of the functions in the system depend on the results of an SOC evaluation of the power battery. Therefore, an accurate SOC estimation of the battery is beneficial to protect the battery, prevent battery overcharge or over discharge, improve battery life and achieve the purpose of energy saving.
At present, SOC estimation has attracted widespread attention and many effective SOC estimation methods have been developed by researchers. Generally, the SOC estimation techniques can be classified into three major categories [5]. The first is to establish the relationship between voltage, current and SOC by directly measuring the voltage and current of the battery [6,7]. The second is a data-driven estimation approach such as neural networks, and both Tong [8] and Yang [9] used a neural network or an improved neural network for lithium battery SOC estimation. The third can be The rest of the paper is organized as follows. Section 2 briefly introduces the equivalent circuit model (ECM) and model parameter identification. Section 3 introduces the C-EKF and C-WLS-EKF algorithms proposed. The performance of the proposed C-WLS-EKF model is verified via experimental data in Section 4. Finally, the conclusion is presented in Section 5.

Equivalent Circuit Model and Parameter Identification
The convergence and accuracy of the battery SOC estimation are largely reliant on the exact identification of the battery model. The mainstream equivalent circuit models (ECMs) include the Rint model, R-C model, PNGV model and Thevenin model [28]. Considering the relationship between estimation accuracy, model complexity and battery dynamic characteristics, the secondorder R-C network model was selected to represent the ECM used in this work, as shown in Figure  1.
where   k  represents the SOC estimation value, N Q is battery capacity, and t  is the sampling interval.

Relationship between SOC and OCV
The accuracy of SOC estimation is affected by many factors such as temperature [29], cells balancing [30,31] and OCV value [27]. In this work, considering the complexity of the model and the limitations of the experimental condition, only SOC is regarded as a non-linear function of OCV. To obtain this nonlinear function, the OCV test was conducted using a 1.6 Ah lithium battery as a case study. The test temperature was maintained at the nominal temperature (25 °C), then the lithium battery was discharged with a 0.2 C constant current. When the battery discharged 10% of its rated capacity, we stopped the discharge process and let it stand for 30 min, then we cycled the process until the discharge was completed. This is recorded in Table 1. where V c represents the terminal voltage of the battery, V oc is the open circuit voltage (OCV), I c denotes the battery's load current during charging and discharging, R 0 represents the Ohm internal resistance, the capacitive resistance branch composed of R s and C s represents the concentration polarization of the battery, whereas the other branch composed of R l and C l denotes electrochemical polarization, and V s and V l are the short-time and long-time transient voltage responses of the battery, respectively.
According to the second-order R-C network model, the battery's discrete state space equation can be established as where ζ(k) represents the SOC estimation value, Q N is battery capacity, and ∆t is the sampling interval.

Relationship between SOC and OCV
The accuracy of SOC estimation is affected by many factors such as temperature [29], cells balancing [30,31] and OCV value [27]. In this work, considering the complexity of the model and the limitations of the experimental condition, only SOC is regarded as a non-linear function of OCV. To obtain this nonlinear function, the OCV test was conducted using a 1.6 Ah lithium battery as a case study. The test temperature was maintained at the nominal temperature (25 • C), then the lithium battery was discharged with a 0.2 C constant current. When the battery discharged 10% of its rated capacity, we stopped the discharge process and let it stand for 30 min, then we cycled the process until the discharge was completed. This is recorded in Table 1. In order to more accurately establish the SOC-OCV relationship, in this experiment, the data in Table 1 are fitted in the form of a sixth polynomial plus logarithm by using the lsqcurvefit function in the MATLAB software, and compared with the eighth-order polynomial fitting implementation used in [32]. The fitting results are shown in Figure 2. Meanwhile, the standard deviation (STD) and goodness of fit (R-square) are utilized to analyze the fitting results, as shown in Table 2. The corresponding results in Table 2 and Figure 2 show that the fitting effect of the sixth-order polynomial plus logarithmic term applied in this paper is better than that of the eighth-order polynomial fitting. Hence, the nonlinear expression of SOC-OCV can be expressed as Energies 2020, 13, x FOR PEER REVIEW 4 of 18 In order to more accurately establish the SOC-OCV relationship, in this experiment, the data in Table 1 are fitted in the form of a sixth polynomial plus logarithm by using the lsqcurvefit function in the MATLAB software, and compared with the eighth-order polynomial fitting implementation used in [32]. The fitting results are shown in Figure 2. Meanwhile, the standard deviation (STD) and goodness of fit (R-square) are utilized to analyze the fitting results, as shown in Table 2. The corresponding results in Table 2 and Figure 2 show that the fitting effect of the sixth-order polynomial plus logarithmic term applied in this paper is better than that of the eighth-order polynomial fitting. Hence, the nonlinear expression of SOC-OCV can be expressed as  

Model Parameter Identification
It is necessary to perform parameter identification of the relevant variables R o , R s , C s , R l , C l in the ECM of the battery when the SOC estimation is conducted. Here, the lithium battery pulse discharge test was implemented to study the battery's dynamic impedance characteristics under different pulse current conditions. Figure 3 presents the battery's voltage response at SOC = 1. Through the excitation response analysis of the experimental data, the least squares (LS) method is utilized to identify the parameters of the ECM under different SOCs.
the ECM of the battery when the SOC estimation is conducted. Here, the lithium battery pulse discharge test was implemented to study the battery's dynamic impedance characteristics under different pulse current conditions. Figure 3 presents the battery's voltage response at SOC = 1. Through the excitation response analysis of the experimental data, the least squares (LS) method is utilized to identify the parameters of the ECM under different SOCs.
According to the experimental data, the voltage changes under pulse currents can be divided into four stages: discharge ( 01 UU ), discharge stop ( During the charge or discharge stop phase, due to the action of s C and l C , the change of the circuit voltage can be regarded as a zero-input response process, and the loop voltage can be expressed as   Table 3.

Parameter Identification Result Verification
As aforementioned, the accuracy of SOC estimation largely depends on the accuracy of the initial parameter identification. Therefore, it is necessary to perform rationality verification on the identification results. The model output can be obtained as According to the experimental data, the voltage changes under pulse currents can be divided into four stages: discharge (U 0 ∼ U 1 ), discharge stop (U 1 ∼ U 2 ), charge (U 2 ∼ U 3 ) and charge stop (U 3 ∼ U 4 ). According to Ohm's law of battery Ohmic resistance R 0 = ∆U/I 0 , the battery charge resistance can be computed as During the charge or discharge stop phase, due to the action of C s and C l , the change of the circuit voltage can be regarded as a zero-input response process, and the loop voltage can be expressed as The LS fitting technique is applied to identify R o , R s , C s , R l , C l , the results are shown in Table 3.

Parameter Identification Result Verification
As aforementioned, the accuracy of SOC estimation largely depends on the accuracy of the initial parameter identification. Therefore, it is necessary to perform rationality verification on the identification results. The model output can be obtained as Energies 2020, 13, 4197 6 of 18 Putting the initial parameters obtained by identification into (5), the result is shown in Figure 4, in which the red line represents the measured voltage, and the blue represents the model output value. It can be seen from Figure 4 that the blue line basically coincides with the red line, which also means that the parameters obtained above are valid.
Energies 2020, 13, x FOR PEER REVIEW 6 of 18 Putting the initial parameters obtained by identification into (5), the result is shown in Figure 4, in which the red line represents the measured voltage, and the blue represents the model output value. It can be seen from Figure 4 that the blue line basically coincides with the red line, which also means that the parameters obtained above are valid.

SOC Estimation Based on EKF with Correntropy
In this section, we focus on the development of C-EKF and C-WLS-EKF. Considering the battery's nonlinear characteristics, a linear state equation and a nonlinear measurement equation are used to represent the battery state.
where k x and Although the EKF has been extensively used for SOC estimation due to its ease of implementation [33][34][35], it still has some flaws such as the filter instability due to the Jacobian matrix [36]. Meanwhile, as pointed out in [37], the EKF is optimal under Gaussian noise conditions. However, in the case of non-Gaussian noise, the model only considers second-order moment information, so its filter is a suboptimal estimate. Therefore, a robust SOC estimation scheme via the C-EKF based on MCC with higher-order information will be developed to achieve better stability and accuracy, and we first review the correntropy in the following subsection.

Correntropy
Correntropy, as a measure of local similarity that is insensitive to outliers and noise, has been successfully applied in fields such as pattern recognition [38], machine learning [39] and signal processing [40]. Correntropy [22,[41][42][43] can be expressed as follows

SOC Estimation Based on EKF with Correntropy
In this section, we focus on the development of C-EKF and C-WLS-EKF. Considering the battery's nonlinear characteristics, a linear state equation and a nonlinear measurement equation are used to represent the battery state. x where x k and y k are the state variable and output variable, and w k ∼ N(0, Q k ) and v k ∼ N(0, R k ) denote the process noise and measurement noise; where Q k represents the covariance of w k , R k represents the covariance of v k , F and B are the state transition coefficient matrix and the gain matrix, respectively, u k is the model input matrix, and g denotes the measurement model. The main difference between EKF and KF is that EKF uses a Jacobian matrix to transform nonlinear equations into linear equations, which can be obtained as Therefore, the new measurement equation can be written as Although the EKF has been extensively used for SOC estimation due to its ease of implementation [33][34][35], it still has some flaws such as the filter instability due to the Jacobian matrix [36]. Meanwhile, as pointed out in [37], the EKF is optimal under Gaussian noise conditions. However, in the case of non-Gaussian noise, the model only considers second-order moment information, so its filter is a suboptimal estimate. Therefore, a robust SOC estimation scheme via the C-EKF based on MCC with higher-order information will be developed to achieve better stability and accuracy, and we first review the correntropy in the following subsection.

Correntropy
Correntropy, as a measure of local similarity that is insensitive to outliers and noise, has been successfully applied in fields such as pattern recognition [38], machine learning [39] and signal processing [40]. Correntropy [22,[41][42][43] can be expressed as follows where X and Y represent random variables, E[·] is the expectation operator, k σ (·, ·) denotes a kernel function with kernel width σ, and f XY (·, ·) is the joint density function of X and Y. Generally, f XY (·, ·) is not obtained in practice and the available samples are limited [22,44]. Hence, the sample estimator of the correntropy is recorded asV The Gaussian kernel function is utilized as the kernel function of correntropy in this paper, Hence, (11) can be rewritten asV We can observe from (12) that the correntropy takes the maximum value when X = Y, which is similar with the minimum error in MSE. Therefore, this paper uses correntropy loss instead of MSE loss to develop a more robust estimation method.

EKF with Correntropy
In this section, the C-EKF filter is mainly introduced by using the concept of correntropy as the loss function in traditional EKF.
Combining (9) and (12), a new loss function is defined as By calculating the partial derivative of the objective function J m with respect to the state variablê x k and setting it equal to zero, we obtain Then, we can from (14) get the state update equation aŝ In the above formula, this paper assumes thatx k ≈x k , so then (6), and we bring it into the denominator part of (15):x where In the follow-up work of this section, this paper will continue to derive on the basis of (15)- (17). Under the assumptionx k ≈x k , (15) can be represented aŝ Energies 2020, 13, 4197 8 of 18 Then, to further introduce the measurement error term, we correct the right side of (18) by factor G σ ( y k −Hx k−1 ) Now, we have Both sides of (20) are divided by (I + G σ ( y k −Hx K ) G σ ( Bu k ) H T H), so we get where To this end, the C-EKF can be summarized in Table 4. Table 4. Correntropy extended Kalman filter (C-EKF).

C-EKF with WLS
As with the traditional EKF, the noise covariance has an impact on the performance of the C-EKF. To further improve the efficiency of the C-EKF, this paper will combine WLS with correntropy to design a more robust filter. Correntropy is used to extract high-order statistical information of an SOC estimation and introduce covariances R −1 k and P −1 k|k−1 through WLS to minimize the variance estimation. Therefore, the new loss function can be expressed as Computing the partial derivative of the J m with respect to the state variablex k and setting it to zero, we have the following: Energies 2020, 13, 4197 9 of 18 Then, after some simplification of (24), we obtain where Adding and subtracting L k H T R −1 k Hx k on the right side of (25) yields Then, the iterative equation can be written aŝ where In the new algorithm, the noise covariance existing in Equations (27) and (29) will affect the gain K k of the system. Finally, the filter, called C-WLS-EKF, uses the correntropy to process the non-Gaussian properties of the data and enhance the digital stability by using WLS, which is summarized in Table 5. Table 5. C-WLS-EKF.

Convergence Analysis of the C-WLS-EKF Algorithm
The Cramér-Rao low bound criterion (CRLB) is usually utilized to verify the convergence of a novel adaptive filtering algorithm [45]. To show the effectiveness of the proposed algorithm theoretically, the CRLB of the proposed C-WLS-EKF algorithm will be derived under nonlinear conditions. Here, we define Y n n 0 +1 = y n 0 +1 , y n 0 +2 , . . . , y n where Y n n 0 +1 represents the output of the measurement equation. The CRLB can be used to evaluate the performance of the filter and ensure the estimation accuracy of the state variable x k . Y n n 0 = H T x n 0 , y n 0 +1 , y n 0 +2 , . . . , y n = y n 0 , y n 0 +1 , y n 0 +2 , . . . , y n According to Van Trees, the amount of Fisher information in an unknown state x k can be expressed by where P n n 0 is the conditional function of Y n n 0 for X n n 0 . P n n 0 = P Y n n 0 X n n 0 X n n 0 = x n 0 , x n 0 +1 , x n 0 +2 , . . . , x n For non-Gaussian noise, P n n 0 can be expressed as [46] P n n 0 = where R n 0 , R n 1 represent the variance of Gaussian noise and short noise, respectively. Taking the logarithm of (35): Then, bringing ln P n n 0 into (33), we can obtain the following This paper addresses the problem that the oscillating system is not suitable for evaluation by using the average of Fisher's information as If the process is ergodic, we usually think that the phase space average is equal to the time average, then According to the definitions of global Lyapunov functions, (39) can be expressed as When the system has chaotic characteristics, the Lyapunov functions define its parameter λ as a positive number. From (40), we have Energies 2020, 13, 4197 11 of 18 According to the result of (41), one can know that the function (40) is convergence, that is, the derived C-WLS-EKF algorithm will converge to an optimal.

Experimental Results
In this section, four sets of simulations are implemented to evaluate the SOC estimation performance of the C-EKF and C-WLS-EKF methods under non-Gaussian noise conditions. The schematic diagram of the lithium battery discharge test system is given in Figure 5. It mainly includes three parts, the first part is the control system which plays the information interaction role, the second is the lithium battery and Thermostat (control temperature and humidity) and the third is the hardware part of the lithium battery charge and discharge test system, which is mainly composed of a signal acquisition unit, regulated power supply, charge and discharge instrument and current transformer. In the following experiments, the original experiment data set is obtained from the lithium battery test system above. Meanwhile, all simulations are implemented on a desktop computer configured with an Intel(R) Core (TM)i5-7400 CPU. To show the robustness of the proposed method, the elements of process noise vector w k , along with the elements of the measurement noise v k , are comprised of Gaussian noise and short noise. Figure 6 shows the non-Gaussian noise applied to the experiment.     Then, the non-Gaussian mixture noise can be expressed by w k = N(µ x , Q) + short noise (42) v k = N(µ y , R) + short noise (43) From Section 3 above, one can observe that the EKF based on correntropy loss will introduce a kernel width σ, the value of which will affect the SOC estimation accuracy to a certain extent. In order to effectively improve the accuracy of SOC estimation, the appropriate kernel width σ should be selected first, and hence the influence of σ for the proposed method is simulated and analyzed.
Firstly, we conduct experiments to investigate the effect of the accuracy of the σ proposed for the C-WLS-EKF method to further obtain an optimal kernel width for the remaining experiments. Here, the initial SOC value is set at 1, and all the parameters except σ remain unchanged. The SOC estimation results under different kernel widths are shown in Figure 7, where one can clearly observe through the partial enlarged view that the SOC estimation accuracy changes with the change of the kernel width σ. In addition, to more intuitively understand the change trend, we conduct 10 independent trails to obtain the average of the estimate errors as shown in Table 6. One can know from the corresponding results that the estimated performance is the best when σ = 0.5. Moreover, the kernel width σ not only affects the estimation accuracy, but also affects the convergence performance of the algorithm. Figure 8 gives the convergence curves of the algorithm when the kernel width σ is set at different values. It can be seen from Figure 8 that there are significant differences in the convergence speed of the algorithm when the kernel width is set at different values. When the kernel function σ is taken as 0.5, the convergence rate of the algorithm is the fastest and the steady-state performance is the highest. Hence, all the kernel widths are taken as 0.5 in this study.
Energies 2020, 13, x FOR PEER REVIEW 13 of 18 width  is set at different values. It can be seen from Figure 8 that there are significant differences in the convergence speed of the algorithm when the kernel width is set at different values. When the kernel function  is taken as 0.5, the convergence rate of the algorithm is the fastest and the steadystate performance is the highest. Hence, all the kernel widths are taken as 0.5 in this study.  Subsequently, the estimation performance comparison experiment between the proposed method and the traditional EKF method was performed. The SOC estimation results and the error curves are shown in Figures 9 and 10, respectively, in which the red line is the experimental data and other lines represent the simulation results. The vertical axis represents the SOC values from 0 to width  is set at different values. It can be seen from Figure 8 that there are significant differences in the convergence speed of the algorithm when the kernel width is set at different values. When the kernel function  is taken as 0.5, the convergence rate of the algorithm is the fastest and the steadystate performance is the highest. Hence, all the kernel widths are taken as 0.5 in this study.  Subsequently, the estimation performance comparison experiment between the proposed method and the traditional EKF method was performed. The SOC estimation results and the error curves are shown in Figures 9 and 10, respectively, in which the red line is the experimental data and other lines represent the simulation results. The vertical axis represents the SOC values from 0 to  Subsequently, the estimation performance comparison experiment between the proposed method and the traditional EKF method was performed. The SOC estimation results and the error curves are shown in Figures 9 and 10, respectively, in which the red line is the experimental data and other lines represent the simulation results. The vertical axis represents the SOC values from 0 to 100%. The horizontal axis represents the simulation time. The battery was discharged by withdrawing a total capacity of 1 C-rate at constant current. One can observe from the simulation results in Figures 9 and 10 that (1) the estimation accuracy is low in the first 500 s because both the EKF and the improved EKF algorithm require a certain convergence time. However, the improved C-EKF and C-WLS-EKF algorithm have a faster convergence speed; (2) after 500 s, all aforementioned algorithms can effectively estimate the SOC; (3) the improved C-WLS-EKF algorithm achieves better estimation accuracy than the traditional EKF and improved C-EKF algorithms under non-Gaussian noise conditions. In addition, to quantitatively compare the performance of each algorithm, the estimation errors and CPU run time of the three algorithms are given in Table 7. According to this result in Table 7, one can know that (1) the estimated error of C-WLS-EKF achieves 0.512%, which is lower than the estimated error of the C-EKF (0.771%) and the estimated error of the EKF (1.361%). Therefore, one can be concluded that the proposed C-WLS-EKF algorithm is effective for SOC estimation under non-Gaussian noise; (2) the proposed C-EKF and C-WLS-EKF algorithms only took a little more run time from the MCC optimization in comparison to the traditional EKF, while the higher estimation precision in this case can be obtained.
Next, this paper investigates the influence of the SOC initial value on the tracking performance of the proposed algorithm. Keeping other parameters unchanged, the initial value range of SOC is set to [0. 2,1], and the experiment result is shown in Figure 11. It can be seen from Figure 11 that as the initialization value decreases, the proposed algorithm requires more iteration time to track the actual state of the battery, which also means that the convergence speed of the algorithm decreases and the estimation error increases. At the same time, it can be observed that the proposed algorithm can achieve an effective estimation after a certain time of convergence.
Energies 2020, 13, x FOR PEER REVIEW 14 of 18 estimation errors and CPU run time of the three algorithms are given in Table 7. According to this result in Table 7, one can know that (1) the estimated error of C-WLS-EKF achieves 0.512%, which is lower than the estimated error of the C-EKF (0.771%) and the estimated error of the EKF (1.361%). Therefore, one can be concluded that the proposed C-WLS-EKF algorithm is effective for SOC estimation under non-Gaussian noise; (2) the proposed C-EKF and C-WLS-EKF algorithms only took a little more run time from the MCC optimization in comparison to the traditional EKF, while the higher estimation precision in this case can be obtained.  Next, this paper investigates the influence of the SOC initial value on the tracking performance of the proposed algorithm. Keeping other parameters unchanged, the initial value range of SOC is set to [0.2,1], and the experiment result is shown in Figure 11. It can be seen from Figure 11 that as the initialization value decreases, the proposed algorithm requires more iteration time to track the actual state of the battery, which also means that the convergence speed of the algorithm decreases and the estimation error increases. At the same time, it can be observed that the proposed algorithm can achieve an effective estimation after a certain time of convergence.   Figure 11. Estimated results under different initial SOCs.
Finally, in order to evaluate the performance of the proposed C-WLS-EKF algorithm, we further conduct a comparative experiment with the data-driven methods (least squares support vector machine (LSSVM) and BP neural network). One can see from Figure 12 that the improved C-WLS-EKF algorithm in this paper has much higher accuracy than the intelligent algorithm. After calculation, the estimation errors of SOC using LSSVM and the BP neural network are 1.84% and 2.49%, respectively, which are higher than the proposed algorithms in this paper. At the same time, for a more intuitive analysis, the SOC estimation error is presented in Figure 13. Moreover, it is worth mentioning that the time complexity of the data-driven methods such as LSSVM and the BP neural network is higher than the proposed C-WLS-EKF-based SOC estimation method. Finally, in order to evaluate the performance of the proposed C-WLS-EKF algorithm, we further conduct a comparative experiment with the data-driven methods (least squares support vector machine (LSSVM) and BP neural network). One can see from Figure 12 that the improved C-WLS-EKF algorithm in this paper has much higher accuracy than the intelligent algorithm. After calculation, the estimation errors of SOC using LSSVM and the BP neural network are 1.84% and 2.49%, respectively, which are higher than the proposed algorithms in this paper. At the same time, for a more intuitive analysis, the SOC estimation error is presented in Figure 13. Moreover, it is worth mentioning that the time complexity of the data-driven methods such as LSSVM and the BP neural network is higher than the proposed C-WLS-EKF-based SOC estimation method.
EKF algorithm in this paper has much higher accuracy than the intelligent algorithm. After calculation, the estimation errors of SOC using LSSVM and the BP neural network are 1.84% and 2.49%, respectively, which are higher than the proposed algorithms in this paper. At the same time, for a more intuitive analysis, the SOC estimation error is presented in Figure 13. Moreover, it is worth mentioning that the time complexity of the data-driven methods such as LSSVM and the BP neural network is higher than the proposed C-WLS-EKF-based SOC estimation method.  calculation, the estimation errors of SOC using LSSVM and the BP neural network are 1.84% and 2.49%, respectively, which are higher than the proposed algorithms in this paper. At the same time, for a more intuitive analysis, the SOC estimation error is presented in Figure 13. Moreover, it is worth mentioning that the time complexity of the data-driven methods such as LSSVM and the BP neural network is higher than the proposed C-WLS-EKF-based SOC estimation method.  Figure 13. Estimation error distribution of the C-WLS-EKF and intelligent algorithm. Figure 13. Estimation error distribution of the C-WLS-EKF and intelligent algorithm.

Conclusions
Two robust SOC estimation methods are introduced in this paper to solve the instability problem of the classic EKF method in the scenario where the system is affected by non-Gaussian noise. Firstly, this paper introduces the establishment process of the C-EKF and effectively solves the problem of SOC estimation under non-Gaussian noise. Secondly, the noise covariance is introduced into the C-EKF model by WLS, then the C-WLS-EKF model is established, and the digital stability and performance of the algorithm are effectively improved by the novel model. Finally, the CRLB criterion is utilized to verify the convergence of the C-WLS-EKF algorithm. The experimental results clarify that the performance of the proposed method under non-Gaussian noise conditions is better than traditional EKF, and the digital stability is effectively improved.
In the future, there are several interesting works that can be further conducted, such as parameters identification based on the robust online method, robust particle filter based on correntropy for SOC estimation, and so on.