function?[F]=F（x1x2）?%目標函數
F=-（x1.^2+2*x2.^2-0.3*cos（3*pi*x1）-0.4*cos（4*pi*x2）+0.7）;
End?
function?[maxxmaxymaxvalue]=antcolony%?蟻群算法求函數最大值的程序
ant=200;%螞蟻數量
times=50;%螞蟻移動次數
rou=0.8;%信息素揮發系數
p0=0.2;%轉移概率常數
lower_1=-1;%設置搜索范圍
upper_1=1;%
lower_2=-1;%
upper_2=1;%
for?i=1:ant??
????X（i1）=（lower_1+（upper_1-lower_1）*rand）;%隨機設置螞蟻的初值位置
????X（i2）=（lower_2+（upper_2-lower_2）*rand）;
????tau（i）=F（X（i1）X（i2））;%第i只螞蟻的信息量?end?%隨機初始每只螞蟻的位置
????step=0.05;%網格劃分單位
????f=‘-（x.^2+2*y.^2-0.3*cos（3*pi*x）-0.4*cos（4*pi*y）+0.7）‘;
????[xy]=meshgrid（lower_1:step:upper_1lower_2:step:upper_2）;
????z=eval（f）;%eval函數，將字符串內的內容執行再賦給對象
????figure（1）;
????mesh（xyz）;%網格圖?
????hold?on;
????plot3（X（:1）X（:2）tau‘k*‘）%螞蟻初始位置
????hold?on;
????text（0.10.8-0.1‘螞蟻的初始分部位置‘）
????xlabel（‘x‘）;
????ylabel（‘y‘）;
????zlabel（‘f（xy）‘）;
????for?t=1:times?%第t次移動
????????lamda=1/t;?%步長系數，隨移動次數增大而減少?
????[tau_best（t）bestindex]=max（tau）;%第t次移動的最優值及其位置

?????????for?i=1:ant?%第i只螞蟻?
????????????p（ti）=（tau（bestindex）-tau（i））/tau（bestindex）;%最優值與第i只螞蟻的值的差比
??????????????%計算狀態轉移概率
????????end??
????????for?i=1:ant
????????????if?p（ti）????????????????temp1=X（i1）+（2*rand-1）*lamda;%移動