initial.m

用于解决TSP问题的蚁群算法

%初始化函数
ant_n=50;
city_n=30;
NC=200; %迭代次数
tobu(ant_n,city_n)=0; %蚂蚁禁忌表
ph_table(city_n,city_n)=0;  %信息素分布表
hu_table(city_n,city_n)=0;  %启发式分布表
a=3;       %信息指数
b=1;       %启发式指数
dispose=0.2; %信息素挥发率
everbest=inf; %历代最优解
adapt_ave(NC)=0; %进化数据
adapt_best(NC)=inf;%历代最优解的变化

%计算城市距离矩阵
for i=1:city_n
    for j=1:city_n
        x1=pos(i,1);
        y1=pos(i,2);
        x2=pos(j,1);
        y2=pos(j,2);
        
        dis_table(i,j)=sqrt((x1-x2)^2+(y1-y2)^2);
    end
end
for i=1:city_n
    dis_table(i,i)=inf;
end

%初始化启发式分布表
for i=1:city_n
    for j=1:city_n
        hu_table(i,j)=100/dis_table(i,j);
    end
end

%初始化信息素
for i=1:city_n
    for j=1:city_n
        ph_table(i,j)=0.5;    %额定信息素取值为0~1
    end
end