[go: up one dir, main page]

CN115267333A - Lithium battery broadband impedance spectrum testing method based on pseudorandom sequence - Google Patents

Lithium battery broadband impedance spectrum testing method based on pseudorandom sequence Download PDF

Info

Publication number
CN115267333A
CN115267333A CN202210904659.4A CN202210904659A CN115267333A CN 115267333 A CN115267333 A CN 115267333A CN 202210904659 A CN202210904659 A CN 202210904659A CN 115267333 A CN115267333 A CN 115267333A
Authority
CN
China
Prior art keywords
impedance
lithium battery
pseudo
frequency
random sequence
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.)
Granted
Application number
CN202210904659.4A
Other languages
Chinese (zh)
Other versions
CN115267333B (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.)
Nanjing Institute of Technology
Original Assignee
Nanjing Institute of Technology
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 Nanjing Institute of Technology filed Critical Nanjing Institute of Technology
Priority to CN202210904659.4A priority Critical patent/CN115267333B/en
Publication of CN115267333A publication Critical patent/CN115267333A/en
Application granted granted Critical
Publication of CN115267333B publication Critical patent/CN115267333B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • 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/385Arrangements for measuring battery or accumulator variables
    • 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
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H21/00Adaptive networks
    • H03H21/0012Digital adaptive filters
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/10Energy storage using batteries

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明公开了一种基于伪随机序列的锂电池宽频阻抗谱测试方法,属于动力锂电池应用技术领域。该锂电池宽频阻抗谱测试方法包括如下步骤:设计一段用于宽频阻抗测量的伪随机序列;通过电池管理系统向锂电池注入伪随机序列,采集注入过程中的锂电池端电压和电流,并且计算锂电池在不同频率下的量测阻抗;在频域内通过自适应双边滤波器对量测阻抗进行滤波处理,得到宽频阻抗。通过对锂电池注入一种伪随机序列,通过双边滤波算法对阻抗测量结果进行滤波处理,降低噪声干扰,从而提高动力锂电池宽频阻抗测量的准确性与稳定性。

Figure 202210904659

The invention discloses a method for testing lithium battery broadband impedance spectrum based on pseudo-random sequence, and belongs to the technical field of power lithium battery application. The lithium battery broadband impedance spectrum testing method includes the following steps: designing a pseudo-random sequence for broadband impedance measurement; injecting the pseudo-random sequence into the lithium battery through the battery management system, collecting the terminal voltage and current of the lithium battery during the injection process, and calculating the lithium battery The measured impedance of the battery at different frequencies; the measured impedance is filtered through an adaptive bilateral filter in the frequency domain to obtain a broadband impedance. By injecting a pseudo-random sequence into the lithium battery, the impedance measurement results are filtered through a bilateral filtering algorithm to reduce noise interference, thereby improving the accuracy and stability of broadband impedance measurement of power lithium batteries.

Figure 202210904659

Description

一种基于伪随机序列的锂电池宽频阻抗谱测试方法A testing method of lithium battery broadband impedance spectroscopy based on pseudo-random sequence

技术领域technical field

本发明属于动力锂电池应用技术领域,具体地,涉及一种基于伪随机序列的锂电池宽频阻抗谱测试方法。The invention belongs to the technical field of power lithium battery application, and in particular relates to a lithium battery broadband impedance spectrum testing method based on a pseudo-random sequence.

背景技术Background technique

锂离子电池具有高功率密度、长循环寿命、无记忆效应和低自放电率等优点,在电力系统、新能源汽车和各类可携带电子器件中受到了大量的关注。为确保锂离子电池能够在安全和可靠的环境下工作,往往需要通过电池管理系统来对电池的状态进行实时监控。其中,电池阻抗与电池的荷电状态、老化状态以及温度等密切相关,通过电化学阻抗谱能够实时对电池的状态进行估计以及故障诊断,但是需要对不同频率下的电池阻抗进行快速、准确的测量。Lithium-ion batteries have the advantages of high power density, long cycle life, no memory effect and low self-discharge rate, and have received a lot of attention in power systems, new energy vehicles and various portable electronic devices. In order to ensure that the lithium-ion battery can work in a safe and reliable environment, it is often necessary to monitor the state of the battery in real time through the battery management system. Among them, the battery impedance is closely related to the state of charge, aging state, and temperature of the battery. Electrochemical impedance spectroscopy can be used to estimate the state of the battery and diagnose faults in real time, but it is necessary to quickly and accurately measure the battery impedance at different frequencies. Measurement.

常见的阻抗谱测试方法是频率响应分析方法,该方法通过正弦信号扫频对宽频阻抗进行测量,能够提供足够准确可靠测量结果,但是前提是需要冗长的时间。通过将正弦信号叠加能够大幅度缩短阻抗谱测试时间,但是正弦叠加信号的发波需要复杂的硬件设计,这极大地限制了该方法的实际用途。采用方波信号测量阻抗谱比较容易实施,但是方波信号包含的频域信息极不充分,因此难以获取完整的阻抗谱信息。相比之下,伪随机二进制序列含有充分的频域谐波,该信号含有有限的信号幅值,能够有效地对硬件设计进行简化,因此特别适用于宽频阻抗的测量。The common impedance spectrum test method is the frequency response analysis method, which measures the broadband impedance by sweeping the sinusoidal signal, which can provide sufficiently accurate and reliable measurement results, but the premise is that it takes a long time. By superimposing sinusoidal signals, the impedance spectrum test time can be greatly shortened, but the wave generation of sinusoidal superposition signals requires complex hardware design, which greatly limits the practical use of this method. Using square wave signal to measure impedance spectrum is relatively easy to implement, but the frequency domain information contained in square wave signal is extremely insufficient, so it is difficult to obtain complete impedance spectrum information. In contrast, the pseudo-random binary sequence contains sufficient harmonics in the frequency domain, and the signal has limited signal amplitude, which can effectively simplify the hardware design, so it is especially suitable for the measurement of broadband impedance.

