A Rest Time-Based Prognostic Framework for State of Health Estimation of Lithium-Ion Batteries with Regeneration Phenomena

State of health (SOH) prognostics is significant for safe and reliable usage of lithium-ion batteries. To accurately predict regeneration phenomena and improve long-term prediction performance of battery SOH, this paper proposes a rest time-based prognostic framework (RTPF) in which the beginning time interval of two adjacent cycles is adopted to reflect the rest time. In this framework, SOH values of regeneration cycles, the number of cycles in regeneration regions and global degradation trends are extracted from raw SOH time series and predicted respectively, and then the three sets of prediction results are integrated to calculate the final overall SOH prediction values. Regeneration phenomena can be found by support vector machine and hyperplane shift (SVM-HS) model by detecting long beginning time intervals. Gaussian process (GP) model is utilized to predict the global degradation trend, and nonlinear models are utilized to predict the regeneration amplitude and the cycle number of each regeneration region. The proposed framework is validated through experimental data from the degradation tests of lithium-ion batteries. The results demonstrate that both the global degradation trend and the regeneration phenomena of the testing batteries can be well predicted. Moreover, compared with the published methods, more accurate SOH prediction results can be obtained under this framework.


Introduction
Lithium-ion batteries have attracted significant attention for many more energy or power demanding applications such as in consumer electronics, aerospace, hybrid and electric vehicles and smart grids applications, due to their superior energy density and high specific energy over other rechargeable battery technologies [1].To reduce the maintenance costs and the chances of failure as well as improve the safety and reliability of systems, battery degradation, prognostics and maintenance optimization have attracted much attention among researchers in the fields of energy, reliability engineering and aerospace engineering [2][3][4].Also, state of health (SOH) prognostics of lithium-ion battery has become a key factor that helps fulfill intelligent battery management [5][6][7].
Traditionally, prognostics can be implemented by using either physics-based approaches or data-driven approaches [8].However, as capacity fade phenomenon is the result of various processes and their interactions, such as chemical side reactions or loss of conductivity [9][10][11], it is difficult to monitor the internal state of a battery in real time and obtain an accurate physics-based model.By Energies 2016, 9, 896 2 of 18 contrast, data-driven approaches can achieve health state assessment and prognostics based only on monitoring parameters and testing data samples.As a result, data-driven techniques draw more and more attention in SOH prognostics.In the literature, statistical, computational and artificial intelligence algorithms, such as autoregressive model [12], particle filter (PF) [13,14], Gaussian process regression [15], Wiener process [16], relevance vector machine (RVM) [17], Bayesian approach [18], support vector machine (SVM) [19] and neural networks [20,21] have been used for battery SOH and remaining useful life (RUL) prognostics in various applications.Capacity fade and impedance increase are the two most used health indicators of batteries.Because capacity is available for monitoring online with low costs, it is widely used in battery health management systems.However, the existence of capacity regeneration phenomena remains a big challenge for prognostics.
Researchers have recognized that reaction products in batteries have a chance to dissipate during the rest time and this phenomenon leads to charge regeneration (or self-recharge).Some researchers have approximately described the law of capacity regeneration by the empirical model of charge regeneration [16,22,23].However, the mechanisms of charge regeneration and capacity regeneration may be quite different.The self-recharge phenomenon just needs about 10 min to become steady, while capacity regeneration takes around 2 h of rest-time to appear [24].Therefore, the short relaxation of 10 min will not lead to the increase of battery capacity.Some researchers have assumed that capacity regeneration has certain regularity over time [24].Liu et al. [15] has predicted SOH regeneration and degradation by using combination Gaussian progress function regression.He et al. [25] has decoupled local regeneration and global degradation with wavelet and predicted them respectively.However, the capacity regeneration phenomenon is strongly related to the rest time which is determined by practical application demands of the battery [24].Therefore, without considering the rest time of the battery, it is hard to predict regeneration phenomenon accurately.Some researchers have treated the capacity regeneration as uncertainty (unpredictable disturbance).Olivares et al. [26] and Orchard et al. [27] have provided some approaches to detect the regeneration phenomenon and predicted SOH by isolating the effects of regeneration.Qin et al. [28] has predicted the global degradation trend directly using the raw SOH data by improving the robustness of a prediction algorithm.Li et al. [29] has utilized the history data of the same type of batteries and indirectly considered the influence of regeneration on battery aging.Although these approaches can capture and predict the global degradation trend of batteries, they have failed to predict the regeneration phenomenon.Therefore, the SOH prediction errors around these related cycles are large.In recent years, although it has been noted that there is a close relationship between the regeneration phenomenon and the rest time [16,24,30], no effective prediction approach that uses the rest time to take regeneration into account has been found.
To address the regeneration phenomenon prediction problem, this paper proposes a rest time-based prognostic framework (RTPF).This framework is composed of extraction phase, prediction phase, and integration phase.In extraction phase, SOH values of regeneration cycles and the number of cycles in regeneration regions and global degradation series are extracted from raw SOH data.After that, the above three series are predicted respectively in prediction phase.In integration phase, the three sets of prediction results are integrated to calculate the final SOH prediction values.Because there is an almost linear relationship between the cycle beginning time interval of two adjacent cycles and the pause time (rest time) during this period, the support vector machine and hyperplane shift (SVM-HS) model is proposed to detect long beginning time intervals to find obvious regeneration phenomena.Also, the Gaussian process (GP) model is used for the prediction of the global degradation trend.Nonlinear models are adopted to predict regeneration amplitude and the number of cycles in regeneration regions, respectively.Finally, the effectiveness of the proposed prognostic framework is validated by using the National Aeronautics and Space Administration (NASA) battery data sets.The results show that the future regeneration phenomena can be well predicted and compared with eight published methods, and more accurate SOH prediction results can be obtained by this framework.
The remainder of this paper is organized as follows: in Section 2, the related classification model and prediction model are introduced.The proposed RTPF is described in detail in Section 3. In Section 4, Energies 2016, 9, 896 3 of 18 the effectiveness of the proposed framework is demonstrated via battery SOH prediction experiments.Finally, the conclusion is provided in Section 5.

