电力系统导纳矩阵与短路计算matlab程序

系统接线示于题图1，已知各元件参数见图上标注：
发电机G-1：＝120MV·A，＝0.23；G-2：＝60MV·A，＝0.14。
变压器T-1：＝120MV·A，＝10.5％；T-2：＝60MV·A，＝10.5％。
线路参数： , 。线路长度L-1：120km，L-2：80km，L-3：70km，L-4: 90km。取＝120MV·A，，试求标幺制下的节点导纳矩阵。另试就下列两种情况分别修改节点导纳矩阵：（1）节点3发生三相短路；（2）线路L-3中间发生三相短路。

matlab程序如下：

// 稳态时导纳矩阵计算
clc;
z=[0.23i        0.105i         inf          inf               inf             inf;
0.105i       -32.3939i   0.1452+0.2903i     inf          0.2178+0.4355i        inf;
inf       0.1452+0.2903i     -43.1965i    0.1271+0.2540i      inf             inf;
inf           inf        0.1271+0.2540i   -40.5022i     0.16335+0.3266i       inf;
inf       0.2178+0.4355i       inf        0.16335+0.3266i   -30.8547i          0.21i;  
inf           inf             inf         inf               0.21i            0.28i];
Y=zeros(6,6);
for(i=1:6)
for(j=1:6)
Y(i,i)=sum(1./z(i,:),2);
if j==i
Y(i,j)=Y(i,i);
else 
Y(i,j)=-1/z(i,j);
end
end
end
Y

结果如下：

第三个节点短路时matlab程序如下：

// 第三节点的短路时matlab程序
clc;
z=[0.23i        0.105i         inf          inf               inf             inf;
0.105i       -32.3939i   0.1452+0.2903i     inf          0.2178+0.4355i        inf;
inf       0.1452+0.2903i     -43.1965i    0.1271+0.2540i      inf             inf;
inf           inf        0.1271+0.2540i   -40.5022i     0.16335+0.3266i       inf;
inf       0.2178+0.4355i       inf        0.16335+0.3266i   -30.8547i          0.21i;  
inf           inf             inf         inf               0.21i            0.28i];
Y=zeros(6,6);
for(i=1:6)
for(j=1:6)
Y(i,i)=sum(1./z(i,:),2);
if j==i
Y(i,j)=Y(i,i);
else 
Y(i,j)=-1/z(i,j);
end
end
end
B=Y;
B(3,:)=[];
B(:,3)=[];
B

结果如下：

三四节点中间短路时matlab程序如下：

// 三四节点中间短路时matlab程序
clc;
z=[0.23i        0.105i         inf          inf               inf             inf;
0.105i       -32.3939i   0.1452+0.2903i     inf          0.2178+0.4355i        inf;
Inf       0.1452+0.2903i     -43.1965i    0.1271+0.2540i      inf             inf;
Inf           inf        0.1271+0.2540i   -40.5022i     0.16335+0.3266i       inf;
Inf       0.2178+0.4355i       inf        0.16335+0.3266i   -30.8547i          0.21i;  
Inf           inf             inf         inf               0.21i            0.28i];
Y=zeros(6,6);
for(i=1:6)
for(j=1:6)
Y(i,i)=sum(1./z(i,:),2);
if j==i
Y(i,j)=Y(i,i);
else 
Y(i,j)=-1/z(i,j);
end
end
end
B=Y;
B(3,4)=0;
B(4,3)=0;
B(3,3)=[B(3,3)+1./z(3,4)-0.0054i];
B(4,4)=[B(4,4)+1./z(3,4)-0.0054i];
B

运算结果如下图：