先新建一個(gè)主程序M文件ACATSP.m

代碼如下：

function ［R_best，L_best，L_ave，Shortest_Route，Shortest_Length］=ACATSP（C，NC_max，m，Alpha，Beta，Rho，Q）

%%=========================================================================

%% 主要符號(hào)說(shuō)明

%% C n個(gè)城市的坐標(biāo)，n×2的矩陣

%% NC_max 蟻群算法MATLAB程序最大迭代次數(shù)

%% m 螞蟻個(gè)數(shù)%% Alpha 表征信息素重要程度的參數(shù)

%% Beta 表征啟發(fā)式因子重要程度的參數(shù)%% Rho 信息素蒸發(fā)系數(shù)

%% Q 表示蟻群算法MATLAB程序信息素增加強(qiáng)度系數(shù)

%% R_best 各代最佳路線

%% L_best 各代最佳路線的長(zhǎng)度

%% 蟻群算法MATLAB程序第一步：變量初始化

n=size（C，1）;%n表示問(wèn)題的規(guī)模（城市個(gè)數(shù)）

D=zeros（n，n）;%D表示完全圖的賦權(quán)鄰接矩陣

for i=1:n

for j=1:n

if i~=jD（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 時(shí)不計(jì)算，應(yīng)該為0，但后面的啟發(fā)因子要取倒數(shù)，用eps（浮點(diǎn)相對(duì)精度）表示

?

%% 蟻群算法MATLAB程序第二步：將m只螞蟻放到n個(gè)城市上

?

%% 蟻群算法MATLAB程序第三步：m只螞蟻按概率函數(shù)選擇下一座城市，完成各自的周游

?

%% 下面計(jì)算蟻群算法MATLAB程序待選城市的概率分布

?

%% 蟻群算法MATLAB程序第四步：記錄本次迭代最佳路線

?

%% 蟻群算法MATLAB程序第五步：更新信息素

?

%% 蟻群算法MATLAB程序第七步：輸出結(jié)果

?

建立一個(gè)子程序DrawRoute.m

代碼如下

?

需要輸入的參數(shù)數(shù)據(jù)有： C： n個(gè)城市的坐標(biāo)，n×2的矩陣NC_max： 蟻群算法MATLAB程序最大迭代次數(shù)M： 螞蟻個(gè)數(shù)Alpha： 表征信息素重要程度的參數(shù)Beta：表征啟發(fā)式因子重要程度的參數(shù)

Rho： 信息素蒸發(fā)系數(shù)Q：表示蟻群算法MATLAB程序信息素增加強(qiáng)度系數(shù)運(yùn)行時(shí)打開(kāi)ACATSP.m 點(diǎn)擊運(yùn)行或輸入ACATSP（C，NC_max，m，Alpha，Beta，Rho，Q）

一個(gè)運(yùn)行實(shí)例：

m=31;Alpha=1;Beta=5;Rho=0.1;NC_max=200;Q=100;

?
? ? ? ? ? ??
? ? ? ? ? ?