微机保护算法综合仿真

波器后,电抗误差明显降低,电阻误差上升。附:AB3.mclearcloseallloadxiti1.mat;k=m';t0=k(:,1);t=t0;tmax=max(t);tmin=min(t);[xy]=size(k);T=0.02;N=round((x*T)/(tmax-tmin));m=k(:,2:7);va=m(:,1);vb=m(:,2);vc=m(:,3);ia=m(:,4);ib=m(:,5);ic=m(:,6);i=(ia-ib);v=(va-vb);%Tukey低通滤波a=1;b=[0134310];i1=filter(b,a,i)/4;v1=filter(b,a,v)/4;i=i1;v=v1;%R-L模型算法[a,b]=size(i);i1=[0i(1:a-1)];i2=[00i(1:a-2)];v1=[0v(1:a-1)];v2=[00v(1:a-2)];r=((i-i1).*(v1+v2)-(v+v1).*(i1-i2))./((i-i1).*(i1+i2)-(i+i1).*(i1-i2));x=pi*((v+v1).*(i1+i2)-(i+i1).*(v1+v2))./((i-i1).*(i1+i2)-(i+i1).*(i1-i2))/N%控制阻抗轨迹显示范围m=[t,r,x];t1=0.08;t2=0.14;a1=t1*N/0.02;a2=t2*N/0.02;m1=m(a1:a2,:);t=m1(:,1);r=m1(:,2);x=m1(:,3);figuresubplot(221);plot(t,r,'k-o');ylabel('r(t)');subplot(223);plot(t,x,'k-o');ylabel('x(t)');xlabel('t/ms');subplot(222);plot(t,r,'k-o');ylabel('r(t)');