TBL.zip_突变量_配电保护_配电网_配电网保护_配电网基于突变量保护的突变量提取

I1=[i1a11a,i1b11b,i1c11c]; I2=[i1a22a,i1b22b,i1c22c]; I3=[i1a33a,i1b33b,i1c33c]; I4=[i1a44a,i1b44b,i1c44c]; t=0:0.0001:0.02;%%%10000HZ x11=I1(1:100:125000,1); x1(1,:)=x11(401:600)-x11(201:400); x12=I1(1:100:125000,2); x1(2,:)=x12(401:600)-x12(201:400); x13=I1(1:100:125000,3); x1(3,:)=x13(401:600)-x13(201:400); x21=I2(1:100:125000,1); x2(1,:)=x21(401:600)-x21(201:400); x22=I2(1:100:125000,2); x2(2,:)=x22(401:600)-x22(201:400); x23=I2(1:100:125000,3); x2(3,:)=x23(401:600)-x23(201:400); x31=I3(1:100:125000,1); x3(1,:)=x31(401:600)-x31(201:400); x32=I3(1:100:125000,2); x3(2,:)=x32(401:600)-x32(201:400); x33=I3(1:100:125000,3); x3(3,:)=x33(401:600)-x33(201:400); x41=I4(1:100:125000,1); x4(1,:)=x41(401:600)-x41(201:400); x42=I4(1:100:125000,2); x4(2,:)=x42(401:600)-x42(201:400); x43=I4(1:100:125000,3); x4(3,:)=x43(401:600)-x43(201:400); xx1=x1(:,1:50) xx2=x2(:,1:50) xx3=x3(:,1:50) xx4=x4(:,1:50) for i=2:3 A=corrcoef(xx1(1,:),xx1(i,:)); c1(i-1)=A(1,2); end A=corrcoef(xx1(2,:),xx1(3,:)); c1(3)=A(1,2) for i=2:3 A=corrcoef(xx2(1,:),xx2(i,:)); c2(i-1)=A(1,2); end A=corrcoef(xx2(2,:),xx2(3,:)); c2(3)=A(1,2) for i=2:3 A=corrcoef(xx3(1,:),xx3(i,:)); c3(i-1)=A(1,2); end A=corrcoef(xx3(2,:),xx3(3,:)); c3(3)=A(1,2) for i=2:3 A=corrcoef(xx4(1,:),xx4(i,:)); c4(i-1)=A(1,2); end A=corrcoef(xx4(2,:),xx4(3,:)); c4(3)=A(1,2) %D=[mean(d1),mean(d2),mean(d3),mean(d3)]%零序电流相关系数 %D1=[c1(1)*c1(2)*c1(3),c2(1)*c2(2)*c2(3),c3(1)*c3(2)*c3(3),c4(1)*c4(2)*c4(3)] D=[min(abs(c1)),min(abs(c2)),min(abs(c3)),min(abs(c4))] t=t+0.04 subplot(411) plot(t(1:50),xx1(1,:),t(1:50),xx1(2,:),t(1:50),xx1(3,:)); subplot(412) plot(t(1:50),xx2(1,:),t(1:50),xx2(2,:),t(1:50),xx2(3,:)); subplot(413) plot(t(1:50),xx3(1,:),t(1:50),xx3(2,:),t(1:50),xx3(3,:)); subplot(414) plot(t(1:50),xx4(1,:),t(1:50),xx4(2,:),t(1:50),xx4(3,:));