遗传算法——matlab代码解析

% 求下列函数的最大值 % % f(x)=10*sin(5x)+7*cos(4x) x∈[0,10] % 将 x 的值用一个10位的二值形式表示为二值问题，一个10位的二值数提供的分辨率是每为 (10-0)/(2^10-1)≈0.01 % 将变量域 [0,10] 离散化为二值域 [0,1023], x=0+10*b/1023, 其中 b 是 [0,1023] 中的一个二值数。 % 2.8 主程序 %遗传算法主程序 %Name:main.m tic; clear,clc; 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(newpop,pc); %交叉 [newpop]=mutation(newpop,pm); %变异 [bestindividual,bestfit]=best(pop,fitvalue); %求出群体中适应值最大的个体及其适应值 y(i)=bestfit; pop_best=bestindividual; x(i)=decodechrom(pop_best,1,chromlength)*10/1023; %转化为十进制解 pop=newpop; %产生子代 end syms t; func = @(t)10*sin(5*t)+7*cos(4*t); %要求的函数 fplot(func,[0 10]); hold on plot(x,y,'r*'); hold off [z, index]=max(y); %计算最大值及其位置 x_best = x(index);%计算最大值对应的x值 toc