A Fast Online State of Health Estimation Method for Lithium-Ion Batteries Based on Incremental Capacity Analysis

: E ﬃ cient and accurate state of health (SoH) estimation is an important challenge for safe and e ﬃ cient management of batteries. This paper proposes a fast and e ﬃ cient online estimation method for lithium-ion batteries based on incremental capacity analysis (ICA), which can estimate SoH through the relationship between SoH and capacity di ﬀ erentiation over voltage ( dQ / dV ) at di ﬀ erent states of charge (SoC). This method estimates SoH using arbitrary dQ / dV over a large range of charging processes, rather than just one or a limited number of incremental capacity peaks, and reduces the SoH estimation time greatly. Speciﬁcally, this method establishes a black box model based on ﬁtting curves ﬁrst, which has a smaller amount of calculation. Then, this paper analyzes the inﬂuence of di ﬀ erent SoC ranges to obtain reasonable ﬁtting curves. Additionally, the selection of a reasonable dV is taken into account to balance the e ﬃ ciency and accuracy of the SoH estimation. Finally, experimental results validate the feasibility and accuracy of the method. The SoH estimation error is within 5% and the mean absolute error is 1.08%. The estimation time of this method is less than six minutes. Compared to traditional methods, this method is easier to obtain e ﬀ ective calculation samples and saves computation time.


Introduction
Developing electric vehicles (EVs) is an important measure to solve energy and environmental problems [1][2][3]. As the power source for electric vehicles, lithium-ion batteries directly affect the performance of EVs [4]. Accurate estimations of battery status, such as the state of charge (SoC), state of power (SoP), and state of health (SoH), are critical to the safe and efficient operation of EVs [5]. Among them, SoH is an important indicator to characterize battery degradation trends and availability. The fast and accurate estimation of SoH is an important challenge for efficient and safe management of batteries [6].
At present, there are three types of SoH estimation methods: the model-based method, the data-driven method, and the incremental capacity analysis (ICA) method. In the model-based method, SoH is usually calculated through the estimated internal resistance of the battery using an equivalent circuit model [7]. To achieve high estimation accuracy, a highly accurate model of the battery is needed. Moreover, the inconsistency of the battery's environment and performance leads to its poor universality [8]. In the data-driven method, computational burden is usually caused [9,10], which reduces the estimated efficiency of the SoH and causes a higher request for CPU performance. ICA is another popular SoH estimation method, which shows a strong correlation with battery aging [11]. The concept of ICA originates from the study of the lithium intercalation process and the corresponding staging phenomenon [12][13][14]. When lithium-ion battery is charged at 0.04 C, its incremental capacity (IC) curve usually has 3-5 peaks, which can be used to estimate SoH [15]. Here, 1 C means that the current is numerically equal to the nominal capacity. For traditional ICA methods, IC peaks are the primary premise to estimate the aging state of battery. However, there are only a limited number of peaks that can be used to estimate SoH [16][17][18][19], that is, the sample range available for SoH estimation is too narrow, and there is no guarantee that the voltage corresponding to the IC peak will be passed in each charging process. Moreover, since IC peaks are the extreme values of the incremental capacity (dQ) during the entire battery charging process, the time to calculate dQ is also the longest, accordingly. For example, if a voltage peak is 3.35 V, experimental results show that when dV = 0.1 V, it takes more than half an hour to calculate the dQ at 3.35 V. Additionally, if the charging time is taken into account, the estimated time of SoH can exceed one hour.
There are two reasons for difficulties in the estimated sample acquisition and the long estimation time of traditional ICA methods, which are the limited number of IC peaks and the excessive capacity calculation time at IC peaks. If an arbitrary capacity differentiation over voltage (dQ/dV) over a large range of charging processes can be used to estimate SoH instead of a few IC peaks, the efficiency and real-time performance of SoH estimation will be improved significantly. Therefore, a SoH estimation method using dQ/dV at different SoCs is proposed. Different from traditional ICA methods, the selection of dQ/dV used to estimate SoH in the proposed method is determined by SoC instead of IC peaks. In this method, the relationship between dQ/dV and SoH under different SoCs is first drawn through experiments. After that, the SoC of the sampling point is determined by calculation. The correspondence between dQ/dV and SoH under this SoC is found by the previously obtained fitting curves. Finally, the corresponding SoH can be estimated by calculating the dQ/dV at the sampling point. The estimated error is limited to 5% when the SoC is within the range of 0.3-0.8, and any dQ/dV in this range can be used to estimate the SoH, which overcomes the disadvantage of traditional ICA methods.

