蚁群算法和简要matlab来源
從1991由意大利學者 M. Dorigo，V. Maniezzo 和 A. Colorni 通過模擬蟻群覓食行為提出了一種基于群體的模擬進化算法——蟻群優化。極大關注，蟻群算法的特點：
??? ① 其原理是一種正反饋機制或稱增強型學習系統; 它通過【最優路徑上螞蟻數量的添加→信息素強度添加→后來螞蟻選擇概率增大→最優路徑上螞蟻數量更大添加】達到終于收斂于最優路徑上L
② 它是一種通用型隨機優化方法, 它吸收了螞蟻的行為特(內在搜索機制) , 它是使用人工螞蟻仿真(也稱螞蟻系統) 來求解問題L但人工螞蟻決不是對實際螞蟻的一種簡單模擬, 它融進了人類的智能L人工螞蟻有一定的記憶; 人工螞蟻不全然是瞎的; 人工螞蟻生活的時空是離散的L
③ 它是一種分布式的優化方法, 不僅適合眼下的串行計算機, 并且適合未來的并行計算機L
④ 它是一種全局優化的方法, 不僅可用于求解單目標優化問題, 并且可用于求解多目標優化問題L
⑤ 它是一種啟示式算法, 計算復雜性為o (Nc*n2*m) , 當中Nc 是迭代次數, m 是螞蟻數目, n 是目的節點數目L
蟻群發現最短路徑的原理和機制[1]
以下用圖 1解釋蟻群發現最短路徑的原理和機制。
如圖 1(a)所看到的。在蟻巢和食物源之間有兩條道路 Nest-A-B-D-Food 和Nest-A-C-D-Food，其長度分別為 4 和 6。單位時間內螞蟻可移動一個單位長度的距離。開始時全部路徑上都沒有外激素。
如圖 1(b)，在 t=0 時刻。20 僅僅螞蟻從蟻巢出發移動到 A。因為路徑上沒有外激素，它們以同樣概率選擇左側或右側道路。因此平均有 10 僅僅螞蟻走左側，另外 10 僅僅走右側。
如圖 1(c)，在 t=4 時刻。第一組先到達食物源的螞蟻將折回。
如圖 1(d)，在 t=5 時刻。兩組螞蟻將在 D 點相遇。
此時 BD 上的外激素數量與 CD 上的同樣。因此返回的 10 僅僅螞蟻中有 5 僅僅選擇 BD 而另 5 僅僅選擇 CD。
如圖 1(e)，在 t=8 時刻，前 5 個螞蟻將返回巢穴，而在 AC、CD 和 AB 上各有 5 個螞蟻。
如圖 1(f)，在 t=9 時刻。前 5 個螞蟻又回到 A 而且再次面對往左還是往右的選擇。
這時。AB 上的軌跡數是 20 而 AC 上是 15。因此將有較為多數的螞蟻選擇往右，從而增強了 AB 上外激素的量。隨著該過程的繼續。兩條道路上外激素數量的差距將越來越大，直至絕大多數螞蟻都選擇了最短的路徑。
正是因為一條道路要比還有一條道路短，因此，在同樣的時間間隔內。短的路線會有很多其它的機會被選擇。
依據仿生學家的研究結果,螞蟻憑借路徑尋優的能力可以找到蟻巢與食物之間的最短路徑,其原理在于:螞蟻在所經過的路徑上留下一種揮發性分泌物(pheromone,下面稱為信息素),信息素隨著時間的推移會逐漸揮發消失.螞蟻在覓食過程中可以感知這樣的物質的存在及其強度,并以此來指導自己的運動方向,傾向于朝著這樣的物質強度高的方向移動,即選擇該路徑的概率與當時這條路徑上該物質的強度成正比.信息素強度越高的路徑,選擇它的螞蟻就越多,則在該路徑上留下的信息素的強度就更大,而強度大的信息素又吸引很多其它的螞蟻,從而形成一種正反饋.通過這樣的正反饋,螞蟻終于可以發現最佳路徑,導致大部分的螞蟻都會走此路徑.
以求解ｎ個城市的ＴＳＰ旅行商問題為例說明ＡＣＡ模型.
設蟻群中螞蟻的數量為ｍ,ｄij (ｉ,ｊ=1,2,…,ｎ)表示城市ｉ和城市ｊ之間的距離,ｂi(ｔ)表示ｔ時刻位于城市ｉ的螞蟻的個數,則有 表示ｔ時刻在城市ｉ,ｊ連線上殘留的信息量.初始時刻,各條路徑上信息量相等,設τｉｊ(0)=Ｃ(Ｃ為常數).螞蟻ｋ(ｋ=1,2,…,ｍ)在運動過程中,依據各條路徑上的信息量決定轉移方向. 表示在ｔ時刻螞蟻ｋ由城市ｉ轉移到城市ｊ的概率.
???? （1）
殘留信息的重要程度;β——啟示信息的重要程度;tabuk——記錄螞蟻ｋ當前所走過的城市,稱為記憶列表,ｋ=1,2,…,ｍ,集合tabuk隨著進化過程作動態調整.經過ｎ個時刻,全部螞蟻都完畢了一次遍歷.此時,計算每一僅僅螞蟻所走過的路徑Ｌｋ,并保存最短路徑Ｌｍｉｎ=ｍｉｎ{Ｌｋ︱ｋ=1,2,…,ｍ}.在螞蟻完畢一次循環以后,各路徑上的信息量進行例如以下調整
τｉｊ(ｔ+1)=(1-ρ)τｉｊ(ｔ)+Δτｉｊ??? (2)
式中ρ∈(0,1),表示信息素τｉｊ(ｔ)隨時間的推移而衰減的程度.所以1-ρ為信息素殘留因子，開始時Δτｉｊ（0）=0，
信息素增量Δτｉｊ可表示???????????? (3)
式中Δτｋｉｊ為螞蟻ｋ在本次循環中在城市ｉ和ｊ之間留下的信息量,它的計算公式依據詳細問題而定.Dorigo曾給出Δτｋｉｊ3種不同的模型,分別稱為Ant-Cycle模型、Ant-Quantity模型、Ant-Density模型，它們的區別就在于信息素的更新機制,即其區別在于Δτｋｉｊ
在Ant-Cycle模型中：
??? （4） 式中，Q表示信息素強度。它在一定程度上影響算法的收斂速度；Lk表示第K僅僅螞蟻在本次循環中所奏路徑的總長度。
在Ant-Quantity模型中：
（5） 式中，Q表示信息素強度。它在一定程度上影響算法的收斂速度；dij表示第K僅僅螞蟻在t和t+1之間經過的（ i, j ）
在Ant-Density模型中：
（6） 差別:式(5)式(6)中利用的是局部信息，即螞蟻完畢一步后更新路徑上的信息素；而式(4)中利用的是總體信息，即螞蟻完畢一個循環后全部路徑上的信息素。經過大量試驗總結研究。採用式(4)性能較好。所以 Ant-Cycle模型是最優的。
以上說明了信息素殘留因子1-ρ、信息啟示式因子α、期望啟示式因子β、信息素強度Q、螞蟻數目M等都是很重要的參數，其選區方式和選區原則直接影響到蟻群算法的全局收斂性和求解效率。我們學習到這樣的“三步走”[2]選擇蟻群算法最優組合參數的有效方法：
（1） 確定螞蟻數目M，依據 城市規模 / 螞蟻數目 ≈1.5的選擇策略來確定螞蟻的總數目。
（2） 參數粗調，即調整數值范圍較大的信息啟示式因子α、期望啟示式因子β、信息素強度Q等參數。已得到較理想的解。
（3） 參數微調。即調整數值范圍較小的信息素殘留因子1-ρ。
2 眼下蟻群算法的應用
盡管對蟻群算法的研究時間不長, 可是初步研究已顯示出它在求解復雜優化問題方面具有非常大的優勢, 特別是1998 年在比利時布魯塞爾專門召開了第一屆螞蟻優化國際研討會后, 如今每兩年召開一次這種螞蟻優化國際研討會。這標志著蟻群算法的研究已經得到了國際上的廣泛支持。使得這種新興的智能進化仿生算法展現出了勃勃生機[3]。
以蟻群算法為代表的群體智能已成為當今分布式人工智能研究的一個熱點,很多源于蜂群和蟻群模型設計的算法已越來越多地被用于企業的運轉模式的研究。
美國五角大樓正在資助關于群體智能系統的研究工作--群體戰略(SWARM STRATEGY),它的一個實戰用途是通過運用成群的空中無人駕駛飛行器和地面車輛來轉移敵人的注意力,讓自己的軍隊在敵人后方不被察覺地安全行進。英國電信公司和美國世界通信公司以電子螞蟻為基礎,對新的電信網絡管理方法進行了試驗。群體智能還被應用于工廠生產計劃的制定和運輸部門的后勤管理。
美國太平洋西南航空公司採用了一種直接源于螞蟻行為研究成果的運輸管理軟件,結果每年至少節約了1000萬美元費用開支。英國聯合利華公司已領先利用群體智能技術改善其一家牙膏廠的運轉狀況。美國通用汽車公司,法國液氣公司,荷蘭公路交通部和美國一些移民事務機構也都採用這樣的技術來改善其運轉的機能。又如美國MCIWorld.com公司一直研究人工螞蟻，并用于管理公司的電話網，對用戶記賬收費等工作。
另外。還設計“人工螞蟻”打算用于因特網的路由管理。鑒于群體智能廣闊的應用前景,美國和歐洲聯盟均于近幾年開始出資資助基于群體智能模擬的相關研究項目, 關在一些院校開設群體智能的相關課程.牛津大學出版社1999年版的E.Bonabeau和M.Dorigo等人編寫的專著《群體智能:從自然到人工系統》(Swarm Intelligence:From Natural to Artificial System),以及2001年出版的J.Kennedy和R.Eberhart編著的《群體智能》(Swarm Intelligence)進一步擴大了群體智能的影響.IEEE進化計算會刊也于2002年8月出版了蟻群優化算特刊。國內也有研究者用螞蟻算法求解全國144個城市的最短回路問題，求得的解同其他方法求到得解一樣精確，這說明螞蟻算法不可是求解組合優化問題的可行方法。并且是一種非常有競爭力的算法。國家自然科學基金"十五"期間學科交叉類優先資助領域中的認知科學及其信息處理的研究內容中也明白列出了群體智能領域的進化,自適應與現場認知主題[4]。并且從1999年開始，差點兒每年都會有幾項相關項目獲得資助。蟻群算法是一種新型的模擬進化算法，其在數據挖掘中的應用正逐步引起人們的關注。眼下。人工蟻群在知識發現的過程中主要用于發掘聚類模型和分類模型。
2．1蟻群算法在數據挖掘中的應用
聚類是將一組對象分成若干個群體，每一個群體構成一個簇，使得簇內的對象盡可能具有最大的相似性。不同簇之間的對象盡可能有最大的相異性。
眼下，聚類方法主要有K均值法，模糊聚類、神經網絡聚類、基于遺傳算法的聚類、小波變換聚類以及將這些算法有效結合而形成的改進方法。隨著蟻群算法研究的興起。人們發如今某些方面採用蟻群模型進行聚類更加接近實際的聚類問題。
將蟻群算法用于聚類分析,靈感源于螞蟻堆積他們的尸體和分類他們的幼體。
基于蟻群算法的聚類方法從原理上可分為兩種：一種是基于蟻堆形成原理來實現數據聚類，還有一種是運用螞蟻覓食的原理，利用信息來實現聚類分析。
而數據是數據挖掘的還有一個重要主題，它是在數據庫對象集合中尋找屬性，并依據分類模式將其劃分為不同類別的過程。分類過程利用歷史數據記錄自己主動推導出對給定數據的分類樹。分類器構造方法有統計學方法、機器學習法、神經網絡、決策樹等。從知識發現的觀點來看，分類規則的表達方式形如if<條件>then<類>規則前件（if 部分）包括一組條件集合，一般由邏輯連接符連接；規則結論（then部分）定義了樣本的預測類，這些樣本的預測屬性滿足規則前件所定義的全部條件[5]。
將蟻群算法引入分類規則的發現。是利用蟻群覓食原理在數據庫中進行搜索，對隨機產生的一組規則進行選擇優化。直到數據庫能被該組規則覆蓋，從而挖掘出隱含在數據庫中的規則。建立最優的分類模型。蟻群算法搜索的初始條件為發現規則的集合為空。且訓練集包括全部的訓練樣本。螞蟻搜索一次要完畢規則生成、規則剪枝、信息素更新三個任務。一次搜索生成一條規則，而且將這條規則增加發現規則集合。同一時候將該條規則所覆蓋的訓練樣本從訓練集中刪除。假設未覆蓋訓練樣本的數目大于用戶定義的閾值。即最大未覆蓋樣本數。就重復運行上述過程，終于算法將得到一組最優分類規則集合[5]。?
最早在這一領域開展工作的是Deneubourg 等[6],他們依據數據對象與其周圍對象的相似性,讓螞蟻隨機地移動、拾起或放下數據對象,以達到聚類數據的目的,這個基本模型已成功地應用于機器人領域。Lumer 等首先改進此算法,提出了LF算法。Wu 等、Ramos等、Yang等[7]從不同角度對LF算法進行了改進,在用蟻群算法進行聚類分析方面取得了一定成效。近幾年，學者在這方面的研究從來沒有間斷過。也取得了一定的研究成果。
2．2 結論?
只是。將蟻群算法運用于數據發掘還存在一些問題，須要進一步研究：
（1）怎樣將現實的挖掘任務轉換成蟻群求解的問題空間，并用適當的方式表達。怎樣定義“人工螞蟻”以及螞蟻間的非直接通信方式（如路徑上的信息素、對象的分布狀態等）的選擇。
（2）怎樣建立正反饋機制，定義啟示函數，遞增地進行問題求解。而且使得到的解與問題定義中現實世界的情況相相應。
（3）基于蟻群的算法要初始化大量的參數。這些參數的選擇會對算法的性能產生較大的影響。但其選取的方法和原則眼下尚無理論上的根據。僅僅能通過多次實驗調優，因此參數的最佳設置原則還有待進一步研究。
（4）蟻群算法的搜索時間較長。怎樣將蟻群算法與遺傳算法、免疫算法等優化算法相結合。改善和提高算法性能。以適應海量數據庫的知識發現。
所以怎樣在數據挖掘中運用蟻群算法高速、高效地獲得高質量的知識越來越受到人們的關注。逐漸成為最近的研究熱點[5]。
?
下面是解放軍信息project大學一個老師編的matlab程序。請尊重原作者勞動，引用時請注明出處。
我經過改動添加了凝視。已經執行過。無誤，
?
function [R_best,L_best,L_ave,Shortest_Route,Shortest_Length]=ACATSP(C,NC_max,m,Alpha,Beta,Rho,Q)
%%-------------------------------------------------------------------------
%% 主要符號說明
%% C n個城市的坐標，n×2的矩陣
%% NC_max 最大迭代次數
%% m 螞蟻個數
%% Alpha 表征信息素重要程度的參數
%% Beta 表征啟示式因子重要程度的參數
%% Rho 信息素蒸發系數
%% Q 信息素添加強度系數
%% R_best 各代最佳路線
%% L_best 各代最佳路線的長度
%%=========================================================================
%%第一步：變量初始化
n=size(C,1);%n表示問題的規模（城市個數）
D=zeros(n,n);%D表示全然圖的賦權鄰接矩陣
for i=1:n
for j=1:n
if i~=j
D(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時不計算，應該為0，但后面的啟示因子要取倒數。用eps（浮點相對精度）表示
end
D(j,i)=D(i,j);?? %對稱矩陣
end
end
Eta=1./D;????????? %Eta為啟示因子，這里設為距離的倒數
Tau=ones(n,n);???? %Tau為信息素矩陣
Tabu=zeros(m,n);?? %存儲并記錄路徑的生成
NC=1;?????????????? %迭代計數器，記錄迭代次數
R_best=zeros(NC_max,n);?????? %各代最佳路線
L_best=inf.*ones(NC_max,1);?? %各代最佳路線的長度
L_ave=zeros(NC_max,1);??????? %各代路線的平均長度
while NC<=NC_max??????? %停止條件之中的一個：達到最大迭代次數，停止
%%第二步：將m僅僅螞蟻放到n個城市上
Randpos=[];?? %隨即存取
for i=1:(ceil(m/n))
Randpos=[Randpos,randperm(n)];
end
Tabu(:,1)=(Randpos(1,1:m))';??? %此句不太理解？
%%第三步：m僅僅螞蟻按概率函數選擇下一座城市。完畢各自的周游
for j=2:n???? %所在城市不計算
for i=1:m????
visited=Tabu(i,1:(j-1)); %記錄已訪問的城市，避免反復訪問
J=zeros(1,(n-j+1));?????? %待訪問的城市
P=J;????????????????????? %待訪問城市的選擇概率分布
Jc=1;
for k=1:n
if length(find(visited==k))==0?? %開始時置0
J(Jc)=k;
Jc=Jc+1;???????????????????????? %訪問的城市個數自加1
end
end
%以下計算待選城市的概率分布
for k=1:length(J)
P(k)=(Tau(visited(end),J(k))^Alpha)*(Eta(visited(end),J(k))^Beta);
end
P=P/(sum(P));
%按概率原則選取下一個城市
Pcum=cumsum(P);???? %cumsum，元素累加即求和
Select=find(Pcum>=rand); %若計算的概率大于原來的就選擇這條路線
to_visit=J(Select(1));
Tabu(i,j)=to_visit;
end
end
if NC>=2
Tabu(1,:)=R_best(NC-1,:);
end
%%第四步：記錄本次迭代最佳路線
L=zeros(m,1);???? %開始距離為0。m*1的列向量
for i=1:m
R=Tabu(i,:);
for j=1:(n-1)
L(i)=L(i)+D(R(j),R(j+1));??? %原距離加上第j個城市到第j+1個城市的距離
end
L(i)=L(i)+D(R(1),R(n));????? %一輪下來后走過的距離
end
L_best(NC)=min(L);?????????? %最佳距離取最小
pos=find(L==L_best(NC));
R_best(NC,:)=Tabu(pos(1),:); %此輪迭代后的最佳路線
L_ave(NC)=mean(L);?????????? %此輪迭代后的平均距離
NC=NC+1????????????????????? %迭代繼續
%%第五步：更新信息素
Delta_Tau=zeros(n,n);??????? %開始時信息素為n*n的0矩陣
for i=1:m
for j=1:(n-1)
Delta_Tau(Tabu(i,j),Tabu(i,j+1))=Delta_Tau(Tabu(i,j),Tabu(i,j+1))+Q/L(i);??????????
%此次循環在路徑（i，j）上的信息素增量
end
Delta_Tau(Tabu(i,n),Tabu(i,1))=Delta_Tau(Tabu(i,n),Tabu(i,1))+Q/L(i);
%此次循環在整個路徑上的信息素增量
end
Tau=(1-Rho).*Tau+Delta_Tau; %考慮信息素揮發。更新后的信息素
%%第六步：禁忌表清零
Tabu=zeros(m,n);???????????? %%直到最大迭代次數
end
%%第七步：輸出結果
Pos=find(L_best==min(L_best)); %找到最佳路徑（非0為真）
Shortest_Route=R_best(Pos(1),:) %最大迭代次數后最佳路徑
Shortest_Length=L_best(Pos(1)) %最大迭代次數后最短距離
subplot(1,2,1)????????????????? %繪制第一個子圖形
DrawRoute(C,Shortest_Route)???? %畫路線圖的子函數
subplot(1,2,2)????????????????? %繪制第二個子圖形
plot(L_best)
hold on???????????????????????? %保持圖形
plot(L_ave,'r')
title('平均距離和最短距離')???? %標題
function DrawRoute(C,R)
%%=========================================================================
%% DrawRoute.m
%% 畫路線圖的子函數
%%-------------------------------------------------------------------------
%% C Coordinate 節點坐標，由一個N×2的矩陣存儲
%% R Route 路線
%%=========================================================================
N=length(R);
scatter(C(:,1),C(:,2));
hold on
plot([C(R(1),1),C(R(N),1)],[C(R(1),2),C(R(N),2)],'g')
hold on
for ii=2:N
plot([C(R(ii-1),1),C(R(ii),1)],[C(R(ii-1),2),C(R(ii),2)],'g')
hold on
end
title('旅行商問題的最優結果 ')
?
總結
以上是生活随笔為你收集整理的蚁群算法和简要matlab来源的全部內容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： FTP 500 OOPS
- 下一篇： mongo api