伪随机序列频域信号的功率含量低,容易受到量测噪声的干扰。通过设计阻抗滤波器能够一定程度上滤除噪声对信号的干扰。传统的滑动均值滤波器虽然简单,但是对于滤除方差较大的噪声性能不足。随着信号功率及信噪比的下降,阻抗滤波结果非常容易发生偏差。由于频率分辨率的限制,低频范围内的可用数据及其稀少,这进一步限制了阻抗测量的精度。目前仍然缺乏有效可靠的手段针对电池宽频阻抗进行在线测量,传统的伪随机序列性能有待进一步提升,电池阻抗测量的精准度及可靠性有待加强。Pseudo-random sequence frequency domain signals have low power content and are easily interfered by measurement noise. The interference of noise to the signal can be filtered out to a certain extent by designing the impedance filter. Although the traditional moving average filter is simple, its performance is insufficient for filtering out noise with large variance. As the signal power and signal-to-noise ratio decrease, the results of impedance filtering are very prone to deviation. Due to limited frequency resolution, the availability of data in the low frequency range is extremely scarce, which further limits the accuracy of impedance measurements. At present, there is still a lack of effective and reliable means for online measurement of battery broadband impedance. The performance of the traditional pseudo-random sequence needs to be further improved, and the accuracy and reliability of battery impedance measurement need to be strengthened.

发明内容Contents of the invention

针对现有技术的不足,本发明提供了一种基于伪随机序列的锂电池宽频阻抗谱测试方法,通过对锂电池注入一种伪随机序列,通过双边滤波算法对阻抗测量结果进行滤波处理,降低噪声干扰,从而提高动力锂电池宽频阻抗测量的准确性与稳定性。Aiming at the deficiencies of the prior art, the present invention provides a lithium battery broadband impedance spectrum testing method based on a pseudo-random sequence. By injecting a pseudo-random sequence into the lithium battery, the impedance measurement result is filtered by a bilateral filtering algorithm, reducing the Noise interference, thereby improving the accuracy and stability of broadband impedance measurement of power lithium batteries.

为了实现上述目的,本发明采用如下技术方案:一种基于伪随机序列的锂电池宽频阻抗谱测试方法,具体包括如下步骤:In order to achieve the above object, the present invention adopts the following technical scheme: a lithium battery broadband impedance spectrum testing method based on a pseudo-random sequence, specifically comprising the following steps:

步骤1、设计一段用于宽频阻抗测量的伪随机序列;Step 1. Design a pseudo-random sequence for broadband impedance measurement;

步骤2、通过电池管理系统向锂电池注入伪随机序列,采集注入过程中的锂电池端电压和电流,并且计算锂电池在不同频率下的量测阻抗;Step 2. Inject a pseudo-random sequence into the lithium battery through the battery management system, collect the terminal voltage and current of the lithium battery during the injection process, and calculate the measured impedance of the lithium battery at different frequencies;

步骤3、在频域内通过自适应双边滤波器对量测阻抗进行滤波处理,得到宽频阻抗。Step 3. Filtering the measured impedance through an adaptive bilateral filter in the frequency domain to obtain broadband impedance.

进一步地,所述伪随机序列为双重伪随机序列或混合伪随机序列。Further, the pseudo-random sequence is a double pseudo-random sequence or a mixed pseudo-random sequence.

进一步地,所述混合伪随机序列SH由第一子序列SA和第二子序列SB叠加而成:SH=0.5(SA+SB),其中第一子序列SA的长度为N1,注入频率为f1,第二子序列SB的长度为N2,注入频率为f2,满足:

Figure BDA0003771819100000021
Further, the mixed pseudo-random sequence SH is formed by superimposing the first subsequence S A and the second subsequence S B : S H =0.5( SA +S B ), wherein the length of the first subsequence S A is N 1 , the injection frequency is f 1 , the length of the second subsequence S B is N 2 , and the injection frequency is f 2 , satisfying:
Figure BDA0003771819100000021

进一步地,所述双重伪随机序列由第三子序列和第四子序列组成,所述第三子序列注入频率为f3和第四子序列的注入频率f4满足:f3=1.5f4Further, the double pseudo-random sequence is composed of a third subsequence and a fourth subsequence, the injection frequency of the third subsequence is f 3 and the injection frequency f 4 of the fourth subsequence satisfies: f 3 =1.5f 4 .

进一步地,所述锂电池在不同频率下的量测阻抗

Figure BDA0003771819100000022
为:Further, the measured impedance of the lithium battery at different frequencies
Figure BDA0003771819100000022
for:

Figure BDA0003771819100000023
Figure BDA0003771819100000023

其中,V(f)为在阻抗频率f下的电压谐波幅值,I(f)为在阻抗频率f下的电流谐波幅值。Among them, V(f) is the voltage harmonic amplitude at impedance frequency f, and I(f) is the current harmonic amplitude at impedance frequency f.

进一步地,所述宽频阻抗

Figure BDA0003771819100000024
为:Further, the broadband impedance
Figure BDA0003771819100000024
for:

Figure BDA0003771819100000025
Figure BDA0003771819100000025

其中,f为阻抗频率,fi-1为滤波窗口的下界,fi+1为滤波窗口的上界,

Figure BDA0003771819100000026
为锂电池在阻抗频率f下的阻抗,g1(f)为关于阻抗频率的核函数,
Figure BDA0003771819100000027
为关于量测阻抗的核函数。Among them, f is the impedance frequency, f i-1 is the lower bound of the filtering window, and f i+1 is the upper bound of the filtering window,
Figure BDA0003771819100000026
is the impedance of the lithium battery at the impedance frequency f, g 1 (f) is the kernel function about the impedance frequency,
Figure BDA0003771819100000027
is the kernel function about the measured impedance.

