三相短路故障计算
1.原始数据：如图表示有五个节点的系统结构接线图，网络参数在图中标出。

在简化的
短路电流计算中，假定发电机电势的标么值等于1，试求：（1）节点1、2、3三相短路时的各支路电流及各节点的电压；（2）在支路1——2的中点发生三相短路时的各支炉电流及各节点的电压；（3）断开支路1——2时，节点发生三相短路时的各支路电流及各节点电压。

其中： 120.25x j =，130.35x j =，230.3x j =，240.015x j =，350.03x j =，''40.02d x j =，''50.02d x j =
2.程序框图：
3.结果显示：
(1)在各节点处发生短路时的值1.2.3节点分别短路时的各值：阻抗矩阵：
0 + 0.1664i 0 + 0.0205i 0 + 0.0208i 0 + 0.0117i 0 + 0.0083i 0 + 0.0205i 0 + 0.0307i 0 + 0.0061i 0 + 0.0175i 0 + 0.0025i 0 + 0.0208i 0 + 0.0061i 0 + 0.0412i 0 + 0.0035i 0 + 0.0165i 0 + 0.0117i 0 + 0.0175i 0 + 0.0035i 0 + 0.0186i 0 + 0.0014i 0 + 0.0083i 0 + 0.0025i 0 + 0.0165i 0 + 0.0014i 0 + 0.0186i
短路点电流：
各节点电压：
各支路电流：
(2) 1-2支路中线短路
阻抗矩阵：
0.1664i 0.0205i 0.0208i 0.0117i 0.0083i 0.0934i 0.0205i 0.0307i 0.0061i 0.0175i 0.0025i 0.0256i 0.0208i 0.0061i 0.0412i 0.0035i 0.0165i 0.0135i 0.0117i 0.0175i 0.0035i 0.0186i 0.0014i 0.0146i 0.0083i 0.0025i 0.0165i 0.0014i 0.0186i 0.0054i 0.0934i 0.0256i 0.0135i 0.0146i 0.0054i 0.1220i
用上述方法可以得其他各值
短路电流I f=08.1957i
各节点电压V:
0.2341 0.7903 0.8898 0.8802 0.9559 0
各支路电流Izl:
0 - 2.5057i 0 +31.5804i 0 +16.1963i 0 + 5.1223i 0 + 4.0107i 0 -30.8900i
(3) 1-2断线
阻抗矩阵：
0 + 0.3935i 0 + 0.0045i 0 + 0.0435i 0 + 0.0026i 0 + 0.0174i 0 + 0.0045i 0 + 0.0318i 0 + 0.0045i 0 + 0.0182i 0 + 0.0018i 0 + 0.0435i 0 + 0.0045i 0 + 0.0435i 0 + 0.0026i 0 + 0.0174i 0 + 0.0026i 0 + 0.0182i 0 + 0.0026i 0 + 0.0190i 0 + 0.0010i 0 + 0.0174i 0 + 0.0018i 0 + 0.0174i 0 + 0.0010i 0 + 0.0190i 同理求取其余各值
短路节点电流If: 0 + 2.5413i
各节点电压V:
0 0.9884 0.8894 0.9934 0.9558
各支路电流Izl:
0 +20.4438i 0 -21.7822i 0 + 0.2723i 0 + 3.8119i 0
4.程序清单
%B=[1 2 0.25i 0 0;1 3 0.35i 0 0;2 3 0.3i 0 0;2 4 0.015i 0 1./0.02i;3 5 0.03i 0 1./0.02i]
%B=[1 6 1./2*0.25i 0 0;1 3 0.35i 0 0;2 3 0.3i 0 0;2 4 0.015i 0 1./0.02i;3 5 0.03i 0 1./0.02i;2 6 1./2*0.25i 0 0]
%B=[1 3 0.35i 0 0;2 3 0.3i 0 0;2 4 0.015i 0 1./0.02i;3 5 0.03i 0 1./0.02i] % ½Úµãµ¼ÄÉÕóµÄÇóÈ¡
clear
n=input('½Úµã¸öÊý:n=');
nl=input('֧·Êý:nl=');
B=input('֧·²ÎÊý:B=');
Yn=zeros(n);%YÕóÇåÁã
for i=1:nl
p=B(i,1);q=B(i,2);
y=B(i,3);s=B(i,4);
t=B(i,5);
if s==0
Yn(p,p)=Yn(p,p)+1./y;
Yn(q,q)=Yn(q,q)+1./y+t;
Yn(p,q)=Yn(p,q)-1./y;
Yn(q,p)=Yn(q,p)-1./y;
end
end%ÇóÈ¡µ¼ÄÉÕó
Zn=inv(Yn); %ÇóÈ¡×迹Õó
f=input('ÊäÈë¶Ì·½ÚµãºÅ;f=');%Çó¶Ì·µãµçÁ÷
disp('×迹¾ØÕóZn:');
disp(Zn);
Zn(f,f)=Zn(f,f);
If(f)=-1./Zn(f,f);
disp('¶Ì·½ÚµãµçÁ÷If:');
disp(If(f));
for i=1:n %Çó¸÷½Úµãµçѹ
V(i)=1-(-Zn(i,f))*If(f);
end
disp('¸÷½ÚµãµçѹV:');
disp(V);
for i=1:nl %Çó֧·µçÁ÷
p=B(i,1); q=B(i,2);
k=-(Zn(p,f)-Zn(q,f))./Zn(p,q);
Izl(i)=k*1./Zn(f,f);
end
disp('¸÷֧·µçÁ÷Izl:');
disp(Izl);。