《2022年遗传算法的MATLAB程序实例 .pdf》由会员分享，可在线阅读，更多相关《2022年遗传算法的MATLAB程序实例 .pdf（5页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、遗传算法的程序实例如求下列函数的最大值f(x)=10*sin(5x)+7*cos(4x) x0,10 一、初始化 (编码 ) initpop.m 函数的功能是实现群体的初始化，popsize 表示群体的大小，chromlength表示染色体的长度(二值数的长度 )，长度大小取决于变量的二进制编码的长度(在本例中取10 位)。代码 :%Name: initpop.m %初始化function pop=initpop(popsize,chromlength) pop=round(rand(popsize,chromlength); % rand随机产生每个单元为 0,1 行数为 popsize，列
2、数为 chromlength的矩阵，% roud对矩阵的每个单元进行圆整。这样产生的初始种群。二、计算目标函数值1、将二进制数转化为十进制数(1) 代码 :%Name: decodebinary.m %产生 2n 2(n-1) . 1 的行向量，然后求和，将二进制转化为十进制function pop2=decodebinary(pop) px,py=size(pop); %求 pop 行和例数for i=1:py pop1(:,i)=2.(py-1).*pop(:,i); py=py-1; end pop2=sum(pop1,2); %求 pop1 的每行之和2、将二进制编码转化为十进制数(2
3、) decodechrom.m函数的功能是将染色体(或二进制编码 )转换为十进制，参数spoint 表示待解码的二进制串的起始位置。 (对于多个变量而言，如有两个变量，采用20 为表示，每个变量10 为，则第一个变量从1开始，另一个变量从11 开始。本例为1)，参数 1ength 表示所截取的长度（本例为10） 。代码 :%Name: decodechrom.m %将二进制编码转换成十进制function pop2=decodechrom(pop,spoint,length) pop1=pop(:,spoint:spoint+length-1); pop2=decodebinary(pop1)
4、;3、计算目标函数值calobjvalue.m函数的功能是实现目标函数的计算，其公式采用本文示例仿真，可根据不同优化问题予以修改。代码 : %Name: calobjvalue.m %实现目标函数的计算function objvalue=calobjvalue(pop) temp1=decodechrom(pop,1,10); %将 pop 每行转化成十进制数名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 5 页 - - - - - - - - - x=temp1*10/
5、1023; %将二值域中的数转化为变量域的数objvalue=10*sin(5*x)+7*cos(4*x);%计算目标函数值三、计算个体的适应值代码 :%Name:calfitvalue.m %计算个体的适应值function fitvalue=calfitvalue(objvalue) global Cmin; Cmin=0; px,py=size(objvalue); for i=1:px if objvalue(i)+Cmin0 temp=Cmin+objvalue(i); else temp=0.0; end fitvalue(i)=temp; end fitvalue=fitvalu
6、e;四 选择复制选择或复制操作是决定哪些个体可以进入下一代。程序中采用赌轮盘选择法选择，这种方法较易实现。根据方程pi=fi/fi=fi/fsum ，选择步骤：1）在第t 代，由（ 1）式计算fsum 和 pi 2）产生0,1 的随机数rand( .) ，求s=rand( .)*fsum 3）求 fi s 中最小的k ，则第k 个个体被选中4）进行N 次 2） 、3）操作，得到N 个个体，成为第t=t+1 代种群代码 : %Name: selection.m %选择复制function newpop=selection(pop,fitvalue) totalfit=sum(fitvalue);
7、%求适应值之和fitvalue=fitvalue/totalfit;%单个个体被选择的概率fitvalue=cumsum(fitvalue); % 累积概率， 如 fitvalue=1 2 3 4 ，则 cumsum(fitvalue)=1 3 6 10 px,py=size(pop); ms=sort(rand(px,1); %从小到大排列，将rand(px,1)产生的一列随机数变成轮盘赌形式的表示方法，由小到大排列fitin=1; %fivalue是一向量， fitin代表向量中元素位，即fitvalue(fitin)代表第 fitin个个体的单个个体被选择的概率newin=1; %同理w
8、hile newin=px if(ms(newin)fitvalue(fitin) %ms(newin)表示的是ms 列向量中第 newin位数值，同理 fitvalue(fitin) newpop(newin,:)=pop(fitin,:); %赋值 , 即将旧种群中的第 fitin个个体保留到下一代名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 5 页 - - - - - - - - - (newpop) newin=newin+1; else fitin=fitin
9、+1; end end五 2.5 交叉交叉 (crossover) ，群体中的每个个体之间都以一定的概率pc 交叉，即两个个体从各自字符串的某一位置（一般是随机确定）开始互相交换，这类似生物进化过程中的基因分裂与重组。例如，假设2 个父代个体x1，x2 为：x1=0100110 x2=1010001 从每个个体的第3 位开始交叉，交又后得到2 个新的子代个体y1，y2 分别为：y10100001 y21010110 这样 2 个子代个体就分别具有了2 个父代个体的某些特征。利用交又我们有可能由父代个体在子代组合成具有更高适合度的个体。事实上交又是遗传算法区别于其它传统优化方法的主要特点之一。代
10、码 : %Name: crossover.m %交叉function newpop=crossover_multiv(pop,pc) global Numv px,py=size(pop); m=py/Numv; for j=1:Numv pop1=ones(px,m); pop2=pop(:,m*(j-1)+1:m*j); %取出相应变量对应的二进制编码段 for i=1:2:px-1 if(randpc) cpoint=round(rand*(m-1); %cpoint为交叉点 pop1(i,:)=pop2(i,1:cpoint) pop2(i+1,cpoint+1:m); pop1(i+
11、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六、变异变异 (mutation) ，基因的突变普遍存在于生物的进化过程中。变异是指父代中的每个个体的每一位都以概率pm 翻转， 即由 “1”变为 “0”，或由 “0”变为 “1”。遗传算法的变异特性可以使求解过程随机地搜索到解可能存在名师资料总结 - - -精品资料欢迎下载 -
12、- - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 5 页 - - - - - - - - - 的整个空间，因此可以在一定程度上求得全局最优解。代码 : %Name: mutation.m %变异function newpop=mutation(pop,pm) px,py=size(pop); newpop=ones(size(pop); for i=1:px if(randpm) mpoint=round(rand*py); if mpointbestfit bestindividual=pop(i,:); bestf
13、it=fitvalue(i); end end八、主程序代码 : %Name:genmain05.m clear clf popsize=20; %群体大小chromlength=10; %字符串长度（个体长度）名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 5 页 - - - - - - - - - pc=0.6; %交叉概率pm=0.001; %变异概率pop=initpop(popsize,chromlength); %随机产生初始群体for i=1:20 %20为
14、迭代次数objvalue=calobjvalue(pop); %计算目标函数fitvalue=calfitvalue(objvalue); %计算群体中每个个体的适应度newpop=selection(pop,fitvalue); %复制newpop=crossover(pop,pc); %交叉newpop=mutation(pop,pc); %变异bestindividual,bestfit=best(pop,fitvalue);%求出群体中适应值最大的个体及其适应值y(i)=max(bestfit); n(i)=i; pop5=bestindividual; x(i)=decodechrom(pop5,1,chromlength)*10/1023; pop=newpop; end fplot(10*sin(5*x)+7*cos(4*x),0 10) hold on plot(x,y,r*) hold off名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 5 页 - - - - - - - - -