24小时热门版块排行榜

返回列表

witzwq

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 696.3
帖子: 210
在线: 46小时
虫号: 2235575
注册: 2013-01-11
性别: MM
专业: 计算机软件

[求助] 遗传算法matlab程序求解

figure(1);
fplot('variable.*sin(10*pi*variable)+2.0',[-1,2]);%画出函数曲线
%定义遗传算法参数
NIND=40;%个体数目（Number of individuals）
MAXGEN=25;%最大遗传代数（Maximum number.of generations)
PRECI=20;%变量的二进制位数(Precision of variables)
GGAP=0.9;%代沟(Generation gap)
trace=zeros(2,MAXGEN);%寻优结果的初始值
FieldD=[20;-1;2;1;0;1;1];%区域描述器(Build field descriptor)格雷编码，算术刻度，包含边界
Chrom=crtbp(NIND,PRECI);%初始种群
gen=0;%代计数器
variable=bs2rv(Chrom,FieldD);%计算初始种群的十进制转换
ObjV=variable.*sin(10*pi*variable)+2.0;%计算目标函数值
while gen FitnV=ranking(-ObjV);%分配适应度值(Assign fitness values)
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);%重插入子代的新种群
gen=gen+1;%代计数器增加
%输出最优解及其序号,并在目标函数图像中标出,Y为最优解,I为种群的序号
[Y,I]=max(ObjV);
hold on;
plot(variable(I),Y,'bo');%b是blue，o是圆圈
trace(1,gen)=max(ObjV);%遗传算法性能跟踪
trace(2,gen)=sum(ObjV)/length(ObjV);
end
variable=bs2rv(Chrom,FieldD);%最优个体的十进制转换
hold on ;
grid on;%显示网格线
plot(variable',ObjV','b*');
figure(2);
plot(trace(1,