Battery Capacity Estimation Based on Incremental Capacity Analysis Considering Charging Current Rate

: Incremental capacity analysis (ICA) is widely used in the battery decay mechanism analysis since the features of battery incremental capacity (IC) curves are closely related to battery aging and maximum available capacity. However, the traditional ICA method to estimate battery capacity mainly focuses on a single charging condition, and the inﬂuence of charging current on IC curves is ignored. In this paper, an adaptive capacity estimation method based on ICA considering the charging current is established. First, the charging experiments using different charging current rates under different battery aging statuses are designed and conducted. Then, the relationship between battery maximum available capacity, IC curve features, and charging current is investigated. Furthermore, the ﬁtting method and data-driven method considering charging current are proposed and compared. Finally, the capacity estimation results prove the accuracy and adaptability of the proposed method.


Introduction
Lithium-ion batteries (LIBs) have been widely applied in electric vehicles because of their high energy/power density and low environmental impact [1].However, due to the harsh working conditions and their inherent poor durability, LIBs must be assisted by an advanced battery management system (BMS), which can ensure the safe and reliable operation of LIBs.Battery capacity is an essential parameter in battery management technology and can indicate the battery health status and remaining useful life [2].However, the nonlinear degradation of battery capacity due to the effect of inside and outside factors brings great challenges to capacity estimation.Some effective battery capacity estimation methods have been proposed, including the state of charge (SOC) based method [3,4], the data-driven based method [5,6], and the incremental capacity analysis (ICA) based method [7,8].Bartlett et al. [9] proposed a dual-nonlinear observer to estimate the SOC and capacity based on a reduced-order electrochemical mode, but the battery aging only considered the loss of cyclable lithium.Cheng et al. [10] proposed an adaptive square root unscented Kalman filter method to improve the SOC estimation accuracy based on a second-order equivalent circuit model.However, the effect of aging and temperature was not taken into account in the model.With the development of artificial intelligence, methods using machine learning [5,11] and deep learning [12,13] have drawn great interest and proved their effectiveness in battery capacity or state of health (SOH) estimation.Hu et al. [5] extracted five characteristics related to the capacity from the charge curve and used k-nearest neighbor regression to build a nonlinear kernel regression model.Wu et al. [12] applied control theory to assist the estimation of SOH, in which the SOH is treated as an internal state.Then, the certain system output data and battery SOH are modeled by the group method of a data handling polynomial neural network.Nevertheless, a great amount of experimental data is essential for the data-driven method, which is hard and time-consuming to acquire.Moreover, the quality of experimental data and the speed of the training process have a great influence on the estimation performance.
Facing the defects of the SOC based method and the data-driven based method, the ICA based method received more and more attention.It has been proved as an effective way to estimate the battery capacity.In addition, it is also considered as a good approach to study the battery fading mechanisms, since there is no need to damage the battery's internal structure.IC curves are obtained by calculating the derivative of the charging/discharging capacity to battery voltage.The voltage plateaus caused by the electrode phase transition is then transferred into valleys and peaks on the incremental capacity (IC) curve, which is more intuitive and identifiable.The reduction, shift, or disappearance of the IC curves are employed to reveal the battery degradation and estimate capacity, which represents different electrochemical processes of the LIBs.However, it is also important to avoid the deformation of the IC curves caused by the polarization effect.Thus, a low charging/discharging current with C/25 is usually adopted in the experiment to ensure the equilibrium state.Overall, ICA is a promising method to estimate battery capacity [14].
The foundation of the capacity estimation based on ICA is to determine the features of IC curves, which have a strong relationship with battery capacity fading [15].A universal capacity model based on the charging curve is proposed in Ref. [15], which established the correlation between SOH and the electrode phase transition.The feasibility and accuracy are verified on batteries of different types and aging statuses.In Ref. [14], the influence of sampling frequency on the IC curves was investigated, and the position of IC peaks was adopted to estimate battery capacity after obtaining the IC curves under different battery aging statuses.Furthermore, Weng et al. [16] employed ICA curves to characterize the aging behavior of LIBs, and further extended the SOH monitoring approach from single cells to battery modules or packs that have parallel-connected cells.
However, the traditional ICA method is limited by the charging condition, such as charging current, charging temperature, or initial SOC, which have a severe influence on ICA [8], while previous ICA studies mainly focused on the standard charging condition (charging from 0% SOC with a 0.5 C current at 25 • C, 1 C = 2.75 A).Therefore, it is necessary to study the IC curve characteristics under different charging conditions to improve the adaptability of the ICA method.This paper will mainly focus on the impact of charging current rate on IC curve features, and further propose an adaptive ICA method to estimate battery capacity.In order to achieve the research goal, we have conducted a series of experiments and introduced two different methods to fit the relationship between the capacity and IC curve features.The final results show a great improvement compared with the traditional method.
The specific work is organized as follows: Section 2 introduces the battery experiments and the IC curve characteristics under different charging current rates.Section 3 proposes two different methods to obtain the relationship between battery capacity and picked IC curve characteristics.Then, the feasibility of the two methods is proved.Section 4 further performs a comparison between the accuracy of the two improved methods and the traditional method by qualitative and quantitative analysis.Some conclusions and summarizations are provided in the final section.

