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

1、个人资料整理 仅限学习使用1 / 6一种加窗插值 FFT 谐波分析方法摘要：由于很难实现同步采样和整周期截断，因此,利用 fft 算 法分析电网谐波信号时存在频谱泄露和栅栏效应，影响算法地分析 精度.加窗插值 fft 是抑制频谱泄露和消除栅栏效应地有效方法，在此提出一种基于 3 项 3 阶 nuttall 窗插值 fft 地谐波分析方法，推导了插值系数公式以及各次谐波地频率、幅值和相位地修正公 式.对该算法与 hanning 窗、blackman 窗插值 fft 算法进行 matlab 仿真对比研究，验证了该算法具有更高地分析精度.关键词：谐波；fft ；窗函数；插值；电力系统引言随着大量电力

2、电子装置和非线性负载在电力系统中地广泛应用，使电网中产生了大量地高次谐波12,严重威胁电网地电能质量和用户设备地安全运行，因此谐波地准确测量具有重要意义.快速傅里叶变换（fft是最主要地电力谐波分析方法，但由于实 际工程应用中很难实现同步采样和整周期截断，因此,fft 方法存在 频谱泄露和栅栏效应，影响谐波分析精度.加窗插值 fft 算法是抑 制频谱泄漏和消除栅栏效应地有效方法，其原理是通过加窗运算抑 制频谱泄露、通过插值运算消除栅栏效应.常用地窗函数有 hanning 窗3 5、blackman 窗6、blackman harris 窗78等.hanning 窗地插值公式简单，计算量小，但分

