State of Health Estimation of Li-ion Batteries with Regeneration Phenomena: a Similar Rest Time-based Prognostic Framework

State of health (SOH) prediction in Li-ion batteries plays an important role in intelligent battery management systems (BMS). However, the existence of capacity regeneration phenomena remains a great challenge for accurately predicting the battery SOH. This paper proposes a novel prognostic framework to predict the regeneration phenomena of the current battery using the data of a historical battery. The global degradation trend and regeneration phenomena (characterized by regeneration amplitude and regeneration cycle number) of the current battery are extracted from its raw SOH time series. Moreover, regeneration information of the historical battery derived from corresponding raw SOH data is utilized in this framework. The global degradation trend and regeneration phenomena of the current battery are predicted, and then the prediction results are integrated together to calculate the overall SOH prediction values. Particle swarm optimization (PSO) is employed to obtain an appropriate regeneration threshold for the historical battery. Gaussian process (GP) model is adopted to predict the global degradation trend, and linear models are utilized to predict the regeneration amplitude and the cycle number of each regeneration region. The proposed framework is validated using experimental data from the degradation tests of Li-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
Due to high efficiency and energy density, wide temperature range, low self-discharge rate and relatively long cycle life, Li-ion batteries have been widely used in various applications, such as hybrid and electric vehicles, trams, satellite, cell phones and laptops [1].To ensure the performance and reliability of Li-ion batteries, the battery management system (BMS) has two core tasks: estimate the state-of-charge (SOC) and the state-of-health (SOH) [2,3].SOC is the percentage of the remaining charge to the battery's current maximum capacity [4], while SOH is the ability of a battery to provide its nominal capacity over its service lifetime [5].Thus far, many robust and accurate approaches to estimate SOC have been studied [2,[6][7][8].However, the development of the SOH estimation methods is still much more challenging problem.
The capacity degradation is the result of various processes and their interactions, such as chemical side reactions or loss of conductivity [1,9,10].Therefore, it is difficult to monitor the internal state of Symmetry 2017, 9, 4 2 of 19 a battery in real time and to establish an accurate physics-based model.In recent years, data-driven approaches have become a popular direction because they can achieve health state prognostics based only on monitoring parameters and testing data samples [11].That is to say, it does not need to explicitly model complex electrochemical reaction and related principles.In the literature, many data-driven methods, such as auto-regressive 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) prediction in various applications.In these applications, because capacity is available for online monitoring with low cost, it becomes the most common used health indicator of battery SOH.However, the existence of capacity regeneration phenomena remains a big challenge for SOH prediction [22].
Researchers have indicated that reaction products in batteries may dissipate during the rest time and lead to charge regeneration.Some researchers have approximately described the law of capacity regeneration by the empirical model of charge regeneration [16,23,24].However, the mechanism of charge regeneration and capacity regeneration may be quite different.Because self-recharge phenomenon just needs about 10 min to become steady, while the appearance of obvious capacity regeneration needs around 2 h of rest-time [25].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 [25].Liu [15] has predicted SOH regeneration and degradation based on a combination of two Gaussian progress function regression models.He [26] has decoupled local regeneration and global degradation with wavelet and predicted them respectively.However, capacity regeneration phenomenon is related to the rest time which is determined by practical application demands [25].It is hard to accurately predict the regeneration phenomenon without considering rest time.Some researchers have treated the capacity regeneration as uncertainty (unpredictable disturbance).Olivares [27] and Orchard [28] have provided some approaches to detect the regeneration phenomenon and predicted SOH by isolating the effects of regeneration.Qin [29] has predicted the global degradation trend directly using the raw SOH data by improving the robustness of a prediction algorithm.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, it has been noted that there is a close relationship between the regeneration phenomenon and the rest time [16,25].There are two possible ways of modeling the relationship between them.One way is to add an extra hardware to record the rest time during each charge-discharge cycle and establish the model of regeneration phenomenon and the rest time directly.The advantage of this approach is that it can predict the regeneration phenomenon even when this type of battery is used in the current scenario for the first time; the drawback is that the added extra hardware will increase the size and the cost of BMS.Another way is to indirectly consider the effects of the rest time by using the SOH data of a historical battery whose status (working time and rest time) is similar to the current battery.The advantage of this method is that it does not need a rest-time-recording hardware; the disadvantage is that it requires the degradation data of a historical battery.By comparing their advantages and disadvantages, designers can make a choice from the two approaches based on the actual demand.
The first approach is described in our previous work [30].In this paper, the second approach called a similar rest time-based prognostic framework (SRTPF) is proposed.This framework is composed of extraction phase and prediction phase.In extraction phase, after giving the regeneration threshold value of the current battery, regeneration amplitude series, regeneration cycle number series and global degradation series of the current battery can be extracted from the raw SOH data.Particle swarm optimization (PSO) algorithm is adopted to determine the threshold value of the historical battery.Regeneration amplitude series and regeneration cycle number series of the historical battery are extracted from the raw SOH data of the historical battery.After that, in prediction phase, the above three series of the current battery are predicted and integrated to calculate the final SOH prediction values.In addition, Gaussian process (GP) model is used for the prediction of the global degradation trend.Linear models are adopted to predict regeneration amplitude series and regeneration cycle number series.Finally, the effectiveness of the proposed prognostic framework is validated using the NASA battery datasets.
The remainder of this paper is organized as follows: In Section 2, the related optimization model and prediction model are introduced.The proposed similar rest time-based prognostic framework is described in detail in Section 3. The effectiveness of the proposed framework is demonstrated via battery SOH prediction using NASA datasets in Section 4. Finally, the conclusion is provided in Section 5.

