请问谁能提供求解TSP的程序

pleaseopen

问题：
  1、能量函数是否唯一？
  2、求解TSP目前较好的方法有哪些？
  3、谁有matlab解TSP的程序呢？
                                                                 
                             谢谢。

gghhjj

蚁群算法 遗传算法 模拟退火算法等等很多
论坛也有部分相关程序，自己搜索一下吧

pleaseopen

神经网络的方法是否可以呢？

gghhjj

用神经网络方法解决tsp问题

1. %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];
   
2. %t0=10;
   
3. % tf=0.7;
   
4. % [f,T]=trp(d,t0,tf)
   
5. function[f,T]=trp(d,t0,tf)
   
6. [m,n]=size(d);
   
7. L=100*n;
   
8. t=t0;
   
9. pi0=1:n;
   
10. min_f=0;
    
11. for k=1:(n-1)
    
12.     min_f=min_f+d(pi0(k),pi0(k+1));
    
13. end
    
14. min_f=min_f+d(pi0(n),pi0(1));
    
15. p_min=pi0;
    
16. while t>tf
    
17.     for k=1:L
    
18.         kk=rand;
    
19.         [d_f,pi_1]=exchange_2(pi0,d);
    
20.         r_r=rand;
    
21.         if d_f<0
    
22.             pi0=pi_1;
    
23.         elseif exp(d_f/t)>r_r
    
24.             pi0=pi_1;
    
25.         else pi0=pi0;
    
26.         end
    
27.     end
    
28.     f_temp=0;
    
29.     for k=1:(n-1)
    
30.         f_temp=f_temp+d(pi0(k),pi0(k+1));
    
31.     end
    
32.     f_temp=f_temp+d(pi0(n),pi0(1));
    
33.     if min_f>f_temp
    
34.         min_f=f_temp;
    
35.         p_min=pi0;
    
36.     end
    
37.     t=0.87*t;
    
38. end
    
39. f=min_f;
    
40. T=p_min;
    
41. function [d_f,pi_r]=exchange_2(pi0,d)
    
42. [m,n]=size(d);
    
43. clear m;
    
44. u=rand;
    
45. u=u*(n-2);
    
46. u=round(u);
    
47. if u<2
    
48.     u=2;
    
49. end
    
50. if u>n-2
    
51.     u=n-2;
    
52. end
    
53. v=rand;
    
54. v=v*(n-u+1);
    
55. v=round(v);
    
56. if v<1
    
57.     v=1;
    
58. end
    
59. v=u+v;
    
60. if v>n
    
61.     v=n;
    
62. end
    
63. pi_1(u)=pi0(v);
    
64. pi_1(v)=pi0(u);
    
65. if u>1
    
66.     for k=1:(u-1)
    
67.         pi_1(k)=pi0(k);
    
68.     end
    
69. end
    
70. if v>(u+1)
    
71.     for k=1:(v-u-1)
    
72.         pi_1(u+k)=pi0(v-k);
    
73.     end
    
74. end
    
75. if v<n
    
76.     for k=(v+1):n
    
77.         pi_1(k)=pi0(k);
    
78.     end
    
79. end
    
80. d_f=0;
    
81. if v<n
    
82.     d_f=d(pi0(u-1),pi0(v))+d(pi0(u),pi0(v+1));
    
83.     for k=(u+1):n
    
84.         d_f=d_f+d(pi0(k),pi0(k-1));
    
85.     end
    
86.     d_f=d_f-d(pi0(u-1),pi0(u))-d(pi0(v),pi0(v+1));
    
87.     for k=(u+1):n
    
88.         d_f=d_f-d(pi0(k-1),pi0(k-1));
    
89.     end
    
90. else
    
91.     d_f=d(pi0(u-1),pi0(v))+d(pi0(u),pi0(1))-d(pi0(u-1),pi0(u))-d(pi0(v),pi0(1));
    
92.     for k=(u+1):n
    
93.         d_f=d_f+d(pi0(k),pi0(k-1));
    
94.     end
    
95.     for k=(u+1):n
    
96.         d_f=d_f-d(pi0(k-1),pi0(k));
    
97.     end
    
98. end
    
99. pi_r=pi_1;

复制代码

gghhjj

hopfield解决TSP问题
http://forum.vibunion.com/forum/thread-30110-1-1.html

pleaseopen

谢谢，正在看
:@D