簡(jiǎn)單蟻群算法求解TSP的源程序[matlab].doc

簡(jiǎn)單蟻群算法求解TSP的源程序（原創(chuàng)）蟻群算法是新興的仿生算法，最初是由意大利學(xué)者Dorigo M于1991年首次提出，由于具有較強(qiáng)的魯棒性，優(yōu)良的分布式計(jì)算機(jī)制和易于與其它方法結(jié)合等優(yōu)點(diǎn)，成為人工智能領(lǐng)域的一個(gè)研究熱點(diǎn)。本程序是實(shí)現(xiàn)簡(jiǎn)單的蟻群算法，TSP問(wèn)題取的是att48，可從http:/www.iwr.uni-heidelberg.de/groups/comopt/software/TSPLIB95獲取，程序運(yùn)行時(shí)間可能會(huì)比較長(zhǎng)，在我的這臺(tái)CPU 1.6G+內(nèi)存256M的機(jī)器上運(yùn)行時(shí)間大概是13分鐘左右。我用的語(yǔ)言是MATLAB 7.1。此程序僅供學(xué)習(xí)所用，如有問(wèn)題請(qǐng)反饋。謝謝。（注：程序沒有計(jì)算最后一個(gè)城市回來(lái)起點(diǎn)城市的距離）function y,val=QACSticload att48 att48;MAXIT=300; % 最大循環(huán)次數(shù)NC=48; % 城市個(gè)數(shù)tao=ones(48,48);% 初始時(shí)刻各邊上的信息最為1rho=0.2; % 揮發(fā)系數(shù)alpha=1;beta=2;Q=100;mant=20; % 螞蟻數(shù)量iter=0; % 記錄迭代次數(shù)for i=1:NC % 計(jì)算各城市間的距離 for j=1:NC distance(i,j)=sqrt(att48(i,2)-att48(j,2)2+(att48(i,3)-att48(j,3)2); endendbestroute=zeros(1,48); % 用來(lái)記錄最優(yōu)路徑routelength=inf; % 用來(lái)記錄當(dāng)前找到的最優(yōu)路徑長(zhǎng)度% for i=1:mant % 確定各螞蟻初始的位置% endfor ite=1:MAXIT for ka=1:mant %考查第K只螞蟻 deltatao=zeros(48,48); % 第K只螞蟻移動(dòng)前各邊上的信息增量為零 routek,lengthk=travel(distance,tao,alpha,beta); if lengthk=chp tabuk=tabuk j; val=val+distance(np,j); break; end endendy=tabuk;function y=isin(x,A) % 判斷數(shù) x 是否在向量 A 中，如在返回 1 ，否則返回 0