微机保护算法仿真示例.pptVIP

算法仿真示例 C4example1：数字滤波器设计 超高速方向继电器新算法的研究(II)——实现与仿真（参阅第八章：综合仿真） (以下引自：电力系统自动化，99/10，p23) 数字滤波器的设计 基于工频电气量原理的继电保护算法，其严重的不足之处在于受谐波分量的影响。数字滤波器的性能好坏直接影响保护装置的性能。 选择以下数字滤波器： H1 (z) =（1 – Z-m/4) 2 H2 (z) =[（1 – Z–7) / （1 – Z–1) ] 2 m为每周期的采样点数，取 m =96； H1(z)是有限种激响应数字滤波器，为有梳状频谱特性的简单数字滤波器，主要用来初步消去整次谐波及直流分量的影响； H2(z)为无限冲激响应数字滤波器，采用递归形式实现，是低通数字滤波器，能较好地消去高频分量的影响。(?！) 分析： H1 (z) =（1 – Z-m/4) 2 结构：相减滤波单元的级联; 阶次 : k = 96/4=24 作用：使用相减滤波单元消除直流分量和 4， 8，12，…l?4 次谐波分量 类型：有限冲激响应数字滤波器 （FIR滤波器） 分析： H2 (z) =[（1 – Z–7) / （1 – Z–1) ] 2 结构：积分滤波单元的级联; 阶次 : k = 8 作用：使用积分滤波单元增加对高次谐波 的衰耗 类型：有限冲激响应数字滤波器 （FIR滤波器；论文误为IIR滤波器） N=96;f0=50;fs=N*f0; a1=1;b1=[1 zeros(1,23) -1]; %设置差分滤波单元的传递函数系数 %k=N/m=24; (谐波阶次 m=4) % 文献给出的积分滤波单元的传递函数 a2=1;b2=[ones(1,9)]; %设置积分滤波单元的传递函数系数 k=8, %另外注意：设置系数b时，b0=1; f=0:1:fs/2; h1=abs(freqz(b1,a1,f,fs)); H1=h1/max(h1); %由传递函数系数确定传递函数的幅频特性 H11=H1.^2; %级联差分滤波器 h2=abs(freqz(b2,a2,f,fs)); H2=h2/max(h2); H22=H2.^2; %级联积分滤波器 h=h1.*h2;H=h./max(h); %计算级联滤波器传递函数的幅频特性 HH=H.^2; Figure %级联差分、积分滤波器幅频特性 subplot(221);plot(f,H11); xlabel(f/Hz);ylabel(H1); subplot(223);plot(f,H22); xlabel(f/Hz);ylabel(H2); subplot(222);plot(f,HH); xlabel(f/Hz);ylabel(H); subplot(224);plot(f,20*log10(HH+10e-10)); axis([0,fs/2,-60,0]);xlabel(f/Hz));ylabel(H/dB); 分析： H2 (z) =[（1 – Z–7) / （1 – Z–1) ] 2 并未针对特定谐波滤波进行设置，最终结 果非常不理想：实际上二次谐波有最大输 出；且3、5、6、7等各次滤波也有很大输 出。 结论： 在N=96时，论文给出的滤波器不能作为工频 选频滤波器。 改进： 考虑积分滤波单元1消除l?3次谐波的影响: m=3 阶次 : k=N/m–1=31 考虑积分滤波单元2消除l?2 次谐波的影响: m=2 阶次 : k=N/m – 1=47 H2 (z) =（1 – Z–31) （1 – Z–47) / （1 – Z–1) 2 仿真m文件作以下修改： N=96;f0=50;fs=N*f0; a1=1;b1=[1 zeros(1,23) -1]; %设置差分滤波单元的传递函数系数 k=24 a2=1;b2=[ones(1,32)]; %设置积分滤波单元的传递函数系数 k=31 a3=1;b3=[ones(1,48)]; %设置积分滤波单元的传递函数系数 k=47 重