A Novel Battery State of Charge Estimation Method Based on a Super-Twisting Sliding Mode Observer

A novel method for Li-ion battery state of charge (SOC) estimation based on a super-twisting sliding mode observer (STSMO) is proposed in this paper. To design the STSMO, the state equation of a second-order RC equivalent circuit model (SRCECM) is derived to represent the dynamic behaviors of the Li-ion battery, and the model parameters are determined by the pulse current discharge approach. The convergence of the STSMO is proven by Lyapunov stability theory. The experiments under three different discharge profiles are conducted on the Li-ion battery. Through comparisons with a conventional sliding mode observer (CSMO) and adaptive extended Kalman filter (AEKF), the superiority of the proposed observer for SOC estimation is validated.


Introduction
Due to tremendous oil consumption and aggravated environmental pollution, electric vehicles (EVs) have exhibited great promise as a new, alternative means of transportation in upcoming decades [1]. EV performance is heavily reliant on the battery characteristics. Among various power batteries for EVs, such as lead-acid and nickel-hydrogen batteries, the Li-ion battery demonstrates its great superiority in terms of high capacity and power density, fast charge capability, long service life, no memory effect, and low self-charge [2].
The state of charge (SOC) reflects the remaining capacity of the battery and plays an important role in the battery management system (BMS) which guarantees the EV's safety and reliable driving. In general, an accurate battery SOC indication is the foundation of other state estimations such as state of power (SOP). It is substantial in maximizing battery energy utilization and in preventing the battery from over-charging/over-discharging. Unfortunately, the SOC cannot be directly measured by sensors. It must be estimated by well-developed methods with the aid of measurable signals such as the voltage and current of the battery.
There have been numerous methods developed to estimate the battery SOC and each one has its own merits in some aspects. The most basic method is a coulomb counting method. It is an open-loop estimator and its accumulated error caused by measurement noise and electromagnetic interference cannot be corrected. Moreover, it also cannot offer the initial SOC value. The open circuit voltage (OCV) method is based on a one-to-one mapping from the OCV to the SOC. It provides the initial SOC value, but the battery needs to be rested for a long time to reach a steady state to obtain the OCV of the battery. The long rest period is impractical in the driving process of EVs. Therefore, the OCV method is not suitable for on-line estimation and is often utilized in combination with parameters are also given in this section. In Section 3, based on the established model, the design methodology of the super-twisting sliding mode observer (STSMO) is elaborated for SOC estimation. Then, the proposed observer is validated by experimental results in Section 4. Section 5 concludes the paper.

The Second-Order RC Equivalent Circuit Model
There have been many attempts to establish an appropriate battery model for simulating the external characteristics of a battery. The common categories are electrochemical models and equivalent circuit models (ECMs). Electrochemical models comprise a set of coupled partial differential equations which describe how the battery's potential is produced and affected by electrochemical reactions inside the battery. It is quite accurate, but its calculation is also quite complex. In general, the electrochemical model is applied to design the battery body. The ECM constructed using basic circuit components such as resistors, capacitors, and voltage sources strikes a balance between complexity and accuracy. Meanwhile, the state-space equation is easy to derive, so that the ECM can be embedded into microprocessors and can output the results in real time. Based on existing research, the higher the RC order, the more accurate the ECM is and the more complex the calculation is. In this paper, the SRCECM is chosen, as shown in Figure 1.
Energies 2018, 11, x FOR PEER REVIEW 3 of 18 model parameters are also given in this section. In Section 3, based on the established model, the design methodology of the super-twisting sliding mode observer (STSMO) is elaborated for SOC estimation. Then, the proposed observer is validated by experimental results in Section 4. Section 5 concludes the paper.

