A Hybrid Method for Remaining Useful Life Estimation of Lithium-Ion Battery with Regeneration Phenomena

The lithium-ion battery has become the primary energy source of many electronic devices. Accurately forecasting the remaining useful life (RUL) of a battery plays an essential role in ensuring reliable operatioin of an electronic system. This paper investigates the lithium-ion battery RUL prediction problem with capacity regeneration phenomena. We aim to reduce the accumulation of the prediction error by integrating different capacity degradation models and thereby improve the prediction accuracy of the long-term RUL. To describe the degradation process more accurately, we decoupled the degradation process into two types: capacity regeneration and normal degradation. Then, we modelled two kinds of degradation processes separately. In the prediction phase, we predicted the battery state of health (SOH) by using the relevance vector machine (RVM) and the gray model (GM) alternately, updated the training dataset according to the prediction results, and then updated the RVM and GM. The RVM and GM correct each other’s prediction results constantly, which reduces the cumulative error of prediction and improves the prediction accuracy of the battery SOH. Experimental results with the National Aeronautics and Space Administration (NASA) battery dataset demonstrated that the proposed method can accurately establish the degradation model and achieve better performance for the RUL estimation as compared with the single RVM or GM methods.


Introduction
With the advantages of high energy/power density, low self-discharge rate, and longevity, lithium-ion batteries have been widely used in consumer electronics, electric vehicles, and even in space systems [1,2]. However, battery performance will degrade with the charging and discharging cycles. In the case of failure, an aged battery may cause significant casualties and economic losses, especially in military communications, navigation, aerospace, and other complex application scenarios [3]. Consequently, the prognostics of the state of health (SOH) of a lithium-ion battery is meaningful for reducing system risks and maintenance costs [4,5].
Generally, the prognostic of a lithium-ion battery is a process of predicting its SOH and remaining useful life (RUL) [6]. The lithium-ion battery RUL is defined as the remaining number of charge-discharge cycles of the battery with a specific output capacity [7]. In recent years, methods for battery prognostics have mainly included a model-based approach, a data-driven approach, and a hybrid approach.
The model-based approach usually includes the electrochemical model [2], the equivalent circuit model [8], and the empirical model. These models usually implement state estimation and prediction However, their iterative strategy could not eliminate the gradually accumulated prediction errors with cycling and aging, which would still affect the accuracy of a long-term prediction.
To obtain a more accurate long-term prediction of RUL, we propose a new hybrid forecasting model using the RVM and the gray model (GM). Our work has improved on the prognostic framework proposed by Qin et al. [30]. In addition, we decoupled capacity regeneration based on rest time. However, this study uses the exponential model to model the regeneration region instead of the linear model, which improves the modeling accuracy. Moreover, the RVM and the GM are used to describe the normal degradation trend. In the prediction phase, we use the RVM and the GM alternately to carry out a multi-step prediction. The prediction results of the RVM and GM model (RVM-GM) update the training dataset and the updated training dataset is used to retain the RVM-GM model. This iterative updating strategy can make the prediction results of the battery's SOH and RUL more accuracy and robust. Finally, we use the NASA battery datasets to verify the prediction performance and effectiveness of the proposed prognostic framework.
The rest of this paper is organized as follows. In Section 2, the algorithm of the RVM and the GM algorithm are introduced. In Section 3, the hybrid model for prognostics of a lithium-ion battery is development. Experimental results and discussion are shown in Section 4 based on the National Aeronautics and Space Administration (NASA) battery dataset. The conclusion is discussed in Section 5.

