CN117909648A - 一种适用多谱线插值补零的谐波测量方法和系统 - Google Patents
一种适用多谱线插值补零的谐波测量方法和系统 Download PDFInfo
- Publication number
- CN117909648A CN117909648A CN202311692858.4A CN202311692858A CN117909648A CN 117909648 A CN117909648 A CN 117909648A CN 202311692858 A CN202311692858 A CN 202311692858A CN 117909648 A CN117909648 A CN 117909648A
- Authority
- CN
- China
- Prior art keywords
- data
- zero padding
- new
- frequency
- interpolation
- 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.)
- Pending
Links
- 238000000691 measurement method Methods 0.000 title claims abstract description 33
- 238000004364 calculation method Methods 0.000 claims abstract description 35
- 238000005070 sampling Methods 0.000 claims abstract description 21
- 238000005259 measurement Methods 0.000 claims abstract description 18
- 230000003595 spectral effect Effects 0.000 claims description 43
- 238000001228 spectrum Methods 0.000 claims description 26
- 238000000034 method Methods 0.000 claims description 22
- 230000006870 function Effects 0.000 claims description 17
- 238000004891 communication Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 7
- 230000003044 adaptive effect Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 description 15
- 238000012937 correction Methods 0.000 description 8
- 230000009471 action Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000002401 inhibitory effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 206010063385 Intellectualisation Diseases 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005067 remediation Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/40—Arrangements for reducing harmonics
Landscapes
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
一种适用多谱线插值补零的谐波测量方法及系统,多谱线插值补零的谐波测量方法包括:对电力系统中电压、电流电气量进行固定长度采样并加窗构建数据窗,然后对数据窗数据进行补零构建新的数据窗,进而进行FFT计算并结合改进多谱线插值拟合计算,最后提取出所测电气量频率成分及其对应的幅值、相位。本发明适用于非整数次幂的FFT,相比常规补零FFT,具有更高的量测精度和更低的运算代价谐波测量方法,且具有很好的工程使用价值。
Description
技术领域
本发明涉及量测方法技术领域,更具体地涉及一种适用多谱线插值补零的谐波测量方法和系统。
背景技术
随着可再生能源技术的越发成熟和电网的智能化程度加快,电力系统面临高比例新能源和高比例电力电子装置的“双高”带来的一系列问题。随着非线性负载的大量接入,给电网带来越来越严重的谐波污染。对谐波参数的准确测量是评估电能质量和进行谐波治理的基础。
对于谐波测量的算法主要有:小波变换法、神经网络法、希尔伯特-黄变换法、FFT(Fast Fourier Transform,快速傅里叶变换)。FFT由于计算量小、易于嵌入式代码实现而被广泛应用于电力系统测量。FFT对数据长度有一定的要求,比如基-2FFT要求数据长度满足是2的整数次幂,基-4FFT要求数据长度满足是4的整数次幂,在国际标准IEC61000-4-7中规定在工频50Hz的电网信号采集分析中,进行一次FFT运算的时间窗必须是10周期,即200ms。不同的采样率,例如6400Hz,每周波采样点数128点,很难保证200ms数据窗是2或4的整数次幂,因此需要非整数次幂的FFT算法。一般通过混合基实现非整数次幂FFT,但是混合基运算要高于普通FFT,且要求数据窗长度为素数,使得灵活应用的场景受限。此外各类加窗FFT与谱线插值校正算法的结合可以使得测量精度普遍达到10-5以上的等级,但对于补零FFT的加窗插值拟合校正,由于补零操作时窗函数的谱线间距减小,而整体频谱又不变,导致未补零的加窗FFT插值拟合公式不再适用于补零加窗FFT。
发明内容
为解决现有技术中存在的不足,本发明提供一种适用多谱线插值补零的谐波测量的方法和系统,为适应非整数次幂数据窗进行FFT提供可靠方法,同时为抑制频谱泄露,针对补零FFT进行加窗并改进插值拟合算法,从而提供一种适合工程应用的高精度谐波测量算法,解决固定时间周期计算数据窗长度非整数次幂无法使用高精度FFT算法的问题。
本发明采用如下的技术方案。
本发明的第一方面提供一种适用多谱线插值补零的谐波测量方法,包括以下步骤:
步骤1,固定频率采集N点电力系统的电压量、电流量序列;
步骤2,对采集的N点电压量、电流量序列进行加窗,得到加窗后N点离散采样数据xw(n);
步骤3,对N点加窗数据补零成N1点数据序列,构成新的数据窗数据xw_new(n),N1为大于N的2的整数次幂整数集合中最小的整数;
步骤4,对步骤3中的新的数据窗数据进行基2FFT计算和三谱线插值拟合计算;
步骤5,提取满足门槛值电气量的频率成分及其幅值、相位。
优选地,步骤1中,对电力系统的电压量、电流量进行离散采样得到x(n)并缓存N点,n=0、1、2…、N-1,采样频率fs,频率分辨率为fs/N,单个信号以如下公式所示:
式中:x(n)为电气量离散采样数据,A0为直流分量,A1为正弦信号的峰值,f1为正弦信号频率,fs为采样频率,为正弦信号初始相位。
优选地,步骤2中,加窗后N点离散采样数据xw(n)以如下公式所示:
xw(n)=x(n)·w5-MSLD(n)
式中:w5-MSLD(n)为5项最快旁瓣衰减速率窗,xw(n)为对离散采样数据进行加窗后的数据,即数据时域的乘积后的结果。
优选地,步骤3中,新的数据窗数据xw_new(n)以如下公式所示:
式中:xw_new(n)为补零后新的数据窗数据,N为原始数据长度。
更进一步地,步骤3中,新的数据窗xw_new(n)的长度N1计算以如下公式所示:
N1=pow2(2,nextpow2(N))
式中:N为原始数据窗长度,nextpow2(N)为取不小于N的最小的2的幂,pow2()为2的幂次方计算。
优选地,步骤4中,采用基2FFT计算后的频谱Xw_new(k)以如下公式所示:
式中:Ak为第k次谐波的幅值,e为自然对数的底数,为第k次谐波的相位,W5-MSLD_new(n)为补零至N1点的混合卷积w5-MSLD_new窗的频谱,k0为真实频谱的谱线位置,Xw_new(k)为基于FFT计算后的频谱。
更进一步地,步骤4中,得到的频谱Xw_new(k)进行改进的三谱线插值拟合,选用Xw_new(k)中极大值附近的三根信息量大的谱线k1、k2、k3,其中谱线k3的信息量最大,谱线k1的信息量在k1、k2、k3最小;
对应的谱线幅值分别X1、X2、X3,X1=|X(k1)|、X2=|X(k2)|、X3=|X(k3)|,幅值大小关系是X1<X2>X3,X1和X3的大小关系取决于真实信号频率处于K2的左侧还是右侧,若真实频率对应谱线k0在k2的左侧,则X1>X3,反之X1<X3;
令α=k0-k2,将X1、X2、X3用于插值拟合后的频谱以如下公式所示:
α=f(β)
式中:β为极大值附近三根谱线组合比值关系,α为真实频率对应谱线k0与极大值谱线k2之间的距离,W5_MSLD_new(1-α)为谱线k3对应的幅值,W5_MSLD_new(-α)为谱线k2对应的幅值,W5_MSLD_new(-1-α)为谱线k1对应的幅值,f(β)为多项式拟合的公式,abs()表示对复数求取模值,。
优选地,步骤5中,从插值拟合后的频谱中提取谐波的幅值、频率和相位以如下公式所示:
式中:k为真实频率所在的谱线位置,Δf为补零后等效的频率分辨率,f为所计算谐波频率,A为所计算谐波幅值,为所计算谐波相位。
更进一步地,步骤5中,当N1大于1000时,函数A以如下公式所示:
A=N1 -1(X1+2X2+X3)g(α)
式中:g(α)为窗函数固定后谱线位置与幅值之间拟合函数,N1 -1为补零后数据窗长度的倒数。
本发明的第二方面提供一种适用多谱线插值补零的谐波测量系统,基于所述的一种适用多谱线插值补零的谐波测量方法,包括:
电力系统电压电流量获取模块,序列加窗模块,加窗补零模块,数据窗计算模块,电压电流频率和幅值相位提取模块;
电力系统电压电流量获取模块用于离散采样获取电力系统中N点电压量、电流量序列;序列加窗模块用于对获取的序列进行加窗处理;加窗补零模块用于加窗数据补零构成新的窗数据;数据窗计算模块用于基2FFT计算和改进多谱线插值拟合计算新的窗数据;电压电流频率和幅值相位提取模块用于提取计算后的窗数据提取电力系统电压量、电流量的频率、幅值、相位。
本发明的第三方面提供一种基于适用多谱线插值补零的谐波测量方法的电子设备,运行所述的一种适用多谱线插值补零的谐波测量方法:
包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序,实现所述一种适用多谱线插值补零的谐波测量方法的步骤。
本发明的第四方面提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现所述的一种适用多谱线插值补零的谐波测量方法。
本发明的有益效果在于,与现有技术相比,本发明提供一种适用多谱线插值补零的谐波测量的方法和系统,为适应非整数次幂数据窗进行FFT提供可靠方法,同时为抑制频谱泄露,针对补零FFT进行加窗并改进插值拟合算法,从而提供一种适合工程应用的高精度谐波测量算法,解决固定时间周期计算数据窗长度非整数次幂无法使用高精度FFT算法的问题;另外本发明适用于非整数次幂的FFT,相比常规补零FFT,具有更高的量测精度和更低的运算代价谐波测量方法,且具有很好的工程使用价值;对于窗函数补零后的拟合校正,本发明给出的公式适用于任何窗函数;本发明可以应用于非整数次幂的FFT算法,不仅计算速度高于一般的FFT计算,还可以灵活应用于不同场景;本发明可以满足对于基-2FFT或者基-4FFT的计算要求,数据长度满足是2的整数次幂,基-4FFT要求数据长度满足是4的整数次幂,在保证计算准确度的情况下,能够提高计算速度,大大减少计算时间。
附图说明
图1是改进多谱线插值补零的电气量特征提取方法流程图;
图2是本发明的一种适用多谱线插值补零的谐波测量系统示意图;
图3是3种窗的频谱响应对比图;
图4是补零前后5-MSLD窗时频域图;
图5基波频率固定时幅值误差;
图6基波频率固定时相位误差;
图7基波频率固定时频率误差。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明的技术方案进行清楚、完整地描述。本申请所描述的实施例仅仅是本发明一部分的实施例,而不是全部实施例。基于本发明精神,本领域普通技术人员在没有作出创造性劳动前提下所获得的有所其它实施例,都属于本发明的保护范围。
本发明的实施例1提供一种适用多谱线插值补零的谐波测量方法,如图1所示,包括以下步骤:
步骤1,固定频率采集N点电力系统的电压量、电流量序列。
在优选但非限制性地实施方式中,步骤1中,对电力系统的电压量、电流量进行离散采样得到x(n)并缓存N点,n=0、1、2…、N-1,采样频率fs,频率分辨率为fs/N,单个信号以如下公式所示:
式中:x(n)为电气量离散采样数据,A0为直流分量,A1为正弦信号的峰值,f1为正弦信号频率,fs为采样频率,为正弦信号初始相位。
步骤2,对采集的N点电压量、电流量序列进行加窗,得到加窗后N点离散采样数据xw(n)。
在优选但非限制性地实施方式中,步骤2中,加窗后N点离散采样数据xw(n)以如下公式所示:
xw(n)=x(n)·w5-MSLD(n)
式中:w5-MSLD(n)为5项最快旁瓣衰减速率窗,xw(n)为对离散采样数据进行加窗后的数据,即数据时域的乘积后的结果。
步骤3,对N点加窗数据补零成N1点数据序列,构成新的数据窗数据xw_new(n),N1为大于N的2的整数次幂整数集合中最小的整数。
在优选但非限制性地实施方式中,步骤3中,新的数据窗数据xw_new(n)以如下公式所示:
式中:xw_new(n)为补零后新的数据窗数据,N为原始数据长度。
更进一步地,步骤3中,新的数据窗xw_new(n)的长度N1计算以如下公式所示:
N1=pow2(2,nextpow2(N))
式中:N为原始数据窗长度,nextpow2(N)为取不小于N的最小的2的幂,pow2()为2的幂次方计算。
步骤4,对步骤3中的新的数据窗数据进行基2FFT计算和三谱线插值拟合计算;
在优选但非限制性地实施方式中,步骤4中,采用基2FFT计算后的频谱Xw_new(k)以如下公式所示:
式中:Ak为第k次谐波的幅值,e为自然对数的底数,为第k次谐波的相位,W5-MSLD_new(n)为补零至N1点的混合卷积w5-MSLD_new窗的频谱,k0为真实频谱的谱线位置,Xw_new(k)为基于FFT计算后的频谱。
更进一步地,步骤4中,得到的频谱Xw_new(k)进行改进的三谱线插值拟合,选用Xw_new(k)中极大值附近的三根信息量大的谱线k1、k2、k3,其中谱线k3的信息量最大,谱线k1的信息量在k1、k2、k3最小;
对应的谱线幅值分别X1、X2、X3,X1=|X(k1)|、X2=|X(k2)|、X3=|X(k3)|,幅值大小关系是X1<X2>X3,X1和X3的大小关系取决于真实信号频率处于K2的左侧还是右侧,若真实频率对应谱线k0在k2的左侧,则X1>X3,反之X1<X3;
令α=k0-k2,将X1、X2、X3用于插值拟合后的频谱以如下公式所示:
α=f(β)
式中:β为极大值附近三根谱线组合比值关系,α为真实频率对应谱线k0与极大值谱线k2之间的距离,W5_MSLD_new(1-α)为谱线k3对应的幅值,W5_MSLD_new(-α)为谱线k2对应的幅值,W5_MSLD_new(-1-α)为谱线k1对应的幅值,f(β)为多项式拟合的公式,abs()表示对复数求取模值,。
步骤5,提取满足门槛值电气量的频率成分及其幅值、相位。
在优选但非限制性地实施方式中,步骤5中,从插值拟合后的频谱中提取谐波的幅值、频率和相位以如下公式所示:
式中:k为真实频率所在的谱线位置,Δf为补零后等效的频率分辨率,f为所计算谐波频率,A为所计算谐波幅值,为所计算谐波相位。
更进一步地,步骤5中,当N1大于1000时,函数A以如下公式所示:
A=N1 -1(X1+2X2+X3)g(α)
式中:g(α)为窗函数固定后谱线位置与幅值之间拟合函数,N1 -1为补零后数据窗长度的倒数。
本发明的实施例2提供一种适用多谱线插值补零的谐波测量系统,基于所述的一种适用多谱线插值补零的谐波测量方法,如图2所示,包括:
电力系统电压电流量获取模块,序列加窗模块,加窗补零模块,数据窗计算模块,电压电流频率和幅值相位提取模块;
电力系统电压电流量获取模块用于离散采样获取电力系统中N点电压量、电流量序列;序列加窗模块用于对获取的序列进行加窗处理;加窗补零模块用于加窗数据补零构成新的窗数据;数据窗计算模块用于基2FFT计算和改进多谱线插值拟合计算新的窗数据;电压电流频率和幅值相位提取模块用于提取计算后的窗数据提取电力系统电压量、电流量的频率、幅值、相位。
本发明的实施例3供一种基于适用多谱线插值补零的谐波测量方法的电子设备,运行所述的一种适用多谱线插值补零的谐波测量方法:
包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序,实现所述一种适用多谱线插值补零的谐波测量方法的步骤。
本发明的实施例4提供一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现所述的一种适用多谱线插值补零的谐波测量方法。
本发明的实施例5通过实施例1中提出的拟合方法分别计算出hanning窗、4项1阶莱夫-文森特窗、5项MSLD窗在2560点补零至4096点时的校正拟合系数,其中:
5项MSLD窗补零三谱线:
α=3.133716434041538β-0.304166041953860β3+0.048340202559174β5-0.008782569614980β7
g(α)=3.052506661907434+0.244228710483206α2+0.010402944634988α4+0.0003186550147148620α6
按照预先设定的门槛提取计算所得的各频率、幅值、相位。
采用仿真对比的方法和电力系统谐波信号进行分析,所施加的仿真谐波模型为:
式中:基波频率f1=50.1Hz,fs为采样频率12800Hz,Ai、为各次幅值和初相位。
谐波模型参数如下表:
表1谐波参数模型表
对上述仿真分别进行加hanning窗即汉明窗、4项1阶莱夫-文森特窗、5项MSLD窗补零FFT计算,采样2560点,并将数据补零至4096点,图3为三种窗函数的频谱响应,可知5项MSLD窗衰减速率最大,其次时4项1阶RV窗,hanning窗衰减速率最差,对于补零FFT,由于信号波形分辨率增加,谱线之间间距减小,相互影响增大,如图4所示,补零后谱线变得更加密集,主瓣内包含更多的谱线,所以衰减速率大的窗更加有利于减少频谱泄露、更加准确计算信号的幅值、频率、相位。
加不同窗的补零FFT根据多谱线插值校正公式计算后,求得谐波信号对应幅值、相位及频率如图5~7所示。
通过本发明的方法对不同窗进行补零多谱线插值校正拟合,并对电力系统谐波信号进行分析,可得到以下结果:
幅值测量结果相对误差如图5所示,5项MSLD窗的效果最好,比4项1阶RV窗的精度提高一个数量级,比Hanning窗的精度提高2~3个数量级;相位测量结果相对误差如图6所示,5项MSLD窗的效果最好,比4项1阶RV窗的精度提高0.5个数量级,比Hanning窗的精度提高3个数量级。频率测量结果相对误差如图7所示,5项MSLD窗的效果最好,比4项1阶RV窗的精度提高0.5个数量级,比Hanning窗的精度提高3个数量级。结果表明这种针对非整数幂的加窗FFT谱线插值校正拟合算法对不同窗函数都是有效的,同时验证旁瓣衰减速率大的窗更有利于提高精度。
本发明的有益效果在于,与现有技术相比,本发明提供一种适用多谱线插值补零的谐波测量的方法和系统,为适应非整数次幂数据窗进行FFT提供可靠方法,同时为抑制频谱泄露,针对补零FFT进行加窗并改进插值拟合算法,从而提供一种适合工程应用的高精度谐波测量算法,解决固定时间周期计算数据窗长度非整数次幂无法使用高精度FFT算法的问题;另外本发明适用于非整数次幂的FFT,相比常规补零FFT,具有更高的量测精度和更低的运算代价谐波测量方法,且具有很好的工程使用价值;对于窗函数补零后的拟合校正,本发明给出的公式适用于任何窗函数;本发明可以应用于非整数次幂的FFT算法,不仅计算速度高于一般的FFT计算,还可以灵活应用于不同场景;本发明可以满足对于基-2FFT或者基-4FFT的计算要求,数据长度满足是2的整数次幂,基-4FFT要求数据长度满足是4的整数次幂,在保证计算准确度的情况下,能够提高计算速度,大大减少计算时间。
在上述实施例中,可以全部或部分地通过软件、硬件、固件或者其任意组合来实现。当使用软件实现时,可以全部或部分地以计算机程序产品的形式实现。所述计算机程序产品包括一个或多个计算机指令。在计算机上加载和执行所述计算机程序指令时,全部或部分地产生按照本发明实施例所述的流程或功能。所述计算机可以是通用计算机、专用计算机、计算机网络、或者其他可编程装置。所述计算机指令可以存储在计算机可读存储介质中,或者从一个计算机可读存储介质向另一个计算机可读存储介质传输,例如,所述计算机指令可以从一个网站站点、计算机、服务器或数据中心通过有线(例如同轴电缆、光纤、数字用户线(DSL))或无线(例如红外、无线、微波等)方式向另一个网站站点、计算机、服务器或数据中心进行传输。所述计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。所述可用介质可以是磁性介质,(例如,软盘、硬盘、磁带)、光介质(例如,DVD)、或者半导体介质(例如固态硬盘Solid State Disk(SSD))等。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。另外,术语“包括”、“包含”或者任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个......”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用相关的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
最后应当说明的是,以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (12)
1.一种适用多谱线插值补零的谐波测量方法,其特征在于,包括以下步骤:
步骤1,固定频率采集N点电力系统的电压量、电流量序列;
步骤2,对采集的N点电压量、电流量序列进行加窗,得到加窗后N点离散采样数据xw(n);
步骤3,对N点加窗数据补零成N1点数据序列,构成新的数据窗数据xw_new(n),N1为大于N的2的整数次幂整数集合中最小的整数;
步骤4,对步骤3中的新的数据窗数据进行基2FFT计算和三谱线插值拟合计算;
步骤5,提取满足门槛值电气量的频率成分及其幅值、相位。
2.根据权利要求1所述一种适用多谱线插值补零的谐波测量方法,其特征在于:
步骤1中,对电力系统的电压量、电流量进行离散采样得到x(n)并缓存N点,n=0、1、2…、N-1,采样频率fs,频率分辨率为fs/N,单个信号以如下公式所示:
式中:x(n)为电气量离散采样数据,A0为直流分量,A1为正弦信号的峰值,f1为正弦信号频率,fs为采样频率,为正弦信号初始相位。
3.根据权利要求1所述一种适用多谱线插值补零的谐波测量方法,其特征在于:
步骤2中,加窗后N点离散采样数据xw(n)以如下公式所示:
xw(n)=x(n)·w5-MSLD(n)
式中:w5-MSLD(n)为5项最快旁瓣衰减速率窗,xw(n)为对离散采样数据进行加窗后的数据,即数据时域的乘积后的结果。
4.根据权利要求1所述一种适用多谱线插值补零的谐波测量方法,其特征在于:
步骤3中,新的数据窗数据xw_new(n)以如下公式所示:
式中:xw_new(n)为补零后新的数据窗数据,N为原始数据长度。
5.根据权利要求4所述一种适用多谱线插值补零的谐波测量方法,其特征在于:
步骤3中,新的数据窗xw_new(n)的长度N1计算以如下公式所示:
N1=pow2(2,nextpow2(N))
式中:N为原始数据窗长度,nextpow2(N)为取不小于N的最小的2的幂,pow2()为2的幂次方计算。
6.根据权利要求1所述一种适用多谱线插值补零的谐波测量方法,其特征在于:
步骤4中,采用基2FFT计算后的频谱Xw_new(k)以如下公式所示:
式中:Ak为第k次谐波的幅值,e为自然对数的底数,为第k次谐波的相位,W5-MSLD_new(n)为补零至N1点的混合卷积w5-MSLD_new窗的频谱,k0为真实频谱的谱线位置,Xw_new(k)为基于FFT计算后的频谱。
7.根据权利要求5所述一种适用多谱线插值补零的谐波测量方法,其特征在于:
步骤4中,得到的频谱Xw_new(k)进行改进的三谱线插值拟合,选用Xw_new(k)中极大值附近的三根信息量大的谱线k1、k2、k3,其中谱线k3的信息量最大,谱线k1的信息量在k1、k2、k3最小;
对应的谱线幅值分别X1、X2、X3,X1=|X(k1)|、X2=|X(k2)|、X3=|X(k3)|,幅值大小关系是X1<X2>X3,X1和X3的大小关系取决于真实信号频率处于K2的左侧还是右侧,若真实频率对应谱线k0在k2的左侧,则X1>X3,反之X1<X3;
令α=k0-k2,将X1、X2、X3用于插值拟合后的频谱以如下公式所示:
α=f(β)
式中:β为极大值附近三根谱线组合比值关系,α为真实频率对应谱线k0与极大值谱线k2之间的距离,W5_MSLD_new(1-α)为谱线k3对应的幅值,W5_MSLD_new(-α)为谱线k2对应的幅值,W5_MSLD_new(-1-α)为谱线k1对应的幅值,f(β)为多项式拟合的公式,abs()表示对复数求取模值,。
8.根据权利要求1所述一种适用多谱线插值补零的谐波测量方法,其特征在于:
步骤5中,从插值拟合后的频谱中提取谐波的幅值、频率和相位以如下公式所示:
式中:k为真实频率所在的谱线位置,Δf为补零后等效的频率分辨率,f为所计算谐波频率,A为所计算谐波幅值,为所计算谐波相位。
9.根据权利要求7所述一种适用多谱线插值补零的谐波测量方法,其特征在于:
当N1大于1000时,函数A以如下公式所示:
A=N1 -1(X1+2X2+X3)g(α)
式中:g(α)为窗函数固定后谱线位置与幅值之间拟合函数,N1 -1为补零后数据窗长度的倒数。
10.根据权利要求1-8任一项权利要求所述方法的一种适用多谱线插值补零的谐波测量系统,包括:电力系统电压电流量获取模块,序列加窗模块,加窗补零模块,数据窗计算模块,电压电流频率和幅值相位提取模块;其特征在于:
电力系统电压电流量获取模块用于离散采样获取电力系统中N点电压量、电流量序列;序列加窗模块用于对获取的序列进行加窗处理;加窗补零模块用于加窗数据补零构成新的窗数据;数据窗计算模块用于基2FFT计算和改进多谱线插值拟合计算新的窗数据;电压电流频率和幅值相位提取模块用于提取计算后的窗数据提取电力系统电压量、电流量的频率、幅值、相位。
11.一种基于适用多谱线插值补零的谐波测量方法的电子设备,运行根据权利要求1-8任意一项权利要求所述的一种适用多谱线插值补零的谐波测量方法,其特征在于:
包括处理器、通信接口、存储器和通信总线,其中,处理器,通信接口,存储器通过通信总线完成相互间的通信;
存储器,用于存放计算机程序;
处理器,用于执行存储器上所存放的程序,实现权利要求1-8任一项权利要求所述的一种适用多谱线插值补零的谐波测量方法的步骤。
12.一种计算机可读存储介质,其上存储有计算机程序,其特征在于:
该程序被处理器执行时实现权利要求1-8任一所述的一种适用多谱线插值补零的谐波测量方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311692858.4A CN117909648A (zh) | 2023-12-11 | 2023-12-11 | 一种适用多谱线插值补零的谐波测量方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311692858.4A CN117909648A (zh) | 2023-12-11 | 2023-12-11 | 一种适用多谱线插值补零的谐波测量方法和系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN117909648A true CN117909648A (zh) | 2024-04-19 |
Family
ID=90684327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311692858.4A Pending CN117909648A (zh) | 2023-12-11 | 2023-12-11 | 一种适用多谱线插值补零的谐波测量方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117909648A (zh) |
-
2023
- 2023-12-11 CN CN202311692858.4A patent/CN117909648A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111984920B (zh) | 次/超同步谐波参数识别方法、装置、设备和介质 | |
CN106018956B (zh) | 一种加窗谱线插值的电力系统频率计算方法 | |
CN111222088B (zh) | 一种改进的平顶自卷积窗加权电力谐波幅值估计方法 | |
CN114977216A (zh) | 振荡信号的参数辨识方法及终端 | |
CN111537776A (zh) | 一种行波波头标定方法、装置、终端以及介质 | |
CN111999635A (zh) | 一种基于4项5阶Nuttall窗的板卡故障信号分析方法及终端 | |
CN110287853B (zh) | 一种基于小波分解的暂态信号去噪方法 | |
CN117909648A (zh) | 一种适用多谱线插值补零的谐波测量方法和系统 | |
CN111523231A (zh) | 一种基于EEMD和Prony方法的次同步振荡分析方法 | |
CN109736793B (zh) | 一种评价储层孔隙结构的方法及设备 | |
Matusiak et al. | Noniterative method for frequency estimation based on interpolated DFT with low-order harmonics elimination | |
Rodrigues et al. | Low-cost embedded measurement system for power quality frequency monitoring | |
CN116150594A (zh) | 一种频谱测试数据中开关类要素特征的识别方法 | |
CN116046968A (zh) | 一种液相色谱工作站数据处理方法、系统及可存储介质 | |
CN113341220B (zh) | 含噪多频衰减实信号频率估计方法 | |
CN115494303A (zh) | 一种emi接收机信号转化方法、装置及存储介质 | |
Nunzi et al. | A procedure for highly reproducible measurements of ADC spectral parameters | |
CN114184838A (zh) | 基于sn互卷积窗的电力系统谐波检测方法、系统及介质 | |
AU2014202396A1 (en) | Adaptive frequency-domain windowing | |
Wolf et al. | Amplitude and frequency estimator for aperiodic multi-frequency noisy vibration signals of a tram gearbox | |
CN117031382B (zh) | 一种基于信噪比和正弦拟合方法的单端行波装置校验方法 | |
CN109948111A (zh) | 一种基于离散数学的emi接收机检波器模型 | |
CN113191317B (zh) | 一种基于极点构造低通滤波器的信号包络提取方法和装置 | |
CN117493769A (zh) | 一种加窗fft算法的噪声畸变信号估计方法及装置 | |
CN114325165A (zh) | 一种电力系统故障信号识别方法 |
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 |