电力系统导纳矩阵与短路计算matlab程序-Toy模板网

这篇具有很好参考价值的文章主要介绍了电力系统导纳矩阵与短路计算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程序
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程序如下：

// 第三节点的短路时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程序如下：

// 三四节点中间短路时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