Non-Invasive Method-Based Estimation of Battery State-of-Health with Dynamical Response Characteristics of Load Surges

: Battery state-of-health (SOH) estimation is an effective approach to evaluate battery reliability and reduce maintenance costs for battery-based backup power supply systems. This paper proposes a novel SOH estimation method for batteries, which only uses the response characteristics of load surges and is, therefore, non-destructive to the estimated battery and its system. The discrete wavelet transform (DWT) method based on multi-resolution analysis (MRA) is used for wavelet energy features extraction, and the fuzzy cerebellar model neural network (FCMNN) is introduced to design the battery SOH estimator. The response voltage signals to load surges are used in the training and detection process of the FCMNN. Compared to conventional methods, the proposed method only exploits characteristics of online response signals to the inrush currents rather than injecting interference signals into the battery. The effectiveness of the proposed method is validated by detailed


Introduction
As a backup power supply, battery energy storage has been widely used in energy storage systems because of its high power density and variable size, which can improve system reliability, reduce maintenance costs, and allow for longer standby times [1].When the power grid is in emergency conditions, these batteries are usually used to power critical loads.Lead-acid batteries, pioneers in the field of energy storage, are the most common type of rechargeable batteries globally [2] and represent the only closed-loop operational battery technology.Despite their relatively lower energy density, they exhibit robust performance, simple control, and lower costs, making them more environmentally friendly compared to other battery types.In consequence, they are primarily used in backup power scenarios for their high safety levels.They are extensively employed in various industries, including data centers, finance, communications, and transportation.Although their growth rate may not match that of lithium-ion batteries, they maintain a leading position in terms of capacity (GWh).This is mainly attributed to the emergence of mild and start-stop hybrid power batteries as the primary growth sector for advanced lead-acid batteries [3].However, the aging and insufficient operation of lead-acid batteries can lead to power outages, shutdowns, and even severe accidents.The health status of the batteries is thus closely related to the security and stability of the system [4].
Battery SOH indicates the health status of the battery, which is usually defined as capacity degradation or internal resistance growth [5][6][7][8].A lead-acid battery is composed of positive and negative plates, separators, and electrolytes.The chemical reactions between the active material of the positive and negative plates and the sulfuric acid solution are the charging and discharging process of the lead-acid battery.The active material of the The decrease in battery SOH is usually caused by the gradual accumulation of lead sulfate in the active material.This part of lead sulfate will not be converted back to active material during the charging process.It inhibits the contact between the electrolyte and the active material.Over time, the small lead sulfate crystals on the positive and negative electrodes will accumulate into large lead sulfate crystals.Therefore, as the battery SOH decreases, the internal resistance of the battery is gradually increasing.
When the available capacity reduces to a preset threshold known as end-of-life value, the invalid battery needs to be replaced.The capacity estimation of batteries in engineering applications mainly requires periodic discharge tests [9].The results obtained from this method are accurate but time-consuming (at least five hours to fully discharge).The measurement time can be shortened by increasing the discharge current, but high currents will lead to incomplete discharge, increasing the error of battery capacity estimation and even causing greater damage to the battery.
Relevant studies have been conducted on estimating battery capacity based on the correlation between battery voltage and capacity.Among them, "coup de fouet" is a particular phenomenon for lead-acid batteries that can be utilized for capacity estimation [10].When a battery is fully charged for a long time, the voltage will drop sharply and rise again within the first few minutes of the battery discharge.The voltage drop observed during the initial moments of discharging a lead-acid battery, following a prolonged charging period is attributed to polarization effects occurring at the electrode surfaces.When the battery ceases charging and transitions to discharge, polarization phenomena alter the chemical state of the electrode surfaces.Upon commencing discharge, this polarization effect briefly causes a voltage reduction.This transient voltage decrease is a result of the battery requiring a certain amount of time to readjust the electrode surface conditions to accommodate the discharge process.The degree of voltage drop is strongly related to the actual capacity of the battery [11].Therefore, the relationship between them can be used to estimate the capacity of the battery.However, this method depends on the battery's state of charge (SOC), which needs to be offline for some time.Another method, the so-called dual-pulses discharge load test, was introduced to estimate the available capacity of the battery in [12], which uses two identical discharge pulse currents.Although this method is user-friendly, it is not suitable for online dynamic estimation and requires two additional pulses to the battery.Therefore, it is feasible to distinguish the SOH of batteries by open-circuit voltage slopes.
A fusion method-based extended Kalman filter (EKF) using a neural network was presented to estimate the SOC, and then, the recursive least squares were used to find the slope of the open circuit voltage versus the SOC curve as the input for the fuzzy system [13].This method has good performance under continuous charging or discharging situations.However, the slope difference is not obvious when the battery pack contains a large number of batteries.The impedance and resistance of batteries increase with the loss of battery life and capacity.Combined with impedance measurement, fuzzy logic data analysis can be used to estimate the SOH of a battery effectively [14,15].While the SOH of batteries can be reflected by the variation trend of internal impedance, it needs more accurate and expensive equipment due to the low impedance of batteries, i.e., less than one-tenth of a milliohm.The measurement of internal impedance will also interfere with excessive fluctuations of discharge current and radiation noise [16].
In addition, estimation models have been developed using data-fitting methods.For instance, the neural network is used to analyze the relationship between the SOH and the battery voltage, current, temperature and other parameters directly measured from batteries and then estimate the capacity of the battery [17][18][19].Sui et al. extracted features from battery voltage using the fuzzy entropy method [20].This method performs well with the change of aging temperature since it can efficiently capture the impact due to aging temperature.With the fractional-order model, a data-driven approach was used to obtain model parameters, and finally, SOC and SOH are estimated by multi-scale double-extended Kalman filtering [21].However, this method is too model-dependent.
In [22], aiming to estimate the degradation of battery capacity, the mathematical model of the battery was established to judge battery capacity.The ideal mathematical model is more suitable for the ideal operating mode.In [23], based on two working states of the battery, the relationship between voltage and battery aging in the discharge state and the relationship between internal resistance and battery aging in the floating charge state were found.Two results were obtained from the corresponding two empirical models weighted as real SOH of batteries.The extended Kalman filter (EKF) was employed to identify the model parameters and estimate the SOH, which reduces the computation of diagnosis and prognosis [24].However, it is just suited in the case that the battery discharge rate is constant.
Another method estimates the remaining capacity of the battery by injecting sinusoidal voltage or current signals into the battery and analyzing its response signals [25,26].However, this method needs an accurate mathematical model; the estimation error will be larger if the model is not accurate.Partial charging curves are commonly used in SOH estimation.In [27], the Kendall coefficients were employed for feature selection.Combined with the linear regression model, the recursive method was used to fit the features with less computation burden.He et al. proposed a new incremental capacity feature from partial constant current charging data [28].The estimated SOH error is within 5% for the most cycle estimation.
Moreover, the sensitive equipment on the DC bus will be affected by the injected harmonic signal, resulting in inaccurate measurement or incorrect operation.A non-linear observer was established to identify bulk capacitance, which converts the SOH estimation of the battery into parameter estimation based on the model describing the dynamic characteristics of the battery with good robustness [29].However, the selection of a battery model has a great influence on the effectiveness and adaptability of the estimation algorithm.
In this paper, a novel method is proposed for battery SOH estimation by using the response signals of the battery to inrush currents, where an inrush current is a large but short-time current triggered by the load.From the response voltage signal during the inrush current, the proposed estimation method first utilizes wavelet transform based on MRA to extract energy features.In addition, a battery SOH estimator is built based on FCMNN.For example, the proposed method can be used in the backup energy system where batteries are infrequently charged or discharged.The inrush current would be offered to activate the motor-based high-voltage breaker.Instead of using the costly charging and capacity curves of the batteries, the FCMNN-based estimator directly uses this fast and simple inrush current to estimate the battery SOH, in turn, does not require any external signal injection, and there is no influence on the system, including the battery.
The rest of this paper is organized as follows: Section 2 introduces the theoretical background and the model-based method for battery SOH estimation.Section 3 presents simulation and experimental results.Finally, concluding remarks are provided in Section 4.