Relevance Vector Machine
The lithium-ion battery RUL prediction is a regression prediction problem [15]. The RVM is a data-driven method based on the Bayesian theory. It includes the benefits of probabilistic predictions, automatic estimation of "nuisance" parameters, and the facility to utilize arbitrary basis functions [31]. It is valuable for the lithium-ion battery RUL prediction. Thus, before performing this method on the problem, we provide a brief introduction of the components of the RVM.
Assuming a dataset can be defined as {x i , t i } M i=1 , where x i ∈ R n is the input variable, n is the dimension of x i , t i ∈ R is the target value, and M is the length of the training data. In terms of the battery RUL prediction, x i is the battery SOH value at the ith cycle, and t i is the battery SOH value to be predicted. In this paper, the percentage of the current capacity to the rated capacity is taken as the battery health status. The regression expression of the RVM is as follows: where y(x, ω) is a nonlinear function which needs to be modelled. w =[w 0 , w 1 , · · · , w M ] T is the corresponding weight vector. x is all the input variables from the training set. K(x,x i ) is the kernel function. ε i ∼ N 0,σ 2 is the independent additive random noise term. The likelihood function of the training set can be expressed as follows: If Equation (3) is directly used to carry out maximum likelihood estimation for ω and σ 2 , the model will be overfitting and not sparse. Therefore, a constraint of zero mean Gaussian prior distribution is added on ω as follows: where α =[α 0 , α 1 , · · · , α N ] is a hyperparameter vector. Then, we can get the posterior Gaussian distribution of ω from the prior distribution according to the Bayesian theory as follows: where µ is the mean of ω, and Σ is the posterior covariance of ω: where A = diag(α 0 , α 1 , · · · , α M ), Φ presents the N × (N + 1) kernel function matrix. The marginal integral of the hyperparameter can be used to calculate the likelihood distribution of the target value t in Equation (3) as follows: The maximum posterior estimate of the weights depends on α and σ 2 , Tipping [31] gives an iterative formula for calculating α new and σ 2 new : where ii is the diagonal element of Σ, µ i is mean of the ith weight of ω. Then, we can train the RVM model by iterating the above process until the results converge. Finally, the maximum a posteriori estimate of α MP and σ 2 MP are obtained. For a newly updated sample x * , the probability distribution is made for the corresponding target t * : Moreover, p t * |t, α MP , σ 2 MP follows a normal distribution N µ * , σ 2 * , in which

Gray Forecasting Theory
From the original data, the GM finds the law of system change by comparing the similarities and differences between the data characteristics. Thus, the GM weakens the randomness of the original data and improves the prediction accuracy. The gray prediction model predicts the growth trend of the data by establishing a differential equation model from the original data [32]. These advantages make it suitable to predict the long-term RUL of the lithium-ion battery. The following is a brief introduction of the GM.
Assuming the original data sequence can be defined as n . In terms of battery RUL prediction, X (0) represents the sequence of the battery SOH value, and n is the number of training dataset or the length of the sequence. First of all, we should accumulate the original data to i , k= 1, 2, · · · , n. Then, we can establish a first-order differential equation based on the gray sequence X (1) as follows: where a and b are the parameters of the gray model to be calculated. According to the discretized Equation (14), a and b can be calculated by the least square method: Therefore, the sequence generation model can be derived by submitting the starting condition 1 into Equation (16), as follows: Under the starting conditionx 1 , the estimation model of the original sequence is as follows:x Finally, if k > n, we can get the prediction of the original sequence by Equations (19) and (20).

Model Development
Although the RVM and the GM can describe the degradation process of the battery's SOH, their long-term prognostic performance is still unsatisfactory. There are two possible reasons. One reason is that the measurement noise will affect modeling accuracy, and another reason is that the capacity regeneration creates uncertainty in the prediction. Hence, our goal is to get a more accurate RUL prediction model by combining the RVM and the GM. Thus, before developing the proposed hybrid model, the prognostic problems of a lithium-ion battery RUL are explained.

