# Prognostics for State of Health of Lithium-Ion Batteries Based on Gaussian Process Regression

Authors: Di Zhou; Hongtao Yin; Ping Fu; Xianhua Song; Wenbin Lu; Lili Yuan; Zuoxian Fu
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Accurate estimation and prediction of the lithium-ion (Li-ion) batteries’ performance has important theoretical and practical significance to make better use of lithium-ion battery and to avoid unnecessary losses. State of health (SOH) estimation is used as a qualitative measure of the capability of a lithium-ion battery to store and deliver energy in a system. To evaluate and predict the SOH of batteries, the Gaussian process regression with neural network (GPRNN) as its variance function is proposed. Experimental results confirm that the proposed method can be effectively applied to Li-ion battery monitoring and prognostics by quantitative comparison with basic GPR, combination LGPFR, combination QGPFR, and the multiscale GPR (SMK-GPR, P-MGPR, and SE-MGPR). The criteria of RMSE and MAPE of the proposed three models are reduced significantly compared to those of other existing methods.


---

## Body

## 1. Introduction

Lithium-ion battery has been widely applied to portable electronic devices, electric and hybrid vehicles, even military electronics, aerospace avionics, and other automotive vehicles because of the properties of lightness, high energy density, high output voltage, long lifetime, and low self-discharge [1, 2]. The safety and reliability of Li-ion battery in actual application are critical factors of the system performance whose failure would not only lead to the performance degradation, repair costs increase, and inconvenience, but also cause catastrophic accidents due to overheating and short circuiting [3]. Therefore, the effective monitoring and prognostics for Li-ion battery have been attracted much more attention and a reliable and effective BMS (battery management system) is indispensable [4, 5]. The precondition of BMS is monitoring and estimation of states which includes estimations of SOC (state of charge), SOH (state of health), and SOL (state of life). The above three states estimations are still the core content in the research of BMS [4, 6–8]. SOH is the ratio of the released capacity with a certain discharge from full charge state to cut-off standard voltage under standard condition and its corresponding initial capacity which is the response of the battery performance degradation and the accidents possibility. In general, SOC describes the short-term changes of the current parameters and SOH describes long-term changes. SOH does not require continuous measurements and only needs periodical measurements in most cases, and the measured periods depend on the applications. Based on the above advantages of SOH, this paper focuses on research of SOH to reflect the performance of Li-ion battery.With the development of Li-ion battery, SOH estimation plays an important role as a quantitative description for the battery performance indicators to store and deliver energy in the system [9, 10]. Generally, there are two mainstream methods to estimate SOH [11]. One is experience-based approach such as cycle counting method, Ah Low and Ah weighted method, and event-oriented cumulative method of aging. The other method is performance-based approach which utilizes different forms of performance models and considers the aging process and stress factors. According to different sources of information used, the performance-based approach is divided into three categories which are mechanism-based method, characteristics-based method, and data-driven method. Due to the complex physical and chemical processes of the battery itself, a lot of the rules are difficult to describe directly through the mechanism. The method which describes the battery performance from the perspective of the test data is called data-driven approach. The acquired data tends to have a strong uncertainty and incompleteness; testing all the possible factors affecting the battery life in the practical application is unrealistic. Although the limitations exist, data-driven method is still an effective prediction method. It is used to excavate the rules of battery performance for SOH prediction through using battery performance test data. Data-driven prediction does not require mechanical knowledge of object systems. It excavates the implied information to predict SOH through a variety of data analysis learning methods. Therefore, it can avoid the difficulty of the model acquisition, so it is a more practical prediction method.There are many common data-driven algorithms such as Particle Filtering (PF), Support Vector Machine (SVM), and Relevant Vector Machine (RVM). PF can combine available information from system measurements and analytic models for SOH prediction [10, 12]. Kalman filter and extended Kalman filter (EKF) are applied to SOH prediction of batteries by using estimation of equivalent circuit model [13, 14]. Unscented filtering and Bayesian filtering are classical methods for battery SOH or SOC estimation. These stochastic filtering methods for prognostics are easier to be obtained and have wider application scope. However, the methods are difficult in containing the environmental interference and loading dynamic characteristics factors and so on, and the dynamic accuracy and adaptability are still limited. Moreover, PF algorithm itself solves the problem of parameters initialization which is time-consuming in the case of more parameters. Some methods are mostly improvement and optimization of PF or KF such as sequential Monte Carlo method referring to PF which has attracted attention to SOH [15]. In recent years, SVM has been used in SOH prediction of Li-ion battery [16]. In addition to basic SVM algorithm, there are some improved SVM algorithms and fusion algorithms. SVM method still has some deficiencies: the kernel function must satisfy the Mercer condition, sparsity is limited, the number of support vectors is sensitive to error boundary, and the loss of function and the penalty factor are difficult to determine. SVM also lacks the expression of uncertainty management. To solve this problem, RVM approaches have been proposed. RVM algorithm has higher calculation accuracy and lower computational complexity compared with SVM and can output the uncertainty information of prediction results [17]. Kim obtained the relevant parameters of exponential model by using EIS test data based on RVM algorithm and gave the prediction outcomes [1].  Oai et al. developed a new technique to predict SOH using a dual-sliding-mode observer. The researchers presented a method, Li-ion battery SOH prediction concept, based on an appropriate SOH definition [18]. However, since the much sparse and dynamic fluctuation of capacity data will lead to forecasting the result of poor stability if RVM is directly used for prediction only, there are many others methods such as neural fuzzy logic [19], networks [20], regressions [21], and distributed active learning [22] that lack specific strategies in dealing with model uncertainties and would lead to prognostics performance deterioration when facing long-term changes in environmental and operating conditions.To consider modeling flexibility and uncertainty representations, Gaussian process regression (GPR) has been investigated for Li-ion battery prognostics, where the degradation trends are learnt from battery datasets with the combination of Gaussian process functions [23]. GPR could provide a confidence measure of output and has been used for SOH predictions as a kind of new machine learning method based on Bayesian theory and statistical learning theory. Meanwhile, it is a nonparametric method which is suitable for high dimension nonlinear regression problems and can give the uncertainty representation of SOH prediction. An improved GPR method has been proposed to capture regeneration phenomena [24]. Li and Xu proposed a new prognostics method for state of health estimation of Li-ion batteries based on a mixture of Gaussian process models and particle filter proposed [25]. However, forecasting accuracy and uncertainty of the above two methods do not meet the requirements of practical application.In this paper, an approach based on GPR with neural network kernel function is presented to predict the SOH of Li-ion batteries. Neural network is ideally suited to the SOH data, since it has smaller misspecification than other stationary covariance functions and draws saturation at different values in the positive and negative directions of variablex [23]. Meanwhile, the log marginal likelihood is much larger than other covariance functions. The proposed method consists of three models: the neural network itself, the sum of the neural network and the Maternard covariance function, and the product of the neural network and the periodic covariance function. Finally, experiments based on the NASA battery datasets are provided to demonstrate the performance of the proposed three prognostics models. The experimental results show that the predictions are near perfect.The remainder of this paper is organized as follows. In Section2, the background knowledge about the SOH of Li-ion batteries is described. Then, the proposed methodology including GPR and its parameters is provided in Section 3. Experiments and analysis are given in Section 4 to demonstrate the performance of proposed prognostics algorithm. Finally, conclusions are drawn in Section 5.