The Second-Order RC Equivalent Circuit Model
There have been many attempts to establish an appropriate battery model for simulating the external characteristics of a battery. The common categories are electrochemical models and equivalent circuit models (ECMs). Electrochemical models comprise a set of coupled partial differential equations which describe how the battery's potential is produced and affected by electrochemical reactions inside the battery. It is quite accurate, but its calculation is also quite complex. In general, the electrochemical model is applied to design the battery body. The ECM constructed using basic circuit components such as resistors, capacitors, and voltage sources strikes a balance between complexity and accuracy. Meanwhile, the state-space equation is easy to derive, so that the ECM can be embedded into microprocessors and can output the results in real time. Based on existing research, the higher the RC order, the more accurate the ECM is and the more complex the calculation is. In this paper, the SRCECM is chosen, as shown in Figure 1. The SRCECM consists of the following: (1) The capacitor Qn represents the total stored energy in the battery. The controlled voltage source Voc(SOC) is used to establish the nonlinear relationship between the SOC and the OCV, as shown in Figure 2. (2) R0 is the ohmic resistance, which indicates the internal energy loss of the battery during the charge/discharge process. (3) The parallel network Rs|| Cs represents the concentration polarization which reflects the effect of the diffusion of the reactant and product in the electrochemical reaction. Rm|| Cm represents the electrochemical polarization, which describes the structure of the double layer of the electrode/solution interface. The SRCECM consists of the following: (1) The capacitor Q n represents the total stored energy in the battery. The controlled voltage source V oc (SOC) is used to establish the nonlinear relationship between the SOC and the OCV, as shown in Figure 2. (2) R 0 is the ohmic resistance, which indicates the internal energy loss of the battery during the charge/discharge process. (3) The parallel network R s C s represents the concentration polarization which reflects the effect of the diffusion of the reactant and product in the electrochemical reaction. R m C m represents the electrochemical polarization, which describes the structure of the double layer of the electrode/solution interface. charge/discharge process. (3) The parallel network Rs|| Cs represents the concentration polarization which reflects the effect of the diffusion of the reactant and product in the electrochemical reaction. Rm|| Cm represents the electrochemical polarization, which describes the structure of the double layer of the electrode/solution interface.  In addition, due to extremely low self-discharge rate of a Li-ion battery, the self-discharge resistance is ignored in the model. V t and I represent the terminal voltage and current, respectively. I > 0 represents discharge and I < 0 represents charge. According to Kirchhoff's law, V t can be expressed as The differential expressions for the SOC, the concentration polarization voltage V s , and the electrochemical polarization voltage V m are . .
During a small range of the SOC indicated by the asterisks in Figure 2, the nonlinear relationship between the OCV and SOC can be considered linear [19]. Therefore, the OCV is expressed as a linear function of the SOC by using the piecewise linearization method: where k 0 and k 1 are the different constants for different SOC ranges. For the fast sampling frequency in SOC estimation, the current change rate . I can be ignored; namely, dI/dt = 0. We substitute Equations (2)-(5) into the derivative of Equation (1) and rearrange the differential equation: where the term ξ is added into Equation (6), which represents the differential error, the modeling error, and uncertain disturbances such as measurement noise. Solving I in Equation (1) and substituting it into Equation (2), as well as rearranging Equations (3) and (4), results in the state-space equations of the SRCECM as follows: where b 1 = k 0 /Q n + R 0 /(R s C s ) + 1/C s + R 0 /(R m C m ) + 1/C m ; b 2 = 1/C s ; b 3 = 1/C m ; a 1 = 1/(R s C s ) + 1/(R m C m ); a 2 = 1/(R 0 Q n ); a 3 = 1/(R s C s ); and a 4 = 1/(R m C m ).
The current and terminal voltage of the battery can be measured, and it is easy to represent the system of the SRCECM as: T is the state vector, and u(t) = I and y(t) = V t are the input and output of the system, respectively. The matrices A, B, and C are defined as where Γ is the constant matrix corresponding to the SOC range.
The main objective of this paper is to design an observer which is able to estimate the battery SOC based on the model given by Equation (7). The prerequisite of the successful application of the observer is that the system given by Equation (8) is observable. The observability proof is as follows: the observability matrix is The full rank of O(M) is verified by using the model parameters (derived from Section 2.2); thus, the system of the SRCECM is observable and the state variables can be estimated by the observer.

