clear?all
clc
r（1410）=0.01938;x（1410）=0.05917;r（1011）=0.04699;x（1011）=0.19797;r（101）=0.05811;x（101）=0.17632;r（142）=0.05403;x（142）=0.22304;r（102）=0.05695;x（102）=0.17388;
r（1014）=0.01938;x（1014）=0.05917;r（1110）=0.04699;x（1110）=0.19797;r（110）=0.05811;x（110）=0.17632;r（214）=0.05403;x（214）=0.22304;r（210）=0.05695;x（210）=0.17388;
r（111）=0.06701;x（111）=0.17103;r（12）=0.01335;x（12）=0.04211;r（212）=0.00000;x（212）=0.23488;r（13）=0.00000;x（13）=0.20452;r（313）=0.00000;x（313）=0.17615;
r（111）=0.06701;x（111）=0.17103;r（21）=0.01335;x（21）=0.04211;r（122）=0.00000;x（122）=0.23488;r（31）=0.00000;x（31）=0.20452;r（133）=0.00000;x（133）=0.17615;
r（14）=0.00000;x（14）=0.53894;r（34）=0.00000;x（34）=0.11001;r（45）=0.03181;x（45）=0.08450;r（126）=0.09498;x（126）=0.19890;r（127）=0.12291;x（127）=0.25581;
r（41）=0.00000;x（41）=0.53894;r（43）=0.00000;x（43）=0.11001;r（54）=0.03181;x（54）=0.08450;r（612）=0.09498;x（612）=0.19890;r（712）=0.12291;x（712）=0.25581;
r（128）=0.06615;x（128）=0.13027;r（49）=0.12711;x（49）=0.27038;r（56）=0.08205;x（56）=0.19207;r（78）=0.22092;x（78）=0.19988;r（89）=0.17093;x（89）=0.34802;
r（812）=0.06615;x（812）=0.13027;r（94）=0.12711;x（94）=0.27038;r（65）=0.08205;x（65）=0.19207;r（87）=0.22092;x（87）=0.19988;r（98）=0.17093;x（98）=0.34802;
?
y（1410）=0.02640;y（1011）=0.02190;y（101）=0.01870;y（142）=0.02460;y（102）=0.01700;y（111）=0.01730;y（12）=0.00640;y（212）=-0.31063;y（13）=-0.10999;y（14）=-0.05936;
y（1014）=0.02640;y（1110）=0.02190;y（110）=0.01870;y（214）=0.02460;y（210）=0.01700;y（111）=0.01730;y（21）=0.00640;y（122）=0.28951;y（31）=0.10757;y（41）=0.05752;
?
for?m=1:14
????t（m）=0;
????for?n=1:14
????????t（m）=t（m）+y（mn）;
????end
????if?m==4
????????t（m）=t（m）+0.190;
????end
end
for?m=1:14
????d（m）=0;
????for?n=1:14
????????if?m==n
????????????
????????elseif?（r（mn）==0）&（x（mn）==0）
???????????d（m）=d（m）+0;?
????????else
???????????d（m）=d（m）+1/（r（mn）+j*x（mn））;????
????????end
????end
end
for?m=1:14
????for?n=1:14
????????if?m==n
???????????Y（mn）=j*t（m）+d（m）;
????????elseif?（r（mn）==0）&（x（mn）==0）
???????????Y（mn）=0;
????????else
???????????Y（mn）=-1/（r（mn）+j*x（mn））;
????????end
????end
end
Y;
G=real（Y）;B=imag（Y）;
%給定節點電壓初值和循環次數k=0
delt（10）=0;delt（11）=0;u（1）=1.0;delt（1）=0;u（2）=1.0;delt（2）=0;delt（12）=0;u（3）=1.0;delt（3）=0;delt（13）=0;u（4）=1.0;delt（4）=0;u（5）=1.0;delt（5）=0;u（6）=1.0;delt（6）=0;u（7）=1.0;delt（7）=0;u（8）=1.0;delt（8）=0;u（9）=1.0;delt（9）=0;
p（10）=0.183;p（11）=-0.942;p（1）=-0.478;q（1）=0.039;p（2）=-0.076;q（2）=-0.016;p（12）=-0.112;p（3）=0;q（3）=0;p（13）=0;p（4）=-0.295;p（5）=-0.09;q（5）=-0.058;p（6）=-0.035;q（6）=-0.018;p（7）=-0.061;q（7）=-0.016;p（8）=-0.135;q（8）=-0.058;p（9）=-0.149;q（9）=-0.05;
k=0;precision=1;N1=13;
%求節點功率的不平衡量??P和??Q
while?precision>0.00001
u（10）=1.0450;u（11）=1.01;u（12）=1.07;u（13）=1.09;u（14）=1.06;delt（14）=0;q（4）=-0.166+（u（4））^2*0.190;
?????for?m=1:N1?????
?????????if?m<=9
?????????????for?n=1:N1+1
pt（n）=u（m）*u（n）*（G（mn）*cos（delt（m）-delt（n））+B（mn）*sin（delt（m）-delt（n）））;