WOA源码

% The Whale Optimization Algorithm
function [Leader_posConvergence_curve]=WOA（SearchAgents_noMaxFEslbubdimfobj）

% initialize position vector and score for the leader
Leader_pos=zeros（1dim）;
Leader_score=inf; %change this to -inf for maximization problems


%Initialize the positions of search agents
Positions=initialization（SearchAgents_nodimublb）;

Convergence_curve=[];
FEs=0;
t=1;
% Main loop
while  FEs < MaxFEs
    for i=1:size（Positions1）
        
        % Return back the search agents that go beyond the boundaries of the search space
        Flag4ub=Positions（i:）>ub;
        Flag4lb=Positions（i:）        Positions（i:）=（Positions（i:）.*（~（Flag4ub+Flag4lb）））+ub.*Flag4ub+lb.*Flag4lb;
        
        % Calculate objective function for each search agent
        fitness=fobj（Positions（i:））;
        FEs=FEs+1;
        % Update the leader
        if fitness for maximization problem
            Leader_score=fi