Matlab实现蚁群算法详解及实例

function [y,val]=QACStic
load att48 att48;
MAXIT=300;       % 最大循环次数
NC=48;           % 城市个数
tao=ones(48,48);% 初始时刻各边上的信息最为1
rho=0.2;         % 挥发系数
alpha=1;
beta=2;
Q=100;
mant=20;         % 蚂蚁数量
iter=0;          % 记录迭代次数
for i=1:NC       % 计算各城市间的距离
    for j=1:NC
        distance(i,j)=sqrt((att48(i,2)-att48(j,2))^2+(att48(i,3)-att48(j,3))^2);
    end
end
bestroute=zeros(1,48); % 用来记录最优路径
routelength=inf;        % 用来记录当前找到的最优路径长度
% for i=1:mant          % 确定各蚂蚁初始的位置
% end
for ite=1:MAXIT
    for ka=1:mant       %考查第K只蚂蚁
        deltatao=zeros(48,48); % 第K只蚂蚁移动前各边上的信息增量为零
        [routek,lengthk]=travel(distance,tao,alpha,beta);
        if lengthk<routelength   % 找到一条更好的路径
            routelength=lengthk;
            bestroute=routek;
        end
        for i=1:NC-1      % 第K只蚂蚁在路径上释放的信息量
            deltatao(routek(i),routek(i+1))=deltatao(routek(i),routek(i+1))+Q/lengthk;