一种加窗插值FFT谐波分析方法.docVIP

一种加窗插值FFT谐波分析方法 摘要：由于很难实现同步采样和整周期截断，因此，利用fft算法分析电网谐波信号时存在频谱泄露和栅栏效应，影响算法的分析精度。加窗插值fft是抑制频谱泄露和消除栅栏效应的有效方法，在此提出一种基于3项3阶nuttall窗插值fft的谐波分析方法，推导了插值系数公式以及各次谐波的频率、幅值和相位的修正公式。对该算法与hanning窗、blackman窗插值fft算法进行matlab仿真对比研究，验证了该算法具有更高的分析精度。 关键词：谐波； fft； 窗函数； 插值； 电力系统 引言 随着大量电力电子装置和非线性负载在电力系统中的广泛应用，使电网中产生了大量的高次谐波［12］，严重威胁电网的电能质量和用户设备的安全运行，因此谐波的准确测量具有重要意义。 快速傅里叶变换(fft)是最主要的电力谐波分析方法，但由于实际工程应用中很难实现同步采样和整周期截断，因此，fft方法存在频谱泄露和栅栏效应，影响谐波分析精度。加窗插值fft算法是抑制频谱泄漏和消除栅栏效应的有效方法，其原理是通过加窗运算抑制频谱泄露、通过插值运算消除栅栏效应。常用的窗函数有hanning窗［35］、blackman窗［6］、blackmanharris窗［78］等。hanning窗的插值公式简单，计算量小，但分析精度较低；blackman和blackmanharris窗插值fft算法的分析精度高，但插值公式复杂，计算量大。 本文首先分析了fft算法频谱泄露的原因和3项3阶nuttall窗函数的特点，然后推导出其插值fft算法的计算公式，其插值系数具有简单的显式表达式，谐波的频率、幅值和相位的修正公式简单明了，易于实现。同hanning窗和blackman窗插值fft算法的仿真对比研究结果表明，所提出算法更加适合于电力系统谐波的精确测量。 1fft频谱泄漏的原因 离散傅里叶变换(dft)变换是针对有限长序列信号进行傅里叶变换的一种数值分析方法。为分析简便，设单一频率信号表达式为:xm(t)=amej(2πfmt+φm)(1)式中am，fm，φm分别为幅值、频率和相位。 以采样频率fs对其进行均匀采样n点，得到离散序列x(n)： x(n)=xm(n)·wr(n)，n=0,1,2,…,n-1(2) 式中wr(n)为矩形窗。 应用dft变换对采样信号x(n)进行分析时，隐含在时域上对其进行周期延拓。在理想同步采样条件下，x(n)周期延拓后的序列与原连续信号x(t)的采样序列完全相同，如图1所示。此时，dft变换能够精确分析原连续信号x(t)的谐波参数。非同步采样时，x(n)周期延拓的序列不再是原连续信号x(t)的采样序列，如图2所示。此时，dft变换将会产生频谱泄漏和栅栏效应，不能精确分析原连续信号x(t)的谐波参数。非同步采样序列x(n)周期延拓后在边界处产生的采样点跳变是产生频谱泄漏和栅栏效应的根本原因。 图1同步采样时的波形拓展图2非同步采样时的波形拓展根据调制定理，信号x(n)经离散时间傅里叶变换（dtft）后的频谱序列为［9］：x(λ)=amwr(λ－λm)ej［－n－1nπ(λ－λm)+φm］ (3) wr(λ)=sin(λπ)/sin(λπ/n) (4)式中：λ=nf/fs为整数；λm=n·fm/fs；wr(λ)是矩形窗的幅度谱。 设λm=km+δm，其中km为正整数，0≤δm＜1，则非同步采样时的幅度频谱如图3所示，可见，此时整数位置上的频谱x(km)与谐波的真实频谱x(λm)不一致，即发生了栅栏效应。为消除dft算法的频谱泄漏和栅栏效应，需要选择合适的窗函数。 2nuttall窗插值fft算法 nuttall窗是一类余弦组合窗［10］，其时域表示为：w(n)=∑m－1m=0(－1)mamcos(2πn·m/n)(5)式中：m为窗函数的项数，n=0,1,2,…,n-1。 当窗函数为3项3阶nuttall窗时：a0=0.375，a1=0.5，a2=0.125，其旁瓣衰减为47 db，旁瓣衰减速度为30 db，适合电力谐波的准确分析。 图3非同步采样时的幅度频谱余弦组合窗的dtft表达式为：w(ejω)=∑m－1m=0(－1)mam2［wr(ej(ω－2πmn))+wr(ej(ω+2πmn))］ ≈∑m－1m=0am2(wr(ω－2πmn)+wr(ω+2πmn))·e－jn－12ω =w(ω)e－jn－12ω (6)式中:wr(ejω)=wr(ω)e－jn－12ω是矩形窗的dtft频谱;wr(ω)=sin(nω/2)/sin(ω/2)是矩形窗的幅度谱;w(ω)为余弦组合窗的幅度谱。 信号x(n)加余弦组合窗的频谱序列为：x(λ)=amw(λ－λm)ej［－n－1nπ(λ－λm)+φm］ (7) w(λ)=∑m－1m=0