微机保护实验报告.docx

The Standardization Office was revised on the afternoon of December 13, 2020
微机保护实验报告
微机继电保护实验报告
or jj=N/4+1:m(2)
USA(jj)=sqrt((UUA(jj)*UUA(jj)+UUA(jj-N/4)*UUA(jj-N/4))/2);
USB(jj)=sqrt((UUB(jj)*UUB(jj)+UUB(jj-N/4)*UUB(jj-N/4))/2);
USC(jj)=sqrt((UUC(jj)*UUC(jj)+UUC(jj-N/4)*UUC(jj-N/4))/2);
end
% %电流
for jj=N/4+1:m(2)
ISA(jj)=sqrt((IIA(jj)*IIA(jj)+IIA(jj-N/4)*IIA(jj-N/4))/2);
ISB(jj)=sqrt((IIB(jj)*IIB(jj)+IIB(jj-N/4)*IIB(jj-N/4))/2);
ISC(jj)=sqrt((IIC(jj)*IIC(jj)+IIC(jj-N/4)*IIC(jj-N/4))/2);
end
%定义测量电压和测量电流
Um=UUA-UUB;
Im=IIA-IIB;
Um1=UUB-UUC;
Im1=IIB-IIC;
Um2=UUC-UUA;
Im2=IIC-IIA;
% %电阻、电抗、相角差
for jj=N/4+1:m(2)
R(jj)=(Um(jj)*Im(jj)+Um(jj-N/4)*Im(jj-N/4))/(Im(jj)*Im(jj)+Im(jj-N/4)*Im(jj-N/4));
X(jj)=(Um(jj-N/4)*Im(jj)-Um(jj)*Im(jj-N/4))/(Im(jj)*Im(jj)+Im(jj-N/4)*Im(jj-N/4));
O(jj)=180/pi*atan((Um(jj-N/4)*Im(jj)-Um(jj)*Im(jj-N/4))/(Um(jj)*Im(jj)+Um(jj-N/4)*Im(jj-N/4)));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
R1(jj)=(Um1(jj)*Im1(jj)+Um1(jj-N/4)*Im1(jj-N/4))/(Im1(jj)*Im1(jj)+Im1(jj-N/4)*Im1(jj-N/4));
X1(jj)=(Um1(jj-N/4)*Im1(jj)-Um1(jj)*Im1(jj-N/4))/(Im1(jj)*Im1(jj)+Im1(jj-N/4)*Im1(jj-N/4));
O1(jj)=180/pi*atan((Um1(jj-N/4)*Im1(jj)-Um1(jj)*Im1(jj-N/4))/(Um1(jj)*Im1(jj)+Um1(jj-N/4)*Im1(jj-N/4)));
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
R2(jj)=(Um2(jj)*Im2(jj)+Um2(jj-N/4)*Im2(jj-N/4))/(Im2(jj)*Im2(jj)+Im2(jj-N/4)*Im2(jj-N/4));
X2(jj)=(Um2(jj-N/4)*Im2(jj)-Um2(jj)*Im2(jj-N/4))/(Im2(jj)*Im2(jj)+Im2(jj-N/4)*Im2(jj-N/4));
O2(jj)=180/pi*atan((Um2(jj-N/4)*Im2(jj)-Um2(jj)*Im2(jj-N/4))/(Um2(jj)*Im2(jj)+Um2(jj-N/4)*Im2(jj-N/4)));
end
%******动作判据*******%
Zset=+i*;
Zm=R+i.*X;
Zm1=R1+i.*X1;
Zm2=R2+i.*X2;
flag=zeros(1,m(2));
flag1=zeros(1,m(2));
flag2=zeros(1,m(2));
for jj=1:m(2)
if abs(Zm(jj)*Zset)<=*abs(Zset)
flag(jj)=1;
end
if abs(Zm1(jj)*Zset)<=*abs(Zset)
flag1(jj)=1;
end
if abs(Zm2(