电力系统

NF=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('请输入短路节点号: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.083i 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 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./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);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 q>m %追加接地树支 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 q>m %追加不接地树支 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: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;

elseif m1==3 Z0=Z end

end %以上程序求出正，负，零序对应的阻抗阵 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); end

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(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,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('故障点处的负序电流标么值为:'); disp(a2);

disp('故障点处的零序电流标么值为:'); disp(a0); for i=1:n1

V1(i)=Vcs(i)-I1(f)*Z1(i,f); end

for i=1:n2

V2(i)=-I2(f)*Z2(i,f); end

for i=1:n0

V0(i)=-I0(f)*Z0(i,f); end

disp('各节点的正序电压标么值为(节点号从小到大排):'); disp( V1);

disp('各节点的负序电压标么值为(节点号从小到大排):'); disp( V2);

disp('各节点的零序电压标么值为(节点号从小到大排):'); disp( V0);

for i=1:3 %求各支路的各序电流标么值 if i==1

B=B1;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('各支路的正序电流标幺值为(顺序同您输入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;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相电流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:nl

disp(Iazl(i)); end

disp('各支路B相电流Ibzl分别为(顺序同你输入B时一样):'); for i=1:nl

disp(Ibzl(i)); end

%本程序的功能是对简单不对称故障一般形式的计算方法 clc;clear;close NF=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+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.0098+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 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=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 i=1:n1

if X1(i,2)~=0; p=X1(i,1);

Y1(p,p)=1./X1(i,2); end end

for 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,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; end

disp('导纳矩阵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); end end

for i=1:nl2

if B2(i,6)==0

p=B2(i,1);q=B2(i,2); else

p=B2(i,2);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; end

disp('导纳矩阵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); end end

for i=1:nl0

if B0(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; end

disp('导纳矩阵Y0='); disp(Y0) Z0=inv(Y0)

disp(Z0)

end %以上程序求出零序对应的阻抗阵

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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); 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(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,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('故障点处的负序电流标么值为;'); disp(a2);

disp('故障点处的零序电流标么值为;'); disp(a0);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%求各节点各序的电压标么值

%注意，负、零序网中节点电压的求法 for i=1:n1

V1(i)=Vcs(i)-I1(f)*Z1(i,f); end

for i=1:n2

V2(i)=-I2(f)*Z2(i,f); end

for i=1:n0

V0(i)=-I0(f)*Z0(i,f); end

disp('各节点的正序电压标么值为(节点号从小到大排)：') disp(V1);

disp('各节点的负序电压标么值为(节点号从小到大排)：') disp(V2);

disp('各节点的零序电压标么值为(节点号从小到大排)：') disp(V0);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%求支路电流的程序与对称短路一样，只不过针对三个序网用不同的参数分别计算 for i=1:3 %求各支路序电流 if i==1

B=B1;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('各支路的正序电流标么值为(顺序同你输入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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%求故障点A,B,C三相各相的总电流，各个节点A,B,C三相各相的总电压，各支路A,B,C三相的总电流

%i=1时处理故障点电流，i=2时处理节点电压，i=3时处理支路电流 for i=1:3

if i==1

s0=a0;s1=a1;s2=a2; elseif i==2

s0=V0;s1=V1;s2=V2; else

s0=I0;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相电流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:nl

disp(Iazl(i)); end

disp('各支路B相电流Ibz1分别为(顺序同你输入B时一样):') for i=1:nl

disp(Ibzl(i)); end

disp('各支路C相电流Icz1分别为(顺序同你输入B时一样):') for i=1:nl

disp(Iczl(i)); end