电力系统短路计算程序

1、对称故障的计算：NF=input(请输入短路点的数目:NF=);n=input(请输入独立节点数:n=);nl=input(请输入支路数:nl=);B=input(请输入由线路参数形成的矩阵:B=);V0=input(请输入由各节点的初电压标幺值形成的列矩阵:V0=);D=input(请输入由短路号，短路点阻抗形成的矩阵:D=);m=0;Z=zeros(n);V=zeros(n);I=zeros(nl);for k1=1:nlp=B(k1,1);q=B(k1,2);if B(k1,6)=0k=1./B(k1,5);else k=B(k1,5);endif p=0if qmZ(q,q)=B(k1

2、,3);m=m+1;elsefor i=1:m,Z(i,m+1)=-Z(i,q);Z(m+1,i)=-Z(q,i);endZ(m+1,m+1)=Z(q,q)+B(k1,3);for i=1:mfor j=1:mZ(i,j)=Z(i,j)-Z(i,m+1)* Z(m+1,j)./ Z(m+1,m+1);endZ(i,m+1)=0;endfor i=1:m+1Z(m+1,i)=0;endendelse if qmfor i=1:mZ(i,q)= Z(i,p)*k; Z(q,i)= Z(p,i)*k;endZ(q,q)=kA2*Z(p,p)+kA2*B(k1,3);m=m+1;elsefor i=1