Experiment Setup
The experiments are performed on a Lishen LR1865EH lithium-ion battery with related battery parameters listed in Table 1. And the battery is produced by Tianjin Lishen Battery Co., Ltd., China. As shown in Figure 1, the test bench is NBT5V10AC16-T battery cycler produced by Ningbo Bayte Measurement and Control Technology Co., Ltd., and it consists of a lithium-ion battery, a programmable thermal chamber, and a host computer for monitoring and processing experimental data. The battery cycler has 16 dependent measurement channels with the voltage limits of 0-5 V and current limits of 0.3-10 A. The errors of the current and voltage sensors are both within 0.1%. incremental capacity (IC) curve usually has 3-5 peaks, which can be used to estimate SoH [15]. Here, 1 C means that the current is numerically equal to the nominal capacity. For traditional ICA methods, IC peaks are the primary premise to estimate the aging state of battery. However, there are only a limited number of peaks that can be used to estimate SoH [16][17][18][19], that is, the sample range available for SoH estimation is too narrow, and there is no guarantee that the voltage corresponding to the IC peak will be passed in each charging process. Moreover, since IC peaks are the extreme values of the incremental capacity (dQ) during the entire battery charging process, the time to calculate dQ is also the longest, accordingly. For example, if a voltage peak is 3.35 V, experimental results show that when dV = 0.1 V, it takes more than half an hour to calculate the dQ at 3.35 V. Additionally, if the charging time is taken into account, the estimated time of SoH can exceed one hour.
There are two reasons for difficulties in the estimated sample acquisition and the long estimation time of traditional ICA methods, which are the limited number of IC peaks and the excessive capacity calculation time at IC peaks. If an arbitrary capacity differentiation over voltage (dQ/dV) over a large range of charging processes can be used to estimate SoH instead of a few IC peaks, the efficiency and real-time performance of SoH estimation will be improved significantly. Therefore, a SoH estimation method using dQ/dV at different SoCs is proposed. Different from traditional ICA methods, the selection of dQ/dV used to estimate SoH in the proposed method is determined by SoC instead of IC peaks. In this method, the relationship between dQ/dV and SoH under different SoCs is first drawn through experiments. After that, the SoC of the sampling point is determined by calculation. The correspondence between dQ/dV and SoH under this SoC is found by the previously obtained fitting curves. Finally, the corresponding SoH can be estimated by calculating the dQ/dV at the sampling point. The estimated error is limited to 5% when the SoC is within the range of 0.3-0.8, and any dQ/dV in this range can be used to estimate the SoH, which overcomes the disadvantage of traditional ICA methods.

Experiment Setup
The experiments are performed on a Lishen LR1865EH lithium-ion battery with related battery parameters listed in Table 1. And the battery is produced by Tianjin Lishen Battery Co., Ltd., China. As shown in Figure 1, the test bench is NBT5V10AC16-T battery cycler produced by Ningbo Bayte Measurement and Control Technology Co., Ltd., and it consists of a lithium-ion battery, a programmable thermal chamber, and a host computer for monitoring and processing experimental data. The battery cycler has 16 dependent measurement channels with the voltage limits of 0-5 V and current limits of 0.3-10 A. The errors of the current and voltage sensors are both within 0.1%.   In actual operation, the battery is retired when the capacity declines to 80%. The relationship between SoH defined by the capacity and dQ/dV at different SoCs is shown in Figure 2. It can be seen that dQ/dV has a satisfactory monotonous relationship with SoH when the SoH is greater than 0.75, which meets the SoH range of actual operation on EVs completely. Therefore, the data of 0-400 cycles is used in the following experiments, and the SoH of 400 cycle is about 0.78.  In actual operation, the battery is retired when the capacity declines to 80%. The relationship between SoH defined by the capacity and dQ/dV at different SoCs is shown in Figure 2. It can be seen that dQ/dV has a satisfactory monotonous relationship with SoH when the SoH is greater than 0.75, which meets the SoH range of actual operation on EVs completely. Therefore, the data of 0-400 cycles is used in the following experiments, and the SoH of 400 cycle is about 0.78.