## 2. The SOH of Lithium-Ion Batteries

There are four definitions of SOH according to the battery characteristics which are shown as follows [26].(1)
From the perspective of remaining power of battery to determine the battery SOH, its definition can be given by(1)SOH=QagedQnew×100%,where Qaged is the current biggest power of battery and Qnew is the initial power.(2)
From the perspective of starting power to define the battery SOH, the expression is given by(2)SOH=CCAocmp-CCAmin⁡CCAnew-CCAmin⁡×100%,where CCAocmp is the real-time starting power and CCAnew is the predicted starting power released by battery when the SOC is 100% and CCAmin⁡ is the minimum starting power needed.(3)
From the perspective of impedance measurement to define the battery SOH, the definition is depicted as follows:(3)SOH=RiR0×100%,where Ri is ith impedance measurement varying with the cycles of charging and discharging and R0 is the initial impedance.(4)
From the perspective of battery capacity power to define the battery SOH, this method can be expressed as(4)SOH=CiC0×100%,where Ci is the ith capacitance value degenerated with cycles and C0 is the initial capacity.In this paper, this definition of SOH using the capacity power is utilized to study the performance of battery.


## 3. Methodology

### 3.1. Gaussian Process Regression

Gaussian process (GP) is used to describe random variables which are scalars or vectors. GP is a stochastic process which also governs the properties of functions and gets Bayesian inferences with the function-space view. GP is a set where any finite random variable has joint Gaussian distribution, and the characters of GP are completely determined by the mean function and the covariance function [28].(5)mx=Efx,kx,x′=Emx-fxmx′-fx′.Consider the simple regression model(6)y=fx+ε,fx=XTw.The training set is {xi,yi∣i=1,…,n}, where xi,i=1,…,n, are the input values, f is a function value, and yi, i=1,…,n, are the output values. ε represents noise which follows an independent, identically distributed Gaussian distribution with zero mean and variance σ2; that is,(7)ε~N0,σ2.Then the prior distribution of the observation y is obtained:(8)y∣X,w~NXTw,kx,x′∣σ2I.It can be used to estimate the predictive values. The joint distribution of the observation y and the predictive values f∗ is derived as follows:(9)yf∗~N0,kX,XkX,x∗kx∗,Xkx∗,x∗,where kX,x∗ represents the covariance of the training data points and predictive values and kx∗,x∗ is the covariance of the predictive values.Computed by Bayes’ rule, the posterior distribution over the weights can be gotten. The posterior distribution to inference in the Bayesian linear model can be shown as(10)f∗∣X,f,x∗~Nkx∗,XkX,X-1f,kx∗,x∗-kx∗,XkX,X-1kX,x∗.


