模拟退火算法 解决旅行商(TSP)问题

1. #include<iostream>
2. #include<string>
3. #include<math.h>
4. #include<cstdlib>
5. using namespace std;
7. // ******************************************************************************************
8. /*
9.         模拟退火算法解决10城市的TSP问题
10.     (1) 初始化：初始温度T(充分大)，初始解状态S(是算法迭代的起点)， 每个T值的迭代次数L 
11.     (2) 对k=1，……，L做第(3)至第6步： 
12. 　　(3) 产生新解S′ 
13. 　　(4) 计算增量Δt′=C(S′)-C(S)，其中C(S)为评价函数 
14. 　　(5) 若Δt′<0则接受S′作为新的当前解，否则以概率exp(-Δt′/T)接受S′作为新的当前解. 
15. 　　(6) 如果满足终止条件则输出当前解作为最优解，结束程序。 
16. 　     　终止条件通常取为连续若干个新解都没有被接受时终止算法。 
17. 　　(7) T逐渐减少，且T->0，然后转第2步。
18.                                                     baixzh  SYSU 2008-11-19 
19. */  
20. // ******************************************************************************************
21.     int i;int j;
22.     struct city//城市节点 坐标表示
23.     {
24.         double x;
25.         double y;
26.         city();
27.          
29.     };
31.     city::city ()
32.     {
33.         x=0.0;
34.         y=0.0;
35.     }
36.      
37.     city c_t[10];//10个城市的tsp问题
38.     double pos[20]=//初始城市坐标数据
39.     {
40.         0.400,  0.4439,//代表x,y
41.         0.2439, 0.1463,
42.         0.1707, 0.2293,
43.         0.2293, 0.7610,
44.         0.5171, 0.9414,
45.         0.8732, 0.6536,
46.         0.6878, 0.5219,
47.         0.8488, 0.3609,
48.         0.6683, 0.2536,
49.         0.6195, 0.2634
51.     };
53.     double dis[10][10];//城市距离矩阵
54.      
55.     int state[10]={7,8,9,1,5,0,6,2,4,3};//给出初始路径
56.     int s1[10],s2[10];//s1原始路径 s2为s1生成的临近路径
57.     double result1,result2;//原路径长度 临近路径长度
60.     void xch(int state[10])//生成临近路径并且计算长度,保存原路径
61.     {   
62.         result1=result2=0.0;
63.         for(i=0;i<10;i++)
64.         {
65.             s1[i]=s2[i]=state[i];
66.         }
67.         for(i=0;i<9;i++)
68.         {
69.             result1=result1+dis[s1[i]][s1[i+1]];
70.         }
71.         result1=result1+dis[s1[0]][s1[9]];
72.         //s1以前的s2改变后的
73.         int k1=rand()%10;int k2=rand()%10;
74.         if(k1==k2)
75.         {
76.             k2=(k1+2)%10;
78.         }
79.         swap(s2[k1],s2[k2]);
80.         for(i=0;i<9;i++)
81.         {
82.             result2=result2+dis[s2[i]][s2[i+1]];
83.         
84.         }
85.         result2=result2+dis[s2[0]][s2[9]];
86.         //最后要返回出发城市
87.     }
88.      
89.     int main()
90.     {   
91.         for(i=0;i<10;i++)
92.         {
93.             c_t[i].x =pos[2*i];
94.             c_t[i].y =pos[2*i+1];
95.         }
96.         //城市节点初始化
97.         memset(dis,0.0,sizeof(dis));
98.         for(i=0;i<10;i++)
99.             for(j=0;j<10;j++)
100.             {
101.                          dis[i][j]=sqrt(pow((c_t[i].x -c_t[j].x),2) +pow((c_t[i].y-c_t[j].y),2)) ; 
102.             }
103.         //距离矩阵初始化
104.         double tmp=1000.0;
105.         int loop_1=1;//外循环计数
106.         int loop_2=0;//内循环计数
107.         while(1)
108.         {       
109.             loop_2=0;
110.             while(1)
111.             {   
112.                 if(loop_2>tmp*3)//内循环结束条件
113.                 {
114.                     break;
115.                 }
116.                 xch(state);
117.                 loop_2++;
118.                 if(result2<result1)
119.                 {
120.                     for(j=0;j<10;j++)
121.                     {
122.                         state[j]=s2[j];
123.                     }
124.                 //如果生成路径长度小于原始路径则替换             
125.                 }
126.                 else
127.                 {
128.                     double d=(result2-result1);
129.                     double e=exp(-d/(0.0001*tmp));
130.                     double r=rand()%100;
131.                     double k=r/100.0;
132.                     if(e>k)
133.                     {
134.                         for(j=0;j<10;j++)
135.                     {
136.                         state[j]=s2[j];
137.                     }
138.                     //如果生成路径长度大于初始路径长度//但是概率e>rand(0,1)也接受新路径
139.                     }
140.                     else
141.                     {
142.                     for(j=0;j<10;j++)
143.                     {
144.                         state[j]=s1[j];
145.                     }
146.                     }//否则不做替换           
147.                 }
148.             }
149.             //输出每次内部循环结束后的路径
150.             xch(state);
151.             cout<<"第"<<loop_1<<"次降温 "<<"温度为"<<tmp<<"     结果为    :"<<endl;
152.             cout<<result1<<"    ";
153.             for(j=0;j<9;j++)
154.             {
155.                 cout<<char('A'+state[j])<<"--";
156.             }
157.             cout<<char('A'+state[9])<<endl;
159.             tmp=tmp*0.95;//温度按照比例下降
160.             loop_1++;
161.             if(tmp<1)//外循环结束条件
162.             {
163.                 break;
164.             }
165.             
166.         } 
168.      
169.         return 0;
170.     }
171.