进一步地,所述关于阻抗频率的核函数g1(f)、关于量测阻抗的核函数

Figure BDA0003771819100000028
分别为:Further, the kernel function g 1 (f) about the impedance frequency, the kernel function about the measured impedance
Figure BDA0003771819100000028
They are:

Figure BDA0003771819100000029
Figure BDA0003771819100000029

Figure BDA00037718191000000210
Figure BDA00037718191000000210

其中,σf为关于频率的标准差,μf为关于频率的均值,μf=f',f'为滤波窗口的中点;

Figure BDA0003771819100000031
关于实部阻抗的标准差,
Figure BDA0003771819100000032
wf为加权系数;
Figure BDA0003771819100000033
为关于虚部阻抗的标准差,
Figure BDA0003771819100000034
Figure BDA0003771819100000035
为关于实部阻抗的均值,
Figure BDA0003771819100000036
Figure BDA0003771819100000037
为锂电池在阻抗频率f下的阻抗的实部,
Figure BDA0003771819100000038
为锂电池在阻抗频率f下的阻抗的虚部,
Figure BDA0003771819100000039
为关于虚部阻抗的均值。Wherein, σ f is the standard deviation about frequency, μ f is the mean value about frequency, μ f =f', and f' is the midpoint of the filtering window;
Figure BDA0003771819100000031
With respect to the standard deviation of the real impedance,
Figure BDA0003771819100000032
w f is the weighting coefficient;
Figure BDA0003771819100000033
is the standard deviation about the imaginary part impedance,
Figure BDA0003771819100000034
Figure BDA0003771819100000035
is the mean value with respect to the real impedance,
Figure BDA0003771819100000036
Figure BDA0003771819100000037
is the real part of the impedance of the lithium battery at the impedance frequency f,
Figure BDA0003771819100000038
is the imaginary part of the impedance of the lithium battery at the impedance frequency f,
Figure BDA0003771819100000039
is the mean value with respect to the imaginary part impedance.

进一步地,所述加权系数wf通过电流谐波幅值I(f)求得:Further, the weighting coefficient w f is obtained by the current harmonic amplitude I(f):

Figure BDA00037718191000000310
Figure BDA00037718191000000310

进一步地,当伪随机序列为双重伪随机序列时,σf=(fi+1-fi-1)/6。Further, when the pseudo-random sequence is a double pseudo-random sequence, σ f =(f i+1 -f i-1 )/6.

进一步地,当伪随机序列为混合伪随机序列时,

Figure BDA00037718191000000311
Further, when the pseudo-random sequence is a mixed pseudo-random sequence,
Figure BDA00037718191000000311

与现有技术相比,本发明具有如下有益效果:本发明基于伪随机序列的锂电池宽频阻抗谱测试方法通过注入一段混合伪随机序列或双重伪随机序列,并采用自适应双边滤波器对测量结果进行滤波处理,从而实现对锂电池的宽频阻抗进行快速测量。此外,本发明所提出的混合伪随机序列或双重伪随机序列能够增强低频范围内信号的信噪比,具有足够高的鲁棒性,相比传统的阻抗测量方案,本发明锂电池宽频阻抗谱测试方法能够快速精确地对锂电池阻抗进行测量,阻抗测量结果精确而稳定,适用于锂电池的各类应用场景。Compared with the prior art, the present invention has the following beneficial effects: the lithium battery broadband impedance spectrum testing method based on the pseudo-random sequence of the present invention injects a mixed pseudo-random sequence or double pseudo-random sequence, and uses an adaptive bilateral filter to measure The result is filtered, so as to realize the fast measurement of the broadband impedance of the lithium battery. In addition, the mixed pseudo-random sequence or double pseudo-random sequence proposed by the present invention can enhance the signal-to-noise ratio of the signal in the low-frequency range, and has sufficiently high robustness. Compared with the traditional impedance measurement scheme, the lithium battery broadband impedance spectrum of the present invention The test method can quickly and accurately measure the impedance of lithium batteries, and the impedance measurement results are accurate and stable, which is suitable for various application scenarios of lithium batteries.

附图说明Description of drawings

附图是用来提供对本发明的进一步理解,并且构成说明书的一部分,与具体实施方式一起用于解释本发明,但并不构成对本发明的限制。The accompanying drawings are used to provide a further understanding of the present invention, and constitute a part of the description, and are used together with the specific implementation to explain the present invention, but do not constitute a limitation to the present invention.

图1是本发明基于伪随机序列的锂电池宽频阻抗谱测试方法的流程图;Fig. 1 is the flow chart of the lithium battery broadband impedance spectrum test method based on the pseudo-random sequence of the present invention;

图2是实施例1在环境温度25℃、荷电状态为80%下的宽频阻抗量测图;Fig. 2 is the broadband impedance measurement figure under the environment temperature 25 ℃, state of charge of embodiment 1 is 80%;

图3是实施例1在环境温度25℃、荷电状态为20%下的宽频阻抗量测图;Fig. 3 is a broadband impedance measurement diagram of embodiment 1 at an ambient temperature of 25°C and a state of charge of 20%;

图4是实施例1在环境温度35℃、荷电状态为50%下的宽频阻抗量测图;Fig. 4 is a broadband impedance measurement diagram of embodiment 1 at an ambient temperature of 35°C and a state of charge of 50%;

图5是实施例1在环境温度15℃、荷电状态为50%下的宽频阻抗量测图;Fig. 5 is a broadband impedance measurement diagram of Example 1 at an ambient temperature of 15°C and a state of charge of 50%;