PSO Algorithm
In 1995, PSO algorithm was first proposed by Kennedy and Eberhart [31].The main thought is to find an optimal solution via fitness function(s).The velocities and positions of particles are updated by [12,32]: where v k i and x k i are the velocity and position of particle i at iteration k, x k i denotes the position of particle i at iteration k, pbest k i is the individual optimal solution of particle i after k iterations, gbest k is the global optimal solution after k iterations, w denotes the inertia weight coefficient, r 1 and r 2 are two random numbers uniformly distributed between 0 and 1, and c 1 and c 2 are positive acceleration constants.

Gaussian Process Model
Gaussian process (GP) model is a flexible nonparametric model.It has been widely applied to multi-step-ahead predictions in time series analysis [26,33].A GP model is completely specified by the mean function m(x) and the covariance function k(x, x ), where m(x) and k(x, x ) are described as follows: The GP regression function is expressed as It is common to utilize a mean function of m(x) = 0, since GP is flexible enough to model the mean arbitrarily well [34].In addition, 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 the length-scale parameter and 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 by where ε i ∼ N 0, σ 2 , i = 1, . . ., N. For new input x * , the prior distribution of the GP is described as where k(x * , X) = [k(x * , x 1 ), k(x * , x 2 ), . . . ,k(x * , x N )].The posterior distribution of x * is where 3. The Proposed Framework

The Regeneration Phenomenon
In this paper, capacity is used to indicate the SOH of Li-ion batteries.It can be obtained by integrating the current over time.The relationship between SOH and capacity can be described as where C i is the capacity value of charge-discharge cycle i and C 0 is the initial value of capacity without degradation.
Figure 1 shows the battery SOH degradation data (Battery No. 05 from NASA PCoE Center, Moffett Field, CA, USA). Figure 1a illustrates the change of SOH over calendar time.It is noted that the value of SOH has increased when the interval between two cycles is relatively large.It demonstrates that capacity regeneration is obvious when the battery has a long rest time.Figure 1b presents the relationship between SOH and the cycle number.It is observed that after obvious regeneration appearance, SOH degradation rate of the next several cycles is much faster than before.Therefore, it is meaningful to analyze the effects of capacity regeneration and global degradation separately.We call the cycle with a significant capacity increase as "regeneration cycle" and the previous one cycle as "the cycle before regeneration".Compared with the cycle before regeneration, the SOH values of the next several cycles have increased.In this article, we call these cycles as "regeneration region" and the difference between the maximum SOH value in a regeneration region and the SOH value at the cycle before regeneration as "regeneration amplitude".Moreover, we call the number of cycles in each regeneration region as "regeneration cycle number" and the cycle except regeneration regions as "global cycle".These definitions are directly used below and illustrated in Figure 1b.

