2025年遗传算法入门新手matlab遗传算法程序.docVIP

%求以下函数的最大值%

%f(x）＝１０*sin(5ｘ)+７*ｃos(４ｘ)x∈[0，10]%

%将x的值用一个10位的二值形式表达为二值问题，一个10位的二值数提供的辨别率是每为（10-0)/(2^１0-1)≈0.0１。％

％将变量域［0,1０]离散化为二值域[0，1023],ｘ＝０+1０＊b/102３,其中b是［０,1023]中的一个二值数。％

％%

%-----－--－-－------－-----------－-－--－-－---－---------－------------------－--------－------－-----－－----------－-----－％

%－------------－-------------－-----------－---－---－－－－---------------－-－-----－－-----------－-----－---－－-－--－----－-％

%编程

％－--－--－---－－---------－-－--------－-----－--－－－－－-

%2．1初始化(编码）

%inｉtpｏp.m函数的功效是实现群体的初始化,popsize表达群体的大小，chｒomlｅnｇtｈ表达染色体的长度(二值数的长度)，

%长度大小取决于变量的二进制编码的长度（在本例中取１0位）。

％遗传算法子程序

%Ｎame:iｎｉｔpop.m

%初始化

funｃtｉonpop=ｉnitｐoｐ(popsize,cｈroｍｌeｎgth)

poｐ=rｏuｎd(rand(pｏpsｉze,chroｍｌｅngｔh));％rand随机产生每个单元为{０,1}行数为popsizｅ,列数为ｃhromｌength的矩阵，

%rｏund对矩阵的每个单元进行圆整。这么产生的初始种群。

%2.2计算目标函数值

%2.2.1将二进制数转化为十进制数（1）

%遗传算法子程序

%Name:dｅcodebinary.m

％产生[2^ｎ２^(ｎ-1）.．.1]的行向量，然后求和,将二进制转化为十进制

functｉonpoｐ2=deｃodｅbinaｒy（pop)

［px,pｙ]=sizｅ(pop);%求poｐ行和列数

ｆｏrｉ=1:pｙ

pop1（：，i)=2.^(py－ｉ）.*pop(:,i）；

enｄ

ｐop2＝suｍ(pｏp1,2);%求ｐop1的每行之和

%２.2.2将二进制编码转化为十进制数（２)

％decodechrom.ｍ函数的功效是将染色体（或二进制编码)转换为十进制,参数ｓpｏiｎt表达待解码的二进制串的起始位置

%(对于多个变量而言,如有两个变量，采取20为表达，每个变量1０为，则第一个变量从1开始，另一个变量从11开始。本例为1),

%参数1engｔh表达所截取的长度（本例为1０）。

%遗传算法子程序

%Name：decｏｄeｃhrom.m

%将二进制编码转换成十进制

ｆｕnctionｐop２=decｏdechrom(pop,spoｉnt,ｌenｇth)

pｏｐ１=pop(：,sｐoinｔ:spoｉｎt+lｅｎｇth-1);

poｐ2＝decoｄeｂinａｒｙ(ｐop1);

％２.２.３计算目标函数值

%ｃａlobjvalue.m函数的功效是实现目标函数的计算，其公式采取本文示例仿真，可依照不一样优化问题予以修改。

%遗传算法子程序

%Name:calobjvａlue.m

%实现目标函数的计算

fｕnctiｏn[oｂｊvalue］=calｏｂｊvalue(pop)

ｔemp1=ｄecodechｒｏm（pop,1,10)；%将pop每行转化成十进制数

ｘ=tｅmｐ1*10/10２3；%将二值域中的数转化为变量域的数

objｖａｌue＝１0*sin(5＊x)+7＊cｏs(4*x)；%计算目标函数值

%2．3计算个体的适应值

％遗传算法子程序

％Naｍe:ｃalfitvalue.m

%计算个体的适应值

ｆunctionｆitｖａｌue=cａlfｉtvalｕe(ｏbjvaｌue)

gｌobalＣmin;

Cmin＝0;

[pｘ,py］＝size(oｂjvaｌuｅ)；

fｏri＝1:px

ifｏbjｖａlue(i)+Cmin０

ｔemp=Cmiｎ＋ｏbｊvalｕe(i);

ｅｌse

temp=0.0;

ｅｎd

fiｔvalue(i)=teｍp;

eｎd

ｆitｖａlue=fｉtvaｌue;

％2.4选择复制

％选择或复制操作是决定哪些个体可以进入下一代。程序中采取赌轮盘选择法选择,这种方法较易实现。

%依照方程ｐi＝fi/∑ｆi=ｆi／ｆsum，选择环节:

%1）在