图6是实施例2在环境温度25℃、荷电状态为80%下的宽频阻抗量测图;Fig. 6 is a broadband impedance measurement diagram of Example 2 at an ambient temperature of 25°C and a state of charge of 80%;

图7是实施例2在环境温度25℃、荷电状态为20%下的宽频阻抗量测图;Fig. 7 is a broadband impedance measurement diagram of Example 2 at an ambient temperature of 25°C and a state of charge of 20%;

图8是实施例2在环境温度35℃、荷电状态为50%下的宽频阻抗量测图;Fig. 8 is a broadband impedance measurement diagram of Example 2 at an ambient temperature of 35°C and a state of charge of 50%;

图9是实施例2在环境温度15℃、荷电状态为50%下的宽频阻抗量测图。FIG. 9 is a measurement diagram of broadband impedance of Example 2 at an ambient temperature of 15° C. and a state of charge of 50%.

具体实施方式Detailed ways

以下结合附图和实施例对本发明的技术方案进行详细说明。应当理解的是,此处所描述的具体实施方式仅用于说明和解释本发明实施例,并不用于限制本发明实施例。The technical solutions of the present invention will be described in detail below in conjunction with the drawings and embodiments. It should be understood that the specific implementation manners described here are only used to illustrate and explain the embodiments of the present invention, and are not intended to limit the embodiments of the present invention.

如图1为本发明基于伪随机序列的锂电池宽频阻抗谱测试方法的流程图,该锂电池宽频阻抗谱测试方法具体包括如下步骤:Figure 1 is a flow chart of the lithium battery broadband impedance spectrum testing method based on a pseudo-random sequence in the present invention, the lithium battery broadband impedance spectrum testing method specifically includes the following steps:

步骤1、设计一段用于宽频阻抗测量的伪随机序列;本发明中伪随机序列为双重伪随机序列或混合伪随机序列。Step 1. Design a pseudo-random sequence for broadband impedance measurement; the pseudo-random sequence in the present invention is a double pseudo-random sequence or a mixed pseudo-random sequence.

混合伪随机序列能够提高频域内低频谐波功率含量,从而能够更好的对电池的低频阻抗进行测量;混合伪随机序列SH由第一子序列SA和第二子序列SB叠加而成:SH=0.5(SA+SB),其中第一子序列SA的长度为N1,注入频率为f1,第二子序列SB的长度为N2,注入频率为f2,满足:

Figure BDA0003771819100000041
The mixed pseudo-random sequence can improve the low-frequency harmonic power content in the frequency domain, so that the low-frequency impedance of the battery can be better measured; the mixed pseudo-random sequence SH is formed by superimposing the first subsequence S A and the second subsequence S B : S H =0.5( SA +S B ) , wherein the length of the first subsequence S A is N 1 , the injection frequency is f 1 , the length of the second subsequence S B is N 2 , and the injection frequency is f 2 , Satisfy:
Figure BDA0003771819100000041

双重伪随机序列能够提高频域内高频谐波功率含量,从而能够更好的对电池的高频阻抗进行测量;双重伪随机序列由第三子序列和第四子序列组成,第三子序列长度为N3,注入频率为f3和第四子序列的长度为N4,注入频率f4满足:f3=1.5f4,N3>N4The double pseudo-random sequence can improve the high-frequency harmonic power content in the frequency domain, so that the high-frequency impedance of the battery can be better measured; the double pseudo-random sequence is composed of the third subsequence and the fourth subsequence, and the length of the third subsequence is N 3 , the injection frequency is f 3 and the length of the fourth subsequence is N 4 , the injection frequency f 4 satisfies: f 3 =1.5f 4 , N 3 >N 4 .

步骤2、通过电池管理系统向锂电池注入伪随机序列,采集注入过程中的锂电池端电压和电流,并且计算锂电池在不同频率下的量测阻抗

Figure BDA0003771819100000042
Step 2. Inject a pseudo-random sequence into the lithium battery through the battery management system, collect the terminal voltage and current of the lithium battery during the injection process, and calculate the measured impedance of the lithium battery at different frequencies
Figure BDA0003771819100000042

Figure BDA0003771819100000043
Figure BDA0003771819100000043

其中,V(f)为在阻抗频率f下的电压谐波幅值,I(f)为在阻抗频率f下的电流谐波幅值。Among them, V(f) is the voltage harmonic amplitude at impedance frequency f, and I(f) is the current harmonic amplitude at impedance frequency f.

步骤3、为了减小量测噪声以及其他外部干扰对阻抗测量的影响,在频域内通过自适应双边滤波器对量测阻抗进行滤波处理,得到宽频阻抗

Figure BDA0003771819100000044
具有计算量小且滤波精度高的特点。Step 3. In order to reduce the impact of measurement noise and other external disturbances on impedance measurement, the measured impedance is filtered through an adaptive bilateral filter in the frequency domain to obtain broadband impedance
Figure BDA0003771819100000044
It has the characteristics of small calculation amount and high filtering precision.

Figure BDA0003771819100000051
Figure BDA0003771819100000051

其中,f为阻抗频率,fi-1为滤波窗口的下界,fi+1为滤波窗口的上界,

Figure BDA0003771819100000052
为锂电池在阻抗频率f下的阻抗,g1(f)为关于阻抗频率的核函数,
Figure BDA0003771819100000053
为关于量测阻抗的核函数。Among them, f is the impedance frequency, f i-1 is the lower bound of the filtering window, and f i+1 is the upper bound of the filtering window,
Figure BDA0003771819100000052
is the impedance of the lithium battery at the impedance frequency f, g 1 (f) is the kernel function about the impedance frequency,
Figure BDA0003771819100000053
is the kernel function about the measured impedance.

