State-of-Charge Estimation of Lithium-Ion Batteries Based on Fractional-Order Square-Root Unscented Kalman Filter

Abstract

The accuracy of the state-of-charge (SOC) estimation of lithium batteries affects the battery life, driving performance, and the safety of electric vehicles. This paper presents a SOC estimation method based on the fractional-order square-root unscented Kalman filter (FSR-UKF). Firstly, a fractional second-order Resistor-Capacitance (RC) circuit model of the lithium battery is derived. The accuracy of the parameterized model is verified, revealing its superiority over integer-order standard descriptions. Then, the FSR-UKF algorithm is developed, combining the advantages of the square-root unscented Kalman filter and the fractional calculus. The effectiveness of the proposed algorithm is proven under a variety of operational conditions in the perspective of the root-mean-squared error, which is shown to be below $1.0\%$. In addition, several experiments illustrate the performance of the FSR-UKF.

1. Introduction

Human society is responding to new energy and environmental challenges posed by climate change. The reduction of greenhouse gases emitted by vehicles is crucial to mitigate the global warming trend, and, in this perspective, the development of sustainable electric vehicles (EVs) is fundamental [1]. With the increasing applicability of EVs, the capacity, safety, health status, and endurance of the power battery became topics needing attention. The lithium-ion batteries possess high energy density, no memory, and long working life, and, consequently, they became popular as energy storage devices in EVs [2]. The battery management system (BMS) is the “intelligent element” of the EV that monitors and controls the battery [3,4,5,6]. It also collects information that is sent to the user in real-time.

The sate-of-charge (SOC) estimation of lithium-ion power batteries is fundamental to the BMS [7,8]. A stable, reliable and accurate estimation method improves the performance of the BMS. The ability to regulate the SOC within a reasonable range can prevent damage and reduction of the lifetime of the battery caused by over-charge and discharge. The design of an efficient SOC estimation method depends on the availability of accurate models able to describe accurately the electrochemical and dynamical behavior of the battery. In the description of the battery, equivalent circuit models (ECMs) are often used [9], as the Rint, Thevenin, PNGV, and high-order RC ECM [10,11]. It has been verified that ECMs with pure capacitance provide limited accuracy. Indeed, the frequency response characteristics of double-layer and pure capacitance models are inconsistent, leading to dispersion phenomena [12,13]. However, the misfitting between pure capacitance ECMs and physical behavior can be solved by adopting constant phase elements (CPE). The combination of ECMs and CPEs leads to the so-called fractional-order (FO) models of the lithium battery [14]. On the other hand, the battery itself has strong nonlinear characteristics and memory, which are well described by the tools of fractional calculus [15]. In recent years, FO ECMs have been widely used to replace traditional ECMs [16,17,18]. For example, in [19], experimental results based on a FO ECM of lithium-ion batteries and ultra-capacitors were analyzed. Other studies addressing FO ECMs of lithium batteries can be found in [16,20].

Stemming from FO ECMs of lithium batteries, a number of FO filtering methods have been proposed for SOC estimation [21,22,23,24]. For example, in [21], the FO extended Kalman filter (FEKF) was adopted to achieve SOC and state of health (SOH) estimation. The experimental results showed that the steady-state errors of the SOC and SOH estimation could be very small in the presence of initial deviation, noise, and disturbances. In [22], a FO unscented Kalman filter (FUKF) was proposed for SOC estimation. The results revealed that the model had good accuracy in estimating the battery voltage, and could estimate the SOC over a wide range of temperatures and aging levels, under distinct operation conditions of the EVs. In [23], the SOC estimation of a lithium-ion battery based on an adaptive FEKF was proposed by combining the advantages of a FO model (FOM) and an adaptive strategy. Moreover, in [24], a FEKF was shown to provide higher accuracy than classical equivalent electric circuits.

The square root unscented Kalman filter (SR-UKF) is an extension of the Kalman filter. The SR-UKF exhibits better performance in solving the nonlinear transfer problem of the mean and covariance, and avoids the error generated by the first-order partial derivative of the Taylor expansion of the nonlinear function. The square root step of the SR-UKF guarantees the positive semi-definiteness of the covariance of the state variables, and improves the stability and precision of the numerical calculations. Therefore, the SR-UKF method has been used in the SOC estimation of lithium batteries [25,26,27,28]. For example, an adaptive SR-UKF algorithm was used in [27] to improve the precision of the SOC estimation. A novel adaptive square root extended Kalman filter combined with the Thevenin ECM was proposed in [28]. The technique could solve the filtering divergence caused by computer rounding errors. A SOC estimation method based on the SR-UKF and the spherical transform was discussed in [29], showing that the scheme could reduce significantly the complexity of the calculations. However, the combination of the SR-UKF [25,26,27,28], for SOC estimation, and the tools of fractional calculus, for describing the dynamics of lithium batteries, have not been fully explored. To take full advantage of the SR-UKF and the FOM of the lithium battery, a new fractional-order square-root unscented Kalman filter (FSR-UKF) is proposed in this paper. Numerical simulations show that the SOC estimation method can yield good accuracy. The FSR-UKF has the following characteristics:

• (1). Inherits the advantages of the SR-UKF. The algorithm can be directly applied to the prediction and estimation of nonlinear systems, ensures the positive semi-definiteness of the state covariance matrix, and improves the stability of the numerical calculations.

• (2). Takes advantage of the tools of fractional calculus to describe the dynamics of lithium batteries.

• (3). Uses the new SOC estimation method, yielding better results than other schemes, namely the EKF, SR-UKF, and FUKF, as shown by tests conducted under three different temperatures and three distinct working conditions.

The paper is organized as follows. Section 2 introduces the main concepts of fractional calculus that are used in the follow-up. Section 3 formulates the fractional second-order RC ECM of the lithium battery and presents the state-space equations. Section 4 discusses the parameter identification and the accuracy of the parameterized battery model. Section 5 and Section 6 analyze the FSR-UKF algorithm and its accuracy, respectively. Section 7 outlines the main conclusions.

2. Preliminaries

The Grünwald-Letnikov (GL), Riemann-Liouville, and Caputo definitions are common in fractional calculus. In this paper, the GL definition is adopted due to its straithforward application to the Kalman filtering [30]. As such we have:

(1)$\begin{array}{c}\hfill {}_a{D}_t^{\alpha }x\left(t\right)=\underset{\Delta T\to 0}{lim}\left(\frac{1}{{(\Delta T)}^{\alpha }}\right)\sum _{j=0}^{[t/\Delta T]}{(-1)}^j\left(\begin{array}{c}\alpha \\ j\end{array}\right)x(t-j\Delta T),\end{array}$

