接上一个遗传算法原理的例子,记录具体是如何编程的。编码方式是二进制编码,选择函数采用轮盘赌方法。
函数与函数图像
这一次我们的函数是
y=10*sin(5*x)+7*abs(x-5)+10
图像如图所示:
染色体长度的确定
因为我只想精确到小数点后两位,区间[0,10],所以应该需要精度为10*10^2=1000,2^10 = 1024,所以染色体长度为10已经足够。
主函数
function main()clear;clc;%种群大小(规模)popsize=100;%二进制编码长度chromlength=10;%交叉概率pc = 0.6;%变异概率pm = 0.001;%初始种群pop = initpop(popsize,chromlength);for i = 1:100 %计算适应度值(函数值) objvalue = cal_objvalue(pop); fitvalue = objvalue; %选择操作,选出优良的父母 newpop = selection(pop,fitvalue); %交叉操作,优良父母交叉产生新后代 newpop = crossover(newpop,pc); %变异操作,新种群跳一跳改变一下 newpop = mutation(newpop,pm); %更新种群 pop = newpop; %寻找最优解 [bestindividual,bestfit] = best(pop,fitvalue); x2 = binary2decimal(bestindividual);%解码最优个体 x1 = binary2decimal(newpop);%解码新种群 y1 = cal_objvalue(newpop);%计算新种群对应的 if mod(i,10) == 0%每迭代十次画一下图 figure; fplot('10*sin(5*x)+7*abs(x-5)+10',[0 10]); hold on; plot(x1,y1,'*'); title(['迭代次数为n=' num2str(i)]); %plot(x1,y1,'*'); endendfprintf('The best X is --->>%5.2f',x2);fprintf('The best Y is --->>%5.2f',bestfit);
初始化种群函数(对种群进行编码)
function pop=initpop(popsize,chromlength)%初始化种群大小%输入变量:%popsize:种群大小%chromlength:染色体长度-->>转化的二进制长度%输出变量:%pop:种群pop = round(rand(popsize,chromlength));%rand(3,4)生成3行4列的0-1之间的随机数% rand(3,4)% % ans =% % 0.8147 0.9134 0.2785 0.9649% 0.9058 0.6324 0.5469 0.1576% 0.1270 0.0975 0.9575 0.9706%round就是四舍五入% round(rand(3,4))=% 1 1 0 1% 1 1 1 0% 0 0 1 1%所以返回的种群就是每行是一个个体,列数是染色体长度
二进制转十进制函数(解码)
function pop2 = binary2decimal(pop)%功能:二进制转化成十进制函数%输入变量:%二进制种群%输出变量%十进制数值[px,py]=size(pop);for i = 1:py%py是染色体的长度 pop1(:,i) = 2.^(py-i).*pop(:,i);end%sum(.,2)对行求和,得到列向量temp = sum(pop1,2);pop2 = temp*10/1023;
适应度函数
function [objvalue] = cal_objvalue(pop)%计算函数目标值%输入变量:二进制数值%输出变量:目标函数值x = binary2decimal(pop);%转化二进制数为x变量的变化域范围的数值objvalue=10*sin(5*x)+7*abs(x-5)+10;
选择函数(为轮盘赌方法)
在开始学习的时候就经常会有一些疑问就是,遗传算法既然是选择优秀的个体(适应度较高的个体),那么应当淘汰一些劣质的个体,执行选择操作之后,新种群应该比原来的种群要小,但是结果运行却是不变的,为什么呢?后来我想明白了就是,如果我们在每次迭代的过程中选择淘汰一些数量,就要保证我们迭代的次数达到后,种群还有一些个体存在,所以我们要适当计算一下每次淘汰的个数是多少。而经典的遗传算法当中,种群大小在执行选择函数之后仍是不变的,因为轮盘赌法在每次迭代过程中会保证有一个个体被选出才会进入下一次循环,而这个数的选出是因为这个数的概率较大(大于随机生成的一个概率),所以新种群中选出的个体有很多是相同的个体,最后达到收敛(种群中所有个体都是一样的,集中在最优秀的地方),下面看算法的操作。
function [newpop] = selection(pop,fitvalue)%如何选择新的个体%输入变量:pop二进制种群,fitvalue:适应度值%输出变量:newpop选择以后的二进制种群%构造轮盘[px,py] = size(pop);totalfit = sum(fitvalue);p_fitvalue = fitvalue/totalfit;index = [];for i = 1:px %px是种群数 pick = rand;while pick == 0pick = rand;end for j = 1:px pick = pick - p_fitvalue(j); if pick < 0 index = [index j]; break; end endfunction [newpop] = selection(pop,fitvalue)%如何选择新的个体%输入变量:pop二进制种群,fitvalue:适应度值%输出变量:newpop选择以后的二进制种群%构造轮盘[px,py] = size(pop);totalfit = sum(fitvalue);p_fitvalue = fitvalue/totalfit;p_fitvalue = cumsum(p_fitvalue);%概率求和排序(这部我没理解为何要求和排序,不求和排序不行吗,评论区求解答[笑哭])ms = sort(rand(px,1));%生成100个(种群数)0-1的数(随机概率)并从小到大排列fitin = 1;newin = 1;while newin<=px if((ms(newin))
交叉函数(染色体相互交叉)
function [newpop] = crossover(pop,pc)%交叉变换%输入变量:pop:二进制的父代种群数,pc:交叉的概率%输出变量:newpop:交叉后的种群数[px,py] = size(pop);newpop = ones(size(pop));for i = 1:2:px-1%相邻的2个染色体(个体)进行交叉,所以是1:2:px-1 if(rand>pc)%如果大于设定的交叉的概率pc,就交叉出新个体,就是原本的父母就丢弃,只要新生儿 cpoint = round(rand*py);%产生随机数,随机选择一个位置点cpoint newpop(i,:) = [pop(i,1:cpoint),pop(i+1,cpoint+1:py)]; newpop(i+1,:) = [pop(i+1,1:cpoint),pop(i,cpoint+1:py)]; else newpop(i,:) = pop(i,:); newpop(i+1,:) = pop(i+1,:); endend
变异函数
function [newpop] = mutation(pop,pm)%关于变异%输入变量:pop:二进制种群,pm:变异概率%输出变量:newpop变异以后的种群[px,py] = size(pop);newpop = ones(size(pop));for i = 1:px if(rand
选择最优个体函数
function [bestindividual bestfit] = best(pop,fitvalue)%求最优适应度函数%输入变量:pop:种群,fitvalue:种群适应度%输出变量:bestindividual:最佳个体,bestfit:最佳适应度值[px,py] = size(pop);bestindividual = pop(1,:);bestfit = fitvalue(1);%其实可以直接max就行,不用冒泡排序之类的for i = 2:px if fitvalue(i)>bestfit bestindividual = pop(i,:); bestfit = fitvalue(i); endend
效果如图
The best X is --->> 0.22
The best Y is --->>52.29