先新建一个主程序 M 文件 ACATSP.m
代码如下 :
function [R_best,L_best,L_ave,Shortest_Route,Shortest_Length]=ACATSP(C,NC_max,m,Alpha,Beta,Rho,Q)
%%================================================================
=========
%% 主要符号说明
%% C n 个城市的坐标, n×2 的矩阵
%% NC_max
蚁群算法 MATLAB程序
最大迭代次数
%% m 蚂蚁个数
%% Alpha 表征信息素重要程度的参数
%% Beta 表征启发式因子重要程度的参数
%% Rho 信息素蒸发系数
%% Q 表示
蚁群算法 MATLAB程序
信息素增加强度系数
%% R_best 各代最佳路线
%% L_best 各代最佳路线的长度
%%================================================================
=========
%% 蚁群算法 MATLAB程序第一步:变量初始化
n=size(C,1);%n 表示问题的规模(城市个数)
D=zeros(n,n);%D 表示完全图的赋权邻接矩阵
for i=1:n
for j=1:n
if i~=j
D(i,j)=((C(i,1)-C(j,1))^2+(C(i,2)-C(j,2))^2)^0.5;
else
D(i,j)=eps; % i = j 时不计算,应该为 0,但后面的启发因子要取
倒数,用 eps (浮点相对精度)表示