33节点的辐射型网络的研究

clear
n=33;%n=input（‘请输入节点数‘）;
b=32;%b=input（‘请输入支路数‘）;
%disp（‘请输入支路阻抗和节点功率矩阵，第一列存支路号，第二列存首节点号，第三列存尾节点号，第四列存支路阻抗，第五列存节点给定功率’）;
Sb=10;
Ub=12.66;
Zb=Ub^2/Sb;
Z=[ 1 0 10.0922/Zb+0.0470i/Zb0.100/Sb+0.060i/Sb
    2 1 20.4930/Zb+0.2511i/Zb0.090/Sb+0.040i/Sb
    3 2 30.3660/Zb+0.1864i/Zb0.120/Sb+0.080i/Sb
    4 3 40.3811/Zb+0.1941i/Zb0.060/Sb+0.030i/Sb
    5 4 50.8190/Zb+0.7070i/Zb0.060/Sb+0.020i/Sb
    6 5 60.1872/Zb+0.6188i/Zb0.200/Sb+0.100i/Sb
    7 6 70.7114/Zb+0.2351i/Zb0.200/Sb+0.100i/Sb
    8 7 81.0300/Zb+0.7400i/Zb0.060/Sb+0.020i/Sb
    9 8 91.0440/Zb+0.7400i/Zb0.060/Sb+0.020i/Sb
    10 9100.1966/Zb+0.0650i/Zb0.045/Sb+0.030i/Sb
    1110110.3744/Zb+0.1238i/Zb0.060/Sb+0.035i/Sb
    1211121.4680/Zb+1.1550i/Zb0.060/Sb+0.035i/Sb
    1312130.5416/Zb+0.7129i/Zb0.120/Sb+0.080i/Sb
    1413140.5910/Zb+0.5260i/Zb0.060/Sb+0.010i/Sb
    1514150.7463/Zb+0.5450i/Zb0.060/Sb+0.020i/Sb
    1615161.2890/Zb+1.7210i/Zb0.