clear
n=40;%n=input（‘請(qǐng)輸入節(jié)點(diǎn)數(shù)‘）;
b=39;%b=input（‘請(qǐng)輸入支路數(shù)‘）;
%disp（‘請(qǐng)輸入支路阻抗和節(jié)點(diǎn)功率矩陣，第一列存支路號(hào)，第二列存首節(jié)點(diǎn)號(hào)，第三列存尾節(jié)點(diǎn)號(hào)，第四列存支路自阻抗第五列存尾節(jié)點(diǎn)給定功率‘）;
Sb=15;
Ub=23;
Zb=Ub^2/Sb;
Z=[1010.169/Zb+i*0.655/Zb0
2120.864/Zb+i*0.751/Zb0
3130.196/Zb+i*0.655/Zb0
4241.374/Zb+i*0.774/Zb0.522/Sb+i*0.174/Sb
5250.864/Zb+i*0.751/Zb0
6360.444/Zb+i*0.439/Zb0
7370.196/Zb+i*0.655/Zb0
8580.864/Zb+i*0.751/Zb0
9590.864/Zb+i*0.751/Zb0
106101.374/Zb+i*0.774/Zb1.08/Sb+i*0.36/Sb
116110.864/Zb+i*0.751/Zb0
127120.279/Zb+i*0.051/Zb0
137130.279/Zb+i*0.051/Zb0
148141.374/Zb+i*0.774/Zb0.36/Sb+i*0.12/Sb
158151.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
169161.374/Zb+i*0.774/Zb0.45/Sb+i*0.15/Sb
179170.864/Zb+i*0.751/Zb0
1811181.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
1911191.374/Zb+i*0.774/Zb0.675/Sb+i*0.225/Sb
2012201.374/Zb+i*0.774/Zb0.45/Sb+i*0.15/Sb
2112210.444/Zb+i*0.439/Zb0
2213220.444/Zb+i*0.439/Zb1.8/Sb+i*0.6/Sb
2313230.444/Zb+i*0.439/Zb0
2417241.374/Zb+i*0.774/Zb0.675/Sb+i*0.225/Sb
2517251.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
2621261.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
2721270.444/Zb+i*0.439/Zb0
2823281.374/Zb+i*0.774/Zb0.855/Sb+i*0.285/Sb
2923290.864/Zb+i*0.751/Zb0
3027301.374/Zb+i*0.774/Zb0.63/Sb+i*0.21/Sb
3127310.864/Zb+i*0.751/Zb0
3229321.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
3327331.374/Zb+i*0.774/Zb0.72/Sb+i*0.24/Sb
3431341.374/Zb+i*0.774/Zb0.765/Sb+i*0.225/Sb
3531350.864/Zb+i*0.751/Zb0
3635361.374/Zb+i*0.774/Zb0.54/Sb+i*0.18/Sb
3735370.864/Zb+i*0.751/Zb0
3837381.374/Zb+i*0.774/Zb0.675/Sb+i*0.225/Sb
3937391.374/Zb+i*0.774/Zb0.522/Sb+i*0.174/Sb];
k=0;
V=ones（n1）;
t=0;
%迭代開(kāi)始處
while?t%算節(jié)點(diǎn)注入電流
?x1=Z（b3）;x=x1-n;
???for?l=1:b
???????j=Z（l3）;
???????ua=V（j+11）;
?????I（j1）=conj（Z（j5）/ua）;
???end
%回推算支路電流
J=zeros（b1）;
l=b;
J（l）=J（l）+I（l）;
for?jj=1:b-1
????l=l-1;
?for?m=l+1:b
?????if?Z（m2）==Z（l3）?
?????J（l）=J（l）+J（m）;
????end?
end
J（l）=J（l）+I（l）;
end
%前推算節(jié)點(diǎn)電壓
???for?l=1:b
????j=Z（l3）+1;
????i=Z（l2）+1;
????V（j1）=V（i1）-Z（l4）*J（l1）;
??end
%收斂判定
t=0;
for?j=2:n
?SS=V（j1）*conj（I（j-11））;
dp=real（SS-Z（j-15））;
dq=imag（SS-Z（j-15））;
S（j-11）=SS;
ddp=abs（dp）;
ddq=abs（dq）;
L1=（ddp<0.000001）&（ddq<0.000001）;
F（j-11）=L1;
if?L1==1
?t=t+1;
end
end
?k=k+1;
end

%輸出結(jié)果與否

?disp（‘輸出直角坐標(biāo)各節(jié)點(diǎn)電壓‘）;
?disp（V）;
?disp（‘顯示迭代次數(shù)‘）;
?disp（k）;
disp（‘顯示收斂節(jié)點(diǎn)情況“1“表示收斂，“0“表示不收斂‘）;
?disp（F）;
?for?j=1:b
?????if?F（j1）==0
?????????disp（‘顯示不收斂節(jié)點(diǎn)號(hào)、計(jì)算功率‘）;
disp（j）;disp（S（j1））;
?????end
?end
for?j=1:n
Vm（j1）=abs（V（j1））;Va（j1）=angle（V（j1））;
end
disp（‘輸出各節(jié)點(diǎn)電壓幅值‘）;
disp（Vm）;
disp（‘輸出各節(jié)點(diǎn)電壓相角‘）;
disp（Va）
?%返回開(kāi)始處

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

?????文件???????2924??2008-06-19?09:15??cthd39-40.m

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

?????????????????2924????????????????????1