Relationship between the dQ/dV and SoH Experiment
A monotonic relationship model between arbitrary dQ/dV and SoH is established by constraining SoC in this method. dQ/dV is defined as:

Relationship between the dQ/dV and SoH Experiment
A monotonic relationship model between arbitrary dQ/dV and SoH is established by constraining SoC in this method. dQ/dV is defined as: where Q(v) represents the capacity Q at voltage v. ∆V is the voltage interval for calculating dQ and it is a fixed value. v 0 represents the voltage at sampling point. v x is a small offset to ensure the overlap of adjacent sampling points, which improves the credibility of each sampling point. The SoH is defined by capacity as follows: where Q 0 represents the total capacity of the new battery and Q 2 represents the total capacity of the current battery.
The definition of SoC is as follows: where Q 1 represents the capacity charging to the sampling point. Based on the above definition, the experiment is carried out as follows: Step 1: For one entire charging process, current values are count up as the total capacity Q 2 , Equation (2) is used to calculate SoH of the charging process. The SoC point is sampled every 0.1 within the interval of 0.1 to 1.
Step 2: For the charging process selected in step 1, the capacity of SoC 0 is calculated in the following equation: The point where the capacity reaches Q 1 is recorded as i. The voltage value is recorded as V i . After that, the capacity between the interval of [V i − dV/2, V i + dV/2] can be accumulated as dQ/dV.
Step 3: Cycle steps 1 and 2 until 400 charging cycles are all processed. After that, the relationship between dQ/dV at SoC 0 and SoH can be obtained.
Step 4: Next, the SoC sampling point is selected in the interval of 0-1 and steps 1, 2, and 3 are executed.
After the above four steps, the relationship between dQ/dV at 10 different SoCs and SoH can be obtained. As can be seen in Figure 3, green points representing SoC = 1 tends to a line parallel to the SoH axis, that is, the fitted curve of SoC = 1 can hardly be applied to the estimation of the SoH, while in the ten kinds of scatter plots, other scatter plots have satisfactory monotonicity, which can prove that the method is feasible over a large charging range.
Moreover, dQ and dV are usually represented by ∆Q and ∆V at the same time interval [11] or the same percent depth of discharge (% DOD) interval [20] in traditional ICA methods. Due to the differential property, dQ/dV is very sensitive to noise. Traditional ICA methods are sensitive to noise and a small measurement error may cause an obvious SoH estimation error. Therefore, traditional ICA methods generally require filtering, while, in this method, the interval of sampling points is determined by SoC instead of the time interval or % DOD interval. Most SoC points used in this paper are located in the voltage platform region where dV changes inconspicuously, relatively, and the influence of the voltage measurement error is also relatively reduced. Corresponding to the selection of sampling points, the new dQ/dV definition is also used in Equation (1). Only dQ is a quantity that varies with SoC. Therefore, the influence of dV measurement errors can be ignored. As can be seen in Figure 3, in the absence of any filtering denoising measure, there are no significant sharp points, that is, the method is relatively insensitive to noise and has satisfactory robustness.
It should be pointed out that in the process of estimating the SoH by the fitting curves in the later stage, in order to obtain the most accurate estimation effect with as few fitting curves as possible, five fitting curves are appropriately increased and some unnecessary curves are deleted. The estimated effect of this method is satisfactory in the SoC interval of 0.3-0.8. Details can be found in Section 4.1.

SoH Estimation Based on a Black Box Model
Through the above experiments and analysis, multiple fitting curves between dQ/dV at different SoCs and SoH can be obtained. Based on the fitting curves, the black box model of the SoH estimation method is established. In this model, the input is dQ/dV and the output is the estimated SoH, as shown in Figure 4. The SoH estimation method is introduced in detail below. Eleven fitting curves of dQ/dV at different SoC (SoC1, SoC2, …, SoC11) and SoH are selected in this method, as shown in Table 2. The estimated SoH is SoH2 and the actual SoH is SoH0. Assuming that the SoC at the sampling point is SoCtemp, as can be seen in Table 2. SoCmin and SoCmax are defined as follows:

SoH Estimation Based on a Black Box Model
Through the above experiments and analysis, multiple fitting curves between dQ/dV at different SoCs and SoH can be obtained. Based on the fitting curves, the black box model of the SoH estimation method is established. In this model, the input is dQ/dV and the output is the estimated SoH, as shown in Figure 4. The SoH estimation method is introduced in detail below.

SoH Estimation Based on a Black Box Model
Through the above experiments and analysis, multiple fitting curves between dQ/dV at different SoCs and SoH can be obtained. Based on the fitting curves, the black box model of the SoH estimation method is established. In this model, the input is dQ/dV and the output is the estimated SoH, as shown in Figure 4. The SoH estimation method is introduced in detail below. Eleven fitting curves of dQ/dV at different SoC (SoC1, SoC2, …, SoC11) and SoH are selected in this method, as shown in Table 2. The estimated SoH is SoH2 and the actual SoH is SoH0. Assuming that the SoC at the sampling point is SoCtemp, as can be seen in Table 2  Eleven fitting curves of dQ/dV at different SoC (SoC 1 , SoC 2 , . . . , SoC 11 ) and SoH are selected in this method, as shown in Table 2. The estimated SoH is SoH 2 and the actual SoH is SoH 0 . Assuming that the SoC at the sampling point is SoC temp , as can be seen in Table 2. SoC min and SoC max are defined as follows:  The corresponding functions y 1 = f 1 (x) and y 2 = f 2 (x) of SoC max and SoC min can be obtained. Here, x is the value of dQ/dV, and f (x) is the calculated SoH value.
In practical application, SoH estimation follows the equation below: As mentioned above, the SoH is about 0.78 at the 400 th cycle. The SoH value is not allowed to exceed 1. Therefore, the following rules can be set for SoH 2 . After the above treatment, the estimated SoH can be obtained.

Influence Analysis of Different SoC Ranges
To select sufficient and appropriate SoC, the figure of SoC, dQ/dV, and SoH is constructed. As shown in Figure 5, when SoC is in the interval of 0.45-0.75, the overall change rate of dQ/dV is not large. Therefore, fewer SoC points are needed in the interval 0.45-0.75. On the contrary, because the three-dimensional surface diagram in the two intervals of 0-0.45 and 0.75-1 are steep, so more SoC points are needed. Moreover, in Figure 6, dV = 0.005 V is selected to amplify the influence of noise, so that the SoC range applicable to the method can be observed more intuitively. In Figure 6, when SoC is 0.1, 0.2, 0.9, and 1, respectively, dQ/dV and SoH are difficult to form a satisfactory linear relationship. Therefore, in the final SoH estimation method, the final selected SoC interval is 0.3-0.8. Within this range, each SoC can be used to estimate the SoH.

Reasonable Selection of dV
To guarantee the estimated efficiency and accuracy of this method, an appropriate dV needs to be selected. A reasonable choice of dV is 0.02 V.
The larger the dV, the better the fitting effect and the more accurate the SoH estimation results. As can be seen from Figure 6, due to the relative amplification of noise, there are many obvious interference points when dV = 0.005 V. With the increase of dV, the noise is reduced relatively and dQ/dV also gets closer to the total capacity Q2 of the current battery. In an extreme case, that is, the value of dV reaches the charging cut-off voltage, dQ/dV is numerically equal to Q2. According to Equation (2), the estimated SoH at this time is the actual SoH. Therefore, as mentioned above, the SoH estimation results become more accurate with increasing dV. However, with the increase of dV, the time required to calculate dQ also increases correspondingly, that is, the estimated efficiency is reduced. Additionally, nonlinear characteristics of the battery state of charge-open circuit voltage (SoC-OCV) curve determines the time taken to calculate dQ during the entire charging process is different. Overall, it is the maximum time to calculate dQ that limits the estimated efficiency. In order to improve the estimation efficiency on the basis of ensuring the accuracy of this method, the root mean square error (RMSE) of the fitting curves with SoC = 0.4 and the maximum time to calculate dQ is obtained by experiment, respectively, when dV = 0.005 V, 0.01 V, 0.02 V, 0.03 V and 0.05 V, as shown in Figure 7.
It can be seen in Figure 7 that the RMSE comes close to 0 gradually with increasing dV. When dV increases to 0.2 V, the declining trend of RSME is obviously slowed down, and there is already no significant difference in the RMSE at dV = 0.03 V and 0.05 V. As for the maximum time, it can be seen from Figure 7 that the maximum time growth trend is linear with increasing dV. Generally, as can be seen in Figure 6, when dV is 0.01 V or less, the noise is relatively amplified and the linearity of dQ and SoH is not satisfactory enough. However, when dV is 0.03 V or more, the maximum required time to calculate dQ is too long. Consequently, dV = 0.02 V is a reasonable option. Based on the above analysis, 11 fitting curves are selected. The third-order fitting R 2 of the fitting curves are shown in Table 2. It can be found that the R 2 of the selected fitting curves are all above 0.95, which indicates that the fitting effect is satisfactory.

