微机继电保护实验报告
项目名称： 微机距离保护算法（2）
姓 名： 田务坤 学 号： K031241629 班 级： K0312416 实验时间： 2015年7月8日 实验地点： 实验楼五楼 实验成绩：
一、 实验目的
1．熟练掌握matlab/simulink的使用方法；
2．通过matlab/simulink输电网络模型与上次实验所编写的继电保护算法程序联合仿真，熟悉继电保护仿真研究的基本方法。
二、实验步骤
1.根据模型参数，计算两段距离保护整定值（距离I段和距离II段）。 2.利用powergui模块分析短路电流、电压中的谐波分量，并根据分析的结果修改实用一中编写的滤波程序。根据仿真数据的特征，对保护算法程序进行修改，要求程序中有距离I段和II段，距离继电器采用方向圆特性。
3.线路不同的位置（10%处，50%处，p处）作短路，观察保护的跳闸情况。
三、实验模型及程序
实验程序
T=VMabc.time;
a=VMabc.signals.values; b=IMabc.signals.values; t=T(:,1)'; UA=a(:,1)'; UB=a(:,2)'; UC=a(:,3)'; IA=b(:,1)'; IB=b(:,2)'; IC=b(:,3)'; Ts=t(1,2)-t(1,1); N=0.02/Ts; m=size(t);
IIA=zeros(1,m(2)); IIB=zeros(1,m(2)); IIC=zeros(1,m(2)); for jj=101:m(2);
IIA(jj)=(IA(jj)-IA(jj-100))/2; IIB(jj)=(IB(jj)-IB(jj-100))/2;
IIC(jj)=(IC(jj)-IC(jj-100))/2; end
plot(t,IIA,'r') figure
plot(t,IIB,'g') figure
plot(t,IIC,'b') figure plot(t,IA) UUA=zeros(1,m(2)); UUB=zeros(1,m(2)); UUC=zeros(1,m(2)); for jj=101:m(2);
UUA(jj)=(UA(jj)-UA(jj-100))/2; UUB(jj)=(UB(jj)-UB(jj-100))/2; UUC(jj)=(UC(jj)-UC(jj-100))/2; end figure
plot(t,UUA,'r') figure
plot(t,UUB,'g') figure
plot(t,UUC,'b') figure plot(t,UUA) USA=zeros(1,m(2)); USB=zeros(1,m(2)); USC=zeros(1,m(2)); for 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
%******?ˉ×÷?D?Y*******% Zset=1.273+i*29.31818; 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)-0.5*Zset)<=0.5*abs(Zset) flag(jj)=1; end