Lithium-Ion Battery RUL Prognostics
Capacity and internal resistance are the two most commonly used health indicators for lithium-ion battery RUL prognostics. The cost of measuring internal resistance online is too high, but the capacity can be calculated easily by integrating the current over time. In order to represent the SOH of a lithium-ion battery for the RUL prognostics in a unified form, capacity is typically converted as follows: where C i and C 0 are the capacity of the cycle i and the rated capacity, respectively. When C i is less than 70% of C 0 , the battery reaches the end-of-life. It is assumed that α is the current charge-discharge cycle and β is the cycle at the end-of-life. By predicting the end-of-life cycle number β, we can get the battery RUL at the α cycle as follows: Figure 1 shows the degradation curve of the #5 lithium-ion battery's SOH with the cycle number and calendar time. The test data were released by the NASA Prognostics Center of Excellence (PCoE) at the Ames Research Center [33]. In Figure 1, DT indicates the rest time of two adjacent discharge cycles. The red triangle shows the capacity regeneration phenomenon. From Figure 1, the capacity regeneration phenomenon affects the state and rate of degradation. When the rest time is long enough, the battery capacity will be regenerated, and the battery's SOH will be increased. Moreover, after a visible regeneration appearance, the SOH degradation rate of the next several cycles is much faster than before. Therefore, it is beneficial to take into account the effect of rest time on the prediction which improves the accuracy of the RUL prediction. Similar to some of the concepts Qin et al. [30] defined, in this paper, the cycle with a significant capacity increase caused by a long rest time is called the "regeneration cycle", the capacity regeneration cycles is called the "regeneration region", the number of cycles in each regeneration region is called the "regeneration cycle number", and the cycle except regeneration regions is called the "normal degradation".
Appl. Sci. 2019, 9, x FOR PEER REVIEW 6 of 16 capacity can be calculated easily by integrating the current over time. In order to represent the SOH of a lithium-ion battery for the RUL prognostics in a unified form, capacity is typically converted as follows: where i C and 0 C are the capacity of the cycle i and the rated capacity, respectively. When i C is less than 70% of 0 C , the battery reaches the end-of-life.
It is assumed that α is the current charge-discharge cycle and β is the cycle at the end-of-life.
By predicting the end-of-life cycle number β , we can get the battery RUL at the α cycle as follows:  Figure 1 shows the degradation curve of the #5 lithium-ion battery's SOH with the cycle number and calendar time. The test data were released by the NASA Prognostics Center of Excellence (PCoE) at the Ames Research Center [33]. In Figure 1, DT indicates the rest time of two adjacent discharge cycles. The red triangle shows the capacity regeneration phenomenon. From Figure 1, the capacity regeneration phenomenon affects the state and rate of degradation. When the rest time is long enough, the battery capacity will be regenerated, and the battery's SOH will be increased. Moreover, after a visible regeneration appearance, the SOH degradation rate of the next several cycles is much faster than before. Therefore, it is beneficial to take into account the effect of rest time on the prediction which improves the accuracy of the RUL prediction. Similar to some of the concepts Qin et al. [30] defined, in this paper, the cycle with a significant capacity increase caused by a long rest time is called the "regeneration cycle", the capacity regeneration cycles is called the "regeneration region", the number of cycles in each regeneration region is called the "regeneration cycle number", and the cycle except regeneration regions is called the "normal degradation".

Hybrid RVM-GM Prediction Method
Our purpose is to get a more accurate battery RUL long-term prediction. Hence, taking into account the effect of the regenerative phenomena,we extract the SOH values of the regeneration cycle, the regeneration cycle number, and normal degradation from the SOH sequence and predict them, respectively.