Reasonable Selection of dV
To guarantee the estimated efficiency and accuracy of this method, an appropriate dV needs to be selected. A reasonable choice of dV is 0.02 V.
The larger the dV, the better the fitting effect and the more accurate the SoH estimation results. As can be seen from Figure 6, due to the relative amplification of noise, there are many obvious interference points when dV = 0.005 V. With the increase of dV, the noise is reduced relatively and dQ/dV also gets closer to the total capacity Q 2 of the current battery. In an extreme case, that is, the value of dV reaches the charging cut-off voltage, dQ/dV is numerically equal to Q 2 . According to Equation (2), the estimated SoH at this time is the actual SoH. Therefore, as mentioned above, the SoH estimation results become more accurate with increasing dV. However, with the increase of dV, the time required to calculate dQ also increases correspondingly, that is, the estimated efficiency is reduced. Additionally, nonlinear characteristics of the battery state of charge-open circuit voltage (SoC-OCV) curve determines the time taken to calculate dQ during the entire charging process is different. Overall, it is the maximum time to calculate dQ that limits the estimated efficiency. In order to improve the estimation efficiency on the basis of ensuring the accuracy of this method, the root mean square error (RMSE) of the fitting curves with SoC = 0.4 and the maximum time to calculate dQ is obtained by experiment, respectively, when dV = 0.005 V, 0.01 V, 0.02 V, 0.03 V and 0.05 V, as shown in Figure 7.

SoH Estimation Results
In order to observe the estimation results of the method comprehensively, the SoH estimation effect on the SoC range of 0.3-0.8 is analyzed. The test SoC point is taken every 0.005 interval for every charging process. A total of 40,000 points are selected. Figure 8A,B show the comparison results and estimated error, respectively. It can be seen that the estimated SoH has a satisfactory agreement It can be seen in Figure 7 that the RMSE comes close to 0 gradually with increasing dV. When dV increases to 0.2 V, the declining trend of RSME is obviously slowed down, and there is already no significant difference in the RMSE at dV = 0.03 V and 0.05 V. As for the maximum time, it can be seen from Figure 7 that the maximum time growth trend is linear with increasing dV. Generally, as can be  Figure 6, when dV is 0.01 V or less, the noise is relatively amplified and the linearity of dQ and SoH is not satisfactory enough. However, when dV is 0.03 V or more, the maximum required time to calculate dQ is too long. Consequently, dV = 0.02 V is a reasonable option.

SoH Estimation Results
In order to observe the estimation results of the method comprehensively, the SoH estimation effect on the SoC range of 0.3-0.8 is analyzed. The test SoC point is taken every 0.005 interval for every charging process. A total of 40,000 points are selected. Figure 8A,B show the comparison results and estimated error, respectively. It can be seen that the estimated SoH has a satisfactory agreement with the actual SoH. Figure 7. The relationship of RSME of fitting curve, the maximum time to calculate dQ and dV.

