C=[1?2;10?3;4?4;1?7;2?4];NC_max=100;m=3;Alpha=1.5;Beta=1.5;Rho=0.9;Q=1;
%%-------------------------------------------------------------------
%%?主要符號(hào)說(shuō)明

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

%%?NC_max?最大迭代次數(shù)

%%?m?螞蟻個(gè)數(shù)

%%?Alpha?表征信息素重要程度的參數(shù)

%%?Beta?表征啟發(fā)式因子重要程度的參數(shù)

%%?Rho?信息素蒸發(fā)系數(shù)

%%?Q?信息素增加強(qiáng)度系數(shù)

%%?R_best?各代最佳路線(xiàn)

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

%%
%%第一步：變量初始化

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

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

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

end

D（ji）=D（ij）;???%對(duì)稱(chēng)矩陣

end

end

Eta=1./D;??????????%Eta為啟發(fā)因子，這里設(shè)為距離的倒數(shù)

Tau=ones（nn）;?????%Tau為信息素矩陣

Tabu=zeros（mn）;???%存儲(chǔ)并記錄路徑的生成

NC=1;???????????????%迭代計(jì)數(shù)器，記錄迭代次數(shù)

R_best=zeros（NC_maxn）;