example5_44.m

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

function ObjVal = objpush(Chrom,sswitch);
% 目标函数的维数
   Dim = 20;
% 初始值
   x0 = [0 0];

% 计算染色体参数
   [Nind,Nvar] = size(Chrom);

% 染色体检查
   % 如果染色体为空，则定义边界矩阵的大小和值
   if Nind == 0
      %返回输出图形的标题文字
      if sswitch == 2
         ObjVal = ['PUSH-CART PROBLEM-' int2str(Dim)];
      %返回全局最小值
      elseif sswitch == 3
         ObjVal = -(1/3 - ((3*Dim-1)/(6*Dim^2)) - (1/(2*Dim^3))*sum((1:Dim-1).^2));
      %定义边界矩阵大小和值
      else   
         %变量的上下边界        
         ObjVal = [0; 5];
         ObjVal = rep(ObjVal,[1 Dim]);
      end
   
   elseif Nvar == Dim
      ObjVal = zeros(Nind,1);
      X = rep(x0,[Nind 1]);
      for irun = 1:Nvar,
         Xsave = X;
         X(:,1) = Xsave(:,2);
         X(:,2) = 2 * X(:,2) - Xsave(:,1) + (1/Dim^2) * Chrom(:,irun);
      end
      X;
      ObjVal = -(X(:,1) - (1/(2*Dim)) * sum((Chrom.^2)')');
   %否则，染色体格式错误
   else
      error('size of matrix Chrom is not correct for function evaluation');
   end