A Dynamic State-of-Charge Estimation Method for Electric Vehicle Lithium-Ion Batteries

: With the increasing environmental concerns, plug-in electric vehicles will eventually become the main transportation tools in future smart cities. As a key component and the main power source, lithium-ion batteries have been an important object of research studies. In order to efficiently control electric vehicle powertrains, the state of charge (SOC) of lithium-ion batteries must be accurately estimated by the battery management system. This paper aims to provide a more accurate dynamic SOC estimation method for lithium-ion batteries. A dynamic Thevenin model with variable parameters affected by the temperature and SOC is established to model the battery. An unscented Kalman particle filter (UPF) algorithm is proposed based on the unscented Kalman filter (UKF) algorithm and the particle filter (PF) algorithm to generate nonlinear particle filter according to the advantages and disadvantages of various commonly used filtering algorithms. The simulation results show that the unscented Kalman particle filter algorithm based on the dynamic Thevenin model can predict the SOC in real time and it also has strong robustness against noises.


Introduction
With the development of renewable energy and the wide concern on global warming, plug-in electric vehicles become the main transportation tool in reducing transport emissions for future smart cities. For plug-in electric vehicles, the energy is mainly provided by lithium-ion batteries, which are charged through an electric power supply network that is preferably powered by clean renewable energy sources. In order to properly control the electric vehicles to meet acceleration, braking, and normal driving needs, the battery management system must properly control the charging/discharging status of the lithium-ion battery pack, and such control is usually based on accurate monitoring and estimation of key battery parameters such as the battery state of charge (SOC). SOC reflects the remaining battery capacity directly [1], which can be translated to the remaining mileage that the vehicle can travel [2]. Therefore, accurate estimation of SOC is an important requirement to manage battery power [3], predict remaining mileage [4], and avoid battery overcharging or over discharging [5].
In recent years, various methods have been proposed to estimate the SOC, including ampere-hour integral method [6], open-circuit voltage method [7], neural network method [8], and model-based estimating method. Open-circuit voltage method is simple in principle, but it cannot achieve online estimation of SOC [7]. The ampere-hour integral method is a common open-loop SOC estimation method [6], but this method requires the initial SOC value, and the error generated in the current measurement process is inevitable. In the actual SOC estimation process, its cumulative error is always increasing, which cannot satisfy the real-time and accurate SOC estimation requirement. Neural

