电力系统谐波的FFT加窗插值算法的研究.pdfVIP

中国电机工程学会第十届青年学术会议 ·吉林 电力系统谐波的 FFT 加窗插值算法研究 刘旭东 叶鹏 （沈阳工业大学电气工程学院 辽宁 沈阳 110178） 摘 要：采用快速傅里叶变换（FFT）进行电力系统谐波分析时很难做到同步采样，故造成频谱泄漏，影响谐波分析的 结果。本文对（FFT）的泄漏原因进行了分析，并研究组合余弦窗对采样数据加权及利用插值对（FFT）的结果进行修正。通 过对卷积窗的谐波理论分析与研究，提出了基于卷积窗的组合余弦窗的加窗算法。文中给出了该算法进行谐波分析的算例， 计算结果表明，该窗的算法具有更高的计算精度。 关键词： 谐波分析； FFT；窗函数 1.引言 近年来，随着直流输电和柔性交流输电技术的采用，电气化铁道的快速发展，化工、冶金、煤炭等工 业部门大量应用电力电子设备等，使得电网的谐波含量大大增加，电网波形畸变越来越严重，对电力系统 的安全、经济运行造成极大的影响。谐波测量是谐波问题研究的主要依据，实时测量电网中的谐波含量， 确切掌握电网中谐波的实际状况，对于防止谐波危害，维护电网的安全运行是十分必要的。因此如何能够 把谐波的危害最大限度地减少，是目前电力电子领域极为关注的问题，而解决这一问题的关键在于定量地 确定谐波的成分、幅值和相位等。这也正是我们进行谐波分析的目的所在。 电力系统的谐波分析，通常都是通过快速傅立叶变换(FFT)实现的。然而 FFT 存在栅栏效应和泄漏现 象，使算出的信号参数即频率、幅值和相位不准，无法满足准确的谐波测量要求。而解决这个问题通常是 通过加窗函数减少泄露现象和通过插值修正减小栅栏效应。对于包含整数倍次谐波的信号而言，如果能够 实现同步于基波的采样，从而使截断数据长度是基波周期的整数倍，就能够使被测的基波和谐波频率恰好 与傅立叶变换频域谱线的离散频点重合，而且此时信号中的基波或谐波在其它频率成分上造成的长范围泄 漏也将为零。然而，要实现同步采样和整倍数截断是比较困难的。无论是 Ignacio 给出的采用插值补点改 变信号抽样率的方法[1] [2] ，还是 Roberto 的利用过采样器进行补点的方法 ，在实现上都是相当复杂的。但 可以通过选择适当的窗函数抑制长范围泄漏，也可以根据所选择的窗函数的形式对频率、相位和幅值进行 插值修正，在一定程度上弥补短范围泄漏造成的误差。很多文献给出了加窗和插值修正算法的研究成果[3-6] 。 本文对（FFT）的泄漏原因进行了分析，并研究组合余弦窗对采样数据加权及利用插值对（FFT）的结果进 行修正[7] 。通过对卷积窗的谐波理论分析与研究，提出了基于卷积窗的组合余弦窗的加窗算法。文中给出 了该算法进行谐波分析的算例，计算结果表明，该窗的算法具有更高的计算精度。 2.算法原理 2.1 频谱泄露和栅栏效应 在实际的谐波测量中，所要处理的信号都是经过采样和 A/D 转换得到的有限长数字序列，这相当于对 原始信号乘以一个矩形窗加以截断。在时域相乘，相当于在频域卷积。因此原信号的单一频率将变成以原 信号频率Ω0 为中心，形状振荡并逐渐衰减的连续谱线 sin(f)/f，这就是说信号频谱泄露到整个频率轴上 了。当采样同步，即时间窗长度正好等于信号周期的整数倍时，离散傅立叶变换得到的各离散频率分量除 了Ω0 点以外都落在 sin(f)/f=0 上，因此，变换的结果不会产生频谱泄露。 在非同步采样时，由于实际信号的各次谐波分量不能正好落在频率分辨点上，而是落在某两个频率分 辨点之间。这样通过 DFT 并不能直接得到各次谐波分量的准确值，而只能以临近的频率分辨点的值来近似 代替。这就是通常所说的栅栏效应。插值算法可以较好地消除栅栏效应，而谐波间的泄露引起的误差则需 要依靠加窗的方法来修正。 加窗插值算法通过加窗减小频谱泄漏，通过插值消除栅栏效应引起的误差。算法中，窗函数的选择非 常重要。通常频谱分析要求窗函数主瓣窄、旁瓣低且跌落速