遗传算法计算无功优化

tic                   %%%%%%%%%%%%%%%
clc;clear all;
%相关原始数据格式说明如下：
%n——节点个数；n1——支路条数；isb——平衡节点号；H——PQ节点个数（为后面形成PVU存储PV节点初始电压用）；pr——误差精度。
%B1——支路参数矩阵，其中第一列和第二列是起始节点编号和终点节点编号第三列、第四列、第五列、第六列分别为：支路电阻、电抗、变压器变比、电纳。（不考虑电导）
%B2——节点参数矩阵，其中第一列和第二列为节点编号和节点类型；第三列到第六列分别为：注入有功、注入无功、电压幅值、电压相位。
%节点类型分类如下：“0”为平衡节点“1”为PQ，“2”为PV节点；“3”为PQ（V）节点，“4”为PI节点。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%数据输入：基准电压10KV，基准功率10MW %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n=33 ;       
n1=32;
isb=1;
H=32;                %%%%%%%%%%%%%18节点加DG    PQV处理
pr=0.0001;          
B1=[1 2 0.00922 0.0047i 1 0;
    2 3 0.00493 0.02511i 1 0;
    3 4 0.0366 0.01864i 1 0;
    4 5 0.03811 0.01941i 1 0;
    5 6 0.0819 0.0707i 1 0;
    6 7 0.01872 0.06188i 1 0;
    7 8 0.07114 0.02351i 1 0;
    8 9 0.103 0.074i 1 0;
    9 10 0.1044 0.074i 1 0;
    10 11 0.01966 0.0065i 1 0;
    11 12 0.03744 0.01238i 1 0;
    12 13 0.1468 0.1155i 1 0;
    13 14 0.05416 0.07129i 1 0;
    14 15 0.05910 0.0526i 1 0;
    15 16 0.07463 0.05450i 1 0;
    16 17 0.1289 0.1721i 1 0;
    17 18 0.0732 0.0574i 1 0;
    2 19 0.0164 0.01565i 1 0;
    19 20 0.15042 0.13554i 1 0;
    20 21 0.04095 0.04784i 1 0;
    21 22 0.07089 0.09373i 1 0;
    3 23 0.04512 0.03083i 1 0;
    23 24 0.08980 0.07091i 1 0;
    24 25 0.08960 0.07011i 1 0;
    6 26 0.0203 0.01034i 1 0;
    26 27 0.02842 0.01447i 1 0;
    27 28 0.1059 0.09337i 1 0;
    28 29 0.08042 0.07006i 1 0;
    29 30 0.05075 0.02585i 1 0;
    30 31 0.09744 0.0963i 1 0;
    31 32 0.03105 0.03619i 1 0;
    32 33 0.03410 0.05302i 1 0];
B2=[1 0 0 0 1.05 0;
    2 1 -0.01 -0.006 1 0;
    3 1 -0.009 -0.004 1 0;
    4 1 -0.012 -0.008 1 0;
    5 1 -0.006 -0.003 1 0;
    6 1 -0.006 -0.002 1 0;
    7 1 -0.02 -0.01 1 0;
    8 1 -0.02 -0.01 1 0;
    9 1 -0.006 -0.002 1 0;
    10 1 -0.006 -0.0035 1 0;
    11 1 -0.0045 -0.003 1 0;
    12 1 -0.006 -0.0035 1 0;
    13 1 -0.006 -0.0035 1 0;
    14 1 -0.012 -0.008 1 0;
    15 1 -0.006 -0.001 1 0;
    16 1 -0.006 -0.002 1 0;
    17 1 -0.006 -0.002 1 0;
    18 3 -0.009+0.1 -0.004 1 0;
    19 1 -0.009 -0.004 1 0;
    20 1 -0.009 -0.004 1 0;
    21 1 -0.009 -0.004 1 0;
    22 1 -0.009 -0.004 1 0;
    23 1 -0.009 -0.005 1 0;
    24 1 -0.042 -0.02 1 0;
    25 1 -0.042 -0.02 1 0;
    26 1 -0.006 -0.0025 1 0;
    27 1 -0.006 -0.0025 1 0;
    28 1 -0.006 -0.002 1 0;
    29 1 -0.012 -0.007 1 0;
    30 1 -0.02 -0.06 1 0;
    31 1 -0.015 -0.007 1 0;
    32 1 -0.021 -0.01 1 0;
    33 4 -0.006 -0.004 1 0];
T=60;%仿真代数               遗传算法
M=40;% 群体规模
pm=0.1;pc=0.8;%交叉变异概率
Ugmax=1.06;Ugmin=1.0;%参数取值范围
Timax=1.1;Timin=0.9;
Qcmax=0.05;Qcmin=-0.05;
Vimax=1.07;Vimin=0.93;
Z1=1;%Z2=1;%%惩罚因子%%%%%%%%%%%%%%%&&&&&&&&&&&&&&&&&&&&&&&*********************具体数值待定
%QGmax=;QGmin=;
Long=5;%单个参数字串长度，总编码长度3L
bestv=-inf;
bval=round（rand（M3*Long））;%初始种群 round 四舍五入取整函数  编码后的初始种群
%%%%%%%%%%%%%bestv=-inf%最优适应度初值%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%