### 3.2. Choice of Parameters

GP is usually parameterized in terms of their covariance functions. There are some most popular covariance functions [23, 29]: the squared exponential covariance function,(11)kx,x′=σf2exp⁡-xi-xj22l2,the constant covariance function,(12)kx,x′=σ02,and the periodic covariance function,(13)kx,x′=σf2-2l2sin2ϖ2πxi-xj.These covariance functions are stationary and are usually used in battery health prognostic with good results. Here in this paper a particular type of neural network is applied in battery health prognostic, which is called neural network covariance function; see (7). Neural network is ideally suited to the SOH data, since it allows saturation at different values in the positive and negative directions of x [23]. For comparison, the predictive distributions for three covariance functions are shown in Figure 1. The 64 data points are produced by a step function with Gaussian noise with standard deviation σ=0.1 which is totally the same as the example taken in [23]. Figure 1 shows the means and 95% confidence intervals for the noisy signal in grey. Particularly, Figure 1(a) expresses the single squared exponential covariance and Figure 1(b) illustrates the sum of two squared exponential covariances. Obviously, this covariance function is more flexible than a single squared exponential, since it has two magnitudes and two length-scale parameters. The predictive distribution looks a little bit better, but they both are not ideal for this scenario. Figure 1(c) shows the neural network covariance function which is ideally suited to this case, since it allows saturation at different values in the positive and negative directions of x. As shown in Figure 1, the predictions are also near perfect.Figure 1
Misspecification example.
(a)
Squared exponential (b)
Sum of two squared exponentials (c)
Neural networkNew kernels can be generated with this kernel; for example, the sum of the neural network covariance function and the Maternard covariance function is a new kernel and the product of the neural network covariance function and the periodic covariance function is also a new kernel. The Maternard covariance function is stationary and nondegenerate and the periodic covariance function is selected to approximate the local regeneration phenomenon of SOH data. Therefore, these two covariance functions are used in combination with the neural network in this paper. According to the covariance function chosen in our paper, the simulation results are better than others’. It turns out that the neural network covariance function can be viewed as a suitable choice of this type of data. The covariance functions used in the proposed paper are as follows [23].The neural network covariance function:(14)k1x,x′=σf2sin-1xTΛ-2x′1+xTΛ-2x1+x′TΛ-2x′.The Maternard covariance function:(15)k2x,x′=σf2fdrdexp⁡-rd.The periodic covariance function:(16)k3x,x′=σf2exp-2l2sin2πx-x′p.In the Maternard covariance function,f is a function of rd,  f1(t)=1,  f3(t)=1+t,  f5t=f3t+t2/3,  rd is the distance between x and x′, and  rd=dx-x′Λ-2x-x′.Typically, the used covariance function has some free parameters which are generally called free hyperparameters. As these parameters vary, the predictions of data are different which will lead to great fitting error. In order to get the suitable hyperparameters, optimization with the maximization of the log-likelihood function is necessary [28, 30]:(17)log⁡py∣X=-12yTK+σ2I-1y-12log⁡K+σ2I-n2log⁡2π.The hyperparameters can be determined from training data. Notice that the length-scale l in the kernel varies, and the signal variance σf2 and the noise variance σ2 would vary. Through experiments, it can be found that when different initial values are set, different fitting curves are obtained and the error between them would be great. According to the proposed issue, an appropriate initial value is chosen.A larger number of training data enable GPR model to get better prognostic prediction. According to the experimental data, the linear mean functionm(x)=ax+b is selected to fit the track of the data. The covariance functions discussed above, that is, the neural network shown as (7), the sum of neural network and Maternard given as kx,x′=k1x,x′+k2x,x′, and the product of neural network and periodic displayed as kx,x′=k1x,x′k3x,x′, are considered. The three covariance functions are denoted as Model I, Model II, and Model III. Then the hyperparameters spaces in the three models are Θ1=a,b,l,sfT, Θ2=a,b,l,sf1,ell,sf2T, and Θ3=a,b,l,sf1,ell,p,sf2T, where a,b are coefficients of the mean function, Λ is l times the unit matrix, and sf, sf1, and sf2 control the variance σf2 in three modes, respectively. These three models are all effective and are designed to fit SOH data and make the prediction more accurate.


## 3.1. Gaussian Process Regression

