clear;clc;
format?compact
%?n=input（‘請(qǐng)輸入節(jié)點(diǎn)數(shù)：n=‘）;
%?pr=input（‘請(qǐng)輸入精度誤差：pr=‘）;
%?Z=input（‘請(qǐng)輸入由各節(jié)點(diǎn)間阻抗形成的距陣：Z=‘）;
%?B=input（‘請(qǐng)輸入由各節(jié)點(diǎn)直接對(duì)地的導(dǎo)納距陣：B=‘）;
%?K=input（‘請(qǐng)輸入由節(jié)點(diǎn)間變壓器變比值：K=‘）;
%?N=input（‘請(qǐng)輸入當(dāng)變壓器兩端節(jié)點(diǎn)對(duì)地等效導(dǎo)納值確定條件：N=‘）;
%數(shù)據(jù)輸入環(huán)節(jié)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n=5;
pr=0.00001;
Z=[0?0.04+0.25i?0.1+0.35i?0?0;
???0.04+0.25i?0?0.08+0.30i?0.015i?0;
???0.1+0.35i?0.08+0.30i?0?0?0.03i;
???0?0.015i?0?0?0;
???0?0?0.03i?0?0];
B=[0.25i?0.5i?0.25i?0?0];
K=[1?1?1?1?1;
???1?1?1?1.05?1;
???1?1?1?1?1.05;
???1?1.05?1?1?1;
???1?1?1.05?1?1;];
N=[0?1?1?0?0];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Y=zeros（n）;C=zeros（n）%（經(jīng)變壓器變換后對(duì)地的等效導(dǎo)納）;
A=zeros（n）;D=zeros（1n）%（經(jīng)變壓器變換后總的對(duì)地導(dǎo)納）;
H=zeros（1n）%（經(jīng)變壓器變換后節(jié)點(diǎn)間互導(dǎo)納）;
for?i=1:n
????for?j=1:n
????????if?Z（ij）~=0
?????????A（ij）=1/Z（ij）
????????end
????????if?Z（ij）==0
??????????A（ij）=0
????????end
????end
end
for?i=1:n
????for?j=1:n
????????H（ij）=A（ij）/K（ij）
????????if?K（ij）==1
???????????C（ij）=0
????????end
????????if?K（ij）~=1
???????????if?N（i）==1
??????????????C（ij）=（1-K（ij））*H（ij）/K（ij）
???????????end
????????????if?N（i）==0
????????????????C（ij）=（K（ij）-1）*H（ij）
???????????end
????????end
????????Q=C（i:）?%（取C矩陣的某一行元素）
????????D（i）=sum（sum（Q））+B（i）
????end
end
for?i=1:n
????for?j=1:n
????????if?i~=j
????????????Y（ij）=-H（ij）
????????end
????????if?i==j
????????????E=H（i:）?%（取H矩陣的某一行元素）
????????????Y（ij）=sum（sum（E））+D（i）
????????end
????end
end???????%求導(dǎo)納距陣
disp（‘導(dǎo)距Y為:‘）;
disp（Y）

????????????
?
???????

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????1715??2010-04-08?23:22??daonaijuzhentongyongchengxu.m

-----------?---------??----------?-----??----

?????????????????1715????????????????????1