电力系统暂态信号分析.docVIP

FFT 在电力系统谐波分析中的应用与研究 王文文 中国矿业大学信息与电气学院，江苏徐州（221008） E-mail:wangww714@ 摘 要：本文首先简单介绍了电力系统谐波的基本概念，利用FFT 结果对电力系统谐波参 数进行分析与研究，主要包括谐波系数的计算和有效值的计算，对电力系统谐波进行分析主 要分为三个步骤：必须对信号进行采样，变换为离散序列；建立数据窗，忽略数据窗前后的 信号波形;用谐波分析方法对数据窗进行分析。最后使用MATLAB 对其仿真，并对FFT 分 析方法和系统分析方法得到的结果进行比较。 关键词：谐波；FFT；MATLAB 中图分类号：TN919.5 1 概述 什么是电力系统谐波？国际上公认的谐波含义是：“谐波是一个周期电气量的正弦波的 分量，其频率为基波频率的整数倍。”在中华人民共和国国家标准中，谐波的定义是：“谐 波是对周期性的交流量进行傅立叶级数分解，得到的频率大于1 的整数倍分量”。从定义中 可以看到： (1)谐波次数必须是个正整数。例如我国电力系统的额定频率是50Hz，则n 次谐波的频 率为n×50Hz;有些国家电力系统的额定频率为60Hz，则其基波频率为60Hz, n 次谐波的频 率为n×60Hz 。 （2)间谐波、次谐波和分数谐波不属于以上定义的谐波范围。 傅里叶级数是一种研究和分析谐波畸变的有效方法。通过傅立叶分解能够对畸变波形的 各种分量进行检查【1】【2】。 2 利用FFT 进行谐波分析 2.1 同步采样 对电力系统谐波进行分析主要分为三个步骤：必须对信号进行采样，变换为离散序列； 建立数据窗，忽略数据窗前后的信号波形;最后用谐波分析方法对数据窗进行分析。 在进行以上三个步骤时，必须满足以下要求:首先满足采样定理，以免引起混叠，这一 点容易做到；其次采样频率必须与信号频率同步。也就是说，当数据窗为l 倍(l 为整数)信号 周期T、数据窗内采样次数为N (N 为整数)、采样频率为 s T 1 时,则当 N T lT s = （1） 时成立，即N 个采样点恰好采集了整数个信号周期时，称采样频率与信号频率同步， 此时称为同步采样。同步采样法又被称为等间隔整周期采样或等周期均匀采样。但是在实际 测量过程中，采样频率和信号频率由于受到很多因素的影响，会随时变化，严格意义上的同 步采样是难以实现的，因为： (1) 信号周期T 的测量不会绝对准确，并且信号周期有时还是变化的； (2)定时器的定时不可能绝对准确，其精度至少要受到最小定时单位(分辨率)的限制； (3) A/D 转换、程序执行等都会产生时间延迟。 -2- 中国科技论文在线 而且还不能采取简单的措施来使得采样频率和信号频率同步，即在采样的过程中必定有 同步误差的存在。正因为有同步误差的存在，给谐波分析带来或多或少的误差，所以如何减 小同步误差对谐波分析的影响一直是一个广为关注的课题[3][4] 。 2.2 FFT 谐波分析方法及参数计算方法 随着计算机和微电子技术的发展，开始采用基于快速傅里叶变换（FFT）的方法来检测 电力谐波。快速傅里叶变换时离散傅里叶变换（DFT）的一种快速算法，能够满足利用计算 机对有限长序列进行频谱分析的要求。 设序列 x(n)的长度为N，其DFT 为： Σ? = = 1 0 ( ) ( ) N n nk N X k x n W k=0,1,…,N-1 （2） 反变换（IDFT）为： Σ? = = ? 1 0 ( ) 1 ( ) N k nk N X k W N x n n=0,1,…,N-1 （3） 其中N j n N W e ? 2 = 。直接计算DFT，乘法次数和加法次数都和N 2成正比，当N 较大时， 运算量太大。快速傅里叶变换利用系数nk N W 的对称性，周期性和可约性的特性，采用蝶形 结构来完成离散傅里叶变换，起结构如图1 所示 K N w 1 ( ) X k ( ) 2 X k 1( ) 2( ) k X w k X kN + ( ) ( ) 1 2 k X w k X kN ? 图 1 信号蝶形流图 当 N = 2L 时，共有L 级蝶形，每级都有2 N 个蝶形运算组成，每对蝶形运算之间的距 离为2m?1。图中，m表示第m级；k表示第k行；把k值表示成L位的二进制数，左移L-m 位，将右边空出的位置补零，即为r 的二进制数。快速傅里叶变换的运算量和( 2 N ) log 2N . 成正比，运算量大大减少。 电流信号经变流器后，在一个周期内对其做等间隔的同步采样。对所获得的N 点时域 离散数据进行快速傅里叶变换，得到信号的频域数据，即可算出基波和各次谐波的幅值，相 位以及谐波电流含量，总畸变率等参数。 2.2.1.谐波系数的计算 通过 FFT运算，假设我们已经求得了被测信号x(n)高达127 次的