DFT Nuttall窗加权谐波分析算法.doc

第30卷 第9期 仪 器 仪 表 学 报 Vol.30 No. 92009年9月 Chinese Journal of Scientific Instrument Sep. 2009Nuttall窗加权谐波分析算法及其在电能计量中的应用*温和，滕召胜，李聪聪，贺静丹，魏双双（湖南大学电气与信息工程学院长沙410082）摘要：采用FFT进行电网信号谐波参数分析时存在频谱泄漏，准确度较低。研究提出了一种基于Nuttall窗函数的频谱相位差校正算法，分析了Nuttall窗函数的时域、频域特性，建立了频谱相位差校正算法的流程，推导了信号谐波参数计算式，给出了算法在三相谐波电能计量中的实现方法。采用Nuttall窗对电压、电流信号进行加权，利用Nuttall窗旁瓣电平低、旁瓣衰减快的特点可有效减少频谱泄漏；运用频谱相位差校正方法进行信号频谱校正与分析，不必求解高次方程即可实现电网信号谐波参数分析。仿真与实际运行结果表明，非同步采样情况下，本文算法可以准确地实现电网信号谐波分析，实现电能的准确计量。关键词：Nuttall窗；谐波分析；频谱泄漏；相位差；电能计量中图分类号：TN911TM714文献标识码：A国家标准学科分类代码：470.4037Harmonic analysis algorithm based on Nuttall window andits application in power measurementWen He, Teng Zhaosheng, Li Congcong, He Jingdan, Wei Shuangshuang（College of Electrical and Information Engineering, Hunan University, Changsha 410082, China）Abstract：The precision of harmonic analysis is affected by spectral leakage when using FFT in electrical signals. A new phase difference revision algorithm of discrete spectrum for Nuttall window weighted signal is proposed in this paper. The characteristics of Nuttall window are studied, and the workflow of the new algorithm is given. The application of the algorithm in power measurement is introduced. Nuttall window has optimal side lobe levels and the side lobes of the Nuttall window decay at a fast rate. Hence, spectral leakage can be reduced when the voltage and current signals are weighted by Nuttall window. The phase difference revision algorithm has the advantages of computational efficiency and easy implementation for Nuttall window weighted signal. Simulation results and power measurement application show the effectiveness and practicability of the algorithm.Key words：Nuttall window; harmonic analysis; spectral leakage; phase difference; power measurement第9期 温和 等：Nuttall窗加权谐波分析算法及其在电能计量中的应用 18251引言电力系统谐波的准确检测为电能计量与电能质量分析、故障诊断、电力谐波抑制与补偿等提供科学依据1-3。快速傅里叶变换（FFT）易于在DSP等嵌入式系统中实现，应用广泛。由于非同步采样和数据截短的影响，采用FFT进行谐波分析会产生频谱泄漏和栅栏效应，降低了电力系统谐波分析的准确度4。收稿日期：2008-08Received Date：2008-08*基金项目：国家自然科学基金（60872128）资助项目电力系统基波频率往往存在波动，即使采用了频率跟踪技术，也难以做到严格的同步采样5。为减少频谱泄漏和栅栏效应的影响，国内外学者研究了加窗插值FFT信号分析方法6-7，分别采用矩形窗5、Hanning窗6、Blackman-Harris窗7-8、Rife-Vincent窗9和矩形卷积窗10等实现插值FFT信号谐波分析，提高了准确度。但采用插值FFT方法进行信号谐波参数分析时，一般需要求解高次方程（一元五次或七次方程），特别是采用高阶余弦组合窗时，计算量较大11-12。为避免复杂运算，丁康等提出了一种基于对称窗的离散频谱相位差校正法13，该方法不需求解高次方程。但是，由于经典窗函数的旁瓣特性不够理想，采用经典窗函数对电压、电流信号进行加权处理后仍存在较大的频谱泄漏14-15，频谱相位差校正法的准确度仍然受到窗函数特性的限制16。为满足嵌入式系统实现和实际测量的要求，在分析了Nuttall窗旁瓣特性的基础上，本文提出一种采用Nuttall窗加权的信号频谱相位差校正算法，给出了其在电能计量中的应用，分析了Nuttall窗函数的时域、频域特性，推导了信号基波及各次谐波分量的参数计算公式。仿真结果表明本文算法能准确地实现基波与谐波参数测量。Nuttall窗加权频谱相位差校正算法在三相多功能谐波电能表中的应用验证了其有效性和准确性。2Nuttall窗函数Nuttall窗是一种余弦组合窗17，其时域表达式为： (1)式中，窗函数的项数为M，且n=1，2，N1。结合傅里叶变换，可得Nuttall窗的频谱函数的表达式为： (2)式中：WR(w)为矩形窗的频谱函数，令w=2k/N，得到Nuttall窗的离散频谱函数为： (3)表1Nuttall窗系数Table 1 Parameters of Nuttall windowbm三项最小旁瓣四项一阶四项三阶四项最小旁瓣b00.424 3800.355 7680.338 9460.363 581 9b10.497 340 60.487 3960.481 9730.489 177 5b20.078 279 30.144 2320.161 0540.136 599 5b30.012 6040.018 0270.010 641 1表1给出了几种典型的Nuttall窗函数系数。图1分别给出了Nuttall窗三项一阶（最小旁瓣）窗、四项一阶窗、四项一阶最小旁瓣窗、四项三阶窗的归一化幅频响应曲线。从图1可见，Nuttall窗四项窗旁瓣电平较低，但其主瓣宽度比三项窗要宽。表2列出了Nuttall窗与经典余弦组合窗的旁瓣特性参数。(a) (b) (c) (d) 图1Nuttall窗幅频响应曲线Fig.1 Frequency response of Nuttall window采用窗函数对信号进行加权时，窗函数的旁瓣电平大小和旁瓣衰减速率影响频谱泄漏量。选用旁瓣电平低且衰减速率大的窗函数能减少频谱泄漏，降低谐波间相互干扰，提高信号谐波参数分析的准确度。参见表2，四项三阶Nuttall窗具有较好的旁瓣特性，旁瓣衰减速率快（30 dB/oct），最大旁瓣电平为83 dB，因此，采用四项三阶Nuttall窗应具有较好的频谱泄漏抑制效果。但是，Nuttall窗的主瓣随着项数的增加而变宽，项数增加的同时也在一定程度上降低了频率分辨率。表2典型窗函数的旁瓣特性Table 2 Comparison of typical window characteristics窗的类型主瓣宽度旁瓣电平（dB）旁瓣衰减速率（dB/oct）Rectangle4/N136Triangle8/N2712Hamming8/N436Hanning窗8/N3218Blackman窗12/N5818Blackman-Harris窗16/N926三项最小旁瓣Nuttall窗16/N716四项一阶Nuttall窗16/N9318四项三阶Nuttall窗16/N8330四项最小旁瓣Nuttall窗16/N9863Nuttall窗加权频谱相位差校正方法采用Nuttall窗对信号进行加权后，可以减少频谱泄漏的影响，对FFT运算后得到的信号频谱进行相位差校正，可准确地求解信号参数。Nuttall窗加权信号频谱相位差校正算法首先在长度为的序列中分别选取1N的序列段1和1+LN+L的序列段2，分别对两段信号进行Nuttall窗加权运算。一般选取N为2的整数次幂（如512或1024），L为信号时域平移长度（0LN/2）。结合Nuttall窗的频谱函数，对长度为的序列（序列段1和序列段2）的离散频谱进行频谱相位差校正，获得信号基波和各次谐波的真实频率、幅值和初相角参数。设包含多项整数次谐波的时域信号表达式为： (4)式中：H为所含谐波的项数；h为整数，代表谐波的次数；f0为基波的频率；Ah和分别为第h项谐波的幅值和初相角。信号x(t)经采样频率为fs的数据采集系统后，得到离散序列： (5)信号x(n)(n=0，1，, N1)被截短为N点长序列段1和序列段2，并采用Nuttall窗函数分别对序列段1和序列段2进行加权处理。FFT运算后，序列段1的离散频谱为： (6)式中：Dkh=khk0(k=fN/fs，k0=f0N/fs)为离散频谱中第h次谐波真实频率处与离散谱线之间的偏差。第h次谐波的频率fh=hf0与第kh根谱线所指频率fh=khfs/N之间存在偏差： (7)利用第kh根谱线及相邻谱线，结合Nuttall窗谱函数，可求解第h次谐波参数。第kh根谱线处相位为： (8)设序列段2的离散频谱为，同理，根据离散傅里叶变换性质得到其第根谱线处的相位为： (9)序列段1和序列段2在第根谱线处形成的相位差为Dh=hh，由于相位差Dh一般在(2，2)的区间，需要对相位差按文献13所示方法进行调整，使Dh处于(，)区间，结合式（8）和式（9），调整后的Dh可写为： (10)参见式(10)，第h次谐波的频率偏差校正量为： (11)根据式(7)和式(11)，离散频谱中第h次谐波真实频率处与离散谱线之间的偏差为： (12)因此，第h次谐波的频率、初相角和幅值分别为： (13) (14) (15)设通过式(14)和式(15)计算后得到的第次谐波的电压、电流的幅值分别为Auh和Aih，电压、电流的初相角分别为和，则第h次谐波的瞬时有功功率、无功功率的计算式为： (16) (17)则第h次谐波的平均有功电能、无功电能为： (18) (19)第9期 温和 等：Nuttall窗加权谐波分析算法及其在电能计量中的应用 18274Nuttall窗加权频谱相位差校正方法4.1含基波和221次谐波的信号仿真实验设包含基波和221次谐波的信号时域表达式为： (20)式中：基波频率f0为50.3 Hz；采样频率为2 500 Hz；h为谐波次数；Ah和分别为第次谐波的幅值和初相角，具体取值见表3。表3基波及谐波参数Table 3 Parameters of the signal1234567(V)207.25.54.63.52.31.6()75604226326658891011121314(V)1.11.30.90.60.40.50.7()2784654662157215161718192021(V)0.850.640.510.390.40.20.1()43231867351948表4给出了选取长度N=512的四项三阶Nuttall窗进行仿真的结果。采用本文算法对微弱幅值频率分量（如第21次谐波）进行检测时，幅值相对误差为0.0002%，相角相对误差为-0.012 1%；对基波和221次谐波检测的幅值相对误差均0.001%，相角相对误差均0.1%。仿真结果表明：非同步采样时，四项三阶Nuttall窗能有效克服频谱泄漏影响，实现对基波和221次谐波分量的幅值、初相角的准确检测。表4Nuttall窗进行谐波分析的仿真结果Table 4 Simulation results with Nuttall window谐波次数1234567幅值相对误差(%)0.000 20.000 6-0.000 3-0.000 9-0.000 30.000 20.000 1初相角相对误差(%)0.002 10.021 80.024 80.033 3-0.004 90.002 10.013 9谐波次数891011121314幅值相对误差(%)-0.000 90.001 40.000 30.000 30.000 1-0.001 10.001 3初相角相对误差(%)0.054 0-0.012 10.003 30.010 80.006 20.079 00.010 5谐波次数15161718192021幅值相对误差(%)0.000 40.000 3-0.000 10.001 80.000 10.000 30.000 2初相角相对误差(%)0.010 00.036 20.028 7-0.020 80.001 60.019 9-0.012 14.2采用不同窗函数时的仿真实验利用式(20)所示信号，分别采用Hanning窗，Blackman窗、Blackman-Harris窗和四项三阶Nuttall窗对信号进行加权，并运用频谱相位差校正算法进行参数计算。图2给出了采用不同窗函数进行仿真得到的基波和221次谐波幅值绝对误差曲线。图2不同窗函数进行仿真的绝对误差Fig.2 Simulation results for different windows由图2可见，分别采用Hanning窗，Blackman窗、Blackman-Harris窗和四项三阶Nuttall窗进行信号谐波幅值检测时，采用四项三阶Nuttall窗的准确度最高，其次为Blackman-Harris窗和Blackman窗。如图2所示，采用四项三阶Nuttall窗得到的幅值仿真检测结果较稳定，幅值相对误差0.001%，存在较大波动。仿真结果表明：Nuttall窗具有良好的旁瓣性能，可有效抑制频谱泄漏，因此，与Hanning窗、Blackman窗、Blackman-Harris窗等经典窗函数相比，基于Nuttall窗函数的频谱相位差校正算法能有效提高谐波参数检测的准确度。4.3白噪声影响时的仿真实验对式(20)所示的信号添加高斯白噪声序列，高斯白噪声的信噪比为60 dB，其余参数不变，选择窗函数长度为512的Nuttall窗和Blackman-Harris窗进行仿真，得到的幅值、初相角相对误差曲线如图3所示。图3白噪声存在时不同窗函数进行仿真的相对误差Fig.3 Simulation results with white noise fordifferent windows参见图3，白噪声对弱幅值分量的影响较大，如采用Blackman-Harris窗时，第1721次谐波分量的幅值检测相对误差在0.1%，1%区间内，相角检测的相对误差约为1%；而采用四项三阶Nuttall窗时，第1721次谐波分量的幅值检测相对误差在0.01%，0.1%区间内，相角检测的相对误差约为0.1%。由图3可见，白噪声存在时，采用四项三阶Nuttall窗进行信号参数分析时，基波和各次谐波分量的幅值、相角检测的准确度均比采用Blackman-Harris窗高12个数量级。5应用与讨论基于Nuttall窗的频谱相位差校正算法已在湖南大学与湖南海兴电器公司合作研制的三相多功能谐波电能表中得到应用。如图4所示，三相多功能谐波电能表采用ADC+DSP+MCU结构形式，ADC选用ADI公司的AD73360，DSP选用ADI公司的ADSP-BF533，MCU选用RENESAS公司的M30624FGPFP。图4三相多功能谐波电能表结构框图Fig.4 Structure sketch of 3 phase power meterAD73360L对三相电网电压、电流信号进行实时采样，ADSP-BF533采用本文算法实现电流、电压的基波和221次谐波参数检测、基波与谐波电能计量，M30624FGPFP实现显示、数据统计、存储、通信以及功能选择等管理功能。三相多功能谐波电能表中，ADC采样频率设置为4kHz，选择长度的四项三阶Nuttall窗，表5给出了第21次谐波参数测量结果。表5第21次谐波参数测量结果Table 5 Measurement results of the 21st harmonicUA(%)UB(%)UC(%)IA(%)IB(%)IC(%)标准表7.027.087.116.757.076.96实测值6.876.696.836.576.416.35误 差0.460.720.593.123.233.15由表5可知，第21次谐波电压幅值的测量误差1%，第21次谐波电流幅值的测量误差3.5%，谐波电压、电流初相角误差5，达到GB/T-14549-93的A类谐波测量仪器标准（谐波电压幅值的测量误差2%，谐波电流幅值的测量误差5%）。仪器于2007年12月通过湖南省电力试验研究院的测试，挂网运行半年来稳定可靠。6结论为提高信号谐波分析准确度，本文提出一种基于Nuttall窗函数的频谱相位差校正算法，并将其应用于电能计量系统中。采用Nuttall窗对信号进行加权处理，并对离散频谱进行相位差校正。在分析了典型Nuttall窗的系数、旁瓣特性的基础上，本文给出了信号基波及各次谐波的频率、幅值和初相角计算公式。本文采用包含基波和221次谐波的信号进行仿真实验，分别选用Hanning窗、Blackman窗、Blackman-Harris窗和四项三阶Nuttall窗函数进行仿真分析，结果表明：与经典窗函数相比，非同步采样时，本文算法能有效抑制谐波间相互干扰和白噪声影响，准确度更高。三相多功能谐波电能表的研制与检验结果证明了Nuttall窗信号加权频谱相位差校正算法的有效性和准确性。参考文献1 COURSES E, SURVEYS T. Dynamics of frequency estimation in the frequency domainJ. IEEE Trans. on Instrumentation and Measurement, 2007,56(6):2111-2118. 2 DOUKOPOULOS X G, LEGOUABLE R. Robust channel estimation via FFT interpolation for multicarrier systemsC. Vehicular Technology Conference, IEEE 65th, 2007:1861-1865. 3 ANTONI J, SCHOUKENS J. A comprehensive study of the bias and variance of frequency-response-function measurements: optimal window selection and overlapping strategiesJ. Automatica, 2007,43(10):1723-1736. 4 DAVID M M, ALIREZA K Z, THOMAS H O. A new technique of measurement of nonstationary harmonicsJ. IEEE Trans. on Power Delivery, 2007,22(1):387-395. 5 JAIN V K, COLLINS W L, DAVIS D C. High-accuracy analog measurements via interpolated FFTJ. IEEE Trans. on Instrumentation and Measurement, 1979,28(2): 113-122. 6 潘文,钱俞寿,周鹗.基于加窗插值FFT的电力谐波测量理论(I)窗函数研究J.电工技术学报,1994,(1):50-54. PAN W, QIAN Y S, ZHOU E. Power harmonics measurement based on windows and interpolated FFT(I) study of windowsJ. Transactions of China Electro Technical Society, 1994,2(1):50-54.7 HARRIS F J. On the use of windows for harmonic analysis with the discrete Fourier transformJ. Proceedings of the IEEE, 1978,66(1):51-83. 8 AGRE D. Interpolation in the frequency domain to improve phase measurementJ. Measurement, 2008, 41(2): 151-159. 9 RIFE 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, 49(2): 197-228. 10 张介秋,陈砚圃,梁昌洪,等.交流电电气参量高精度测量的加权算法J.仪器仪表学报,2005,26(2): 172-176. ZHANG J Q, CHEN Y P, LIANG C H, et al. Weighted algorithm for high-accuracy measurement of AC electrical parametersJ. Chinese Journal of Scientific Instrument, 2005,26(2):172-176.11 庞浩,李东霞,俎云霄,等.应用FFT进行电力系统谐波分析的改进算法J.中国电机工程学报,2003,23(6):50-54. PANG H, LI D X, ZU Y X, et al. An improved algorithm for harmonic analysis of power system using FFT techniqueJ. Proceedings of the CSEE, 2003,23(6): 0-54. 12 AGRE D. Power measurement in the non-coherent samplingJ. Measurement, 2008,41(3):230-23513 谢明,张晓飞,丁康.频谱分析中用于相位和频率校正的相位差校正法J.振动工程学报,1999,12(4):454-459.XIE M, ZHANG X F, DING K. A phase difference correction method for phase and frequency correction in spectral analysisJ. Journal of Vibration Engineering, 1999,12(4):454-459. 14 沈国峰,王祁,张世平,等.电力交流参数测量中非均匀采样问题的误差分析J.仪器仪表学报,2003,24(4): 512-514.SHEN G F, WANG Q, ZHANG S, et al. Error analysis on non-uniformly sampling of power system AC parameter measurementJ. Chinese Journal of Scientific Instrument, 2003,24(4):512-514.15 戴