Gaussian process (GP) is used to describe random variables which are scalars or vectors. GP is a stochastic process which also governs the properties of functions and gets Bayesian inferences with the function-space view. GP is a set where any finite random variable has joint Gaussian distribution, and the characters of GP are completely determined by the mean function and the covariance function [28].(5)mx=Efx,kx,x′=Emx-fxmx′-fx′.Consider the simple regression model(6)y=fx+ε,fx=XTw.The training set is {xi,yi∣i=1,…,n}, where xi,i=1,…,n, are the input values, f is a function value, and yi, i=1,…,n, are the output values. ε represents noise which follows an independent, identically distributed Gaussian distribution with zero mean and variance σ2; that is,(7)ε~N0,σ2.Then the prior distribution of the observation y is obtained:(8)y∣X,w~NXTw,kx,x′∣σ2I.It can be used to estimate the predictive values. The joint distribution of the observation y and the predictive values f∗ is derived as follows:(9)yf∗~N0,kX,XkX,x∗kx∗,Xkx∗,x∗,where kX,x∗ represents the covariance of the training data points and predictive values and kx∗,x∗ is the covariance of the predictive values.Computed by Bayes’ rule, the posterior distribution over the weights can be gotten. The posterior distribution to inference in the Bayesian linear model can be shown as(10)f∗∣X,f,x∗~Nkx∗,XkX,X-1f,kx∗,x∗-kx∗,XkX,X-1kX,x∗.


## 3.2. Choice of Parameters

GP is usually parameterized in terms of their covariance functions. There are some most popular covariance functions [23, 29]: the squared exponential covariance function,(11)kx,x′=σf2exp⁡-xi-xj22l2,the constant covariance function,(12)kx,x′=σ02,and the periodic covariance function,(13)kx,x′=σf2-2l2sin2ϖ2πxi-xj.These covariance functions are stationary and are usually used in battery health prognostic with good results. Here in this paper a particular type of neural network is applied in battery health prognostic, which is called neural network covariance function; see (7). Neural network is ideally suited to the SOH data, since it allows saturation at different values in the positive and negative directions of x [23]. For comparison, the predictive distributions for three covariance functions are shown in Figure 1. The 64 data points are produced by a step function with Gaussian noise with standard deviation σ=0.1 which is totally the same as the example taken in [23]. Figure 1 shows the means and 95% confidence intervals for the noisy signal in grey. Particularly, Figure 1(a) expresses the single squared exponential covariance and Figure 1(b) illustrates the sum of two squared exponential covariances. Obviously, this covariance function is more flexible than a single squared exponential, since it has two magnitudes and two length-scale parameters. The predictive distribution looks a little bit better, but they both are not ideal for this scenario. Figure 1(c) shows the neural network covariance function which is ideally suited to this case, since it allows saturation at different values in the positive and negative directions of x. As shown in Figure 1, the predictions are also near perfect.Figure 1
Misspecification example.
(a)
Squared exponential (b)
Sum of two squared exponentials (c)
Neural networkNew kernels can be generated with this kernel; for example, the sum of the neural network covariance function and the Maternard covariance function is a new kernel and the product of the neural network covariance function and the periodic covariance function is also a new kernel. The Maternard covariance function is stationary and nondegenerate and the periodic covariance function is selected to approximate the local regeneration phenomenon of SOH data. Therefore, these two covariance functions are used in combination with the neural network in this paper. According to the covariance function chosen in our paper, the simulation results are better than others’. It turns out that the neural network covariance function can be viewed as a suitable choice of this type of data. The covariance functions used in the proposed paper are as follows [23].The neural network covariance function:(14)k1x,x′=σf2sin-1xTΛ-2x′1+xTΛ-2x1+x′TΛ-2x′.The Maternard covariance function:(15)k2x,x′=σf2fdrdexp⁡-rd.The periodic covariance function:(16)k3x,x′=σf2exp-2l2sin2πx-x′p.In the Maternard covariance function,f is a function of rd,  f1(t)=1,  f3(t)=1+t,  f5t=f3t+t2/3,  rd is the distance between x and x′, and  rd=dx-x′Λ-2x-x′.Typically, the used covariance function has some free parameters which are generally called free hyperparameters. As these parameters vary, the predictions of data are different which will lead to great fitting error. In order to get the suitable hyperparameters, optimization with the maximization of the log-likelihood function is necessary [28, 30]:(17)log⁡py∣X=-12yTK+σ2I-1y-12log⁡K+σ2I-n2log⁡2π.The hyperparameters can be determined from training data. Notice that the length-scale l in the kernel varies, and the signal variance σf2 and the noise variance σ2 would vary. Through experiments, it can be found that when different initial values are set, different fitting curves are obtained and the error between them would be great. According to the proposed issue, an appropriate initial value is chosen.A larger number of training data enable GPR model to get better prognostic prediction. According to the experimental data, the linear mean functionm(x)=ax+b is selected to fit the track of the data. The covariance functions discussed above, that is, the neural network shown as (7), the sum of neural network and Maternard given as kx,x′=k1x,x′+k2x,x′, and the product of neural network and periodic displayed as kx,x′=k1x,x′k3x,x′, are considered. The three covariance functions are denoted as Model I, Model II, and Model III. Then the hyperparameters spaces in the three models are Θ1=a,b,l,sfT, Θ2=a,b,l,sf1,ell,sf2T, and Θ3=a,b,l,sf1,ell,p,sf2T, where a,b are coefficients of the mean function, Λ is l times the unit matrix, and sf, sf1, and sf2 control the variance σf2 in three modes, respectively. These three models are all effective and are designed to fit SOH data and make the prediction more accurate.


