CN110275114B - Storage battery internal resistance on-line monitoring method based on combined filtering algorithm - Google Patents

Storage battery internal resistance on-line monitoring method based on combined filtering algorithm Download PDF

Info

Publication number
CN110275114B
CN110275114B CN201910662857.2A CN201910662857A CN110275114B CN 110275114 B CN110275114 B CN 110275114B CN 201910662857 A CN201910662857 A CN 201910662857A CN 110275114 B CN110275114 B CN 110275114B
Authority
CN
China
Prior art keywords
storage battery
signal
internal resistance
formula
max
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910662857.2A
Other languages
Chinese (zh)
Other versions
CN110275114A (en
Inventor
刘林
丁笑迎
杨元浩
官晓鹏
王康
张溱旼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong Zhengchen Polytron Technologies Co ltd
Original Assignee
Shandong Zhengchen Polytron Technologies Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong Zhengchen Polytron Technologies Co ltd filed Critical Shandong Zhengchen Polytron Technologies Co ltd
Priority to CN201910662857.2A priority Critical patent/CN110275114B/en
Publication of CN110275114A publication Critical patent/CN110275114A/en
Application granted granted Critical
Publication of CN110275114B publication Critical patent/CN110275114B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/367Software therefor, e.g. for battery testing using modelling or look-up tables
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/36Arrangements for testing, measuring or monitoring the electrical condition of accumulators or electric batteries, e.g. capacity or state of charge [SoC]
    • G01R31/389Measuring internal impedance, internal conductance or related variables

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Tests Of Electric Status Of Batteries (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

The invention discloses an online monitoring method for internal resistance of a storage battery based on a combined filtering algorithm. The online monitoring method for the internal resistance of the storage battery filters transient pulse noise components in alternating current response signals by adopting a mathematical morphology alternating mixing filtering method, and extracts characteristic quantities of approximate signals by utilizing a wavelet threshold denoising algorithm, thereby not only filtering high-frequency noise and transient pulse interference in the alternating current response signals, but also overcoming the influence of line resistance on a detection result, and remarkably improving the online measurement precision of the internal resistance of the storage battery by an alternating current injection method.

Description

Storage battery internal resistance on-line monitoring method based on combined filtering algorithm
Technical Field
The invention relates to a storage battery internal resistance online monitoring method based on a combined filtering algorithm, in particular to a storage battery internal resistance online monitoring method based on a combined filtering algorithm, which is suitable for UPS power storage batteries and electric automobile power batteries in the field of green buildings.
Background
At present, the evaluation of the health state of the storage battery generally takes the change of the internal resistance of the storage battery as an important basis of the capacity and the residual life of the storage battery, and when the internal resistance value is increased by 20 percent compared with the standard internal resistance value and the capacity of the storage battery is less than 80 percent of the rated capacity, the storage battery generally has a fault and needs to be replaced immediately to ensure the normal operation of a direct current power supply system.
At present, a direct current instantaneous large current discharge method and an alternating current impedance method are generally adopted for internal resistance detection. Patent No. CN105092977A 'method and circuit for measuring internal resistance of accumulator, method and system for detecting health status', adopts direct current instantaneous heavy current discharge method to calculate internal resistance of accumulator by measuring instantaneous change of terminal voltage of accumulator and volt-ampere characteristic curve of accumulator when load is switched on or removed. The large load current is instantly connected, so that the influence of capacitance effect between the electrode plates of the storage battery on internal resistance measurement can be effectively avoided, but the service life and safety of the storage battery are adversely affected by the instant large current, and the storage battery is not suitable for the storage batteries with medium and small capacities. Patent No. CN104950180A 'EIS-based accurate measurement method of internal resistance of single storage battery' applies sinusoidal excitation current with fixed frequency to the storage battery by AC impedance method, and AC voltage drop signals are formed at two ends of the storage battery under the action of equivalent internal resistance of the storage battery, and the signals are used as AC response signals of the system, and the internal resistance value of the storage battery can be calculated by voltammetry after wavelet transformation. The method is not limited by the capacity of the storage battery, has a wide application range and supports online monitoring. Because the amplitude of the excitation signal is small, the generated ripple has small influence on a direct current load and a storage battery, but is influenced by interference such as high-frequency noise and transient pulse generated by loads such as an electromagnetic switch, a charger and an inverter, so that an alternating current response signal is subjected to wavelet transformation to cause local distortion, and the measurement precision is reduced. The simple wavelet transformation noise reduction method cannot meet the precision requirement of internal resistance measurement when the storage battery outputs power under the complex working condition. Therefore, how to effectively avoid the distortion influence of the high-strength pulse interference on the signal processing and improve the measurement precision of the internal resistance of the storage battery is a problem to be solved urgently at present. .
Disclosure of Invention
In order to overcome the defects of the technical problems, the invention provides an online monitoring method for the internal resistance of a storage battery based on a combined filtering algorithm.
The invention relates to a storage battery internal resistance on-line monitoring method based on a combined filtering algorithm1A sinusoidal excitation signal with amplitude A; a contrast sampling resistor R is connected in series on a loop between the excitation source and the storage battery1(ii) a The method is characterized in that the online monitoring method of the internal resistance of the storage battery comprises the following steps: firstly, collecting alternating current response signals at two ends of a storage battery and two ends of a comparison resistor, then respectively determining the scale expansion range of structural elements of two pairs of alternating current response signals, and then carrying out multi-scale alternate mixed filtering based on open-close operation on the alternating current response signals by utilizing the scale expansion range of the structural elements so as to eliminate transient pulse noise components in the alternating current response signals; then, after wavelet decomposition, threshold filtering processing of soft-hard compromise is carried out, effective values of the storage battery and the comparison resistance response signals are calculated, and finally, the internal resistance of the storage battery is calculated by adopting a comparison method.
The invention discloses a storage battery internal resistance online monitoring method based on a combined filtering algorithm, which is characterized by comprising the following steps of: the online monitoring method for the internal resistance of the storage battery is realized by the following steps:
a) AC signal acquisition at the same sampling frequency f under the condition that the excitation source generates a continuous and stable sinusoidal excitation signal2Collecting alternating current response signals f at two ends of storage batterybsw1(n) and a comparative sampling resistor R1Ac response signal f at both endsdsw2(N), counting the number of sampling points of the two alternating current response signals to be N;
b) determining the scale expansion range of the structural element, and collecting the AC discrete signal fbsw1(n)、fdsw2(n) performing the processes of the steps b-1) to b-3), respectively:
b-1), determining a maximum position sequence, determining a maximum in the sampling signal by searching and comparing the numerical values of sampling points in the discrete sampling signal, and determining a maximum position point set according to the sampling sequence:
Post-max={postmaxi|i=1,2,...,Nnum-max} (1)
in the formula, Post-maxIs a set of maximum location points, post-maxiIs the abscissa of the maximum value with the number i, Nost-maxThe number of maximum points;
b-2) calculating the distance between adjacent maximum values, and calculating the horizontal distance set D of the adjacent maximum values according to the abscissa of each maximum value obtained in the step b-1)p
Dp={dpm|dpm=post-max(m+1)-post-maxm,m=1,2,...,Nnum-max-1} (2)
In the formula (d)pmIs the horizontal distance between the maximum value with the serial number m +1 and the maximum value with the serial number m;
b-3) calculating the scale of the structural element, firstly, calculating the maximum value K of the scale of the structural element by using a formula (3)maxAnd a minimum value Kmin
Figure BDA0002139113460000031
In the formula, ceil [ ] is rounded up, floor [ ] is rounded down; from this, the set sg of dimensions of the structural elements is:
sg={sgi|sgi=Kmin,Kmin+1,…,Kmax} (4)
c) multi-scale alternate hybrid filtering, setting discrete sampling signal fbsw1(n)、fdsw2(n) after the processing of step b), the obtained scale sets of the structural elements are sg-sg respectivelyi|sgi=Kmin,Kmin+1,...,Kmax}、sg′={sg′i|gg′i=K′min,K′min+1,...,K′max}; adopting the multi-scale alternative mixed filtering algorithm based on the open-close operation given by the formula (5) to the sampling signal fbsw1(n)、fdsw2(n) filtering:
Figure BDA0002139113460000032
in the formula (f)bs(n)、fds(n) are each fbsw1(n)、fdsw2(n) represents a closed operation,
Figure BDA0002139113460000034
represents an on operation; wherein, the open operation mainly eliminates the positive pulse interference, and the close operation mainly eliminates the negative pulse interference;
d) wavelet decomposition with DB10Wavelet basis function pair morphological filtering signal fbs(n)、fds(n) performing 4-layer wavelet decomposition, wherein the expression after decomposition is as follows:
Figure BDA0002139113460000033
in the formula (I), the compound is shown in the specification,
Figure BDA0002139113460000041
for the i-th layer detail signal,
Figure BDA0002139113460000042
approximating the signal for the ith layer;
then using formula (7) to respectively process four layers of detail signals
Figure BDA0002139113460000043
Carrying out threshold filtering processing of soft and hard compromise:
Figure BDA0002139113460000044
in the formula, σ ∈ (0, 1)]In order to be the weight coefficient,
Figure BDA00021391134600000412
is the standard deviation of the detail signal of layer i, sgn [ [ alpha ] ]]Is a sign function, i is 1, 2, 3, 4, N' is Di(n) totality of discrete pointsCounting; di(n) is
Figure BDA0002139113460000045
Or
Figure BDA0002139113460000046
e) Signal reconstruction, setting four layers of detail signals
Figure BDA0002139113460000047
The signals processed by the soft and hard compromise threshold function of formula (7) are respectively marked as Wi、Wi', using WiAnd approximate signal
Figure BDA0002139113460000048
Reconstructing to obtain a response voltage signal U of the storage battery1(n) using Wi' sum approximation signal
Figure BDA0002139113460000049
Reconstructing to obtain a contrast voltage signal U2(n), i ═ 1, 2, 3, 4; then, the response voltage signal and the comparison voltage signal are subjected to effective value calculation using equation (8):
Figure BDA00021391134600000410
in the formula, N is the number of sampling points;
f) calculating internal resistance of the battery by using the effective value U of the battery response voltage signalp1Comparing the effective value U of the voltage signalp2Calculating equivalent internal resistance R of the storage battery by adopting a comparison methodac
Figure BDA00021391134600000411
In the formula, R1The resistance value of the comparison resistor is shown.
The method for monitoring the internal resistance of the storage battery on line based on the combined filtering algorithm comprises the steps of obtaining the internal resistance of the storage battery, judging whether the current internal resistance increment of the storage battery exceeds 20% of the standard internal resistance of the storage battery, giving a warning and a prompt for replacing the storage battery if the current internal resistance increment of the storage battery exceeds 20% of the standard internal resistance of the storage battery, and giving the internal resistance of the storage battery and the percentage exceeding the standard internal resistance if the current internal resistance increment of the.
The invention has the beneficial effects that: the invention relates to an on-line monitoring method for internal resistance of a storage battery, which comprises the steps of firstly collecting alternating current response signals at two ends of the storage battery and a comparison resistor under the condition of adding alternating current and sine excitation signals, then filtering transient pulse noise components in the alternating current response signals by adopting a mathematical morphology alternating mixing filtering method, extracting characteristic quantities of approximate signals by utilizing a wavelet threshold noise reduction algorithm, and finally calculating the internal resistance of the storage battery by adopting a comparison method.
Drawings
FIG. 1 is a circuit diagram of the online monitoring of internal resistance of a storage battery according to the present invention;
FIG. 2 is a signal processing flow diagram of the combined filtering algorithm of the present invention;
FIG. 3 is a flow chart of the method for on-line monitoring the internal resistance of the storage battery according to the present invention;
in fig. 4, a is a diagram of an ac response signal at two ends of an unprocessed storage battery, b is a diagram of a waveform of an ac signal at two ends of a storage battery after being processed by a simple wavelet threshold filtering process, and c is a diagram of a waveform of an ac signal at two ends of a storage battery after being processed by a combined filtering algorithm according to the present invention.
Detailed Description
The invention is further described with reference to the following figures and examples.
As shown in figure 1, the circuit diagram for on-line monitoring of internal resistance of the storage battery is provided, the lead-acid storage battery is a battery with internal resistance to be detected, the alternating current excitation signal generator is used for injecting alternating current sinusoidal signals into two ends of the storage battery, and a comparison sampling resistor R is connected in series on a line between the alternating current excitation signal generator and the storage battery1AC excitation signal generationThe equivalent resistance of the line between the device and the accumulator is R2. While injecting the AC excitation signal, the AC response signal collected from the two ends of the storage battery is fbsw1(n) the AC response signal collected from both ends of the comparison sampling resistor is fdsw2(n) of (a). Since the internal resistance of the battery is in milliohm level under normal conditions, the specific sampling resistance R is1Also in the milliohm range, e.g. 20 milliohms. The processing flow of the collected alternating current response signals at the two ends of the storage battery and the comparison resistor is as follows: firstly, according to the frequency and peak value self-setting structural element parameters of a sampling signal, a mathematical morphology alternating hybrid filter is used for conducting pre-filtering processing on the sampling signal, a wavelet threshold noise reduction function is used for conducting characteristic quantity extraction on an alternating current response signal of a monitoring system, and online accurate measurement on internal resistance of a lead-acid storage battery is achieved.
As shown in fig. 2, a signal processing flow chart of the combined filtering algorithm of the present invention is given, and fig. 3 is a flow chart of the online monitoring method for internal resistance of storage battery of the present invention, which is specifically realized by the following steps:
a) AC signal acquisition at the same sampling frequency f under the condition that the excitation source generates a continuous and stable sinusoidal excitation signal2Collecting alternating current response signals f at two ends of storage batterybsw1(n) and a comparative sampling resistor R1Ac response signal f at both endsdsw2(n), wherein n is the number of sampling points;
b) determining the scale expansion range of the structural element, and collecting the AC discrete signal fbsw1(n)、fdsw2(n) performing the treatments of the steps b-1) to b-3), respectively:
b-1), determining a maximum position sequence, determining a maximum in the sampling signal by searching and comparing the numerical values of sampling points in the discrete sampling signal, and determining a maximum position point set according to the sampling sequence:
Post-max={post-maxi|i=1,2,...,Nnum-max} (1)
in the formula, Post-maxIs a set of maximum location points, post-maxiIs a serial number iAbscissa of maximum value of, Nost-maxThe number of maximum points;
b-2) calculating the distance between adjacent maximum values, and calculating the horizontal distance set D of the adjacent maximum values according to the abscissa of each maximum value obtained in the step b-1)p
Dp={dpm|dpm=post-max(m+1)-post-maxm,m=1,2,...,Nnum-max-1} (2)
In the formula (d)pmIs the horizontal distance between the maximum value with the serial number m +1 and the maximum value with the serial number m;
b-3) calculating the scale of the structural element, firstly, calculating the maximum value K of the scale of the structural element by using a formula (3)maxAnd a minimum value Kmin
Figure BDA0002139113460000061
In the formula, ceil [ ] is rounded up, floor [ ] is rounded down; from this, the set sg of dimensions of the structural elements is:
sg={sgi|sgi=Kmin,Kmin+1,...,Kmax} (4)
c) multi-scale alternate hybrid filtering, setting discrete sampling signal fbsw1(n)、fdsw2(n) after the processing of step b), the obtained scale sets of the structural elements are sg-sg respectivelyi|sgi=Kmin,Kmin+1,...,Kmax}、sg′={sg′i|gg′i=K′min,K′min+1,...,K′max}; adopting the multi-scale alternative mixed filtering algorithm based on the open-close operation given by the formula (5) to the sampling signal fbsw1(n)、fdsw2(n) filtering:
Figure BDA0002139113460000071
in the formula (f)bs(n)、fds(n) are each fbsw1(n)、fdsw2(n) represents a closed operation,
Figure BDA00021391134600000713
represents an on operation; wherein, the open operation mainly eliminates the positive pulse interference, and the close operation mainly eliminates the negative pulse interference;
d) wavelet decomposition with DB10Wavelet basis function pair morphological filtering signal fbs(n)、fds(n) performing 4-layer wavelet decomposition, wherein the expression after decomposition is as follows:
Figure BDA0002139113460000072
in the formula (I), the compound is shown in the specification,
Figure BDA0002139113460000073
for the i-th layer detail signal,
Figure BDA0002139113460000074
approximating the signal for the ith layer;
then using formula (7) to respectively process four layers of detail signals
Figure BDA0002139113460000075
Carrying out threshold filtering processing of soft and hard compromise:
Figure BDA0002139113460000076
in the formula, σ ∈ (0, 1)]In order to be the weight coefficient,
Figure BDA0002139113460000077
is the standard deviation of the detail signal at layer j, sgn [ [ deg. ] ]]Is a sign function, i is 1, 2, 3, 4; di(n) is
Figure BDA0002139113460000078
Or
Figure BDA0002139113460000079
e) Signal reconstruction, setting four layers of detail signals
Figure BDA00021391134600000710
The signals processed by the soft and hard compromise threshold function of formula (7) are respectively marked as Wi、Wi', using WiAnd approximate signal
Figure BDA00021391134600000711
Reconstructing to obtain a response voltage signal U of the storage battery1(n) using Wi' sum approximation signal
Figure BDA00021391134600000712
Reconstructing to obtain a contrast voltage signal U2(n), i ═ 1, 2, 3, 4; then, the response voltage signal and the comparison voltage signal are subjected to effective value calculation using equation (8):
Figure BDA0002139113460000081
in the formula, N1Is U1(N) number of discrete voltage values contained, N2Is U2(n) the number of discrete voltage values contained;
f) calculating internal resistance of the battery by using the effective value U of the battery response voltage signalp1Comparing the effective value U of the voltage signalp2Calculating equivalent internal resistance R of the storage battery by adopting a comparison methodac
Figure BDA0002139113460000082
In the formula, R1The resistance value of the comparison resistor is shown.
After the internal resistance of the storage battery is obtained, whether the current internal resistance increment of the storage battery exceeds 20% of the standard internal resistance is judged, if yes, an alarm and a prompt for replacing the storage battery are given, and if not, the internal resistance of the storage battery and the percentage exceeding the standard internal resistance are given.
As shown in fig. 4, a is an unprocessed ac response signal at two ends of the storage battery, b is a waveform diagram of the ac signal at two ends of the storage battery after being processed only by the simple wavelet threshold filtering, and c is a waveform diagram of the ac signal at two ends of the storage battery after being processed by the combined filtering algorithm of the present invention. By comparing the waveform diagrams in the graph b and the graph a, it is found that the simple wavelet threshold filtering can only filter small positive and negative impulse noise in the signal, and is difficult to filter large transient impulse noise components in the alternating current response signal. By comparing the graph c with the graph a, the combined filtering algorithm can filter out high-frequency noise and transient impulse noise in the alternating current response signal, and is favorable for calculating more accurate internal resistance of the storage battery.
The storage battery internal resistance online monitoring method based on the combined filtering algorithm not only ensures the filtering capability of the morphological filter on transient pulse interference, but also ensures the filtering capability of the wavelet threshold algorithm on high-frequency noise, and fully exerts the respective characteristics of the morphological filter and the wavelet threshold algorithm. Compared with the wavelet threshold denoising algorithm for directly processing the alternating current response signal, the combined method reduces the influence of the transient pulse component on the peak value of the alternating current response signal and effectively improves the signal-to-noise ratio of the filtering signal. Meanwhile, the method effectively avoids measurement deviation caused by factors such as electromagnetic interference and the like under the background of a complex environment, and has important significance for detecting the internal resistance of the storage battery.

