点容易做到；其次采样频率必须与信号频率同步。也就是说，当数据窗为 l 倍(l 为整数)信号
1
周期 T、数据窗内采样次数为 N (N 为整数)、采样频率为 时,则当
Ts
lT = N Ts
（1）
时成立，即 N 个采样点恰好采集了整数个信号周期时，称采样频率与信号频率同步， 此时称为同步采样。同步采样法又被称为等间隔整周期采样或等周期均匀采样。但是在实际
表 1 摘取了部分谐波系数来说明 FFT 分析方法的准确性以及和系统分析方法的比较。
谐波次数
0 1 5 10 15 20 25 30
频率（HZ）
0 50 250 500 750 1000 1250 1500
表 1 部分谐波系数分析表
谐波系数
谐波系数
FFT 分析
0.0016
0.0016
0.9600
0.9600
t
图 4 FFT 结果模值
从图 3 中可以看出由于含有高次谐波，所以波形不平滑，在实际系统中这种信号时不允 许出现的，这里只是用它来验证结果是否合理。
我们用 MATLAB 对输入信号进行 64 点离散采样后，得到了 64 个样本点，我们将以 64 点样本点为采样数据，分别用 MATLAB 和系统分析对其进行 FFT 变换从而来实现谐波分析。
从结果可知，经 FFT 分析得到的谐波系数和真实值完全相等，所以 FFT 分析方法可用
于谐波的分析。
谐波有效值 1．1001
表 2 三种方法得到有效值比较 FFT 分析谐波有效值 1.1001
系统分析谐波有效值 1.1023
谐波有效值即：根据原图形信号的各次谐波的幅值的平方和开方后所得的真实有效值； 而 FFT 分析谐波有效值即：在原图形信号抽样 64 点后，对抽样信号做 FFT，应用 FFT 结果，根据公式（8）（9）可计算信号出有效值； 系统分析谐波有效值即：在原图形信号抽样 64 点后，不做 FFT，而是直接根据系统分 析谐波的公式（此处未列出），计算信号有效值。 可见，从利用 FFT 结果频谱分析出来的有效值和谐波真实的有效值相比较，是正确的， 而系统分析谐波则存在误差，所以 FFT 在系统谐波分析中得到广泛的应用。
Abstract This article describes the basic definition of Power System Harmonic(PSH) , then analyzes and researches on parameters of PSH consisting essentially of the calculation of harmonic coefficient and effective values through the Fast Fourier Transform (FFT) way.PSH is divided into three steps: the signal must be sampled, transformed into discrete sequence; establish the data window and ignore the signal wave after the data window; harmonic analysis method to analyze the data window though the harmonic analysis method . Finally, uses MATLAB software to simulate it. At last the paper compares the results which were got from the analysis of the FFT method and system analysis. Key words: harmonic;FF;MATLAB
可得电压与电流的有效值为：
U m=
127
∑ A2 Um
(k
)
k =1
2
（8）
127
∑ A2 Im
(k
)
I m=
k =1
2
（9）
其中，m=1，2，3；U m , I m 分别为被测电压电流u1(n), u2 (n),u3(n), i1(n), i2 (n),i3(n),
的有效值。由此可以滤出信号中的直流分量，而又含有基波至 127 次的有效值。
（2)间谐波、次谐波和分数谐波不属于以上定义的谐波范围。 傅里叶级数是一种研究和分析谐波畸变的有效方法。通过傅立叶分解能够对畸变波形的 各种分量进行检查【1】【2】。
2 利用 FFT 进行谐波分析
2.1 同步采样
对电力系统谐波进行分析主要分为三个步骤：必须对信号进行采样，变换为离散序列；
建立数据窗，忽略数据窗前后的信号波形;最后用谐波分析方法对数据窗进行分析。 在进行以上三个步骤时，必须满足以下要求:首先满足采样定理，以免引起混叠，这一
AX (k) = X (k) =
X
2 r
(k)
+
X
2 i
(k
)
k=0,1,2,…,127
（5）
在三相交流电采集系统中，对于每一条被测线路各有三个电压输入和三个电流输入：
u1(n), u2 (n),u3(n), i1(n), i2 (n),i3(n), 根据上式可得电压和电流的振幅：
•
AUm (k) = Um (k ) =
-4-
中国科技论文在线

