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

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