## 4. Experiments and Analysis

### 4.1. Raw Data Description of Li-Ion Batteries

To comprehensively illustrate the efficiency of the proposed method, the dataset about the aging of Li-ion batteries is selected which is obtained from the data repository of the NASA Ames Prognostics Center of Excellence (PCoE) [31].The dataset of the NASA Ames PCoE consists of 36 Li-ion batteries’ data: No. 5, No. 6, No. 7, No. 18, and Nos. 25–56. Because the lengths of SOHs extracted from the 36 batteries are different, Li-ion batteries of No. 5, No. 6, and No. 7 are the most commonly used in related literatures [24, 25, 27] because they are much longer. They were run through 3 different operational profiles (charge, discharge, and impedance) at room temperature. Charging was carried out in a constant current (CC) mode at 1.5 A until the battery voltage reached 4.2 V and then continued in a constant voltage (CV) mode until the charge current dropped to 20 mA. Discharge was carried out at a constant current (CC) level of 2 A until the battery voltage fell to 2.7 V, 2.5 V, 2.2 V, and 2.5 V for batteries Nos. 5, 6, 7, and 18, respectively. Impedance measurement was carried out through an electrochemical impedance spectroscopy (EIS) frequency sweep from 0.1 Hz to 5 kHz. Repeated charge and discharge cycles result in accelerated aging of the batteries while impedance measurements provide insight into the internal battery parameters that change as aging progresses. The experiments were stopped when the batteries reached end-of-life (EOL) criteria, which was a 30% fade in rated capacity (from 2 Ahr to 1.4 Ahr). This dataset can be used for the prediction of SOC, SOH, and the remaining useful life (RUL) of batteries [23].In this paper, the data from batteries No. 5, No. 6, and No. 7 is processed using (4) to implement the experiments including training and testing. Figure 2 shows the SOH curves of these three batteries. The total charge/discharge cycles are all 168. Figure 2 indicates that the SOH shows an obvious global degradation trend and local regeneration phenomenon.Figure 2
SOHs of Li-ion batteries No. 5, No. 6, and No. 7.


### 4.2. Training with GPR

For No. 5, No. 6, and No. 7, the first 100 SOH points are utilized to train the proposed GPR model and the remaining 68 points are used for prediction. For the training data from the above three batteries, the mean function is chosen as the linear functionm(x)=ax+b and the covariance function is chosen as three functions discussed above; according to these three models, the algorithm for the training and prediction based on GPR is produced as follows:(a)
Determine the training datasetsi,xii=1N, where i is the number of charge/discharge cycles and xi is the corresponding value of SOH at the ith charge/discharge cycle. For batteries No. 5, No. 6, and No. 7, N equals 100.(b)
Initialize the hyperparameters included in the mean function and covariance function for the different training sets. Here, the initial hyperparameters for Model I areΘ1=a,b,l,sfT=0.1,0,rand,2T, for Model II are Θ2=a,b,l,sf1,ell,sf2T=0.1,0,rand,2,0.9,2T, and for Model III are Θ3=a,b,l,sf1,ell,p,sf2T=0.5,1,rand,2,0.9,2,2T, wherein signal “rand” which is the value of parameter “l” in these three models is a random number between 0 and 1.To compare with other methods, the hyperparameter in basic GPR is selected asΘ′=l,sfT=1,1T and in combination LGPFR is Θ′′=a,b,sf1,l,sf2,ell,ωT=0,0,0.1,1,0.2,1,5T which is the same as the hyperparameters used in [24].(c)
Optimize the hyperparameters with the maximization of the log-likelihood function.(d)
Set the prediction stepsm, for No. 5, No. 6, and No. 7, m=68.(e)
Apply the training datai,xii=1N and the testing data i,xii=N+1N+m to Models I, II, and III with the optimal hyperparameters. Then the prediction outputs will be derived.


### 4.3. Prediction and Comparison