where ${}_a{D}_t^{\alpha }$ stands for the FO operator, a and t represent the lower and upper bounds of the operator, respectively, $x\left(t\right)$ is a given function, $\Delta T$ denotes the sampling time, $[t/\Delta T]$ represents the memory length and

(2)$\begin{array}{c}\hfill \left(\begin{array}{c}\alpha \\ j\end{array}\right)=\frac{\mathsf{\Gamma }(\alpha +1)}{\mathsf{\Gamma }(j+1)\mathsf{\Gamma }(\alpha -j+1)},\end{array}$

where $\mathsf{\Gamma }\left(\alpha \right)$ is the Gamma function, defined by:

(3)$\begin{array}{c}\hfill \mathsf{\Gamma }\left(\alpha \right)={\int }_0^{+\infty }{\xi }^{\alpha -1}{e}^{-\xi }d\xi .\end{array}$

The continuous-time GL derivative can be discretized with a fixed memory length L [15], yielding:

(4)$\begin{array}{c}\hfill {}_a{D}_t^{\alpha }{x}_{k+1}\approx \frac{1}{{(\Delta T)}^{\alpha }}{x}_{k+1}+\left(\frac{1}{{(\Delta T)}^{\alpha }}\right)\sum _{j=1}^{L+1}{(-1)}^j\left(\begin{array}{c}\alpha \\ j\end{array}\right)x_{k+1-j}.\end{array}$

3. Fractional Order Modeling of Lithium-Ion Batteries

The electrochemical impedance spectroscopy (EIS) is an effective tool for determining the internal characteristics of lithium-ion batteries [31]. The process starts by applying to the device a low amplitude sine wave AC current $x\left(t\right)$, with different frequencies, and measuring the output AC signal $y\left(t\right)$, which has the same frequency as the input, but distinct amplitude, and phase. From the measurements, the transfer function $G\left(\omega \right)=Y\left(\omega \right)/X\left(\omega \right)$ can be obtained. The image of $G\left(\omega \right)$ in the complex plane represents the EIS characteristic of the lithium battery.

Figure 1 shows the Nyquist diagram of the battery obtained by the EIS procedure. We note the existence of three parts, namely the low-, intermediate- and high-frequency regions. The low-frequency region consists of a straight line with a fixed slope, representing the solid diffusion process of lithium ions inside the battery. We can use an RC branch to describe this part. The mid-frequency region is a semicircular area, which represents the charge transfer process. This part can be modeled by another RC branch. In the high-frequency region, the intersection of the impedance spectrum with the real axis corresponds to the ohmic internal resistance of the battery.

The second-order RC-ECM has been widely used to describe the behavior of lithium batteries. However, an ideal capacitor does not fully characterize the phenomena captured by the EIS. Instead of ideal capacitors, CPEs can be used. The complex impedance of a CPE is given by:

(5)$\begin{array}{c}\hfill Z\left(s\right)=\frac{1}{C{s}^{\alpha }},\end{array}$

where Z denotes the complex impedance, C is a constant representing the main capacitance effect, $\alpha$ denotes the FO, and s stands for the Laplace variable.

Hereafter, we will adopt the fractional second-order RC-ECM shown in Figure 2, where $OCV$ represents the open-circuit voltage, $R_0$ is the ohmic internal resistance, $R_1$ and $R_2$, and $C_1$ and $C_2$ denote the resistance and the capacitance in the two parallel branches, respectively, $V_0$ stands for the terminal voltage of the battery, and I corresponds to the current flowing through the circuit. The complex impedance values of the two CPEs are represented by $Z_{CPE1}\left(s\right)$ and $Z_{CPE2}\left(s\right)$, that is:

(6)$\begin{array}{c}\hfill Z_{CPE1}\left(s\right)=\frac{1}{C_1{s}^{\alpha }},Z_{CPE2}\left(s\right)=\frac{1}{C_2{s}^{\beta }},\end{array}$

where $\alpha ,\beta \in (0,1)$ are the FO.

If we denote by $V_{CPE1}\left(t\right)$ and $V_{CPE2}\left(t\right)$ the voltages on the two parallel associations, respectively, then the equations of the dynamics are:

(7)$\begin{array}{c}\hfill \left\{\begin{array}{c}\hfill D^{\alpha }{V}_{CPE1}\left(t\right)=-\frac{V_{CPE1}\left(t\right)}{R_1{C}_1}+\frac{I\left(t\right)}{C_1},\\ \hfill D^{\beta }{V}_{CPE2}\left(t\right)=-\frac{V_{CPE2}\left(t\right)}{R_2{C}_2}+\frac{I\left(t\right)}{C_2}.\end{array}\right.\end{array}$

Let us use $Q_n$ to denote the nominal capacity of the lithium battery, and $\eta$ to represent the Coulomb efficiency. The $\mathit{SOC}$ of the battery can be given by:

(8)$\begin{array}{c}\hfill \frac{d\left(\mathit{SOC}\right(t\left)\right)}{dt}=-\frac{\eta }{Q_n}I\left(t\right).\end{array}$

It follows from Kirchhoff’s voltage law that the terminals voltage $V_0\left(t\right)$ is given by:

(9)$\begin{array}{c}\hfill V_0\left(t\right)=OCV\left(\mathit{SOC}\left(t\right)\right)-I\left(t\right)R_0-V_{CPE1}\left(t\right)-V_{CPE2}\left(t\right),\end{array}$

where $OCV$ is a nonlinear function of the $\mathit{SOC}$. The relationship between the $OCV$ and the $\mathit{SOC}$ is usually expressed as:

(10)$\begin{array}{c}\hfill g\left[\theta \left(t\right)\right]=OCV\left(\mathit{SOC}\left(t\right)\right)=\sum _{k=0}^4{d}_k{\mathit{SOC}}^k\left(t\right),\end{array}$

where $d_k$$(k=0,\dots ,4)$ are the coefficients of a polynomial. Equation (10) has been used in the literature to describe the $OCV$-$\mathit{SOC}$ relationship [21,32].

Further, the state space equation of the lithium-ion battery can be written as,

(11)$\begin{array}{c}\hfill \left\{\begin{array}{c}\hfill D^{\eta }x\left(t\right)=Ax\left(t\right)+Bu\left(t\right),\\ \hfill y\left(t\right)=Cx\left(t\right)+Du\left(t\right),\end{array}\right.\end{array}$

where $\eta ={[\alpha ,\beta ,1]}^T$ is the incommensurate vector, $x\left(t\right)$ is the state vector, $D^{\eta }x\left(t\right)={[D^{\alpha }{V}_{CPE1}\left(t\right),D^{\beta }{V}_{CPE2}\left(t\right),D^1\mathit{SOC}\left(t\right)]}^T$, $u\left(t\right)$ denotes the system input (battery current $I\left(t\right)$) and $y\left(t\right)$ represents the system output (battery terminal voltage $V_0\left(t\right)$). The matrices $A,B,C$ and D are obtained as:

(12)$\begin{array}{c}\hfill A=\left[\begin{array}{ccc}-\frac{1}{R_1{C}_1}& 0& 0\\ 0& -\frac{1}{R_2{C}_2}& 0\\ 0& 0& 0\end{array}\right],B=\left[\begin{array}{c}\frac{1}{C_1}\\ \frac{1}{C_2}\\ -\frac{\eta }{Q_n}\end{array}\right],\end{array}$

(13)$\begin{array}{c}\hfill C=\left[\begin{array}{ccc}-1& -1& \frac{g\left[\theta \right(t\left)\right]}{\mathit{SOC}\left(t\right)}\end{array}\right],\end{array}$

(14)$\begin{array}{c}\hfill D=-R_0.\end{array}$

With Equation (4) in mind, we can write model (11) in the discrete time:

(15)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill x_{k+1}& =A_1{x}_k+B_1{u}_k-\sum _{j=2}^{L+1}{(-1)}^j{\gamma }_j^{\eta }{x}_{k+1-j},\hfill \\ & :=f(x_k,u_k)-\sum _{j=2}^{L+1}{(-1)}^j{\gamma }_j^{\eta }{x}_{k+1-j},\hfill \\ \hfill y_k& =g\left({\theta }_k\right)-V_{CPE1k}-V_{CPE2k}-R_0{I}_k,\hfill \\ & :=h\left(x_k\right).\hfill \end{array}\right.\end{array}$

where $f(x_k,u_k)$ stands for the system model, and $h\left(x_k\right)$ is the system measurement model. The matrices $A_1$, $B_1$ and ${\gamma }_j^{\eta }$ are as follows:

(16)$\begin{array}{c}\hfill A_1=diag({(\Delta T)}^{\alpha },{(\Delta T)}^{\beta },(\Delta T))A+diag(\alpha ,\beta ,1),\end{array}$

(17)$\begin{array}{c}\hfill B_1=diag({(\Delta T)}^{\alpha },{(\Delta T)}^{\beta },(\Delta T))B,\end{array}$

(18)$\begin{array}{c}\hfill {\gamma }_j^{\eta }=diag\left(\begin{array}{c}\left(\begin{array}{c}\alpha \\ j\end{array}\right),\left(\begin{array}{c}\beta \\ j\end{array}\right),\left(\begin{array}{c}1\\ j\end{array}\right)\end{array}\right).\end{array}$

4. Model Parameter Identification and Validation

4.1. Description of the Experimental Data

The current and voltage data of the battery INR 18650-20R, with a capacity of 2000 mAh, are provided by the CALCE Battery Research Group. We consider four current and voltage datasets for testing, which correspond to typical vehicle operating conditions [33]: Dynamic Stress Test (DST), Federal Urban Driving Schedule (FUDS), US06 Highway Driving Schedule (US06), and Beijing Dynamic Stress Test (BJDST). Their current and voltage profiles at the temperature of 25 ${}^{\circ }$C are shown in Figure 3. The current and voltage data at 0 ${}^{\circ }$C and 45 ${}^{\circ }$C, to be used in Section 5, can be found at the CALCE Battery Research Group (https://web.calce.umd.edu/batteries/index.html, accessed on 1 May 2020), being omitted here due to space constraints.

4.2. Parameter Identification

The current and voltage data under the DST conditions are used to identify the parameters of the model. The current and voltage data for the other three working conditions are used to verify the accuracy of the SOC estimation algorithm. Figure 3 shows that the DST is composed of several small cycles, each with a duration of 350 s [33]. Indeed, the DST consists of a variety of current steps with different amplitudes and lengths, and takes into account regenerative charging. In order to reduce the complexity of the analysis, we select one cycle for parameter identification. Figure 4 depicts the current and voltage data of a cycle, obtained from the 25 ${}^{\circ }$C DST dataset (data point ranging from 9430 to 9780), and with SOC values ranging from 0.549 to 0.575.

Online identification methods are more practical than their offline versions. An online method can simultaneously identify the parameters and estimate the SOC with a limited number of preliminary experiments. Meanwhile, the time-varying parameters of the lithium battery can be tackled efficiently. However, as pointed out in [22], FOMs adopt frequently off-line parameter identification methods, due to the additional complexity posed by the fractional calculus. Therefore, considering the complexity and computational burden involved, a particle swarm optimization (PSO) algorithm is employed to identify the parameters off-line. The PSO has the advantages of simplicity, fast convergence speed, and a limited number of variables.

Herein, the 12 parameters to be identified are $d_0$, $d_1$, $d_2$, $d_3$, $d_4$, $R_0$, $R_1$, $C_1$, $R_2$, $C_2$, $\alpha$ and $\beta$. In the process of parameter identification, we aim to minimize the root-mean-squared error (RMSE) between the measured and the estimated voltages. We define the objective function E as:

(19)$\begin{array}{c}\hfill minE=\sum _{k=1}^M{[V_o\left(k\right)-{\widehat{V}}_o\left(k\right)]}^2,\end{array}$

where $V_o\left(k\right)$ and ${\widehat{V}}_o\left(k\right)$ represent the measured and estimated voltages, respectively, and M is the number of sampling points in the time response.

Table 1 shows the results of the parameter identification of the FOM. Since we intend to compare the FOM with the integer-order model, we perform also the identification of the parameters of the integer second-order RC-ECM under the same conditions used for the FOM. The number of parameters to be identified for the integer-order model are merely 10, since $\alpha$ and $\beta$ are related just with the FOM. Table 2 lists the parameter identification results of the integer-order model.

4.3. Model Accuracy Verification

The accuracy of the parameters identification method is important for characterizing the quality of the battery model. Figure 5 and Figure 6 show the accuracy for the FOM identification. We verify that (i) the maximum error is 33.23 mV, (ii) the error can be kept within 30 mV most of the time, and (iii) the RMSE is 11.6 mV. Therefore, the parameters identified are accurate, and characterize adequately the battery. Figure 7 and Figure 8 illustrate the accuracy of the parameters identification method of the integer-order model. We can see that (i) the maximum error is 47.75 mV, (ii) the RMSE is 16.7 mV, and (iii) the error is kept within 40 mV most of the time.

From Figure 5, Figure 6, Figure 7 and Figure 8, it is clear that the terminal voltage estimation error of the FOM is much smaller than the one exhibited by the integer-order model. In fact, we can say that it is the main reason why the SOC estimation of the FSR-UKF is better than that of the conventional SR-UKF.

4.4. Model Parameters Sensitivity Analysis

The sensitivity analysis assesses the uncertainty in the output of a model with respect to the uncertainty in its input. After completing the parameter optimization phase, the sensitivity analysis is carried out to test the robustness of the values of the parameters. Herein, we analyze the 12 parameters of the lithium-ion battery FOM. We apply random noise in four ranges, namely [(${10}^{-5}$, ${10}^{-4}$), (${10}^{-4}$, ${10}^{-3}$), (${10}^{-3}$, ${10}^{-2}$), (${10}^{-2}$, ${10}^{-1}$)] to the parameters, and we assess their sensitivity in the perspective of the output voltage error. Three groups of experiments are carried out. The results are shown in Figure 9, Figure 10 and Figure 11. We verify that for the first three ranges of the noise, the model exhibits strong adaptability, and the output voltage error can be kept within 0.06 V. However, when the random noise perturbation is in the range $({10}^{-2},{10}^{-1})$, their impact on the voltage error becomes serious. Therefore, we conclude that random noise below ${10}^{-2}$ has little effect on the accuracy of the model.

5. SOC Estimation

The integer-order UKF (IUKF) [34] and the FUKF [22,35] have been applied in the SOC estimation of lithium batteries. Similar to the evolution verified from the IUKF towards the FUKF, we can also think of the generalization of the integer-order SR-UKF to obtain the FSR-UKF. In other words, can the accuracy of the SOC estimation be improved using the FSR-UKF? In this perspective, a FSR-UKF method is proposed.

The FSR-UKF and the FUKF are based on FOMs of lithium-ion batteries. Note that the UKF usually suffers from the ill-conditioning of the covariance matrix, because the Cholesky decomposition requires a positive definite matrix during the iterations stage. However, this requirement may result in unsatisfied with an improper selection of the initial values. In contrast, the SR-UKF does not need decomposition and reconstruction of the covariance matrix, avoiding the re-factorization of the covariance matrix at each time step. Therefore, the positive semi-definiteness of the covariance matrix can be guaranteed. Herein, the FSR-UKF method is proposed by combining the advantages of the FOM and the SR-UKF.

Consider the following FO system:

(20)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill D^{\eta }{x}_{k+1}& =f(x_k,u_k)+{\omega }_k,\hfill \\ \hfill x_{k+1}& =D^{\eta }{x}_{k+1}-\sum _{j=2}^{L+1}{(-1)}^j{\gamma }_j^{\eta }{x}_{k+1-j},\hfill \\ \hfill y_k& =h\left(x_k\right)+{\upsilon }_k,\hfill \end{array}\right.\end{array}$

where $x_k$, $u_k$ and $y_k$ represent the state, input and output system variables, respectively. The variables ${\omega }_k$ and ${\upsilon }_k$ represent the Gaussian process and measurement noise, respectively. The symbols Q and H denote the covariance matrices of ${\omega }_k$ and ${\upsilon }_k$, respectively.

The FSR-UKF algorithmic flow graph is shown in Figure 12. The detailed FSR-UKF steps are presented as follows:

• (1). Initialization

  Specify the initial state $x_o$, the matrices Q and H, and the initial state estimation error covariance $P_o$. Set the Cholesky factor of the covariance $P_o$ as $S_o$, and define $S=chol\left(P_o\right)$ ($chol$ represents the Cholesky decomposition).

• (2). Time updating

  • (a). Calculate $2n+1$ sigma sampling points:

    (21)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill {\rho }_{k-1|k-1}& =\sqrt{n+\lambda }{S}_{k-1|k-1},\hfill \\ \hfill x_{0,k-1|k-1}& ={\widehat{x}}_{k-1|k-1},\hfill \\ \hfill x_{i,k-1|k-1}& ={\widehat{x}}_{k-1|k-1}+{\rho }_{k-1|k-1}^j,i=1,2,\dots ,n,j=1,2,3,\hfill \\ \hfill x_{i,k-1|k-1}& ={\widehat{x}}_{k-1|k-1}-{\rho }_{k-1|k-1}^j,i=n+1,n+2,\dots ,2n,j=1,2,3,\hfill \end{array}\right.\end{array}$

    where n denotes the dimension of the system, and ${\rho }_{k-1|k-1}^j$ is the jth column of ${\rho }_{k-1|k-1}$. The weight of each sampling point is calculated as follows:

    (22)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill {\omega }_m^o& =\frac{\lambda }{n+\lambda },\hfill \\ \hfill {\omega }_c^o& =\frac{\lambda }{n+\lambda }+(1-{\alpha }_1^2+{\beta }_1),\hfill \\ \hfill {\omega }_m^i& ={\omega }_c^i=\frac{1}{2(n+\lambda )},i=1,2,\dots ,2n,\hfill \end{array}\right.\end{array}$

    where ${\omega }_m^i$ and ${\omega }_c^i$ are the weights of the mean of the sampling points and the covariance, respectively. The parameter $\lambda ={\alpha }_1^2(n+k)-n$ is a scaling factor that can be used to reduce the total prediction error of the system. The non-negative factor ${\alpha }_1$ generally takes a small positive value to reduce the influence of the higher-order moments, and k is the turning parameter. The non-negative weight coefficient ${\beta }_1$ is used to reduce the peak error of the estimated state and to improve the accuracy of the covariance.

  • (b). Propagate the sigma sampling points using the nonlinear function $f(·)$:

    (23)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill {\varphi }_{i,k|k-1}& =f(x_{i,k-1|k-1},u_{k-1}),i=0,1,\dots 2n,\hfill \\ \hfill D^{\eta }{\widehat{x}}_{k|k-1}& =\sum _{j=0}^{2n}{\omega }_m^i{\varphi }_{i,k|k-1}.\hfill \end{array}\right.\end{array}$

  • (c). Update the prior states estimation:

    The square root mean and covariance propagation update is:

    (24)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill S_{k|k-1}^{\Delta \Delta }& =QR\left[\left(\sqrt{{\omega }_c^i}({\varphi }_{1:2n,k|k-1}-D^{\eta }{\widehat{x}}_{k|k-1})\sqrt{Q}\right)\right],\hfill \\ \hfill S_{k|k-1}^{\Delta \Delta }& =Cholupdate(S_{k|k-1}^{\Delta \Delta },{\varphi }_{0,k|k-1}-D^{\eta }{\widehat{x}}_{k|k-1},{\omega }_c^0),\hfill \\ \hfill S_{k|k-1}^{x\Delta }& =\sum _{j=0}^{2n}(x_{i,k-1|k-1}-{\widehat{x}}_{k-1|k-1})\times ({\varphi }_{1:2n,k|k-1}-D^{\eta }{\widehat{x}}_{k|k-1}),\hfill \end{array}\right.\end{array}$

    where $QR$ denotes the $QR$ decomposition, and $Cholupdate$ represents the Cholesky factor updating.

    (25)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill {\widehat{x}}_{k|k-1}& =D^{\eta }{\widehat{x}}_{k|k-1}-\sum _{j=1,}^k{(-1)}^j{\gamma }_j{\widehat{x}}_{k-j|k-j},\hfill \\ \hfill S_{k|k-1}& =S_{k|k-1}^{\Delta \Delta }+{\gamma }_1{S}_{k|k-1}^{x\Delta }+S_{k|k-1}^{\Delta x}{\gamma }_1+\sum _{j=1}^k{\gamma }_j{S}_{k-j|k-j}{\gamma }_j.\hfill \end{array}\right.\end{array}$

• (3). Observation updating

  • (a). Update the sigma points:

    (26)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill x_{0,k|k-1}& ={\widehat{x}}_{k|k-1},\hfill \\ \hfill x_{i,k|k-1}& ={\widehat{x}}_{k|k-1}+\sqrt{n+\lambda }{S}_{k|k-1},i=1,2,\dots ,n,\hfill \\ \hfill x_{i,k|k-1}& ={\widehat{x}}_{k|k-1}-\sqrt{n+\lambda }{S}_{k|k-1},i=n+1,n+2,\dots ,2n,\hfill \end{array}\right.\end{array}$

    where the weight remains equal to the one in Equation (22).

  • (b). Propagate the sigma sampling points using the nonlinear measurement function $h(·)$:

    (27)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill {\theta }_{i,k|k-1}& =h\left(x_{i,k|k-1}\right),i=0,1,\dots ,2n,\hfill \\ \hfill {\widehat{y}}_{k|k-1}& =\sum _{j=0}^{2n}{\omega }_m^i{\theta }_{i,k|k-1},\hfill \end{array}\right.\end{array}$

  • (c). Calculate the observation-error covariance matrix:

    (28)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill S_{k|k-1}^{yy}& =QR\left[\left(\sqrt{{\omega }_c^i}({\theta }_{1:2n,k|k-1}-{\widehat{y}}_{k|k-1})\sqrt{H}\right)\right],\hfill \\ \hfill S_{k|k-1}^{yy}& =Cholupdate(S_{k|k-1}^{yy},{\theta }_{0,k|k-1}-{\widehat{y}}_{k|k-1},{\omega }_c^0),\hfill \end{array}\right.\end{array}$

  • (d). Compute the cross covariance matrix:

    (29)$\begin{array}{c}\hfill \left\{\begin{array}{c}\hfill P_{k|k-1}^{xy}=\sum _{j=0}^{2n}{\omega }_c^i({\varphi }_{i,k|k-1}-D^{\eta }{\widehat{x}}_{k|k-1}),\\ \hfill \times ({\theta }_{i,k|k-1}-{\widehat{y}}_{k|k-1}),i=0,1,\dots ,2n,\end{array}\right.\end{array}$

  • (e). Update the posterior states estimation:

    (30)$\begin{array}{c}\hfill \left\{\begin{array}{cc}\hfill K& =P_{k|k-1}^{xy}{\left(S_{k|k-1}^{yy}{\left(S_{k|k-1}^{yy}\right)}^T\right)}^{-1},\hfill \\ \hfill {\widehat{x}}_{k|k}& ={\widehat{x}}_{k|k-1}+K(y_k-{\widehat{y}}_{k|k-1}),\hfill \\ \hfill U& =K{S}_{k|k-1}^{yy},\hfill \\ \hfill S_{k|k}& =Cholupdate(S_{k|k-1},U,-1).\hfill \end{array}\right.\end{array}$

We adopt an empirical trial and error method to find the best covariance matrix. The process noise covariance matrix Q represents the statistical characteristics of the model. If the value of Q increases, then the gain matrix and the correction weight also increase, improving the dynamic performance and stability of the algorithm. If the noise covariance H increases, then the gain matrix decreases, weakening the correction weight and compromising the transient and steady-state responses. Therefore, the effects of Q and H are opposite and their values should be selected according to the actual application. In what concerns the error covariance matrix $P_o$, its size determines the amount of overshoot of the state estimate. A large $P_o$ leads to overshooting, but the estimated value evolves fast to the correct one. A small $P_o$ may result in a state estimate value that does not reach the correct one. Hereafter, we select $Q=diag\left([{10}^{-8},{10}^{-8},{10}^{-8}]\right)$, $H={10}^{-2}$ and $P_o=diag\left([{10}^{-3},{10}^{-3},{10}^{-3}]\right)$.

6. Simulation Verification and Discussion

In Section 4 we analyzed the current-voltage curves under four operating conditions at the temperature of 25 ${}^{\circ }$C. The current-voltage curve under DST conditions serves for parameter identification, while the other three are used to verify the accuracy of the SOC estimation.

Figure 13a shows the SOC estimation results under three operating conditions at 25 ${}^{\circ }$C. In all cases, the initial SOC values are selected equal to 0.7. Meanwhile, we choose the EKF, SR-UKF, and FUKF for comparison, in order to verify the FSR-UKF performance. The light green line is represented as the FSR-UKF. The red, blue, and magenta lines stand for EKF, SR-UKF, and FUKF, respectively. Also, the four algorithms are compared with the reference value represented by a black line. The closer to the reference, the higher the estimation accuracy of the algorithm. To allow good visualization of the experimental results, we use a partial magnification of the plot. Figure 13b shows the SOC estimation error associated to the EKF, SR-UKF, FUKF and FSR-UKF under the three operating conditions at 25 ${}^{\circ }$C. In order to improve the visualization of the magnitude of the estimation error, a 0 division line is added. It follows from Figure 13b that the estimation error of the EKF is the largest, followed by the ones of the SR-UKF and FUKF. On the other hand, the FSR-UKF maintains a high SOC estimation accuracy under all operating conditions.

The working environment of batteries poses some challenges, since frequently they have to operate in extreme temperatures. Therefore, the accuracy of the battery SOC estimation at extreme temperatures must be assessed. In order to verify the robustness of the FSR-UKF estimation method, we study the SOC both at low (0 ${}^{\circ }$C) and at high (45 ${}^{\circ }$C) temperatures.

6.1. Verification of the Algorithm at Low Temperatures 0 ${}^{\circ }$C

Figure 14a illustrates the SOC estimation results at 0 ${}^{\circ }$C when using the EKF, SR-UKF, FUKF and FSR-UKF, respectively. The initial SOC values under the three operating conditions are $0.7,0.75,0.7$, respectively. Figure 14a is partially magnified in order to highlight the differences. According to Figure 14a, one can see that the light green line, representing the FSR-UKF, is closer to the reference value, which also shows that the estimation accuracy of the FSR-UKF is higher. Figure 14b depicts the SOC estimation error. We verify that they are similar to those obtained at 25 ${}^{\circ }$C, and, consequently, the EKF, SR-UKF, and FUKF lead to higher estimation errors than the FSR-UKF. Most of the time, the estimation error of the FSR-UKF is kept within 0.005. The estimation errors of the other three algorithms are more than 0.005 most of the time. The FSR-UKF has small fluctuation, which shows that it can maintain good estimation accuracy and has certain stability even at low temperatures.

6.2. Verification of the Algorithm at High Temperatures 45 ${}^{\circ }$C

Figure 15a portraits the SOC estimation results under three operating conditions at 45 ${}^{\circ }$C. The initial values of the SOC are equal to 0.8. Figure 15a is also partially magnified so that we can more clearly observe which line is closer to the reference value represented by the black line. The SOC estimation errors are shown in Figure 15b. From Figure 15b, under the three working conditions, the SOC estimation error of our algorithm does not exceed 0.005 most of the time, but the estimation errors of the other three algorithms are much more than 0.005 and even 0.01. So we can verify that the FSR-UKF is better than the EKF, SR-UKF, and FUKF. The fluctuation of the FSR-UKF is smaller than that of the other three algorithms. This further verifies that our algorithm based on the FSR-UKF has certain stability at high temperatures. Therefore, both at low and high-temperature conditions, the proposed algorithm exhibits superior estimation accuracy.

Table 3, Table 4 and Table 5 list the RMSE values of the EKF, SR-UKF, FUKF and FSR-UKF estimation methods at the three temperatures (25 ${}^{\circ }$C, 0 ${}^{\circ }$C, 45 ${}^{\circ }$C). One can verify that the accuracy of the FSR-UKF is always superior to the one yielded by the other methods for all operating conditions (FUDS, US06, BJDST), and that its RMSE can always be kept within 1.0%.

From Figure 13, Figure 14 and Figure 15, we observe that the FUKF and FSR-UKF, which are based on FOMs, have smaller SOC estimation errors than the EKF and SR-UKF, which are based on integer-order models. The results prove that the FO filter methods surpass their integer-order versions. The results in Table 3, Table 4 and Table 5 show that the FSR-UKF exhibits always the smallest RMSE, which proves the effectiveness and robustness of the proposed strategy.

7. Conclusions

In this paper, a new FSR-UKF was proposed to estimate the SOC of lithium-ion batteries. The new strategy inherits the advantages of the SR-UKF, which does not require the reconstruction of the covariance matrix and, therefore, does not suffer from the ill-conditioning of that matrix, as occurs with the UKF. Moreover, the proposed method extends the SR-UKF to the FO domain, and takes full advantage of the fractional calculus to model the battery accurately. Simulation experiments were carried out at extreme temperatures (0 ${}^{\circ }$C and 45 ${}^{\circ }$C) for testing the effectiveness of the method. It was verified that the RMSE of the FSR-UKF can always be kept within 1%, demonstrating that the new scheme has good accuracy. In addition, tests with the EKF, SR-UKF, and FUKF showed the superiority of the proposed algorithm.

Author Contributions

Methodology, L.C.; software, X.W.; validation and data curation, X.D.; writing—original draft preparation, X.W.; writing—review and editing, A.M.L. and J.A.T.M.; supervision and project administration, P.L. We would like to dedicate this paper to J.A.T.M., who unfortunately passed away just before the paper was submitted for publication. J.A.T.M. played an essential role in the research described here and he is greatly missed. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Funds of China (No. 62073114, No. 11971032), Ministry of Education China Mobile Research Fund (MCM20180404), Chongqing Basic Research and Frontier Exploration Project (cstc2018jcyjAX0167) and Chongqing Outstanding Youth Fund Project (cstc2021jcyj-jqx0001).

Data Availability Statement

The data that support the findings of this study are available which has included references to the CALCE article that describes the experiments conducted for generating the data.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Nyquist diagram of a lithium battery.

Figure 2. The ECM of a lithium battery.

Figure 3. Current and voltage profiles of the DST (a), FUDS (b), US06 (c) and BJDST (d).

Figure 4. Current and voltage data of a complete DST cycle.

Figure 5. Accuracy verification of the FOM.

Figure 6. Output terminals voltage error with the FOM.

Figure 7. Accuracy verification of the integer-order model.

Figure 8. Output terminals voltage error with the integer-order model.

Figure 9. Sensitivity analysis of the model parameters (Group one).

Figure 10. Sensitivity analysis of the model parameters (Group two).

Figure 11. Sensitivity analysis of the model parameters (Group three).

Figure 12. The FSR-UKF algorithmic flow graph.

Figure 13. (a) The SOC estimation results of four algorithms under three operating conditions at 25 [Forumla omitted. See PDF.]C. (b) The SOC estimation errors of four algorithms under three operating conditions at 25 [Forumla omitted. See PDF.]C.

Figure 14. (a) The SOC estimation results of four algorithms under three operating conditions at 0 [Forumla omitted. See PDF.]C. (b) The SOC estimation errors of four algorithms under three operating conditions at 0 [Forumla omitted. See PDF.]C.

Figure 15. (a) The SOC estimation results of four algorithms under three operating conditions at 45 [Forumla omitted. See PDF.]C. (b) The SOC estimation errors of four algorithms under three operating conditions at 45 [Forumla omitted. See PDF.]C.

Table 1

The results of the FOM parameter identification.

${\mathit{d}}_0$	${\mathit{d}}_1$	${\mathit{d}}_2$	${\mathit{d}}_3$	${\mathit{d}}_4$	${\mathit{R}}_0$
$2.4877$	$1.8243$	$0.6608$	$1.1131$	$-3.2348$	$0.0687$
${\mathit{R}}_{\mathbf{1}}$	${\mathit{C}}_{\mathbf{1}}$	${\mathit{R}}_{\mathbf{2}}$	${\mathit{C}}_{\mathbf{2}}$	$\mathit{\alpha }$	$\mathit{\beta }$
0.5975	264.25	1.2679	448.54	0.4325	0.4380

Table 2

The results of the integer-order model parameter identification.

${\mathit{d}}_0$	${\mathit{d}}_1$	${\mathit{d}}_2$	${\mathit{d}}_3$	${\mathit{d}}_4$
$4.1498$	$-1.1940$	$0.0907$	$0.6634$	$-0.6488$
${\mathit{R}}_{\mathbf{0}}$	${\mathit{R}}_{\mathbf{1}}$	${\mathit{C}}_{\mathbf{1}}$	${\mathit{R}}_{\mathbf{2}}$	${\mathit{C}}_{\mathbf{2}}$
0.0788	0.8513	92.65	0.5337	75.34

Table 3

The RMSE of four algorithms under three operating conditions at 25 ${}^{\circ }$C.

RMSE	EKF	SR-UKF	FUKF	FSR-UKF
FUDS	1.09%	0.66%	0.64%	0.19%
US06	0.89%	0.34%	0.31%	0.17%
BJDST	0.92%	0.74%	0.41%	0.19%

Table 4

The RMSE of four algorithms under three operating conditions at 0 ${}^{\circ }$C.

RMSE	EKF	SR-UKF	FUKF	FSR-UKF
FUDS	1.07%	0.77%	0.49%	0.27%
US06	0.95%	0.85%	0.34%	0.16%
BJDST	0.50%	0.42%	0.40%	0.21%

Table 5

The RMSE of four algorithms under three operating conditions at 45 ${}^{\circ }$C.

RMSE	EKF	SR-UKF	FUKF	FSR-UKF
FUDS	0.93%	0.61%	0.61%	0.22%
US06	0.93%	0.82%	0.70%	0.23%
BJDST	1.11%	0.65%	0.45%	0.17%

References

1. Saw, L.H.; Ye, Y.; Tay, A.A.O. Integration issues of lithium-ion battery into electric vehicles battery pack. J. Clean. Prod.; 2016; 113, pp. 1032-1045. [DOI: https://dx.doi.org/10.1016/j.jclepro.2015.11.011]

2. Diouf, B.; Pode, R. Potential of lithium-ion batteries in renewable energy. Renew. Energy; 2015; 76, pp. 375-380. [DOI: https://dx.doi.org/10.1016/j.renene.2014.11.058]

3. Xiong, R.; Li, L.; Tian, J. Towards a smarter battery management system: A critical review on battery state of health monitoring methods. J. Power Sources; 2018; 405, pp. 18-29. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2018.10.019]

4. Xiong, R.; Zhang, Y.; Wang, J.; He, H.; Peng, S.; Pecht, M. Lithium-ion battery health prognosis based on a real battery management system used in electric vehicles. IEEE Trans. Veh. Technol.; 2019; 68, pp. 4110-4121. [DOI: https://dx.doi.org/10.1109/TVT.2018.2864688]

5. Li, P.; Zhang, Z.; Xiong, Q.; Ding, B.; Hou, J.; Luo, D.; Rong, Y.; Li, S. State-of-health estimation and remaining useful life prediction for the lithium-ion battery based on a variant long short term memory neural network. J. Power Sources; 2020; 459, 228069. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2020.228069]

6. Li, P.; Zhang, Z.; Grosu, R.; Deng, Z.; Hou, J.; Rong, Y.; Wu, R. An end-to-end neural network framework for state-of-health estimation and remaining useful life prediction of electric vehicle lithium batteries. Renew. Sustain. Energy Rev.; 2022; 156, 111843. [DOI: https://dx.doi.org/10.1016/j.rser.2021.111843]

7. Cheng, K.W.E.; Divakar, B.P.; Wu, H.; Ding, K.; Ho, H.F. Battery management system (BMS) and SOC development for electrical vehicles. IEEE Trans. Veh. Technol.; 2011; 60, pp. 76-88. [DOI: https://dx.doi.org/10.1109/TVT.2010.2089647]

8. Zou, C.; Hu, X.; Wei, Z.; Tang, X. Electrothermal dynamics conscious Lithium-ion battery cell-level charging management via state monitored predictive control. Energy; 2017; 141, pp. 250-259. [DOI: https://dx.doi.org/10.1016/j.energy.2017.09.048]

9. Hu, X.; Li, S.; Peng, H. A comparative study of equivalent circuit models for Li-ion batteries. J. Power Sources; 2012; 198, pp. 359-367. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2011.10.013]

10. Seaman, A.; Dao, T.S.; Mcphee, J. A survey of mathematics-based equivalent-circuit and electrochemical battery models for hybrid and electric vehicle simulation. J. Power Sources; 2014; 256, pp. 410-423. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2014.01.057]

11. Nejad, S.; Gladwin, D.T.; Stone, D.A. A systematic review of lumped parameter equivalent circuit models for real-time estimation of lithium-ion battery state. J. Power Sources; 2016; 316, pp. 183-196.

12. Huang, J.; Li, Z.; Liaw, B.Y.; Zhang, J. Graphical analysis of electrochemical impedance spectroscopy data in Bode and Nyquist representations. J. Power Sources; 2016; 309, pp. 82-98. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2016.01.073]

13. Zhou, X.; Huang, J.; Pan, Z.; Ouyang, M. Impedance characterization of lithium-ion batteries aging under high-temperature cycling: Importance of electrolyte-phase diffusion. J. Power Sources; 2019; 426, pp. 216-222. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2019.04.040]

14. Yang, Q.; Xu, J.; Cao, B.; Li, X. A simplified fractional order impedance model and parameter identification method for lithium-ion batteries. PLoS ONE; 2017; 12, e0172424. [DOI: https://dx.doi.org/10.1371/journal.pone.0172424]

15. Zou, C.; Zhang, L.; Hu, X.; Wang, Z.; Wik, T.; Pecht, M. A review of fractional-order techniques applied to lithium-ion batteries, lead-acid batteries, and supercapacitors. J. Power Sources; 2018; 390, pp. 286-296. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2018.04.033]

16. Wang, B.; Li, S.E.; Peng, H.; Liu, Z. Fractional-order modeling and parameter identification for lithium-ion batteries. J. Power Sources; 2015; 293, pp. 151-161. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2015.05.059]

17. Wang, Y.; Chen, Y.; Liao, X. State-of-art survey of fractional order modeling and estimation methods for lithium-ion batteries. Fract. Calc. Appl. Anal.; 2019; 22, pp. 1449-1479. [DOI: https://dx.doi.org/10.1515/fca-2019-0076]

18. Tian, J.; Xiong, R.; Shen, W.; Wang, J.; Yang, R. Online simultaneous identification of parameters and order of a fractional order battery model. J. Clean. Prod.; 2019; 247, 119147. [DOI: https://dx.doi.org/10.1016/j.jclepro.2019.119147]

19. Wang, Y.; Li, M.; Chen, Z. Experimental study of fractional-order models for lithium-ion battery and ultra-capacitor: Modeling, system identification, and validation. Appl. Energy; 2020; 278, 115736. [DOI: https://dx.doi.org/10.1016/j.apenergy.2020.115736]

20. Yu, M.; Li, Y.; Podlubny, I.; Gong, F.; Zhang, C. Fractional-order modeling of lithium-ion batteries using additive noise assisted modeling and correlative information criterion. J. Adv. Res.; 2020; 25, pp. 49-56. [DOI: https://dx.doi.org/10.1016/j.jare.2020.06.003]

21. Hu, X.; Yuan, H.; Zou, C.; Li, Z.; Zhang, L. Co-estimation of state of charge and state of health for lithium-ion batteries based on fractional-order calculus. IEEE Trans. Veh. Technol.; 2018; 67, pp. 10319-10329. [DOI: https://dx.doi.org/10.1109/TVT.2018.2865664]

22. Xiong, R.; Tian, J.; Shen, W.; Sun, F. A novel fractional order model for state of charge estimation in lithium ion batteries. IEEE Trans. Veh. Technol.; 2019; 68, pp. 4130-4139. [DOI: https://dx.doi.org/10.1109/TVT.2018.2880085]

23. Zhu, Q.; Xu, M.; Liu, W.; Zheng, M. A state of charge estimation method for lithium-ion batteries based on fractional order adaptive extended Kalman filter. Energy; 2019; 187, 115880. [DOI: https://dx.doi.org/10.1016/j.energy.2019.115880]

24. Mawonou, K.; Eddahech, A.; Dumur, D.; Beauvois, D.; Godoy, E. Improved state of charge estimation for Li-ion batteries using fractional order extended Kalman filter. J. Power Sources; 2019; 435, 226710. [DOI: https://dx.doi.org/10.1016/j.jpowsour.2019.226710]

25. Gholizade-Narm, H.; Charkhgard, M. Lithium-ion battery state of charge estimation based on square-root unscented Kalman filter. IET Power Electron.; 2013; 6, pp. 1833-1841. [DOI: https://dx.doi.org/10.1049/iet-pel.2012.0706]

26. Aung, H.; Low, K.S. Temperature dependent state-of-charge estimation of lithium ion battery using dual spherical unscented Kalman filter. IET Power Electron.; 2015; 8, pp. 2026-2033. [DOI: https://dx.doi.org/10.1049/iet-pel.2014.0863]

27. Wang, K.; Feng, X.; Ren, J.; Duan, C.; Li, L. State of charge (SOC) estimation of lithium-ion battery based on adaptive square root unscented Kalman filter. Int. J. Electrochem. Sci.; 2020; 15, pp. 9499-9516.

28. Jiang, C.; Wang, S.; Wu, B.; Fernandez, C.; Xiong, X.; Coffie-Ken, J. A state-of-charge estimation method of the power lithium-ion battery in complex conditions based on adaptive square root extended Kalman filter. Energy; 2021; 219, 119603. [DOI: https://dx.doi.org/10.1016/j.energy.2020.119603]

29. Aung, H.; Soon Low, K.; Ting Goh, S. State-of-Charge estimation of lithium-ion battery using square root spherical unscented Kalman filter (Sqrt-UKFST) in nanosatellite. IEEE Trans. Power Electron.; 2015; 30, pp. 4774-4783. [DOI: https://dx.doi.org/10.1109/TPEL.2014.2361755]

30. Monje, C.A.; Chen, Y.; Vinagre, B.M.; Xue, D.; Feliu-Batlle, V. Fractional-Order Systems and Controls: Fundamentals and Applications; Springer: London, UK, 2010.

31. Waag, W.; Kaebitz, S.; Sauer, D.U. Application-specific parameterization of reduced order equivalent circuit battery models for improved accuracy at dynamic load. Measurement; 2013; 46, pp. 4085-4093. [DOI: https://dx.doi.org/10.1016/j.measurement.2013.07.025]

32. Wang, B.; Liu, Z.; Li, S.E.; Moura, S.J.; Peng, H. State-of-Charge estimation for lithium-ion batteries based on a nonlinear fractional model. IEEE Trans. Control. Syst. Technol.; 2017; 25, pp. 3-11. [DOI: https://dx.doi.org/10.1109/TCST.2016.2557221]

33. Zheng, F.; Xing, Y.; Jiang, J.; Sun, B.; Kim, J.; Pecht, M. Influence of different open circuit voltage tests on state of charge online estimation for lithium-ion batteries. Appl. Energy; 2016; 183, pp. 513-525. [DOI: https://dx.doi.org/10.1016/j.apenergy.2016.09.010]

34. Li, B.; Bei, S. Estimation algorithm research for lithium battery SOC in electric vehicles based on adaptive unscented Kalman filter. Neural Comput. Appl.; 2018; 31, pp. 8171-8183. [DOI: https://dx.doi.org/10.1007/s00521-018-3901-7]

35. Chen, Y.; Huang, D.; Zhu, Q.; Liu, W.; Liu, C.; Xiong, N. A new state of charge estimation algorithm for lithium-ion batteries based on the fractional unscented Kalman filter. Energies; 2017; 10, 1313. [DOI: https://dx.doi.org/10.3390/en10091313]