Support Vector Machine
Let S = {(x i , y i )|i = 1, 2, ..., N} be a two-class dataset, where x i ∈ R m are the data points and y i ∈ {1, −1} are the corresponding labels.These samples can be used to determine a decision function to separate two classes, which can be expressed by a linear boundary: where w (w ∈ R m ) and b (b ∈ R) are decision coefficients obtained by learning from the dataset S.And x * is a new sample that need to be classified.Hard margin SVM attempts to find the max-margin hyperplane w T x + b = 0.In this way, the points in the positive class have a positive margin and points in the negative class have a negative margin.The optimization problem can be expressed by: min subject to: Introducing Lagrange multipliers α i corresponding to the constraints in Equation (2b), this optimization problem can be solved by using its dual problem counterpart [31].
subject to: Because Equation ( 3) is a standard quadratic programming, the optimization problem can be easily solved, and then the decision coefficients can be obtained as follows: where (x s , y s ) is any support vector which is defined by the condition α i > 0.
For the cases that two classes in the dataset are not linearly separable, the constraints are relaxed by introducing a slack variable ξ i for each data point.Then, the soft margin SVM optimization problem is as: Energies 2016, 9, 896 4 of 18 subject to: where C is the penalty parameter, and ξ i is the slack variable.The goal of large margin can be gotten by minimizing the term ||w|| 2 in the objective function, while the goal of small margin violations can be obtained by minimizing the term ∑ N i=1 ξ i .The parameter C is defined by users' large-margin favor or small-violations favor.Lagrange multipliers α i are introduced for Equation (6b) and β i are for Equation (6c).Then β i can be removed by KKT condition β i = C − α i .The final dual expression of soft margin SVM is as follows: subject to: Decision coefficients can also be obtained by Equations ( 4) and (5), where (x s , y s ) is any support vector which is defined by the condition 0 < α i < C.Moreover, it can be seen from Equations ( 3) and (7) that the only difference between the dual problem of hard margin SVM and the dual problem of soft margin SVM is the constraint C on α i .This means that if the dataset is linearly separable, hard margin SVM can be obtained by setting a large penalty parameter C in the soft margin SVM.
For any test vector x * , the decision function of both hard margin SVM and soft margin SVM is finally given by: The Vapnik-Chervonenkis (VC) dimension is usually used to measure the complexity of a model.The VC dimension of SVM is bounded as follows [32]: where R is the smallest diameter which encloses all training sample, d is the dimension of given training samples and M is the margin.For general d dimension classification models, the VC dimension is (d + 1).For example, the VC dimension of perceptron learning algorithm (PLA), one of the simplest classifiers in machine learning, is shown as follows [33]: An advantage of SVM-HS method is that the complexity can be bounded by R 2 /M 2 and does not explicitly depend on the dimension d.

