State of Health Prediction of Lithium-Ion Batteries Based on the Discharge Voltage and Temperature

Accurate state of health (SOH) prediction of lithium-ion batteries is essential for battery health management. In this paper, a novel method of predicting the SOH of lithium-ion batteries based on the voltage and temperature in the discharging process is proposed to achieve the accurate prediction. Both the equal voltage discharge time and the temperature change during the discharge process are regarded as health indicators (HIs), and then, the Pearson and Spearman relational analysis methods are applied to evaluate the relevance between HIs and SOH. On this basis, we modify the relevance vector machine (RVM) to a multiple kernel relevance vector machine (MKRVM) by combining Gaussian with sigmoid function to improve the accuracy of SOH prediction. The particle swarm optimization (PSO) is used to find the optimal weight and kernel function parameters of MKRVM. The aging data from NASA Ames Prognostics Center of Excellence are used to verify the effectiveness and accuracy of the proposed method in numerical simulations, whose results show that the MKRVM method has higher SOH prediction accuracy of lithium-ion batteries than the relevant methods.


Introduction
Lithium-ion batteries have been widely used in aviation [1], electronic technology, and other fields [2] due to the high energy density, long service life, and low self-discharge rate. However, the capacity of the lithium-ion battery will degenerate after repeated charge and discharge. When the capacity of the lithium-ion battery degenerates a threshold, the lithium-ion battery should be replaced to ensure system safety and avoid the occurrence of serious accidents. Therefore, we need to understand and take it as a guide to provide security for the system operation, which makes it of great significance to study the state of health (SOH) of lithium-ion batteries [3].
The prediction of SOH is still facing great challenges. The prediction methods of SOH for lithium-ion batteries are mainly divided into three types: model-based method [4], data-driven method [5,6], and fusion method [7]. The model-based method is mainly to establish a mathematical model of describing the failure of the battery through the prior knowledge in the battery degradation e.g., equivalent circuit model [8], empirical model [9], electrochemical model [10], etc. The data-driven method uses machine learning techniques (such as Gaussian process regression (GPR) [11,12], support vector regression (SVR) [13], relevance vector machine (RVM) [14,15], artificial neural network (ANN), deep learning [16,17], and so on) to predict the SOH based on the battery degradation data. The prediction can be divided into direct prediction and indirect prediction. Direct prediction is to predict the SOH of lithium-ion batteries by taking capacity or internal resistance as a health indicator (HI). For example, Qu et al. [18] proposed complete ensemble empirical mode decomposition with adaptive noise based on capacity and combined long short-term memory with particle swarm optimization (PSO) as well as an attention mechanism for battery degradation process, we utilize the global sigmoid kernel function and the local Gaussian kernel function to improve the performance of the RVM model to achieve higher accuracy of SOH prediction. Moreover, the datasets provided by NASA are used to verify the effectiveness and accuracy of the method proposed in this paper.
The remainder of this paper is organized as follows. In Section 2, the lithium-ion battery datasets, as well as the feature extraction process and correlation analysis are described. The method of MKRVM is introduced in Section 3. The experimental result verification and analysis are provided in Section 4. Finally, the conclusion is discussed in Section 5.

Feature Extraction and Analysis Based on Discharge Curves
In this section, the features, which can reflect the battery degradation during the discharging process, are extracted for SOH prediction.

The Definition of SOH
SOH is a significant indicator that can be used to measure the health of the lithiumion battery, which leads to the accurate prediction of SOH being able to guide prognosis management. The SOH of lithium-ion batteries can be defined based on the capacity, internal resistance, cycle number, and so on. In this paper, SOH is defined by using the battery capacity as follows [34]: where C i is the capacity at the ith cycle, and C 0 is the initial capacity.

