本文程式碼基於謝菲爾德大學遺傳演算法工具箱。

問題：求在0~4範圍內的最大值

figure(1);
% fplot('variable.*sin(10*pi*variable)+2.0',[-1,2]); %畫出函式曲線
fplot('-variable.*(variable - 4)',[0,5]);   %畫出函式曲線
%定義遺傳演算法引數
NIND=40;        %個體數目(Number of individuals)
MAXGEN=25;      %最大遺傳代數(Maximum number of generations)
PRECI=20;       %變量的二進位制位數(Precision of variables)
GGAP=0.9
 
;       %代溝(Generation gap)
trace=zeros(2, MAXGEN);                        %尋優結果的初始值
FieldD=[20;0;5;1;0;1;1];                      %區域描述器(Build field descriptor)表示染色體個體到實際的數字的對映關係
Chrom=crtbp(NIND, PRECI);                      %初始種群
gen=0;                                         %代計數器
variable=bs2rv(Chrom, FieldD);             %計
 
算初始種群的十進位制轉換
ObjV=-variable.*(variable - 4);        %計算目標函式值
while gen<MAXGEN
   FitnV=ranking(-ObjV);     %分配適應度值(Assign fitness values)         
   SelCh=select('sus', Chrom, FitnV, GGAP);               %選擇
   SelCh=recombin('xovsp', SelCh, 0.7);                   %重組
   SelCh=mut(SelCh);                                      %變
 
異
   variable=bs2rv(SelCh, FieldD); %子代個體的十進位制轉換
   ObjVSel=-variable.*(variable - 4);%計運算元代的目標函式值
   [Chrom ObjV]=reins(Chrom, SelCh, 1, 1, ObjV, ObjVSel);
                                          %重插入子代的新種群
   variable=bs2rv(Chrom, FieldD);
   gen=gen+1;                            %代計數器增加
   %輸出最優解及其序號，並在目標函式影象中標出，Y為最優解,I為種群的序號
   [Y, I]=max(ObjV);hold on;
   plot(variable(I), Y, 'bo');
   trace(1, gen)=max(ObjV);             %遺傳演算法效能跟蹤
   trace(2, gen)=sum(ObjV)/length(ObjV);
end
variable=bs2rv(Chrom, FieldD);            %最優個體的十進位制轉換
hold on, grid;
plot(variable,ObjV,'b*');
figure(2);
plot(trace(1,:));
hold on;
plot(trace(2,:),'-.');grid
legend('解的變化','種群均值的變化')

步驟：

1. 初始化種群
2. 分配適應度值FitnV
3. 選擇、重組、變異
4. 子代插入父代形成新種群
5. 回到2，迴圈直至最大代數
6. 得到最優種群，即最優解

本文程式碼大部分來自《MATLAB_遺傳演算法工具箱與應用》雷英傑 一書