Battery Experiments
A kind of commercial NCM/graphite cylindrical battery with a rated capacity of 2750 mAh and a nominal capacity of 2900 mAh is used in this study.The battery charge and discharge cut-off voltage are 4.2 V and 2.5 V, respectively.The recommended charging temperature ranges from 0 • C to 45 • C. The experiment setups include a battery charge and discharge system, a temperature measurement station, an environmental chamber, and a host computer, as shown in the following Figure 1.
World Electr.Veh.J. 2021, 12, x FOR PEER REVIEW 3 of 13 A kind of commercial NCM/graphite cylindrical battery with a rated capacity of 2750 mAh and a nominal capacity of 2900 mAh is used in this study.The battery charge and discharge cut-off voltage are 4.2 V and 2.5 V, respectively.The recommended charging temperature ranges from 0 °C to 45 °C.The experiment setups include a battery charge and discharge system, a temperature measurement station, an environmental chamber, and a host computer, as shown in the following Figure 1.To obtain the relationship between IC features, capacity fading, and charging current rate, some charging experiments are designed and conducted.The specific experimental steps are shown in Table 1.Considering the 0.5 C charging current is relatively large for vehicular application, this study focuses on charging current rates smaller than 0.5 C, and includes the charging experiments of 0.1 C, 0.2 C, 0.3 C, 0.4 C, and 0.5 C. It should be noted that the charging experiment is performed periodically during the battery accelerated aging test.To obtain the relationship between IC features, capacity fading, and charging current rate, some charging experiments are designed and conducted.The specific experimental steps are shown in Table 1.Considering the 0.5 C charging current is relatively large for vehicular application, this study focuses on charging current rates smaller than 0.5 C, and includes the charging experiments of 0.1 C, 0.2 C, 0.3 C, 0.4 C, and 0.5 C. It should be noted that the charging experiment is performed periodically during the battery accelerated aging test.Traditionally, the numerical derivative method is adopted to acquire IC curves.According to the definition of ICA, the charging capacity divided by terminal voltage change within the equal voltage interval (EVI) or equal time interval (ETI) is the IC value.The acquisition of IC curves through the numerical derivative method can adopt the following equation.
EVI : However, the IC curve will be severely affected by the selected voltage or time interval.Therefore, it is necessary to adopt a proper filter to obtain a smooth IC curve.In our previous study [8], a fast acquisition method for IC curves is proposed, and in this paper, we still use the Kalman filter (KF) method to smooth IC curves.
Based on the KF method mentioned above and the experimental data obtained in Section 2.1.,IC curves with the different charging current rates at different aging cycles are obtained, as shown in Figure 2.
However, the IC curve will be severely affected by the selected voltage or time interval.Therefore, it is necessary to adopt a proper filter to obtain a smooth IC curve.In our previous study [8], a fast acquisition method for IC curves is proposed, and in this paper, we still use the Kalman filter (KF) method to smooth IC curves.
Based on the KF method mentioned above and the experimental data obtained in Section 2.1.,IC curves with the different charging current rates at different aging cycles are obtained, as shown in Figure 2. Overall, the shape of the IC curves with the same charging current is basically the same under different aging cycles.However, the location and intensity of curves change under different charging current rates.With the decrease in the charging current, the IC curves have a moving trend to low voltage.For cycle 0, every peak or valley moves about 0.1 V left as the charging current decreases from 0.5 C to 0.1 C. Additionally, the fourth peak is shown around 4.15 V when the current is less than 0.3 C.Moreover, when the battery SOH is high and the charging current is low, as shown in Figure 2a, b, a smaller peak is generated at 3.6 V. Therefore, the charging/discharging process with extremely low current is essential for the analysis of the battery degradation mechanism in detail.In addition, it can be seen that with the decrease in the charging current, the height of all peaks increases, and the height of the valleys decreases.Taking the second peak around 3.7 V as an example, for 0, 200, 400, and 600 cycles, the height of the peak increases by 0.43, 0.67, 0.92, and 1.04, respectively, when the current decreases from 0.5 C to 0.1 C. It is also noteworthy that, the more serious the aging status of the battery, the greater the IC curves of the battery affected by the charging current.As the SOH decreases, the fluctuation of the curves becomes flat.When the highest peak of IC curves under 0.1 C is 6.5 Ah/V for cycle 0, it only reaches 4.8 Ah/V for cycle 600.Thus, taking the coupling effect of aging and charging current to IC curve characteristics into consideration is helpful for expanding the adaptability of the ICA method in capacity estimation.Overall, the shape of the IC curves with the same charging current is basically the same under different aging cycles.However, the location and intensity of curves change under different charging current rates.With the decrease in the charging current, the IC curves have a moving trend to low voltage.For cycle 0, every peak or valley moves about 0.1 V left as the charging current decreases from 0.5 C to 0.1 C. Additionally, the fourth peak is shown around 4.15 V when the current is less than 0.3 C.Moreover, when the battery SOH is high and the charging current is low, as shown in Figure 2a, b, a smaller peak is generated at 3.6 V. Therefore, the charging/discharging process with extremely low current is essential for the analysis of the battery degradation mechanism in detail.In addition, it can be seen that with the decrease in the charging current, the height of all peaks increases, and the height of the valleys decreases.Taking the second peak around 3.7 V as an example, for 0, 200, 400, and 600 cycles, the height of the peak increases by 0.43, 0.67, 0.92, and 1.04, respectively, when the current decreases from 0.5 C to 0.1 C. It is also noteworthy that, the more serious the aging status of the battery, the greater the IC curves of the battery affected by the charging current.As the SOH decreases, the fluctuation of the curves becomes flat.When the highest peak of IC curves under 0.1 C is 6.5 Ah/V for cycle 0, it only reaches 4.8 Ah/V for cycle 600.Thus, taking the coupling effect of aging and charging current to IC curve characteristics into consideration is helpful for expanding the adaptability of the ICA method in capacity estimation.