Experiment Data
The effectiveness of the MKRVM prediction model is evaluated based on the data from the data repository of the NASA Ames Prognostics Center of Excellence [35]. This datasets are collected from a battery prognostic test-bed at NASA 18650 comprising commercially available lithium-ion rechargeable batteries whose nominal capacity is 2 Ah and nominal voltage is 3.7 V [36,37]. In order to verify the effectiveness of the MKRVM, a set of batteries #5, #6, #7, and #18, and another set of batteries #54 and #55 in three different operating conditions are used in this paper. Lithium-ion batteries are run through three different operational profiles (charge, discharge, and impedance). Batteries #5, #6, #7, and #18 are tested at room temperature (24 • C). First, charging is carried out in a constant current mode at 1.5 A until the battery voltage reached 4.2 V and then continued in a constant voltage mode until the charge current dropped to 20 mA. After that, the discharge process is carried out at a constant current of 2 A until the battery voltage falls to 2.7, 2.5, 2.2, and 2.5 V for batteries #5, #6, #7, and #18, respectively. Thus, batteries #5, #6, and #7 are cycled with different depth of discharge (DoD), batteries #6 and #18 are cycled with the same DoD. However, the cycles corresponding to the failure threshold of battery #18 are much shorter than those of battery #6 due to the different types of batteries. Meanwhile, batteries #54 and #55 discharged at a constant current of 2A stop at 2.2 V and 2.5 V in a low-temperature environment (4 • C). Repeated charge and discharge cycles result in accelerated aging of batteries, and the capacity with cycle numbers are shown in Figure 1. In Figure 1b, the first point and the last point have been removed from the capacity data due to a big shift for batteries #54 and #55. For batteries #54 and #55, we can observe from Figure 1b that the battery with many local noises due to being under low temperature, batteries #54 and #55, are used to verify the extrapolation performance of the method proposed in this paper.

Feature Extraction Based on Discharging Curves
Feature extraction can be regarded as an information simplification process to extract meaningful signals from the original data to reduce the input size without affecting the performance of the model. There are many factors affecting battery degradation, such as temperature, current, voltage, cyclic life, etc. Since the current is constant in the discharging process, the voltage and temperature are considered in this paper. Meanwhile, considering the battery will not start charging until the cut-off voltage, part of the data of the low state of charge (SOC) area is discarded. In addition, when the battery is open circuit and not in a working state, which with the existence of self-discharge, the battery may not be in a high SOC region during working. In addition, using partial data can reduce the calculation time and solve information redundancy. Taking the data of battery #5 as an example, the feature extraction is introduced in Figure 2, in which Figure 2a shows the voltage changing of battery #5 at different cycle numbers and Figure 2b shows the temperature changes in different cycles. It can be observed from Figure 2 that the time to reach a specific voltage is different for different cycles, and the voltage variation between 3.0 and 3.8V covers the main information of the discharge process [38].
(1) TIEDVD: From Figure 2a, we can see that the time of different cycles in the same voltage interval is different. As the cycle of battery charge and discharge increases, the voltage drop rate is different. Taking the 150th cycle as an example, t1(150) corresponds to the time when the battery voltage reaches 3V, t2(150) corresponds to the time when the battery voltage reaches 3.8V, and TIEDVD is equal to t1-t2. The result is shown in Figure 3a,we can see that TIEDVD is consistent with the capacity change trend and local fluctuations. The TIEDVD and TIEDTD, which are selected from voltage and temperature, are used as inputs of the MKRVM model to obtain the SOH prediction of lithium-ion batteries.

Feature Extraction Based on Discharging Curves
Feature extraction can be regarded as an information simplification process to extract meaningful signals from the original data to reduce the input size without affecting the performance of the model. There are many factors affecting battery degradation, such as temperature, current, voltage, cyclic life, etc. Since the current is constant in the discharging process, the voltage and temperature are considered in this paper. Meanwhile, considering the battery will not start charging until the cut-off voltage, part of the data of the low state of charge (SOC) area is discarded. In addition, when the battery is open circuit and not in a working state, which with the existence of self-discharge, the battery may not be in a high SOC region during working. In addition, using partial data can reduce the calculation time and solve information redundancy. Taking the data of battery #5 as an example, the feature extraction is introduced in Figure 2, in which Figure 2a shows the voltage changing of battery #5 at different cycle numbers and Figure 2b shows the temperature changes in different cycles. It can be observed from Figure 2 that the time to reach a specific voltage is different for different cycles, and the voltage variation between 3.0 and 3.8 V covers the main information of the discharge process [38].
(1) TIEDVD: From Figure 2a, we can see that the time of different cycles in the same voltage interval is different. As the cycle of battery charge and discharge increases, the voltage drop rate is different. Taking the 150th cycle as an example, t1(150) corresponds to the time when the battery voltage reaches 3V, t2(150) corresponds to the time when the battery voltage reaches 3.8V, and TIEDVD is equal to t1-t2. The result is shown in Figure 3a,we can see that TIEDVD is consistent with the capacity change trend and local fluctuations. (2) TIEDTD: During the time of voltage change, we also analyze the temperature change during this period. T1(150) is the point corresponding to t1(150), T2(150) is the point corresponding to t2(150), and TIEDTD is equal to T2-T1. From Figure 2b, we can see that temperature change simultaneously is not the same. So, we extract the temperature change within the same time as other HI. The normalized HIs are shown in Figure 3b, while TIEDTD and the capacity change trend show the opposite trend in this paper.
The TIEDVD and TIEDTD, which are selected from voltage and temperature, are used as inputs of the MKRVM model to obtain the SOH prediction of lithium-ion batteries.

