遗传算法MATLAB实现(3):多元函数优化举例

多峰的Shubert为：

　　

求f(x,y)在[-10,10]x[-10,10]上的最大值。

 

MATLAB代码：

    fun_mutv函数为：

function my=fun_mutv(x,y)
t1=zeros(size(x));
t2=t1;
for i=1:5
    t1=t1+i*cos((i+1)*x+i);
    t2=t2+i*cos((i+1)*y+i);
end
my=t1.*t2;

　　

opt_minmax=1;   %优化目标类型：1最大化 0最小化
num_ppu=60;     %种群规模，个体个数。
num_gen=100;    %最大遗传代数
num_v=2;        %变量个数
len_ch=20;      %基因长度
gap=0.9;        %代沟
sub=-10;        %变量取值下限
up=10;          %变量取值上限
cd_gray=1;      %是否选择格雷码编码方式 1是，0否
sc_log=0;           %是否选择对数标度：1是，0否
trace=zeros(num_gen,2);     %遗传迭代性能跟踪器
%区域描述器，rep为矩阵复制函数
fieldd=[rep([len_ch],[1,num_v]);rep([sub;up],[1,num_v]);rep([1-cd_gray;sc_log;1;1],[1,num_v])];
chrom=crtbp(num_ppu,len_ch*num_v);  %初始化生产种群
k_gen=0;
x=bs2rv(chrom,fieldd);      %翻译初始化种群为10进制
fun_v=fun_mutv(x(:,1),x(:,2));  %计算目标函数值
[tx,ty]=meshgrid(-10:1:10); 
mesh(tx,ty,fun_mutv(tx,ty));
xlabel('x');ylabel('y');zlabel('z');
title('多元函数优化结果');
hold on;
while k_gen 
   

　　

转载于:https://www.cnblogs.com/zxhyxiao/p/9401230.html