基于遗传算法的配电网重构研究（Matlab代码实现）

function [U,ploss]=powerflow()
n=69;
b=68;
Sb=10;       %基准功率
Ub=12.66;    %基准电压
Zb=Ub^2/Sb;  %基准阻抗
s1=13; s2=69; s3=55; s4=18; s5=61;
Q1=0; Q2=0; Q3=0;
%输入系统参数
Z=[1, 1, 2,0.005/Zb+0.0012i/Zb,0
    2, 2, 3,0.005/Zb+0.0012i/Zb,0
    3, 3, 4,0.0015/Zb+0.0036i/Zb,0
    4, 4, 5,0.0251/Zb+0.02941i/Zb,0
    5, 5, 6,0.366/Zb+0.1869i/Zb,0.0026/Sb+0.0022i/Sb
    6, 6, 7,0.3811/Zb+0.1964i/Zb,0.0404/Sb+0.03i/Sb
    7, 7, 8,0.0922/Zb+0.047i/Zb,0.075/Sb+0.054i/Sb
    8, 8, 9,0.0493/Zb+0.0251i/Zb,0.030/Sb+0.022i/Sb
    9, 9, 10,0.8090/Zb+0.2707i/Zb,0.028/Sb+0.019i/Sb
    10, 10,11,0.1872/Zb+0.0691i/Zb,0.145/Sb+0.104i/Sb
    11,11,12,0.7114/Zb+0.2351i/Zb,0.145/Sb+0.104i/Sb
    12,12,13,1.03/Zb+0.34i/Zb,0.008/Sb+0.0055i/Sb
    13,13,14,1.044/Zb+0.345i/Zb,0.008/Sb+0.0055i/Sb
    14,14,15,1.058/Zb+0.3496i/Zb,0
    15,15,16,0.1966/Zb+0.065i/Zb,0.0455/Sb+0.030i/Sb
    16,16,17,0.3744/Zb+0.1238i/Zb,0.060/Sb+0.035i/Sb
    17,17,18,0.0047/Zb+0.0016i/Zb,0.060/Sb+0.035i/Sb
    18,18,19,0.3276/Zb+0.1083i/Zb,0
    19,19,20,0.2106/Zb+0.0696i/Zb,0.001/Sb+0.0006i/Sb
    20,20,21,0.3416/Zb+0.1129i/Zb,0.114/Sb+0.081i/Sb
    21,21,22,0.014/Zb+0.0046i/Zb,0.0053/Sb+0.0035i/Sb
    22, 22,23,0.1591/Zb+0.0526i/Zb,0
    23,23,24,0.3463/Zb+0.11451i/Zb,0.028/Sb+0.020i/Sb
    24,24,25,0.7488/Zb+0.2745i/Zb,0
    25, 25,26,0.3089/Zb+0.1021i/Zb,0.014/Sb+0.010i/Sb
    26,26,27,0.1732/Zb+0.5728i/Zb,0.014/Sb+0.010i/Sb
    27,3,20,0.0044/Zb+0.0108i/Zb,0.026/Sb+0.0186i/Sb
    28,28,29,0.064/Zb+0.1565i/Zb,0.026/Sb+0.0186i/Sb
    29,29,30,0.3978/Zb+0.1315i/Zb,0
    30,30,31,0.0702/Zb+0.0232i/Zb,0
    31,32,33,0.3510/Zb+0.1160i/Zb,0
    32,32,33,0.8390/Zb+0.2816i/Zb,0.014/Sb+0.010i/Sb
    33, 33,34,1.7080/Zb+0.5645i/Zb,0.0194/Sb+0.014i/Sb
    34, 34,35,1.4740/Zb+0.4673i/Zb,0.0060/Sb+0.004i/Sb
    35,3,59,0.0044/Zb+0.0108i/Zb,0.0260/Sb+0.01855i/Sb
    36,59,60,0.064/Zb+0.15650i/Zb,0.0260/Sb+0.01855i/Sb
    37,60,61,0.10530/Zb+0.123i/Zb,0
    38,61,62,0.0304/Zb+0.355i/Zb,0.024/Sb+0.017i/Sb
    39,62,63,0.018/Zb+0.021i/Zb,0.024/Sb+0.017i/Sb
    40,63,64,0.7283/Zb+0.8509i/Zb,0.0012/Sb+0.001i/Sb
    41,64,65,0.310/Zb+0.3623i/Zb,0
    42,65,66,0.041/Zb+0.0478i/Zb,0.006/Sb+0.0043i/Sb
    43,66,67,0.0092/Zb+0.0116i/Zb,0
    44,67,68,0.1089/Zb+0.1373i/Zb,0.03922/Sb+0.0263i/Sb
    45,68,69,0.0009/Zb+0.0012i/Zb,0.03922/Sb+0.0263i/Sb
    46,4,36,0.0034/Zb+0.0034i/Zb,0
    47,36,37,0.0851/Zb+0.2033i/Zb,0.079/Sb+0.0564i/Sb
    48,37,38,0.2898/Zb+0.7091i/Zb,0.3847/Sb+0.2745i/Sb
    49,38,39,0.0822/Zb+0.2011i/Zb,0.3847/Sb+0.2745i/Sb
    50,8,40,0.0928/Zb+0.0473i/Zb,0.0405/Sb+0.0283i/Sb
    51,40,41,0.3391/Zb+0.1114i/Zb,0.0036/Sb+0.0027i/Sb
    52,9,42,0.1740/Zb+0.0886i/Zb,0.00435/Sb+0.0035i/Sb
    53,42,43,0.2030/Zb+0.1034i/Zb,0.0264/Sb+0.019i/Sb
    54,43,44,0.2842/Zb+0.1447i/Zb,0.024/Sb+0.0172i/Sb
    55,44,45,0.2813/Zb+0.1433i/Zb,0
    56,45,46,1.59/Zb+0.5337i/Zb,0
    57,46,47,0.7837/Zb+0.2630i/Zb,0
    58,47,48,0.3042/Zb+0.1006i/Zb,0.1/Sb+0.072i/Sb
    59,48,49,0.3861/Zb+0.1172i/Zb,0
    60,49,50,0.5075/Zb+0.2585i/Zb,1.244/Sb+0.888i/Sb
    61,50,51,0.0974/Zb+0.0496i/Zb,0.032/Sb+0.023i/Sb
    62,51,52,0.1450/Zb+0.0738i/Zb,0
    63,52,53,0.7105/Zb+0.3619i/Zb,0.227/Sb+0.162i/Sb
    64,53,54,1.041/Zb+0.5302i/Zb,0.059/Sb+0.042i/Sb
    65,11,55,0.2012/Zb+0.0611i/Zb,0.018/Sb+0.013i/Sb
    66,55,56,0.0047/Zb+0.0014i/Zb,0.018/Sb+0.013i/Sb
    67,12,57,0.7394/Zb+0.2444i/Zb,0.028/Sb+0.020i/Sb
    68,57,58,0.0047/Zb+0.0016i/Zb,0.028/Sb+0.020i/Sb
    69,11,66,0.5/Zb+0.5/Zb,0.018/Sb+0.013i/Sb
    70,13,20,0.5/Zb+0.5/Zb,0.114/Sb+0.081i/Sb
    71,15,69,1/Zb+1/Zb,0.018/Sb+0.013i/Sb
    72,27,54,1/Zb+1/Zb,0.024/Sb+0.0172i/Sb
    73,39,48,2/Zb+2/Zb,0.3847/Sb+0.2745i/Sb
    ]; 
