clear;clc;
format?compact
%?n=input（‘請輸入節點數：n=‘）;
%?pr=input（‘請輸入精度誤差：pr=‘）;
%?Z=input（‘請輸入由各節點間阻抗形成的距陣：Z=‘）;
%?B=input（‘請輸入由各節點直接對地的導納距陣：B=‘）;
%?K=input（‘請輸入由節點間變壓器變比值：K=‘）;
%?N=input（‘請輸入當變壓器兩端節點對地等效導納值確定條件：N=‘）;
%數據輸入環節%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
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）%（經變壓器變換后對地的等效導納）;
A=zeros（n）;D=zeros（1n）%（經變壓器變換后總的對地導納）;
H=zeros（1n）%（經變壓器變換后節點間互導納）;
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???????%求導納距陣
disp（‘導距Y為:‘）;
disp（Y）

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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

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

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

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