|
觉得有几处错误。
for j=1:Numv
pop1=ones(px,m);
pop2=pop(:,m*(j-1)+1:m*j); %取出相应变量对应的二进制编码段
for i=1:2:px-1
if(rand<pc)
cpoint=round(rand*(m-1)); %cpoint为交叉点
pop1(i,:)=[pop2(i,1:cpoint) pop2(i+1,cpoint+1:m)];
pop1(i+1,:)=[pop2(i+1,1:cpoint) pop2(i,cpoint+1:m)];
else
pop1(i,:)=pop2(i,1:m);
pop1(i+1,:)=pop2(i+1,1:m);
end
end
newpop(:,m*(j-1)+1:m*j)=pop1; %将交叉后的一个参数的编码放入新种群中
end
觉得红色字体这两处互换才对
%Name:genmain05.m
clear
clf
popsize=20; %群体大小
chromlength=10; %字符串长度(个体长度)
pc=0.6; %交叉概率
pm=0.001; %变异概率
pop=initpop(popsize,chromlength); %随机产生初始群体
for i=1:20 %20为迭代次数
[objvalue]=calobjvalue(pop); %计算目标函数
fitvalue=calfitvalue(objvalue); %计算群体中每个个体的适应度
[newpop]=selection(pop,fitvalue); %复制
[newpop]=crossover(pop,pc); %交叉
[newpop]=mutation(pop,pc); %变异
———这两处应该换成newpop吧? |
|