参考文献
[1] 吴竞昌.供电系统谐波[M]. 北京:中国电力出版社，1998 [2」中华人民共和国国家标准.电能质量公用电网谐波[M]. GB/T 14549-9 [3] 陈文进，江道灼.电力系统数据采样的几种同步方法〔J〕.继电器，2003, 31(11):2s-29 [4] 毛筱，肖雁鸿，龚理专等.FFT 应用于谐波测量中频谱泄漏的分析和处理[J]电工技术杂志，2001, 2: 3-4 [5] 陈怀琛，吴大正，高西全，MATLAB 及在电子信息课程中的应用[第 2 版]，北京，电子工业出版社， 2003 [6] 薛年喜，MATLAB 在数字信号处理中的应用，清华大学出版社，2003
随着计算机和微电子技术的发展，开始采用基于快速傅里叶变换（FFT）的方法来检测 电力谐波。快速傅里叶变换时离散傅里叶变换（DFT）的一种快速算法，能够满足利用计算 机对有限长序列进行频谱分析的要求。
设序列 x(n)的长度为 N，其 DFT 为：
N −1
∑ X (k) = x(n)WNnk n=0
反变换（IDFT）为：
-1-
中国科技论文在线

而且还不能采取简单的措施来使得采样频率和信号频率同步，即在采样的过程中必定有 同步误差的存在。正因为有同步误差的存在，给谐波分析带来或多或少的误差，所以如何减 小同步误差对谐波分析的影响一直是一个广为关注的课题[3][4] 。
2.2 FFT 谐波分析方法及参数计算方法
0.0800
0.0800
0.0640
0.0640
0.0600
0.0600
0.0040
0.0040
0.0016
0.0016
0.0006
0.0006
谐波系数 系统分析
0.0010 0.9600 0.0802 0.0656 0.0638 0.0092 0.0061 0.0022
绝对误差
0.0006 0.0000 0.0002 0.0016 0.0038 0.0052 0.0045 0.0016
1 2
U
2 mr
(k
)
+
U
2 mi
(k
)
（6）
•
AIm (k ) = Im (k ) =
1 2
I
2 mr
(k
)
+
I
2 mi
(k
)
（7）
其中，m=1，2，3； AUm (k) , AIm (k ) 分别为三相交流电压和电流第 k 次的谐波分量的
幅值[5][6] 。
2.2.2.有效值的计算 非正弦周期电压或电流的有效值等于各次谐波电压电流有效值的平方和的开方，由此，

原始信号图像
1.5
1
0.5
采样后的信号图像 1.5
1
0.5
0 0
-0.5 -0.5
-1 -1
-1.5 0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 t
图 3 原始图形
-1.5 0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02
1 概述
什么是电力系统谐波？国际上公认的谐波含义是：“谐波是一个周期电气量的正弦波的 分量，其频率为基波频率的整数倍。”在中华人民共和国国家标准中，谐波的定义是：“谐 波是对周期性的交流量进行傅立叶级数分解，得到的频率大于 1 的整数倍分量”。从定义中 可以看到：
(1)谐波次数必须是个正整数。例如我国电力系统的额定频率是 50Hz，则 n 次谐波的频 率为 n×50Hz;有些国家电力系统的额定频率为 60Hz，则其基波频率为 60Hz, n 次谐波的频 率为 n×60Hz 。
中国科技论文在线

FFT 在电力系统谐波分析中的应用与研究
王文文
中国矿业大学信息与电气学院，江苏 徐州（221008）
E-mail:wangww714@
摘 要：本文首先简单介绍了电力系统谐波的基本概念，利用 FFT 结果对电力系统谐波参 数进行分析与研究，主要包括谐波系数的计算和有效值的计算，对电力系统谐波进行分析主 要分为三个步骤：必须对信号进行采样，变换为离散序列；建立数据窗，忽略数据窗前后的 信号波形;用谐波分析方法对数据窗进行分析。最后使用 MATLAB 对其仿真，并对 FFT 分 析方法和系统分析方法得到的结果进行比较。 关键词：谐波；FFT；MATLAB 中图分类号：TN919.5