example5_42.m

MATLAB7辅助控制系统设计与仿真 第5章

function ObjVal = objharv(Chrom,sswitch);
% 目标函数的维数
   Dim = 20;
% 参数设置
   a = 1.1;
   x0 = 100;
   xend = x0;
   XENDWEIGHT = 0.4/(Dim^0.6);
   
% 种群的大小
   [Nind,Nvar] = size(Chrom);

% 种群检查
   if Nind == 0
      %返回输出图形的标题文字
      if sswitch == 2
         ObjVal = ['HARVEST PROBLEM-' int2str(Dim)];
      %返回全局最小值
      elseif sswitch == 3
         ObjVal = -sqrt(x0*(a^Dim-1)^2/(a^(Dim-1)*(a-1)));
      %定义边界矩阵大小和值
      else   
         % 变量的边界        
         ObjVal1 = [0; 10*Dim];
         ObjVal = rep(ObjVal1,[1 Dim]);
      end
   elseif Nvar == Dim
      ObjVal = zeros(Nind,1);
      X = rep(x0,[Nind 1]);
      for irun = 1:Nvar,
         X = a*X - Chrom(:,irun);
      End
X;
ObjVal = -(sum(sqrt(Chrom)')' - XENDWEIGHT * abs(X-x0));
%否则，染色体格式错误
else
error('size of matrix Chrom is not correct for function evaluation');
   end