Feature Analysis Based on Pearson and Spearman
It can be seen from Figure 3 that the TIEDVD and the capacity degradation trend are the same. While the temperature changes within the same time, as the cycle changes, it shows the opposite trend to the capacity, and the fluctuation part corresponds to it. All the tested batteries show a similar change to that of battery #5, as shown in Figure 3. To verify the effectiveness of the extracted HIs, the correlation analysis between HIs and SOH is presented as follows.
Pearson Correlation Analysis: Assuming that there are two variables and , the Pearson correlation coefficient between the two variables can be calculated by the following formula: . (2) In this paper, is the SOH, and is the HI. The value of the correlation coefficient is between −1 and +1. The greater the absolute value of the correlation coefficient, the stronger the correlation. The correlation analysis between the HIs and SOH is shown in Table 1. It can be seen from Table 1 that the absolute values of the Pearson correlation coefficient of HIs proposed in this paper are above 0.9, and it can be seen that HIs are strongly correlated with SOH.

Feature Analysis Based on Pearson and Spearman
It can be seen from Figure 3 that the TIEDVD and the capacity degradation trend are the same. While the temperature changes within the same time, as the cycle changes, it shows the opposite trend to the capacity, and the fluctuation part corresponds to it. All the tested batteries show a similar change to that of battery #5, as shown in Figure 3. To verify the effectiveness of the extracted HIs, the correlation analysis between HIs and SOH is presented as follows.
Pearson Correlation Analysis: Assuming that there are two variables and , the Pearson correlation coefficient between the two variables can be calculated by the following formula: . (2) In this paper, is the SOH, and is the HI. The value of the correlation coefficient is between −1 and +1. The greater the absolute value of the correlation coefficient, the stronger the correlation. The correlation analysis between the HIs and SOH is shown in Table 1. It can be seen from Table 1 that the absolute values of the Pearson correlation coefficient of HIs proposed in this paper are above 0.9, and it can be seen that HIs are strongly correlated with SOH.

Feature Analysis Based on Pearson and Spearman
It can be seen from Figure 3 that the TIEDVD and the capacity degradation trend are the same. While the temperature changes within the same time, as the cycle changes, it shows the opposite trend to the capacity, and the fluctuation part corresponds to it. All the tested batteries show a similar change to that of battery #5, as shown in Figure 3. To verify the effectiveness of the extracted HIs, the correlation analysis between HIs and SOH is presented as follows.
Pearson Correlation Analysis: Assuming that there are two variables X and Y, the Pearson correlation coefficient between the two variables can be calculated by the following formula: In this paper, X is the SOH, and Y is the HI. The value of the correlation coefficient is between −1 and +1. The greater the absolute value of the correlation coefficient, the stronger the correlation. The correlation analysis between the HIs and SOH is shown in Table 1. It can be seen from Table 1 that the absolute values of the Pearson correlation Electronics 2021, 10, 1497 6 of 18 coefficient of HIs proposed in this paper are above 0.9, and it can be seen that HIs are strongly correlated with SOH. Spearman Correlation Analysis: Spearman uses the rank of two variables to make a linear correlation analysis, and it belongs to a non-parametric statistical method that has a wider range of applications [39,40]. Spearman's calculation formulas are as follows: n is the number of variables, and ρ(a, b) is the Spearman rank correlation coefficient. The value of the Spearman correlation coefficient is also from −1 to +1, when the absolute value of correlation is close to 1, the greater the correlation between variables X and Y.
From Table 1, one can see that the absolute value of Pearson and Spearman correlation coefficients of the TIEDVD and TIEDTD and SOH are all above 0.9. Thus, the HIs chosen in this paper are highly correlated with SOH and are effective, so that it is reasonable for building an RVM model and achieving accurate SOH predictions.