Capacity Estimation Based on the Fitting Method
Based on the analysis in Section 2.2, the relationship between IC curve characteristics and capacity under different current rates can be fitted for capacity estimation.In this paper, we adopt the height of the IC peaks to estimate the battery capacity.The peaks are associated with the phase transition process in the positive and negative electrodes, as discussed in Refs.[17][18][19].Basically, capacity fading is mainly caused by the loss of lithium inventory (LLI), the loss of active material (LAM), and resistance increase during the phase transition process.However, which aging factors dominate depends on the specific using conditions of the battery.The area under each peak represents the charging capacity during the corresponding phase change process.The decrease in the peak value means LAM of the negative electrode.If a certain peak loses more proportion than the others, it indicates that LLI also happens.Additionally, the moving trend towards high voltage means that the plateau voltage corresponding to each phase transition process has a significant increase, indicating that the internal resistance has increased.Here, because peak II around 3.7 V and peak III around 3.9 V have a strong correlation with the battery capacity under all current rates, as shown in Figure 3, these two features will be employed for the ICA based estimation method.As analyzed before, it can be found that LAM and resistance increase happened during both the phase transition process of peak II and peak III.In addition, LLI also happened in peak II after cycle 400.

Capacity Estimation Based on the Fitting Method
Based on the analysis in Section 2.2., the relationship between IC curve characteristics and capacity under different current rates can be fitted for capacity estimation.In this paper, we adopt the height of the IC peaks to estimate the battery capacity.The peaks are associated with the phase transition process in the positive and negative electrodes, as discussed in Refs.[17][18][19].Basically, capacity fading is mainly caused by the loss of lithium inventory (LLI), the loss of active material (LAM), and resistance increase during the phase transition process.However, which aging factors dominate depends on the specific using conditions of the battery.The area under each peak represents the charging capacity during the corresponding phase change process.The decrease in the peak value means LAM of the negative electrode.If a certain peak loses more proportion than the others, it indicates that LLI also happens.Additionally, the moving trend towards high voltage means that the plateau voltage corresponding to each phase transition process has a significant increase, indicating that the internal resistance has increased.Here, because peak II around 3.7 V and peak Ⅲ around 3.9 V have a strong correlation with the battery capacity under all current rates, as shown in Figure 3, these two features will be employed for the ICA based estimation method.As analyzed before, it can be found that LAM and resistance increase happened during both the phase transition process of peak II and peak III.In addition, LLI also happened in peak II after cycle 400.The linear relationship between the battery capacity and the features can be established with the following equation: where i represents the ith IC feature,  represents the estimated capacity according to the ith feature. represents the identified height of the ith feature.The corresponding coefficients  and  with different currents are listed in Table 2.The linear relationship between the battery capacity and the features can be established with the following equation: where i represents the ith IC feature, C F i represents the estimated capacity according to the ith feature.H F i represents the identified height of the ith feature.The corresponding coefficients α i and β i with different currents are listed in Table 2. From Table 2, we can observe that almost all R 2 of the coefficients of peak II and peak III are greater than 0.95, which indicates the high correlation between the coefficients and capacity.In addition, it can be found that for peak II, the value of α 1 gradually drops as the charging current increases, while that of β 1 shows an increasing trend.The coefficients of peak III show a complete opposite change rule.From Figure 4 below, we can have a more intuitive understanding of the regular.
From Table 2, we can observe that almost all R 2 of the coefficients of peak Ⅱ and peak Ⅲ are greater than 0.95, which indicates the high correlation between the coefficients and capacity.In addition, it can be found that for peak Ⅱ, the value of  gradually drops as the charging current increases, while that of  shows an increasing trend.The coefficients of peak Ⅲ show a complete opposite change rule.From Figure 4 below, we can have a more intuitive understanding of the regular.The relationship can be fitted by nonlinear functions shown in the Equations ( 3 Finally, we can obtain the estimated capacity with the following steps: first, two groups of coefficients  and  corresponding to two characteristics under different charging current rates are calculated using Equations ( 3) and ( 4).Secondly, the estimated capacity of two characteristics under a specific charging current can be acquired by substituting the coefficients into Equation (2).Lastly, the average of two results will be considered as the final estimated capacity.

Capacity Estimation Based on the Data-Driven Method
In the previous method, there are two error sources: one comes from the fitting process of the linear relationship between the battery capacity and the features, the other comes from the fitting process of the coefficients and the charging current.In order to further improve the precision of the estimation, we adopted a data-driven Gaussian process regression (GPR) [20] method to acquire the relationship among the capacity, IC curve features, and charging current directly.GPR is a machine learning regression method based on statistical mathematics theory.It is a probabilistic machine learning method that can avoid local optima during parameter optimization, and can also solve nonlinear or high-dimension problems.The procedure of the GPR method is explained as follows: First, consider the training set where xi and yi represent the input and output data.
Pretend that there is a function f (•), which maps the input value xn to output value yn, as written in Equation ( 6): The relationship can be fitted by nonlinear functions shown in the Equations ( 3) and ( 4). (3) Finally, we can obtain the estimated capacity with the following steps: first, two groups of coefficients α i and β i corresponding to two characteristics under different charging current rates are calculated using Equations ( 3) and ( 4).Secondly, the estimated capacity of two characteristics under a specific charging current can be acquired by substituting the coefficients into Equation (2).Lastly, the average of two results will be considered as the final estimated capacity.

Capacity Estimation Based on the Data-Driven Method
In the previous method, there are two error sources: one comes from the fitting process of the linear relationship between the battery capacity and the features, the other comes from the fitting process of the coefficients and the charging current.In order to further improve the precision of the estimation, we adopted a data-driven Gaussian process regression (GPR) [20] method to acquire the relationship among the capacity, IC curve features, and charging current directly.GPR is a machine learning regression method based on statistical mathematics theory.It is a probabilistic machine learning method that can avoid local optima during parameter optimization, and can also solve nonlinear or high-dimension problems.The procedure of the GPR method is explained as follows: First, consider the training set where x i and y i represent the input and output data.
Pretend that there is a function f (•), which maps the input value x n to output value y n , as written in Equation ( 6): where the value of noise ξ n follows the Gaussian process ξ n ~N(0, σ n 2 ).Here, we use squared exponential kernel functions as the kernel of the GPR model, as written in Equation (7): where l = [l 0 , l 1 , . . .l d ] denotes hyperparameter in the kernel function.Substituting Equation ( 7) into Equation ( 6), we can obtain Equation (8).
Therefore, the distribution of output y follows: Using maximum likelihood estimation to find the hyperparameters of the kernel function of the formula, we finally determine the kernel function: Therefore, the combined distribution of the training output and the test output can be expressed as: where the mean and covariance of the predicted distribution are as follows: The mean value of the predicted distribution in Equation ( 11) is actually the estimation value of the output.
Here, the training inputs of the GPR model are the charging current and the heights of peak II and peak III, and capacity obtained from the experiment is used as the test output, as shown in Figure 5.The training output showed a high precision with a training error under 0.8% in most conditions.Only in the situation of high aging cycles under 0.1 C is the error relatively high, ranging from 0.8% to 1.4%.Thus, it is feasible to use the model for the capacity estimation.
,    ,    , (8) Therefore, the distribution of output y follows: |,   , Using maximum likelihood estimation to find the hyperparameters of the kernel function of the formula, we finally determine the kernel function: Therefore, the combined distribution of the training output and the test output can be expressed as: where the mean and covariance of the predicted distribution are as follows: The mean value of the predicted distribution in Equation ( 11) is actually the estimation value of the output.
Here, the training inputs of the GPR model are the charging current and the heights of peak II and peak III, and capacity obtained from the experiment is used as the test output, as shown in Figure 5.The training output showed a high precision with a training error under 0.8% in most conditions.Only in the situation of high aging cycles under 0.1 C is the error relatively high, ranging from 0.8% to 1.4%.Thus, it is feasible to use the model for the capacity estimation.

Comparison of the Estimation Accuracy between the Fitting and Data-Driven Methods
In order to validate the effectiveness and further compare the accuracy of the two proposed methods, we also performed the same experiments on another battery.The estimation results are shown in Figure 6, and we can figure out that the capacity estimation error of the data-driven method is much lower than the fitting method in most conditions, which indicates better estimation accuracy.The maximum estimation error difference between the two methods is up to 2%.Additionally, it can be observed that the estimation errors of the two methods both increase with the battery aging in general.From the aspect of charging current, the error of the two methods shows a downward trend with increasing of the current rate.It is noteworthy that the GPR method has more stable performance in the whole lifecycle.For example, the error of the fitting method fluctuates from 2.0% to

Comparison of the Estimation Accuracy between the Fitting and Data-Driven Methods
In order to validate the effectiveness and further compare the accuracy of the two proposed methods, we also performed the same experiments on another battery.The estimation results are shown in Figure 6, and we can figure out that the capacity estimation error of the data-driven method is much lower than the fitting method in most conditions, which indicates better estimation accuracy.The maximum estimation error difference between the two methods is up to 2%.Additionally, it can be observed that the estimation errors of the two methods both increase with the battery aging in general.From the aspect of charging current, the error of the two methods shows a downward trend with increasing of the current rate.It is noteworthy that the GPR method has more stable performance in the whole lifecycle.For example, the error of the fitting method fluctuates from 2.0% to 0.09%, while that of the data-driven method is much lower, ranging from 0.7% to 0.2% under a 0.3 C charging current.
In order to have a quantitative analysis of the two methods, we further observed average error and the root mean squared error (RMSE) of the two methods during the whole lifecycle, as shown in Table 3.We can tell that both the average error and RMSE of the data-driven method are better than that of the fitting method.While the average error of the fitting method remains above 1% and achieves 3.9% under 0.1 C, the data-driven method shows the improvement of precision up to 1%.Similar results can be obtained in terms of RMSE.Thus, the data-driven method is effective in the corrected ICA method.In order to have a quantitative analysis of the two methods, we further observed average error and the root mean squared error (RMSE) of the two methods during the whole lifecycle, as shown in Table 3.We can tell that both the average error and RMSE of the data-driven method are better than that of the fitting method.While the average error of the fitting method remains above 1% and achieves 3.9% under 0.1 C, the data-driven method shows the improvement of precision up to 1%.Similar results can be obtained in terms of RMSE.Thus, the data-driven method is effective in the corrected ICA method.

Comparison of the Accuracy between the GPR Method and Traditional Method
Figure 7 shows an intuitive comparison of the data-driven method and traditional method in the whole life cycle of the battery.Here, the traditional method employs the data under 0.5 C in Table 2 to estimate capacity without considering the charging current.It is confirmed that the proposed method is effective.Compared with the traditional method, the precision of the improved method considering charging current has an improvement of up to 15.175%.In addition, the smaller the charging current is, the bigger the improvement will be.For the method ignoring charging current, the error ranges from around 19.5% under 0.1 C to 0.39% under 0.5 C.Although the accuracy of the proposed method also reduces slightly with the decrease in charging current, the variation of the error is still within an acceptable range, not more than 5.6%.It is worth noting that the estimation accuracy no longer has a linear relationship with the charging current compared with the traditional method, which means the robustness to current interference has improved.During the whole life cycle, the estimation precision also remains at a high level; thus, the proposed method also has good adaptability to capacity attenuation.
World Electr.Veh.J. 2021, 12, x FOR PEER REVIEW 9 of 13 Figure 7 shows an intuitive comparison of the data-driven method and traditional method in the whole life cycle of the battery.Here, the traditional method employs the data under 0.5 C in Table 2 to estimate capacity without considering the charging current.It is confirmed that the proposed method is effective.Compared with the traditional method, the precision of the improved method considering charging current has an improvement of up to 15.175%.In addition, the smaller the charging current is, the bigger the improvement will be.For the method ignoring charging current, the error ranges from around 19.5% under 0.1 C to 0.39% under 0.5 C.Although the accuracy of the proposed method also reduces slightly with the decrease in charging current, the variation of the error is still within an acceptable range, not more than 5.6%.It is worth noting that the estimation accuracy no longer has a linear relationship with the charging current compared with the traditional method, which means the robustness to current interference has improved.During the whole life cycle, the estimation precision also remains at a high level; thus, the proposed method also has good adaptability to capacity attenuation.

The Influence of Charging Current Offset on Capacity Estimation
For both of the two methods mentioned above, the charging current is used as a vital input.However, the measurement error may exist in the actual test conditions and could

The Influence of Charging Current Offset on Capacity Estimation
For both of the two methods mentioned above, the charging current is used as a vital input.However, the measurement error may exist in the actual test conditions and could lead to the deviation of the estimation results.For the fitting method, different charging current rates correspond to different IC curves, which means the characteristics and corresponding coefficients will also change with it.For the data-driven method, the deviation of the charging current will also affect the estimation error.Thus, the robustness ability against the interference of charging current is of great importance.
In order to validate the robustness of the two methods, interference verification was conducted.Here, we impose a ±10 mA offset on the actual charging current and investigate the effects on the capacity estimation error of the two methods.Figure 8 demonstrates the interference error of the fitting method caused by the current offset, which indicates the increased deviation from the original estimation results shown in Figure 6.
against the interference of charging current is of great importance.
In order to validate the robustness of the two methods, interference verification was conducted.Here, we impose a ± 10 mA offset on the actual charging current and investigate the effects on the capacity estimation error of the two methods.Figure 8 demonstrates the interference error of the fitting method caused by the current offset, which indicates the increased deviation from the original estimation results shown in Figure 6.
It is noteworthy that the interference error is falling with the increase in the charging current.In the high charging current conditions from 0.3 C to 0.5 C, the interference error maintains in an acceptable range within 1%.However, the error under low charging current is relatively high.Additionally, while the error has a downward trend with the battery aging under the current rate from 0.1 C to 0.3 C, an opposite trend is observed under 0.4 C and 0.5 C. Furthermore, the interference error of −10 mA offset is always lower than +10 mA offset, and the gap narrows with increasing of the charging current rate.It is noteworthy that the interference error is falling with the increase in the charging current.In the high charging current conditions from 0.3 C to 0.5 C, the interference error maintains in an acceptable range within 1%.However, the error under low charging current is relatively high.Additionally, while the error has a downward trend with the battery aging under the current rate from 0.1 C to 0.3 C, an opposite trend is observed under 0.4 C and 0.5 C. Furthermore, the interference error of −10 mA offset is always lower than +10 mA offset, and the gap narrows with increasing of the charging current rate.
Figure 9 shows the interference error of the data-driven method caused by the current offset.Compared with the fitting method, this method performs much better in resistance to the charging current deviation.First, the interference error remains steady under different charging current rates, ranging from 0.02% to 0.09%, indicating better robustness under different charging conditions.The interference error of both offsets overlaps under most conditions.Thus, the data-driven method also performs more stably under different charging current offsets.Overall, the data-driven method has a better and more stable performance than the fitting method.
Figure 9 shows the interference error of the data-driven method caused by the current offset.Compared with the fitting method, this method performs much better in resistance to the charging current deviation.First, the interference error remains steady under different charging current rates, ranging from 0.02% to 0.09%, indicating better robustness under different charging conditions.The interference error of both offsets overlaps under most conditions.Thus, the data-driven method also performs more stably under different charging current offsets.Overall, the data-driven method has a better and more stable performance than the fitting method.

Conclusions
This paper investigates the capacity estimation using the ICA method considering the charging current rate.First, the characteristics of IC curves with different charging current rates under different aging statuses are investigated.Then, relationships between the charging current and the coefficients of different features are established to construct the fitting method.Furthermore, a data-driven GPR method is introduced to reduce the estimation error caused by the two-step fitting process of the first method.In order to validate the robustness against the measurement error of charging current, the interference test also is performed on two methods.Finally, the accuracy and adaptability of different methods are validated and compared.

Conclusions
This paper investigates the capacity estimation using the ICA method considering the charging current rate.First, the characteristics of IC curves with different charging current rates under different aging statuses are investigated.Then, relationships between the charging current and the coefficients of different features are established to construct the fitting method.Furthermore, a data-driven GPR method is introduced to reduce the estimation error caused by the two-step fitting process of the first method.In order to validate the robustness against the measurement error of charging current, the interference test also is performed on two methods.Finally, the accuracy and adaptability of different methods are validated and compared.
Some conclusions and summarizations are drawn below: • Charging current has a significant impact on IC curves.With the decrease in the charging current, the IC curves have a moving trend to low voltage, and the height of all peaks increases, and the height of valleys decreases.In addition, the more serious the aging status of the battery, the more the IC curves of the battery will be affected by the charging current.

Figure 3 .
Figure 3.The correlation between the features of IC curves and battery capacity: (a) peak II; (b) peak III.

Figure 3 .
Figure 3.The correlation between the features of IC curves and battery capacity: (a) peak II; (b) peak III.

Figure 4 .
Figure 4.The correlation between coefficients and current rate: (a) peak II; (b) peak III.

Figure 4 .
Figure 4.The correlation between coefficients and current rate: (a) peak II; (b) peak III.

Figure 5 .
Figure 5. Capacity estimation mode based on GPR.

Figure 5 .
Figure 5. Capacity estimation mode based on GPR.

Figure 6 .
Figure 6.Accuracy comparison of two methods considering charging current.

Figure 6 .
Figure 6.Accuracy comparison of two methods considering charging current.

Figure 7 .
Figure 7. Accuracy comparison of the traditional method and data-driven method.

Figure 7 .
Figure 7. Accuracy comparison of the traditional method and data-driven method.

Figure 8 .
Figure 8. Interference error generated by ± 10 mA offset of the fitting method.Figure 8. Interference error generated by ±10 mA offset of the fitting method.

Figure 8 .
Figure 8. Interference error generated by ± 10 mA offset of the fitting method.Figure 8. Interference error generated by ±10 mA offset of the fitting method.

Figure 9 .
Figure 9. Interference error generated by ± 10 mA offset of the data-driven method.

Figure 9 .
Figure 9. Interference error generated by ±10 mA offset of the data-driven method.

Table 1 .
Charging experiments at different charging current rates.

Table 1 .
Charging experiments at different charging current rates.

Table 2 .
The coefficients  and  with different currents.

Table 2 .
The coefficients α i and β i with different currents.

Table 3 .
Comparison of the average error and RMSE of the two methods.

Table 3 .
Comparison of the average error and RMSE of the two methods.