In this part, the prediction results of five GPR models, basic GPR, combination LGPFR, and Models I, II, and III for batteries No. 5, No. 6, and No. 7, are shown in Figures3–5, respectively. As shown in Figure 3, the red lines are the actual SOH value, the blue scatter plots are the prediction result, and the grey regions represent the 95% confidence intervals. Taking Figure 3, for example, for the method of basic GPR shown in Figure 3(a), the mean prediction output is further away from the actual SOH for the subsequent cycles. Meanwhile, the 95% confidence intervals increase significantly. The above two aspects demonstrate the poor prediction effect of basic GPR. For the method of combination LGPFR shown in Figure 3(b), although the prediction values are close to the actual SOH data points, the 95% confidence intervals are still wide. The wide confidence intervals indicate the high uncertainty of the prediction result. Experiments have shown the limitations of Li-ion batteries’ SOH prediction based on the basic GPR model and the combination LGPFR model. The simulations with Models I, II, and III for No. 5 are shown in Figures 3(c)–3(e); it is clear that both the point prediction results and the uncertainty representations are improved greatly.Figure 3
Prediction results of five GPR models for battery No. 5.
(a)
Basic GPR (b)
Combination LGPFR (c)
Model I (d)
Model II (e)
Model IIIFigure 4
Prediction results of five GPR models for battery No. 6.
(a)
Basic GPR (b)
Combination LGPFR (c)
Model I (d)
Model II (e)
Model IIIFigure 5
Prediction results of five GPR models for battery No. 7.
(a)
Basic GPR (b)
Combination LGPFR (c)
Model I (d)
Model II (e)
Model IIIThe concrete prediction errors comparisons of the five GPR models for three batteries are indicated in Table1. Here, two criteria of root mean square error (RMSE) and the mean absolute percentage error (MAPE) are introduced to evaluate the prediction performance, which are shown as follows:(18)RMSE=∑i=1myi-y^i2m,MAPE=1m∑i=1myi-y^iyi,where m is the prediction steps, yi represents the actual SOH, and y^i represents the predicted SOH. Table 1 shows the prognostic RMSEs (%) and MAPEs (%) of different SOH prediction methods, respectively. Compared with the other popular published methods in [24, 27], both RMSE and MAPE values have been reduced greatly. It is found that, for three batteries, the prognostics accuracies of SMK-GPR, SE-MGRP, and P-MGPR in [27] are improved evidently compared with the combination LGPFR and combination QGPFR in [24]. It is also seen that, for the proposed Models I, II, and III, the maximum prediction errors of batteries Nos. 5, 6, and 7 are averagely less by 50% than SE-MGRP and P-MGPR models, respectively. Furthermore, the MAPE and RMSE values are less than 1% for three batteries using Model II. In conclusion, it is found that the proposed method shows a much better prediction performance than the other methods. The small values of MAPE and RMSE indicate that the proposed approach can be suitable to perform accurate SOH prediction and satisfy the engineering requirement.Table 1
Prediction errors comparison of different methods for batteries Nos. 5, 6, and 7.
(a) Battery number Error criteria Basic GPR LGPFR QGPFR Combination LGPFR Combination QGPFR SMK-GPR 5 MAPE (%) 12.10 23.0 1.90 1.60 2.10 1.65 RMSE (%) 13.03 1.71 1.50 1.36 1.80 1.38 6 MAPE (%) 27.0 10.30 7.70 10.20 29.0 10.60 RMSE (%) 22.51 6.90 5.12 6.86 20.44 7.08 7 MAPE (%) 19.20 1.90 5.40 1.70 2.60 1.91 RMSE (%) 20.70 1.59 5.52 1.73 2.69 1.88  (b) Battery number Error criteria P-MGPR SE-MGPR Model I Model II Model III 5 MAPE (%) 1.55 1.38 0.80 0.94 0.91 RMSE (%) 1.36 1.20 0.74 0.83 0.83 6 MAPE (%) 2.96 2.93 1.00 0.99 2.18 RMSE (%) 2.12 2.11 0.82 0.81 1.71 7 MAPE (%) 1.09 1.02 0.74 0.73 1.12 RMSE (%) 1.14 1.07 0.82 0.81 1.71 Note. Some experimental results are obtained from [24, 27]; LGPFR: Gaussian process functional regression with linear mean function [24]; QGPFR: Gaussian process functional regression with quadratic polynomial mean function [24]; Combination LGPFR: LGPFR with combination of squared exponential covariance function and periodic covariance function [24]; Combination QGPFR: QGPFR with combination of squared exponential covariance function and periodic covariance function [24]; SMK-GPR: the GPR method with spectral mixture kernels [27]; SE-MGPR: multiscale GPR methods with squared exponential function [27]; P-MGPR: multiscale GPR methods with periodic covariance function [27].The prediction results for No. 6 and No. 7 are similar to No. 5’s. The prediction results of five GPR models are shown in Figures4 and 5.From the above three experiments for different Li-ion batteries, the efficiency of the proposed method is verified.


