




已阅读5页,还剩2页未读, 继续免费阅读
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
一种加窗插值FFT谐波分析方法摘要:由于很难实现同步采样和整周期截断,因此,利用fft算法分析电网谐波信号时存在频谱泄露和栅栏效应,影响算法的分析精度。加窗插值fft是抑制频谱泄露和消除栅栏效应的有效方法,在此提出一种基于3项3阶nuttall窗插值fft的谐波分析方法,推导了插值系数公式以及各次谐波的频率、幅值和相位的修正公式。对该算法与hanning窗、blackman窗插值fft算法进行matlab仿真对比研究,验证了该算法具有更高的分析精度。关键词:谐波; fft; 窗函数; 插值; 电力系统引言随着大量电力电子装置和非线性负载在电力系统中的广泛应用,使电网中产生了大量的高次谐波12,严重威胁电网的电能质量和用户设备的安全运行,因此谐波的准确测量具有重要意义。快速傅里叶变换(fft)是最主要的电力谐波分析方法,但由于实际工程应用中很难实现同步采样和整周期截断,因此,fft方法存在频谱泄露和栅栏效应,影响谐波分析精度。加窗插值fft算法是抑制频谱泄漏和消除栅栏效应的有效方法,其原理是通过加窗运算抑制频谱泄露、通过插值运算消除栅栏效应。常用的窗函数有hanning窗35、blackman窗6、blackmanharris窗78等。hanning窗的插值公式简单,计算量小,但分析精度较低;blackman和blackmanharris窗插值fft算法的分析精度高,但插值公式复杂,计算量大。本文首先分析了fft算法频谱泄露的原因和3项3阶nuttall窗函数的特点,然后推导出其插值fft算法的计算公式,其插值系数具有简单的显式表达式,谐波的频率、幅值和相位的修正公式简单明了,易于实现。同hanning窗和blackman窗插值fft算法的仿真对比研究结果表明,所提出算法更加适合于电力系统谐波的精确测量。1fft频谱泄漏的原因离散傅里叶变换(dft)变换是针对有限长序列信号进行傅里叶变换的一种数值分析方法。为分析简便,设单一频率信号表达式为:xm(t)=amej(2fmt+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)ejn1n(m)+m (3)wr()=sin()/sin(/n) (4)式中:=nf/fs为整数;m=nfm/fs;wr()是矩形窗的幅度谱。设m=km+m,其中km为正整数,0m1,则非同步采样时的幅度频谱如图3所示,可见,此时整数位置上的频谱x(km)与谐波的真实频谱x(m)不一致,即发生了栅栏效应。为消除dft算法的频谱泄漏和栅栏效应,需要选择合适的窗函数。2nuttall窗插值fft算法nuttall窗是一类余弦组合窗10,其时域表示为:w(n)=m1m=0(1)mamcos(2nm/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)=m1m=0(1)mam2wr(ej(2mn)+wr(ej(+2mn)m1m=0am2(wr(2mn)+wr(+2mn)ejn12=w()ejn12 (6)式中:wr(ej)=wr()ejn12是矩形窗的dtft频谱;wr()=sin(n/2)/sin(/2)是矩形窗的幅度谱;w()为余弦组合窗的幅度谱。信号x(n)加余弦组合窗的频谱序列为:x()=amw(m)ejn1n(m)+m (7)w()=m1m=0am2wr(m)+wr(+m) (8)定义x(n)幅度谱线上的两个相邻峰值之比为:m=|x(km+1)|/|x(km)|(9)由于n比较大,而且0m1,对于3项3阶nuttall窗,可得:m=(2+m)/(3m) (10)即:m=(3m2)/(1+m)(11)根据m可以估计出谐波的频率、幅值和相位:fm=(km+m)fs/n(12)am=1n|x(km)|m(12m)(42m)1.5sin(m)(13)m=anglex(km)mn1n(14)3仿真结果为验证提出的3项3阶nuttall窗插值fft算法的有效性,将该算法与hanning窗及blackman窗插值fft算法进行matlab仿真对比研究。设信号表达式为:x(t)=9m=1amcos(2mf1t+m)(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算法的谐波频率、幅值和相位的计算公式。提出的谐波分析方法具有插值修正公式简单,计算量小,分析精度高的优点。仿真结果表明,该算法适合电力系统谐波的准确测量。表3相位参数的仿真结果谐波次数设定值/()估计值 /()hanningblackmannuttall12020.020 720.019 120.017 524042.030 341.186 039.900 635049.976 050.020 050.053 447069.977 370.030 970.064 05100100.052 8100.069 8100.086 76130130.060 4130.043 1130.014 07200200.005 1200.020 4200.035 58250249.903 4249.991 4250.058 29290290.061 4290.066 9290.069 6误差标准差0.637 60.363 20.052 0参考文献1李红,杨善水.电力系统谐波检测的现状与发展j.现代电子技术,2004,27(9):8184.2李圣清,朱英浩,周有庆,等.电网谐波检测方法的综述j.高电压技术,2004,30(3):3942.3祁才君,陈隆道,王小海.应用插值fft算法精确估计电网谐波参数j.浙江大学学报:工学版,2003,37(1):112116.4熊杰锋,王柏林,孙艳.电力系统间谐波和谐波分析的海宁窗插值算法j.自动化仪表,2010,31(4):2526.5高索丹,巴鹏.hanning窗在电力系统谐波分析中的应用j.自动化技术与仪表,2008,27(11):124127.6周俊,王小海,祁才君.基于blackman窗函数的插值fft在电网谐波信号分析中的应用j.浙江大学学报:理学版,2006,33(6):650653.7zhang f, geng z, yuan w. the algorithm of interpolating windowed fft for harmonic analysis of electric power system j. ieee transactions on power delivery, 2001, 16(2): 160164.8赵文春,马伟明,胡安.电机测试中谐波分析的高精度fft算法j.中国电机工程学报,2001,21(12):8387.9grandke t. interpolation algorithms for discrete fourier transforms of weighted signals j. ieee transactions on instrumen
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论