2012年月电能质量第三讲.ppt

第三章 电能质量分析数学方法 四川大学电气信息学院 肖先勇 2012年3月 许多数学分析方法已在电能质量领域得到应用，主要用于特性分析、检测、分类、定位、评估等方面。 主要方法包括：时域、频域和基于变换的方法等三类： 时域法应用最广泛，用时域仿真程序分析电能质量扰动，如EMPT、EMTDC、NETOMAC、Matlab、PSPICE等； 频域法主要用于谐波分析，包括：频谱分布、谐波潮流计算等； 基于变换（数学）的方法主要有：傅立叶变换法、短时傅立叶变换法、矢量变换法、小波变换法、S变换法，以及人工神经网络法等。 一、 非正弦周期信号分解为傅里叶三角级数 任意周期性电压和电流等都可用一个周期函数表示： 其中：T——基本周期。 非正弦周期函数满足狄里赫利条件时可分解为傅立叶级数，而在电力系统中的光滑函数通常都能满足这个条件。 傅里叶的三角级数形式为： 也可写为： 式中：w1——周期函数的角频率，h——谐波次数。 给定离散时间序列x1,x2,……,xN-1，设该序列绝对可和，即满足 则 (n＝0，1，…,N-1 ) 被称为序列 的离散傅里叶变换。 (n＝0，1，…,N-1) 被称为序列 的逆离散傅里叶变换。 式中，n相当于对时间域的离散化，k相当于频率域的离散化，且它们都是以N点为周期的。而离散傅里叶序列 是以2p为周期的，且具有共轭对称性。 由式 可见，要求出N点需要次复数乘法，N(N-1)次复数加法。 众所周知，实现1次复数乘法需要4次实数相乘和2次实数相加，实现1次复数加法则需要2次实数相加。当N很大时，计算量相当可观。例如，若N=1024，则需要104.8576万次复数实数相乘。所需时间过长，难以“实时”计算——这是算法本身决定的。 由离散傅里叶变换式系数的共轭对称性，即 ，可见 ，幅频特性是与纵坐标轴对称的。由 的周期性可见： 即幅频特性为周期性的偶函数，如图： 当采样点数为N时，仅给出N/2个频谱分量的数值。例如，选取每周期128个采样点时，只能得到64个及以下的谐波幅值。 因此，可对采样定理作如下解释： 采样频率至少是原信号最高频率的2倍以上(fs2fc)，采样才能正确地表述原信号的信息。将最高频率的2倍频率2fc称为奈奎斯特频率。 当采样频率低于内奎斯特频率时，原信号中高于fs/2的频谱分量将会低于fs/2的频率中再现，即出现频谱混叠，使频谱分析出现误差，如图： 为防止频谱混叠，可先用带宽低于fs/2的低通滤波器滤除高于fs/2的分量，然后作离散傅立叶变换。这样原信号中低于fs/2的频率分量能得到准确刻画，但会失去高于fs/2的频率分量。 例如：对于方波信号，若不经低通滤波而对其采样作离散傅立叶变换，则会因频率混叠引入误差；若经低通滤波，如使其仅含7次以下的谐波分量，对其采样作16点以上的离散傅立叶变换频谱分析，不会出现混叠，但已先在方波中舍去了频率高于7次的谐波分量。 利用W因子的周期性和对称性导出高效的快速算法——快速傅立叶变换算法(FFT)。 FFT使N点DFT的乘法计算量由N2次降到 次。如：当N=1024时，计算量为5120次，仅为原来的4.88%。因此，人们公认，FFT是数字信号处理史上的一个转折点，也可称之为一个里程碑。 对于离散傅里叶逆变换（IDFT）式: 将以上思路推广到点N=2M的DFT的一般情况，不难看出，第m次分解的结果是由2m个N/2m 点的DFT两两组成共2m-1个N/2m-1点的DFT。由于N=2M，通过M=log2N次分解后，最终达到了N/2个两点DFT的运算，从而构成了由x(n)到X(k)的M级运算过程。其迭代过程如图： 傅里叶变换是从时域到频域的相互转换的工具。从物理意义上讲，其实质是把f(t) 波形分解成许多不同频率的正弦波的叠加和。这样就可以把对原函数f(t)的研究转化为其权系数即傅里