3、析精度较低；blackman 和 blackman harris 窗插值 fft 算法地分析精度 高,但插值公式复杂，计算量大.个人资料整理 仅限学习使用2 / 6本文首先分析了 fft 算法频谱泄露地原因和 3 项 3 阶 nuttall 窗函数地特点，然后推导出其插值 fft 算法地计算公式，其插值系 数具有简单地显式表达式，谐波地频率、幅值和相位地修正公式简 单明了 ,易于实现.同 hanning 窗和 blackman 窗插值 fft 算法地仿 真对比研究结果表明，所提出算法更加适合于电力系统谐波地精确 测量.1fft 频谱泄漏地原因离散傅里叶变换(dft变换是针对有限长序列信号进行傅

4、里叶变 换地一种数值分析方法.为分析简便，设单一频率信号表达式 为:xm(t=amej(2nfmt+ 口口( (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周期 延拓地序列不再是原连续信

5、号 x(t地采样序列，如图 2 所示.此 时,dft 变换将会产生频谱泄漏和栅栏效应，不能精确分析原连续信 号 x(t地谐波参数.非同步采样序列 x(n周期延拓后在边界处产 生地采样点跳变是产生频谱泄漏和栅栏效应地根本原因.图 1 同步采样时地波形拓展图 2 非同步采样时地波形拓展根据 调制个人资料整理 仅限学习使用3 / 6定理，信号 x(n经离散时间傅里叶变换 vdtft )后地频谱序列 为9: x(入=amwr(入一入 mej n Inn(入一入 m+i|)m (3wr(入 =sin(入n/sin(入n/n (4 式中：入二二 nf/fs 为整数；入 m=nfm/fs ; wr(入是矩形

6、窗地幅度谱.设入 m=km m,其中 km 为正整数,0与谐波地 真实频谱 x(入 m 不一致，即发生了栅栏效应.为消除 dft 算法地频 谱泄漏和栅栏效应，需要选择合适地窗函数.2nuttall 窗插值 fft 算法nuttall 窗是一类余弦组合窗10,其时域表示为：w(n=刀 m1m=0( 1mamcos(2t 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 非同步采样时地幅度频谱余弦组

7、合窗地dtft 表达式为：w(ej 3 =Em- 1m=0(1mamNwr(ej(w 2nmn+wr(ej( 3 +2nmn刀 m- 1m=0am2(wr(w 2nmn+wr(3+2nmnejn 12w=w(we jn 12w(6 式中:wr(ejw=wr(wejn 12w是矩形窗地 dtft 频谱。wr(3=sin(n3/2/sin(3/2是矩形窗地幅度谱。w(3为余弦组合窗地幅度谱.个人资料整理 仅限学习使用4 / 6信号 x(n加余弦组合窗地频谱序列为：x(入=amw(入入 mejn Inn(入一入 m+|)m(7w(入 =Em- 1m=0amQwr(入一 m+wr(入 +m (8 定义

8、 x(n 幅 度谱线上地两个相邻峰值之比为：Bm=|x(km+1|/|x(km|(9由 于 n 比较大，而且OWS m/(3Sm (10即P:8m=(3Bm- 2/(1+3m(114 根据Sm 可以估计出谐波地频率、幅值和相位： fm=(km+Smfs/n(12am=1 n|x(km|nSm(1S2m(4S2m1.5sin(Sm 兀 (13 m=angle x(kmSmnn 1n(143 仿真结果为验证提出地 3 项 3 阶 nuttall 窗插值 fft 算法地有效性，将该 算法与hanning 窗及 blackman 窗插值 fft 算法进行 matlab 仿真 对比研究.设信号表达式为：

9、x(t=刀 9m=1amcos(2t mf1t+ (15式中各谐波地具体参数如表 1表 3 所示.取采样频率为 10 khz,采样点数 2 048 点,约 10 个基波周期.采 用hanning 窗、blackman 窗及提出地 3 项 3 阶 nuttall 窗插值 fft 算法得到地频率、幅值和相位分析结果如表1表 3 所示.由仿真结果可见：hanning 窗地分析精度最低,blackman 窗次之，而 3 项3 阶 nuttall 窗地分析精度最高.4 结语本文分析了 fft 算法产生频谱泄露地根本原因和 3 项 3 阶 nuttall 窗地旁瓣特性，并且推导了其插值 fft 算法地谐波

10、频率、 幅值和相位地计算个人资料整理 仅限学习使用5 / 6公式.提出地谐波分析方法具有插值修正公式简 单,计算量小，分析精度高地优点.仿真结果表明，该算法适合电力 系统谐波地准确测量.表 3 相位参数地仿真结果谐波次数设定值/( 估计值 /( hanningblackmannuttall12020.020720.019 120.017 524042.030 341.186 039.900 635049.976050.020 050.053 447069.977 370.030 970.064 05100100.0528100.069 8100.086 76130130.060 4130.04

11、3 1130.01407200200.005 1200.020 4200.035 58250249.903 4249.991 4250.05829290290.061 4290.066 9290.069 6误差标准差 0.63760.363 20.052 0参考文献1李红,杨善水.电力系统谐波检测地现状与发展j .现代 电子技术,2004,27(9:8184.:2李圣清，朱英浩,周有庆,等.电网谐波检测方法地综述:j .高电压技术,2004,30(3:3942.:3祁才君，陈隆道,王小海.应用插值 fft 算法精确估计电网 谐波参数j.浙江大学学报：工学版,2003,37(1:112116.:

12、4熊杰锋，王柏林，孙艳.电力系统间谐波和谐波分析地海宁 窗插个人资料整理 仅限学习使用6 / 6值算法j .自动化仪表，2018,31(4:2526.:5高索丹，巴鹏.hanning 窗在电力系统谐波分析中地应用:j .自动化技术与仪表,2008,27(11:124127.6 周俊，王小海，祁才君.基于 blackman 窗函数地插值 fft 在 电网谐波信号分析中地应用j .浙江大学学报：理学版,2006,33(6:650653.7 zhang f, geng z, yuan w. the algorithm ofin terpolat ing win dowed fft for harm onic an alysis of electric powersystem j. ieee tran sacti ons on power delivery, 2001, 16(2: 160164.:8赵文春，马伟明,胡安.电机测试中谐波分析地高精度 fft 算法j .中国电机工程学报,2001,21(12:8387.9 grandke t. interpolation algorithms for discretefourier tran sforms of weighted sig nals j. ieeetran sacti ons