本发明中:

Figure BDA0003771819100000054
In the present invention:
Figure BDA0003771819100000054

Figure BDA0003771819100000055
Figure BDA0003771819100000055

其中,σf为关于频率的标准差,当伪随机序列为双重伪随机序列时,σf=(fi+1-fi-1)/6;当伪随机序列为混合伪随机序列时,

Figure BDA0003771819100000056
μf为关于频率的均值,μf=f',f'为滤波窗口的中点;
Figure BDA0003771819100000057
关于实部阻抗的标准差,
Figure BDA0003771819100000058
wf为加权系数,加权系数wf通过电流谐波幅值I(f)求得
Figure BDA0003771819100000059
Figure BDA00037718191000000510
为关于虚部阻抗的标准差,
Figure BDA00037718191000000511
Figure BDA00037718191000000512
Figure BDA00037718191000000513
为关于实部阻抗的均值,
Figure BDA00037718191000000514
Figure BDA00037718191000000515
为锂电池在阻抗频率f下的阻抗的实部,
Figure BDA00037718191000000516
为锂电池在阻抗频率f下的阻抗的虚部,
Figure BDA00037718191000000517
为关于虚部阻抗的均值。Among them, σ f is the standard deviation about the frequency, when the pseudo-random sequence is a double pseudo-random sequence, σ f =(f i+1 -f i-1 )/6; when the pseudo-random sequence is a mixed pseudo-random sequence,
Figure BDA0003771819100000056
μ f is the mean value about the frequency, μ f =f', f' is the midpoint of the filter window;
Figure BDA0003771819100000057
With respect to the standard deviation of the real impedance,
Figure BDA0003771819100000058
w f is the weighting coefficient, and the weighting coefficient w f is obtained by the current harmonic amplitude I(f)
Figure BDA0003771819100000059
Figure BDA00037718191000000510
is the standard deviation about the imaginary part impedance,
Figure BDA00037718191000000511
Figure BDA00037718191000000512
Figure BDA00037718191000000513
is the mean value with respect to the real impedance,
Figure BDA00037718191000000514
Figure BDA00037718191000000515
is the real part of the impedance of the lithium battery at the impedance frequency f,
Figure BDA00037718191000000516
is the imaginary part of the impedance of the lithium battery at the impedance frequency f,
Figure BDA00037718191000000517
is the mean value with respect to the imaginary part impedance.

实施例1Example 1

在本实施例中,以额定容量1.5Ah、额定电压3.6V的18650圆柱形磷酸铁锂电池为例,基于混合伪随机序列的锂电池宽频阻抗谱测试方法,得到宽频阻抗:In this embodiment, taking a 18650 cylindrical lithium iron phosphate battery with a rated capacity of 1.5Ah and a rated voltage of 3.6V as an example, the broadband impedance is obtained based on the mixed pseudo-random sequence lithium battery broadband impedance spectrum test method:

步骤1:设计一段用于宽频阻抗测量的混合伪随机序列,本实施例中第一子序列SA的长度N1为8191位,注入频率f1为7000Hz,伪随机序列类别为最大长度序列;第二序列SB的长度N2为63位,注入频率f2为53.8Hz,伪随机序列类别为最大长度序列。Step 1: Design a mixed pseudo-random sequence for broadband impedance measurement. In this embodiment, the length N1 of the first subsequence S A is 8191 bits, the injection frequency f1 is 7000 Hz, and the category of the pseudo-random sequence is the maximum length sequence; The length N 2 of the second sequence S B is 63 bits, the injection frequency f 2 is 53.8 Hz, and the category of the pseudo-random sequence is a maximum-length sequence.

步骤2:通过电池管理系统向18650圆柱形磷酸铁锂电池中注入步骤1所设计的混合伪随机序列,在信号注入过程中采集电池两端的端电压以及电流,计算锂电池在不同频率下的量测阻抗。本实施例中,信号采样频率为70kHz。Step 2: Inject the mixed pseudo-random sequence designed in step 1 into the 18650 cylindrical lithium iron phosphate battery through the battery management system, collect the terminal voltage and current at both ends of the battery during the signal injection, and calculate the lithium battery at different frequencies Measure impedance. In this embodiment, the signal sampling frequency is 70kHz.

步骤3、在频域内通过自适应双边滤波器对量测阻抗进行滤波处理,得到宽频阻抗。Step 3. Filtering the measured impedance through an adaptive bilateral filter in the frequency domain to obtain broadband impedance.

