Documentation of objharv

Global Index (all files) (short | long) | Local contents | Local Index (files in subdir) (short | long)

Function Synopsis

ObjVal = objharv(Chrom, option);

Help text

 OBJective function for HARVest problem

 This function implements the HARVEST PROBLEM.

 Syntax:  ObjVal = objharv(Chrom, option)

 Input parameters:
    Chrom     - Matrix containing the chromosomes of the current
                population. Each row corresponds to one individual's
                string representation.
                if Chrom == [], then speziell values will be returned
    option    - if Chrom == [] and
                option == 1 (or []) return boundaries
                option == 2 return title
                option == 3 return value of global minimum

 Output parameters:
    ObjVal    - Column vector containing the objective values of the
                individuals in the current population.
                if called with Chrom == [], then ObjVal contains
                option == 1, matrix with the boundaries of the function
                option == 2, text for the title of the graphic output
                option == 3, value of global minimum
                
 See also: objdopi, objlinq, objpush

Listing of function objharv

% Author:     Hartmut Pohlheim
% History:    18.02.94     file created (copy of vallinq.m)
%             01.03.94     name changed in obj*
%             17.02.95     direct Dim removed and function cleaned

function ObjVal = objharv(Chrom, option);

   
% Compute population parameters
   [Nind, Nvar] = size(Chrom);

% Check size of Chrom and do the appropriate thing
   % if Chrom is [], then define size of boundary-matrix and values
   if Nind == 0
      % Default dimension of objective function
      Dim = 20;
      % return text of title for graphic output
      if option == 2
         ObjVal = ['HARVEST PROBLEM'];
      % return value of global minimum
      elseif option == 3
         ObjVal = -sqrt(x0*(a^Dim-1)^2/(a^(Dim-1)*(a-1)));
      % define size of boundary-matrix and values
      else   
         % lower and upper bound, identical for all n variables        
         ObjVal = repmat([0; 10*Dim], [1 Dim]);
      end
   % compute values of function
   else
      % values from MICHALEWICZ
      a = 1.1; x0 = 100; xend = x0; XENDWEIGHT = 0.4/(Nvar^0.6);
      % Start computation of objective function
      ObjVal = zeros(Nind,1);
      X = repmat(x0, [Nind 1]);
      for irun = 1:Nvar,
         X = a * X - Chrom(:,irun);
      end
      ObjVal = -(sum(sqrt(Chrom)')' - XENDWEIGHT * abs(X - x0));
   end   


% End of function