凯塞窗四谱线插值FFT的电力谐波分析方法

　　摘要：针对在采用快速傅里叶变换分析电力谐波时，对信号非同步采样以及非整周期截断的情况下，会产生频谱泄露和栅栏效应，无法有效检测谐波的各参数值等问题，提出了加窗插值算法，能够提高计算精度抑制频谱泄露，Kaiser窗的频带内能量主要集中在主瓣中，拥有绝佳的旁瓣抑制性能，该算法充分利用峰值谱线频点最近的四条谱线进行加权运算提高运算精度，运用MATLAB多项式拟合推导出实用的插值修正公式。仿真结果表明，四谱线插值算法得到的21次谐波的幅值相对误差0。00001，相位相对误差0。000001，基波频率计算相对误差为6。21011，并且能够有效消除频谱泄露和栅栏效应的影响，计算精度更高。
　　关键词：谐波分析；Kaiser窗；四谱线插值；FFT
　　DOI：10。15938j。jhust。2018。02。010
　　中图分类号：TM935
　　文献标志码：A
　　文章编号：10072683（2018）02005305
　　Abstract：WhentheFastFourierTransform（FFT）methodisusedtoanalyzetheharmonicsinapowernetwork，thefastFouriertransform（FFT）canproducespectralleakageandfenceeffectunderthesituationofnonsynchronizedsamplingandnonintegralperiodtruncation，whichleadstoagreaterrorintheresultsofcalculation。Thewindowfunctionsandinterpolationalgorithmcanreducespectralleakageandimprovetheaccuracyofcomputation。ByMATLABdatafitting，practicalcorrectionformulascorrespondingtocommonlyusedwindowfunctionsareobtainedbyusingweightedarithmeticoffourpillarsspectrallineneartheharmonicfrequencypoint。Thesimulationandtestresultsshowthattherelativeerrorsofamplitudefor21orderharmonicsarenomorethan0。00001，therelativeerrorsofphrasecollectionfor21orderharmonicsarenomorethan0。000001，andtherelativeerrorsoffundamentalwavefrequencyare6。21011。Thealgorithmcaneliminatetheinfluenceofspectralleakageandfenceeffecteffectively，improvingtheaccuracyofharmonicanalysis。
　　Keywords：harmonicanalysis；Kaiserwindow；fourspectrumlineinterpolation；FFT
　　0引言
　　在供配电网络中，存在着含有大量的非线性、冲击性负荷的电力电子设备，导致电网谐波污染越来越严重。为了维护电网的安全稳定运行，获得实时的电网质量参数，对谐波参数进行高精度的检测显得尤其重要。常用的谐波分析方法是快速傅里叶变换（FFT），在非同步采样时会产生频谱泄露和栅栏效应，影响谐波分析的精度。采用性能良好的窗函数可以有效抑制频谱泄露，对频域信号进行插值运算可减少栅栏效应引起的误差〔12〕。常用的窗函数有三角窗、Hanning窗〔3〕、BlackmanHarris窗〔4〕、RifeVincent窗〔5〕、Nuttall窗〔67〕等。在插值算法中，常用的有纹紫卟逯邓惴〔89〕，双谱线插值算法〔1014〕，三次样条插值算法〔1516〕，三谱线插值算法等高精度插值算法〔1719〕。
　　传统的窗函数无法随意改变主旁瓣能量之间的比重，而Kaiser窗较为灵活〔20〕。本文对Kaiser窗的频谱特性进行分析，在双谱线插值和三谱线插值的基础上，提出了基于Kaiser窗四谱线插值FFT的电力谐波分析方法，同时为了计算的方便快捷，运用MATLAB多项式拟合推导出实用的插值修正公式，并对本文提出的算法进行了验证。
　　1Kaiser窗特性
　　Kaiser窗由零阶贝塞尔函数构成，可自由选择主瓣宽度和旁瓣高度之间的比重，频带内能量主要集中在主瓣中，其时域表达式为：
　　图1（a）、（b）给出了〔0，4，8，11〕时Kaiser窗的时域和频域特性。由时域图形可见，当0时，Kaiser窗变为矩形窗函数，其旁瓣渐进衰减速率为6dBoct，当8时，旁瓣峰值为58。3dB，旁瓣渐进衰减速率为12dBoct，随着值的增大，旁瓣峰值继续下降，旁瓣渐进衰减速率增加。本文算法选择20。
　　2基于Kaiser窗的四谱线插值算法
　　3仿真验证及结果分析
　　3。1简单信号分析为验证基于Kaiser窗的四谱线插值算法的精确性，采用已被众多文献验证的无噪声信号模型，其表达式为：
　　本文选取Hanning窗、Kaiser窗双谱线、三谱线和四谱线（20，N1024）进行仿真。
　　从表1数据可以看出，相比于Hanning窗、Kaiser双谱线、Kaiser三谱线，本文Kaiser窗四谱线插值算法具有良好的抑制频谱泄露的作用，并且精度更高。
　　3。2利用复杂信号对算法进行的仿真验证
　　为了便于比较分析，本文采用了已被众多文献验证的包含21次谐波的复杂信号模型：
　　依次用Kaiser窗双谱线插值算法、三谱线插值算法和本文提出的四谱线插值算法对谐波信号进行处理（双谱线插值采用文〔10〕的修正公式及结果，三谱线插值采用文〔17〕提出的插值算法），仿真结果如表3、表4所示。其中表3为基波与各次谐波幅值的相对误差，最后一项为基波频率的相对误差；表4为基波与各次谐波相位的相对误差。
　　从表3、表4仿真结果可以看出，在其他条件相同且20的情况下，本文四谱线插值算法得到的21次谐波的幅值相对误差0。00001，相位相对误差0。000001，由于本文算法流程，只求一次基波频率，谐波频率根据基波频率的倍数求得，所以基波频率的精度非常重要，由表3可得双谱线、三谱线和四谱线基波频率计算相对误差分别为1。4109、41010、6。21011，四谱线插值算法比双谱线和三谱线精确度更高。可见本文提出的Kaiser窗四谱线插值算法与双谱线和三谱线相比，由于提高了修正公式的精确度，从而显著的提高了谐波参数检测准确性。
　　4结论
　　由于Kaiser窗拥有绝佳的频谱泄露抑制效果，并且，四谱线插值较高的运算精度，本文提出了基于Kaiser窗四谱线插值FFT的电力谐波分析方法。并且为了计算方便，根据MATLAB多项式拟合求出了实用的插值修正公式。通过对简单信号和复杂信号在不同插值算法的仿真对比，结果表明本文所提算法在对采样数据进行非整周期截断情况下，与Kaiser窗的双谱线和三谱线插值算法相比，具有更高的计算准确度。
　　参考文献：
　　〔1〕许珉，杨阳，陈飞，等。基于Nuttall（I）窗的插值FFT算法〔J〕。电力系统保护与控制，2011，39（23）：4447。
　　〔2〕曾博，唐求，卿柏元，等，基于Nuttall自卷积窗的改进FFT谱分析方法〔J〕。电工技术学报，2014，29（7）：5965。
　　〔3〕BARROSJ，DIEGORI。OntheUseoftheHanningWindowforHarmonicAnalysisintheStandardFramework〔J〕。IEEETransonPowerDelivery，2006，21（1）：538539。
　　〔4〕HARRISFJ。OntheUseofWindowsforHarmonicAnalysiswiththeDiscreteFourierTransform〔J〕。ProceedingsoftheIEEE，1978，66（1）：5183。
　　〔5〕RIFEDC，VINCENTGA。UseoftheDiscreteFourierTransformintheMeasurementofFrequenciesandLevelsofTones〔J〕。TheBellSystemTechnicalJournal，1970，49（2）：197228。
　　〔6〕NUTTALLAH。SomeWindowswithaVeryGoodSideLobeBehavior〔J〕。IEEETrans。onAcoustics，SpeechandSignalProcessing，1981，29（1）：8491。
　　〔7〕卿柏元，滕召胜，高云鹏，等。基于Nuttall窗双谱线插值FFT的电力谐波分析方法〔J〕。中国电机工程学报，2008，28（25）：153158。
　　〔8〕王楠，肖先勇，曾颂崎，等。基于插值线性调频Z变换的间谐波分析方法〔J〕。电网技术，2007，31（18）：4347。
　　〔9〕肖先勇，王楠，刘亚梅。基于多项式逼近的单峰谱线插值算法在间谐波分析中的应用〔J〕。电网技术，2008，32（18）：5761。
　　〔10〕高云鹏，腾召胜，卿柏元。基于Kaiser窗双谱线插值FFT的谐波分析方法〔J〕。仪器仪表学报，2010，31（2）：287292。
　　〔11〕高云鹏，滕召胜，温和，等。凯塞窗差值FFT的电力谐波分析与应用〔J〕。中国电机工程学报，2010，30（4）：4348。
　　〔12〕曾博，滕邵胜，高云鹏，等。基于RifeVincent窗的高嗜范鹊缌巢嗔考扑惴椒〔J〕。电工技术报，2009，24（8）：154159。
　　〔13〕庞浩，李东霞，俎云霄，等。应用FFT进行电力系统谐波分析的改进算法〔J〕。中国电机工程学报，2003，23（6）：5054。
　　〔14〕薛蕙，杨仁刚。基于FFT的高精度谐波检测算法〔J〕。中国电机工程学报，2002，22（12）：106110。
　　〔15〕王继东，杜旭浩，杨帆。基于三次样条插值信号重构的微网谐波及间谐波分析算法〔J〕。电网技术，2012，36（11）：711。
　　〔16〕杨晓萍，刘普森，钟彦儒。基于经验模式分解的有源滤波器谐波检测〔J〕。电工技术学报，2009，24（5）：197202。
　　〔17〕陈立伟，刘欢，冯驰，等。Kaiser窗三谱线插值电力谐波分析〔J〕。应用科技，2014（2）：511。
　　〔18〕牛胜锁，梁志瑞，张建华，等。基于三谱线插值FFT的电力谐波分析算法〔J〕。中国电机工程学报，2012，32（16）：130136。
　　〔19〕王宏民，苏燕涛，杨新华。Nuttall窗的光伏系统间谐波分析方法〔J〕。哈尔滨理工大学学报，2015（4）：667075。
　　〔20〕黄知超，延红艳，杨升振，等。凯瑟窗改进谐波法测电力电容器介质损耗因数〔J〕。哈尔滨理工大学学报，2012（5）：131825。