Similar Rest Time-Based Prognostics Strategy
It has been noticed that long time rest will lead to battery capacity regeneration [16,25,35].Our aim is to decouple the global degradation and regeneration from battery SOH, predict them, and integrate their results together to get the overall SOH prediction results.In this article, the influences of rest time to capacity regeneration are characterized by regeneration amplitude and regeneration cycle number.As shown in Figure 2,  is the mapping between the rest time and the regeneration cycle number, and  is the mapping between the rest time and the regeneration amplitude.The different working conditions and the individual difference of the current and historical batteries result in the differences of the mappings of the two batteries, which are distinguished by the subscripts c and h, respectively.When the rest time of the current battery is unavailable, the relationship of the regeneration of the current and historical batteries can be associated by mappings  The schematic of similar rest time-based prognostic framework is shown in Figure 3.This framework contains extraction phase and prediction phase.In the extraction phase, SOH series of current battery and SOH series of historical battery can be obtained by capacity monitoring during the usage of the current and historical batteries.The regeneration threshold of the current battery is given by considering the magnitude of noises.Then, the global degradation trend and regeneration phenomena (regeneration amplitude and regeneration cycle number) of the current battery can be decoupled.Based on the similarity between the current and historical batteries, the threshold of the historical battery can be determined automatically by PSO.After that, the regeneration phenomena of the historical battery can be decoupled with this threshold.In the prediction phase, the global

Similar Rest Time-Based Prognostics Strategy
It has been noticed that long time rest will lead to battery capacity regeneration [16,25,35].Our aim is to decouple the global degradation and regeneration from battery SOH, predict them, and integrate their results together to get the overall SOH prediction results.In this article, the influences of rest time to capacity regeneration are characterized by regeneration amplitude and regeneration cycle number.As shown in Figure 2, φ is the mapping between the rest time and the regeneration cycle number, and ϕ is the mapping between the rest time and the regeneration amplitude.The different working conditions and the individual difference of the current and historical batteries result in the differences of the mappings of the two batteries, which are distinguished by the subscripts c and h, respectively.When the rest time of the current battery is unavailable, the relationship of the regeneration of the current and historical batteries can be associated by mappings

Similar Rest Time-Based Prognostics Strategy
It has been noticed that long time rest will lead to battery capacity regeneration [16,25,35].Our aim is to decouple the global degradation and regeneration from battery SOH, predict them, and integrate their results together to get the overall SOH prediction results.In this article, the influences of rest time to capacity regeneration are characterized by regeneration amplitude and regeneration cycle number.As shown in Figure 2,  is the mapping between the rest time and the regeneration cycle number, and  is the mapping between the rest time and the regeneration amplitude.The different working conditions and the individual difference of the current and historical batteries result in the differences of the mappings of the two batteries, which are distinguished by the subscripts c and h, respectively.When the rest time of the current battery is unavailable, the relationship of the regeneration of the current and historical batteries can be associated by mappings  The schematic of similar rest time-based prognostic framework is shown in Figure 3.This framework contains extraction phase and prediction phase.In the extraction phase, SOH series of current battery and SOH series of historical battery can be obtained by capacity monitoring during the usage of the current and historical batteries.The regeneration threshold of the current battery is given by considering the magnitude of noises.Then, the global degradation trend and regeneration phenomena (regeneration amplitude and regeneration cycle number) of the current battery can be decoupled.Based on the similarity between the current and historical batteries, the threshold of the historical battery can be determined automatically by PSO.After that, the regeneration phenomena of the historical battery can be decoupled with this threshold.In the prediction phase, the global The schematic of similar rest time-based prognostic framework is shown in Figure 3.This framework contains extraction phase and prediction phase.In the extraction phase, SOH series of current battery and SOH series of historical battery can be obtained by capacity monitoring during the usage of the current and historical batteries.The regeneration threshold of the current battery is given by considering the magnitude of noises.Then, the global degradation trend and regeneration phenomena (regeneration amplitude and regeneration cycle number) of the current battery can be decoupled.Based on the similarity between the current and historical batteries, the threshold of the historical battery can be determined automatically by PSO.After that, the regeneration phenomena of the historical battery can be decoupled with this threshold.In the prediction phase, the global SOH degradation series, regeneration amplitude series and regeneration cycle number series are predicted, and then integrated together to get the overall prediction results.The extraction phase is described in detail in Figure 4, where H k  can be obtained by capacity monitoring during the usage of the battery.In addition, cycle number series are denoted as The extraction phase is described in detail in Figure 4, where {H ce (k)} n k=1 is the SOH time series of the current battery from cycle 1 to cycle n, which is used to extract the global degradation and regeneration information.To avoid the influence of noises, we consider that capacity regeneration phenomenon has happened only when the value of the increased capacity between two cycles are no less than Th c .In other words, Th c is the regeneration threshold of the current battery.Because the purpose of setting the threshold is to separate regeneration phenomena from noises, when the noise is stronger (or weaker), the value of Th c should be larger (or lower) accordingly.We suggest that Th c can be set as three to five times of the median of the SOH fluctuations caused by noises.
When the SOH of cycle (k + 1) is much better than that of cycle k, there may be significant capacity regeneration during this period.Then the cycle before regeneration can be put into set It is assumed that ce  C contains q elements, which means that obvious regeneration of the current battery has appeared q times.We sort the elements of ce  C in ascending order and get a new series The algorithm to extract regeneration region sets and regeneration cycle number series is as Algorithm 1. {DH ce (k)} n−1 k=1 is the forward difference series of SOH, where DH ce (k) is the SOH increment between cycle k and k + 1 and it is computed as When the SOH of cycle (k + 1) is much better than that of cycle k, there may be significant capacity regeneration during this period.Then the cycle before regeneration can be put into set C ce .
It is assumed that C ce contains q elements, which means that obvious regeneration of the current battery has appeared q times.We sort the elements of C ce in ascending order and get a new series {c ce (k)} q k=1 .The values of regeneration from cycle c ce (k) to cycle (c ce (k) + 1) are described by The algorithm to extract regeneration region sets and regeneration cycle number series is as Algorithm 1. while(c (i) + j ≤ v) and (H(c (i) + j) − H(c (i)) ≥ 0) do 5: put cycle number c (i) + j into set S i 6: j = j + 1 7:

Algorithm 1 Extraction of regeneration region sets and regeneration cycle
end while 8: end for 9: for i = 1 to w − 1 do 10: for j = i + 1 to w do 11: end for 13: series {NUM(w)} w k=1 ; The proposed Algorithm 1 can be called repeatedly in the extraction phase for different inputs.For example, when the algorithm is called at the first time, the input will be w = q, {c (k All the cycle number of global degradation is denoted by set C " ce , which can be expressed as By sorting the g elements of C " ce in ascending order we can get series c ce (k) . With series , the global degradation of SOH can be described by Similarly, elements in SOH forward difference series of the historical battery {DH he (k are shown as: To obtain similar regeneration information from the historical battery, PSO algorithm is adopted to identify the value of the regeneration threshold Th h from the historical battery.Moreover, the range of Th h is given as Th h ∈ [0.1 * Th c , 10 * Th c ].The fitness functions are described by where q is the number of regeneration regions of the current battery and p is the number of regeneration regions of the historical battery within the first n cycles.NU M ce (k) is the regeneration cycle number of the current battery in the kth regeneration region.NU M he (k) is the regeneration cycle number of the historical battery in the kth regeneration region.p and {NUM he (k)} p k=1 can be obtained by extracting regeneration information of the historical battery with PSO algorithm and Algorithm 1. x k i is the particle of PSO algorithm.If the minimum value of F 1 can be obtained from different particles, the particle which makes the value of F 2 minimum is the optimal solution.
The set C he and set C hp denote the cycle before regeneration of the historical battery and they are expressed as follows: It is assumed that the number of the element of set C he and set C hp are p e and p p , respectively.
After sorting the elements of C he and C hp in ascending order, we can obtain series c he (k)  In this article, the working time profiles of the current and historical batteries are similar.Moreover, it is assumed that the relationship of capacity regeneration between two batteries is linear.This means that  In this article, the working time profiles of the current and historical batteries are similar.Moreover, it is assumed that the relationship of capacity regeneration between two batteries is linear.This means that φ c • φ −1 h and ϕ c • ϕ −1 h are constants.The regeneration cycle number series NU M cp (k) p p k=1 can be predicted by linear model and the elements are described by where [ ] is the symbol of floor function.Similarly, the elements of DH cp (k) can be predicted as In addition, Gaussian process model described in Section 2.2 is used in global degradation prediction.
After all prediction processes, the predicted SOH values of each regeneration region are inserted into the global regeneration series.Then the final SOH prediction values of all cycles are estimated.
The SOH series of the ith regeneration region is . Based on the assumption that SOH linearly degrades in the regeneration regions, the value of the element RH (i) cp (k) can be calculated by Series is updated by

Battery Data Set
The Li-ion battery capacity degradation data adopted to conduct the experiments of SOH prediction are obtained from the data repository of the NASA Ame's Prognostics Center of Excellence (PCoE) [36].In the battery prognostics testbed, 18,650 sized Li-ion batteries were run through three operational profiles (charge, discharge, and impedance) at ambient temperature of 4, 24 and 43 • C. Repeated charge and discharge of the batteries resulted in accelerated aging.Charging was performed using a 1.5 A constant current (CC) until the battery voltage of 4.2 V was reached, then in a constant voltage (CV) mode until the current dropped to 20 mA.Discharge was carried out at three constant levels of current (1 A, 2 A, and 4 A) until the discharge voltage fell to 2.2, 2.5 and 2.7 V, respectively.The experimental parameters of the batteries are summarized in Table 1.Detailed information about the testbed can be found in [24].Based on Equation ( 12), SOH can be easily calculated by the experimental capacity data.Figure 6 illustrates the trajectory of SOH for these batteries.It can be noted that the SOH time series presents obvious regeneration phenomena after long time rest.Moreover, the SOH degradation rate in regeneration regions is much faster than in other cycles.In addition, although the degradation rate varies with different end of discharge voltage, after similar long time rest, their capacity regeneration phenomena are similar.Therefore, it is meaningful to extract the similar regeneration information by using the available data of another battery to establish the prediction model of the current battery.

Prediction and Comparison
In this article, the rules of the batteries No. 05, No. 06 and No. 07 are alternatively changed to validate the effectiveness of the proposed method.Specifically, in the first group, battery No. 05 is considered as the current battery and battery No. 07 is assumed as the historical battery.In the second group, battery No. 06 is considered the current battery and battery No. 05 is assumed as the historical battery.In the third group, battery No. 07 is considered as the current battery and battery

Prediction and Comparison
In this article, the rules of the batteries No. 05, No. 06 and No. 07 are alternatively changed to validate the effectiveness of the proposed method.Specifically, in the first group, battery No. 05 is considered as the current battery and battery No. 07 is assumed as the historical battery.In the second group, battery No. 06 is considered the current battery and battery No. 05 is assumed as the historical battery.In the third group, battery No. 07 is considered as the current battery and battery No. 06 is assumed as the historical battery.
For the first group, SOH values of the first 100 cycles of battery No. 05 and all SOH values of battery No. 07 are utilized as the inputs of this framework (n = 100, m = 68).In addition, the regeneration threshold is set at 0.1% (Th c = 0.1%).In the PSO algorithm, the parameters of fitness function F 1 are given as a = 0.8 and b = 0.2.The regeneration influences of the current and historical batteries are shown in Figure 7. From Figure 7, in the 100 training samples, it is noted that the regeneration regions of the historical battery and the current battery are all accurately extracted.It indicates that under this proposed framework, the value of the regeneration threshold determined by PSO algorithm is appropriate.
The regeneration cycle number and the regeneration amplitude of the current battery are predicted by Equations ( 26) and ( 27), respectively.The predicted cycle number and the actual regeneration cycle number of each regeneration region are shown in Figure 8, while the predicted regeneration amplitude and the actual regeneration amplitude are shown in Figure 9.It is noted that the maximum value of prediction errors of the regeneration cycle number is equal to one cycle and the predicted values of the regeneration amplitude are rather close to the corresponding actual values.It demonstrates that it is effective to use linear models to predict regeneration cycle number and regeneration amplitude in this proposed framework.From Figure 7, in the 100 training samples, it is noted that the regeneration regions of the historical battery and the current battery are all accurately extracted.It indicates that under this proposed framework, the value of the regeneration threshold determined by PSO algorithm is appropriate.
The regeneration cycle number and the regeneration amplitude of the current battery are predicted by Equations ( 26) and ( 27), respectively.The predicted cycle number and the actual regeneration cycle number of each regeneration region are shown in Figure 8, while the predicted regeneration amplitude and the actual regeneration amplitude are shown in Figure 9.It is noted that the maximum value of prediction errors of the regeneration cycle number is equal to one cycle and the predicted values of the regeneration amplitude are rather close to the corresponding actual values.It demonstrates that it is effective to use linear models to predict regeneration cycle number and regeneration amplitude in this proposed framework.
regeneration cycle number of each regeneration region are shown in Figure 8, while the predicted regeneration amplitude and the actual regeneration amplitude are shown in Figure 9.It is noted that the maximum value of prediction errors of the regeneration cycle number is equal to one cycle and the predicted values of the regeneration amplitude are rather close to the corresponding actual values.It demonstrates that it is effective to use linear models to predict regeneration cycle number and regeneration amplitude in this proposed framework.The prediction steps of the other two groups are the same as the first one.To quantitatively describe the effectiveness of the proposed method, mean absolute percentage error (MAPE) and root mean square error (RMSE) are taken as the criteria to evaluate the performance of different prognostic methods.The prediction steps of the other two groups are the same as the first one.To quantitatively describe the effectiveness of the proposed method, mean absolute percentage error (MAPE) and root mean square error (RMSE) are taken as the criteria to evaluate the performance of different prognostic methods.
Table 2 presents the comparison results of the performance of nine methods for batteries No. 05, No. 06 and No. 07.The prediction errors of the published methods are from References [15,26,29].It is obvious that the prediction performance of the proposed SRTPF is much better than that of the eight published methods.More specifically, the lowest prediction MAPE on batteries No. 05, No. 06, No. 07 among the eight published methods are 0.82%, 2.28% and 1.02%, respectively, while the counterparts of this proposed framework are 0.62%, 1.41% and 0.76%, respectively.Similarly, for prediction RMSE, we can conclude that the proposed SRTPF can obtain the lowest values for all the three batteries.There are two possible reasons for the excellent prediction performance of SRTPF.The first reason is that this framework can well predict the regeneration information.The second reason is that global SOH degradation trend without the effects of regeneration becomes smoother and easier to be predicted.

Prediction with Various Operation Conditions
To further validate the effectiveness of the SRTPF method, SOH degradation of batteries under different temperatures and operation currents are predicted in this subsection.The BMS has thermal management module like fan and electric heater to do the cooling control and heating control [37].That is to say, under a certain scenario, the temperature around the battery is almost constant.Similarly, the discharge current of the battery can be considered as a constant under a certain scenario.For example, Han [38] has pointed out that a typical driving cycle for a pure EV can be similarly expressed by constant current cycles.Therefore, the key point lies in the validation of this method under various operating ambient temperatures and with different discharge current levels.The following two scenarios are selected to validate the effectiveness of SRTPF.In each scenario, the two batteries are taken as the current and historical batteries alternatively, and the prediction results are illustrated in Figure 12.From the figure, it is noted that the predicted SOH values and the actual values are very close.It means that the proposed method can be well used in different temperatures and operation currents.levels.The following two scenarios are selected to validate the effectiveness of SRTPF.
Scenario 1: For batteries No. 30 and No. 32, the ambient temperature is 43 °C and the discharge current is 4 A.
Scenario 2: For batteries No. 47 and No. 48, the ambient temperature is 4 °C and the discharge current is 1 A.
In each scenario, the two batteries are taken as the current and historical batteries alternatively, and the prediction results are illustrated in Figure 12.From the figure, it is noted that the predicted SOH values and the actual values are very close.It means that the proposed method can be well used in different temperatures and operation currents.

Prediction with Two or More Historical Batteries
The proposed similar rest time-based prognostic framework can realize the prediction of capacity regeneration of the current battery depending only on the degradation SOH data of one historical battery.This is meaningful when the available data sets of historical batteries are not many.It is noteworthy that the first fitness function of PSO (Equation ( 21)) is used for finding the regeneration threshold of the historical battery to obtain the most similar regenerate regions of the current battery.The second fitness function of PSO (Equation ( 22)) is to ensure that the threshold difference of the current and historical batteries is the minimum.Therefore, in the same scenario, when there are two or more available historical batteries have similar rest time to the current battery, the second fitness function F 2 can also be used to select a more similar degradation trajectory to the current battery from them.Then, the capacity data set of this historical battery can be adopted to predict the SOH of the current battery under this framework.
To validate the effectiveness of the historical battery selection based on similarity of degradation trajectory, battery No. 07 is considered as the current battery and batteries No. 05 and No. 06 are assumed as the historical batteries.The values of the parameters are set to the same selected values in Section 4.2.As shown in Table 3, when battery No. 05 is considered as the historical battery, the value of fitness function F 2 is 0; when battery No. 06 is considered as the historical battery, the value of fitness function F 2 is 0.75%.Therefore, compared with battery No. 06, the regeneration phenomenon of battery No. 05 is more similar to that of battery No. 07.The prediction results of battery No. 07 are illustrated in Figure 13.It is noted that the predicted SOH degradation curve of battery No. 07 by using the data of battery No. 05 is closer to the actual curve than by using the data of battery No. 06.As shown in Table 2, when battery No. 05 is considered as the historical battery, the prediction MAPE and RMSE are 0.30% and 0.28%, which are respectively less than 0.76% and 0.83% calculated by considering battery No. 06 as the historical battery.The prediction results can demonstrate that the proposed method can automatically choose a more similar trajectory to predict the regeneration phenomena of the current battery.
degradation trajectory, battery No. 07 is considered as the current battery and batteries No. 05 and No. 06 are assumed as the historical batteries.The values of the parameters are set to the same selected values in Section 4.2.As shown in Table 3, when battery No. 05 is considered as the historical battery, the value of fitness function 2 F is 0; when battery No. 06 is considered as the historical battery, the value of fitness function 2 F is 0.75%.Therefore, compared with battery No. 06, the regeneration phenomenon of battery No. 05 is more similar to that of battery No. 07.The prediction results of battery No. 07 are illustrated in Figure 13.It is noted that the predicted SOH degradation curve of battery No. 07 by using the data of battery No. 05 is closer to the actual curve than by using the data of battery No. 06.As shown in Table 2, when battery No. 05 is considered as the historical battery, the prediction MAPE and RMSE are 0.30% and 0.28%, which are respectively less than 0.76% and 0.83% calculated by considering battery No. 06 as the historical battery.The prediction results can demonstrate that the proposed method can automatically choose a more similar trajectory to predict the regeneration phenomena of the current battery.

Conclusions and Future Work
To obtain more accurate prediction results of battery SOH with capacity regeneration phenomenon, we have proposed a similar rest time-based prognostic framework, in which the regeneration information of a historical battery is used.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.PSO is utilized to determine the threshold of capacity regeneration of the historical battery.Moreover, GP model is used to predict the global SOH degradation trend; linear models are utilized to predict the regeneration cycle number and the regeneration amplitude.
The results of case studies demonstrate that both the global degradation trend and the regeneration phenomena of the testing batteries can be well predicted under different temperatures and discharge currents.Moreover, this proposed framework can obtain more accurate prediction results of battery SOH than the published methods.In addition, when two or more historical batteries are available, the fitness function of PSO can be used to choose a more similar trajectory from them to predict the regeneration phenomena of the current battery.
In future research, additional evaluations need to be conducted on other Li-ion batteries under more complex operating conditions.In addition, more effort should be focused on improving the adaptability of the proposed method.

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

Figure 2 .
Figure 2. Mappings between current and historical batteries.

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

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

Figure 2 .
Figure 2. Mappings between current and historical batteries.

Figure 2 .
Figure 2. Mappings between current and historical batteries.
Symmetry 2016, 9, 4 6 of 20SOH degradation series, regeneration amplitude series and regeneration cycle number series are predicted, and then integrated together to get the overall prediction results.

Figure 3 .
Figure 3. Schematic of similar rest time-based prognostic framework.
the SOH time series of the current battery from cycle 1 to cycle n, which is used to extract the global degradation and regeneration information.To avoid the influence of noises, we consider that capacity regeneration phenomenon has happened only when the value of the increased capacity between two cycles are no less than c T h .In other words, c T h is the regeneration threshold of the current battery.Because the purpose of setting the threshold is to separate regeneration phenomena from noises, when the noise is stronger (or weaker), the value of c T h should be larger (or lower) accordingly.We suggest that c T h can be set as three to five times of the median of the SOH fluctuations caused by noises.Similarly, time series of the historical battery from cycle 1 to cycle (n + m).The data of the first n cycles are mainly used to obtain the value of h Th and they are denoted as k  .The following m cycles of the historical battery are used as the reference to predict the regeneration information of the current battery, denoted as

Figure 3 .
Figure 3. Schematic of similar rest time-based prognostic framework.
Similarly, {H h (k)} n+m k=1 is the SOH time series of the historical battery from cycle 1 to cycle (n + m).The data of the first n cycles are mainly used to obtain the value of Th h and they are denoted as {H he (k)} n k=1 .The following m cycles of the historical battery are used as the reference to predict the regeneration information of the current battery, denoted as H hp (k) m k=1 .Series {H ce (k)} n k=1 , {H he (k)} n k=1 and H hp (k) m k=1 can be obtained by capacity monitoring during the usage of the battery.In addition, cycle number series are denoted as {c e (k)} n k=1 = {1, 2, . . . ,n} and c p (k) m k=1 = {n + 1, n + 2, . . ., n + m}.

Figure 4 .
Figure 4. Detailed description of extraction phase.

Figure 4 .
Figure 4. Detailed description of extraction phase.

.
Thus, regeneration amplitude series DH he (k) p e k=1 and DH hp (k) p p k=1 can be written asDH he (k) = DH he c he (k) (24)DH hp (k) = DH hp c hp (k)(25)By calling Algorithm 1 in the extraction phrase, the regeneration cycle number series {NUM he (k)} p e k=1 and NU M hp (k) prediction phase is shown in Figure 5. Based on the regeneration cycle number prediction algorithm, the series NU M cp (k) p p k=1 can be predicted.In addition, based on the regeneration amplitude prediction algorithm, the series DH cp (k) p p k=1 can be obtained.Taking global SOH degradation series H ce (k) g k=1 as the training samples, the m − ∑ p p l=1 NU M cp (l) steps' global SOH prediction values can be calculated and written as H

Figure 5 .
Figure 5. Detailed description of prediction phase.
can be predicted by linear model and the elements are described by

Figure 5 .
Figure 5. Detailed description of prediction phase.

Figure 6 .
Figure 6.State of health (SOH) values over time: (a) batteries No. 05, No. 06 and No. 07 with calendar time; (b) batteries No. 05, No. 06 and No. 07 with cycle number; (c) batteries No. 30 and No. 32 with calendar time; (d) batteries No. 30 and No. 32 with cycle number; (e) batteries No. 47 and No. 48 with calendar time; and (f) batteries No. 47 and No. 48 with cycle number.

Figure 6 .
Figure 6.State of health (SOH) values over time: (a) batteries No. 05, No. 06 and No. 07 with calendar time; (b) batteries No. 05, No. 06 and No. 07 with cycle number; (c) batteries No. 30 and No. 32 with calendar time; (d) batteries No. 30 and No. 32 with cycle number; (e) batteries No. 47 and No. 48 with calendar time; and (f) batteries No. 47 and No. 48 with cycle number.

Figure 8 .
Figure 8.The actual cycle number and predicted cycle number of each regeneration region.Figure 8.The actual cycle number and predicted cycle number of each regeneration region.

Figure 8 . 20 Figure 9 .
Figure 8.The actual cycle number and predicted cycle number of each regeneration region.Figure 8.The actual cycle number and predicted cycle number of each regeneration region.Symmetry 2016, 9, 4 14 of 20

Figure 9 .
Figure 9.The actual regeneration amplitude and predicted regeneration amplitude.

Figure 10 .
Figure 10.Prediction results of global SOH degradation.

Figure 10 .
Figure 10.Prediction results of global SOH degradation.

Scenario 1 :
For batteries No. 30 and No. 32, the ambient temperature is 43 • C and the discharge current is 4 A. Scenario 2: For batteries No. 47 and No. 48, the ambient temperature is 4 • C and the discharge current is 1 A.

Figure 12 .
Figure 12.Prediction results of SOH: (a) battery No. 30 is the current battery and battery No. 32 is the historical battery; (b) battery No. 32 is the current battery and battery No. 30 is the historical battery; (c) battery No. 47 is the current battery and battery No. 48 is the historical battery; and (d) battery No. 48 is the current battery and battery No. 47 is the historical battery.

Figure 12 .
Figure 12.Prediction results of SOH: (a) battery No. 30 is the current battery and battery No. 32 is the historical battery; (b) battery No. 32 is the current battery and battery No. 30 is the historical battery; (c) battery No. 47 is the current battery and battery No. 48 is the historical battery; and (d) battery No. 48 is the current battery and battery No. 47 is the historical battery.

Table 1 .
Discharging parameters of the batteries.

Table 2 .
Comparison of different prognostic methods for batteries No. 05, No. 06 and No. 07.

Table 3 .
Prediction errors with different historical batteries.

Table 3 .
Prediction errors with different historical batteries.