Methodologies
This section mainly describes the method proposed in this paper. The RVM method has been widely used in SOH prediction [41]. Compared with SVM, RVM can provide probability interpretation under the Bayesian framework, and it can also achieve a better generalization performance with fewer kernel functions. So, we use RVM and describe the MKRVM method in this section.

MKRVM
..N denote a set of input vectors, and their corresponding targets where N is the length of the training data. In this paper, x i represents HI, and SOH is denoted as t i . The RVM learns the relationship between the target variable t and the input variable x to accurately predictt throughx. We assume the relationship between the input vector x and the true output vector t can be written as: In Equation (4), ε i is the independent errors following ε i ∼ N 0, σ 2 . In Equation (5), K(x, x i ) is the kernel function, and w i = (w 1 , w 2 , . . . , w N ) T is the corresponding weight, and then, the likelihood can be written as: Electronics 2021, 10, 1497 To prevent the over-fitting problem caused by maximizing the likelihood, in which a zero-mean Gaussian prior distribution N over w is defined as: With η, a vector of N + 1 hyperparameters, which is also the target of our optimization. Then, the posterior distribution over the weights is given the data from Bayes' rule: where the posterior covariance matrix and mean vector are: Then, to obtain the value of η and σ 2 , whose most-probable (MP) values are calculated with an iterative re-evaluation algorithm: With a new inputx, the corresponding output distribution can be given as follows: where the mean and the variance are respectively: It can be seen from Equation (14) that the selection of the kernel function is very important for the establishment of the model. Compared with a single kernel function, a multi-kernel functions has higher flexibility. The generalization ability of a single kernel function is weak. Combining multiple different kernel functions to get better characteristics is the basic idea of a mixed kernel function. Therefore, the Gaussian kernel function and sigmoid function as the hybrid kernel functions of MKRVM, and the expression of the mixed kernel function is: where K j is the basic kernel function, and m j is the weight coefficient, that is, the combination coefficient. In this paper, a mixed kernel function of Gaussian and sigmoid is used. The Gaussian kernel function has strong nonlinear processing capabilities, and it is often used as a kernel function: Among them, d is the bandwidth, which is used to control the local scope of the Gaussian kernel function.
The sigmoid kernel function from the neural network has been proved to have good classification performance [42], and it can be expressed as: where t is a scalar, which is an amplitude adjustment parameter of the input data, and c is a displacement parameter that controls the mapping threshold. Among them, the combination coefficient m and the kernel function parameters d, c are to be optimized and determined. In order to determine the optimal kernel function weights and kernel function parameters, PSO is applied to optimize them, so as to reduce the model generalization error and prevent over-fitting.

PSO
PSO is a global random search algorithm based on swarm intelligence by simulating the migration and clustering behavior of birds in the process of foraging. It has been widely used in function optimization, neural network training, fuzzy system control, and other application areas.
The procedures of the PSO algorithm are as follows: Step 1. Initialize the particle swarm position and speed.
Step 2. Calculate the fitness value of each particle.
Step 3. Compare the current position of each particle with the historical optimal position fitness, and the optimal position is regarded as the individual optimal position. Step 4. Compare the optimal position of each particle with the optimal individual position in the group, and the optimal position is obtained as the global optimal position. Step 5. Update particle speed and position.
x k+1 where υ k i , x k i , and p k i are the particle velocity, current location, and best location in the history of particle i when the iteration is k, respectively. c 1 and c 2 are accelerating factors, λ 1 and λ 2 are random numbers between 0 and 1.
With this algorithm, the optimal parameters of MKRVM can be found automatically, and the prediction accuracy of the model can be further improved. The PSO algorithm flow chart is shown in Figure 4.

