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

ys=filter(bs,a,y); yc=filter(bc,a,y); ym=2*abs(complex(ys,yc))./N; %付立叶算法计算原始信号幅值 y1=y(1,1:N); %截取一周波原始信号 I=sum(y1)/N; tt1 = T* I/(y(1)-y(N+1)); A=(y(1)-y(N+1))/(1-exp(-T/tt1)); y0=A*exp(-t./tt1); %计算衰减分量 y1=y-y0; %对原始信号修正 figure plot(t,y,t,y1,t,y2,t,exp(-t/t1)); legend(‘y(原始信号)’,‘y0’,‘y-y1(修正信号), exp(-t/t1)) ys2=filter(bs,a,y2); yc2=filter(bc,a,y2); ym2=2*abs(complex(ys2,yc2))./N; %付立叶算法计算修正信号幅值 figure plot(t,ym,t,ym2); C4example4：滤序算法 f=50;N=24;t=0:0.02/N:0.03; va=sin(2*pi*f*t); [x,y]=size(va); % y:序列长度 vb=sin(2*pi*f*t-2*pi/3); vc=sin(2*pi*f*t+2*pi/3); ％ 输入正序电压 va1=[0,0,0,0, va(1:y-4)]; ％ 延迟NTs/6 vb1=[0,0,0,0, vb(1:y-4)]; vc1=[0,0,0,0, vc(1:y-4)]; v1=[(vb-va)-(vc1-va1)]/3; v2=[(vc-va)-(vb1-va1)]/3; figure subplot(221) plot(t,va,t,vb,t,vc) legend(va,vb,vc) title( 三相正序电压) subplot(223) plot(t,v1,-o,t,v2,-*) legend(正序电压v1,负序电压v2) ％类似写出输入负序电压时的m语句 滤序算法仿真 C4example4：任意角移相算法 f=50;N=24;t=0:0.02/N:0.04; v1=sin(2*pi*f*t); [x,y]=size(v1); v2=[0,v1(1:y-1)]; p=3.7319;q=3.8637; v0=p*v1-q*v2; figure plot(t,v1,t,v0) legend(‘v1’,‘v0’)；grid 任意角移相算法仿真 人有了知识，就会具备各种分析能力， 明辨是非的能力。 所以我们要勤恳读书，广泛阅读， 古人说“书中自有黄金屋。 ”通过阅读科技书籍，我们能丰富知识， 培养逻辑思维能力； 通过阅读文学作品，我们能提高文学鉴赏水平， 培养文学情趣； 通过阅读报刊，我们能增长见识，扩大自己的知识面。 有许多书籍还能培养我们的道德情操， 给我们巨大的精神力量， 鼓舞我们前进。 * 算法仿真示例 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 作用：使用积分滤波单元增加对高