系統(tǒng)接線示于題圖1，已知各元件參數(shù)見圖上標(biāo)注：
發(fā)電機(jī)G-1：＝120MV·A，＝0.23；G-2：＝60MV·A，＝0.14。
變壓器T-1：＝120MV·A，＝10.5％；T-2：＝60MV·A，＝10.5％。
線路參數(shù)： , 。線路長(zhǎng)度L-1：120km，L-2：80km，L-3：70km，L-4: 90km。取＝120MV·A，，試求標(biāo)幺制下的節(jié)點(diǎn)導(dǎo)納矩陣。另試就下列兩種情況分別修改節(jié)點(diǎn)導(dǎo)納矩陣：（1）節(jié)點(diǎn)3發(fā)生三相短路；（2）線路L-3中間發(fā)生三相短路。

matlab程序如下：

// 穩(wěn)態(tài)時(shí)導(dǎo)納矩陣計(jì)算
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

結(jié)果如下：

第三個(gè)節(jié)點(diǎn)短路時(shí)matlab程序如下：

// 第三節(jié)點(diǎn)的短路時(shí)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

結(jié)果如下：

三四節(jié)點(diǎn)中間短路時(shí)matlab程序如下：

// 三四節(jié)點(diǎn)中間短路時(shí)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

運(yùn)算結(jié)果如下圖：
文章來源地址http://www.zghlxwxcb.cn/news/detail-511599.html

到了這里，關(guān)于電力系統(tǒng)導(dǎo)納矩陣與短路計(jì)算matlab程序的文章就介紹完了。如果您還想了解更多內(nèi)容，請(qǐng)?jiān)谟疑辖撬阉鱐OY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章，希望大家以后多多支持TOY模板網(wǎng)！