电力系统短路课程设计.docVIP

2 短路计算原理框图 利用MATLAB程序编程计算故障点的短路故障电流及故障点相电压。两相短路故障计算的原理框图如图2.1所示。 图2.1两相短路故障计算的原理框图 3 MATLAB程序设计 NF=input(请输入短路点数：NF=); for i1=1:NF clear; n1=input(请输入正序节点数:n1=); n2=input(请输入负序节点数:n2=); n0=input(请输入零序节点数:n0=); nl1=input(请输入正序支路数:nl1=); nl2=input(请输入负序支路数:nl2=); nl0=input(请输入零序支路数:nl0=); f=input(请输入短路号:f=); zf=input(请输入接地阻抗:zf=); zg=input(请输入接地阻抗:zg=); B1=input(请输入正序支路参数:B1=); B2=input(请输入负序支路参数:B2=); B0=input(请输入零序序支路参数:B0=); 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(nl0,1); I1=zeros(nl1,1);I2=zeros(nl2,1);I0=zeros(nl0,1); for m1=1:3 m=0; if m1==1 nl=nl1;B=B1;Z=Z1;n=n1; elseif m1==2 nl=nl2;B=B2;Z=Z2;n=n2; elseif m1==3 nl=nl0;B=B0;Z=Z0;n=n0; end m=0;Z=zeros(n);V=zeros(n);I=zeros(nl); for K1=1:nl p=B(K1,1);q=B(K1,2); if B(K1,6)==0 k=1./B(K1,5); else k=B(K1,5); end if p==0 if 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); end Z(m+1,m+1)=Z(q,q)+B(K1,3); for i=1:m for j=1:m Z(i,j)=Z(i,j)-Z(i,m+1)*Z(m+1,j)./Z(m+1,m+1); end Z(i,m+1)=0; end for i=1:m+1 Z(m+1,i)=0; end end else if qm %追加不接地树枝 for i=1:m Z(i,q)=Z(i,p)*k; Z(q,i)=Z(p,i)*k; end Z(q,q)=k^2*Z(p,p)+k^2*B(K1,3); m=m+1; else for 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); end Z(m+1,m+1)=k^2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k^2*B(K1,3); for i=1