电力系统谐波分析的分块FFT+算法研究.pdfVIP

·620 · 2008 全国博士生学术论坛电气工程论文集 电力系统谐波分析的分块 FFT 算法 梅 永 1，2 王柏林 1 1 河海大学电气工程学院 江苏 南京 210098 2 南京信息工程大学电子与信息工程学院 江苏 南京 210044 【摘 要】 为了得到良好的精度，FFT （快速傅里叶变换）算法要求时域采样点数N 足够大，但 N 越 大计算量也越大。针对电力系统谐波分析中频域点数远小于时域点数的特点，提出了一种 适合于电力系统谐波分析的分块 FFT 算法，并研究了加窗分块 FFT 算法在电力系统谐波 分析中的应用。研究结果表明，分块 FFT 算法与标准 FFT 算法计算结果完全相同，但前 者的计算量和占用内存的数量都比标准 FFT 少得多。 【关键词】 谐波分析 加窗 FFT 分块 FFT 在谐波分析 中，最常用的算法是 FFT[1 ，2] 。为了提高精度，FFT 算法要求时域采样点数 N 足够 大，N 越大 FFT 的维数就越高。研究 电力系统的谐波分析不难发现：N 维 FFT 求出的 N 个频域系数大 部分是不必要的 有用的仅仅是其 中的几十个系数。 频 域 点 数 （ H ）远 小 于 时 域 点 数 （N ） 的 FFT 已 经 引 起 了 研 究 者 的 关 注 ， 文献 [3]论 证 了 当 H log2 N 时直接 DFT 比标准 FFT 计算量少。还有一个 ‘修剪算法 ’[4 ，5] ，它是在标准 FFT 的计算 ‘树’上修剪掉无用的分支 ，从而减少计算量，不过对于数千点的 FFT ，这种修剪 的规则比较复杂 。本 文提出一种适合于电力系统谐波分析 的分块 FFT 算法， 它规则简单 ，计算结果与标准 FFT 完全相同， 而计算量和 占用内存 的数量都比标准 FFT 少得多。 1 DFT 的矩阵表示 可以证明，N 维 DFT 可以表示成矩 阵形式[6] 1 X = Fy （1） N 式中 ⎡ 1 1 1 1 ⎤ ⎢ 2π 2π 2π ⎥ ⎢ 1 e− j N e−j N 2 e− j N (N −1) ⎥ ⎢ ⎥ F = ⎢ − j 2π 2 −j 2π 4 −j 2π 2(N −1) ⎥ （2 ） ⎢ 1 e N e N e N ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ − 2π ( −1) − 2π 2(