SoH Estimation Results
In order to observe the estimation results of the method comprehensively, the SoH estimation effect on the SoC range of 0.3-0.8 is analyzed. The test SoC point is taken every 0.005 interval for every charging process. A total of 40,000 points are selected. Figure 8A,B show the comparison results and estimated error, respectively. It can be seen that the estimated SoH has a satisfactory agreement with the actual SoH. As mentioned above, a total of 11 fitting curves belonging to different SoC are selected. The more the fitting curves, the more accurate the estimation results. In Figure 8B, because there are more fitting curves in regions a and b, the overall errors of these two regions are smaller than that of other regions. This proves that the error of SoH estimation decreases as the number of fitting curves increases. On As mentioned above, a total of 11 fitting curves belonging to different SoC are selected. The more the fitting curves, the more accurate the estimation results. In Figure 8B, because there are more fitting curves in regions a and b, the overall errors of these two regions are smaller than that of other regions. This proves that the error of SoH estimation decreases as the number of fitting curves increases. On the whole, the SoH error estimated by this method is within 5% and the MAE is 1.08%, which has satisfactory accuracy.

Influence of SoC Estimation Error on SoH Estimation Results
In this method, the sampling points are selected according to SoC. The magnitude of the SoC estimation error directly affects the estimation effect of SoH. Therefore, a comparative experiment is designed under different SOC estimation errors. As can be seen in Figure 9, when the SoC estimation error is within 1%, the SoH estimation results are not much different from the SoC error of 0 and the MAE is only 1.17%. Therefore, when an SoC estimation method can limit the error to 1%, the influence of the SoC estimation error on the SoH estimation is negligible and the SoH estimation effect is satisfactory. In addition, in Figure 9, when the error of the SoC is 2%, the estimation results of the SoH have shown a few obvious fluctuation points and MAE has risen to 1.36%. Combined with Figure 10B, there are a few obvious errors at both ends of the figure. As for the cause of the above phenomenon, it is the degree of difference in the fitting curves of different SoC regions. As can be seen from the analysis in Section 4.1 and Figure 10, in the SoC interval of 0.45-0.75, the difference in the different fitting curves is not large, and the relationship between dQ/dV and SoH under different SoCs is quite similar, that is, the accuracy of the SoC estimation in this interval has relatively inconspicuous influence on SoH estimation, while in the SoC region at both ends there are a few obvious fluctuation points, that is, when the error of the SoC is amplified to 2%, the estimation accuracy of the SoH estimation method has begun to be affected by the SoC error in the region of the SoC from 0.3-0.45. Fortunately, all the commonly used SoC estimation methods limit the SoC error to 1% or less [21][22][23]. Therefore, in this paper, the selection of the SoC estimation methods has not been discussed. In the following work, the impact of different SoC methods on the SoH estimation method will be discussed in detail. estimation error directly affects the estimation effect of SoH. Therefore, a comparative experiment is designed under different SOC estimation errors. As can be seen in Figure 9, when the SoC estimation error is within 1%, the SoH estimation results are not much different from the SoC error of 0 and the MAE is only 1.17%. Therefore, when an SoC estimation method can limit the error to 1%, the influence of the SoC estimation error on the SoH estimation is negligible and the SoH estimation effect is satisfactory. In addition, in Figure 9, when the error of the SoC is 2%, the estimation results of the SoH have shown a few obvious fluctuation points and MAE has risen to 1.36%. Combined with Figure 10B, there are a few obvious errors at both ends of the figure. As for the cause of the above phenomenon, it is the degree of difference in the fitting curves of different SoC regions. As can be seen from the analysis in Section 4.1 and Figure 10, in the SoC interval of 0.45-0.75, the difference in the different fitting curves is not large, and the relationship between dQ/dV and SoH under different SoCs is quite similar, that is, the accuracy of the SoC estimation in this interval has relatively inconspicuous influence on SoH estimation, while in the SoC region at both ends there are a few obvious fluctuation points, that is, when the error of the SoC is amplified to 2%, the estimation accuracy of the SoH estimation method has begun to be affected by the SoC error in the region of the SoC from 0.3-0.45. Fortunately, all the commonly used SoC estimation methods limit the SoC error to 1% or less [21][22][23]. Therefore, in this paper, the selection of the SoC estimation methods has not been discussed. In the following work, the impact of different SoC methods on the SoH estimation method will be discussed in detail.  estimation error directly affects the estimation effect of SoH. Therefore, a comparative experiment is designed under different SOC estimation errors. As can be seen in Figure 9, when the SoC estimation error is within 1%, the SoH estimation results are not much different from the SoC error of 0 and the MAE is only 1.17%. Therefore, when an SoC estimation method can limit the error to 1%, the influence of the SoC estimation error on the SoH estimation is negligible and the SoH estimation effect is satisfactory. In addition, in Figure 9, when the error of the SoC is 2%, the estimation results of the SoH have shown a few obvious fluctuation points and MAE has risen to 1.36%. Combined with Figure 10B, there are a few obvious errors at both ends of the figure. As for the cause of the above phenomenon, it is the degree of difference in the fitting curves of different SoC regions. As can be seen from the analysis in Section 4.1 and Figure 10, in the SoC interval of 0.45-0.75, the difference in the different fitting curves is not large, and the relationship between dQ/dV and SoH under different SoCs is quite similar, that is, the accuracy of the SoC estimation in this interval has relatively inconspicuous influence on SoH estimation, while in the SoC region at both ends there are a few obvious fluctuation points, that is, when the error of the SoC is amplified to 2%, the estimation accuracy of the SoH estimation method has begun to be affected by the SoC error in the region of the SoC from 0.3-0.45. Fortunately, all the commonly used SoC estimation methods limit the SoC error to 1% or less [21][22][23]. Therefore, in this paper, the selection of the SoC estimation methods has not been discussed. In the following work, the impact of different SoC methods on the SoH estimation method will be discussed in detail.