## 4.1. Raw Data Description of Li-Ion Batteries

To comprehensively illustrate the efficiency of the proposed method, the dataset about the aging of Li-ion batteries is selected which is obtained from the data repository of the NASA Ames Prognostics Center of Excellence (PCoE) [31].The dataset of the NASA Ames PCoE consists of 36 Li-ion batteries’ data: No. 5, No. 6, No. 7, No. 18, and Nos. 25–56. Because the lengths of SOHs extracted from the 36 batteries are different, Li-ion batteries of No. 5, No. 6, and No. 7 are the most commonly used in related literatures [24, 25, 27] because they are much longer. They were run through 3 different operational profiles (charge, discharge, and impedance) at room temperature. Charging was carried out in a constant current (CC) mode at 1.5 A until the battery voltage reached 4.2 V and then continued in a constant voltage (CV) mode until the charge current dropped to 20 mA. Discharge was carried out at a constant current (CC) level of 2 A until the battery voltage fell to 2.7 V, 2.5 V, 2.2 V, and 2.5 V for batteries Nos. 5, 6, 7, and 18, respectively. Impedance measurement was carried out through an electrochemical impedance spectroscopy (EIS) frequency sweep from 0.1 Hz to 5 kHz. Repeated charge and discharge cycles result in accelerated aging of the batteries while impedance measurements provide insight into the internal battery parameters that change as aging progresses. The experiments were stopped when the batteries reached end-of-life (EOL) criteria, which was a 30% fade in rated capacity (from 2 Ahr to 1.4 Ahr). This dataset can be used for the prediction of SOC, SOH, and the remaining useful life (RUL) of batteries [23].In this paper, the data from batteries No. 5, No. 6, and No. 7 is processed using (4) to implement the experiments including training and testing. Figure 2 shows the SOH curves of these three batteries. The total charge/discharge cycles are all 168. Figure 2 indicates that the SOH shows an obvious global degradation trend and local regeneration phenomenon.Figure 2
SOHs of Li-ion batteries No. 5, No. 6, and No. 7.


## 4.2. Training with GPR

For No. 5, No. 6, and No. 7, the first 100 SOH points are utilized to train the proposed GPR model and the remaining 68 points are used for prediction. For the training data from the above three batteries, the mean function is chosen as the linear functionm(x)=ax+b and the covariance function is chosen as three functions discussed above; according to these three models, the algorithm for the training and prediction based on GPR is produced as follows:(a)
Determine the training datasetsi,xii=1N, where i is the number of charge/discharge cycles and xi is the corresponding value of SOH at the ith charge/discharge cycle. For batteries No. 5, No. 6, and No. 7, N equals 100.(b)
Initialize the hyperparameters included in the mean function and covariance function for the different training sets. Here, the initial hyperparameters for Model I areΘ1=a,b,l,sfT=0.1,0,rand,2T, for Model II are Θ2=a,b,l,sf1,ell,sf2T=0.1,0,rand,2,0.9,2T, and for Model III are Θ3=a,b,l,sf1,ell,p,sf2T=0.5,1,rand,2,0.9,2,2T, wherein signal “rand” which is the value of parameter “l” in these three models is a random number between 0 and 1.To compare with other methods, the hyperparameter in basic GPR is selected asΘ′=l,sfT=1,1T and in combination LGPFR is Θ′′=a,b,sf1,l,sf2,ell,ωT=0,0,0.1,1,0.2,1,5T which is the same as the hyperparameters used in [24].(c)
Optimize the hyperparameters with the maximization of the log-likelihood function.(d)
Set the prediction stepsm, for No. 5, No. 6, and No. 7, m=68.(e)
Apply the training datai,xii=1N and the testing data i,xii=N+1N+m to Models I, II, and III with the optimal hyperparameters. Then the prediction outputs will be derived.


## 4.3. Prediction and Comparison