3、:mZ(i,m+1)=k*Z(i,p)- Z(i,q);Z(m+1,i)=k*Z(p,i)- Z(q,i);endZ(m+1,m+1)= k*Z(p,p)+Z(q,q)-2*k*Z(p,q)+ kA2*B(k1,3);for i=1:mfor j=1:mZ(i,j)= Z(i,j)- Z(i,m+1)* Z(m+1,j)./ Z(m+1,m+1);endZ(i,m+1)=0;endfor i=1:m+1Z(m+1,i)=0;endendendendfor k=1:NFI(D(k,1),D(k,1)=V0(D(k,1),1)./(Z(D(k,1),D(k,1)+D(k,2);ft=num2str

4、(D(k,1);ts1=(点短路时);ts2=(电流的标幺值);dn=strcat(ft,ts1,ts2);disp(dn);disp(I(D(k,1),D(k,1);endfor i=1:nV(i,i)=V0(i,1)-I(D(k,1),D(k,1)*Z(i,D(k,1);ts3=(点短路时);ts4=（电压的标幺值）;dn1=strcat（i,ts3,ts4）;disp（dnl）;disp（V（i,i）;end简单不对称故障的计算：NF=input（请输入短路点的数目:NF=）;for i1=1:NFclear;n1= input（请输入正序节点数:n1=）;n2= input（请输入负序

5、节点数:n2=）;n0= input（请输入零序节点数:n0=）;nl1= input（请输入正序支路数:nl1=）;nl2= input（请输入负序支路数:nl2=）;nl0= input（请输入零序支路数:nl0=）;Lf= input（请输入短路类型:Lf=）;f= input（请输入短路节点号:f=）;If= input（请输入纵横故障标志:If =）;zf= input（请输入接地阻抗:Zf=）;zg= input（请输入接地阻抗:zg=）;B1= input（请输入正序支路参数:B1=）;B2= input（请输入负序支路参数:B2=）;B0= input（请输入零序支路参数:B0

6、=）;A1= input（请输入不对称系数:A1=）;Vcs= input（请输入由各节点的初始电压标幺值形成的列矩阵:Vcs=）;Z1=zeros（n1）;Z2=zeros（n2）;Z0=zeros（n0）;Y1=zeros（n1）;Y2=zeros（n2）;Y0=zeros（n0）;V1=zeros（n1,1）;V2=zeros（n2,1）;V0=zeros（n0,1）;I1=zeros（nl1,1）; I2=zeros（nl2,1）; I0=zeros（nl0,1）;for m1=1:3m=0;if m1=1nl=nl1;B=B1;Z=Z1;n=n1;else if m1=2nl=nl2

7、;B=B2;Z=Z2;n=n2;else if m1=3nl=nl0;B=B0;Z=Z0;n=n0;endm=0;Z=zeros(n);V=zeros(n);I=zeros(nl);for k1=1:nlp=B(k1,1);q=B(k1,2);if B(k1,6)=0k=1./B(k1,5);else k=B(k1,5);endif p=0if qm%追加接地树支Z(q,q)=B(k1,3);m=m+1;else%追加接地连支for i=1:m,Z(i,m+1)=-Z(i,q);Z(m+1,i)=-Z(q,i);endZ(m+1,m+1)=Z(q,q)+B(k1,3);for i=1:mfor

8、 j=1:mZ(i,j)= Z(i,j)- Z(i,m+1)* Z(m+1,j)./ Z(m+1,m+1);endZ(i,m+1)=0;endfor i=1:m+1Z(m+1,i)=0;endendelse if qm%追加不接地树支for i=1:mZ(i,q)=Z(i,p)*k; Z(q,i)=Z(p,i)*k;endZ(q,q)=kA2*Z(p,p)+ k人2*B(k1,3);m=m+1;elsefor i=1:m%追加不接地连支Z(i,m+1)=k*Z(i,p)-Z(i,q);Z(m+1,i)=k*Z(p,i)-Z(q,i);endZ(m+1,m+1)=kA2*Z(p,p)+Z(q,q

9、)-2*k*Z(p,q)+kA2*B(k1,3);for i=1:mfor j=1:mZ(i,j)= Z(i,j)- Z(i,m+1)* Z(m+1,j)./ Z(m+1,m+1);endZ(i,m+1)=0;endfor i=1:m+1Z(m+1,i)=0;endendendendif m1=1Z1=Z;elseif m1=2Z2=Z;else if m1=3Z0=Z;endend%以上程序求出正，负，零序对应的阻抗阵if If=0%求出横向故障的阻抗阵Z1(f,f)=Z1(f,f);Z2(f,f)=Z2(f,f);Z0(f,f)=Z0(f,f);else%求出纵向故障的阻抗阵Z1(f,f)

10、= Z1(f,f)+ Z1(1,1)-2* Z1(f,1);Z2(f,f)= Z2(f,f)+ Z2(1,1)-2* Z2(f,1);Z0(f,f)= Z0(f,f)+ Z0(1,1)-2* Z0(f,1);endif Lf=1%求出正负零序短路节点的电流标幺值Z4=Z2(f,f)+Z0(f,f)+3*zf;k2=1;k0=1;I1(f)=Vcs (f) ./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=2Z4=Z2(f,f)+2*zf;k2=-1;k0=0;I1(f)=Vcs(f

11、)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);else if Lf=3Z4=zf+(Z2(f,f)+zf)*(Z0(f,f)+zf+3*zg)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);k2=-(Z0(f,f)+zf+3*zg)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);k0=-(Z2(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf+3*zg);I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I

12、1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=4Z4=zf+(Z2(f,f)+zf)*(Z0(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf)k2=-(Z0(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf);k0=-(Z2(f,f)+zf)./(Z2(f,f)+Z0(f,f)+2*zf);I1(f)=Vcs(f)./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseif Lf=5Z4=Z2(f,f)+Z0(f,f)+3*zf;k

13、2=1;k0=1;I1(f)=Vcs (f) ./(Z1(f,f)+Z4);I2(f)=k2*I1(f);I0(f)=k0*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);enddisp(故障点处的正序电流标幺值为：);disp(al);disp(故障点处的负序电流标幺值为：);disp(a2);disp(故障点处的零序电流标幺值为：);disp(a0);for i=1:nl%求各节点各序的电压标幺值V1(i)=Vcs(i)-I1(f)*Z1(i,f);endfor i=1:n2V2(i)=-I2(f)*Z2(i,f);endfor i=1:n0V0(i)=-I0(f)*Z0

14、(i,f);enddisp(各节点的正序电压标幺值为(节点号从小到大排)：);disp(V1);disp(各节点的负序电压标幺值为(节点号从小到大排)：);disp(V2);disp(各节点的零序电压标幺值为(节点号从小到大排)：);disp(V0);for i=1:3%求各节点的各序电流标幺值if i=1B=B1;V=V1;nl=nl1;elseif i=2B=B2;V=V2;nl=nl2;elseif i=3B=B0;V=V0;nl=nl0;endfor j=1:nlif B(j,6)=0k=B(j,5);else k=1./B(j,5);endp=B(j,1);q=B(j,2);if p

15、=0e=0;b=B(j,3);I=(e-V(q)./k)./b;elseI(j)=(V(p)-V(q)./k)./B(j,3);endendif i=1I1=I;elseif i=2I2=I;elseif i=3I0=I;endenddisp(各支路的正序电流标幺值为(顺序同您输入B时一样)：);for i=1:nldisp(I1(i);enddisp(各支路的负序电流标幺值为(顺序同您输入B时一样)：);for i=1:nldisp(I2(i);enddisp(各支路的零序电流标幺值为(顺序同您输入B时一样)：);for i=1:nldisp(I0(i);endfor i=1:3%求故障点A

16、，B，C三相各相的总电流，各个节点A，B，C三相各相 的总电压，各支路A，B，C三相各相的总电流if i=1s0=a0;s1=a1;s2=a2;elseif i=2s0=V0;s1=V1;s2=V2;elses0=I0；s1=I1;s2=I2;endfor j=1:3CA=s0*A1(j,1)+s1*A1(j,2)+s2*A1(j,3);if j=1A=CA;else if j=2B=CA;else if j=3C=CA;endendif i=1Iadl=A;Ibdl=B;Icdl=C;else if i=2Vajy=A;Vbjy=B;Vcjy=C;elseIazl=A;Ibzl=B;Iczl

17、=C;endenddisp(故障点处A相电流Iadl为：);disp(Iadl);disp(故障点处B相电流Ibdl为：);disp(Ibdl);disp(故障点处C相电流Icdl为：);disp(Icdl);disp(各节点处A相电压Vajy分别为(节点号从小到大排)：);disp(Vajy);disp(各节点处B相电压Vbjy分别为(节点号从小到大排)：);disp(Vbjy);disp(各节点处C相电压Vcjy分别为(节点号从小到大排)：);disp(Vcjy);disp(各支路A相电流Iazl分别为(顺序同您输入B时一样)：);for i=1:nldisp(Iazl(i);enddis

18、p(各支路B相电流Ibzl分别为(顺序同您输入B时一样)：);for i=1:nldisp(Ibzl(i);enddisp(各支路C相电流Iczl分别为(顺序同您输入B时一样)：);for i=1:nldisp(Iczl(i);endend请输入短路点的数目:NF=1请输入正序节点数:n1=1请输入负序节点数:n2=1请输入零序节点数:n0=1请输入正序支路数:nl1=2请输入负序支路数:nl2=2请输入零序支路数:nl0=2请输入短路类型:Lf=1请输入短路节点号:f=1请输入纵横故障标志:If =0请输入接地阻抗:zf=0请输入接地阻抗:zg=0 TOC o 1-5 h z 请输入正序支路

19、参数:B1=0 1 1.055i 0 1 0;011.055i010请输入负序支路参数:B2=0 1 0.555i 0 1 0;010.555i010请输入零序支路参数:B0=0 1 0.555i 0 1 0;010.555i010请输 入 不 对 称系数:A1=1 1 1;1 -1./2-i*sqrt(3)/2-1./2+i*sqrt(3)/2;1-1./2+i*sqrt(3)/2 -1./2-i*sqrt(3)/2请输入由各节点的初始电压标幺值形成的列矩阵:Vcs=1.08;1.08;1故障点处的正序电流标幺值为：0 - 0.9977i故障点处的负序电流标幺值为：0 - 0.9977i故障点处的零序电