hilloptions.m

爬山算法 可以解决许多函数优化方面的问题,可以与其他算法,如蚂蚁算法,粒子群算法等融合使用,有其一定的研究意义

function options=hilloptions(varargin)
%Syntax: options=hilloptions(varargin)
%_____________________________________
%
% Options definition for hill.
%
% options is the options struct:
%          'space': the 2-column matrix with the boundaries of the
%                   initial solution
%         'climbs': the number of climbs
%           'prec': the number of the decimal digits arithmetic precision
%           'line': the number of additional steps if a climb is made
%        'Display': displays the progress (or part of it) on the screen
%      'TimeLimit': time in seconds before the optimization stops
%           'Goal': a value to be reached
%
%
% Example:
%    options = hilloptions('space',[-ones(8,1) ones(8,1)],'climbs',8);
%    options = hilloptions(options,'line',1000,'Goal',0);
%
% Alexandros Leontitsis
% Department of Education
% University of Ioannina
% 45110- Dourouti
% Ioannina
% Greece
% 
% University e-mail: [e-mail]
% Lifetime e-mail: [e-mail]
% Homepage: http://www.geocities.com/CapeCanaveral/Lab/1421
% 
% 8 Dec, 2004.


options.space = [0 1];
options.MaxIter = 200;
options.prec = 4;
options.line = 99;
options.Display = 'Final';
options.TimeLimit = 30;
options.Goal = -inf;

if (nargin == 0) && (nargout == 0)
    fprintf('          space: [ real matrix  | {[0 1]} ]\n');
    fprintf('        MaxIter: [ integer positive scalar | {200} ]\n');
    fprintf('           prec: [ integer positive scalar | {4} ]\n');
    fprintf('           line: [ integer positive scalar | {99} ]\n');
    fprintf('        Display: [ integer non-negative scalar | {"Final"} ]\n');
    fprintf('      TimeLimit: [ real positive scalar | {30} ]\n');
    fprintf('           Goal: [ real scalar | {-inf} ]\n');
end

if nargin > 1
    if rem(nargin,2)==0
        for i=1:nargin/2
            switch varargin{i*2-1}
                case 'space'
                    % varargin{i*2} must be a 2 dimensinal matrix
                    if ndims(varargin{i*2})>3
                        error(['"' varargin{i*2-1} '" must be a 2 dimensinal matrix.']);
                    end
                    % varargin{i*2-1} must contain exactly 2 columns
                    if size(varargin{i*2},2)~=2
                        error(['"' varargin{i*2-1} '" must contain exactly 2 columns.']);
                    end
                case {'MaxIter','prec','TimeLimit'}
                    %  varargin{i*2} must be a scalar
                    if sum(size(varargin{i*2}))>2
                        error([varargin{i*2-1} ' must be a scalar.']);
                    end
                    % varargin{i*2} must be positive
                    if varargin{i*2}<=0
                        error([varargin{i*2-1} ' must be positive.']);
                    end
                case {'line'}
                    %  varargin{i*2} must be a scalar
                    if sum(size(varargin{i*2}))>2
                        error([varargin{i*2-1} ' must be a scalar.']);
                    end
                    % varargin{i*2} must be non-negative
                    if varargin{i*2}<0
                        error([varargin{i*2-1} ' must be non-negative.']);
                    end
                case {'Display'}
                    if isa(varargin{i*2},'numeric')==1
                        %  varargin{i*2} must be a scalar
                        if sum(size(varargin{i*2}))>2
                            error([varargin{i*2-1} ' must be a scalar.']);
                        end
                        % varargin{i*2} must be positive
                        if varargin{i*2}<=0
                            error([varargin{i*2-1} ' must be positive.']);
                        end
                        % varargin{i*2} must be an integer
                        if varargin{i*2}~=round(varargin{i*2})
                            error([varargin{i*2-1} ' must be an integer.']);
                        end
                    elseif strcmp(varargin{i*2},'Final')==0
                        error('The only non-numeric value of "Show" is "Final".');
                    end
                case 'Goal'
                    %  varargin{i*2} must be a scalar
                    if sum(size(varargin{i*2}))>2
                        error([varargin{i*2-1} ' must be a scalar.']);
                    end
                otherwise
                    error(['The field "' varargin{i*2-1} '" does not exist in the struct hilloptions.']);
            end
            options = setfield(options,varargin{i*2-1},varargin{i*2});
        end
    else
        error('The number of input arguments must be even.');
    end
end