Model Parameter Identification
The experimental platform shown in Figure 3 was built to identify the model parameters and estimate the SOC. The Li-ion battery used in the experiment comprises a graphite cathode and a ternary composite anode of nickel, cobalt, and aluminum. It has a nominal capacity of 2.5 Ah and a nominal voltage of 3.6 V. The dimensions of the Li-ion battery are 18 mm diameter and 65 mm height. This kind of battery is commonly referred to as being of the "18650" format. The battery test system is used to charge or discharge the battery, which has a built-in current integrator that can be used to calculate the SOC reference. Its measurement ranges are 0~5 V and 0~5 A, and the measurement accuracy is ±5 mV/±5 mA. The upper software controls the charging or discharging mode and records the voltage and current data. The sampling interval was 10 s, and the experiment was carried out at room temperature.
Model parameter identification is used to make the model output as close to the battery output as possible utilizing the input and output of the battery with the aid of some mathematical algorithms. Currently, the off-line identification method is more widely applied, since on-line identification is sensitive to the initial values of coefficients and has stability problems due to the cross interference between the model parameters and SOC [30]. A more robust off-line identification method is to fit rest periods of experimental data with exponential functions; therefore, based on the above principle, the pulse current discharge profile as an off-line approach was applied to identify the  Figure 4a, the pulse current profile consists of many sequences of discharge and rest. Initially, the battery was fully charged (SOC = 100%) and rested for one hour to return to the equilibrium state. Then, the battery was discharged by 0.5 A, corresponding to a 0.2 C discharge ratio for 15 min. After the discharge process, the battery was rested for 1 h to allow the terminal voltage to reach the OCV. Each pulse current discharged approximately 5% of the actual capacity of the battery, which is equivalent to about 5% of the SOC. The sequence of discharge and rest was iterated until the terminal voltage was less than the cut-off voltage of 2.5 V, meaning that the battery was fully discharged (SOC = 0%). The last discharge data was abandoned because of the lack of the last resting period data. It can be seen that there were 18 transient voltage responses generated to determine 18 groups of model parameters for 18 different SOCs. Incidentally, the nonlinear relationship between the OCV and the SOC was obtained, as shown in Figure 2. The specific parameter identification process is as follows.
e system of the SRCECM is observable and the state variables can be estimated by the observer.