Hyperplane Shift
When the dataset is unbalanced, the majority class is likely to have more representatives.The noise of majority class may cause the SVM to move the separating hyperplane towards the minority class.The described effects are illustrated in Figure 1.One approach to solving this problem is to shift the hyperplane.SVM can be trained in the traditional way, but the decision function given by Equation ( 8) is replaced by: where ( 1,1) p ∈ − is the offset.The shift range is within the positive and negative margin bounds.
The hyperplane is shifted towards positive data points when 0 p < , and towards negative data points when 0 p > .Also, the shifted hyperplane gives a larger margin for the minority class, and a smaller margin for the majority class.Moreover, the sum of positive and negative margin is unchanged, which means that the large margin principle of SVM is still effective.Because the offset parameter p is constant, the complexity of SVM-HS is the same as that of SVM.If the dataset is linearly separable, using hard margin SVM or soft margin SVM with large penalty parameter, the hyperplane shift has no effect on the accuracy of classification of all the data points.If the dataset is not linearly separable, using soft margin SVM with large penalty parameter, the shifted hyperplane is in the area of overlap.In this situation, the classification accuracy may be increased or decreased, because solving the optimization problem (Equation ( 6)) is to minimize the term  One approach to solving this problem is to shift the hyperplane.SVM can be trained in the traditional way, but the decision function given by Equation ( 8) is replaced by: where p ∈ (−1, 1) is the offset.The shift range is within the positive and negative margin bounds.
The hyperplane is shifted towards positive data points when p < 0, and towards negative data points when p > 0. Also, the shifted hyperplane gives a larger margin for the minority class, and a smaller margin for the majority class.Moreover, the sum of positive and negative margin is unchanged, which means that the large margin principle of SVM is still effective.Because the offset parameter p is constant, the complexity of SVM-HS is the same as that of SVM.If the dataset is linearly separable, using hard margin SVM or soft margin SVM with large penalty parameter, the hyperplane shift has no effect on the accuracy of classification of all the data points.If the dataset is not linearly separable, using soft margin SVM with large penalty parameter, the shifted hyperplane is in the area of overlap.In this situation, the classification accuracy may be increased or decreased, because solving the optimization problem (Equation ( 6)) is to minimize the term One approach to solving this problem is to shift the hyperplane.SVM can be trained in the traditional way, but the decision function given by Equation ( 8) is replaced by: where is the offset.The shift range is within the positive and negative margin bounds.
The hyperplane is shifted towards positive data points when 0 p < , and towards negative data points when 0 p > .Also, the shifted hyperplane gives a larger margin for the minority class, and a smaller margin for the majority class.Moreover, the sum of positive and negative margin is unchanged, which means that the large margin principle of SVM is still effective.Because the offset parameter p is constant, the complexity of SVM-HS is the same as that of SVM.If the dataset is linearly separable, using hard margin SVM or soft margin SVM with large penalty parameter, the hyperplane shift has no effect on the accuracy of classification of all the data points.If the dataset is not linearly separable, using soft margin SVM with large penalty parameter, the shifted hyperplane is in the area of overlap.In this situation, the classification accuracy may be increased or decreased, because solving the optimization problem (Equation ( 6)) is to minimize the term

Gaussian Process Model
GP model is a flexible nonparametric model, which has been widely applied to multi-step-ahead predictions in time series analysis [25,34].A GP model is completely specified by its mean function and covariance function.The mean function m (x) and the covariance function k (x, x ) of a real process f (x) are as: The regression function is denoted by: It is common to utilize a mean function of m (x) = 0, since the GP is flexible enough to model the mean arbitrarily well [35].Also, the squared exponential covariance function applied in SOH prognostics is as follows: where σ 2 f controls the vertical scale of this function, and l is a length-scale parameter that governs the speed of the correlation decrease as the input data distance increases.
For training dataset D = {(x i , y i )|i = 1, 2, ..., N}, the target y i is described as: where For a new test input x * , the prior distribution for GP is as follows: where k (x * , X) = [k (x * , x 1 ) , k (x * , x 2 ) , ..., k (x * , x N )].The posterior distribution for test input x * is: where:

