电力系统.docx

1、电力系统%本程序的功能是对简单不对故障一般形式的计算方法clc;clear;closeNF=1;%input(请输入短路点的数目:NF=);for i1=1;NF clear; n1=9;%input(请输入正序节点数:n1=); n2=9;%input(请输入负序节点数:n2=); n0=9;%input(请输入零序节点数:n0=); nl1=9;%input(请输入正序支路数:nl1=); nl2=9;%input(请输入负序支路数:nl2=); nl0=9;%input(请输入零序支路数:nl0=); Lf=2;%input(请输入短路类型:Lf=); f=7;%input(请输入短路节

2、点号:f=); If=0;%input(请输入纵横故障标志;If=); zf=0;%input(请输入接地阻抗:zf=); zg=0;%input(请输入接地阻抗:zg=); B1=1 4 0.0476i 0 1 0;2 7 0.0474i 0 1 0;3 8 0.0484i 0 1 0;4 5 0.008+0.07i 0.1774i 1 0;4 6 0.014+0.075i 0.158i 1 0;5 7 0.026+0.133i 0.302i 1 0;6 8 0.032+0.14i 0.358i 1 0;7 9 0.007+0.0594i 0.15i 1 0;8 9 0.0098+0.083

3、i 0.21i 1 0;%input(请输入正序支路参数:B1=); B2=1 4 0.0476i 0 1 0;2 7 0.0474i 0 1 0;3 8 0.0484i 0 1 0;4 5 0.008+0.07i 0.1774i 1 0;4 6 0.014+0.075i 0.158i 1 0;5 7 0.026+0.133i 0.302i 1 0;6 8 0.032+0.14i 0.358i 1 0;7 9 0.007+0.0594i 0.15i 1 0;8 9 0.0098+0.083i 0.21i 1 0;%input(请输入负序支路参数:B2=); B0=1 4 0.1428i 0 1

4、 0;2 7 0.1422i 0 1 0;3 8 0.1452i 0 1 0;4 5 0.024+0.21i 0.5322i 1 0;4 6 0.042+0.225i 0.474i 1 0;5 7 0.078+0.399i 0.906i 1 0;6 8 0.096+0.42i 1.074i 1 0;7 9 0.021+0.1782i 0.45i 1 0;8 9 0.0294+0.249i 0.63i 1 0;%input(请输入零序支路参数:B0=); 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

5、./2-i*sqrt(3)/2;%input(请输入不对称系数:A1=); Vcs=1.0500;1.0430;1.0430;1.0376;1.0090;1.0268;1.0431; 1.0499;1.0340;%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)

6、;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;

7、 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)=k2*Z(p,p)+k2*B(k1,

8、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)=k2*Z(p,p)+Z(q,q)-2*k*Z(p,q)+k2*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 end end if m1=1 Z1=Z; elseif m1=2 Z2=Z; e

9、lseif m1=3 Z0=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)=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

10、)=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=2 Z4=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);elseif Lf=3 Z4=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)+

11、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);elseif Lf=4 Z4=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)./

12、(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=5 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);enddisp(故障点处的正序电流标么值为:);disp(a1);disp(故障点处的负序电流标么值为:);disp(

13、a2);disp(故障点处的零序电流标么值为:);disp(a0);for i=1:n1 V1(i)=Vcs(i)-I1(f)*Z1(i,f);endfor i=1:n2 V2(i)=-I2(f)*Z2(i,f);endfor i=1:n0 V0(i)=-I0(f)*Z0(i,f);enddisp(各节点的正序电压标么值为(节点号从小到大排):);disp( V1);disp(各节点的负序电压标么值为(节点号从小到大排):);disp( V2);disp(各节点的零序电压标么值为(节点号从小到大排):);disp( V0);for i=1:3 %求各支路的各序电流标么值 if i=1 B=B1