. Model Parameter Identification
The experimental platform shown in Figure 3 was built to identify the model parameters an timate the SOC. The Li-ion battery used in the experiment comprises a graphite cathode and nary composite anode of nickel, cobalt, and aluminum. It has a nominal capacity of 2.5 Ah and minal voltage of 3.6 V. The dimensions of the Li-ion battery are 18 mm diameter and 65 mm ight. This kind of battery is commonly referred to as being of the "18650" format. The battery te stem is used to charge or discharge the battery, which has a built-in current integrator that can b ed to calculate the SOC reference. Its measurement ranges are 0~5 V and 0~5 A, and th easurement accuracy is ±5 mV/±5 mA. The upper software controls the charging or dischargin ode and records the voltage and current data. The sampling interval was 10 s, and the experimen s carried out at room temperature. Model parameter identification is used to make the model output as close to the battery outpu possible utilizing the input and output of the battery with the aid of some mathematic orithms. Currently, the off-line identification method is more widely applied, since on-lin entification is sensitive to the initial values of coefficients and has stability problems due to th ss interference between the model parameters and SOC [30]. A more robust off-line identificatio ethod is to fit rest periods of experimental data with exponential functions; therefore, based on th ove principle, the pulse current discharge profile as an off-line approach was applied to identif e model parameters. As shown in Figure 4a, the pulse current profile consists of many sequence discharge and rest. Initially, the battery was fully charged (SOC = 100%) and rested for one hour t turn to the equilibrium state. Then, the battery was discharged by 0.5 A, corresponding to a 0.2 scharge ratio for 15 min. After the discharge process, the battery was rested for 1 h to allow th minal voltage to reach the OCV. Each pulse current discharged approximately 5% of the actu pacity of the battery, which is equivalent to about 5% of the SOC. The sequence of discharge an st was iterated until the terminal voltage was less than the cut-off voltage of 2.5 V, meaning th During each rest period of the battery, the terminal voltage response can be expressed using Equation (11): where τ s = R s C s and τ m = R m C m are the time constants of concentration polarization and electrochemical polarization, respectively. As an example, the seventh terminal voltage (actual value) response is indicated by a red square in Figure 4b where the framed section is magnified in Figure 5. The terminal voltage response profile was fitted based on Equation (11) with the Cftool Matlab fitting tool; thereby, the parameters of V s , V m , τ s , and τ m were solved. At the time of t 1 , the battery finishes the discharge of a certain sequence, and the polarization voltages are Energies 2018, 11, 1211 7 of 21 R s , C s , R m , and C m can be derived as follows: In addition, the ohmic resistance is given by The identified model parameters for the different SOCs are listed in Table 1 with the curve fitting errors represented by root mean square errors (RMSE). The SRCECM has a very small deviation in the fitting profile. For the last rest sequence, the identified parameters cause the largest model error because the polarization effect of the battery is the most serious. The fitted voltage of the model is compared with the actual voltage in Figures 4b and 6, showing the voltage error of the model. The voltage error is quite small in the range SOC ∈ [0. 1 1]. The largest error also exists around the cut-off voltage region. The model RMSE is 0.0037. To test the model's flexibility, another larger constant current (0.8 C) profile was applied to both the battery and the established model. The model voltage compared with the actual voltage of the battery is shown in Figure 7a, in which the model's voltage error has similar characteristics as in Figure 4b. The voltage error is shown in Figure 7b. The discharge end area also has large modeling error while the other error is quite small. The model RMSE is 0.0156. value) response is indicated by a red square in Figure 4b where the framed section is magnified in Figure 5. The terminal voltage response profile was fitted based on Equation (11) with the Cftool Matlab fitting tool; thereby, the parameters of , , τ , and τ were solved. At the time of , the battery finishes the discharge of a certain sequence, and the polarization voltages are , , , and can be derived as follows: In addition, the ohmic resistance is given by The identified model parameters for the different SOCs are listed in Table 1 with the curve fitting errors represented by root mean square errors (RMSE). The SRCECM has a very small deviation in the fitting profile. For the last rest sequence, the identified parameters cause the largest model error because the polarization effect of the battery is the most serious. The fitted voltage of the model is compared with the actual voltage in Figures 4b and 6, showing the voltage error of the model. The voltage error is quite small in the range SOC ∈ [0. 1 1]. The largest error also exists around the cut-off voltage region. The model RMSE is 0.0037. To test the model's flexibility, another larger constant current (0.8 C) profile was applied to both the battery and the established model. The model voltage compared with the actual voltage of the battery is shown in Figure 7a, in which the model's voltage error has similar characteristics as in Figure 4b. The voltage error is shown in Figure 7b. The discharge end area also has large modeling error while the other error is quite small. The model RMSE is 0.0156.
In addition, the ohmic resistance is given by The identified model parameters for the different SOCs are listed in Table 1 with the curve fitting errors represented by root mean square errors (RMSE). The SRCECM has a very small deviation in the fitting profile. For the last rest sequence, the identified parameters cause the largest model error because the polarization effect of the battery is the most serious. The fitted voltage of the model is compared with the actual voltage in Figures 4b and 6, showing the voltage error of the model. The voltage error is quite small in the range SOC ∈ [0. 1 1]. The largest error also exists around the cut-off voltage region. The model RMSE is 0.0037. To test the model's flexibility, another larger constant current (0.8 C) profile was applied to both the battery and the established model. The model voltage compared with the actual voltage of the battery is shown in Figure 7a, in which the model's voltage error has similar characteristics as in Figure 4b. The voltage error is shown in Figure 7b. The discharge end area also has large modeling error while the other error is quite small. The model RMSE is 0.0156.

Super-Twisting Sliding Mode Observer
In this paper, the ST algorithm was applied to design a sliding mode observer for the SOC estimation of the battery. Considering the following system of the SRCECM, and defining as the sliding mode variable and = − , the proposed SMO is formulated as