SOH Prediction Based on the PSO-MKRVM
The flow chart of SOH prediction based on MKRVM is shown in Figure 5. First, we extract the corresponding HIs from the lithium-ion battery discharge data, perform correlation analysis, and then use the extracted HIs as the inputs of the MKRVM model for training; finally, the effectiveness of the proposed method is verified by the test sets. Here, the following experiments are based on Matlab 2019a on a laptop with an Inter Pentium 2.6 GHz (i7-9750H). factors, and are random numbers between 0 and 1. With this algorithm, the optimal parameters of MKRVM can b and the prediction accuracy of the model can be further improve flow chart is shown in Figure 4.  The flow chart of SOH prediction based on MKRVM is shown in Figure 5. First, we extract the corresponding HIs from the lithium-ion battery discharge data, perform correlation analysis, and then use the extracted HIs as the inputs of the MKRVM model for training; finally, the effectiveness of the proposed method is verified by the test sets. Here, the following experiments are based on Matlab 2019a on a laptop with an Inter Pentium 2.6GHz (i7-9750H).

SOH Prediction Based on the PSO-MKRVM
In Figure 5, the details of battery SOH prognostic are described as follows: (1) Feature extraction: TIEDVD and TIEDTD are extracted from the discharge data of the lithium-ion battery, .
(2) Calculation of the SOH of the battery: .
(3) Division of the dataset: Divide the dataset into the training set and test set, set the prediction starting point to ; then, the training set input is , the corresponding target variable is , and the test set can be described as , .

Evaluation Indicators
In order to evaluate the performance of the MKRVM method, root mean square error (RMSE) and mean absolute percentage error (MAPE) act as evaluation indicators in this paper, whose calculation formulas are given as follows:

Evaluation Indicators
In order to evaluate the performance of the MKRVM method, root mean square error (RMSE) and mean absolute percentage error (MAPE) act as evaluation indicators in this paper, whose calculation formulas are given as follows: whereŷ i denotes the predicted SOH, and y i denotes the real SOH.

Numerical Experiments and Results Analysis
To verify the method proposed in this paper, batteries #5, #6, #7, and #18 are tested to prove the performance of MKRVM and the results with different models at the same starting points are analyzed, while the experiments at different starting points are also carried out. In addition, batteries #54 and #55 are tested with the MKRVM method at different starting points to verify the MKRVM for SOH predictions and analyze the results.

SOH Prediction at Same Starting Points for Different HIs
In order to present the advantages of the multi-kernel functions, we take the battery #5 as an example, the starting point is cycle 80, and the prediction results of a single TIEDVD, TIEDTD, and the combination of TIEDVD and TIEDTD are compared and shown in Figure 6, while the corresponding RMSE and MAPE are shown in Table 2. It can be seen from Figure 6 and Table 2 that the prediction result of the two HIs combined together is the best, and the RMSE and MAPE are also the smallest, with the values of 0.004 and 0.2932, respectively. The prediction result of TIEDTD is the worst, and the RMSE is 0.0159.
TIEDVD only considers the impact of voltage on capacity degradation, while TIEDTD only considers the influence of temperature on battery degradation. However, the influencing factors of battery degradation are related to both voltage and temperature, so TIEDVD combined with TIEDTD can not only accurately predict the trend but also capture the capacity recovery. Therefore, TIEDVD and TIEDTD are regarded as the inputs of MKRVM to improve the prediction accuracy of SOH in this paper. so TIEDVD combined with TIEDTD can not only accurately p capture the capacity recovery. Therefore, TIEDVD and TIEDTD of MKRVM to improve the prediction accuracy of SOH in this p

SOH Prediction at Same Starting Points for Different Models
In order to verify the effectiveness of the MKRVM method els at the same starting points and compared them with SVR an widely used to predict the SOH of lithium-ion batteries. Figur results of SOH by using the proposed method, RVM with Gaus with Gaussian kernel function, and SVR with Gaussian functio are denoted as M1, M2, M3, and M4, respectively. In Figure 7, for batteries #5, #6, and #7, while the start point is cycle 60 for b ison experiments, M2 is considered to illustrate the affection of