The Regeneration Phenomenon in Battery Degradation
In this paper, capacity is adopted to indicate the SOH of lithium-ion batteries, and it can be calculated by integrating current over time from full charge to full discharge.The relationship between SOH and capacity is described as: where C i is the capacity value of cycle i and C 0 is the initial capacity value with no degradation.Many articles have mentioned that capacity regeneration can improve the health state of the battery and slow down battery aging, and it has been confirmed by controlled experiments by Eddahech et al. [24].However, the mechanism of capacity regeneration is still a much-debated topic.Agubra and Fergus [36] suggested that the loss or consumption of recyclable lithium ions due to the Energies 2016, 9, 896 7 of 18 growth of the solid-electrolyte interface (SEI) layer is a main cause of the reduction in the reversible capacity of the lithium ion battery.During the rest time, the acidic impurities generated at the positive electrode can destroy the SEI layer at the negative and liberate lithium ions [37,38].Eddahech et al. [39] pointed out that recovery phenomenon may be related to charge redistribution when there is no charge or discharge force acting on them.Baghdadi et al. [40] argued that a rearrangement of the lithium within the active material crystalline structure occurs thanks to solid diffusion during the battery rest.Thus, less disorder and more accessibility to lithium intercalation sites promotes the battery capacity, allowing more ampere-hours to be stored.
Figure 3 shows the battery SOH degradation data (Battery No. 5 from NASA Ames Prognostics Center of Excellence (PCoE) Center, Moffett Field, CA, USA). Figure 3a illustrates the curve of SOH and calendar time.The amplitude of SOH has an increase when the beginning time interval between two cycles is relatively long.It demonstrates that the regeneration phenomenon is obvious when the rest time is long enough.Figure 3b shows the relationship between SOH and cycle number.Moreover, we can see that after obvious regeneration appearance, SOH degradation rate at the next several cycles is much faster than before.Therefore, it is necessary to analyze the effects of regeneration and global degradation separately.We call the cycle with a significant capacity increase caused by a long rest time as "regeneration cycle" and the previous one cycle as "the cycle before regeneration".Compared with the cycle before regeneration, the capacity values of the next several cycles increase.We call these cycles as "regeneration region", the number of cycles in each regeneration region as "regeneration cycle number" and the cycle except regeneration regions as "global cycle".These definitions will be directly used below and they are illustrated in Figure 3b.
Energies 2016, 9, 896 7 of 18 Many articles have mentioned that capacity regeneration can improve the health state of the battery and slow down battery aging, and it has been confirmed by controlled experiments by Eddahech et al. [24].However, the mechanism of capacity regeneration is still a much-debated topic.Agubra and Fergus [36] suggested that the loss or consumption of recyclable lithium ions due to the growth of the solid-electrolyte interface (SEI) layer is a main cause of the reduction in the reversible capacity of the lithium ion battery.During the rest time, the acidic impurities generated at the positive electrode can destroy the SEI layer at the negative and liberate lithium ions [37,38].Eddahech et al. [39] pointed out that recovery phenomenon may be related to charge redistribution when there is no charge or discharge force acting on them.Baghdadi et al. [40] argued that a rearrangement of the lithium within the active material crystalline structure occurs thanks to solid diffusion during the battery rest.Thus, less disorder and more accessibility to lithium intercalation sites promotes the battery capacity, allowing more ampere-hours to be stored.
Figure 3 shows the battery SOH degradation data (Battery No. 5 from NASA Ames Prognostics Center of Excellence (PCoE) Center, Moffett Field, CA, USA). Figure 3a illustrates the curve of SOH and calendar time.The amplitude of SOH has an increase when the beginning time interval between two cycles is relatively long.It demonstrates that the regeneration phenomenon is obvious when the rest time is long enough.Figure 3b shows the relationship between SOH and cycle number.Moreover, we can see that after obvious regeneration appearance, SOH degradation rate at the next several cycles is much faster than before.Therefore, it is necessary to analyze the effects of regeneration and global degradation separately.We call the cycle with a significant capacity increase caused by a long rest time as "regeneration cycle" and the previous one cycle as "the cycle before regeneration".Compared with the cycle before regeneration, the capacity values of the next several cycles increase.We call these cycles as "regeneration region", the number of cycles in each regeneration region as "regeneration cycle number" and the cycle except regeneration regions as "global cycle".These definitions will be directly used below and they are illustrated in Figure 3b.