如图2为环境温度为25℃、荷电状态为80%下的宽频阻抗量测图,可以看出测量阻抗能够准确地和真实值拟合,说明实施例1所采用的方法在环境温度为25℃、荷电状态为80%下具有较高的阻抗测量精度;如图3为环境温度为25℃、荷电状态为20%下的宽频阻抗量测图,可以看出测量阻抗能够准确地和真实值拟合,说明实施例1所采用的方法在环境温度为25℃、荷电状态为20%下具有较高的阻抗测量精度。因此,实施例1所采用的方法适用于不同荷电状态的电池阻抗测量。如图4为环境温度为35℃、荷电状态为50%下的宽频阻抗量测图,可以看出测量阻抗能够准确地和真实值拟合,说明实施例1所采用的方法在环境温度为35℃、荷电状态为50%下具有较高的阻抗测量精度;如图5为环境温度为15℃,荷电状态为50%下的宽频阻抗量测图,可以看出测量阻抗能够准确地和真实值拟合,说明实施例1所采用的方法在环境温度为15℃、荷电状态为50%下具有较高的阻抗测量精度;由图4和5的结合可以看出:实施例1所采用的方法适用于不同环境温度的电池阻抗测量。As shown in Figure 2, the ambient temperature is 25°C and the state of charge is 80% of the broadband impedance measurement diagram, it can be seen that the measured impedance can be accurately fitted with the real value, indicating that the method used in Example 1 is at an ambient temperature of It has a high impedance measurement accuracy at 25°C and a state of charge of 80%. Figure 3 is a broadband impedance measurement diagram at an ambient temperature of 25°C and a state of charge of 20%. It can be seen that the measured impedance can be accurately measured. Fitting with the real value shows that the method adopted in Example 1 has a relatively high impedance measurement accuracy at an ambient temperature of 25° C. and a state of charge of 20%. Therefore, the method adopted in Embodiment 1 is applicable to the measurement of battery impedance in different states of charge. As shown in Figure 4, the ambient temperature is 35 ° C, the state of charge is 50% of the broadband impedance measurement diagram, it can be seen that the measured impedance can be accurately fitted with the real value, indicating that the method used in Example 1 is at an ambient temperature of It has high impedance measurement accuracy at 35°C and 50% state of charge; as shown in Figure 5, the broadband impedance measurement diagram at an ambient temperature of 15°C and a state of charge of 50% shows that the measured impedance can be accurately measured Fitting with the real value shows that the method adopted in embodiment 1 has higher impedance measurement accuracy at ambient temperature of 15°C and state of charge of 50%; as can be seen from the combination of Fig. 4 and 5: embodiment 1 The adopted method is suitable for battery impedance measurements at different ambient temperatures.

实施例2Example 2

在本实施例中,以额定容量1.5Ah、额定电压3.6V的18650圆柱形磷酸铁锂电池为例,基于双重伪随机序列的锂电池宽频阻抗谱测试方法,得到宽频阻抗:In this embodiment, taking a 18650 cylindrical lithium iron phosphate battery with a rated capacity of 1.5Ah and a rated voltage of 3.6V as an example, the broadband impedance is obtained based on the double pseudo-random sequence lithium battery broadband impedance spectrum test method:

步骤1:设计一段用于宽频阻抗测量的双重伪随机序列,本实施例中第三子序列的长度为8191位,注入频率f3为7000Hz,第四子序列的长度为1023位,注入频率f4为4667Hz。Step 1: Design a double pseudo-random sequence for broadband impedance measurement. In this embodiment, the length of the third subsequence is 8191 bits, the injection frequency f3 is 7000 Hz, the length of the fourth subsequence is 1023 bits, and the injection frequency f 4 is 4667Hz.

步骤2:通过电池管理系统向18650圆柱形磷酸铁锂电池中注入步骤1所设计的双重伪随机序列,在信号注入过程中采集电池两端的端电压以及电流,计算锂电池在不同频率下的量测阻抗。本实施例中,信号采样频率为70kHz。Step 2: Inject the double pseudo-random sequence designed in step 1 into the 18650 cylindrical lithium iron phosphate battery through the battery management system, collect the terminal voltage and current at both ends of the battery during the signal injection, and calculate the lithium battery at different frequencies Measure impedance. In this embodiment, the signal sampling frequency is 70kHz.

步骤3、在频域内通过自适应双边滤波器对量测阻抗进行滤波处理,得到宽频阻抗。Step 3. Filtering the measured impedance through an adaptive bilateral filter in the frequency domain to obtain broadband impedance.

如图6为环境温度为25℃、荷电状态为80%下的宽频阻抗量测图,可以看出测量阻抗能够准确地和真实值拟合,说明实施例2所采用的方法在环境温度为25℃、荷电状态为80%下具有较高的阻抗测量精度;如图7为环境温度为25℃、荷电状态为20%下的宽频阻抗量测图,可以看出测量阻抗能够准确地和真实值拟合,说明实施例2所采用的方法在环境温度为25℃、荷电状态为20%下具有较高的阻抗测量精度。将图6和图7结合可以看出:实施例2所采用的方法适用于不同荷电状态的电池阻抗测量。如图8为环境温度为35℃、荷电状态为50%下的宽频阻抗量测图,可以看出测量阻抗能够准确地和真实值拟合,说明实施例2所采用的方法在环境温度为35℃、荷电状态为50%下具有较高的阻抗测量精度;如图9为环境温度为15℃,荷电状态为50%下的宽频阻抗量测图,可以看出测量阻抗能够准确地和真实值拟合,说明实施例2所采用的方法在环境温度为15℃、荷电状态为50%下具有较高的阻抗测量精度。将图8和图9结合,说明实施例2所采用的方法适用于不同环境温度的电池阻抗测量。As shown in Figure 6, the ambient temperature is 25°C and the state of charge is 80% of the broadband impedance measurement diagram, it can be seen that the measured impedance can be accurately fitted with the real value, indicating that the method used in Example 2 is used when the ambient temperature is It has high impedance measurement accuracy at 25°C and a state of charge of 80%. Figure 7 is a broadband impedance measurement diagram at an ambient temperature of 25°C and a state of charge of 20%. It can be seen that the measured impedance can be accurately measured. Fitting with the real value shows that the method adopted in Example 2 has a relatively high impedance measurement accuracy at an ambient temperature of 25° C. and a state of charge of 20%. Combining Fig. 6 and Fig. 7, it can be seen that the method adopted in Embodiment 2 is applicable to the measurement of battery impedance in different states of charge. As shown in Figure 8, the ambient temperature is 35°C and the state of charge is 50% of the broadband impedance measurement diagram, it can be seen that the measured impedance can be accurately fitted with the real value, indicating that the method used in Example 2 is at an ambient temperature of It has high impedance measurement accuracy at 35°C and 50% state of charge; as shown in Figure 9, it is a broadband impedance measurement diagram at an ambient temperature of 15°C and a state of charge of 50%. It can be seen that the measured impedance can be accurately measured Fitting with the real value shows that the method adopted in Example 2 has a relatively high impedance measurement accuracy at an ambient temperature of 15° C. and a state of charge of 50%. Combining FIG. 8 and FIG. 9 , it is illustrated that the method adopted in Embodiment 2 is applicable to battery impedance measurement at different ambient temperatures.

