简单的遗传算法MATLAB实现精品毕业设计完整版文档格式.docx
- 文档编号:15293668
- 上传时间:2022-10-29
- 格式:DOCX
- 页数:16
- 大小:35KB
简单的遗传算法MATLAB实现精品毕业设计完整版文档格式.docx
《简单的遗传算法MATLAB实现精品毕业设计完整版文档格式.docx》由会员分享,可在线阅读,更多相关《简单的遗传算法MATLAB实现精品毕业设计完整版文档格式.docx(16页珍藏版)》请在冰豆网上搜索。
对于函数优化问题,一般来说,有两种编码方式,一是实数编码,一是二进制编码,两者各有优缺点,二进制编码具有稳定性高、种群多样性大等优点,但是需要的存储空间大,需要解码过程并且难以理解;
而实数编码直接用实数表示基因,容易理解并且不要解码过程,但是容易过早收敛,从而陷入局部最优。
本文以最常用的二进制编码为例,说明遗传编码的过程。
从遗传算法求解的过程来看,需要处理好两个空间的问题,一个是编码空间,另一个是解空间,如下图所示
从解空间到编码空间的映射过程成为编码过程;
从编码空间到解空间的映射过程成为解码过程。
下面就以求解一个简单的一维函数f(x)=-(x-1)^2+4,x的取值范围为[-1,3]最大值为例,来说明编码及解码过程。
编码:
在编码之前需要确定求解的精度,在这里,我们设定求解的精度为小数点后四位,即1e-4。
这样可以将每个自变量xi的解空间划分为个等分。
以上面这个函数为例,即可以将x的解空间划分为(3-(-1))*1e+4=40000个等分。
使ni满足,这里ni表示使上式成立的最小整数,即表示自变量xi的基因串的长度。
因为215<
40000<
216
,这里ni取16。
例如0000110110000101就表示一个解空间中的基因串。
表示所有自变量x=(x1,x2,…,xk)的二进制串的总长度称为一个染色体(Chromosome)的长度或者一个个体(Individual)的长度,。
编码过程一般在实现遗传算法之前需要指定。
解码:
解码即将编码空间中的基因串翻译成解空间中的自变量的实际值的过程。
对于二进制编码而言,每个二进制基因串都可以这样翻译成一个十进制实数值,。
例如基因串0000110110000101,可以翻译为,这里二进制基因串转变成十进制是从左至右进行的。
1.2初始化种群
在开始遗传算法迭代过程之前,需要对种群进行初始化。
设种群大小为pop_size,每个染色体或个体的长度为chromo_size,种群的大小决定了种群的多样性,而染色体的长度则是由前述的编码过程决定的。
一般随机生成初始种群,但是如果知道种群的实际分布,也可以按照此分布来生成初始种群。
假设生成的初始种群为(v1,v2,…,vpop_size)。
1.3选择操作
选择操作即从前代种群中选择个体到下一代种群的过程。
一般根据个体适应度的分布来选择个体。
以初始种群(v1,v2,…,vpop_size)为例,假设每个个体的适应度为(fitness(v1),fitness(v2),…,fitness(vpop_size)),一般适应度可以按照解码的过程进行计算。
以轮盘赌的方式选择个体,如下图
随机转动一下轮盘,当轮盘停止转动时,若指针指向某个个体,则该个体被选中。
很明显,具有较高适应度的个体比具有较低适应度的个体更有机会被选中。
但是这种选择具有随机性,在选择的过程中可能会丢失掉比较好的个体,所以可以使用精英机制,将前代最优个体直接选到下一代中。
轮盘赌选择具体算法如下(这里假定种群中个体是按照适应度从小到大进行排列的,如果不是,可以按照某种排序算法对种群个体进行重排):
SelectionAlgorithm
varpop,pop_new;
/*pop为前代种群,pop_new为下一代种群*/
varfitness_value,fitness_table;
/*fitness_value为种群的适应度,fitness_table为种群累积适应度*/
fori=1:
pop_size
r=rand*fitness_table(pop_size);
/*随机生成一个随机数,在0和总适应度之间,因为fitness_table(pop_size)为最后一个个体的累积适应度,即为总适应度*/
first=1;
last=pop_size;
mid=round((last+first)/2);
idx=-1;
/*下面按照排中法选择个体*/
while(first<
=last)&
&
(idx==-1)
ifr>
fitness_table(mid)
first=mid;
elseifr<
last=mid;
else
idx=mid;
break;
endif
if(last-first)==1
idx=last;
endwhile
forj=1:
chromo_size
pop_new(i,j)=pop(idx,j);
endfor
/*是否精英选择*/
ifelitism
p=pop_size-1;
p=pop_size;
p
pop(i,j)=pop_new(i,j);
/*若是精英选择,则只将pop_new前pop_size-1个个体赋给pop,最后一个为前代最优个体保留*/
1.3交叉操作
交叉操作是对任意两个个体进行的(在这里我们实现的算法是直接对相邻的两个个体进行的)。
随机选择两个个体,如下图所示
然后随机生成一个实数0<
=r<
=1,如果r<
cross_rate,0<
cross_rate<
1为交叉概率,则对这两个个体进行交叉,否则则不进行。
如果需要进行交叉,再随机选择交叉位置(rand*chromo_size),如果等于0或者1,将不进行交叉。
否则将交叉位置以后的二进制串进行对换(包括交叉位置)。
(注意:
有时候还可以进行多点交叉,但是这里只讨论单点交叉的情况)
单点交叉具体算法如下:
Crossoveralgorithm
2:
if(rand<
cross_rate)/*cross_rate为交叉概率*/
cross_pos=round(rand*chromo_size);
/*交叉位置*/
ifor(cross_pos==0,cross_pos==1)
continue;
/*若交叉位置为0或1,则不进行交叉*/
forj=cross_pos:
pop(i,j)<
->
pop(i+1,j);
/*交换*/
1.4变异操作
变异操作是对单个个体进行的。
首先生成一个随机实数0<
mutate_rate,则对此个体进行变异操作,0<
mutate_rate<
1为变异概率,一般为一个比较小的实数。
对每一个个体,进行变异操作,如下图所示
如个体需要进行变异操作,首先需要确定变异位置(rand*chromo_size),若为0则不进行变异,否则则对该位置的二进制数字进行变异:
1变成0,0变成1.(当然也可以选择多点进行变异)
单点变异的具体算法描述如下:
Mutationalgorithm
ifrand<
mutate_rate/*mutate_rate为变异概率*/
mutate_pos=round(rand*chromo_size);
/*变异位置*/
ifmutate_pos==0
/*若变异位置为0,则不进行变异*/
pop(i,mutate_pos)=1-pop(i,mutate_pos);
/*将变异位置上的数字至反*/
1.5遗传算法流程
遗传算法计算流程图如下图所示
1.6MATLAB程序实现
初始化:
%初始化种群
%pop_size:
种群大小
%chromo_size:
染色体长度
functioninitilize(pop_size,chromo_size)
globalpop;
pop(i,j)=round(rand);
end
cleari;
clearj;
计算适应度:
(该函数应该根据具体问题进行修改,这里优化的函数是前述的一维函数)
%计算种群个体适应度,对不同的优化目标,此处需要改写
functionfitness(pop_size,chromo_size)
globalfitness_value;
globalG;
fitness_value(i)=0.;
ifpop(i,j)==1
fitness_value(i)=fitness_value(i)+2^(j-1);
fitness_value(i)=-1+fitness_value(i)*(3.-(-1.))/(2^chromo_size-1);
fitness_value(i)=-(fitness_value(i)-1).^2+4;
对个体按照适应度大小进行排序:
%对个体按适应度大小进行排序,并且保存最佳个体
functionrank(pop_size,chromo_size)
globalfitness_table;
globalfitness_avg;
globalbest_fitness;
globalbest_individual;
globalbest_generation;
pop_size
fitness_table(i)=0.;
min=
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 简单 遗传 算法 MATLAB 实现 精品 毕业设计 完整版
![提示](https://static.bdocx.com/images/bang_tan.gif)