Rest Time-Based Prognostics Strategy
Our motivation is to decouple global degradation trend and regeneration, predict them respectively and then integrate them to get the overall prediction results.The proposed framework is outlined in Figure 4.This framework is composed of extraction phase, prediction phase and integration phase.
In extraction phase, it is assumed that the battery is fully charged in the initial state, and the cycle starts from discharge.
( ) T k is a series of the calendar time of the beginning of each cycle.
( ) H k is a time series of SOH from cycle 1 to cycle n, and it can be obtained by capacity

Rest Time-Based Prognostics Strategy
Our motivation is to decouple global degradation trend and regeneration, predict them respectively and then integrate them to get the overall prediction results.The proposed framework is outlined in Figure 4.This framework is composed of extraction phase, prediction phase and integration phase.
In extraction phase, it is assumed that the battery is fully charged in the initial state, and the cycle starts from discharge.{T E (k)} n k=1 is a series of the calendar time of the beginning of each cycle.{H E (k)} n k=1 is a time series of SOH from cycle 1 to cycle n, and it can be obtained by capacity monitoring during the usage of the battery.Also, cycle number series is denoted as {c E (k)} n k=1 = {1, 2, ..., n}.  ( ) DT k is expressed as: ( ) ( ) ( ) Moreover, ( ) DT k can also be described by: where tD(k) is the total discharge time of cycle k.tC(k) is the total charge time of cycle k.tP(k) is the total pause time (rest time) of cycle k.Although the battery may pause several times in a cycle, the lasting time of battery charge and discharge always has little change under certain working conditions.As a result, from Equation ( 23), the beginning time interval DTE(k) is almost linear with the sum of pause time of cycle k.Moreover, compared with the rest time, the beginning time interval calculated by the recorded beginning time of each cycle is much easier to be obtained.Therefore, in this paper, the beginning time interval instead of the rest time is adopted as the input of this framework.{DT E (k)} n−1 k=1 is the forward difference series of the beginning time of each cycle.DT E (k) is expressed as: Moreover, DT E (k) can also be described by: where t D (k) is the total discharge time of cycle k.t C (k) is the total charge time of cycle k.t P (k) is the total pause time (rest time) of cycle k.Although the battery may pause several times in a cycle, the lasting time of battery charge and discharge always has little change under certain working conditions.As a result, from Equation ( 23), the beginning time interval DT E (k) is almost linear with the sum of pause time of cycle k.Moreover, compared with the rest time, the beginning time interval calculated by the recorded beginning time of each cycle is much easier to be obtained.Therefore, in this paper, the beginning time interval instead of the rest time is adopted as the input of this framework.{DH E (k)} n−1 k=1 is the forward difference series of SOH, where DH E (k) is the SOH increment between cycle k and k + 1 and is shown as: Energies 2016, 9, 896 9 of 18 If the SOH of cycle (k + 1) is much better than that of cycle k, there may be significant regeneration, measurement errors, unexpected fluctuations and other factors during this period.{y E (k)} n−1 k=1 is to distinguish whether obvious SOH regeneration has appeared in each cycle.
where y E (k) = −1 means that the value of capacity has a large increase from cycle k to cycle (k + 1), while y E (k) = 1 means the opposite.Th is the given threshold to distinguish between the two states.Generally, regeneration cycles account for a very small number of the total cycles, which makes the dataset {(DT E (k) , y E (k))|k = 1, 2, ..., n − 1} unbalanced.Moreover, in the SOH degradation data, there is noise besides regeneration, and they are linearly non-separable.In this situation, some common methods, such as PLA, are not applicable.That is because PLA can only solve linearly separable problems.However, SVM-HS model can well balance the accuracy and generalization with the linearly non-separable samples.The hyperplane shift parameter p can be determined according to the unbalanced degree of the two-class data.Significantly, in this case, d = 1 and R 2 /M 2 ≥ 1, so the VC dimension of SVM-HS equals to 2 and it is the same as that of PLA.This mean that the complexity of SVM-HS and that of PLA is the same.
Then, the trained SVM-HS model is used to update {y E (k)} n−1 k=1 .This is to ensure that the capacity value increase is indeed caused by regeneration.Put the cycle before regeneration into set C E which can be described by: It is assumed that C E contains q elements.It means that obvious regeneration of the battery has appeared q times.We sort the elements of C E in ascending order and then get a new series {c E (k)} q k=1 .The time interval between the beginning of cycle c E (k) and cycle (c E (k) + 1) is expressed as: The regeneration capacity amplitude from cycle c E (k) to cycle (c E (k) + 1) is described by: The algorithm to extract regeneration region sets and regeneration cycle number series are as Algorithm 1.
Algorithm 1. Extraction of regeneration region sets and regeneration cycle number series.