Equivalent Circuit Model of the Battery
Various battery models have been proposed [30][31][32][33], most of which are the equivalent circuit models.The Ceraolo model is widely used to demonstrate electrochemistry, as shown in Figure 1.Here, E m (θ, SOC) is the open circuit voltage of the battery when switch S 1 and switch S 2 are all off, which is related to the SOC of the battery and the electrolyte temperature θ; Ip is the parasitic branch current, I p ≈ 0; U o is the battery terminal voltage; I o is the battery discharge current; R 0 is the internal resistance of the battery; R x − C x (x = 1, 2, . .., n) are the parallel blocks modeling the transient behavior of the battery, which represent the activity of the battery associated with the SOH.

Equivalent Circuit Model of the Battery
Various battery models have been proposed [30][31][32][33], most of which are the equivalent circuit models.The Ceraolo model is widely used to demonstrate electrochemistry, as shown in Figure 1.Here, Em(θ, SOC) is the open circuit voltage of the battery when switch S1 and switch S2 are all off, which is related to the SOC of the battery and the electrolyte temperature θ; Ip is the parasitic branch current, Ip ≈ 0; Uo is the battery terminal voltage; Io is the battery discharge current; R0 is the internal resistance of the battery; Rx − Cx (x = 1, 2, …, n) are the parallel blocks modeling the transient behavior of the battery, which represent the activity of the battery associated with the SOH.

AC Charger
Lead-acid battery equivalent circuit network [24] Load and Charger The internal impedance is a common performance indicator for the SOH estimation.According to Kirchhoff's voltage law, the battery terminal voltage will drop suddenly when an inrush current happens.The parameters of the battery, such as the internal resistance, have obvious changes, accompanied by the increase in the number of charging and discharging.Hence, the battery terminal voltage drops at different aging levels are different.

