IEEE33节点

%%IEEE33节点潮流计算程序%%
countnum=0;
BranchData=[
    1 2 0.0922 0.0470;
    2 3 0.4930 0.2511;
    3 4 0.3660 0.1864;
    4 5 0.3811 0.1941;
    5 6 0.8190 0.7070;
    6 7 0.1872 0.6188;
    7 8 0.7114 0.2351;
    8 9 1.0300 0.7400;
    9 10 1.0440 0.7400;
    10 11 0.1966 0.0650;
    11 12 0.3744 0.1238;
    12 13 1.4680 1.1550;
    13 14 0.5416 0.7129;
    14 15 0.5910 0.5260;
    15 16 0.7463 0.5450;
    16 17 1.2890 1.7210;
    17 18 0.7320 0.5
    2 19 0.1640 0.1565;
    19 20 1.5042 1.3554;
    20 21 0.4095 0.4784;
    21 22 0.7089 0.9373;
    3 23 0.4512 0.3083;
    23 24 0.8980 0.7091;
    24 25 0.8960 0.7011;
    6 26 0.2030 0.1034;
    26 27 0.2842 0.1447;
    27 28 1.0590 0.9337;
    28 29 0.8042 0.7006;
    29 30 0.5075 0.2585;
    30 31 0.9744 0.9630;
    31 32 0.3105 0.3619;
    32 33 0.3410 0.5302;];%支路阻抗
NodeData=[
    2 100.00 60.00;
    3 90.00 40.00;
    4 120.00 80.00;
    5 60.00 30.00;
    6 60.00 20.00;
    7 200.00 100.00;
    8 200.00 100.00;
    9 60.00 20.00;
    10 60.00 20.00;
    11 45.00 30.00;
    12 60.00 35.00;
    13 60.00 35.00;
    14 120.00 80.00;
    15 300.00 200.00;
    16 60.00 20.00;
    17 60.00 20.00;
    18 90.00 40.00;
    19 90.00 40.00;
    20 90.00 40.00;
    21 90.00 40.00;
    22 90.00 40.00;
    23 90.00 50.00;
    24 420.00 200.00;
    25 300.00 200.00;
    26 60.00 25.00;
    27 60.00 25.00;
    28 60.00 20.00;
    29 120.00 70.00;
    30 200.00 600.00;
    31 150.00 70.00;
    32 210.00 100.00;
    33 300.00 40.00;];%节点，负荷
pv=[15;25];pvg=2;hk=0;Qx=[-100 -500;-100 -500];pvv=[1;1];%PV节点
pii=29;piig=1;piii=0.0633;Qk=0;                         %PI节点
pqv=33;pqvg=1;pqx=[0.8 0.5 50];                        %PQ（V）节点
dv=0;
UB=12.66;%电压基准kV
SB=10;%功率基准MVA
ZB=UB^2/SB;%阻抗基准0hm
BranchData（:[34]）=BranchData（:[34]）/ZB;%阻抗标幺化
NodeData（:[23]）=NodeData（:[23]）/SB/1000;%功率标幺化
Qx=Qx/SB/1000;
NN=33;  %节点数
A0=zeros（NN）;
for n=1:NN-1
    A0（BranchData（n1）BranchData（n2））=1;
end %形成A0

X=zeros（pvg）;
for n=1:pvg  %形成PV节点的电抗矩阵
    s=pv（n）;
    while s>1
        X（nn）=X（nn）+BranchData（s-14）;%求自导纳
        s=BranchData（s-11）;
    end;
end;
for m=1:pvg
    for