Super-Twisting Sliding Mode Observer
In this paper, the ST algorithm was applied to design a sliding mode observer for the SOC estimation of the battery. Considering the following system of the SRCECM, and defining V t as the sliding mode variable and e 1 = V t −V t , the proposed SMO is formulated as whereV t , SÔC,V s , andV m are the estimated states; and r 1 , r 2 , and r 3 are the coefficients of the . v 1 (e 1 ) = λ 1 sign(e 1 ) (16) and λ 0 and λ 1 are the gains of the ST algorithm. sign(σ) is the symbol function defined as follows: The ST algorithm is composed of two parts: The first part is a continuous function v 0 (e 1 ), which can be regarded as a proportional term. While the system error appears, it has quick response speed to reduce the error for the following reasons: when |e 1 | > 1, |e 1 | 1/2 is large, which increases the proportional gain λ 0 |e 1 | 1/2 . When 0 < |e 1 | < 1, the STSMO achieves large gain with small bias for |e 1 | 1/2 > |e 1 |. The second part is the time integral v 1 , which helps to eliminate the system error for the observed system.
We define other observation errors as Subtracting Equation (15) from Equation (7), the state estimation error system is given as follows: . . . .

Stability Proof
Assumption 1. The systems of the equivalent circuit model (Equation (7)) and the observer system (Equation (15)) are bounded input, bounded state (BIBS) since this is a physical system and the battery current and voltage are both bounded according to the battery operation conditions |I| ≤ 2.5 A, 2.5 V ≤ |V t | ≤ 4.2 V, which are controlled by the battery test system. Theorem 1. The term ξ represents the error system uncertainties, and is assumed to satisfy the following condition: . ξ ≤ δ, ξ ≤ χ and ξ = 0 for e 1 = 0 where the bounds δ and χ are positive constants which can be determined by the largest model error from experiments or experts' experience. Therefore, the state estimation error system (Equations (19)- (22)) is asymptotically stable.
Proof of Theorem 1. The proof is divided into two parts. First, Equation (19) is proven to converge to zero in finite time. Second, the resulting reduced-order dynamics of Equations (20)- (22) are proven to be exponentially stable [28].
The system given by Equation (19) can be rewritten as ξ. Equation (24) can be understood as the sliding mode dynamic property along the system trajectory. e 1 is the sliding mode variable and the relative order of the system given by Equation (19) is 1. Therefore, the ST algorithm can make . e 1 = e 1 = 0 in finite time with λ 0 > 0 and λ 1 > 0 , while e 1 = 0 can be made for the CSM.
Based on Assumption 1, the state does not go to infinity in finite time for the reason that the input I is bounded. Moreover, V t is bounded and all the states of the observer are also bounded for a finite time. Consequently, the observation error e 1 is also bounded. It follows from Equations (20) Then, the following Lyapunov function V is introduced: which can be rewritten in quadratic form where ζ T = |e 1 | 1 2 sign(e 1 ) ϕ and P = 1 2 Energies 2018, 11, 1211

of 21
Note that V is differentiable everywhere except {e 1 = 0} and is radially unbounded. The system state e 1 does not stay at {e 1 = 0} until the system converges to the origin . e 1 = e 1 = 0. With d|e 1 |/dt = . e 1 sign(e 1 ), the derivative of ζ is Taking the derivative of Equation (27) along the trajectory, where Because V = ζ T Pζ is a quadratic positive-definite function, there exists where λ min (P) and λ max (P) are the minimum eigenvalue and maximum eigenvalue of the matrix P, respectively; · 2 is a 2-norm of Euclidean space R 2 ; and ζ 2 2 = ζ 2 1 + ζ 2 2 = |e 1 | + ϕ 2 . It can be further obtained as follows: Therefore, Using Equation (32) and the fact that it follows that By the comparison principle, it follows that V and therefore ζ converge to zero in finite time. Thus, Theorem 1 is proven.
Consider a candidate Lyapunov function for the system in Equation (38): where there exists a symmetric positive-definite matrix P r which satisfies W T P + PW = −I 3×3 and I 3×3 is an identity matrix. The time derivative of V(e r ) along the trajectories of the system in Equation (38) is presented as e T r P r e r + e T r P r . e r = e T r W T P r e r + e T r P r We r = −e T r W T P r + P r W e r = − e r 2 .
According to Equation (41), . V r is always negative. Therefore, the resulting reduced-order dynamics in Equation (38) are proven to be exponentially stable.