1:
for i = 1 to q do 2: set j = 1 3: put cycle number c E (i) + j into set S i 5: j = j + 1 6: end while 7: end for 8: for i = 1 to q − 1 do 9: for j = i + 1 to q do 10: end for 12: NU M E (i) = |S i | 13: end for 14: NU M E (q) = S q Energies 2016, 9, 896 10 of 18 All the cycle numbers of global degradation are denoted by set C E , which can be expressed as: We sort the g elements of C E into ascending order and get a new series {c E (k)} g k=1 .With series {c E (k)} g k=1 , the global degradation of SOH can be described by: In prediction phase, {T p (k)} m k=1 consists of the beginning time of cycle (n + 1) to cycle (n + m).It is used for SOH prediction and can be acquired from the working plan or history information of used relevant batteries.Cycle number series is c p (k) is the forward difference series of the beginning time of each cycle, where DT p (k) is as: Using series {DT p (k)} m−1 k=0 as the input, the trained SVM-HS model is adopted to predict the value of {y p (k)} n−1 k=0 .Put the cycles before regeneration into set C p which can be described by: With the assumption that obvious regeneration phenomena have happened p times, there are p elements in C p .Sort them into ascending order and then get the series {c p (k)} p k=1 .Long time intervals that lead to obvious regeneration compose the series {DT p (k)} p k=1 , where DT p (k) is as: Based on the regeneration amplitude prediction algorithm, the series {DH p (k)} p k=1 can be predicted.Also, based on the regeneration cycle number prediction algorithm, the series {NUM p (k)} .The capacity regeneration phenomenon is strongly related to the rest time but it is not affected by aging, because the phenomenon is reproduced in the same way even when the battery becomes aged [24].The relationship between rest time and the amplitude of capacity regeneration is a nonlinear correlation.The battery can be at rest for a very long time, but capacity regeneration amplitude is limited.Therefore, as shown in Figure 5, we assume that the increase rate of capacity regeneration amplitude declines as the rest time increases.
The capacity regeneration phenomenon is strongly related to the rest time but it is not affected by aging, because the phenomenon is reproduced in the same way even when the battery becomes aged [24].The relationship between rest time and the amplitude of capacity regeneration is a nonlinear correlation.The battery can be at rest for a very long time, but capacity regeneration amplitude is limited.Therefore, as shown in Figure 5, we assume that the increase rate of capacity regeneration amplitude declines as the rest time increases.The mapping function adopted in this article is a hyperbolic tangent function which is widely used in machine learning.It is described as: The mapping function adopted in this article is a hyperbolic tangent function which is widely used in machine learning.It is described as: where c is a length-scale parameter.The relationship between regeneration amplitude and DT (k) can be expressed by: Therefore, the elements of {DH p (k)} p k=1 can be predicted by: Similarly, we assume that the increase rate of regeneration cycle number declines with the increase of rest time and take hyperbolic tangent function as the mapping function.Thus, the regeneration cycle number series {NUM p (k)} p k=1 can be predicted by: In this paper, Gaussian process model recommended in Section 2.2 is used in global degradation prediction.
In the integration phase, the predicted SOH values of each regeneration region are added into the global regeneration series.Then, the SOH prediction values of all cycles are obtained.The SOH series of the ith regeneration region is as {RH

k=1
. Based on the assumption that SOH linearly degrades in the regeneration regions, the value of element RH (i) p (k) can be calculated by: Series {H (i) is updated by: For different sizes and chemistries of lithium-ion batteries, the capacity regeneration phenomenon may or may not appear.If it appears, the proposed framework can be used for decoupling the global degradation and local regeneration; and if capacity regeneration phenomenon doesn't appear, the value of discriminate function (Equation ( 25)) always equals to 1 and the proposed framework only contains prediction method of global degradation.In other words, {H (p) Therefore, the proposed framework is applicable for different batteries with or without the capacity regeneration phenomenon.

Battery Data Set
The lithium-ion battery capacity degradation data adopted to conduct the experiments of SOH prediction are obtained from the data repository of the NASA Ames PCoE [41].A set of commercial 18,650 lithium-ion batteries were run at room temperature through three operational profiles (charge, Energies 2016, 9, 896 12 of 18 discharge, and impedance).Charging was carried out at a constant level of 1.5 A until the battery voltage is up to 4.2 V and then it was continued with the voltage of 4.2 V until the charge current fell to 20 mA.Discharge was carried out in a constant mode at 2 A until the battery voltage dropped to 2.7, 2.5, and 2.2 V for batteries Nos. 5, 6, and 7, respectively.Impedance was measured by an electrochemical impedance spectroscopy (EIS) with frequency from 0.1 to 5 kHz.Repeated charge and discharge of the batteries resulted in their accelerated aging.The experiments were stopped when the batteries reached end-of-life criteria, a 30% fade in rated capacity (from 2 to 1.4 Ah).
Based on Equation ( 21), SOH can be easily derived from the experimental capacity data.

Prediction and Comparison
The SOH values and the beginning time of the first 100 cycles are adopted as the inputs of the method given in Section 3 to extract the regeneration information.In this paper, the threshold is set at 0.2% (Th = 0.2%) and the hyperplane shift parameter is set at −0.5 (p = −0.5).
The extraction results of regeneration information are illustrated in Figure 7. From the figure, we can see that all the cycles with a remarkable increase of SOH value are extracted from the 100 training samples.It indicates that for these three sets of degradation data, regeneration is the main factor of capacity fluctuations.Moreover, all the extracted cycles before regeneration of the three batteries are the same.This is because, in the experiments, the three batteries have the same beginning time of every cycle.However, for these three batteries, regeneration amplitude and regeneration cycle number of each regeneration region are different.This is because of individual differences and their different working conditions.

Prediction and Comparison
The SOH values and the beginning time of the first 100 cycles are adopted as the inputs of the method given in Section 3 to extract the regeneration information.In this paper, the threshold is set at 0.2% (Th = 0.2%) and the hyperplane shift parameter is set at −0.5 (p = −0.5).
The extraction results of regeneration information are illustrated in Figure 7. From the figure, we can see that all the cycles with a remarkable increase of SOH value are extracted from the 100 training samples.It indicates that for these three sets of degradation data, regeneration is the main factor of capacity fluctuations.Moreover, all the extracted cycles before regeneration of the three batteries are the same.This is because, in the experiments, the three batteries have the same beginning time of every cycle.However, for these three batteries, regeneration amplitude and regeneration cycle number of each regeneration region are different.This is because of individual differences and their different working conditions.
Using the intervals of the beginning time of two adjacent cycles which can be known from battery working plan as the input, the cycle before regeneration series c p (k) are shown in Table 1.
we can see that all the cycles with a remarkable increase of SOH value are extracted from the 100 training samples.It indicates that for these three sets of degradation data, regeneration is the main factor of capacity fluctuations.Moreover, all the extracted cycles before regeneration of the three batteries are the same.This is because, in the experiments, the three batteries have the same beginning time of every cycle.However, for these three batteries, regeneration amplitude and regeneration cycle number of each regeneration region are different.This is because of individual differences and their different working conditions.The regeneration cycle number of batteries Nos. 5, 6, and 7 is predicted based on Equation (37).The results are shown in Table 2.The regeneration amplitude of batteries Nos. 5, 6, and 7 predicted by Equation ( 36) is shown in Table 3.The curves of global SOH degradation of three batteries are illustrated in Figure 8.Compared with the curves in Figure 6, they are smooth and have fewer fluctuations.GP model is used to predict the global degradation of SOH values of the batteries Nos. 5, 6, and 7, and the initial hyper-parameters Θ = l, σ f , σ n are set by experience, as l = g, σ f = max (DH t (k 1 )) − min (DH t (k 2 )) (k 1 = 1, 2, ..., t; k 2 = 1, 2, ..., t), σ n = 1.5max (|DH t (k) − DH t (k + 1)|) (k = 1, 2, ..., t − 1).The particle swarm optimization algorithm is used to determine the hyper-parameters of GP model for the three Moreover, different types of discharge/charge profiles may lead to the different degrees of capacity regeneration.If there are regeneration phenomena during battery degradation, this framework can be employed to decouple global degradation trend and regeneration.In this article, the prediction methods of global degradation, regeneration amplitude and regeneration cycle number are all data-driven approaches which realize prediction (extrapolation) by learning from training samples.Therefore, for different capacity regeneration phenomena, the parameters of these three prediction methods will be adopted from corresponding training samples to suit various degradation trends.Therefore, for different discharge/charge profiles, the parameters of prediction models (global degradation prediction model, regeneration cycle number prediction model, and regeneration amplitude prediction model) would be changed, while the structure of this framework would not be changed.

Conclusions
To obtain more accurate long-term prediction results of battery SOH with capacity regeneration phenomenon, we have proposed a RTPF.It consists of extraction phase, prediction phase and integration phase.Based on this framework, regeneration phenomenon and global degradation are decoupled.The two parts are predicted separately and finally integrated together to obtain the overall SOH degradation results.SVM-HS model is proposed to classify the unbalanced dataset of beginning time intervals of two adjacent cycles.Moreover, GP model is used to predict the global degradation trend of battery SOH, and nonlinear models are utilized to predict the regeneration amplitude and the cycle number of each regeneration region, respectively.
Case studies are carried out to evaluate the performance of the proposed framework.The results show that: (1) in extraction phase, SVM-HS model can extract the beginning time intervals of two adjacent cycles that cause obvious capacity regeneration; in prediction phase, it can successfully predict all the cycles before regeneration; (2) even around the regeneration cycles, this framework also has a high prediction accuracy; and (3) using the previous capacity data of the current battery, compared with the published methods, this model can achieve more accurate SOH prediction results.

Figure 1 .
Figure 1.The effects of noise of majority class: (a) without noise; and (b) with noise.

Figure 2 .
Figure 2. Hyperplane shifted with data points: (a) two classes are linearly separable; and (b) two classes are not linearly separable.

Figure 1 .
Figure 1.The effects of noise of majority class: (a) without noise; and (b) with noise.

Figure 1 .
Figure 1.The effects of noise of majority class: (a) without noise; and (b) with noise.

Figure 2 .
Figure 2. Hyperplane shifted with data points: (a) two classes are linearly separable; and (b) two classes are not linearly separable.Figure 2. Hyperplane shifted with data points: (a) two classes are linearly separable; and (b) two classes are not linearly separable.

Figure 2 .
Figure 2. Hyperplane shifted with data points: (a) two classes are linearly separable; and (b) two classes are not linearly separable.Figure 2. Hyperplane shifted with data points: (a) two classes are linearly separable; and (b) two classes are not linearly separable.

Figure 3 .
Figure 3.The experimental degradation data: (a) with calendar time; and (b) with cycle number.

Figure 3 .
Figure 3.The experimental degradation data: (a) with calendar time; and (b) with cycle number.
usage of the battery.Also, cycle number series is denoted as

Figure 4 .
Figure 4. Framework of rest time-based prognostic strategy.

Figure 4 .
Figure 4. Framework of rest time-based prognostic strategy.
p k=1 can be predicted.Taking global SOH degradation series {H E (k)} g k=1 as the training samples, the m − ∑ p l=1 NU M p (l) steps' SOH prediction values of global degradation can be obtained and written as {H (0) p (k)} m−∑ p l=1 NUM p (l) k=1 steps' SOH prediction values of global degradation can be obtained and written as ( )

Figure 5 .
Figure 5.The assumed relationship between rest time and regeneration amplitude.

Figure 5 .
Figure 5.The assumed relationship between rest time and regeneration amplitude.
with the trained SVM-HS model obtained from the extraction phrase.Although the corresponding series c p (k) p k=1 and DT p (k) p k=1 are determined by different SVM-HS models, all the predicted cycles before regeneration of batteries Nos. 5, 6, and 7 are the same.This is because the rest time of the three batteries is similar.The prediction results of c p (k) p k=1 and DT p (k) p k=1

Table 1 .
Cycle before regeneration and corresponding beginning time intervals.