Claims (2)

1. A storage battery internal resistance on-line monitoring method based on a combined filtering algorithm is characterized in that excitation sources are connected to two ends of a storage battery, and the frequency generated by the excitation sources is f1A sinusoidal excitation signal with amplitude A; a contrast sampling resistor R is connected in series on a loop between the excitation source and the storage battery1(ii) a The method is characterized in that the online monitoring method of the internal resistance of the storage battery comprises the following steps: firstly, alternating current response signals at two ends of a storage battery and two ends of a comparison resistor are collected, then the scale expansion range of structural elements of two pairs of alternating current response signals is respectively determined, and then the scale expansion range of the structural elements is utilized to conduct opening-based operation on the alternating current response signals-closed-operation multiscale alternating mixing filtering to eliminate transient impulse noise components in the ac response signal; then, after wavelet decomposition, threshold filtering processing of soft-hard compromise is carried out, effective values of the storage battery and the comparison resistance response signals are calculated, and finally, the internal resistance of the storage battery is calculated by adopting a comparison method;
the online monitoring method for the internal resistance of the storage battery is realized by the following steps:
a) AC signal acquisition at the same sampling frequency f under the condition that the excitation source generates a continuous and stable sinusoidal excitation signal2Collecting alternating current response signals f at two ends of storage batterybsw1(n) and a comparative sampling resistor R1Ac response signal f at both endsdsw2(N), counting the number of sampling points of the two alternating current response signals to be N;
b) determining the scale expansion range of the structural element, and collecting the AC discrete signal fbsw1(n)、fdsw2(n) performing the treatments of the steps b-1) to b-3), respectively:
b-1), determining a maximum position sequence, determining a maximum in the sampling signal by searching and comparing the numerical values of sampling points in the discrete sampling signal, and determining a maximum position point set according to the sampling sequence:
Post-max={post-maxi|i=1,2,...,Nnum-max} (1)
in the formula, Post-maxIs a set of maximum location points, post-maxiIs the abscissa of the maximum value with the number i, Nnum-maxThe number of maximum points;
b-2) calculating the distance between adjacent maximum values, and calculating the horizontal distance set D of the adjacent maximum values according to the abscissa of each maximum value obtained in the step b-1)p
Dp={dpm|dpm=post-max(m+1)-post-maxm,m=1,2,...,Nnum-max-1} (2)
In the formula (d)pmIs the horizontal distance between the maximum value with the serial number m +1 and the maximum value with the serial number m;
b-3) solving structureThe element scale is obtained by first calculating the maximum value K of the structural element scale by using formula (3)maxAnd a minimum value Kmin
Figure FDA0003010540250000021
In the formula, ceil [ ] is rounded up, floor [ ] is rounded down; from this, the set sg of dimensions of the structural elements is:
sg={sgi|sgi=Kmin,Kmin+1,...,Kmax} (4)
c) multi-scale alternate hybrid filtering, setting discrete sampling signal fbsw1(n)、fdsw2(n) after the processing of step b), the obtained scale sets of the structural elements are sg-sg respectivelyi|sgi=Kmin,Kmin+1,...,Kmax}、sg′={sg′i|sg′i=K′min,K′min+1,...,K′max}; adopting the multi-scale alternative mixed filtering algorithm based on the open-close operation given by the formula (5) to the sampling signal fbsw1(n)、fdsw2(n) filtering:
Figure FDA0003010540250000022
in the formula (f)bs(n)、fds(n) are each fbsw1(n)、fdsw2(n) represents a closed operation,
Figure FDA0003010540250000027
represents an on operation; wherein, the open operation mainly eliminates the positive pulse interference, and the close operation mainly eliminates the negative pulse interference;
d) wavelet decomposition with DB10Wavelet basis function pair morphological filtering signal fbs(n)、fds(n) performing 4-layer wavelet decomposition, wherein the expression after decomposition is as follows:
Figure FDA0003010540250000023
in the formula (I), the compound is shown in the specification,
Figure FDA0003010540250000024
for the i-th layer detail signal,
Figure FDA0003010540250000025
approximating the signal for the ith layer;
then using formula (7) to respectively process four layers of detail signals
Figure FDA0003010540250000026
Carrying out threshold filtering processing of soft and hard compromise:
Figure FDA0003010540250000031
in the formula, σ ∈ (0, 1)]In order to be the weight coefficient,
Figure FDA0003010540250000032
is the standard deviation of the detail signal of layer i, sgn [ [ alpha ] ]]Is a sign function, i is 1, 2, 3, 4, N' is Di(n) total number of discrete points; di(n) is
Figure FDA0003010540250000033
Or
Figure FDA0003010540250000034
e) Signal reconstruction, setting four layers of detail signals
Figure FDA0003010540250000035
The signals processed by the soft and hard compromise threshold function of formula (7) are respectively marked as Wi、Wi', using WiAnd approximate signal
Figure FDA0003010540250000036
Reconstructing to obtain a response voltage signal U of the storage battery1(n) using Wi' sum approximation signal
Figure FDA0003010540250000037
Reconstructing to obtain a contrast voltage signal U2(n), i ═ 1, 2, 3, 4; then, the response voltage signal and the comparison voltage signal are subjected to effective value calculation using equation (8):
Figure FDA0003010540250000038
in the formula, N is the number of sampling points;
f) calculating internal resistance of the battery by using the effective value U of the battery response voltage signalp1Comparing the effective value U of the voltage signalp2Calculating equivalent internal resistance R of the storage battery by adopting a comparison methodac
Figure FDA0003010540250000039
In the formula, R1The resistance value of the comparison resistor is shown.
2. The combined filtering algorithm-based storage battery internal resistance online monitoring method according to claim 1, characterized in that: after the internal resistance of the storage battery is obtained, whether the current internal resistance increment of the storage battery exceeds 20% of the standard internal resistance is judged, if yes, an alarm and a prompt for replacing the storage battery are given, and if not, the internal resistance of the storage battery and the percentage exceeding the standard internal resistance are given.
CN201910662857.2A 2019-07-22 2019-07-22 Storage battery internal resistance on-line monitoring method based on combined filtering algorithm Active CN110275114B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910662857.2A CN110275114B (en) 2019-07-22 2019-07-22 Storage battery internal resistance on-line monitoring method based on combined filtering algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910662857.2A CN110275114B (en) 2019-07-22 2019-07-22 Storage battery internal resistance on-line monitoring method based on combined filtering algorithm

