c65.m

MATLAB中常用的程序

function ObjVal = objlinq(Chrom,sswitch);

% 目标函数的维数
   Dim = 45;

%参数值
   x0 = 100;              % 初始点
   var = 1;               
   Para = [   1    1    1    1    1    16180.3399];             
   s = Para(var,1); r = Para(var,2); q = Para(var,3);
   a = Para(var,4); b = Para(var,5); GlobalMinimum = Para(var,6);

% 种群大小
   [Nind,Nvar] = size(Chrom);

% 检验染色体大小
   % 如果染色体为空，则定义边界矩阵大小和值
   if Nind == 0
      %返回输出图形的标题文字
      if sswitch == 2
         ObjVal = ['Linear-quadratic problem (dis)-' int2str(Dim)];
      %返回全局最小值
      elseif sswitch == 3
         ObjVal = GlobalMinimum;
      % 定义边界矩阵大小和值
      else   
         %上界和下界        
         ObjVal1 = [-100 -70 -50; 20 20 20];
         ObjVal = [ObjVal1 rep([-30;20],[1 Dim-3])];
      end
   % 计算适值函数值
   elseif Nvar == Dim
      ObjVal = zeros(Nind,1);
      X = zeros(Nind,Nvar+1);
      X(:,1) = rep(x0,[Nind 1]);
      for irun = 1:Nvar,
         X(:,irun+1) = a*X(:,irun) + b*Chrom(:,irun);
      end
      ObjVal = q * X(:,Nvar+1).^2 + sum((s * X(:,1:Nvar).^2 + r * Chrom.^2)')';
   % 否则，染色体格式错误
   else
      error('size of matrix Chrom is not correct for function evaluation');
   end