SOH Prediction at Same Starting Points for Different Models
In order to verify the effectiveness of the MKRVM method, we set up different models at the same starting points and compared them with SVR and GPR models, which are widely used to predict the SOH of lithium-ion batteries. Figure 7 shows the prediction results of SOH by using the proposed method, RVM with Gaussian kernel function, GPR with Gaussian kernel function, and SVR with Gaussian function for four batteries, which are denoted as M1, M2, M3, and M4, respectively. In Figure 7, the start point is cycle 80 for batteries #5, #6, and #7, while the start point is cycle 60 for battery #18. In the comparison experiments, M2 is considered to illustrate the affection of the kernel function in the proposed method, and M3 and M4 are also tested to analyze the affection of RVM in the proposed model. Figure 7 shows the prediction results of different models for the four batteries, from which one can see that all the four models can predict the trend of SOH well, indicating the advantage of machine learning. However, the results of M2, M3, and M4 deviate from the true value; since multi-kernel functions can learn better than the single kernel function, M3 and M4 can provide better results than M2, which indicates that RVM can be improved a lot. It can be observed from Figure 7 that the prediction results by using MKRVM are close to the actual values and more accurate than those by using other methods, which indicates that MKRVM effectively provides the accurate prediction. Therefore, the results demonstrate that the multi-kernel functions of the proposed method effectively obtain the improved prediction performance. MKRVM, which is the improved RVM, the RMSE is the best. It can be observed tha MKRVM method can improve the performance of SOH prediction.

SOH Prediction at Different Starting Points for the Same Model
In order to verify the generalization ability of the MKRVM, the SOH predictio The RMSE and MAPE of SOH prediction results of the four methods are shown in Figure 8. It can be observed from Figure 8 that M1 shows the lowest RMSE and MAPE, which are 0.004 and 0.293, respectively. The comparisons between different models indicate that the proposed method has the best SOH prediction performance among the four methods. It is clear that the MKRVM is highly accurate and effective for SOH prediction. The error comparisons between M1 and M2 indicate that the Gaussian combined with sigmoid multi-kernel functions have better performance than the Gaussian kernel function, which indicates that the proposed method is valid for SOH prediction. The RMSE of M1 is the smallest, and the RMSE of M2 is the largest. Although the RMSEs of M3 and M4 are less than M2, the M3 and M4 models are not sparse, and the M4 model cannot provide probabilistic prediction results. So, we consider improving the M2. The error of M2 is not satisfying; this is because for M2 without a global sigmoid kernel function, and with the MKRVM, which is the improved RVM, the RMSE is the best. It can be observed that the MKRVM method can improve the performance of SOH prediction.

SOH Prediction at Different Starting Points for the Same Model
In order to verify the generalization ability of the MKRVM, the SOH prediction batteries #5, #6, #7, and #18 at different starting points are tested, whose results are p sented in Figure 9. In numerical experiments, the starting points for batteries #5, #6, #7 are 60, 80, and 100, respectively, while the starting points for batteries #18, #54, and are 40, 60, and 80, respectively. Moreover, the SOH prediction errors are shown in Fig  9b,d,f,h; it is clear that the errors of SOH prediction are very low, which indicates that accuracy of the MKRVM method is high. It can be observed that the proposed method predict the SOH accurately in general. For instance, the prediction of battery #5 is sim to the actual value, which demonstrates that the MKRVM method can effectively ob more accurate prediction results. It is also shown that the predictions at the later t point are close to actual values; the reason is that the earlier the forecast starting point be set up, the more historical information can be tested to build the model. In addit the late period predictions are more important than the early ones for battery health m agement. The SOH prediction in the later period can provide us with guidance on w to replace the battery. However, the prediction errors from different starting points very small, and it can be concluded that the proposed method in this paper shows hig reliability and strong generalization ability, which can be attributed to the multi-ker functions and the extracted HIs.
To verify the effectiveness of the MKRVM method for the SOH prediction of batte at different temperatures, batteries #54 and #55 that operate at low temperature are furt used for SOH prediction. Batteries #54 and #55 have capacity recovery phenomena dur the degradation process and present various local fluctuations. Figure 9 presents S results and errors on the two batteries; it can be observed from Figure 9 that the pre tions of the proposed method are close to the true value, which can capture the local f tuations and especially at the change points. This indicates that the proposed method provide effective SOH prediction when there are many local regenerations in the batt capacity. Moreover, the SOH errors are kept in a certain range, which demonstrates t the MKRVM method is able to improve the prediction accuracy of the proposed meth Table 3 shows the RMSE and MAPE of different starting points. The MAPE of batt #5 is the smallest among the four batteries. The RMSE of different starting points of MKRVM approach is controlled below 0.01, and we obtain the minimum RMSE MAPE, which are 0.0027 and 0.29%, and the maximum MAPE is less than 0.4%. Ot

