ieee14节点的潮流计算

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）））;