Hybrid RVM-GM Prediction Method
Our purpose is to get a more accurate battery RUL long-term prediction. Hence, taking into account the effect of the regenerative phenomena, we extract the SOH values of the regeneration cycle, the regeneration cycle number, and normal degradation from the SOH sequence and predict them, respectively. Figure 2 shows the RVM and GM prediction curves of the #5 battery starting prediction at cycle 80. Starting prediction at cycle 80 means the length of the training data is 80, which is a long-term prediction for the battery's RUL. Although the RVM and GM can describe future trends of the degradation, the long-term prediction performance is still unsatisfactory. The SOH prediction results of the RVM are accurate over the first few cycles, but the accuracy will lose or even diverge after a few cycles. This problem occurs because the relevance vectors are easily affected by data noise and the coefficients matrix of the trained RVM cannot be updated during the prediction process. The GM identifies the different degree of the development trend among the system factors to predict the trend of the system, but disregards the uncertainty among the factors in the system. To predict the normal degradation trends more accurately, we propose a hybrid model for fusing the RVM and GM. Figure 3 illustrates the framework of the hybrid RVM-GM prediction method. This framework includes the modeling phase and the prediction phase. The execution steps are as follows:  Figure 2 shows the RVM and GM prediction curves of the #5 battery starting prediction at cycle 80. Starting prediction at cycle 80 means the length of the training data is 80, which is a long-term prediction for the battery's RUL. Although the RVM and GM can describe future trends of the degradation, the long-term prediction performance is still unsatisfactory. The SOH prediction results of the RVM are accurate over the first few cycles, but the accuracy will lose or even diverge after a few cycles. This problem occurs because the relevance vectors are easily affected by data noise and the coefficients matrix of the trained RVM cannot be updated during the prediction process. The GM identifies the different degree of the development trend among the system factors to predict the trend of the system, but disregards the uncertainty among the factors in the system. To predict the normal degradation trends more accurately, we propose a hybrid model for fusing the RVM and GM. Figure. 3 illustrates the framework of the hybrid RVM-GM prediction method. This framework includes the modeling phase and the prediction phase. The execution steps are as follows: Step 1 decoupling the training data into two parts of normal degradation and regeneration.
Step 2 modeling of the capacity regeneration based on rest time.
Step 3 modeling of the normal degradation by the RVM and the GM.
Step 4 predicting the normal degradation and regeneration region, respectively. In particular, the normal degradation is predicted by the alternating iterative method of the RVM-GM.
Step 5 adding the predicted SOH results of each regeneration region into the normal degradation prediction sequence and getting the final SOH prediction results. Step 1 decoupling the training data into two parts of normal degradation and regeneration.
Step 2 modeling of the capacity regeneration based on rest time.
Step 3 modeling of the normal degradation by the RVM and the GM.
Step 4 predicting the normal degradation and regeneration region, respectively. In particular, the normal degradation is predicted by the alternating iterative method of the RVM-GM.
Step 5 adding the predicted SOH results of each regeneration region into the normal degradation prediction sequence and getting the final SOH prediction results.  Figure 2 shows the RVM and GM prediction curves of the #5 battery starting prediction at cycle 80. Starting prediction at cycle 80 means the length of the training data is 80, which is a long-term prediction for the battery's RUL. Although the RVM and GM can describe future trends of the degradation, the long-term prediction performance is still unsatisfactory. The SOH prediction results of the RVM are accurate over the first few cycles, but the accuracy will lose or even diverge after a few cycles. This problem occurs because the relevance vectors are easily affected by data noise and the coefficients matrix of the trained RVM cannot be updated during the prediction process. The GM identifies the different degree of the development trend among the system factors to predict the trend of the system, but disregards the uncertainty among the factors in the system. To predict the normal degradation trends more accurately, we propose a hybrid model for fusing the RVM and GM. Figure. 3 illustrates the framework of the hybrid RVM-GM prediction method. This framework includes the modeling phase and the prediction phase. The execution steps are as follows: Step 1 decoupling the training data into two parts of normal degradation and regeneration.
Step 2 modeling of the capacity regeneration based on rest time.
Step 3 modeling of the normal degradation by the RVM and the GM.
Step 4 predicting the normal degradation and regeneration region, respectively. In particular, the normal degradation is predicted by the alternating iterative method of the RVM-GM.
Step 5 adding the predicted SOH results of each regeneration region into the normal degradation prediction sequence and getting the final SOH prediction results. We assume that the battery is fully charged and discharged in each cycle. Thus, we can get the time sequence of the SOH {H k } α k=1 and the start time sequence of two adjacent discharge cycles from cycle 1 to cycle α, where α is the length of training data.
In the modelling phase, we decouple normal degradation and regeneration from the training data, using s k to distinguish whether there is a visible SOH regeneration at the kth cycle.
where DH k = H k+1 − H k is the forward difference of the SOH between the kth and the (k + 1)th cycle, DH th is the threshold for capacity regeneration. s k = 1 means normal degradation, and s k = −1 means that there is a visible SOH regeneration from cycle kth to (k + 1)th. The set of cycle indices before a regeneration phenomenon can be represented as A = {k|s k = −1, k ∈ {1, 2, · · · , α − 1}}. The number of elements in the set A is N a , which indicates the number of regeneration phenomena in the training data. The cycle index of the regeneration region can be extracted according to A, which can be represented as B = k + j H k+ j > H k , k ∈ A, j ∈ N + . Continuous cycles in B means a regeneration region. According to the division, the cycles of each regeneration region can be determined. As a result, the cycle index of the normal degenerate can be extracted, which can be represented as E = {k|k ∈ {1, 2, · · · , α}, k B}.
Based on the decoupled data, the RVM and the GM for normal degradation and the rest time-based model for regeneration phenomenon are trained. Moreover, the training data set for the RVM is R RVM = (H k , H k+1 ) k, k + 1 ∈ E , and the original sequences training data set for the GM is To predict the regeneration phenomenon, we need to model the regeneration amplitude and the regeneration cycle number, respectively. Figure 4 shows the data points and fitting curves of the rest time and regeneration amplitude of the experimental data of four groups of batteries, indicating a nonlinear relationship. Furthermore, it can be seen that after a visible regeneration, the SOH degradation rate at the next several cycles will be much faster than before. Qin et al. [30] assumed that the growth rate of the regeneration amplitude decreased with an increase in rest time and that the degradation of the capacity regeneration region was a uniform linear process, and therefore the regeneration amplitude and the regeneration cycle number are accurately predicted. However, the modeling of the degradation process in the regeneration region was not accurate. By observing the data of [33], we find that the degradation process of the regeneration zone was decreasing exponentially. Therefore, by obtaining the regeneration amplitude and the regeneration quantity prediction model, the regeneration area is modelled using the exponential model.
To predict the regeneration phenomenon, we need to model the regeneration amplitude and the regeneration cycle number, respectively. Figure 4 shows the data points and fitting curves of the rest time and regeneration amplitude of the experimental data of four groups of batteries, indicating a nonlinear relationship. Furthermore, it can be seen that after a visible regeneration, the SOH degradation rate at the next several cycles will be much faster than before. Qin et al. [30] assumed that the growth rate of the regeneration amplitude decreased with an increase in rest time and that the degradation of the capacity regeneration region was a uniform linear process, and therefore the regeneration amplitude and the regeneration cycle number are accurately predicted. However, the modeling of the degradation process in the regeneration region was not accurate. By observing the data of [33], we find that the degradation process of the regeneration zone was decreasing exponentially. Therefore, by obtaining the regeneration amplitude and the regeneration quantity prediction model, the regeneration area is modelled using the exponential model.  The rest time is defined as the forward difference of the beginning time of each cycle, which can be represented as: The rest times dataset for estimating regeneration phenomena is where s is a scale parameter, φ(•) is a hyperbolic tangent function, NUM e i is the regeneration cycle number to be estimated, N a is the number of times that regeneration occurs in the training data, and NUM j is the cycle number of the jth regeneration in the training data.
When using the exponential degradation model to model the regeneration region, it is assumed that j is the cycle number before the ith regeneration and the health state sequence of the regeneration region can be expressed as follows: where a and b are undetermined parameters.
In the prediction phase, we use the RVM and the GM alternately to carry out a multi-step prediction. We suppose that the plan for future use of the battery is known, which means {T k } n k=α+1 is known.
The rest time sequence DT p k n−1 k=α+1 used for prediction is obtained by Equation (26). Because the SOH of the prediction process is unknown, we use the rest time threshold DT th to indicate the occurrence of regeneration.
The regeneration amplitude set DH  (25) and (26), where M is the number of times that regeneration occurs in the prediction process. Our proposed hybrid prediction algorithm is summarized as Algorithm 1.
During the prediction process, all the new prediction results are added to the original sequence to update the GM. However, in contrast to the iterative updating strategy proposed by [15], which adds all the prediction results to the training dataset of the RVM, all we need is the original dataset and the latest prediction results. Because the GM cannot carry out an accurate local prediction, adding all the GM prediction results in the training set of the RVM will reduce the prediction accuracy.

