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

1、凯塞窗四谱线插值FFT的电力谐波分析方法摘 要：针对在采用快速傅里叶变换分析电力谐波时，对信号非同步采样以及非整周期截断的情况下，会产生频谱泄露和栅栏效应，无法有效检测谐波的各参数值等问题，提出了加窗插值算法，可以进步计算精度抑制频谱泄露，Kaiser窗的频带内能量主要集中在主瓣中，拥有绝佳的旁瓣抑制性能，该算法充分利用峰值谱线频点最近的四条谱线进展加权运算进步运算精度，运用MATLAB多项式拟合推导出实用的插值修正公式。仿真结果说明，四谱线插值算法得到的21次谐波的幅值相对误差0.00001%，相位相对误差0.000001%，基波频率计算相对误差为6.2×10-11%，并且可以有效

2、消除频谱泄露和栅栏效应的影响，计算精度更高。关键词：谐波分析；Kaiser窗；四谱线插值；FFTDOI：10.15938/j.jhust.2021.02.010中图分类号： TM935文献标志码： A文章编号： 1007-2683202102-0053-05Abstract：When the Fast Fourier Transform FFTmethod is used to analyze the harmonics in a power network， the fast Fourier transform FFTcan produce spectral leakage and fenc

3、e effect under the situation of non-synchronized sampling and non-integral period truncation， which leads to a great error in the results of calculation. The window functions and interpolation algorithm can reduce spectral leakage and improve the accuracy of computation. By MATLAB data fitting， prac

4、tical correction formulas corresponding to commonly used window functions are obtained by using weighted arithmetic of four pillars spectral line near the harmonic frequency point. The simulation and test results show that the relative errors of amplitude for 21 order harmonics are no more than 0.00

5、001%， the relative errors of phrase collection for 21 order harmonics are no more than 0.000001%， and the relative errors of fundamental wave frequency are 6.2×10-11%. The algorithm can eliminate the influence of spectral leakage and fence effect effectively， improving the accuracy of harmonic

6、analysis.Keywords：harmonic analysis； Kaiser window； four-spectrum-line interpolation； FFT0 引 言在供配电网络中，存在着含有大量的非线性、冲击性负荷的电力电子设备，导致电网谐波污染越来越严重。为了维护电网的平安稳定运行，获得实时的电网质量参数，对谐波参数进展高精度的检测显得尤其重要。常用的谐波分析方法是快速傅里叶变换FFT，在非同步采样时会产生频谱泄露和栅栏效应，影响谐波分析的精度。采用性能良好的窗函数可以有效抑制频谱泄露，对频域信号进展插值运算可减少栅栏效应引起的误差1-2。常用的窗函数有三角窗、Han

7、ning窗3、Blackman-Harris窗4、Rife-Vincent窗5、Nuttall窗6-7等。在插值算法中，常用的有纹紫卟逯邓惴8-9，双谱线插值算法10-14，三次样条插值算法15-16，三谱线插值算法等高精度插值算法17-19。传统的窗函数无法随意改变主旁瓣能量之间的比重，而Kaiser窗较为灵敏20。本文对Kaiser窗的频谱特性进展分析，在双谱线插值和三谱线插值的根底上，提出了基于Kaiser窗四谱线插值FFT的电力谐波分析方法，同时为了计算的方便快捷，运用MATLAB多项式拟合推导出实用的插值修正公式，并对本文提出的算法进展了验证。1 Kaiser窗特性Kaiser窗由零

8、阶贝塞尔函数构成，可自由选择主瓣宽度和旁瓣高度之间的比重，频带内能量主要集中在主瓣中，其时域表达式为：图1a、b给出了=0，4，8，11时Kaiser窗的时域和频域特性。由时域图形可见，当=0时，Kaiser窗变为矩形窗函数，其旁瓣渐进衰减速率为6dB/oct，当=8时，旁瓣峰值为-58.3dB，旁瓣渐进衰减速率为12dB/oct，随着值的增大，旁瓣峰值继续下降，旁瓣渐进衰减速率增加。本文算法选择=20。2 基于Kaiser窗的四谱线插值算法3 仿真验证及结果分析3.1 简单信号分析 为验证基于Kaiser窗的四谱线插值算法的准确性，采用已被众多文献验证的无噪声信号模型，其表达式为：本文选取H

9、anning窗、Kaiser窗双谱线、三谱线和四谱线=20，N=1024进展仿真。从表1数据可以看出，相比于Hanning窗、Kaiser双谱线、Kaiser三谱线，本文Kaiser窗四谱线插值算法具有良好的抑制频谱泄露的作用，并且精度更高。3.2 利用复杂信号对算法进展的仿真验证为了便于比较分析，本文采用了已被众多文献验证的包含21次谐波的复杂信号模型：依次用Kaiser窗双谱线插值算法、三谱线插值算法和本文提出的四谱线插值算法对谐波信号进展处理双谱线插值采用文10的修正公式及结果，三谱线插值采用文17提出的插值算法 ，仿真结果如表3、表4所示。其中表3为基波与各次谐波幅值的相对误差，最后一

