程序简化代码如下:
NIND=40; % 定义个体数目
MAXGEN=25; % 定义最大遗传代数
PRECI=20; % 编码长度
GGAP=0.9; % 代沟
FieldD=[20;-1;2;1;0;1;1]; % 区域描述器
Chrom=crtbp(NIND, PRECI); % 产生初始种群
gen=0; % 代计数器
variable=bs2rv(Chrom, FieldD); % 初始种群的十进制转换
ObjV=variable.*sin(10*pi*variable)+2.0; % 计算初始种群目标函数值
while gen<MAXGEN+1
FitnV=ranking(-ObjV); % 分配适应度值
SelCh=select('sus', Chrom, FitnV, GGAP); % 选择
SelCh=recombin('xovsp', SelCh, 0.7); % 重组
SelCh=mut(SelCh); % 变异
variable=bs2rv(SelCh, FieldD); % 子代个体的十进制转换
ObjVSel=variable.*sin(10*pi*variable)+2.0; % 计算子代的目标函数值
[Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); % 重插入子代到种
群
Y=max(ObjV) % 当前种群最优解
gen=gen+1; % 代计数器增加
end
第 2 页 / 共 74 页