蚁群算法的原理及其代码

%%来自别人的研究生数学建模的一片论文，侵权即删
clear
clc
t0 = clock;
% P01=xlsread('data.xlsx','E70:F72');
% P01=[72.3 82.3 85.6 82.3 85.6 85.6 83.4 85.6;
% 64.7 66 71.2 67.3 58.2 60.8 59.5 55.6]';
P01=[1 1 1 2 2 3 4 4 5 6 7 7 9 9 10 11 13 13;
4 10 16 2 14 7 4 10 16 13 11 8 15 5 11 16 14 11]';
%P01是点的（x,y)的坐标信息,找出了最短路径，还是TSP问题
%-------------------------------------------------------------------------%
n=size(P01,1);
D=zeros(n,n);
for i=1:n
for j=1:n
if i~=j
D(i,j)=sqrt( sum( ( P01(i,:)-P01(j,:) ).^2 ) );
else
D(i,j)=1e-4;
end
end
end
%-------------------------------------------------------------------------%
%%初始化参数
m = 100; % 蚂蚁数量
alpha = 1; % 信息素重要程度因子
beta = 5; % 启发函数重要程度因子
vol = 0.2; % 信息素挥发(volatilization)因子
Q = 10; % 常系数
Heu_F = 1./D; % 启发函数(heuristic function)
Tau = ones(n,n); % 信息素矩阵
Table = zeros(m,n); % 路径记录表
iter = 1; % 迭代次数初值
iter_max = 200; % 最大