ga.m

请大家看一看, 我编的这个用遗传算法求 f(x)=xsin(10pi*x)+2.0 x为-1到2区间的值

%遗传算法主程序
popsize=100;        %初始化时种群的数目;
chromlength=22;     % 染色体的长度;
pc=0.42;           %交叉概率;
pm=0.01;           %变异概率;
%X=ones(popsize,chromlength);
%初始化种群
inipop=inipop1(popsize,chromlength);
decData1=bin2dec(inipop);                     %把二进制数转换成十进制数;   
X=decData1;
a=-1;
b=2;
decdata=bin2decData(X,a,b,chromlength);      %把十进制数转换到对应域,并输出;
decdata;
fitness=decdata.*sin(decdata.*10*pi)+2.0;   %计算适应度;
fitness;
newpops=selection1(inipop,fitness)
newpopc=crossover1(newpops,pc);              %交叉
newpopc;
newpopm=mutation1(newpopc,pm);              %变异
newpopm
for i=1:300
    decData1=bin2Dec(newpopm);  
  i
X=decData1;
a=-1;
b=2;
decdata=bin2decData(X,a,b,chromlength);
decdata;
fitness=decdata.*sin(decdata.*10*pi)+2.0;
newpops=selection1(newpopm,fitness);
newpopc=crossover1(newpops,pc);
newpopc;
newpopm=mutation1(newpopc,pm);
newpopm;
%fitness
m=max(fitness)
end;