Coefficients of the STSMO
According to the Schur lemma, matrices P and Q will be positive when satisfying the condition λ 0 > 0 and λ 1 > 0 large enough, as for the example λ 1 > χ and λ 2 0 > λ 1 where χ is explained in Equation (23), and the system errors will then converge to zero in finite time [31]. In order to make the ST algorithm convenient for application, the coefficients λ 0 and λ 1 are further discussed in [32,33]. Reference [33] gives the following sufficient conditions: where L is the Lipschitz's constant of the sliding mode variable, L > 0. Condition (42) results from a very crude estimation. In practice, one way to choose λ 0 and λ 1 is to take with properly chosen values of µ 0 and µ 1 . In particular, both µ 0 = 1, µ 1 = 1.1 and µ 0 = 4, µ 1 = 0.5 are valid choices. The other coefficients r 1 , r 2 , and r 3 of the STSMO can be determined by experiments after many rounds of debugging.

Experimental Results for SOC Estimation
The configuration of the STSMO for the SOC estimation is illustrated in Figure 8. The upper software executes the discharge command and the Li-ion battery is discharged according to the established current mode. The actual current and voltage of the battery are sampled constantly every 10 s and fed into the model to obtain the model parameters. Then, the STSMO updates the SOC based on the observation error utilizing all the parameters and sampling data.
The coefficients of STSMO were selected as λ 0 = 8 × 10 −8 , λ 1 = 7 × 10 −15 , r 1 = 5 × 10 −3 , r 2 = −7 × 10 −4 , and r 3 = −3 × 10 −4 . To further illustrate the advantages of the proposed observer, the SOC estimation results are firstly compared with the CSMO to verify the necessity of HOSM application. Then the superiority of the STSMO is proven by comparing with the well-established algorithm, namely, the AEKF proposed in [12]. The initial values of the process noise variance matrix R k , the measurement noise variance matrix Q k , and the system covariance matrix P k are R k0 = 10 −6 I, Q k0 = 10 −6 I, and P k0 = 10 −6 I, where I is an identity matrix. In addition, three different discharge profiles are utilized to evaluate the SOC estimation effect of the STSMO. The evaluation indices include SOC estimation accuracy, computation time, convergence speed, and robustness against noise disturbances. Besides the 0.5 C pulse current discharge (PCD) profile, other two discharge profiles were 0.2 C constant current discharge (CCD), and dynamic stress test (DST) discharge (shown in Figure 9). The DST discharge profile is a simplified version of the United States Advanced battery Consortium (USABC) Electric Vehicle Battery Test Procedures Manual. Before loading each discharge current, the battery is fully charged (the initial SOC = 100%) and rested for 1 h to reach the steady state.

Experimental Results for SOC Estimation
The configuration of the STSMO for the SOC estimation is illustrated in Figure 8. The upper software executes the discharge command and the Li-ion battery is discharged according to the established current mode. The actual current and voltage of the battery are sampled constantly every 10 s and fed into the model to obtain the model parameters. Then, the STSMO updates the SOC based on the observation error utilizing all the parameters and sampling data. are firstly compared with the CSMO to verify the necessity of HOSM application. Then the superiority of the STSMO is proven by comparing with the well-established algorithm, namely, the AEKF proposed in [12]. The initial values of the process noise variance matrix Rk, the measurement noise variance matrix Qk, and the system covariance matrix Pk are k0 =10 P I , where I is an identity matrix. In addition, three different discharge profiles are utilized to evaluate the SOC estimation effect of the STSMO. The evaluation indices include SOC estimation accuracy, computation time, convergence speed, and robustness against noise disturbances. Besides the 0.5 C pulse current discharge (PCD) profile, other two discharge profiles were 0.2 C constant current discharge (CCD), and dynamic stress test (DST) discharge (shown in Figure 9). The DST discharge profile is a simplified version of the United States Advanced battery Consortium (USABC) Electric Vehicle Battery Test Procedures Manual. Before loading each discharge current, the battery is fully charged (the initial SOC = 100%) and rested for 1 h to reach the steady state.   are firstly compared with the CSMO to verify the necessity of HOSM application. Then the superiority of the STSMO is proven by comparing with the well-established algorithm, namely, the AEKF proposed in [12]. The initial values of the process noise variance matrix Rk, the measurement noise variance matrix Qk, and the system covariance matrix Pk are  Figure 9). The DST discharge profile is a simplified version of the United States Advanced battery Consortium (USABC) Electric Vehicle Battery Test Procedures Manual. Before loading each discharge current, the battery is fully charged (the initial SOC = 100%) and rested for 1 h to reach the steady state.

