模擬退火算法求解TSP問(wèn)題Matlab源碼

《模擬退火算法求解TSP問(wèn)題Matlab源碼》由會(huì)員分享，可在線(xiàn)閱讀，更多相關(guān)《模擬退火算法求解TSP問(wèn)題Matlab源碼（3頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

1、function [f,T]=TSPSA(d,t0,tf) %TSP問(wèn)題（貨郎擔(dān)問(wèn)題，旅行商問(wèn)題）的模擬退火算法通用malab源程序 % f目標(biāo)最優(yōu)值,T最優(yōu)路線(xiàn),d距離矩陣,t0初始溫度,tf結(jié)束溫度 [m,n]=size(d); L=100*n; t=t0; pi0=1:n; min_f=0; for k=1：n-1 min_f=min_f+d(pi0(k),pi0(k+1)); end min_f=min_f+d(pi0(n),pi0(1)); p_min=pi0; while t>tf for k=1:L; kk=rand; [d_f,pi_1]=exc

2、hange_2(pi0,d); r_r=rand; if d_f<0 pi0=pi_1; elseif exp(d_f/t)>r_r pi0=pi_1; else pi0=pi0; end end f_temp=0; for k=1：n-1 f_temp=f_temp+d(pi0(k),pi0(k+1)); end f_temp=f_temp+d(pi0(n),pi0(1)); if min_f>f_temp min_f=f_temp; p_min=pi0; end t=0.87*t; end f=min_f; T=p_min; %aiwa要調(diào)用

3、的子程序,用于產(chǎn)生新解 function [d_f,pi_r]=exchange_2(pi0,d) [m,n]=size(d); clear m; u=rand; u=u*(n-2); u=round(u); if u<2 u=2; end if u>n-2 u=n-2; end v=rand; v=v*(n-u+1); v=round(v); if v<1 v=1; end v=u+v; if v>n v=n; end pi_1(u)=pi0(v); pi_1(v)=pi0(u); if u>1 for k=1：u-1 pi_1(k)=p

4、i0(k); end end if v>(u+1) for k=1：v-u-1 pi_1(u+k)=pi0(v-k); end end if v