Publications (2)

Publication Number Publication Date
CN110275114A CN110275114A (en) 2019-09-24
CN110275114B true CN110275114B (en) 2021-06-25

Family

ID=67965065

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910662857.2A Active CN110275114B (en) 2019-07-22 2019-07-22 Storage battery internal resistance on-line monitoring method based on combined filtering algorithm

Country Status (1)

Country Link
CN (1) CN110275114B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110988232A (en) * 2019-11-28 2020-04-10 南通乐尔环保科技有限公司 Chromatographic baseline noise reduction method and device
CN112269067A (en) * 2020-10-22 2021-01-26 国网天津市电力公司 Capacitor dielectric loss angle on-line monitoring method based on morphology and 3 sigma criterion
CN112732684B (en) * 2020-12-10 2023-04-14 中国西安卫星测控中心 Denoising method and device for current data of satellite storage battery
CN114397581A (en) * 2021-12-09 2022-04-26 国网天津市电力公司 New energy automobile battery SOC anti-interference assessment method oriented to charging monitoring data of direct-current charging pile

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101404408A (en) * 2008-11-13 2009-04-08 昆明理工大学 Distribution network cable-wire mixed line failure route selection method by utilizing relative energy
CN102435897A (en) * 2011-09-14 2012-05-02 浙江省电力试验研究院 Method for filtering online monitoring data of arrester leakage current based on morphological principle
CN102508013A (en) * 2011-11-25 2012-06-20 国网技术学院 Current alarm signal processing method for improved mathematical morphology filter
CN102736000A (en) * 2012-07-06 2012-10-17 桂林电子科技大学 On-line detection method and device for power capacitor
CN102944818A (en) * 2012-11-27 2013-02-27 广东电网公司佛山供电局 Geographic information system (GIS) terminal cable fault on-line location method
CN103454537A (en) * 2013-09-16 2013-12-18 国家电网公司 Wind power generation low-voltage ride-through detection equipment and method based on wavelet analysis
CN103617630A (en) * 2013-12-18 2014-03-05 国网重庆潼南县供电有限责任公司 Infrared image dividing method and system for power system equipment based on wavelet analysis
CN104367316A (en) * 2014-11-13 2015-02-25 重庆邮电大学 Electrocardiosignal denoising method based on morphological filtering and lifting wavelet transformation
CN104714076A (en) * 2015-03-24 2015-06-17 清华大学 Method and device for smoothing current signals
CN104950180A (en) * 2015-06-02 2015-09-30 国电南京自动化股份有限公司 Method for precisely measuring internal resistance of single storage batteries on basis of EIS (electrochemical impedance spectroscopy)
CN105913835A (en) * 2016-06-15 2016-08-31 华北电力大学 Self-adaptive filtering method based on mathematical morphology
CN106019076A (en) * 2016-05-16 2016-10-12 西安理工大学 Fault distance detecting method for high-voltage DC transmission line
CN106054028A (en) * 2016-07-29 2016-10-26 咸亨国际(杭州)电气制造有限公司 Cable fault automatic ranging method based on time domain characteristics and wavelet analysis
CN106780148A (en) * 2016-12-30 2017-05-31 华北电力大学 Adaptive morphology filter method based on composite construction element
CN106845334A (en) * 2016-11-29 2017-06-13 国网辽宁省电力有限公司沈阳供电公司 A kind of innovative noise extracting method based on mathematical morphology
KR20180012891A (en) * 2016-07-27 2018-02-07 주식회사 라이프사이언스테크놀로지 Method for remove noise of Bio signal

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11029347B2 (en) * 2017-04-26 2021-06-08 Nokomis, Inc Electronics equipment testing apparatus and method utilizing unintended RF emission features

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101404408A (en) * 2008-11-13 2009-04-08 昆明理工大学 Distribution network cable-wire mixed line failure route selection method by utilizing relative energy
CN102435897A (en) * 2011-09-14 2012-05-02 浙江省电力试验研究院 Method for filtering online monitoring data of arrester leakage current based on morphological principle
CN102508013A (en) * 2011-11-25 2012-06-20 国网技术学院 Current alarm signal processing method for improved mathematical morphology filter
CN102736000A (en) * 2012-07-06 2012-10-17 桂林电子科技大学 On-line detection method and device for power capacitor
CN102944818A (en) * 2012-11-27 2013-02-27 广东电网公司佛山供电局 Geographic information system (GIS) terminal cable fault on-line location method
CN103454537A (en) * 2013-09-16 2013-12-18 国家电网公司 Wind power generation low-voltage ride-through detection equipment and method based on wavelet analysis
CN103617630A (en) * 2013-12-18 2014-03-05 国网重庆潼南县供电有限责任公司 Infrared image dividing method and system for power system equipment based on wavelet analysis
CN104367316A (en) * 2014-11-13 2015-02-25 重庆邮电大学 Electrocardiosignal denoising method based on morphological filtering and lifting wavelet transformation
CN104714076A (en) * 2015-03-24 2015-06-17 清华大学 Method and device for smoothing current signals
CN104950180A (en) * 2015-06-02 2015-09-30 国电南京自动化股份有限公司 Method for precisely measuring internal resistance of single storage batteries on basis of EIS (electrochemical impedance spectroscopy)
CN106019076A (en) * 2016-05-16 2016-10-12 西安理工大学 Fault distance detecting method for high-voltage DC transmission line
CN105913835A (en) * 2016-06-15 2016-08-31 华北电力大学 Self-adaptive filtering method based on mathematical morphology
KR20180012891A (en) * 2016-07-27 2018-02-07 주식회사 라이프사이언스테크놀로지 Method for remove noise of Bio signal
CN106054028A (en) * 2016-07-29 2016-10-26 咸亨国际(杭州)电气制造有限公司 Cable fault automatic ranging method based on time domain characteristics and wavelet analysis
CN106845334A (en) * 2016-11-29 2017-06-13 国网辽宁省电力有限公司沈阳供电公司 A kind of innovative noise extracting method based on mathematical morphology
CN106780148A (en) * 2016-12-30 2017-05-31 华北电力大学 Adaptive morphology filter method based on composite construction element

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
融合形态学-复小波变换的PD干扰信号抑制方法;王友元 等;《仪器仪表学报》;20130315;707-713 *