Comparison between STSMO and CSMO
The CSMO is designed as where α is the switching gain and α = 2 × 10 −3 4 × 10 −4 −7 × 10 −5 −7 × 10 −5 T was chosen. As an example, the results of the 0.2 C CCD profile are shown in Figure 10. It can be found the SOC estimation of the CSMO has high track accuracy and proper chattering. Unfortunately, when the initial SOC deviates from its true value, e.g., 30% and 60% error, the SOC estimation errors are reduced to 5% after 930 s and 3530 s, respectively, and the convergence speed is very slow, as shown in Figure 11. In Figure 12, with relatively large gain, the CSMO achieves faster SOC estimation convergence, but significant chattering is generated on the SOC estimation results due to its discontinuous switching property. The error of Figure 12 shows that the chattering width is approximate 0.5%. As shown in Figure 13, the STSMO can eliminate the chattering significantly with reasonable convergence speed. This is due to the fact that the STSMO transfers the chattering into the HOSM surface and thereby attenuates the chattering. Therefore, compared with the CSMO, the STSMO demonstrates its superiority in terms of chattering for SOC estimation.
ors are reduced to 5% after 930 s and 3530 s, respectively, and the convergence speed is ver w, as shown in Figure 11. In Figure 12, with relatively large gain, the CSMO achieves faster SOC imation convergence, but significant chattering is generated on the SOC estimation results due t discontinuous switching property. The error of Figure 12 shows that the chattering width i proximate 0.5%. As shown in Figure 13, the STSMO can eliminate the chattering significantl th reasonable convergence speed. This is due to the fact that the STSMO transfers the chatterin o the HOSM surface and thereby attenuates the chattering. Therefore, compared with the CSMO STSMO demonstrates its superiority in terms of chattering for SOC estimation.  w, as shown in Figure 11. In Figure 12, with relatively large gain, the CSMO achieves faster SO timation convergence, but significant chattering is generated on the SOC estimation results due t discontinuous switching property. The error of Figure 12 shows that the chattering width proximate 0.5%. As shown in Figure 13, the STSMO can eliminate the chattering significantl th reasonable convergence speed. This is due to the fact that the STSMO transfers the chatterin to the HOSM surface and thereby attenuates the chattering. Therefore, compared with the CSMO e STSMO demonstrates its superiority in terms of chattering for SOC estimation.   Table  ts the RMSE of the SOC estimation at different discharge profiles. It can be seen that the RMSE o e STSMO is always smaller than that of AEKF for the three different discharge profiles.
In addition, Table 3 shows the computation time for SOC estimation by the two methods a  In addition, Table 3 shows the computation time for SOC estimation by the two methods a ferent discharge profiles. The computation time was obtained by the tic/toc program of MATLA

Comparison between STSMO and AEKF
• SOC Estimation with True Initial SOC Value and No Noise Figures 14-16 show the comparison of SOC estimation results between the AEKF and STSMO under three different discharge profiles. It can be found that both methods can track the SOC reference accurately. However, according to the error in Figures 14-16, the SOC estimation errors of the STSMO are closer to zeros while the errors of the AEKF are just bounded in by about ±5%. To quantificationally illustrate the improved accuracy of the STSMO-based SOC estimation, Table 2 lists the RMSE of the SOC estimation at different discharge profiles. It can be seen that the RMSE of the STSMO is always smaller than that of AEKF for the three different discharge profiles.
In addition, Table 3 shows the computation time for SOC estimation by the two methods at different discharge profiles. The computation time was obtained by the tic/toc program of MATLAB (2015a, MathWorks, Natick, MA, USA). The computer's CPU was an i5-4590 and the basic frequency was 3.3 GHz. After running the program many times, the maximum value of the time was chosen and presented in Table 3. It is obvious that the computation time of the STSMO is less than that of the AEKF. The reason for this is that the STSMO makes a direct state correction while the AEKF executes a priori estimation and then corrects for the system states. Therefore, the STSMO has shown outstanding capability to track the SOC reference accurately.

