五个城市的TSP问题MATLAB程序.doc

  1. 1、本文档共4页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
五个城市的TSP问题,距离矩阵为: 0 10 15 6 2 10 0 8 13 9 15 8 0 20 15 6 13 20 0 5 2 9 15 5 0 设初始温度为10,结束温度为0.7,则程序如下: d=[0, 10, 15, 6, 2;10, 0, 8, 13, 9;15, 8, 0, 20, 15;6, 13, 20, 0, 5;2, 9, 15, 5, 0]; t0=10; tf=0.7; [f,T]=trp(d,t0,tf) 得结果如下: f=43 T=1 4 5 2 3 MATLAB 的程序如下: function [f,T]=trp(d,t0,tf) % f为目标函数最优值,T为最优路线,d为距离矩阵,t0为初始温度,tf为结束温度 [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 ttf for k=1:L kk=rand; [d_f,pi_1]=exchange_2(pi0,d); r_r=rand; if d_f0 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_ff_temp min_f=f_temp; p_min=pi0; end t=0.87*t; end f=min_f; T=p_min; %下面的函数产生新解 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 u2 u=2; end if un-2 u=n-2; end v=rand; v=v*(n-u+1); v=round(v); if v1 v=1; end v=u+v; if vn v=n; end pi_1(u)=pi0(v); pi_1(v)=pi0(u); if u1 for k=1:u-1 pi_1(k)=pi0(k); end end if vu+1 for k=1:(v-u-1) pi_1(u+k)=pi0(v-k); end end if vn for k=(v+1):n pi_1(k)=pi0(k); end end d_f=0; if vn d_f=d(pi0(u-1),pi0(v))+d(pi0(u),pi0(v+1)); for k=(u+1):n d_f=d_f+d(pi0(k),pi0(k-1)); end d_f=d_f-d(pi0(u-1),pi0(u))-d(pi0(v),pi0(v+1)); for k=(u+1):n d_f=d_f-d(pi0(k-1),pi0(k)); end else d_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):n d_f=d_f+d(pi0(k),pi0(k-1)); end for k=(u+1):n d_f=d_f-d(pi0(k-1),pi0(k)); end end pi_r=pi_1;

文档评论(0)

xingyuxiaxiang + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档