clc;clear all;
close all;
errcode = loadlibrary('epanet2.dll','epanet2.h');%加载EPANET文件
errcode = calllib('epanet2','ENopen','pump.inp','pump.rpt','');%打开inp文件
jc = [2 3 4 5 6];%节点
nodevalue = 0;
nVar = length(jc); % x的长度
EN_FLOW = 0;
time = 0;%初始化工况时间,可以是任意值
number = 0;%number为水力步长进行到的阶段数
tstep = 1;%初始化水力分析的步数,可以是任意非零值
totaltime = 86400;%模拟时长
errcode = calllib('epanet2','ENsettimeparam',0,totaltime);%设置模拟历时
timestep = 3600;%时间步长
nPop = 10; % 种群规模大小
maxIt = 500; % 最大迭代次数
nPc = 0.8; % 子代规模的比例
nC = round(nPop * nPc/2)*2;% 子代规模的大小 round()得出最接近的整数,同时进行偶数化操作
nMu = 0.01; % 变异概率
% 用结构体方式来存储变量,使得x,y有一个绑定的关系
template.x = []; % 模板
template.y = []; % 模板
% repmat函数主要是用于快速的产生一个大的矩阵
%%
遗传算法优化
于 2023-09-20 09:37:51 首次发布