Also Published As

Publication number Publication date
CN110275114A (en) 2019-09-24

Similar Documents

Publication Publication Date Title
CN110275114B (en) Storage battery internal resistance on-line monitoring method based on combined filtering algorithm
CN102269798B (en) Method and device for detecting faulty battery
CN114430080B (en) Power battery cell abnormal self-discharge identification method based on operation data
CN111505532A (en) Online detection method for early internal short circuit of series lithium battery pack based on SOC correlation coefficient
CN112924884B (en) Quantitative diagnosis method for short circuit in battery based on incremental capacity curve peak area
JP5803053B2 (en) Power storage device state estimation method and state estimation device
CN113311239B (en) Electric automobile insulation resistance detection circuit and method
CN106405424B (en) Method and device for metering residual electric quantity of lithium ion battery
CN108828448A (en) Battery charge state estimation on line method based on charging voltage curve fusion Kalman filtering
CN111123136A (en) Online measuring device and online measuring method for health state of lithium ion battery
Chen et al. An approach for state of charge estimation of Li-ion battery based on Thevenin equivalent circuit model
CN110554328B (en) HHT-based storage battery internal resistance measurement method and system
Huang et al. An approach to measurements of electrical characteristics of lithium-ion battery with open-circuit voltage function
Deng et al. An improved algorithm of SOC testing based on open-circuit voltage-ampere hour method
CN115840121A (en) Insulation detection circuit of energy storage high-voltage battery system and detection method thereof
CN115923592A (en) Battery state detection method and device for electric vehicle
Tang et al. An aging-and load-insensitive method for quantitatively detecting the battery internal-short-circuit resistance
CN112379271B (en) Passivation-considered capacity detection method for carbon-coated lithium thionyl chloride battery
CN115469226A (en) Real-time safety early warning method for electric vehicle power battery based on operation big data
CN113376435B (en) Insulation monitoring circuit and method and electric automobile
CN114114024A (en) Storage battery consistency screening device and screening method
CN111308366A (en) Storage battery pack multi-parameter monitoring system and using method thereof
CN216526140U (en) Insulation detection system for ship power battery pack
Hossain et al. Performance analysis of filter sensing board for measuring the battery online impedance
Grcić et al. Electric vehicle charging station fault detection: a machine learning approach

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant