微机保护的算法.pptVIP

三、基于信号为周期函数模型的算法前几种算法按输入信号为正弦函数来考虑，但电力系统的故障信号为非纯正弦，包含直流、基波、高次谐波分量。采用前述几种算法时，需要严格的预滤波。预滤波与算法做统一考虑时，前述算法所用的计算时间就较长。傅氏算法：基本思路来自傅里叶分解，本身具有滤波作用。傅氏算法依据：将傅氏级数与继电保护信号进行积分变换，求出保护所需基频分量或倍频分量实部、虚部的系数，进而求出其幅值和相位。1.周期函数的傅氏级数及各次谐波的关系一个周期函数i(t)可由傅里叶级数展开为：an、bn－各次谐波分量的正弦项、余弦项系数；?1－基波角频率；n－谐波次数。根据傅立叶级数性质与三角函数的正交性，有：将基波分量改写为：可见：基波分量的有效值和相角为：其它各次谐波分量的求法与此相同。用复数表示为:2.全周波傅氏算法将a1、b1离散化得（）式中，N—每工频周期采样点数；ik—第k次采样值。基波分量的模值(幅值)、相位角为：一般将a1、b1记作基波分量的实部IR1、虚部II1。基波分量的向量表达式为：有些微机保护中，需用到谐波分量。由全周傅氏算法，可得其n次倍频分量的正弦、余弦系数，离散形式：采样频率确定后，cosnk?/N、sinnk?/N可先离线计算出来。k123456789101112sink300.50.86610.8660.50-0.5-0.866-1-0.866-0.50cosk300.8660.50-0.5-0.866-1-0.866-0.500.50.8661全周傅氏算法系数表（fs=600Hz,N=12，求取基波）上表中，全周傅氏算法的系数仅有三组值：0.5、0.866、1，计算基波分量的正弦、余弦系数的表达式为:（fs=600Hz）为求出故障分量，i1~i12必须是故障后的采样值，即所需数据窗为一个周波（20ms）。3.半周波傅氏算法对半周波傅氏算法，其n次倍频分量的正弦、余弦系数，其离散形式：当fs=600Hz,n=1时：为求出故障分量，所需数据窗为半个周波（10ms）。4.一阶差分后半周波傅氏算法半周傅氏算法不能滤除恒定直流分量和偶次谐波分量，可先作一次差分运算后再应用半周傅氏算法，提高计算精度。当fs=600Hz,n=1时：一阶差分后，信号中的直流分量被完全滤除。5.傅氏算法的简化计算用离散值累加和代替连续积分，结果会产生误差。傅氏算法计算信号中的一个频率分量，需2N次乘、(2N-1)次加，计算量大，可采用递推简化算法。设信号i(t)每周采样次数N，Ts＝2?/N，则一个基波周期后的t=mTs(mN)采样时刻的正弦、余弦系数：此时只需２次乘、4次加减，与Ｎ无关，运算量减小。6.傅氏算法的滤波特性全周傅氏算法能滤除恒定直流分量和各整数倍次的谐波分量。半周傅氏算法不能滤除恒定直流和偶次谐波分量.两种滤波器不能完全滤掉按指数衰减的非周期分量、非整次高频分量，但对它们有一定的抑制能力。不同初相角情况下，傅氏算法对基波、直流、整次谐波分量的频率响应相同，但对非整次谐波分量的频率响应有较大区别。全周傅氏算法的幅频特性不同初相角时全周傅氏的滤波特性半周傅氏算法的幅频特性不同初相角时半周傅氏的滤波特性全周傅氏算法的Matlab仿真设输入信号为：x(t)＝sin(?t)+0.3sin(3?t)+0.1sin(5?t)时，每周采样点数N=60，利用傅氏算法计算输入信号的基波、三次、五次谐波幅值。Matlab的M文件如下：N=60;t=0:0.02/N:0.06;y=sin(2*pi*50*t)+0.3*sin(2*3*pi*50*t)+0.1*sin(2*5*pi*50*t);i=1:N;bs1=sin(i*2*pi/N);bc1=cos(i*2*pi/N);bs3=sin(i*2*3*pi/N);bc3=cos(i*2*3*pi/N);bs5=sin(i*2*5*pi/N);bc5=cos(i*2*5*pi/N);a=1ys1=filter(bs1,a,y);yc1=filter(bc1,a,y);ys3=filter(bs3,a,y);yc3=filter(bc3,a,y);ys5=filter(bs5,a,y);yc5=filter(bc5,a,y);ym1=2*sqrt(ys1.^2+yc1.^2)/N;ym3=2*sqrt(