Dynamical Response Characteristics Analysis
The voltage response to pulsed load for lead acid battery at different SOHs and SOCs are shown in Figures 2 and 3 by using the SIMULINK toolbox in MATLAB.As can be seen from Figure 2, affected by inrush current, minimum and maximum voltage drops occur under 100% and 60% SOHs, respectively.It confirms that the battery aging will bring more evidently changes subjected to a pulse load than a new battery.The internal impedance is a common performance indicator for the SOH estimation.According to Kirchhoff's voltage law, the battery terminal voltage will drop suddenly when an inrush current happens.The parameters of the battery, such as the internal resistance, have obvious changes, accompanied by the increase in the number of charging and discharging.Hence, the battery terminal voltage drops at different aging levels are different.

Dynamical Response Characteristics Analysis
The voltage response to pulsed load for lead acid battery at different SOHs and SOCs are shown in Figures 2 and 3 by using the SIMULINK toolbox in MATLAB.As can be seen from Figure 2, affected by inrush current, minimum and maximum voltage drops occur under 100% and 60% SOHs, respectively.It confirms that the battery aging will bring more evidently changes subjected to a pulse load than a new battery.
(1) Effect of battery SOH According to the equivalent circuit in Figure 1 and the simulation results in Figure 2, the relation between the response signal of an inrush current and internal parameters can be deduced as follows: (a) t < t 0 : switch S 1 is on, and switch S 1 is off, then the battery is in the charging situation.
When the actual direction of the charging current is opposite to the reference direction, that is, I o < 0, then U o > E m (t 0− ), u cx < 0 (x = 1, 2, . .., n), then the terminal voltage of the battery is greater than its open circuit voltage.
(b) t = t 0 : switch S 1 is on, then the load is connected to the battery.If the load current is larger than the charger can offer, the battery will switch to the discharging state.On the other hand, the voltages of parallel R-C blocks would not change at the moment because the voltage of the capacitor u cx cannot change instantaneously.Therefore, we have the following equation: Energies 2024, 17, 583 5 of 17 Since E m (t 0+ ) = E m (t 0− ), the terminal voltage changed at this moment is Equation ( 3) indicates that ∆U o (t 0 ) is related to the internal resistant R 0 of the battery.As shown in Figure 2, for the same inrush current, the greater R 0 is, the deeper the voltage drop is.
(c) t 0 < t < t 1 : switch S 1 and switch S 2 are both on, U o will decrease in accordance with Equation (4): where and C m is the capacity of the battery.Therefore, the response voltage signal from t 0 to t 1 contains the main SOH information of the battery.(d) t = t 1 : switch S 1 is off, that is when the inrush current disappears.
For E m (t (e) t > t 1 : switch S 1 and switch S 2 are all off, the battery is in the recovery process: According to Equation ( 11), the recovery process only depends on the dynamical parameter Energies 2024, 17, x FOR PEER REVIEW 4 of 18

Dynamical Response Characteristics of Inrush Current
2.1.1.Equivalent Circuit Model of the Battery Various battery models have been proposed [30][31][32][33], most of which are the equivalent circuit models.The Ceraolo model is widely used to demonstrate electrochemistry, as shown in Figure 1.Here, Em(θ, SOC) is the open circuit voltage of the battery when switch S1 and switch S2 are all off, which is related to the SOC of the battery and the electrolyte temperature θ; Ip is the parasitic branch current, Ip ≈ 0; Uo is the battery terminal voltage; Io is the battery discharge current; R0 is the internal resistance of the battery; Rx − Cx (x = 1, 2, …, n) are the parallel blocks modeling the transient behavior of the battery, which represent the activity of the battery associated with the SOH.

AC Charger
Lead-acid battery equivalent circuit network [24] Load and Charger The internal impedance is a common performance indicator for the SOH estimation.According to Kirchhoff's voltage law, the battery terminal voltage will drop suddenly when an inrush current happens.The parameters of the battery, such as the internal resistance, have obvious changes, accompanied by the increase in the number of charging and discharging.Hence, the battery terminal voltage drops at different aging levels are different.

Dynamical Response Characteristics Analysis
The voltage response to pulsed load for lead acid battery at different SOHs and SOCs are shown in Figures 2 and 3 by using the SIMULINK toolbox in MATLAB.As can be seen from Figure 2, affected by inrush current, minimum and maximum voltage drops occur under 100% and 60% SOHs, respectively.It confirms that the battery aging will bring more evidently changes subjected to a pulse load than a new battery.(1) Effect of battery SOH According to the equivalent circuit in Figure 1 and the simulation results in Figure 2, the relation between the response signal of an inrush current and internal parameters can be deduced as follows: (a) t < t0: switch S1 is on, and switch S1 is off, then the battery is in the charging situation.
When the actual direction of the charging current is opposite to the reference direction, that is, Io < 0, then Uo > Em(t0−), ucx < 0 (x = 1, 2, …, n), then the terminal voltage of the battery is greater than its open circuit voltage.
(b) t = t0: switch S1 is on, then the load is connected to the battery.If the load current is larger than the charger can offer, the battery will switch to the discharging state.On the other hand, the voltages of parallel R-C blocks would not change at the moment because the voltage of the capacitor ucx cannot change instantaneously.Therefore, we have the following equation: Since Em(t0+) = Em(t0−), the terminal voltage changed at this moment is Equation ( 3) indicates that ΔUo(t0) is related to the internal resistant R0 of the battery.As shown in Figure 2, for the same inrush current, the greater R0 is, the deeper the voltage drop is.
(c) t0 < t < t1: switch S1 and switch S2 are both on, Uo will decrease in accordance with Equation (4): where 0 ( ) (2) Effect of Battery SOC Similarly, the variation of the battery SOC also influences the steepness of the terminal voltage drops.In this case, we supposed that the AC source is offline.Then, the voltage response to the same inrush current at different SOCs and the same SOH is shown in Figure 3.It is obvious that the smaller the SOC of the battery is, the larger the battery terminal voltage drops.When the battery SOC is at a low value, there is a higher presence of non-active substances in the battery, reducing the availability of active materials.Simultaneously, the concentration of solutes in the electrolyte decreases, resulting in a lower rate of electrochemical reactions.This subsequently causes an increase in internal resistance within the battery, influencing the voltage response and causing a more significant voltage drop.A battery with a low SOC will express greater battery voltage change to sustain high currents.It illustrates that battery voltage drops are related not only to SOH but also to SOC.However, there are some different characteristics between these two kinds of response signals, which can be seen from the following formulas deduced in a similar way above.
(a) t < t 0 : switch S 1 and switch S 2 are all off, for I o = 0, then The initial value of E m (t 0− ) is decided by the SOC under the same temperature.(b) t = t 0 : switch S 1 turns on at this time, then the drop voltage will be According to Equation ( 13), the voltage drops at the moment when the inrush current is injected, independent of the initial voltage E m (t 0− ).Therefore, it is not influenced by the SOC.
(c) t 0 < t < t 1 : switch S 1 is on, then (d) t ≥ t 1 : switch S 1 turns off, the step-up voltage is the same as Equation (10).And the recovery voltage is also the same as Equation (11).
(3) Discussions of Dynamical Response Characteristics According to (10), the step voltage at the moment when the inrush current disappeared only depends on the internal resistant R 0 .Therefore, for different SOCs but the same SOH of a battery, the values of step voltage are equal under the same step current.On the contrary, the value of step voltage will be bigger for an aging battery with a bigger internal resistant R 0 .
In summary, the voltage response to the pulsed load is related to the SOH and the SOC of the battery.However, the relationship between SOH and the voltage response signal is not explicit, and the voltage signal is susceptible to interference.Therefore, in this paper, a feature extraction method based on wavelet transform is adopted.Then, an estimator with a moderate computational complexity is designed using a fuzzy cerebellar model neural network.same SOH of a battery, the values of step voltage are equal under the same step current.On the contrary, the value of step voltage will be bigger for an aging battery with a bigger internal resistant R0.

Battery SOH Estimation Method
In summary, the voltage response to the pulsed load is related to the SOH and the SOC of the battery.However, the relationship between SOH and the voltage response signal is not explicit, and the voltage signal is susceptible to interference.Therefore, in this paper, a feature extraction method based on wavelet transform is adopted.Then, an estimator with a moderate computational complexity is designed using a fuzzy cerebellar model neural network.

Flowchart of Estimation Process
Aiming to estimate the real-time SOH of the battery using only the voltage response signal of the inrush current, as shown in Figure 4

Feature Extraction Based on DWT
There are various signal decomposition methods, such as DWT, ensemble empirical mode decomposition (EEMD) [34], and improved complete ensemble empirical mode decomposition with adaptive noise (ICEEMDAN) [35], etc.In this study, the voltage re-

Feature Extraction Based on DWT
There are various signal decomposition methods, such as DWT, ensemble empirical mode decomposition (EEMD) [34], and improved complete ensemble empirical mode decomposition with adaptive noise (ICEEMDAN) [35], etc.In this study, the voltage response of the battery contains rich information on SOH and SOC; it is affected by the current's magnitude and duration, showing a non-stationary nature.Wavelet analysis is used since it can obtain the information in both time and frequency domains from raw degradation signal and is well suited for processing non-stationary signals [32].
With the scaled and shifted wavelets, the multi-resolution signal analysis method based on DWT can avoid the numerical complexity and be easily implemented for online applications [33].Given an arbitrary signal f (n) of finite energy, its multi-resolution decomposition at level M is defined by the formula: where a m,k is the scaling coefficient, d m,k is the wavelet coefficient.φ(n) and ψ(n) are the scaling function and wavelet function, respectively.
The parameters a m,k and d m,k can be calculated as: Here, h(m) is the low-pass filter, and g(m) is the high-pass filter.DWT based on MRA can be used to process data samples.For a signal V * , it can be decomposed into approximate coefficient a j and detailed coefficient d j at different scales by using low-pass filter h(n) and high-pass filter g(n), and a j can be further decomposed into a j+1 and d j+1 , etc. Coefficient reconstruction is a reverse decomposition process.The energy content of every detailed coefficient can be calculated as: where q j is the length of wavelet coefficients, and c j is the length of scaling coefficient at scale j.The number of decomposed layers is in accordance with the following rule: where N j,d is the number of detail coefficients of the j-th layer after wavelet decomposition, L is the supporting length of the wavelet function, and j is determined as the final scale.
Here, three different SOHs are chosen, and the voltage curves under the same SOC illustrate the process of feature selection.The information of these three voltage curves is shown in Figure 5a.According to the decomposition principle shown in Equation ( 15), the number of decomposition layers is finally determined to be 14 layers.

Battery SOH Estimator Based on FCMNN
FCMNN can be thought of as a learning mechanism imitating the cerebellu human being.It is often referred to as an associative NN, where only a small sub dressed by the input vector determines instantaneous output.Compared with bac agation neural networks (BPNNs) and deep learning algorithms, FCMNN has vantages of faster learning speed, better generalization and less computation res which is increasingly being employed for classification and forecast [33,36,37].Th Figure 5b-d show the reconstructed curves of the 1st, 10th, and 14th layer detail coefficients.These layers contain information about the battery's internal resistance, open circuit voltage and rebound voltage, and battery polarization, respectively.The energy of each level of detail coefficient is obtained and used as an effective classification feature.It is found that the reconstruction curves of layers 1-8, 9-11, and 12-14 are similar to those in Figure 5b, Figure 5c, and Figure 5d, respectively.The detailed energy component of each layer is calculated using Equation (18).From the above analysis, only three layers need to be selected as features.

Battery SOH Estimator Based on FCMNN
FCMNN can be thought of as a learning mechanism imitating the cerebellum of a human being.It is often referred to as an associative NN, where only a small subset addressed by the input vector determines instantaneous output.Compared with backpropagation neural networks (BPNNs) and deep learning algorithms, FCMNN has the advantages of faster learning speed, better generalization and less computation resources, which is increasingly being employed for classification and forecast [33,36,37].Therefore, FCMNN is proposed as a battery SOH estimation model in this paper.
Figure 6 shows the architecture of the battery SOH estimator based on FCMNN and its training method.The structure of the battery SOH estimator based on FCMNN has five spaces [33].

The cost function
( ) ( ) ( ) Layer 3. Receptive-field space Rs: It is formed by multiple input regions.
( ) Layer 3. Receptive-field space R s : It is formed by multiple input regions.
where b jk is the general basis function associated with the jth layer and kth block of association memory space.Layer 4. Weight memory space W s : The weight memory space connects the receptivefield space with the output space.Where w jko denotes the connecting weight of o-th (o = 1, 2, . .., n o ) output associated with the k-th block of j-th receptive field.n o is the number of battery SOH estimator outputs.
Therefore, we have the last layer, i.e., the output space.Where the output of the proposed estimator can be represented by:

Experiment Platform
Figure 7 shows an experimental platform constructed for subsequent experiments: a battery, a battery charge-discharge testing system, and a computer.The battery has a capacity of 20 AH and a rated voltage of 12 V.The pulsed current and the duration used in the experiment are listed in Table 1.
Energies 2024, 17, x FOR PEER REVIEW 11 of 18 Therefore, we have the last layer, i.e., the output space.Where the output of the proposed estimator can be represented by: ( ) ( )

Experiment Platform
Figure 7 shows an experimental platform constructed for subsequent experiments: a battery, a battery charge-discharge testing system, and a computer.The battery has a capacity of 20 AH and a rated voltage of 12 V.The pulsed current and the duration used in the experiment are listed in Table 1.

Battery charge and discharge test system
Name: Battery charge and discharge test system 4CH Model: CBTS-50A-060V-4CH Company: China Electronics Technology Group Corporation Features: The device has 4 test channels, which can test multiple cells of the battery pack at the same time; the voltage range is 0-60V, the maximum current is 50A, the maximum power of the channel is 3kW, and the test accuracy is ±0.1% FSR; the current hardware response time is <100ms; It has test functions such as constant current charging, constant voltage charging, constant current to constant voltage charging, and constant current discharge.

Simulation Case Study
In this section, two examples were modeled in Simulink of MATLAB for simulation study.

Auxiliary Supply System of Transformer Substation:
A classical auxiliary DC supply system with the proposed online dynamical battery SOH estimator is shown in Figure 8.The rectifier converts the 220 V AC voltage to 220 V DC voltage.The battery bank paralleled with the DC bus is used as a backup supply during an outage.This system has two kinds of loads: the constant load and the pulse load.Due to high inrush currents that the rectifier cannot offer, most or all of the inrush currents are supplied by a battery bank to activate the motor-based high-voltage breaker when the AC grid is normal or during outages.Therefore, the response signal can be used to estimate the real-time battery bank SOH with the proposed method.The main parameters of the auxiliary supply system in this simulation case are de scribed in Table 2.When switch S1 is on, it means that the AC grid is working well; other wise, the grid will fail.Table 3 lists different operated values and durations of inrush cur rents.Usually, the battery is regarded to be at the end of life when the capacity reaches 80% of the rated capacity [38].The SOC of the battery varies from 28~100%, and the range of SOH is 70~100%.Two estimation models are built by FCMNN for load surges that may occur when the AC source is normal or faulty.Table 4 lists six intervals of the SOH of the battery.The main parameters of the auxiliary supply system in this simulation case are described in Table 2.When switch S 1 is on, it means that the AC grid is working well; otherwise, the grid will fail.Table 3 lists different operated values and durations of inrush currents.Usually, the battery is regarded to be at the end of life when the capacity reaches 80% of the rated capacity [38].The SOC of the battery varies from 28~100%, and the range of SOH is 70~100%.Two estimation models are built by FCMNN for load surges that may occur when the AC source is normal or faulty.Table 4 lists six intervals of the SOH of the battery.By using the feature extraction method, the battery voltage response signal is decomposed into 11 layers when the AC power supply is online, and the battery voltage response signal is decomposed into 13 layers when the AC power supply is offline.This difference is caused by different lengths of the impulse current time selected in Table 3.The time segment of the AC power supply online inrush current selection is short, and the time segment of the AC power supply offline inrush current selection is long.
According to Section 3, it can be noted that the detailed energy of decomposition can be roughly divided into three categories.Therefore, suitable energy can be selected as the characteristic quantity among different classes.Figure 9a shows a three-dimensional view of the detailed energy of the 1st, 2nd, and 9th layers when the AC power supply is online.Figure 9b shows the three-dimensional view of the AC power offline, and the detailed energy of the 4th, 5th, and 6th layers are selected.There were 300 data sets obtained in the normal operation of the AC source, and suppose that SOC is 100%.To analyze the performance of FCMNN, the samples were divided into three parts: training sets, test sets and validation sets, respectively.Without the loss of generality, 30 sets of data were randomly selected as the training set from 40 sets of data under each battery SOH interval, and the remaining 10 sets of data were used as the test sets.Thus, 180 sets of data were obtained as the training set data and 60 data sets were used as the test set data.Another 60 sets of data were used to validate the effectiveness of the well-trained SOH estimation model.
A total of 1020 sets of data were collected in AC source offline under different SOCs.Similarly, 120 sets of data were randomly selected as the training set from the 170 sets of data under each capacity interval, and 30 sets of data were used as the test sets.Thus, 720 sets of data were obtained as the training set data, 180 sets of data were used as the test set data, and the remaining 120 sets were for validation.Table 5 shows the sample partition and the classification accuracy of the testing set using an estimator based on a different kind of neural network.There were 300 data sets obtained in the normal operation of the AC source, and suppose that SOC is 100%.To analyze the performance of FCMNN, the samples were divided into three parts: training sets, test sets and validation sets, respectively.Without the loss of generality, 30 sets of data were randomly selected as the training set from 40 sets of data under each battery SOH interval, and the remaining 10 sets of data were used as the test sets.Thus, 180 sets of data were obtained as the training set data and 60 data sets were used as the test set data.Another 60 sets of data were used to validate the effectiveness of the well-trained SOH estimation model.
A total of 1020 sets of data were collected in AC source offline under different SOCs.Similarly, 120 sets of data were randomly selected as the training set from the 170 sets of data under each capacity interval, and 30 sets of data were used as the test sets.Thus, 720 sets of data were obtained as the training set data, 180 sets of data were used as the test set data, and the remaining 120 sets were for validation.Table 5 shows the sample partition and the classification accuracy of the testing set using an estimator based on a different kind of neural network.
A Back-Propagation Neural Network (BPNN) and Support Vector Machine (SVM) are used to compare with FCMNN in the paper.As listed in Table 5, the proposed method performs particularly well when the AC source is online.The FCMNN-based SOH estimator has better capabilities than the other two methods when the AC source is offline.The battery capacity and voltage used in the experiment did not reach 100 AH or 220 V.This can be regarded as a group of batteries connected in a series.In this case, the current is identical for each battery, so the voltage response should be the same.Therefore, for security purposes, a battery with a capacity of 20 AH and a rated voltage of 12 V was used in the experiment.Therefore, in order to verify the validity of the experiment, the battery capacity and voltage adjustment in the simulation are consistent with the experiment, listed in Table 6.A total of 864 sets of data were obtained under AC power offline.For the sample data, the training set is 600 sets of data, the test set is 120 sets of data, and the verification set is 144 sets of data.
The set simulation scenes are shown in Table 7.The final validation set results and test set results are shown in Figure 10.It can be seen that using the FCMNN-based estimator, the battery SOH can still be distinguished with a high classification accuracy of 97.5%.battery capacity and voltage adjustment in the simulation are consistent with the experiment, listed in Table 6.A total of 864 sets of data were obtained under AC power offline.For the sample data, the training set is 600 sets of data, the test set is 120 sets of data, and the verification set is 144 sets of data.
The set simulation scenes are shown in Table 7.The final validation set results and test set results are shown in Figure 10.It can be seen that using the FCMNN-based estimator, the battery SOH can still be distinguished with a high classification accuracy of 97.5%.

Experiment Verification
Figure 11 shows the voltage and current curves of the experimental battery under different SOHs.There were 420 data samples, including 380 samples for training and 40 samples for testing.Table 8 shows the SOH range of the experiment setting and the corre-

Experiment Verification
Figure 11 shows the voltage and current curves of the experimental battery under different SOHs.There were 420 data samples, including 380 samples for training and 40 samples for testing.Table 8 shows the SOH range of the experiment setting and the corresponding label.Using the proposed feature extraction method, the end condition was reached when the number of decomposition layers was selected to be 11 layers.Figure 12 shows the detailed energy of each layer obtained after decomposing 11 layers of a certain curve.It can be seen from Figure 12 that the detail energies of layers 6-9 have a strong linear correlation, so one of them was chosen.Here, the energy of the 8th layer was selected.Then, two layers with larger energy were chosen as features, namely the 10th and the 11th layers.The three features of each curve are shown in the three-dimensional graph, as shown in Figure 13.These four categories can be easily distinguished from the figure.If the current difference is large or only one current is used, it is easier to distinguish.The three features obtained from the experimental data were fed into the FCMNN,   It can be seen from Figure 12 that the detail energies of layers 6-9 have a strong linear correlation, so one of them was chosen.Here, the energy of the 8th layer was selected.Then, two layers with larger energy were chosen as features, namely the 10th and the 11th layers.The three features of each curve are shown in the three-dimensional graph, as shown in Figure 13.These four categories can be easily distinguished from the figure.If the current difference is large or only one current is used, it is easier to distinguish.The three features obtained from the experimental data were fed into the FCMNN, SVM, and BPNN models, and the results obtained are shown in Table 9.Table 9 shows It can be seen from Figure 12 that the detail energies of layers 6-9 have a strong linear correlation, so one of them was chosen.Here, the energy of the 8th layer was selected.Then, two layers with larger energy were chosen as features, namely the 10th and the 11th layers.The three features of each curve are shown in the three-dimensional graph, as shown in Figure 13.These four categories can be easily distinguished from the figure.If the current difference is large or only one current is used, it is easier to distinguish.
The three features obtained from the experimental data were fed into the FCMNN, SVM, and BPNN models, and the results obtained are shown in Table 9.Table 9 shows that the proposed battery estimation method based on FCMNN can obtain higher classification accuracy than SVM and BPNN algorithms.This is because the associative memory and receptive-field space within FCMNN enhance the network's capability to process uncertain information, adaptively learn, and generalize effectively.These improvements enable the estimator to better approximate complex non-linear mapping relationships.Compared with the traditional discharging method [9], both methods can achieve high accuracy; however, the traditional discharging method requires more than 5 h while the proposed method has the potential of online measurements.The three features obtained from the experimental data were fed into the FCMNN SVM, and BPNN models, and the results obtained are shown in Table 9.Table 9 shows that the proposed battery estimation method based on FCMNN can obtain higher classification accuracy than SVM and BPNN algorithms.This is because the associative memory and receptive-field space within FCMNN enhance the networkʹs capability to process uncertain information, adaptively learn, and generalize effectively.These improvements enable the estimator to better approximate complex non-linear mapping relationships.Com-  Moreover, the traditional method needs the battery to be fully discharged; otherwise, it would cause large measurement errors.For the proposed method based on FCMNN, no additional injections and interference are needed during the measurements, which is non-destructive to the battery and the entire system.Furthermore, due to the similarity in the operating principles of different types of batteries, the proposed voltage response characteristics are also applicable to these batteries, e.g., lithium-ion batteries.They are gradually becoming the preferred energy storage solution for various applications due to their high energy density and long lifespan characteristics; it is important to estimate their health state.Therefore, there is potential for applying the proposed FCMNN-based SOH estimator due to the load surges.

Conclusions
This paper presents a novel battery SOH estimation method based on the response characteristics of the battery terminal voltage to pulsed loads.DWT-based MRA is employed to extract effective features from the dynamical response voltage data to inrush currents.The SOH estimator based on FCMNN is then utilized to estimate the SOH of the battery dynamically online.The simulation and experiment results show that the proposed approach has better performance.The advantage of the proposed method is that it can avoid interference from rectifier equipment and the other sensitive loads connected to the DC bus.Meanwhile, no additional injected high currents are needed.The proposed method is suitable for the scenario where the battery bank is used as a backup power supply and is not charging and discharging frequently, such as the online UPS system and the auxiliary.The DC supply system was deployed in transformer substations with high voltage circuit breakers.It can also be used in other applications with load surges, such as the starting process of electric vehicles.

Figure 1 .
Figure 1.Battery equivalent circuit with load and charger.

Figure 2 .
Figure 2. Battery voltage drop affected by inrush current under different SOHs.

Figure 1 .
Figure 1.Battery equivalent circuit with load and charger.

Figure 1 .
Figure 1.Battery equivalent circuit with load and charger.

Figure 2 .
Figure 2. Battery voltage drop affected by inrush current under different SOHs.

Figure 3 .
Figure 3. Battery voltage drop affected by inrush current under different SOCs.
Based on FCMNN 2.2.1.Flowchart of Estimation Process Aiming to estimate the real-time SOH of the battery using only the voltage response signal of the inrush current, as shown in Figure 4, DWT-based MRA is introduced for feature extraction and FCMNN for the battery SOH estimation is designed by applying the fuzzy rule to machine learning algorithms.Raw data are divided into a training set and a testing set.The data sample base contains battery terminal voltage signals collected from different values of inrush currents under different SOCs.After the training, the proposed method can be used for the online battery SOH estimation with real-time sampled terminal voltage signals of the battery.

Figure 4 .
Figure 4.The flowchart of the proposed battery SOH estimation.

Figure 4 .
Figure 4.The flowchart of the proposed battery SOH estimation.

Figure 5 .
Figure 5. (a) Voltage curve of different SOHs.(b) The 1st layer detail coefficient reconstruction curve.(c) The 10th level detail coefficient reconstruction curve.(d) The 14th layer detail coefficient reconstruction curve.

Figure 6 .Layer 1 .
Figure 6.The architecture of the battery SOH estimator is based on FCMNN and its learning process.Layer 1. Input space Es: The input data is the energy feature data Ei, i = 1, 2, •••, ne, which are the features extracted from V * .i represents the i-th input feature, and ne is the number of input signals to the estimator.Layer 2. Association memory space As: It consists of nl layers and nb blocks.rijk denotes the output of the k-th layer and j-th block of the receptive-field basis function for the i-th input.Each block has two adjustable parameters, mean mijk, and variance vijk.The Gaussian function is employed as the receptive-field basis function in this paper.

Figure 6 .
Figure 6.The architecture of the battery SOH estimator is based on FCMNN and its learning process.Layer 1. Input space E s : The input data is the energy feature data E i , i = 1, 2, • • • , n e , which are the features extracted from V * .i represents the i-th input feature, and n e is the number of input signals to the estimator.Layer 2. Association memory space A s : It consists of n l layers and n b blocks.r ijk denotes the output of the k-th layer and j-th block of the receptive-field basis function for the i-th input.Each block has two adjustable parameters, mean m ijk , and variance v ijk .The Gaussian function is employed as the receptive-field basis function in this paper.

Figure 7 .Table 1 .
Figure 7. Experiment platform with equivalent pulsed load.(a) Battery charge and discharge test system.(b) Experimental battery and computer.

Figure 7 .
Figure 7. Experiment platform with equivalent pulsed load.(a) Battery and discharge test system.(b) Experimental battery and computer.

18 Figure 8 .
Figure 8. Architecture of an auxiliary supply system.

Figure 8 .
Figure 8. Architecture of an auxiliary supply system.

Figure 9 .
Figure 9. Three-dimensional feature distribution for (a) AC source online and (b) AC source offline.

Figure 9 .
Figure 9. Three-dimensional feature distribution for (a) AC source online and (b) AC source offline.

Figure 11 .
Figure 11.Experimental results of voltage drop affected by the inrush current under different SOHs.

Figure 13 .
Figure 13.Three-dimensional feature distribution of experiment data samples.

Figure 11 .
Figure 11.Experimental results of voltage drop affected by the inrush current under different SOHs.

Figure 11 .
Figure 11.Experimental results of voltage drop affected by the inrush current under different SOHs.

Figure 13 .
Figure 13.Three-dimensional feature distribution of experiment data samples.

Figure 12 .
Figure 12.Detail energy for different wavelet decomposition layers.

Figure 12 .
Figure 12.Detail energy for different wavelet decomposition layers.

Figure 13 .
Figure 13.Three-dimensional feature distribution of experiment data samples.

Figure 13 .
Figure 13.Three-dimensional feature distribution of experiment data samples.

Table 1 .
Current Changing of Circuit Breaker used in the experiment.

Table 2 .
Main Parameters of Auxiliary DC Supply System.

Table 3 .
Current Changing of Circuit Breaker in auxiliary DC supply system.

Table 4 .
Battery SOH Range Setting in auxiliary DC supply system.

Table 2 .
Main Parameters of Auxiliary DC Supply System.

Table 3 .
Current Changing of Circuit Breaker in auxiliary DC supply system.

Table 4 .
Battery SOH Range Setting in auxiliary DC supply system.

Table 5 .
Test Set Accuracy in the Simulation.

Table 5 .
Test Set Accuracy in the Simulation.

Table 6 .
Main Parameters of System for Example II.

Table 7 .
Current Changing of Circuit Breaker for Example II.

Table 6 .
Main Parameters of System for Example II.

Table 7 .
Current Changing of Circuit Breaker for Example II.

Table 8 .
Experiment Battery SOH Range Setting.

Table 9 .
The Set Results in the Experiment.