function?[R_bestL_bestL_aveShortest_RouteShortest_Length]=ACATSP（CNC_maxmAlphaBetaRhoQ）
%%-------------------------------------------------------------------------
%%?主要符號說明
%%?C?n個城市的坐標，n×2的矩陣
%%?NC_max?最大迭代次數
%%?m?螞蟻個數
%%?Alpha?表征信息素重要程度的參數
%%?Beta?表征啟發式因子重要程度的參數
%%?Rho?信息素蒸發系數
%%?Q?信息素增加強度系數
%%?R_best?各代最佳路線
%%?L_best?各代最佳路線的長度
%%?=========================================================================
%%第一步：變量初始化
n=size（C1）;%n表示問題的規模（城市個數）
D=zeros（nn）;%D表示完全圖的賦權鄰接矩陣
for?i=1:n
for?j=1:n
if?i~=j
D（ij）=（（C（i1）-C（j1））^2+（C（i2）-C（j2））^2）^0.5;
else
D（ij）=eps;??????%i=j時不計算，應該為0，但后面的啟發因子要取倒數，用eps（浮點相對精度）表示
end
D（ji）=D（ij）;???%對稱矩陣
end
end
Eta=1./D;??????????%Eta為啟發因子，這里設為距離的倒數
Tau=ones（nn）;?????%Tau為信息素矩陣
Tabu=zeros（mn）;???%存儲并記錄路徑的生成
NC=1;???????????????%迭代計數器，記錄迭代次數
R_best=zeros（NC_maxn）;???????%各代最佳路線
L_best=inf.*ones（NC_max1）;???%各代最佳路線的長度
L_ave=zeros（NC_max1）;????????%各代路線的平均長度
?
while?NC<=NC_max????????%停止條件之一：達到最大迭代次數，停止

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????4462??2015-09-19?20:12??ACATSP.m
?????文件?????????600??2015-09-19?15:15??coords.mat
?????文件????????1058??2015-09-19?18:00??Run.m