SOH Prediction at Different Starting Points for the Same Model
In order to verify the generalization ability of the MKRVM, the SOH prediction of batteries #5, #6, #7, and #18 at different starting points are tested, whose results are presented in Figure 9. In numerical experiments, the starting points for batteries #5, #6, and #7 are 60, 80, and 100, respectively, while the starting points for batteries #18, #54, and #55 are 40, 60, and 80, respectively. Moreover, the SOH prediction errors are shown in Figure 9b,d,f,h; it is clear that the errors of SOH prediction are very low, which indicates that the accuracy of the MKRVM method is high. It can be observed that the proposed method can predict the SOH accurately in general. For instance, the prediction of battery #5 is similar to the actual value, which demonstrates that the MKRVM method can effectively obtain more accurate prediction results. It is also shown that the predictions at the later time point are close to actual values; the reason is that the earlier the forecast starting point can be set up, the more historical information can be tested to build the model. In addition, the late period predictions are more important than the early ones for battery health management. The SOH prediction in the later period can provide us with guidance on when to replace the battery. However, the prediction errors from different starting points are very small, and it can be concluded that the proposed method in this paper shows higher reliability and strong generalization ability, which can be attributed to the multi-kernel functions and the extracted HIs.
To verify the effectiveness of the MKRVM method for the SOH prediction of batteries at different temperatures, batteries #54 and #55 that operate at low temperature are further used for SOH prediction. Batteries #54 and #55 have capacity recovery phenomena during the degradation process and present various local fluctuations. Figure 9 presents SOH results and errors on the two batteries; it can be observed from Figure 9 that the predictions of the proposed method are close to the true value, which can capture the local fluctuations and especially at the change points. This indicates that the proposed method can provide effective SOH prediction when there are many local regenerations in the battery capacity. Moreover, the SOH errors are kept in a certain range, which demonstrates that the MKRVM method is able to improve the prediction accuracy of the proposed method. smaller prediction error with the proposed method, which indicates the ability of th multi-kernel functions. In addition, it is also noted that the performance of the extracte HIs are valid. Although batteries #54 and #55 have many local fluctuations, the predictio results can follow this fluctuation, which proves that the proposed method is still effectiv for SOH prediction when the temperature changes or has many local fluctuations.      Table 3 shows the RMSE and MAPE of different starting points. The MAPE of battery #5 is the smallest among the four batteries. The RMSE of different starting points of the MKRVM approach is controlled below 0.01, and we obtain the minimum RMSE and MAPE, which are 0.0027 and 0.29%, and the maximum MAPE is less than 0.4%. Other batteries also follow a similar trend. The largest RMSE and MAPE of the proposed method among four batteries are 0.0115 and 1.457%, respectively. In addition, the prediction results change with the different starting points, and the error range of the method proposed in this paper is relatively small. It can be concluded that the accuracy of the model proposed in this paper is high. Moreover, the RMSE and MAPE results of batteries #54 and #55 are also shown in Table 3. It can be observed that the results of SOH prediction in low temperature have a smaller prediction error with the proposed method, which indicates the ability of the multi-kernel functions. In addition, it is also noted that the performance of the extracted HIs are valid. Although batteries #54 and #55 have many local fluctuations, the prediction results can follow this fluctuation, which proves that the proposed method is still effective for SOH prediction when the temperature changes or has many local fluctuations.

Conclusions
We propose a method to predict the SOH of lithium-ion batteries based on changes in discharge voltage and temperature in this paper. The voltage and temperature extracted from the discharge process can be used for online prediction. Then, the MKRVM model is improved based on the basic RVM model in terms of the kernel function design in order to improve the SOH estimate accuracy, and the PSO algorithm is used to optimize the parameters of the MKRVM. The experimental results of different models of batteries #5, #6, #7, and #18 show that the method proposed in this paper has higher accuracy than the other methods, the experimental results of batteries #5-#18 and batteries #54 and #55 at different starting points show that the method presented in this paper has high accuracy.
The further works can be considered as follows: (1) The study of long-term SOH prediction for other kinds of batteries based on MKRVM due to that long-term SOH prediction is of great significance to the preventive management of batteries; (2) The design of more accurate SOH prediction methods by optimizing the current model; (3) The study of an adaptive kernel function of MKRVM for lithium-ion batteries by reducing the calculation time. (4) We will utilize the random datasets, which are continuously cycled with randomly generated current profiles, to verify the effectiveness of the proposed method in our future work.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.