Future Discussions
This paper proposes a SoH estimation method based on incremental capacity analysis, which has satisfactory accuracy and estimated efficiency. However, there are still many aspects about the method that can be optimized and improved, including the following two aspects: (1) Filtering and denoising. The selection of dV is directly related to the estimated efficiency of the SoH estimation. A smaller dV means a shorter SoH estimation time. However, with the decrease of dV, the noise is relatively amplified. If the influence of noise can be reduced or eliminated, then a smaller dV can be selected, and the estimated efficiency of this method can be further

Future Discussions
This paper proposes a SoH estimation method based on incremental capacity analysis, which has satisfactory accuracy and estimated efficiency. However, there are still many aspects about the method that can be optimized and improved, including the following two aspects: (1) Filtering and denoising. The selection of dV is directly related to the estimated efficiency of the SoH estimation. A smaller dV means a shorter SoH estimation time. However, with the decrease of dV, the noise is relatively amplified. If the influence of noise can be reduced or eliminated, then a smaller dV can be selected, and the estimated efficiency of this method can be further improved. A Kalman filter, or other method, will be considered in future works to eliminate noise to the maximum extent. (2) Further research on the SoH estimation of HEVs. In actual operation, due to fuel economy considerations, the energy management strategy of HEVs tends to maintain SoC at 0.5 [24][25][26]. This method has the best estimation effect when the SoC = 0.5 or so, which means only 2-3 fitting curves are needed for HEVs to obtain a very satisfactory SoH estimation effect. Therefore, the next research direction is to apply the SoH estimation method to HEVs.
In addition, due to the long time of battery aging experiments, this paper only verifies the feasibility of the method at 25 • C. In the next stage of work, we will carry out the aging test of battery cells and battery packs at different temperatures. The universality and practicability of the proposed method will be further studied.

Conclusions
A fast and efficient online SoH estimation method based on ICA is proposed in this paper. The relationship between dQ/dV and SoH at different SoCs is analyzed and a black box model for SoH estimation is established based on this relationship. In this model, only an arbitrary dQ/dV is required to estimate the corresponding SoH. The influence of different SoC ranges on SoH estimation is analyzed, and the estimated accuracy and efficiency are balanced by choosing the reasonable dV. The maximum time to estimate SoH is less than six minutes and the estimated error is within 5%. Additionally, the mean absolute error (MAE) is 1.08%. Compared with traditional ICA methods, this method removes the dependence on IC peaks, so that most of the dQ/dV in the whole charging process can be used to estimate the SoH, which shows a strong adaptability.