www.gusucode.com > 溷沌优化工具箱 - chaoshcacs源码程序 > randInit.m

    function [ants,optIndex,sigma] = randInit(bounds,antnum)
ant = struct('pos',[],'fitness',0);
for i = 1:antnum
    ant.pos =(bounds(:,1) + rand(size(bounds,1),1).*(bounds(:,2)-bounds(:,1)))';
    ant.fitness = evalCandidate(ant.pos);
    ants(i) = ant;
    if i==1
        optIndex = 1;
    elseif ants(optIndex).fitness<ants(i).fitness
        optIndex = i;    
    end
end
sigma =3* max(bounds(:,2) - bounds(:,1));