用模拟退火求TSP如何？

慢跑20

一直EXCEL中为33个点之间，任意两个点之间的距离（没有数字的是无穷）。 求经过所有点之间的最短路线问题，。
, j0 ]" ]) [, ^- d9 h2 t! I( [4 G
A=xlsread('33点矩阵s上三角');
>> for i=1:33;
       for j = 1:33
                   if i<j
                 temp(i,j)=A(i,j);
; L; e( R0 c+ n( I' x* G! Y9 `                 A(j,i)=temp(i,j);
                   end   
                 if isnan(A(i,j))
                 A(i,j)=inf;
7 l5 B$ S: |# K. k, f                 end
         
    end
' w3 U+ q. ^. i( b* q, q4 y$ T end
+ G. B2 V" Y; H  H' g

这样A为邻接矩阵了。然后运行百度的代码：
$ {# Z; k- [- Z  Q# P
  u: s3 H3 S5 [1 l, g7 e& i
9 J5 H) D9 [8 q4 pfunction [f,T]=TSPSA(d,t0,tf)
%TSP问题（货郎担问题，旅行商问题）的模拟退火算法通用malab源程序
2 V# ]0 W- k, W% s# _7 @% f目标最优值,T最优路线,d距离矩阵,t0初始温度,tf结束温度
% ?! O5 T& [/ v, I9 C* _; g[m,n]=size(d);
L=100*n;
t=t0;
- ]- w7 @( H. }) fpi0=1:n;
min_f=0;
for k=1:n-1
min_f=min_f+d(pi0(k),pi0(k+1));
end
1 J6 k/ ~' a: I' Qmin_f=min_f+d(pi0(n),pi0(1));
p_min=pi0;
while t>tf
1 `/ u+ ?1 C  A% V4 `) t  m. U3 vfor k=1:L;
, D9 V' g- K6 h) Q8 c7 b  ?kk=rand;
1 q$ K$ N& J4 i9 R6 u6 {9 D[d_f,pi_1]=exchange_2(pi0,d);
r_r=rand;
' b- j6 L5 s+ i6 W4 ]% R& p: \) bif d_f<0
, b/ t' G' B0 y6 u8 ?pi0=pi_1;
% c& E! \8 H  I6 R. g% t# jelseif exp(d_f/t)>r_r
; M: |" u8 _0 R1 {0 O( qpi0=pi_1;
else
% A1 Y+ s/ K8 {& k9 W1 _pi0=pi0;
end
end
! f! y2 {8 o  A4 Q3 \f_temp=0;
for k=1:n-1
! m; N  m/ B  Pf_temp=f_temp+d(pi0(k),pi0(k+1));
; o* n/ z1 f2 d/ Uend
f_temp=f_temp+d(pi0(n),pi0(1));
+ k1 C; y' _0 I+ gif min_f>f_temp
min_f=f_temp;
  P) H# A/ `2 k- I# a, ~7 x) ip_min=pi0;
end
4 |6 ^; W# H/ [t=0.87*t;
: K  E1 C4 `: Wend
f=min_f;
T=p_min;
%aiwa要调用的子程序,用于产生新解
0 [3 _1 N3 H. @( H' Efunction [d_f,pi_r]=exchange_2(pi0,d)
1 ^, ]5 t* c- u4 w5 ~1 W+ g[m,n]=size(d);
clear m;
u=rand;
u=u*(n-2);
5 U3 N$ _& P! M9 o8 p: ?u=round(u);
if u<2
u=2;
end
if u>n-2
u=n-2;
end
8 E# Y6 |; @( r/ _v=rand;
/ i: [) s1 A; P/ Q+ O; Pv=v*(n-u+1);
# K( W) \5 O, T# N# c$ kv=round(v);
* ?% `/ A) V- {& X& Yif v<1
5 E- Y8 x- @; E% t) W  V" ov=1;
end
v=u+v;
if v>n
v=n;
end
pi_1(u)=pi0(v);
: E& P' i9 u! t0 J* g# upi_1(v)=pi0(u);
1 ?% a- F3 `2 }if u>1
for k=1:u-1
+ Q; y: |' C6 U! p. ^# cpi_1(k)=pi0(k);
end
end
if v>(u+1)
for k=1:v-u-1
pi_1(u+k)=pi0(v-k);
0 N6 K1 z; Q# h  [' E# o0 v4 \end
end
" P0 Q+ P3 r/ ^, f! o: R8 Qif v<n
! }6 B* s# }4 Z. F& x* Gfor k=(v+1):n
( j. W% O8 `$ x# p+ P3 h' A6 Spi_1(k)=pi0(k);
8 U. z  O1 I0 Q. J) @9 E4 i# ^" nend
- B9 K5 ?: K. J7 \end
d_f=0;
3 V( N' x& i- T; Bif v<n
( {( q0 G, c5 r9 n. I7 |d_f=d(pi0(u-1),pi0(v))+d(pi0(u),pi0(v+1));
$ W5 m3 E8 {* D8 K; w9 E$ Ifor k=(u+1):n
) k* I+ L: A7 [+ {d_f=d_f+d(pi0(k),pi0(k-1));
8 j9 t8 O- z" pend
d_f=d_f-d(pi0(u-1),pi0(u))-d(pi0(v),pi0(v+1));
for k=(u+1):n
  R6 D) ~5 O. ~d_f=d_f-d(pi0(k-1),pi0(k));
( B0 P# y& n% u* l) e9 Tend
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));
" ~; u4 Z" W$ R0 S+ @6 H6 mend
for k=(u+1):n
d_f=d_f-d(pi0(k-1),pi0(k));
end
end
pi_r=pi_1;
( Q2 f  N) M# i1 M
得到：
  [f,T]=TSPSA(A,0,99)
2 [1 V& V  f& G* k! k& g
1 p( s/ C8 V( d, [f =
* L9 L3 C1 s1 _$ J
   Inf

) N# P% x  u- f" _+ [0 T
T =
+ P: `+ v$ e8 T: X  q/ I* c
. b% q6 @4 P9 L# F" |% o0 H  Columns 1 through 18

     1     2     3     4     5     6     7     8     9    10    11    12    13    14    15    16    17    18

  Columns 19 through 33

, N! A: X. C2 W) ^4 Y0 G) w- i    19    20    21    22    23    24    25    26    27    28    29    30    31    32    33

这个初始，结束温度是自己随便设定的？？
得到的这个F是无穷？？？难道？？？  T又是什么意思呢？？

冰E柠檬

在利用模拟退火算法进行优化之前，必须首先选取一个优化的起始点，优化起始点可以随机选取也可根据经验选取.
, B6 D) z! T$ V) \. [F是目标最优值,T为最优路线。
$ e* y+ \+ S) F6 L! S5 kF是无穷可能说明所求目标值没有路径到达没有最优解，或者是不是由于哪个步骤出现了什么问题（比如数据，或者邻接矩阵那里等等。。。）导致没有最优解。。。T表示的是最优解的起点到终点的最优路线。。。。。

magic2728

模拟退火算法是一种模拟工业退火时的一种智能算法，是利用这种自然现象带来的启发帮助算法设计。
应该来说，这个算法的各个参数的调整是需要一些经验的。初始温度和结束温度确实需要自己设定才是，而且设定好坏直接决定你的解的好坏。f是无穷表明当前的得到的最优解的路径长度是无穷，也就是说，你现在的路径里，存在两个不相通的目标点；T是路径，可以看出这个路径就是从第一个点顺次走下去，所以应该是你的温度值设定不当，导致退货过程不佳而造成的，应该调整温度来重新测试。