In this part, the prediction results of five GPR models, basic GPR, combination LGPFR, and Models I, II, and III for batteries No. 5, No. 6, and No. 7, are shown in Figures3–5, respectively. As shown in Figure 3, the red lines are the actual SOH value, the blue scatter plots are the prediction result, and the grey regions represent the 95% confidence intervals. Taking Figure 3, for example, for the method of basic GPR shown in Figure 3(a), the mean prediction output is further away from the actual SOH for the subsequent cycles. Meanwhile, the 95% confidence intervals increase significantly. The above two aspects demonstrate the poor prediction effect of basic GPR. For the method of combination LGPFR shown in Figure 3(b), although the prediction values are close to the actual SOH data points, the 95% confidence intervals are still wide. The wide confidence intervals indicate the high uncertainty of the prediction result. Experiments have shown the limitations of Li-ion batteries’ SOH prediction based on the basic GPR model and the combination LGPFR model. The simulations with Models I, II, and III for No. 5 are shown in Figures 3(c)–3(e); it is clear that both the point prediction results and the uncertainty representations are improved greatly.Figure 3
Prediction results of five GPR models for battery No. 5.
(a)
Basic GPR (b)
Combination LGPFR (c)
Model I (d)
Model II (e)
Model IIIFigure 4
Prediction results of five GPR models for battery No. 6.
(a)
Basic GPR (b)
Combination LGPFR (c)
Model I (d)
Model II (e)
Model IIIFigure 5
Prediction results of five GPR models for battery No. 7.
(a)
Basic GPR (b)
Combination LGPFR (c)
Model I (d)
Model II (e)
Model IIIThe concrete prediction errors comparisons of the five GPR models for three batteries are indicated in Table1. Here, two criteria of root mean square error (RMSE) and the mean absolute percentage error (MAPE) are introduced to evaluate the prediction performance, which are shown as follows:(18)RMSE=∑i=1myi-y^i2m,MAPE=1m∑i=1myi-y^iyi,where m is the prediction steps, yi represents the actual SOH, and y^i represents the predicted SOH. Table 1 shows the prognostic RMSEs (%) and MAPEs (%) of different SOH prediction methods, respectively. Compared with the other popular published methods in [24, 27], both RMSE and MAPE values have been reduced greatly. It is found that, for three batteries, the prognostics accuracies of SMK-GPR, SE-MGRP, and P-MGPR in [27] are improved evidently compared with the combination LGPFR and combination QGPFR in [24]. It is also seen that, for the proposed Models I, II, and III, the maximum prediction errors of batteries Nos. 5, 6, and 7 are averagely less by 50% than SE-MGRP and P-MGPR models, respectively. Furthermore, the MAPE and RMSE values are less than 1% for three batteries using Model II. In conclusion, it is found that the proposed method shows a much better prediction performance than the other methods. The small values of MAPE and RMSE indicate that the proposed approach can be suitable to perform accurate SOH prediction and satisfy the engineering requirement.Table 1
Prediction errors comparison of different methods for batteries Nos. 5, 6, and 7.
(a) Battery number Error criteria Basic GPR LGPFR QGPFR Combination LGPFR Combination QGPFR SMK-GPR 5 MAPE (%) 12.10 23.0 1.90 1.60 2.10 1.65 RMSE (%) 13.03 1.71 1.50 1.36 1.80 1.38 6 MAPE (%) 27.0 10.30 7.70 10.20 29.0 10.60 RMSE (%) 22.51 6.90 5.12 6.86 20.44 7.08 7 MAPE (%) 19.20 1.90 5.40 1.70 2.60 1.91 RMSE (%) 20.70 1.59 5.52 1.73 2.69 1.88  (b) Battery number Error criteria P-MGPR SE-MGPR Model I Model II Model III 5 MAPE (%) 1.55 1.38 0.80 0.94 0.91 RMSE (%) 1.36 1.20 0.74 0.83 0.83 6 MAPE (%) 2.96 2.93 1.00 0.99 2.18 RMSE (%) 2.12 2.11 0.82 0.81 1.71 7 MAPE (%) 1.09 1.02 0.74 0.73 1.12 RMSE (%) 1.14 1.07 0.82 0.81 1.71 Note. Some experimental results are obtained from [24, 27]; LGPFR: Gaussian process functional regression with linear mean function [24]; QGPFR: Gaussian process functional regression with quadratic polynomial mean function [24]; Combination LGPFR: LGPFR with combination of squared exponential covariance function and periodic covariance function [24]; Combination QGPFR: QGPFR with combination of squared exponential covariance function and periodic covariance function [24]; SMK-GPR: the GPR method with spectral mixture kernels [27]; SE-MGPR: multiscale GPR methods with squared exponential function [27]; P-MGPR: multiscale GPR methods with periodic covariance function [27].The prediction results for No. 6 and No. 7 are similar to No. 5’s. The prediction results of five GPR models are shown in Figures4 and 5.From the above three experiments for different Li-ion batteries, the efficiency of the proposed method is verified.


## 5. Conclusions

In this paper, a battery SOH estimation approach based on GPRNN is presented. In order to improve the long-term prediction accuracy and uncertainty based on basic GPR and combination GPFR and multiscale GPR, three GPR models based on neural network are used. Compared with the other mixture methods of GPR and multiscale GPR algorithms, the proposed GPRNN models are simpler; therefore, they satisfy the need of real-time prognostics. Moreover, the experimental simulations show that the proposed models have higher prediction accuracy and lower prediction uncertainty for Li-ion batteries with the complicated phenomenon. To further verify the efficiency of the proposed algorithm, using dataset with greater number of degradation data samples in real industrial applications is necessary. Finally, exploring the methodology with better prognostic performances will be the future work.



---

*Source: 8358025-2018-04-01.xml*