《模拟退火算法求解TSP问题Matlab源码》由会员分享，可在线阅读，更多相关《模拟退火算法求解TSP问题Matlab源码（3页珍藏版）》请在金锄头文库上搜索。
1、function f,T=TSPSA(d,t0,tf)%TSP问题（货郎担问题，旅行商问题）的模拟退火算法通用malab源程序% f目标最优值,T最优路线,d距离矩阵,t0初始温度,tf结束温度m,n=size(d);L=100*n;t=t0;pi0=1:n;min_f=0;for k=1：n-1min_f=min_f+d(pi0(k),pi0(k+1);endmin_f=min_f+d(pi0(n),pi0(1);p_min=pi0;while ttffor k=1:L;kk=rand;d_f,pi_1=exchange_2(pi0,d);r_r=rand; if d_fr_rpi0=pi_
2、1;elsepi0=pi0;endendf_temp=0;for k=1：n-1f_temp=f_temp+d(pi0(k),pi0(k+1);endf_temp=f_temp+d(pi0(n),pi0(1);if min_ff_tempmin_f=f_temp;p_min=pi0;endt=0.87*t;endf=min_f;T=p_min;%aiwa要调用的子程序,用于产生新解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 un-2u=n-2;endv=rand;v=
3、v*(n-u+1);v=round(v);if vnv=n;endpi_1(u)=pi0(v);pi_1(v)=pi0(u);if u1for k=1：u-1pi_1(k)=pi0(k);endendif v(u+1)for k=1：v-u-1pi_1(u+k)=pi0(v-k);endendif vnfor k=(v+1):npi_1(k)=pi0(k);endendd_f=0;if vnd_f=d(pi0(u-1),pi0(v)+d(pi0(u),pi0(v+1);for k=(u+1):nd_f=d_f+d(pi0(k),pi0(k-1);endd_f=d_f-d(pi0(u-1),pi0(u)-d(pi0(v),pi0(v+1);for k=(u+1):nd_f=d_f-d(pi0(k-1),pi0(k);endelsed_f=d(pi0(u-1),pi0(v)+d(pi0(u),pi0(1)-d(pi0(u-1),pi0(u)-d(pi0(v),pi0(1);for k=(u+1):nd_f=d_f+d(pi0(k),pi0(k-1);endfor k=(u+1):nd_f=d_f-d(pi0(k-1),pi0(k);endendpi_r=pi_1;
