CN106814246A - 基于Kaiser窗FFT单峰插值的电能计量方法 - Google Patents
基于Kaiser窗FFT单峰插值的电能计量方法 Download PDFInfo
- Publication number
- CN106814246A CN106814246A CN201611157854.6A CN201611157854A CN106814246A CN 106814246 A CN106814246 A CN 106814246A CN 201611157854 A CN201611157854 A CN 201611157854A CN 106814246 A CN106814246 A CN 106814246A
- Authority
- CN
- China
- Prior art keywords
- frequency
- amplitude
- unimodal
- formula
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R22/00—Arrangements for measuring time integral of electric power or current, e.g. electricity meters
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
Abstract
本发明公开了一种基于Kaiser窗FFT单峰插值的电能计量方法,包括:基于Kaiser窗FFT单峰插值修正公式、FFT的频移性和窗函数Kaiser的频谱在基波、谐波和间谐波峰值频率点左、右两侧最大谱线和次最大谱线的幅值,利用多项式拟合公式,得到信号基波、谐波和间谐波的通用幅值、相位和频率插值修正公式;根据通用幅值、相位和频率插值修正公式计算被测电压信号、电流信号的基波、谐波和间谐波的幅值、相位和频率参量;进而计算基波电能、谐波电能、间谐波电能和总电能。本发明具有如下优点:单峰插值修正算法,利用Kaiser窗函数形状参数的灵活性,实现动态信号的高精度谐波参量的插值修正,具有较高的电能计量计算准确度。
Description
技术领域
本发明涉及供电系统电能计量领域以及数字化电能计量领域,具体涉及一种基于Kaiser窗FFT单峰插值的电能计量方法。
背景技术
随着越来越多的电气化铁道、轧钢机、电弧炉、电动汽车充电桩、提升机、智能家电、电梯等非线性动态电力负荷的投入使用,在供电电网中产生了大量电力谐波和间谐波。大量电力谐波、间谐波、直流衰减、电压跌落等干扰成分,造成了供电电网电流、电压波形的畸变。这些干扰的存在,不仅严重影响供电电能质量,危及供电系统的安全稳定运行,甚至导致严重的电力事故,造成巨大的经济损失,而且影响电能计量的准确性。
电能是电网公司、售电公司等进行经济核算的依据,其能否被正确、准确地计量,直接关系到电力供需双方的经济效益。传统的电能计量方式,是直接按电能定义,将电压、电流信号的瞬时采样值相乘、再进行时间积分来计算电能的所谓全电能计量方式。而风力发电机和光伏电站的并网发电,以及电弧炉、轧钢机、电力机车、充电桩等非线性动态电力负荷的大量投运,在供电电网中,不仅作为负荷要消耗功率,而且作为谐波源、污染源,造成供电电流、电压波形畸变的同时,会向供电电网注入大量谐波及间谐波功率等。若按传统的全电能计量方式,谐波源用户虽然会向电网注入大量谐波干扰,却还少交电费;而广大的普通用电户不仅受到谐波的干扰和危害,而且还要多交电费;同时,这种全电能计量方式也没有考虑因非同步采样而导致的频谱泄露和栅栏效应,故采用这种电能计量方式,无疑可能带来较大误差。
而另一种所谓基波电能计量方式,虽然可解决普通电力用户的用电公平性问题,使普通电力用户不必为谐波电能多交电费,但谐波源用户仍仅需支付基波电能电费,而其向供电电网注入的电力谐波污染却未受到任何处罚。
可见,不论是采用全电能计量方式还是使用基波电能计量方式,均未能解决如何对谐波源用户向供电电网注入谐波、间谐波等所造成电能计量不准确的问题,更谈不上如何通过惩罚措施尽可能削弱谐波源用户向供电系统注入谐波等干扰。
因此,针对具体、典型动态非线性负荷造成供电系统电流、电压波形畸变,如何更准确地提取出基波、谐波和间谐波的幅值、相位和频率等参数,是高精度电能计量的一个前提。常用的方法是对被测信号施加窗函数,傅里叶变换、插值修正等运算。具体是在对被测信号进行FFT前,施加合适的窗函数,以抑制因非同步采样和非整周期截断造成的频谱泄漏;再借助插值修正,补偿栅栏效应导致的误差。论文“基于多项式逼近的单峰谱线插值算法在间谐波分析中的应用”(肖先勇等.电网技术,2008),给出了基于Hanning窗、Hamming窗、Blackman-Harris窗、Rife-Vincent窗的单峰谱线插值修正公式;论文“应用FFT进行电力系统谐波分析的改进算法”(庞浩等.中国电机工程学报,2003),给出了基于矩形窗、Triangle窗、Hanning窗、Hamming窗、Blackman窗的双峰插值修正公式;论文“基于三谱线插值FFT的电力谐波分析算法”(牛胜锁等.中国电机工程学报,2012)、论文“基于五项莱夫-文森特窗的三谱线插值FFT谐波分析”(朱俊伟等.电气技术,2015)、论文“基于Nuttall窗的三峰插值谐波算法分析”(翟瑞淼等.电力系统保护与控制,2015),分别给出了基于Hanning窗、Hamming窗、Blackman窗、Blackman-Harris窗、Rife-Vincent窗、Nuttall窗的三峰谱线插值修正公式;论文“基于四谱线插值FFT的电网谐波检测方法”(郝柱等.电力系统保护与控制,2014)和论文“基于Nuttall窗四谱线插值FFT的电力谐波分析”(李得民等.电力系统保护与控制,2016),分别给出了基于Triangle窗、Hanning窗、Hamming窗、Blackman、Nuttall窗的四峰谱线插值修正公式。
但这些已有算法仍存在的问题是:1)由于窗函数的频谱性能是固定不变的,故在对动态信号进行加窗插值修正算法时,计算精度受到窗函数固定旁瓣性能的限制,导致被测信号中的谐波成分由于受到频谱泄露的影响的无法精确提取,从而使电能计量的计算准确度降低;2)单峰插值修正公式较复杂、易受频谱泄漏影响,计算准确度较差;双峰插值算法对峰值频点附近蕴含谱线信息利用不充分;三峰插值算法未考虑峰值频点左右对称谱线的全面信息;四峰插值算法计算精度最高但算法计算量大大增加,难以满足实时电能计量等更高要求;3)现有插值修正算法中,普遍存在偶次谐波电能计算准确度较差的问题。
有鉴于此,有必要提供一种基于Kaiser窗FFT单峰插值的电能计量方法,以满足实际应用需要。
发明内容
本发明旨在至少解决上述技术问题之一。
为此,本发明的目的在于提出一种基于Kaiser窗FFT单峰插值的电能计量方法,利用Kaiser窗函数形状参数的灵活性,通过设置合适的参数获得具有优良的主、旁瓣性能的Kaiser窗函数,以实现单峰插值算法下动态被测电压、电流信号的高精度电能计量的目的。
为了实现上述目的,本发明所采用的技术方案是:一种基于Kaiser窗FFT单峰插值的电能计量方法,其特征在于,包括以下步骤:
S1:基于Kaiser窗FFT单峰插值修正公式、FFT的频移性和窗函数Kaiser的频谱在基波、谐波和间谐波峰值频率点左、右两侧最大谱线和次最大谱线的幅值,利用多项式拟合,得到信号基波、谐波和间谐波的通用幅值、相位和频率插值修正公式;
S2:根据所述通用幅值、相位和频率修正公式计算被测电压信号、电流信号的基波、谐波和间谐波的幅值、相位和频率参量;
S3:根据所述被测电压信号、电流信号的基波、谐波和间谐波的幅值、相位和频率参量分别计算出基波电能、谐波电能、间谐波电能和总电能。
根据本发明实施例的基于Kaiser窗FFT单峰插值修正的电能计量方法,Kaiser窗函数是一组可调的窗函数,通过自由选择形状参数β来调整主瓣宽度和旁瓣高度的比重,β值越大,窗函数频谱的主瓣宽度越宽,旁瓣峰值越小。本专利采用形状参数β=22的Kaiser窗函数,其主瓣宽度是28π/N,旁瓣峰值为-172dB,旁瓣衰减速率14dB/otc,具有较好的主、旁瓣性能。
另外,根据本发明上述实施例的基于Kaiser窗FFT单峰插值修正的电能计量方法,还可以具有如下附加的技术特征:
进一步地,步骤S1进一步包括:
S101:基于Kaiser窗FFT单峰插值修正公式、FFT的频移性和窗函数Kaiser的频谱在基波、谐波和间谐波峰值频率点左、右两侧最大谱线和次最大谱线的幅值表达出单峰谱线幅值比值参数γ的公式;
S102:对比值参数γ进行多项式拟合逼近,拟合出基于Kaiser窗FFT单峰谱线的通用频率修正系数δ的公式;
S103:根据被测信号的幅值公式拟合出通用的幅值修正系数g(τ)的公式;
S104:根据通用的幅值修正系数g(τ)的公式和通用频率修正系数δ的公式得到通用幅值修正公式、通用相位修正公式和通用频率修正公式。
进一步地,步骤S2进一步包括:
S201:对电压、电流信号分别施加Kaiser窗函数,并进行FFT变换以得到所述Kaiser窗函数的频谱;
S202:利用峰值检测技术,分别检测出被测电压、电流信号的基波、谐波以及间谐波峰值点左、右两侧最大谱线和次最大谱线的幅值;
S203:利用单峰谱线幅值比值参数比值公式,分别计算出电压、电流信号的基波、谐波以及间谐波峰值点处的单峰谱线幅值比值参数γ的值;
S204:根据电压、电流信号的基波、谐波以及间谐波峰值点处的单峰谱线幅值比值参数γ的值,并根据所述Kaiser窗FFT单峰插值修正公式计算被测信号基波、谐波以及间谐波峰值点处的频率修正系数δ的值;
S205:根据修正系数δ的值得到被测信号基波、谐波以及间谐波峰值点处的幅值修正系数g(τ)的值;
S206:根据修正系数δ的值和修正系数g(τ)的值插值修正出被测电压、电流信号的基波、谐波和间谐波的幅值、相位和频率参量。
进一步地,通用频率修正系数δ的计算公式为:
δ=-4.62432199+11.5068053·γ-13.5107196·γ2+12.9060092·γ3-8.82476649·γ4+3.97411947·γ5-1.05034078·γ6+0.12321496·γ7。
进一步地,通用幅值修正系数g(τ)的计算公式为:
g(τ)=7.53791319+1.60993789·τ2+0.179661988·τ4+0.0143863034·τ6;
式中,当基波频率小于等于50Hz时,τ=-(1-δ);当基波频率大于50Hz时,τ=δ。
进一步地,通用幅值的修正公式为:
A=yi·g(τ)/N;
其中,如果y1对应的谱线为最大谱线,则i取1;否则,i取2;N为采样数据长度。
与现有技术相比,本发明具有如下的有益效果:本发明的基于Kaiser窗FFT单峰插值修正的电能计量方法,利用Kaiser窗函数形状参数的灵活性,实现动态信号的高精度谐波参量的插值修正,具有较高的电能计量计算准确度。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1是本发明的基于Kaiser窗FFT单峰插值修正的电能计量方法的流程图;
图2是本发明一个实施例的电网基波频率大于或小于50Hz的情况下峰值点左、右两侧最大谱线和次最大谱线的示意图。
图3是Kaiser窗函数在形状参数β分别为0、4、8、22时的频谱特性。随着β的增大,旁瓣峰值越小,但主瓣宽度变宽。当β=22时,具有较低的旁瓣峰值和较快的旁瓣衰减速率,同时具有适中的主瓣宽度。
具体实施方式
下面通过参考附图描述的实施例是示例性的,仅用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性。
参照下面的描述和附图,将清楚本发明的实施例的这些和其他方面。在这些描述和附图中,具体公开了本发明的实施例中的一些特定实施方式,来表示实施本发明的实施例的原理的一些方式,但是应当理解,本发明的实施例的范围不受此限制。相反,本发明的实施例包括落入所附加权利要求书的精神和内涵范围内的所有变化、修改和等同物。
图1是本发明一个实施例的基于Kaiser窗FFT单峰插值修正的电能计量方法的流程图。如图1所示,一种基于Kaiser窗FFT单峰插值修正的电能计量方法,包括以下步骤:
S1:基于Kaiser窗FFT单峰插值修正公式、FFT的频移性和窗函数Kaiser的频谱在基波、谐波和间谐波峰值频率点左、右两侧最大谱线和次最大谱线的幅值得到信号基波、谐波和间谐波的通用幅值、相位和频率插值修正公式。
在本发明的一个实施例中,步骤S1进一步包括:
S101:基于Kaiser窗FFT单峰插值修正公式、FFT的频移性和窗函数Kaiser的频谱在基波、谐波和间谐波峰值频率点左、右两侧最大谱线和次最大谱线的幅值表达出单峰谱线幅值比值参数γ的公式。
S102:对比值参数γ进行多项式拟合逼近,拟合出基于Kaiser窗FFT单峰谱线的通用频率修正系数δ的公式。
S103:根据被测信号的幅值公式拟合出通用的幅值修正系数g(τ)的公式。
S104:根据通用的幅值修正系数g(τ)的公式和通用频率修正系数δ的公式得到通用幅值修正公式、通用相位修正公式和通用频率修正公式。
S2:根据所述通用幅值、相位和频率修正公式计算被测电压信号、电流信号的基波、谐波和间谐波的幅值、相位和频率参量。
在本发明的一个实施例中,步骤S2进一步包括:
S201:对电压、电流信号分别施加Kaiser窗函数,并进行FFT变换以得到所述Kaiser窗函数的频谱;
S202:利用峰值检测技术,分别检测出被测电压、电流信号的基波、谐波以及间谐波峰值点左、右两侧最大谱线和次最大谱线的幅值;
S203:利用单峰谱线幅值比值参数比值公式,分别计算出电压、电流信号的基波、谐波以及间谐波峰值点处的单峰谱线幅值比值参数γ的值;
S204:根据电压、电流信号的基波、谐波以及间谐波峰值点处的单峰谱线幅值比值参数γ的值,并根据所述Kaiser窗FFT单峰插值修正公式计算被测信号基波、谐波以及间谐波峰值点处的频率修正系数δ的值;
S205:根据修正系数δ的值以及基波频率的判定结果得到被测信号基波、谐波以及间谐波峰值点处的幅值修正系数g(τ)的值;
S206:根据修正系数δ的值和修正系数g(τ)的值插值修正出被测电压、电流信号的基波、谐波和间谐波的幅值、相位和频率。
S3:根据所述被测电压信号、电流信号的基波、谐波和间谐波的幅值、相位和频率参量分别计算出基波电能、谐波电能、间谐波电能和总电能。
在本发明的一个实施例中,通用频率修正系数δ的计算公式为:
δ=-4.62432199+11.5068053·γ-13.5107196·γ2+12.9060092·γ3-8.82476649·γ4+3.97411947·γ5-1.05034078·γ6+0.12321496·γ7
在本发明的一个实施例中,通用幅值修正系数g(τ)的计算公式为:
g(τ)=7.53791319+1.60993789·τ2+0.17966199·τ4+0.01438630·τ6
式中,当基波频率小于等于50Hz时,τ=-(1-δ);当基波频率大于50Hz时,τ=δ。
进一步地,通用幅值的修正公式为:
A=yi·g(τ)/N;
其中,如果y1对应的谱线为最大谱线,则i取1;否则,i取2;N为采样数据长度。
为使本领域技术人员进一步理解本发明,将通过以下实施例进行进一步地说明。
(1)基于Kaiser窗FFT单峰插值修正公式
Kaiser窗函数的时域表达式:
式中,I0(β)是第1类变形零阶贝塞尔函数;β是Kaiser函数的形状参数;
式中,α为Kaiser窗的主瓣值和旁瓣值之间的差值,通过设置不同的β值,可以改变Kaiser窗函数的主瓣和旁瓣性能。
通用频率修正系数公式δ,
δ=-4.62432199+11.5068053·γ-13.5107196·γ2+12.9060092·γ3-8.82476649·γ4+3.97411947·γ5-1.05034078·γ6+0.12321496·γ7 (1)
通用幅值修正系数公式g(τ),
g(τ)=7.53791319+1.60993789·τ2+0.17966199·τ4+0.01438630·τ6 (2)
式中,当基波频率小于等于50Hz时,τ=-(1-δ);当基波频率大于50Hz时,τ=δ。
通用幅值修正公式:
A=yi·g(τ)/N (3)
通用相位修正公式:
通用频率修正公式:
f0=k0Δf=(ki-1+δ)fs/N (5)
式中,如果y1对应的谱线为最大谱线,则i取1;否则,i取2。
(2)推导过程
假设含有多次谐波(间谐波)的信号,以采样频率fs进行采样后,得到离散序列x(n):
式中,m表示谐波次数;Am和θm分别表示m次谐波的幅值和相位;f0表示基波频率。
将式(6)改用欧拉公式表示为:
Kaiser窗函数的时域形式为w(n),其离散频谱为则x(n)加窗后的FFT表达式为:
如果忽略负频点-mf0处谱峰的旁瓣影响,在正频点mf0附近的连续频谱函数为:
对式(9)进行采样,可得到其离散傅里叶变换的表达式为:
式中,Δf=fs/N为离散频率间隔;N为采样数据长度。
以基波频率f0=k0Δf为例。由于电网频率存在波动性,且一般不是准确的50Hz,如此,以固定频率对被测信号进行采样,就会发生栅栏效应,即,其最高离散谱线的频率点很难正好是频谱曲线峰值的频率,亦即k0一般不是整数。设频谱曲线峰值点左、右两侧最大谱线和次最大谱线分别为第k1条和第k2条谱线(k1≤k0≤k2),对应幅值分别为如图2所示。引入参数δ=k0-k1,得到δ∈[0,1]。令两条谱线幅值比值参数比值γ::
根据多项式拟合逼近公式polyfit(γ,δ,7),求出反函数δ=f-1(γ),即求出通用频率修正系数δ的公式,见公式(1)。
幅值A的修正公式为:
进一步优化,令τ=k0-ki,当0≤δ≤0.5,τ=δ;当0.5≤δ≤1时,τ=-(1-δ),于是,幅值的修正公式变为:
令g(τ)=2·N/|W[2π(-τ)/N]|,当N较大时通过polyfit(τ,g(τ),7)程序进行多项式逼近,求出幅值修正系数g(τ)的公式,见式(2)。将式(2)带入式(13),于是得到通用幅值修正公式,见公式(3)。
由公式(3)可看出,被测信号实际峰值是利用峰值点处最大谱线的幅值经加权修正获得的,故该算法称为单峰插值修正算法。
(3)电能计算公式
有了被测信号实际基波分量幅值、相位和频率的修正式,就可足够准确地计算出被测实际信号的基波、各次谐波和间谐波分量的幅值、相位和频率,进而就可足够准确地计算出基波、各次谐波和间谐波电能。由三角函数的正交性可知,不同整数次电压电流谐波分量不产生有功电能,故基波、谐波、间谐波的电能计算式为:
式中,Wn为第n次谐波或间谐波的有功电能,当n取1时,W1为基波电能;当n取大于1的整数时,Wn为谐波电能;当n取非整数时,Wn为间谐波电能;Un、In分别为电压信号、电流信号的第n次谐波或间谐波的幅值,当n取1时为基波幅值;αn、βn分别为电压信号、电流信号的第n次谐波或间谐波的初相位,当n取1时为基波相位;T为电网电压电流的基波周期;k为时间窗的个数,按照IEC61000-4-7标准,一般取10个基波周期。
计及基波、谐波、间谐波的总电能计算公式为:
W=CjWj-CposxWposx+CoppxWoppx-CposjxWposjx+CoppjxWoppjx (16)
其中,Wj为基波电能;Wposx为正向谐波电能;Woppx为反向谐波电能;Wposjx正向间谐波电能;Woppjx反向间谐波电能;Cj为基波电能加权系数;Cposx为正向谐波电能加权系数;Coppx为反向谐波电能加权系数;Cposjx为正向间谐波电能加权系数;Coppjx为反向间谐波电能加权系数。建议加权系数取值:Cj=1,0<Cposx、Cposjx<1,Coppx、Coppjx>1。
(4)电能仿真
供电电网电压、电流仿真模型如表1所示,其中含有多次谐波成分,具体为
式中,f0为基波频率,通常在49.7Hz~50.3Hz之间波动,采样频率fs=4000Hz,按照IEC61000-4-7标准,谐波测量的频谱分析时间窗被统一规定为200ms。
表1电压、电流信号的基波及谐波参数
分别采用Hanning窗、Hamming窗、Blackman窗、Blackman-Harris窗以及Kaiser窗(β=22)基于单峰插值算法进行电能计量仿真计算,窗函数的长度均选为800点,基波频率为50.1Hz,仿真计算结果给出在表2中。
表2电能的仿真计算结果
由仿真结果可以看出,基于Kaiser窗FFT单峰插值修正的电能计量算法的计算准确度相比其他四种算法高一个甚至两个数量级,同时偶次谐波电能的准确度也得到大大提高,说明本发明提出的基于Kaiser窗FFT单峰插值修正的电能计量算法的确具有较高的计算准确度。
另外,本发明实施例的基于Kaiser窗FFT单峰插值修正的电能计量方法的其它构成以及作用对于本领域的技术人员而言都是已知的,为了减少冗余,不做赘述。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例,本领域的普通技术人员可以理解:在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由权利要求及其等同限定。
Claims (6)
1.一种基于Kaiser窗FFT单峰插值的电能计量方法,其特征在于,包括以下步骤:
S1:基于Kaiser窗FFT单峰插值修正公式、FFT的频移性和窗函数Kaiser的频谱在基波、谐波和间谐波峰值频率点左、右两侧最大谱线和次最大谱线的幅值,利用多项式拟合,得到信号基波、谐波和间谐波的通用幅值、相位和频率插值修正公式;
S2:根据所述通用幅值、相位和频率修正公式计算被测电压信号、电流信号的基波、谐波和间谐波的幅值、相位和频率参量;
S3:根据所述被测电压信号、电流信号的基波、谐波和间谐波的幅值、相位和频率参量分别计算出基波电能、谐波电能、间谐波电能和总电能。
2.根据权利要求1所述的基于Kaiser窗FFT单峰插值的电能计量方法,其特征在于,步骤S1进一步包括:
S101:基于Kaiser窗FFT单峰插值修正公式、FFT的频移性和窗函数Kaiser的频谱在基波、谐波和间谐波峰值频率点左、右两侧最大谱线和次最大谱线的幅值表达出单峰谱线幅值比值参数γ的公式;
S102:对比值参数γ进行多项式拟合逼近,拟合出基于Kaiser窗FFT单峰谱线的通用频率修正系数δ的公式;
S103:根据被测信号的幅值公式拟合出通用的幅值修正系数g(τ)的公式;
S104:根据通用的幅值修正系数g(τ)的公式和通用频率修正系数δ的公式得到通用幅值修正公式、通用相位修正公式和通用频率修正公式。
3.根据权利要求2所述的基于Kaiser窗FFT单峰插值修正的电能计量方法,其特征在于,步骤S2进一步包括:
S201:对电压、电流信号分别施加Kaiser窗函数,并进行FFT变换以得到所述Kaiser窗函数的频谱;
S202:利用峰值检测技术,分别检测出被测电压、电流信号的基波、谐波以及间谐波峰值点左、右两侧最大谱线和次最大谱线的幅值;
S203:利用单峰谱线幅值比值参数比值公式,分别计算出电压、电流信号的基波、谐波以及间谐波峰值点处的单峰谱线幅值比值参数γ的值;
S204:根据电压、电流信号的基波、谐波以及间谐波峰值点处的单峰谱线幅值比值参数γ的值,并根据所述Kaiser窗FFT单峰插值修正公式计算被测信号基波、谐波以及间谐波峰值点处的频率修正系数δ的值;
S205:根据修正系数δ的值得到被测信号基波、谐波以及间谐波峰值点处的幅值修正系数g(τ)的值;
S206:根据修正系数δ的值和修正系数g(τ)的值插值修正出被测电压、电流信号的基波、谐波和间谐波的幅值、相位和频率参量。
4.根据权利要求2所述的基于Kaiser窗FFT单峰插值修正的电能计量方法,其特征在于,通用频率修正系数δ的计算公式为:
δ=-4.62432199+11.5068053·γ-13.5107196·γ2+12.9060092·γ3-8.82476649·γ4+3.97411947·γ5-1.05034078·γ6+0.12321496·γ7。
5.根据权利要求4所述的基于Kaiser窗FFT单峰插值修正的电能计量方法,其特征在于,通用幅值修正系数g(τ)的计算公式为:
g(τ)=7.53791319+1.60993789·τ2+0.179661988·τ4+0.0143863034·τ6;
式中,当基波频率小于等于50Hz时,τ=-(1-δ);当基波频率大于50Hz时,τ=δ。
6.根据权利要求5所述的基于Kaiser窗FFT单峰插值修正的电能计量方法,其特征在于,通用幅值的修正公式为:
A=yi·g(τ)/N;
其中,如果y1对应的谱线为最大谱线,则i取1;否则,i取2;N为采样数据长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611157854.6A CN106814246A (zh) | 2016-12-15 | 2016-12-15 | 基于Kaiser窗FFT单峰插值的电能计量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611157854.6A CN106814246A (zh) | 2016-12-15 | 2016-12-15 | 基于Kaiser窗FFT单峰插值的电能计量方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106814246A true CN106814246A (zh) | 2017-06-09 |
Family
ID=59109725
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611157854.6A Pending CN106814246A (zh) | 2016-12-15 | 2016-12-15 | 基于Kaiser窗FFT单峰插值的电能计量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106814246A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107945147A (zh) * | 2017-12-04 | 2018-04-20 | 云南电网有限责任公司普洱供电局 | 一种基于Kaiser窗FFT四峰插值修正的图像融合方法 |
CN107991591A (zh) * | 2017-12-04 | 2018-05-04 | 云南电网有限责任公司普洱供电局 | 一种基于Kaiser窗FFT单峰插值修正的图像融合方法 |
CN109581045A (zh) * | 2018-12-24 | 2019-04-05 | 中国船舶重工集团公司第七〇九研究所 | 一种满足iec标准框架的间谐波功率计量方法 |
CN113702702A (zh) * | 2021-08-23 | 2021-11-26 | 中国电力科学研究院有限公司 | 一种改进频谱叠加算法的多间谐波识别方法 |
CN114167198A (zh) * | 2021-10-18 | 2022-03-11 | 国网山东省电力公司平原县供电公司 | 一种同期线损数据的测量方法及测量平台 |
CN114290939A (zh) * | 2021-12-29 | 2022-04-08 | 国网新疆电力有限公司经济技术研究院 | 直流充电桩电压误差补偿方法、系统、设备及可读存储介质 |
-
2016
- 2016-12-15 CN CN201611157854.6A patent/CN106814246A/zh active Pending
Non-Patent Citations (3)
Title |
---|
徐文佳 等: ""计及间谐波的电能计量方法"", 《电网技术》 * |
肖先勇 等: ""基于多项式逼近的单峰谱线插值算法在间谐波分析中的应用"", 《电网技术》 * |
高云鹏 等: ""基于 Kaiser 窗双谱线插值 FFT 的谐波分析方法"", 《仪器仪表学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107945147A (zh) * | 2017-12-04 | 2018-04-20 | 云南电网有限责任公司普洱供电局 | 一种基于Kaiser窗FFT四峰插值修正的图像融合方法 |
CN107991591A (zh) * | 2017-12-04 | 2018-05-04 | 云南电网有限责任公司普洱供电局 | 一种基于Kaiser窗FFT单峰插值修正的图像融合方法 |
CN109581045A (zh) * | 2018-12-24 | 2019-04-05 | 中国船舶重工集团公司第七〇九研究所 | 一种满足iec标准框架的间谐波功率计量方法 |
CN109581045B (zh) * | 2018-12-24 | 2020-12-11 | 中国船舶重工集团公司第七一九研究所 | 一种满足iec标准框架的间谐波功率计量方法 |
CN113702702A (zh) * | 2021-08-23 | 2021-11-26 | 中国电力科学研究院有限公司 | 一种改进频谱叠加算法的多间谐波识别方法 |
CN113702702B (zh) * | 2021-08-23 | 2024-01-23 | 中国电力科学研究院有限公司 | 一种改进频谱叠加算法的多间谐波识别方法 |
CN114167198A (zh) * | 2021-10-18 | 2022-03-11 | 国网山东省电力公司平原县供电公司 | 一种同期线损数据的测量方法及测量平台 |
CN114167198B (zh) * | 2021-10-18 | 2024-03-01 | 国网山东省电力公司平原县供电公司 | 一种同期线损数据的测量方法及测量平台 |
CN114290939A (zh) * | 2021-12-29 | 2022-04-08 | 国网新疆电力有限公司经济技术研究院 | 直流充电桩电压误差补偿方法、系统、设备及可读存储介质 |
CN114290939B (zh) * | 2021-12-29 | 2024-04-23 | 国网新疆电力有限公司经济技术研究院 | 直流充电桩电压误差补偿方法、系统、设备及可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106814246A (zh) | 基于Kaiser窗FFT单峰插值的电能计量方法 | |
CN106814230A (zh) | 基于Kaiser窗FFT四峰插值的电能计量方法 | |
CN106526312A (zh) | 基于r‑v (ⅲ)窗fft双峰插值的电能计量方法 | |
CN201352236Y (zh) | 基于Kaiser窗双谱线插值FFT的谐波电能表 | |
CN203149027U (zh) | 一种基于能量算子和频谱校正的电压波动与闪变检测装置 | |
CN107179439A (zh) | 基于五项余弦窗fft三峰插值的电能计量方法 | |
CN109030941A (zh) | Hanning自乘卷积窗FFT三谱线插值谐波分析方法 | |
CN102662106A (zh) | 谐波电网电能计量方法 | |
CN104360235B (zh) | 含dg复杂配电网的电压暂降源定位方法 | |
Li et al. | Research of harmonic distortion power for harmonic source detection | |
CN104897961A (zh) | 基于互乘法窗函数的三谱线插值fft谐波分析方法及系统 | |
CN103809148B (zh) | 一种测量电能表动态特性的方法及系统 | |
CN102938130A (zh) | 一种基于PQView的电能质量性能指标综合评估方法 | |
CN107543962A (zh) | 主导间谐波频谱分布的计算方法 | |
CN103472333A (zh) | 风电并网电能质量综合性能检测方法 | |
CN108535613A (zh) | 一种基于组合窗函数的电压闪变参数检测方法 | |
CN103399204A (zh) | 一种基于Rife-Vincent(II)窗插值FFT的谐波与间谐波检测方法 | |
CN106154040A (zh) | 一种用于计算电网单点等值时等效谐波阻抗方法 | |
CN102185314A (zh) | 电能质量综合评判方法 | |
CN105006820A (zh) | 一种电能质量对输电线路综合能耗影响的定量分析方法 | |
CN103543317A (zh) | 用于测量iec标准电压闪变的方法 | |
CN106405287B (zh) | 适合日常电能质量监测系统应用的主导谐波源追踪方法 | |
CN106771552B (zh) | 一种畸变功率计量方法 | |
CN106849792A (zh) | 电机设备及集群系统的能耗计算和节能措施评估方法 | |
CN106709623A (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170609 |