Experiment Battery Dataset
The effectiveness of our proposed method is verified by the lithium-ion battery experimental data provided by the NASA PCoE [33]. This experimental dataset which consisted of 18,650 commercial lithium-ion battery sets was obtained by charging and discharging in cycles at room temperature. Details of the experiment can be found in reference [33]. We used the experimental data of batteries #5, #6, #7, and #18. The batteries reached the end-of-life when their capacity was reduced to 30% of the rated capacity. Figure 5 shows the SOH degradation curves of batteries #5, #6, #7, and #18. We can see that each battery has some visible SOH regeneration, and the SOH degradation rate at the next several cycles of regeneration will be much faster than before. In our experiment, the initial capacity value is set to 2 Ah and the end-of-life capacity is set to 1.4 Ah. The battery data length for batteries #5, #6, and #7 is 168, and #18 is 132. The actual life cycle of batteries #5, #6, and #18 is 124, 109, and 121, respectively. For battery #7, the capacity decreases to 1.432 Ah at the end. When we set its end-of-life capacity threshold at 72%, then its actual life cycle is 146.  Figure 6 shows the results of normal degradation and capacity regeneration decoupling of the first 100 cycle data of each battery. We can see that the starting time of the capacity regeneration in batteries #5, #6, and #7 is the same, because the starting time of each charge-discharge cycle is the same during the three battery tests. However, due to the influence of initial conditions and other uncertainties in the test process, the regeneration amplitude and quantity of those three batteries are different.  Figure 6 shows the results of normal degradation and capacity regeneration decoupling of the first 100 cycle data of each battery. We can see that the starting time of the capacity regeneration in batteries #5, #6, and #7 is the same, because the starting time of each charge-discharge cycle is the same during the three battery tests. However, due to the influence of initial conditions and other uncertainties in the test process, the regeneration amplitude and quantity of those three batteries are different. Figure 6 shows the results of normal degradation and capacity regeneration decoupling of the first 100 cycle data of each battery. We can see that the starting time of the capacity regeneration in batteries #5, #6, and #7 is the same, because the starting time of each charge-discharge cycle is the same during the three battery tests. However, due to the influence of initial conditions and other uncertainties in the test process, the regeneration amplitude and quantity of those three batteries are different. The absolute error of the RUL prediction results and the root mean square error of the SOH prediction results are used to quantify and evaluate the prediction results, which are defined as follows: The absolute error of the RUL prediction results and the root mean square error of the SOH prediction results are used to quantify and evaluate the prediction results, which are defined as follows:

Experiment and Discussion
The prediction comparison for the RVM, the GM, and the hybrid method for batteries #5, #6, #7, and #18 are shown in Table 1 and Figure 7. The prediction start cycle means predicting the RUL at the current cycle, and the length of the training data is equal to the started prediction cycle. It can be seen from Table 1 that the RUL prediction errors of our hybrid are less than the single RVM and the single GM method at a different started predicton cycle. Moreover, the AE and RMSE are also smaller than the single model. The error message illustrates that our method has high accuracy over the prediction period. The performance of the hybrid model is still satisfied with the expansion of the prediction horizon. More specifically, the prognostic AE started prediction at cycle 60 of batteries #5, #6, #7, and #18 among the single RVM methods 54, 15, 74, and 1, respectively. But, the counterparts of our hybrid model are 4, 15, 1, and 4, respectively. Although the AE on battery #18 is slightly larger, we can see that the prognostic of the proposed method is more accurate and stable. Similarly, as compared with the single GM method, the prognostic results on four batteries are the smallest.    To compare more broadly with other types of prediction methods, we investigated many other model-based and data-driven approaches. The RVM-KF method proposed in references [15] is based on the empirical model, RVM and Kalman filter for lithium-ion battery prognostic without considering the influence of the regeneration phenomenon, although they did not give the experimental results of battery #18, which was heavily affected by the regeneration phenomenon. The proposed hybrid method not only accurately predicts the local regeneration phenomenon but also has higher accuracy in predicting the battery RUL for short-term and long-term. Table 2 shows the comparison of prediction results of the hybrid method with the RVM-KF method. It should be noted that the RVM-KF algorithm has better prediction performance than the proposed hybrid method started prediction at cycle 80 of battery #6. One possible reason is that the hybird method does not accurately predict the abnormal regeneration, and the end-of-life is in this inaccurate prediction. The rest time-based prognostic method proposed in reference [30] uses the data-driven method of the SVM and the Gaussian process model for lithium-ion battery prognostic.However, they only gave short-term prediction results for started prediction at the start cycle 100. The present experiment is started at different stages of battery life, and the prediction performance is better in both the short-and long-term prognostic.
Moreover, as compared with the electrochemical model [2] and the equivalent circuit model [8,25,34] methods that provide high-precision short-term SOH prediction, the proposed hybrid method not only provides very closed short-term SOH prediction accuracy, but also provides effective long-term RUL prediction results.

Conclusions
To improve the long-term RUL prediction performance of the lithium-ion battery with capacity regeneration phenomenon, a hybrid prediction method based on the RVM and the GM algorithm is explored. The proposed method decouples regeneration phenomenon and normal degradation and predicts them separately. After obtaining the amplitude of regeneration and the regeneration cycle number, an exponential model is proposed to describe the degraded region. Moreover, an iterative updating method is proposed to predict the normal degradation. The RVM and the GM are alternately used to carry out a multi-step prediction. The RVM-GM model prediction results update the training dataset, and the updated training dataset retrains the RVM-GM model. Finally, we add the prediction results of each regeneration region into the normal degradation prediction sequence to get the overall SOH prediction results.
As compared with a pure RVM and a pure GM, our hybrid method has higher accuracy and stability. Moreover, the RVM and the GM in the proposed hybrid method only rely on historical data, and it can construct a prognostic model without any prior knowledge about the complex battery degradation. Thereby, this method can predict the RUL of different types of batteries with historical data accumulation.
The limitations of our hybrid method in real applications are that the training data should have typical data reflecting the regenerated phenomena, and the future battery use plan should be known to determine the rest time in the prediction phase. In addition, other factors such as temperature and hysteresis can also affect the accuracy of the prediction. Therefore, our future work will investigate how to build the degradation model in the case of insufficient data and will consider the influence of unknown rest time, temperature, data lag, and other factors on the prediction uncertainty.
Author Contributions: Y.W. and L.Z. conceived the method and structured the whole paper; Y.W. and J.C. performed the experiments, analyzed the data, wrote and revised the paper.
Funding: This research was funded by the National Natural Science Foundation of China (Nos. 61633008, 61773132 and 61803115).