• SOC Estimation with Different Initial SOCs
A fast convergence speed helps to improve the capability to track the SOC reference at any initial SOC (e.g., 70% and 40%). Figure 17 presents the SOC convergence process of the STSMO and AEKF at different initial SOCs for the DST discharge profile. Table 3 gives the corresponding convergence time, which is the time that the SOC error is reduced to 5%. According to Figure 17 and Table 4, it is clear that the convergence speed of the STSMO-based SOC estimation is faster than that of the AEKF-based SOC estimation, and the SOC convergence time of the AEKF is more than 3 times that of the STSMO. Therefore, compared with the AEKF, the STSMO has faster convergence speed at certain initial SOC deviations. and Table 4, it is clear that the convergence speed of the STSMO-based SOC estimation is faster than that of the AEKF-based SOC estimation, and the SOC convergence time of the AEKF is more than 3 times that of the STSMO. Therefore, compared with the AEKF, the STSMO has faster convergence speed at certain initial SOC deviations.  In an on-line system, it is difficult to obtain accurate voltage/current values due to sensor precision and electromagnetic interference. To evaluate the robustness of the STSMO, random normally distributed noise was added to the measurement signals, i.e., voltage and current. The mean value of the noise was zero and the standard deviations (SD) SD = 0.01 and SD = 0.03 were studied. The SOC estimation results based on the STSMO and AEKF with two kinds of different noises are shown in Figure 18, and the error range is summarized in Table 5.
Obviously, the maximum errors of both the STSMO and AEKF increase with the increase of the SD, but the values of the STSMO increase far less than those of the AEKF. The STSMO can also  In an on-line system, it is difficult to obtain accurate voltage/current values due to sensor precision and electromagnetic interference. To evaluate the robustness of the STSMO, random normally distributed noise was added to the measurement signals, i.e., voltage and current. The mean value of the noise was zero and the standard deviations (SD) SD = 0.01 and SD = 0.03 were studied. The SOC estimation results based on the STSMO and AEKF with two kinds of different noises are shown in Figure 18, and the error range is summarized in Table 5. track the SOC reference with small error, bounded by about ±2%. In addition, it can be found that the maximum of the AEKF becomes unexpectedly smaller when the SD is from 0 to 0.01. The reason for this is that it hardly works for relatively accurate measurement signals such as laboratory data, although AEKF has strong robustness against low levels of noise. On the whole, the STSMO has stronger robustness against noise disturbance.
(a) (b)    Obviously, the maximum errors of both the STSMO and AEKF increase with the increase of the SD, but the values of the STSMO increase far less than those of the AEKF. The STSMO can also track the SOC reference with small error, bounded by about ±2%. In addition, it can be found that the maximum of the AEKF becomes unexpectedly smaller when the SD is from 0 to 0.01. The reason for this is that it hardly works for relatively accurate measurement signals such as laboratory data, although AEKF has strong robustness against low levels of noise. On the whole, the STSMO has stronger robustness against noise disturbance.

Conclusions
On the basis of the established SRCECM, the STSMO was presented to estimate battery SOC. Lyapunov stability analysis was utilized to prove the convergence of the error system of the observer. Three different discharge profiles-PCD, CCD, and DST-were carried out to evaluate the performance of the SOC estimation. The experimental comparisons with the CSMO showed that the STSMO can solve the chattering problem and improve the convergence speed, and the superiority of the HOSM application is thereby verified. Furthermore, the comparisons with the AEKF indicate that the STSMO has better performance in improving the estimation accuracy, reducing the computation time, accelerating the convergence speed and enhancing the robustness against noise disturbances. The proposed STSMO is suitable for on-line SOC estimation.
Author Contributions: All authors have worked on this manuscript together and all authors have read and approved this manuscript.