Battery Thevenin Model
According to the principle of Thevenin equivalence, a Thevenin model whose parameters are affected by temperature T and SOC is established, as shown in Figure 1 [19]: According to Figure 1, the open-circuit voltage U oc can be calculated as follows: Energies 2020, 13, 121 3 of 16 where U 1 represents terminal voltage, R pol and C pol represent the polarization resistance and the polarization capacitance, respectively; the terminal voltage of C pol is the polarization voltage U pol . R ohm and R i are the ohmic internal resistance and the other internal resistance, respectively, and I 1 is the battery charging or discharging current. Note that C pol , R pol , R ohm , and R i are all variable parameters affected by battery SOC and temperature T. A single lithium-ion battery with a rated capacity of 9 Ah is selected, and the parameters are measured by the hybrid impulse characteristics experiment (HPPC) [20] in order to better fit the model parameters.   (1) where U1 represents terminal voltage, Rpol and Cpol represent the polarization resistance and the polarization capacitance, respectively; the terminal voltage of Cpol is the polarization voltage Upol. Rohm and Ri are the ohmic internal resistance and the other internal resistance, respectively, and I1 is the battery charging or discharging current. Note that Cpol, Rpol, Rohm, and Ri are all variable parameters affected by battery SOC and temperature T. A single lithium-ion battery with a rated capacity of 9 Ah is selected, and the parameters are measured by the hybrid impulse characteristics experiment (HPPC) [20] in order to better fit the model parameters.  Figure 2 shows the relationship among Rohm, SOC, and T measured by multiple HPPC experiments. Rohm changes simultaneously with both T and SOC when SOC is less than 0.3. This Rohm changes mainly with T when SOC is greater than 0.3, and it is slightly higher when SOC is less than 0.7. Therefore, the Rohm can be characterized when SOC is considered separately over the three different intervals, [0, 0.3], (0.3, 0.7], and (0.7, 1]. Similar to the above analysis, the relationship among Rpol, SOC, and T can be measured by HPPC experiments, as shown in Figure 3. We can see that Rpol in Figure 3 has similar variation patterns to those seen in Figure 2 Figure 2 shows the relationship among R ohm , SOC, and T measured by multiple HPPC experiments. R ohm changes simultaneously with both T and SOC when SOC is less than 0.3. This R ohm changes mainly with T when SOC is greater than 0.3, and it is slightly higher when SOC is less than 0.7. Therefore, the R ohm can be characterized when SOC is considered separately over the three different intervals, [ where U1 represents terminal voltage, Rpol and Cpol represent the polarization resistance and the polarization capacitance, respectively; the terminal voltage of Cpol is the polarization voltage Upol. Rohm and Ri are the ohmic internal resistance and the other internal resistance, respectively, and I1 is the battery charging or discharging current. Note that Cpol, Rpol, Rohm, and Ri are all variable parameters affected by battery SOC and temperature T. A single lithium-ion battery with a rated capacity of 9 Ah is selected, and the parameters are measured by the hybrid impulse characteristics experiment (HPPC) [20] in order to better fit the model parameters.  Figure 2 shows the relationship among Rohm, SOC, and T measured by multiple HPPC experiments. Rohm changes simultaneously with both T and SOC when SOC is less than 0.3. This Rohm changes mainly with T when SOC is greater than 0.3, and it is slightly higher when SOC is less than 0.7. Therefore, the Rohm can be characterized when SOC is considered separately over the three different intervals, [0, 0.3], (0.3, 0.7], and (0.7, 1]. Similar to the above analysis, the relationship among Rpol, SOC, and T can be measured by HPPC experiments, as shown in Figure 3. We can see that Rpol in Figure 3 has similar variation patterns to those seen in Figure 2 Similar to the above analysis, the relationship among R pol , SOC, and T can be measured by HPPC experiments, as shown in Figure 3. We can see that R pol in Figure 3 has similar variation patterns to those seen in Figure 2

Identification of Battery Model Parameters
In order to obtain the actual relationship between Rohm, Rpol, SOC, and T, a discrete data-fitting method can be used to fit the measured resistance values at different temperature points under different SOCs, and the fitting method is as follows.

Identification of Parameters Rohm and Rpol
The above analysis shows that Rohm and Rpol change with T and SOC simultaneously when the SOC is no more than 0.3. In order to establish a more accurate battery model while still having a simple battery model, the SOC and T can be used as input for the equivalent internal resistance model by the two-level factorial design of experimental (DOE) fitting method. Central composite design

Identification of Battery Model Parameters
In order to obtain the actual relationship between R ohm , R pol , SOC, and T, a discrete data-fitting method can be used to fit the measured resistance values at different temperature points under different SOCs, and the fitting method is as follows.

Identification of Parameters R ohm and R pol
The above analysis shows that R ohm and R pol change with T and SOC simultaneously when the SOC is no more than 0.3. In order to establish a more accurate battery model while still having a simple battery model, the SOC and T can be used as input for the equivalent internal resistance model by the two-level factorial design of experimental (DOE) fitting method. Central composite design (CCD) is the most common DOE method because of its simple design and strong system correlations. In this paper, the DOE based on CCD is applied to the internal resistance model, and the internal resistance is fitted by using T and SOC as inputs. The fitted equivalent internal resistance model includes the first-order, second-order, and cross terms, and the fitting principle is shown in Figure 4.

Identification of Parameters Rohm and Rpol
The above analysis shows that Rohm and Rpol change with T and SOC simultaneously when the SOC is no more than 0.3. In order to establish a more accurate battery model while still having a simple battery model, the SOC and T can be used as input for the equivalent internal resistance model by the two-level factorial design of experimental (DOE) fitting method. Central composite design (CCD) is the most common DOE method because of its simple design and strong system correlations. In this paper, the DOE based on CCD is applied to the internal resistance model, and the internal resistance is fitted by using T and SOC as inputs. The fitted equivalent internal resistance model includes the first-order, second-order, and cross terms, and the fitting principle is shown in Figure 4. The temperature T is the main factor that affects the internal resistance change, and the internal resistance can be fitted by the least squares method using the temperature T as the input when the SOC is greater than 0.3.
Through the above analysis, Rohm and Rpol can be fitted as shown in Equation (2). The corresponding piecewise fitted equation is as follows: Each parameter in Table 1 is the average of the same type and the same batch of multiple battery cells for many times. In practical applications, there should be a confidence interval of 95% confidence level because there will be a large number of batteries, and the same kind of battery cell parameters  The temperature T is the main factor that affects the internal resistance change, and the internal resistance can be fitted by the least squares method using the temperature T as the input when the SOC is greater than 0.3.
Through the above analysis, R ohm and R pol can be fitted as shown in Equation (2). The corresponding piecewise fitted equation is as follows: Each parameter in Table 1 is the average of the same type and the same batch of multiple battery cells for many times. In practical applications, there should be a confidence interval of 95% confidence level because there will be a large number of batteries, and the same kind of battery cell parameters are basically different, then the piecewise fitted relationship diagram of R ohm against SOC and T is shown in Figure 5. are basically different, then the piecewise fitted relationship diagram of Rohm against SOC and T is shown in Figure 5. Following a similar analysis as above, the piecewise fitted equation of Rpol and the relation between Rpol and SOC, T can be obtained as follows:  In practical application, the parameters in Table 2 should also have confidence intervals of 95% confidence level, and the relationship between Rpol and SOC, T obtained by piecewise fitting is given as Figure 6.
Rohm and Rpol are fitted by SOC and T in Figures 5 and 6 following the values measured by the HPPC test in Figures 2 and 3, respectively. Through correlation curve analysis, it is also possible to  Following a similar analysis as above, the piecewise fitted equation of R pol and the relation between R pol and SOC, T can be obtained as follows: In practical application, the parameters in Table 2 should also have confidence intervals of 95% confidence level, and the relationship between R pol and SOC, T obtained by piecewise fitting is given as Figure 6.

Identification of Parameter Cpol
Parameter τ can be identified by the charge−discharge experiments, because the polarization voltage will decay 95% after 3τ time when the battery stops discharging. Table 3 gives the measured values of τ at different typical temperatures. R ohm and R pol are fitted by SOC and T in Figures 5 and 6 following the values measured by the HPPC test in Figures 2 and 3, respectively. Through correlation curve analysis, it is also possible to verify that R ohm and R pol will change with the temperature T and battery SOC when SOC is no more than 0.3; and change mainly with temperature T when SOC is greater than 0.3. In other words, the fitted R ohm and R pol have higher accuracy. Compared with R ohm and R pol , the other internal resistance R i can be ignored in the battery model since it is relatively much smaller.

Identification of Parameter C pol
Parameter τ can be identified by the charge−discharge experiments, because the polarization voltage will decay 95% after 3τ time when the battery stops discharging. Table 3 gives the measured values of τ at different typical temperatures. The fitted equation of the parameter τ is obtained through the analysis of the measurement data, where the parameters in the equation are identified in Table 4, and the parameters in Table 4 should have a confidence level of 95% confidence level in practical application. The equation of C pol is given as follows:

SOC Estimation
To characterize the influence of temperature on SOC estimation, the SOC state equation can be expressed as follows: where SOC(t) and SOC(t − 1) in the formula represent SOC values of time t and time t − 1 respectively; λ represents a coefficient which changes under different temperatures(see Figure 7); i(τ) is the charging−discharging current; T denotes temperature; and Q full represents the nominal capacity of battery at normal temperature. In order to establish a more accurate battery model while keeping a simple model structure, the charge discharge efficiency is taken as 1 [21]. From Figure 7, λ can be identified through data fitting as follows: Its parameter values are shown in Table 5.
From Figure 7,  can be identified through data fitting as follows: Its parameter values are shown in Table 5.  According to the variable parameter Thevenin model and Equation (1), the following equation is obtained: The discrete state-space equations can be obtained by discretizing Equations (6)   According to the variable parameter Thevenin model and Equation (1), the following equation is obtained: The discrete state-space equations can be obtained by discretizing Equations (6) and (8) as follows: where x k represents the kth calculated SOC, ∆t represents the sampling period; w k and v k represent the system state noise and measurement noise, respectively; w k and v k both have zero-means and their covariances are the uncorrelated Gaussian white noise of Q and R, respectively. The observation variable y k characterizes the terminal voltage U 1 calculated by the kth calculation of the battery model, and the input variable u k characterizes the battery charge−discharge current I 1 during the kth calculation. According to the state space equation, the process for achieving the state of lithium battery SOC estimation by unscented Kalman particle filter algorithm is given as follows: (1) Using (0, 1) uniform distribution to generate N initial values of SOC, the N initial state particles and covariance are obtained as follows: Energies 2020, 13, 121 8 of 16 P 0 is the initial value of the covariance of each particle. Q 0 and R 0 are initial values of system and observation noises.
(2) The proposed state estimation principle is shown in Figure 8, and is further explained below in steps (2.1) and (2.2).  (14) 0 P is the initial value of the covariance of each particle. 0 Q and 0 R are initial values of system and observation noises.
(2) The proposed state estimation principle is shown in Figure 8, and is further explained below in steps (2.1) and (2.2). The state variables are as follows: UKF update procedure： Update and based on each sample point using UKF algorithm Figure 8. Principle of UPF (unscented Kalman particle filter) state estimation.
The state variables are as follows: The variance of state is as follows: (2.1) The specific estimation process of UPF algorithm is given as follows: An unscented Kalman filter algorithm is used for each sample point to calculate x (i)+ k and P (i)+ k . The transmission and estimation of system noise and observation noise are considered in this algorithm. The specific steps are given as follows: (2.1.1) Particle Sigma Point Sampling: At this point, the sampling point set becomes {X (i)+ j,k−1 , j = 1, 2 . . . .2M, M = n + q + r}, where q and r are the dimensions of Q and R, respectively. X j,z,k−1 are the first n dimensional column vector components, n + 1 to n + q dimensional column vector components, and n + q + 1 to n + r dimensional column vector components of the sample point respectively.
(2.1.2) Particle and estimation error covariance time update process: Energies 2020, 13, 121 9 of 16 (2.1.3) Particle measurement update process: (2.2) Resampling phase The posterior probability p i of particle x (i)+ k is calculated from the measurement result y k : This can be normalized as follows.
Thus, a set of particles x

Battery Model Validation
A single LiFePO 4 battery with rated capacity of 9 Ah was selected to conduct charge−discharge experiments at different temperatures, and the experimental test platform physical map under different conditions is shown in Figure 9.

Battery Model Validation
A single LiFePO4 battery with rated capacity of 9 Ah was selected to conduct charge−discharge experiments at different temperatures, and the experimental test platform physical map under different conditions is shown in Figure 9. This paper selected 20 cycles of dynamic stress test (DST) for verification. Figure 10a shows the DST current curve. Taking the temperature 25 °C as an example, the measurement voltage and the simulation are shown in Figure 10b. This paper selected 20 cycles of dynamic stress test (DST) for verification. Figure 10a shows the DST current curve. Taking the temperature 25 • C as an example, the measurement voltage and the simulation are shown in Figure 10b. This paper uses the mean absolute error to evaluate the fit of the battery model. The mean absolute error defined in this paper is shown as (32). The mean absolute errors are listed in Table 6.
where MAE represents the mean absolute error, and ek represents the absolute error at time k. It can be seen that the battery model established by the HPPC experimental method still has the advantages of high precision, small error by the analysis of the terminal voltage change of the lithium battery in the DST condition experiments in Figure 10, and mean absolute error of terminal voltage in Table 6, which fully verify the accuracy of the battery model.

Validation of SOC Estimation Accuracy
In order to verify the accuracy of the proposed UPF-based SOC estimation, a single LiFePO4 battery was selected, where the rated capacity is 9 Ah, and its 3 A current pulse discharge waveform is shown in Figure 11. In order to compare our UPF and existing EKF-, UKF-, and PF-based SOC estimation methods in terms of the estimation error, robustness, and other aspects by MATLAB simulations, we set the number of particles to 80 in the simulation.
Through analyzing Figures 12 and 13, it is found that the newly proposed UPF-based SOC estimation method has higher estimation accuracy and strong robustness against noises. The This paper uses the mean absolute error to evaluate the fit of the battery model. The mean absolute error defined in this paper is shown as (32). The mean absolute errors are listed in Table 6.
where MAE represents the mean absolute error, and e k represents the absolute error at time k. It can be seen that the battery model established by the HPPC experimental method still has the advantages of high precision, small error by the analysis of the terminal voltage change of the lithium battery in the DST condition experiments in Figure 10, and mean absolute error of terminal voltage in Table 6, which fully verify the accuracy of the battery model.

Validation of SOC Estimation Accuracy
In order to verify the accuracy of the proposed UPF-based SOC estimation, a single LiFePO 4 battery was selected, where the rated capacity is 9 Ah, and its 3 A current pulse discharge waveform Energies 2020, 13, 121 11 of 16 is shown in Figure 11. In order to compare our UPF and existing EKF-, UKF-, and PF-based SOC estimation methods in terms of the estimation error, robustness, and other aspects by MATLAB simulations, we set the number of particles to 80 in the simulation. battery was selected, where the rated capacity is 9 Ah, and its 3 A current pulse discharge waveform is shown in Figure 11. In order to compare our UPF and existing EKF-, UKF-, and PF-based SOC estimation methods in terms of the estimation error, robustness, and other aspects by MATLAB simulations, we set the number of particles to 80 in the simulation.
Through analyzing Figures 12 and 13, it is found that the newly proposed UPF-based SOC estimation method has higher estimation accuracy and strong robustness against noises. The maximum estimation errors at 25 °C of the estimation methods are compared in Table 7.  Through analyzing Figures 12 and 13, it is found that the newly proposed UPF-based SOC estimation method has higher estimation accuracy and strong robustness against noises. The maximum estimation errors at 25 • C of the estimation methods are compared in Table 7.   In order to verify whether the newly proposed UPF-based SOC estimation method still has higher estimation accuracy under low temperature, a comparison of different SOC estimation methods at −10 °C is provided in Figure 14.   In order to verify whether the newly proposed UPF-based SOC estimation method still has higher estimation accuracy under low temperature, a comparison of different SOC estimation methods at −10 °C is provided in Figure 14.  In order to verify whether the newly proposed UPF-based SOC estimation method still has higher estimation accuracy under low temperature, a comparison of different SOC estimation methods at −10 • C is provided in Figure 14. In order to verify whether the newly proposed UPF-based SOC estimation method still has higher estimation accuracy under low temperature, a comparison of different SOC estimation methods at −10 °C is provided in Figure 14. Based on Figure 15 and the maximum estimation error of each estimation algorithm in Table 8, it can be found that the newly proposed UPF-based SOC estimation method still has the best accuracy at low temperature. Based on Figure 15 and the maximum estimation error of each estimation algorithm in Table 8, it can be found that the newly proposed UPF-based SOC estimation method still has the best accuracy at low temperature.  In order to verify the robustness of the established UPF algorithm when the model parameters are under disturbances, the open-circuit voltage Uoc was disturbed under the condition of 3 A pulsed current discharge in Figure 11, and the SOC estimation error curve of each algorithm was obtained, as shown in Figure 16.  Based on Figure 17 and the maximum estimation error of each estimation algorithm in Table 9, it can be found that the newly proposed UPF-based SOC estimation method is more robust to interference caused by model parameter disturbances.
The above verification is based on the normal charging and discharging conditions in the laboratory, however, in the actual driving process, the working conditions are random. In order to verify the effectiveness of the proposed method in the dynamic process, the comparison verification based on the dynamic conditions was added. The federal urban driving schedule (FUDS) [22] was selected to verify the proposed method. The current and voltage profiles are shown in Figure 18.  In order to verify the robustness of the established UPF algorithm when the model parameters are under disturbances, the open-circuit voltage U oc was disturbed under the condition of 3 A pulsed current discharge in Figure 11, and the SOC estimation error curve of each algorithm was obtained, as shown in Figure 16.  In order to verify the robustness of the established UPF algorithm when the model parameters are under disturbances, the open-circuit voltage Uoc was disturbed under the condition of 3 A pulsed current discharge in Figure 11, and the SOC estimation error curve of each algorithm was obtained, as shown in Figure 16.  Based on Figure 17 and the maximum estimation error of each estimation algorithm in Table 9, it can be found that the newly proposed UPF-based SOC estimation method is more robust to Based on Figure 17 and the maximum estimation error of each estimation algorithm in Table 9, it can be found that the newly proposed UPF-based SOC estimation method is more robust to interference caused by model parameter disturbances. The FUDS test was run at 25 °C . Figures 19 and 20 represent the SOC and error curves estimated by the proposed and other methods. The maximum estimation errors at 25 °C of the estimation methods are compared in Table 10. It can be seen from the analysis that the proposed method has higher accuracy than other methods in dynamic operation conditions.
To verify the accuracy of the proposed method at low temperature, the FUDS test was run at −10 °C . The results are shown in Figures 21 and 22 and Table 11.
From the above results, it can be seen that the proposed method still has high accuracy under the condition of low-temperature operation. According to the above analysis, the proposed method still has higher estimation accuracy than other methods under dynamic conditions.
By analyzing the estimation curves and error curves of each algorithm in each state in the above simulation experiment, it can be found that the UPF-based SOC estimation has higher estimation accuracy than other commonly used methods, and it also has better anti-interference performance for the disturbance caused by model parameters.  The above verification is based on the normal charging and discharging conditions in the laboratory, however, in the actual driving process, the working conditions are random. In order to verify the effectiveness of the proposed method in the dynamic process, the comparison verification based on the dynamic conditions was added. The federal urban driving schedule (FUDS) [22] was selected to verify the proposed method. The current and voltage profiles are shown in Figure 18. The FUDS test was run at 25 °C . Figures 19 and 20 represent the SOC and error curves estimated by the proposed and other methods. The maximum estimation errors at 25 °C of the estimation methods are compared in Table 10. It can be seen from the analysis that the proposed method has higher accuracy than other methods in dynamic operation conditions.
To verify the accuracy of the proposed method at low temperature, the FUDS test was run at −10 °C . The results are shown in Figures 21 and 22 and Table 11.
From the above results, it can be seen that the proposed method still has high accuracy under the condition of low-temperature operation. According to the above analysis, the proposed method still has higher estimation accuracy than other methods under dynamic conditions. By analyzing the estimation curves and error curves of each algorithm in each state in the above simulation experiment, it can be found that the UPF-based SOC estimation has higher estimation accuracy than other commonly used methods, and it also has better anti-interference performance for the disturbance caused by model parameters. The FUDS test was run at 25 • C. Figures 19 and 20 represent the SOC and error curves estimated by the proposed and other methods. The maximum estimation errors at 25 • C of the estimation methods are compared in Table 10. It can be seen from the analysis that the proposed method has higher accuracy than other methods in dynamic operation conditions.        To verify the accuracy of the proposed method at low temperature, the FUDS test was run at −10 • C. The results are shown in Figures 21 and 22 and Table 11.

Conclusions
In order to obtain a SOC estimation method with high accuracy and wide temperature application range, the UPF algorithm and the variable parameter battery Thevenin model are adopted in this paper to propose a new SOC estimation method. This new method has the following features.
The parameters of the battery model are identified by using the parameters fitting method of CCD-based DOE and least squares method to achieve 95% confidence level according to the variations of the internal resistance. The obtained variable-parameter Thevenin model can also accurately characterize the influence of temperature, T, and SOC, thus ensure the accuracy of the battery model.
The parameter fitting method used here can obtain a more accurate battery model, based on less measured data, and reduce the workload while ensuring the accuracy of the battery model. Compared with other SOC estimation methods, the UPF-based SOC estimation method can estimate battery SOC in real time and has stronger anti-interference performance. From the above results, it can be seen that the proposed method still has high accuracy under the condition of low-temperature operation. According to the above analysis, the proposed method still has higher estimation accuracy than other methods under dynamic conditions. By analyzing the estimation curves and error curves of each algorithm in each state in the above simulation experiment, it can be found that the UPF-based SOC estimation has higher estimation accuracy than other commonly used methods, and it also has better anti-interference performance for the disturbance caused by model parameters.

Conclusions
In order to obtain a SOC estimation method with high accuracy and wide temperature application range, the UPF algorithm and the variable parameter battery Thevenin model are adopted in this paper to propose a new SOC estimation method. This new method has the following features.
The parameters of the battery model are identified by using the parameters fitting method of CCD-based DOE and least squares method to achieve 95% confidence level according to the variations of the internal resistance. The obtained variable-parameter Thevenin model can also accurately characterize the influence of temperature, T, and SOC, thus ensure the accuracy of the battery model.
The parameter fitting method used here can obtain a more accurate battery model, based on less measured data, and reduce the workload while ensuring the accuracy of the battery model. Compared with other SOC estimation methods, the UPF-based SOC estimation method can estimate battery SOC in real time and has stronger anti-interference performance.