以上仅是本发明的优选实施方式,本发明的保护范围并不仅局限于上述实施例,凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理前提下的若干改进和润饰,应视为本发明的保护范围。The above are only preferred implementations of the present invention, and the protection scope of the present invention is not limited to the above-mentioned embodiments, and all technical solutions under the idea of the present invention belong to the protection scope of the present invention. It should be pointed out that for those skilled in the art, some improvements and modifications without departing from the principle of the present invention should be regarded as the protection scope of the present invention.

Claims (10)

1. A method for testing a broadband impedance spectrum of a lithium battery based on a pseudorandom sequence is characterized by comprising the following steps:
step 1, designing a section of pseudo-random sequence for broadband impedance measurement;
step 2, injecting a pseudorandom sequence into the lithium battery through a battery management system, collecting the terminal voltage and current of the lithium battery in the injection process, and calculating the measurement impedance of the lithium battery under different frequencies;
and 3, filtering the measured impedance through a self-adaptive bilateral filter in a frequency domain to obtain the broadband impedance.
2. The method for testing the broadband impedance spectrum of the lithium battery based on the pseudorandom sequence as claimed in claim 1, wherein the pseudorandom sequence is a dual pseudorandom sequence or a mixed pseudorandom sequence.
3. The method as claimed in claim 2, wherein the pseudo-random sequence S is a mixed pseudo-random sequence SHFrom a first subsequence SAAnd a second subsequence SBAnd (3) stacking: sH=0.5(SA+SB) Wherein the first subsequence SAHas a length of N1Injection frequency of f1A second subsequence SBIs of length N2Injection frequency of f2And satisfies the following conditions:
Figure FDA0003771819090000011
4. the method as claimed in claim 2, wherein the dual pseudo-random sequence comprises a third subsequence and a fourth subsequence, and the injection frequency of the third subsequence is f3And injection frequency f of the fourth subsequence4Satisfies the following conditions: f. of3=1.5f4
5. The method as claimed in claim 1, wherein the measured impedance of the lithium battery at different frequencies is measured by the pseudo-random sequence-based method for testing the broadband impedance spectrum of the lithium battery
Figure FDA0003771819090000012
Comprises the following steps:
Figure FDA0003771819090000013
where V (f) is the voltage harmonic amplitude at the impedance frequency f and I (f) is the current harmonic amplitude at the impedance frequency f.
6. The method as claimed in claim 1, wherein the method for testing the broadband impedance spectrum of the lithium battery based on the pseudorandom sequence comprises
Figure FDA0003771819090000014
Comprises the following steps:
Figure FDA0003771819090000015
wherein f is the impedance frequency, fi-1To the lower bound of the filter window, fi+1In order to be the upper bound of the filtering window,
Figure FDA0003771819090000016
is the impedance, g, of the lithium battery at an impedance frequency, f1(f) Is a kernel function with respect to the frequency of the impedance,
Figure FDA0003771819090000021
is a kernel function with respect to the measured impedance.
7. The method as claimed in claim 6, wherein the pseudo-random sequence-based test method for broadband impedance spectrum of lithium battery is characterized in thatKernel function g of impedance frequency1(f) Kernel function for measuring impedance
Figure FDA0003771819090000022
Respectively as follows:
Figure FDA0003771819090000023
Figure FDA0003771819090000024
wherein σfAs standard deviation with respect to frequency, μfIs a mean value with respect to frequency, muf= f ', f' is the midpoint of the filtering window;
Figure FDA0003771819090000025
with respect to the standard deviation of the real part impedance,
Figure FDA0003771819090000026
wfis a weighting coefficient;
Figure FDA0003771819090000027
as regards the standard deviation of the imaginary impedance,
Figure FDA0003771819090000028
Figure FDA0003771819090000029
as an average value with respect to the real part impedance,
Figure FDA00037718190900000210
Figure FDA00037718190900000211
is the real part of the impedance of the lithium battery at the impedance frequency f,
Figure FDA00037718190900000212
being the imaginary part of the impedance of the lithium battery at the impedance frequency f,
Figure FDA00037718190900000213
as a mean value with respect to the imaginary impedance.
8. The method as claimed in claim 7, wherein the weighting factor w is a function of the impedance spectrum of the lithium batteryfThe current harmonic amplitude I (f) is used to obtain:
Figure FDA00037718190900000214
9. the method as claimed in claim 7, wherein when the pseudo-random sequence is a dual pseudo-random sequence, σ is determined by the pseudo-random sequencef=(fi+1-fi-1)/6。
10. The method as claimed in claim 7, wherein, when the pseudo-random sequence is a mixed pseudo-random sequence,
Figure FDA00037718190900000215
CN202210904659.4A 2022-07-29 2022-07-29 A lithium battery broadband impedance spectrum testing method based on pseudo-random sequence Active CN115267333B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210904659.4A CN115267333B (en) 2022-07-29 2022-07-29 A lithium battery broadband impedance spectrum testing method based on pseudo-random sequence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210904659.4A CN115267333B (en) 2022-07-29 2022-07-29 A lithium battery broadband impedance spectrum testing method based on pseudo-random sequence

Publications (2)

Publication Number Publication Date
CN115267333A true CN115267333A (en) 2022-11-01
CN115267333B CN115267333B (en) 2025-04-25

Family

ID=83771169

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210904659.4A Active CN115267333B (en) 2022-07-29 2022-07-29 A lithium battery broadband impedance spectrum testing method based on pseudo-random sequence

