智能优化作业.docVIP

【问题】在－5=Xi=5,i=1,2区间内，求解??????? f(x1,x2)=-20*exp(-0.2*sqrt(0.5*(x1.^2x2.^2)))-exp(0.5*(cos(2*pi*x1)+ cos(2*pi*x2))) 22.71282的最小值。?【分析】种群大小10，最大代数1000，变异率0.1交叉率0.3?【运行结果】?p =? 0.0000 -0.0000 0.0055?【程序清单】??? ％源函数的matlab代码?????? function [eval]=f(sol)????????numv=size(sol,2);????????x=sol(1:numv);????????eval=-20*exp(-0.2*sqrt(sum(x.^2)/numv)))-exp(sum(cos(2*pi*x))/numv) 22.71282;??%适应度函数的matlab代码?????? function [sol,eval]=fitness(sol,options)????????numv=size(sol,2)-1;?????? x=sol(1:numv);???????eval=f(x);???????eval=-eval;??%遗传算法的matlab代码?????? bounds=ones(2,1)*[-5 5];?????? [p,endPop,bestSols,trace]=ga(bounds,fitness)? %初始化function pop=initpop(popsize,chromlength)?pop=round(rand(popsize,chromlength));%计算目标函数值function pop2=decodebinary(pop)[px,py]=size(pop); %求pop行和列数for i=1:pypop1(:,i)=2.^(py-i).*pop(:,i);endpop2=sum(pop1,2); %求pop1的每行之和%将二进制编码转化为十进制数function pop2=decodechrom(pop,spoint,length)pop1=pop(:,spoint:spoint+length-1);pop2=decodebinary(pop1);% 计算目标函数值function [objvalue]=calobjvalue(pop)temp1=decodechrom(pop,1,10); %将pop每行转化成十进制数x=temp1*10/1023; %将二值域 中的数转化为变量域 的数objvalue=10*sin(5*x)+7*cos(4*x); %计算目标函数值%计算个体的适应值function fitvalue=calfitvalue(objvalue)global Cmin;Cmin=0;[px,py]=size(objvalue);for i=1:pxif objvalue(i)+Cmin0temp=Cmin+objvalue(i);elsetemp=0.0;endfitvalue(i)=temp;endfitvalue=fitvalue;%选择复制function [newpop]=selection(pop,fitvalue)totalfit=sum(fitvalue); %求适应值之和fitvalue=fitvalue/totalfit; %单个个体被选择的概率fitvalue=cumsum(fitvalue); [px,py]=size(pop);ms=sort(rand(px,1)); %从小到大排列fitin=1;newin=1;while newin=pxif(ms(newin))fitvalue(fitin)newpop(newin)=pop(fitin);newin=newin+1;elsefitin=fitin+1;endend% 交叉function [newpop]=crossover(pop,pc)[px,py]=size(pop);newpop=ones(size(pop));for i=1:2:px-1if(randpc)cpoint=round(rand*py);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)];elsenewpop(i,:)=pop(i);newpop(i+1,:)=pop(i+1