电力系统分析实验指导书(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)=Z(i,p)*k-Z(i,q); Z(m+1,i)=Z(p,i)*k-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

34

end end end

for k=1:NF

I(D(k,1),D(k,1))=V0(D(k,1),1)./(Z(D(k,1), D(k,1))+D(k,2)); ft=num2str(D(k,1)); ts1=('点短路时');

ts2=('电流的标么值If='); dn=strcat(ft,ts1,ts2); disp(dn);

disp(I(D(k,1),D(k,1))); for i=1:n

V(i,i)=V0(i,1)-I(D(k,1),D(k,1))*Z(i,D(k,1));%求各支路的短路电流标么值 end for i=1:nl

if B(i,6)==0 k=B(i,5); else k=1./B(i,5); end

p=B(i,1);q=B(i,2); if p==0

e=0;b=B(i,3);

I(i,i)=(e-V(q,q)./k)./b; else

I(i,i)=(V(p,p)-V(q,q)./k)./B(i,3); end end

disp('各节点的电压标么值U为（节点号从小到大排）：'); for i=1:n disp(V(i,i)); end

disp('各之路短路电流的标么值I为（顺序同您输入B时一样）：'); for i=1:nl disp(I(i,i)); end end 四 实验数据

所示电力系统，负荷全部略去，简化后的各电抗标么值注于等值网络中。试计算f点三相短路时的短路电流及网络中的电流分布。

35

B=[0 1 0.2i 0 1 0;0 2 4i 0 1 0;1 3 0.51i 0 1 0;2 3 0.59i 0 1 0;3 4 1.43i 0 1 0]; V0=[1;1;1;1]; D=[4,0];

五、实验过程

1.手算图1所示网络的等值电抗。

2.认真阅读源程序，建立matlab的M程序，输入短路电流计算程序。 3.运行程序，输入实验数据：

矩阵B的每行是由下列参数构成的：

? 某支路的首端号P；

? 某支路末端号Q，且P

? 变比折算到哪一侧的标志（如果支路的首端P处于高压侧则请输入“1”， 否则请输入“0”）。 矩阵D：由短路点号，短路点阻抗

列矩阵V0：由各节点的初电压标么值形成 4.分析计算结果，检验程序的正确性。 六、要求

认真书写实验报告并上交实验数据

请输入短路点的数目：NF=1 请输入节点数:n=4 请输入支路数:nl=5

请输入由各支路参数形成的矩阵:B=[0 1 0.2i 0 1 0;0 2 4i 0 1 0;1 3 0.51i 0 1 0;2 3 0.59i 0 1 0;3 4 1.43i 0 1 0];

请输入由各节点的初电压标么值形成的列矩阵: V0=[1;1;1;1];

36

请输入由短路号，短路点阻抗组成的矩阵D=[4,0];

37