五個城市的TSP問題，距離矩陣為：

01015621008139158020156132005291550設初始溫度為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));endmin_f=min_f+d(pi0(n),pi0(1));p_min=pi0;while t>tf for k=1:L kk=rand; [d_f,pi_1]=exchange_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;endf=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 u<2 u=2;endif u>n-2 u=n-2;endv=rand;v=v*(n-u+1);v=round(v);if v<1 v=1;endv=u+v;if v>n v=n;endpi_1(u)=pi0(v);pi_1(v)=pi0(u);if u>1 for k=1:u-1 pi_1(k)=pi0(k); endendif v>u+1 for k=1:(v-u-1) pi_1(u+k)=pi0(v-k); endendif v

總結

以上是生活随笔為你收集整理的tsp matlab,五个城市的TSP问题MATLAB程序.doc的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。