14、;V=V1;nl=nl1; elseif i=2 B=B2;V=V2;nl=nl2; elseif i=3 B=B0;V=V0;nl=nl0; end for j=1:nl if B(j,6)=0 k=B(j,5); else k=1./B(j,5); end p=B(j,1);q=B(j,2); if p=0 e=0;b=B(j,3); I(j)=(e-V(q)./k)./b; else I(j)=(V(p)-V(q)./k)./B(j,3); end end if i=1 I1=I; elseif i=2 I2=I; elseif i=3 I0=I; end end disp(各支路的正序

15、电流标幺值为(顺序同您输入B时一样):); for i=1:nl disp(I1(i); end disp(各支路负序电流的标幺值为(顺序同您输入B时一样):); for i=1:nl disp(I2(i); end disp(各支路零序电流的标幺值为(顺序同您输入B时一样):); for i=1:nl disp(I0(i); end for i=1:3 %求故障点A,B,C三相各相的总电流，各个节点A,B,C三相各相的总电压，各支路A,B,C三相各相的总电流 if i=1 s0=a0;s1=a1;s2=a2; elseif i=2 s0=V0;s1=V1;s2=V2; else s0=I0;

16、s1=I1;s2=I2; end for j=1:3 CA=s0*A1(j,1)+s1*A1(j,2)+s2*A1(j,3); if j=1 A=CA; elseif j=2 B=CA; elseif j=3 C=CA; end end if i=1 Iadl=A;Ibdl=B;Icdl=C; elseif i=2 Vajy=A;Vbjy=B;Vcjy=C; else Iazl=A;Ibzl=B;Iczl=C; end end disp(故障点处A相电流Iadl为:); disp(Iadl); disp(故障点处B相电流Ibdl为:); disp(Ibdl); disp(故障点处C相电流Icd

17、l为:); disp(Icdl); disp(各节点处A相电压Vajy分别为(节点号从小到大排):); disp(Vajy); disp(各节点处B相电压Vbjy分别为(节点号从小到大排):); disp(Vbjy); disp(各节点处C相电压Vcjy分别为(节点号从小到大排):); disp(Vcjy); disp(各支路A相电流Iazl分别为(顺序同你输入B时一样):); for i=1:nl disp(Iazl(i); end disp(各支路B相电流Ibzl分别为(顺序同你输入B时一样):); for i=1:nl disp(Ibzl(i); end%本程序的功能是对简单不对称故障一

18、般形式的计算方法clc;clear;closeNF=1;n1=9;n2=9;n0=9;nl1=9;nl2=9;nl0=9;Lf=2;f=7;If=0;zf=0;zg=0;%各序网的支路参数B1=1 4 0.0476i 0 1 0;2 7 0.0474i 0 1 0;3 8 0.0484i 0 1 0;4 5 0.008+0.07i 0.1774i 1 0;4 6 0.014+0.075i 0.158i 1 0;5 7 0.026+0.133i 0.302i 1 0;6 8 0.032+0.14i 0.358i 1 0;7 9 0.007+0.0594i 0.15i 1 0;8 9 0.0098

19、+0.083i 0.21i 1 0;X1=1 0.0424i; 2 0.065i; 3 0.0933i;4 0; 5 0.7428i; 6 1.0541i;7 0;8 0; 9 0.9438i;B2=1 4 0.0476i 0 1 0;2 7 0.0474i 0 1 0;3 8 0.0484i 0 1 0;4 5 0.008+0.07i 0.1774i 1 0;4 6 0.014+0.075i 0.158i 1 0;5 7 0.026+0.133i 0.302i 1 0;6 8 0.032+0.14i 0.358i 1 0;7 9 0.007+0.0594i 0.15i 1 0;8 9 0.0

20、098+0.083i 0.21i 1 0;X2=1 0.0424i; 2 0.065i; 3 0.0933i;4 0; 5 0.7428i; 6 1.0541i;7 0;8 0; 9 0.9438i;B0=1 4 0.1428i 0 1 0;2 7 0.1422i 0 1 0;3 8 0.1452i 0 1 0;4 5 0.024+0.21i 0.5322i 1 0;4 6 0.042+0.225i 0.474i 1 0;5 7 0.078+0.399i 0.906i 1 0;6 8 0.096+0.42i 1.074i 1 0;7 9 0.021+0.1782i 0.45i 1 0;8 9

21、0.0294+0.249i 0.63i 1 0;X0=1 0; 2 0; 3 0;4 0; 5 0.7428i; 6 1.0541i;7 0;8 0; 9 0.9438i;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.0500;1.0430;1.0430;1.0376;1.0090;1.0268;1.0431; 1.0499;1.0340; Z1=zeros(n1);Z2=zeros(n2);Z0=zeros(n0); Y1=zeros(n1);Y2=ze

22、ros(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 i=1:n1 if X1(i,2)=0; p=X1(i,1); Y1(p,p)=1./X1(i,2); endendfor i=1:nl1 if B1(i,6)=0 p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y1(p,q)=Y1(p,q)-1./(B1(i,

23、3)*B1(i,5); Y1(q,p)=Y1(p,q); Y1(q,q)=Y1(q,q)+1./(B1(i,3)*B1(i,5)2)+B1(i,4)./2; Y1(p,p)=Y1(p,p)+1./B1(i,3)+B1(i,4)./2;enddisp(导纳矩阵Y1=);disp(Y1)Z1=inv(Y1)disp(Z1)end %以上程序求出正序对应的阻抗阵for i=1:n2 if X2(i,2)=0; p=X2(i,1); Y2(p,p)=1./X2(i,2); endendfor i=1:nl2 if B2(i,6)=0 p=B2(i,1);q=B2(i,2); else p=B2(i,2

24、);q=B2(i,1); end Y2(p,q)=Y2(p,q)-1./(B2(i,3)*B2(i,5); Y2(q,p)=Y2(p,q); Y2(q,q)=Y2(q,q)+1./(B2(i,3)*B2(i,5)2)+B2(i,4)./2; Y2(p,p)=Y2(p,p)+1./B2(i,3)+B2(i,4)./2;enddisp(导纳矩阵Y2=);disp(Y2)Z2=inv(Y2)disp(Z2)end %以上程序求出负序对应的阻抗阵for i=1:n0 if X0(i,2)=0; p=X0(i,1); Y0(p,p)=1./X0(i,2); endendfor i=1:nl0 if B0

25、(i,6)=0 p=B0(i,1);q=B0(i,2); else p=B0(i,2);q=B0(i,1); end Y0(p,q)=Y0(p,q)-1./(B0(i,3)*B0(i,5); Y0(q,p)=Y0(p,q); Y0(q,q)=Y0(q,q)+1./(B0(i,3)*B0(i,5)2)+B0(i,4)./2; Y0(p,p)=Y0(p,p)+1./B0(i,3)+B0(i,4)./2;enddisp(导纳矩阵Y0=);disp(Y0)Z0=inv(Y0)disp(Z0)end %以上程序求出零序对应的阻抗阵%if If=0 %求出横向故障的故障端口输入阻抗 Z1(f,f)=Z1(

26、f,f) Z2(f,f)=Z2(f,f) Z0(f,f)=Z0(f,f)else %求出纵向故障的故障端口输入阻抗 Z1(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);end% %求出正负零序短路节点的电流标么值 %Z4即正序等效定则中的附加阻抗 if 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(

27、f);a1=I1(f); a2=I2(f);a0=I0(f); elseif Lf=2 Z4=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); elseif Lf=3 Z4=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,

28、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); elseif Lf=4 Z4=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=5 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); end disp(故障点处的正序电流标么值为;); disp(a1); disp(故障点处的负序电