10、项为基波频率的相对误差；表4为基波与各次谐波相位的相对误差。从表3、表4仿真结果可以看出，在其他条件一样且=20的情况下，本文四谱线插值算法得到的21次谐波的幅值相对误差0.00001%，相位相对误差0.000001%，由于本文算法流程，只求一次基波频率，谐波频率根据基波频率的倍数求得，所以基波频率的精度非常重要，由表3可得双谱线、三谱线和四谱线基波频率计算相对误差分别为-1.4×10-9%、4×10-10%、6.2×10-11%，四谱线插值算法比双谱线和三谱线准确度更高。可见本文提出的Kaiser窗四谱线插值算法与双谱线和三谱线相比，由于进步了修正公式的准确度，

11、从而显著的进步了谐波参数检测准确性。4 结 论由于Kaiser窗拥有绝佳的频谱泄露抑制效果，并且，四谱线插值较高的运算精度，本文提出了基于Kaiser窗四谱线插值FFT的电力谐波分析方法。并且为了计算方便，根据MATLAB多项式拟合求出了实用的插值修正公式。通过对简单信号和复杂信号在不同插值算法的仿真比照，结果说明本文所提算法在对采样数据进展非整周期截断情况下，与Kaiser窗的双谱线和三谱线插值算法相比，具有更高的计算准确度。参 考 文 献：1许珉， 杨阳， 陈飞， 等. 基于 NuttallI窗的插值 FFT算法J. 电力系统保护与控制， 2021， 3923： 44-47.2曾博， 唐求

12、， 卿柏元， 等， 基于 Nuttall 自卷积窗的改进FFT 谱分析方法J. 电工技术学报， 2021， 297： 59-65.3BARROS J， DIEGO R I. On the Use of the Hanning Window for Harmonic Analysis in the Standard Frame-workJ.IEEE Trans on Power Delivery， 2006， 211： 538-539.4HARRIS F J. On the Use of Windows for Harmonic Analysis with the Discrete Fourie

13、r Transform J.Proceedings of the IEEE， 1978， 661： 51-83.5RIFE D C，VINCENT G A.Use of the Discrete Fourier Transform in the Measurement of Frequencies and Levels of TonesJ.The Bell System Technical Journal，1970，492：197-228.6NUTTALL A H.Some Windows with a Very Good Side Lobe Behavior J.IEEE Trans. on

14、 Acoustics，Speech and Signal Processing， 1981， 291：84-91.7卿柏元，滕召胜，高云鹏，等.基于 Nuttall 窗双谱线插值 FFT 的电力谐波分析方法J.中国电机工程学报，2021，2825：153-158.8王楠，肖先勇，曾颂崎，等.基于插值线性调频 Z 变换的间谐波分析方法J.电网技术，2007，3118：43-47.9肖先勇，王楠，刘亚梅.基于多项式逼近的单峰谱线插值算法在间谐波分析中的应用J.电网技术，2021，3218：57-61.10高云鹏， 腾召胜， 卿柏元. 基于 Kaiser 窗双谱线插值FFT 的谐波分析方法J. 仪器

15、仪表学报， 2021， 312： 287-292.11高云鹏， 滕召胜， 温和， 等. 凯塞窗差值 FFT 的电力谐波分析与应用J. 中国电机工程学报， 2021， 304：43-48.12曾博， 滕邵胜， 高云鹏， 等. 基于 Rife-Vincent 窗的高嗜范鹊缌巢嗔考扑惴椒J. 电工技术报，2021， 248： 154-159.13庞浩，李东霞，俎云霄，等. 应用 FFT 进展电力系统谐波分析的改进算法J. 中国电机工程学报，2003，23 6： 50-54.14薛蕙，杨仁刚. 基于 FFT 的高精度谐波检测算法J.中国电机工程学报，2002，22 12： 106-110.15王继东， 杜旭浩， 杨帆. 基于三次样条插值信号重构的微网谐涉及间谐波分析算法J. 电网技术， 2021，3611： 7-11.16杨晓萍， 刘普森， 钟彦儒. 基于经历形式分解的有源滤波器谐波检测J. 电工技术学报， 2021， 245：197-202.17陈立伟，刘欢，冯驰，等.Kaiser 窗三谱线插值电力谐波分析J.应用科技，20212：5-11.18牛胜锁，