%%?I.?清空環境變量
clear?all
clc

%%?II.?導入數據
load?citys_data.mat

%%?III.?計算城市間相互距離
n?=?size（citys1）;
D?=?zeros（nn）;
for?i?=?1:n
????for?j?=?1:n
????????if?i?~=?j
????????????D（ij）?=?sqrt（sum（（citys（i:）?-?citys（j:））.^2））;
????????else
????????????D（ij）?=?1e-4;??????
????????end
????end????
end

%%?IV.?初始化參數
m?=?50;??????????????????????????????%?螞蟻數量
alpha?=?1;???????????????????????????%?信息素重要程度因子
beta?=?5;????????????????????????????%?啟發函數重要程度因子
rho?=?0.1;???????????????????????????%?信息素揮發因子
Q?=?1;???????????????????????????????%?常系數
Eta?=?1./D;??????????????????????????%?啟發函數
Tau?=?ones（nn）;?????????????????????%?信息素矩陣
Table?=?zeros（mn）;??????????????????%?路徑記錄表
iter?=?1;????????????????????????????%?迭代次數初值
iter_max?=?200;??????????????????????%?最大迭代次數?
Route_best?=?zeros（iter_maxn）;??????%?各代最佳路徑???????
Length_best?=?zeros（iter_max1）;?????%?各代最佳路徑的長度??
Length_ave?=?zeros（iter_max1）;??????%?各代路徑的平均長度??

%%?V.?迭代尋找最佳路徑
while?iter?<=?iter_max
?????%?隨機產生各個螞蟻的起點城市
??????start?=?zeros（m1）;
??????for?i?=?1:m
??????????temp?=?randperm（n）;
??????????start（i）?=?temp（1）;
??????end
??????Table（:1）?=?start;?
??????citys_index?=?1:n;
??????%?逐個螞蟻路徑選擇
??????for?i?=?1:m
??????????%?逐個城市路徑選擇
?????????for?j?=?2:n
?????????????tabu?=?Table（i1:（j?-?1））;???????????%?已訪問的城市集合（禁忌表）
?????????????allow_index?=?~ismember（citys_indextabu）;
?????????????allow?=?citys_index（allow_index）;??%?待訪問的城市集合
?????????????P?=?allow;
?????????????%?計算城市間轉移概率
?????????????for?k?=?1:length（allow）
?????????????????P（k）?=?Tau（tabu（end）allow（k））^alpha?*?Eta（tabu（end）allow（k））^beta;
?????????????end
?????????????P?=?P/sum（P）;
?????????????%?輪盤賭法選擇下一個訪問城市
?????????????Pc?=?cumsum（P）;?????
????????????target_index?=?find（Pc?>=?rand）;?
????????????target?=?allow（target_index（1））;
????????????Table（ij）?=?target;
?????????end
??????end
??????%?計算各個螞蟻的路徑距離
??????Length?=?zeros（m1）;
??????for?i?=?1:m
??????????Route?=?Table（i:）;
??????????for?j?=?1:（n?-?1）
??????????????Length