for i=1:73
    if Z(i,1)==s1
        Z(i,5)=0;
    end
end
for i=1:73
    if Z(i,1)==s2
        Z(i,5)=0;
    end
end
for i=1:73
    if Z(i,1)==s3
        Z(i,5)=0;
    end
end
for i=1:73
    if Z(i,1)==s4
        Z(i,5)=0;
    end
end
for i=1:73
    if Z(i,1)==s5
        Z(i,5)=0;
    end
end
%计算潮流
v=ones(1,n);   %初始化节点电压
k=0;
L1=0;
while  k<10&&L1<1     %收敛条件
    v_before=v;
    %计算节点注入电流，b=32为系统支路数
    I=zeros(b,1);   %生成一个b行1列的零矩阵
    for h=1:b
        I(h)=conj(0.5*Z(h,5)/v(1,h+1)); %求节点注入电流，conj为求共轭
    end
    J=zeros(b,1);
    J=J+I;    %支路电流，J矩阵中存储的是各节点注入电流，J有32行1列
    for h=b-1:-1:1
        nn=Z(h,3);   %取系统中末端节点编号
        for j=b:-1:1
            if Z(j,2)==nn    %说明他们之间首末端点是相连的
                J(h)=J(h)+J(j);
            end
        end
    end
   %回推节点电压
    for h=1:b
        m=Z(h,2);mm=Z(h,3);
        v(1,mm+1)=v(1,m+1)-Z(h,4)*J(h);
    end
    %收敛条件
    for h=1:b
        dat_v(h)=abs(v(h))-abs(v_before(h));
    end
    if max(abs(dat_v))<1.0e-4
        L1=1;    %不能用dat_v（i）就指最后的元素第i个
    end
    k=k+1;
end
%计算网损
ws=zeros(b,1);
for i=1:1:b;
    ws(i,1)=J(i,1)*conj(J(i,1))*real(Z(i,4))*Sb;
end
ploss=zeros(n,1);
ploss=sum(ws);
U=zeros(n,1);
for i=1:n;
    U(i)=real(v(i));
end
ploss
end