search.m

来自「一个ACO实验程序」· M 代码 · 共 42 行

%解搜索函数

%初始化禁忌表
for i=1:ant_n
    for j=1:city_n
        tobu(i,j)=0;
    end
end
%随机放置蚂蚁初始位置
for i=1:ant_n
    city=round((city_n-1)*rand)+1;  %随机产生初始位置城市
    tobu(i,1)=city;
end

for i=1:ant_n  %对每一只蚂蚁进行解构造
    n=1;
        
    for j=1:city_n-1  %每只蚂蚁进行city_n-1次探索
        
        for k=1:city_n      %临时池供选择用
            temp_pool(k)=(hu_table(tobu(i,n),k)^b)*(ph_table(tobu(i,n),k)^a);
        end
        for k=1:n           %将已经访问过的城市的选择概率置0
            temp_pool(tobu(i,k))=0;
        end

        ada_sum=0;
        for k=1:city_n
            ada_sum=ada_sum+temp_pool(k);
        end
        r=rand*ada_sum;  %随机产生一个数
        ada_temp=0;      %初始化累加值为0
        h=0;
        while(ada_temp<r)
            h=h+1;
            ada_temp=ada_temp+temp_pool(h);        
        end
        %退出循环时的h值即为被选择的个体序号
        n=n+1;
        tobu(i,n)=h;    %当前访问的城市加入禁忌表
    end
end