clc;
clear;
%配電網(wǎng)編號(hào)從零開(kāi)始，電源點(diǎn)為0，電源點(diǎn)為平衡點(diǎn)
z=input（‘請(qǐng)輸入配電網(wǎng)參數(shù)z=‘）;????%第一列為支路首節(jié)點(diǎn)號(hào)，第二列為支路末節(jié)點(diǎn)號(hào)，第三列為支路參數(shù)，第四列為末節(jié)點(diǎn)負(fù)荷第五列為末節(jié)點(diǎn)類型（1為PQ節(jié)點(diǎn)，2為PV節(jié)點(diǎn)）
n=input（‘請(qǐng)輸入配電網(wǎng)節(jié)點(diǎn)數(shù)n=‘）;
e=input（‘請(qǐng)輸入收斂精度e=‘）;
Sb=1;
Vb=10.5;
Zb=Vb^2/Sb;
Z（:3）=Z（:3）/Zb;
Z（:4）=Z（:4）/Sb/1000;???????????%修改配電網(wǎng)參數(shù)，使之統(tǒng)一至標(biāo)幺值
V=ones（n-11）;???????????????????%為各節(jié)點(diǎn)賦電壓初值（電源點(diǎn)除外）

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%開(kāi)始計(jì)算
for?i=1:20???????????????????????%i為迭代次數(shù)，最大迭代次數(shù)為20，超過(guò)20則收斂性較差
????I=（z（:4）.‘）*（diag（1./V））;???%求各節(jié)點(diǎn)注入電流
????I=I‘;
????I1=I;
???for?k=n-1:-1:1????????????????%前推支路電流及功率分布
??????A=（find（z（:1）==k））;???????%判斷矩陣Z第1列是否存在值為k的元素，并把行值賦給矩陣A，亦即找出網(wǎng)絡(luò)的末節(jié)點(diǎn)
??????if?size（A）==0??????????????%k為末節(jié)點(diǎn)
?????????I（k1）;
??????else?
?????????I（k1）=sum（I（A1））+I（k1）;%利用KCL求支路電流