Country Status (1)

Country Link
CN (1) CN115267333B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2385477Y (en) * 1999-07-30 2000-06-28 河南省华康通信研究所 All-digitizing high voltage transmission line digital carrier device
US6975689B1 (en) * 2000-03-30 2005-12-13 Mcdonald James Douglas Digital modulation signal receiver with adaptive channel equalization employing discrete fourier transforms
US20080107352A1 (en) * 2006-11-03 2008-05-08 Siemens Corporate Research, Inc. System and Method for Structure Enhancement and Noise Reduction in Medical Images
CN107067367A (en) * 2016-09-08 2017-08-18 南京工程学院 A kind of Image Super-resolution Reconstruction processing method
CN107478993A (en) * 2017-07-11 2017-12-15 中国科学技术大学 A kind of bilateral insulation resistance monitoring device of electric automobile power battery
CN108510459A (en) * 2018-04-08 2018-09-07 哈尔滨理工大学 One kind is based on wavelet adaptive threshold and bilateral filtering image noise reduction algorithm
CN109035156A (en) * 2018-06-15 2018-12-18 浙江医院 Medicine CT image denoising method based on DNST
CN109497988A (en) * 2018-12-13 2019-03-22 复旦大学 Multi-modal electro-physiological signals monitor system for domestic sleeping monitoring field
CN114720890A (en) * 2022-03-23 2022-07-08 四川大学 Lithium battery internal temperature estimation method based on broadband impedance spectrum

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2385477Y (en) * 1999-07-30 2000-06-28 河南省华康通信研究所 All-digitizing high voltage transmission line digital carrier device
US6975689B1 (en) * 2000-03-30 2005-12-13 Mcdonald James Douglas Digital modulation signal receiver with adaptive channel equalization employing discrete fourier transforms
US20080107352A1 (en) * 2006-11-03 2008-05-08 Siemens Corporate Research, Inc. System and Method for Structure Enhancement and Noise Reduction in Medical Images
CN107067367A (en) * 2016-09-08 2017-08-18 南京工程学院 A kind of Image Super-resolution Reconstruction processing method
CN107478993A (en) * 2017-07-11 2017-12-15 中国科学技术大学 A kind of bilateral insulation resistance monitoring device of electric automobile power battery
CN108510459A (en) * 2018-04-08 2018-09-07 哈尔滨理工大学 One kind is based on wavelet adaptive threshold and bilateral filtering image noise reduction algorithm
CN109035156A (en) * 2018-06-15 2018-12-18 浙江医院 Medicine CT image denoising method based on DNST
CN109497988A (en) * 2018-12-13 2019-03-22 复旦大学 Multi-modal electro-physiological signals monitor system for domestic sleeping monitoring field
CN114720890A (en) * 2022-03-23 2022-07-08 四川大学 Lithium battery internal temperature estimation method based on broadband impedance spectrum

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李伟恒;黄秋安;杨维明;杨昌平;张久俊;: "基于伪随机二进制序列的阻抗谱快速重构及其在电化学能源领域的应用", 电化学, no. 03, 28 June 2020 (2020-06-28), pages 57 - 75 *

Also Published As

Publication number Publication date
CN115267333B (en) 2025-04-25

Similar Documents

Publication Publication Date Title
CN111736085B (en) Lithium ion battery health state estimation method based on electrochemical impedance spectrum
US10663525B2 (en) Method and device for monitoring state of charge and state of health of lithium-ion battery
CN108845264B (en) Wavelet-based battery health state estimation method
CN113659245B (en) An electrochemical device heating method, electrochemical device and electrical equipment
CN111458642A (en) A kind of non-destructive testing method for lithium evolution of lithium ion battery
Du et al. Hybrid pseudorandom sequence for broadband impedance measurements of lithium-ion batteries
CN112946489B (en) A Fast Capacity Assessment Method Based on Low Frequency EIS
CN107656215A (en) A kind of battery functi on method for diagnosing status based on constant current mode impedance spectrum
Song et al. Anomalous diffusion models in frequency-domain characterization of lithium-ion capacitors
CN114720890B (en) Lithium battery internal temperature estimation method based on broadband impedance spectrum
CN111123136A (en) On-line measuring device and on-line measuring method of lithium-ion battery state of health
WO2016106501A1 (en) Equivalent circuit model of battery
CN117347886B (en) Lithium ion battery impedance analysis lithium separation detection method based on DEIS technology
CN109991556A (en) A diagnostic method for short-term failure of lithium iron phosphate power battery
CN110275114A (en) Accumulator internal resistance on-line monitoring method based on combined filter algorithm
CN114675195A (en) On-line monitoring system and method for lithium battery core temperature
Shu et al. Research on the Electrochemical Impedance Spectroscopy Evolution of Sodium-Ion Batteries in Different States
Geng et al. Short binary sequence-based fast electrochemical impedance spectroscopy measurement for lithium-ion batteries
CN118625149A (en) An equivalent circuit model of a lithium-ion battery and an online parameter identification method thereof
CN202330560U (en) Storage battery impedance detection device excited by SPWM (Sinusoidal Pulse Width Modulation) signals generated by single-chip microcomputer
CN115480181A (en) A Method for Analyzing Battery Health State Based on Reversible and Irreversible Heat
CN115267333A (en) Lithium battery broadband impedance spectrum testing method based on pseudorandom sequence
Erol Equivalent Circuit Model for Electrochemical Impedance Spectroscopy of Commercial 18650 Lithium‐Ion Cell Under Over‐Discharge and Overcharge Conditions
CN114397588A (en) Lithium ion battery capacity prediction method utilizing echelon
CN210863986U (en) A battery internal resistance measurement system

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