電力系統不對稱故障計算的Matlab算法程序電力系統不對稱故障計算的Matlab算法程序電力系統不對稱故障計算的Matlab算法程序xxx公司電力系統不對稱故障計算的Matlab算法程序文件編號：文件日期：修訂次數：第1.0次更改批準審核制定方案設計，管理制度NF=1;fori1=1:NFclear;n1=9;n2=9;n0=9;nl1=15;nl2=15;nl0=15;Lf=input('請輸入短路類型（單相接地短路輸入‘1’，兩相短路輸入‘2’，兩相短路接地輸入‘f=7;If=0;zf=0;zg=0;B1=[01010;02010;03010;1400;05010;06010;2700;08010;3900;45011;46011;57011;69011;78011;89011];B2=[01010;02010;03010;1400;05010;06010;2700;08010;3900;45011;46011;57011;69011;78011;89011];B0=[01010;02010;03010;1400;05010;06010;2700;08010;3900;45*3011;46*3011;57*3011;69*3011;78*3011;89*3011];A1=[111;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=[;;;;;;;;];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);form1=1:3m=0;ifm1==1nl=nl1;B=B1;Z=Z1;n=n1;elseifm1==2nl=nl2;B=B2;Z=Z2;n=n2;elseifm1==3nl=nl0;B=B0;Z=Z0;n=n0;endm=0;Z=zeros(n);V=zeros(n);I=zeros(nl);fork1=1:nlp=B(k1,1);q=B(k1,2);ifB(k1,6)==0k=1./B(k1,5);elsek=B(k1,5);endifp==0ifq>mZ(q,q)=B(k1,3);m=m+1;elsefori=1:mZ(i,m+1)=-Z(i,q);Z(m+1,i)=-Z(q,i);endZ(m+1,m+1)=Z(q,q)+B(k1,3);fori=1:mforj=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;endfori=1:m+1Z(m+1,i)=0;endendelseifq>mfori=1:mZ(i,q)=Z(i,p)*k;Z(q,i)=Z(p,i)*k;endZ(q,q)=k^2*Z(p,p)+k^2*B(k1,3);m=m+1;elsefori=1: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^2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k^2*B(k1,3);fori=1:mforj=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;endfori=1:m+1Z(m+1,i)=0;endendendendifm1==1Z1=Z;elseifm1==2Z2=Z;elseifm1==3Z0=Z;endendifIf==0Z1(f,f)=Z1(f,f);Z2(f,f)=Z2(f,f);Z0(f,f)=Z0(f,f);elseZ1(f,f)=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);endifLf==1Z4=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);elseifLf==2Z4=Z2(f,f)+2*zf;k2=-1;k0=0;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);elseifLf==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*I1(f);a1=I1(f);a2=I2(f);a0=I0(f);elseifLf==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);elseifLf==5Z4=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);enddisp('故障點處的正序電流標幺值為：');disp(a1);disp('故障點處的負序電流標幺值為：');disp(a2);disp('故障點處的零序電流標幺值為：');disp(a0);fori=1:n1V1(i)=Vcs(i)-I1(f)*Z1(i,f);endfori=1:n2V2(i)=-I2(f)*Z2(i,f);endfori=1:n0V0(i)=-I0(f)*Z0(i,f);enddisp('各節點的正序電壓標幺值為（節點號從小到大排）：');disp(V1);disp('各節點的負序電壓標幺值為（節點號從小到大排）：');disp(V2);disp('各節點的零序電壓標幺值為（節點號從小到大排）：');disp(V0);fori=1:3ifi==1B=B1;V=V1;nl=nl1;elseifi==2B=B2;V=V2;nl=nl2;elseifi==3B=B0;V=V0;nl=nl0;endforj=1:nlifB(j,6)==0k=B(j,5);elsek=1./B(j,5);endp=B(j,1);q=B(j,2);ifp==0e=0;b=B(j,3);I(j)=(e-V(q)./k)./b;elseI(j)=(V(p)-V(q)./k)./B(j,3);endendifi==1I1=I;elseifi==2I2=I;elseifi==3I0=I;endenddisp('各支路的正序電流標幺值為（順序同您輸入B時一樣）;');fori=1:nldisp(I1(i));enddisp('各支路的負序電流標幺值為（順序同您輸入B時一樣）;');fori=1:nldisp(I2(i));enddisp('各支路的零序電流標幺值為（順序同您輸入B時一樣）;');fori=1:nldisp(I0(i));endfori=1:3ifi==1s0=a0;s1=a1;s2=a2;elseifi==2s0=V0;s1=V1;s2=V2;elses0=I0;s1=I1;s2=I2;endforj=1:3CA=s0*A1(j,1)+s1*A1(j,2)+s2*A1(j,3);ifj==1A=CA;elseifj==2B=CA;elseifj==3C=CA;endendifi==1Iadl=A;Ibdl=B;Icdl=C;elseifi==2Vajy=A;Vbjy=B;Vcjy=C;elseIazl=A;Ibzl=B;Iczl=C;endenddisp('故障點處A相電流Iadl為：');disp(Iadl);disp('故障點處B相電流Ibdl為：');disp(Ibdl);disp('故障點處C相電流Icdl為：');disp(Icdl);disp('各節點處A相電壓Vajy分別為(節點號從小到大排）：');disp(Vajy);disp('各節