模拟退火算法（TSP问题）
生活随笔
收集整理的這篇文章主要介紹了
模拟退火算法（TSP问题）
小編覺(jué)得挺不錯(cuò)的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
模擬退火算法解決TSP問(wèn)題
算法思想
模擬退火算法（Simulate Anneal，SA）是一種通用概率演算法，用來(lái)在一個(gè)大的搜尋空間內(nèi)找尋命題的最優(yōu)解
模擬退火算法來(lái)源于固體退火原理，將固體加溫至充分高，再讓其徐徐冷卻，加溫時(shí)，固體內(nèi)部粒子隨溫升變?yōu)闊o(wú)序狀，內(nèi)能增大，而徐徐冷卻時(shí)粒子漸趨有序，在每個(gè)溫度都達(dá)到平衡態(tài)，最后在常溫時(shí)達(dá)到基態(tài)，內(nèi)能減為最小。根據(jù)Metropolis準(zhǔn)則，粒子在溫度T時(shí)趨于平衡的概率為e(-ΔE/(kT))，其中E為溫度T時(shí)的內(nèi)能，ΔE為其改變量，k為Boltzmann常數(shù)。用固體退火模擬組合優(yōu)化問(wèn)題，將內(nèi)能E模擬為目標(biāo)函數(shù)值f，溫度T演化成控制參數(shù)t，即得到解組合優(yōu)化問(wèn)題的模擬退火算法：由初始解i和控制參數(shù)初值t開(kāi)始，對(duì)當(dāng)前解重復(fù)“產(chǎn)生新解→計(jì)算目標(biāo)函數(shù)差→接受或舍棄”的迭代，并逐步衰減t值，算法終止時(shí)的當(dāng)前解即為所得近似最優(yōu)解，這是基于蒙特卡羅迭代求解法的一種啟發(fā)式隨機(jī)搜索過(guò)程。退火過(guò)程由冷卻進(jìn)度表(Cooling Schedule)控制，包括控制參數(shù)的初值t及其衰減因子Δt、每個(gè)t值時(shí)的迭代次數(shù)L和停止條件S。
設(shè)計(jì)思路
初始化溫度T，初始解狀態(tài)S，每個(gè)溫度t下的迭代次數(shù)L；
當(dāng)k = 1，2，……，L時(shí)，進(jìn)行3~6；
對(duì)當(dāng)前解進(jìn)行變換得到新解S’（例如對(duì)某些解中的元素進(jìn)行互換，置換）；
計(jì)算增量Δt′=C(S′)-C(S)，其中C(S)為評(píng)價(jià)函數(shù)；
若Δt′<0則接受S′作為新的當(dāng)前解，否則以概率exp(-Δt′/(KT))接受S′作為新的當(dāng)前解(k為玻爾茲曼常數(shù)，數(shù)值為：K=1.3806505(24) × 10^-23 J/K）；
如果滿足終止條件則輸出當(dāng)前解作為最優(yōu)解，結(jié)束程序；
減小T，轉(zhuǎn)到第2步，直到T小于初始設(shè)定的閾值。
程序代碼（TSP問(wèn)題）
第一種
#include <stdlib.h> #include <stdio.h> #include <iostream> #include<time.h> #include<math.h> #define N 109 #define big 9999999 using namespace std;double nowDist(double d[N][N],int temp[N],int m){//計(jì)算當(dāng)前temp路線組合距離double dist=0;for(int i=1;i<m;i++)dist+=d[temp[i]][temp[i+1]];dist+=d[temp[m]][temp[1]];return dist; }int main(){double d[N][N];//距離矩陣int temp[N];//臨時(shí)地點(diǎn)排列順序表double dist=0;//最小距離int n,m;//n為輸入行數(shù),m為地點(diǎn)個(gè)數(shù)int res[N];double mostdist=big;//數(shù)據(jù)輸入freopen("input.txt","r",stdin);cin>>m>>n;while(n-->0){int i,j;double k;cin>>i>>j>>k;d[i][j]=d[j][i]=k;}for(int i=1;i<=m;i++){d[i][i]=big;temp[i]=i;if(i!=m) dist+=d[i][i+1];}dist+=d[m][1];cout<<"初始距離: "<<dist<<endl;int count=40;while(count--){//取count次迭代的最小結(jié)果double T=1;//初始溫度double a=0.999;//退火率double low=1e-30;//最低溫度long ct=0;//迭代步數(shù)long ctj=0;//有效迭代步數(shù)srand(time(0));//隨機(jī)數(shù)種子long all=300000;//最大迭代次數(shù)while(T>low){int t1=rand()%m+1,t2=rand()%m+1;//隨機(jī)獲得兩個(gè)位置if(t1!=t2){swap(temp[t1],temp[t2]);//交換兩位置double ndist=nowDist(d,temp,m);//交換后的路線總距離double df=(ndist-dist);if(df<0){//當(dāng)前解更優(yōu)dist=ndist;ctj++;}else if(exp(-df/T)>(rand()%100)/100.0){//大于dist=ndist;ctj++;}else{swap(temp[t1],temp[t2]);}T*=a;//降溫}ct++;if(ct>all)break;}cout<<"當(dāng)前最短距離:"<<dist<<endl;if(dist<mostdist){mostdist=dist;for(int i=1;i<=m;i++)res[i]=temp[i];}}cout<<"最短距離:"<<mostdist<<endl<<"路線:";for(int i=1;i<=m;i++){cout<<res[i]<<" ";}return 0; }第二種
#include <stdlib.h> #include <stdio.h> #include <iostream> #include<time.h> #include<math.h>#include<fstream> #include<algorithm> #include<memory.h> using namespace std;const int num = 1000;//city number const int width = 100; const int height = 100;typedef struct node {int x;int y; }city; city citys[num];//citys double dic[num][num];//distance from two citys; bool visit[num];//visited int N;//real citys int seq[num];//最優(yōu)路徑序列 double answer;//最優(yōu)路徑長(zhǎng)度 const int tempterature = 1000;//初始溫度 const double u = 0.998;//成功降溫因子 const double v = 0.999;//失敗降溫因子 int k = 100;//對(duì)每個(gè)溫度迭代次數(shù) void init() {//set N&&x-yN = 51;citys[0].x = 37; citys[0].y = 52;citys[1].x = 49; citys[1].y = 49;citys[2].x = 52; citys[2].y = 64;citys[3].x = 20; citys[3].y = 26;citys[4].x = 40; citys[4].y = 30;citys[5].x = 21; citys[5].y = 47;citys[6].x = 17; citys[6].y = 63;citys[7].x = 31; citys[7].y = 62;citys[8].x = 52; citys[8].y = 33;citys[9].x = 51; citys[9].y = 21;citys[10].x = 42; citys[10].y = 41;citys[11].x = 31; citys[11].y = 32;citys[12].x = 5; citys[12].y = 25;citys[13].x = 12; citys[13].y = 42;citys[14].x = 36; citys[14].y = 16;citys[15].x = 52; citys[15].y = 41;citys[16].x = 27; citys[16].y = 23;citys[17].x = 17; citys[17].y = 33;citys[18].x = 13; citys[18].y = 13;citys[19].x = 57; citys[19].y = 58;citys[20].x = 62; citys[20].y = 42;citys[21].x = 42; citys[21].y = 57;citys[22].x = 16; citys[22].y = 57;citys[23].x = 8; citys[23].y = 52;citys[24].x = 7; citys[24].y = 38;citys[25].x = 27; citys[25].y = 68;citys[26].x = 30; citys[26].y = 48;citys[27].x = 43; citys[27].y = 67;citys[28].x = 58; citys[28].y = 48;citys[29].x = 58; citys[29].y = 27;citys[30].x = 37; citys[30].y = 69;citys[31].x = 38; citys[31].y = 46;citys[32].x = 46; citys[32].y = 10;citys[33].x = 61; citys[33].y = 33;citys[34].x = 62; citys[34].y = 63;citys[35].x = 63; citys[35].y = 69;citys[36].x = 32; citys[36].y = 22;citys[37].x = 45; citys[37].y = 35;citys[38].x = 59; citys[38].y = 15;citys[39].x = 5; citys[39].y = 6;citys[40].x = 10; citys[40].y = 17;citys[41].x = 21; citys[41].y = 10;citys[42].x = 5; citys[42].y = 64;citys[43].x = 30; citys[43].y = 15;citys[44].x = 39; citys[44].y = 10;citys[45].x = 32; citys[45].y = 39;citys[46].x = 25; citys[46].y = 32;citys[47].x = 25; citys[47].y = 55;citys[48].x = 48; citys[48].y = 28;citys[49].x = 56; citys[49].y = 37;citys[50].x = 30; citys[50].y = 40; } void set_dic() {//set distancefor (int i = 0; i<N; ++i) {for (int j = 0; j<N; ++j) {dic[i][j] = sqrt(pow(citys[i].x - citys[j].x, 2) + pow(citys[i].y - citys[j].y, 2));}} } double dic_two_point(city a, city b) {return sqrt(pow(a.x - b.x, 2) + pow(a.y - b.y, 2)); } double count_energy(int* conf) {double temp = 0;for (int i = 1; i<N; ++i) {temp += dic_two_point(citys[conf[i]], citys[conf[i - 1]]);}temp += dic_two_point(citys[conf[0]], citys[conf[N - 1]]);return temp; } bool metro(double f1, double f2, double t) {if (f2 < f1)return true;//else// return false;double p = exp(-(f2 - f1) / t);int bignum = 1e9;if (rand() % bignum<p*bignum)return true;return false; } void generate(int* s) {//隨機(jī)產(chǎn)生一組新解bool v[num];memset(v, false, sizeof(v));for (int i = 0; i<N; ++i) {s[i] = rand() % N;while (v[s[i]]) {s[i] = rand() % N;}v[s[i]] = true;} } void generate1(int* s) {//隨機(jī)交換序列中的一組城市順序int ti = rand() % N;int tj = ti;while (ti == tj)tj = rand() % N;for (int i = 0; i<N; ++i)s[i] = seq[i];swap(s[ti], s[tj]); } void generate2(int* s) {//隨機(jī)交換序列中的兩組城市順序int ti = rand() % N;int tj = ti;int tk = ti;while (ti == tj)tj = rand() % N;while (ti == tj || tj == tk || ti == tk)tk = rand() % N;for (int i = 0; i<N; ++i)s[i] = seq[i];swap(s[ti], s[tj]);swap(s[tk], s[tj]); } void generate3(int* s) {//隨機(jī)選序列中的三個(gè)城市互相交換順序int ti = rand() % N;int tj = ti;int tm = rand() % N;int tn = ti;while (ti == tj)tj = rand() % N;while (tm == tn)tn = rand() % N;for (int i = 0; i<N; ++i)s[i] = seq[i];swap(s[ti], s[tj]);swap(s[tm], s[tn]); } void generate0(int* s) {//以上三種交換方式等概率選擇int temp = rand() % 3;if (temp == 0)generate1(s);else if (temp == 1)generate2(s);else if (temp == 2)generate3(s); } void moni() {double t = tempterature;int seq_t[num];for (int i = 0; i<N; ++i) {//初始化當(dāng)前序列seq[i] = seq_t[i] = i;}double new_energy = 1, old_energy = 0;while (t>1e-9&&fabs(new_energy - old_energy)>1e-9) {//溫度作為控制變量int t_k = k;int seq_tt[num];while (t_k--&&fabs(new_energy - old_energy)>1e-9) {//迭代次數(shù)作為控制變量generate1(seq_tt);new_energy = count_energy(seq_tt);//newold_energy = count_energy(seq_t);//oldif (metro(old_energy, new_energy, t))for (int i = 0; i < N; ++i)seq_t[i] = seq_tt[i];}new_energy = count_energy(seq_t);//newold_energy = answer;//oldif (metro(old_energy, new_energy, t)) {for (int i = 0; i < N; ++i)seq[i] = seq_t[i];answer = count_energy(seq);t *= u;//接受新?tīng)顟B(tài)降溫因子0.98}elset *= v;//不接受新?tīng)顟B(tài)降溫因子0.99}answer = count_energy(seq); } void output() {cout << "the best road is : \n";for (int i = 0; i < N; ++i) {cout << seq[i];if (i == N - 1)cout << endl;elsecout << " -> ";}cout << "the length of the road is " << answer << endl; } void test() {ifstream ifile("data.txt");if (!ifile) {cout << "open field\n";return;}while (!ifile.eof()) {int te = 0;ifile >> te;ifile >> citys[te - 1].x >> citys[te - 1].y;} } int main() {srand(time(0));int t;while (cin >> t) {//僅作為重啟算法開(kāi)關(guān)使用，無(wú)意義init();//使用程序內(nèi)置數(shù)據(jù)使用init()函數(shù)，//test();//使用文件讀取數(shù)據(jù)使用test()函數(shù)，set_dic();//計(jì)算每個(gè)城市之間的距離moni();//退火output();//輸出}return 0; }測(cè)試?yán)?/strong>
第一種，通過(guò)以下代碼生成數(shù)據(jù)集文件
第二種
1 37 52 2 49 49 3 52 64 4 20 26 5 40 30 6 21 47 7 17 63 8 31 62 9 52 33 10 51 21 11 42 41 12 31 32 13 5 25 14 12 42 15 36 16 16 52 41 17 27 23 18 17 33 19 13 13 20 57 58 21 62 42 22 42 57 23 16 57 24 8 52 25 7 38 26 27 68 27 30 48 28 43 67 29 58 48 30 58 27 31 37 69 32 38 46 33 46 10 34 61 33 35 62 63 36 63 69 37 32 22 38 45 35 39 59 15 40 5 6 41 10 17 42 21 10 43 5 64 44 30 15 45 39 10 46 32 39 47 25 32 48 25 55 49 48 28 50 56 37 51 30 40運(yùn)行結(jié)果
第一種
第二種
the best road is : 13 -> 24 -> 22 -> 23 -> 50 -> 4 -> 48 -> 9 -> 32 -> 44 -> 38 -> 29 -> 33 -> 20 -> 28 -> 34 -> 35 -> 2 -> 21 -> 27 -> 30 -> 19 -> 1 -> 49 -> 37 -> 8 -> 15 -> 31 -> 10 -> 0 -> 7 -> 25 -> 42 -> 6 -> 47 -> 26 -> 5 -> 45 -> 36 -> 14 -> 43 -> 41 -> 18 -> 39 -> 40 -> 12 -> 16 -> 3 -> 46 -> 11 -> 17分析
模擬退火算法具有以下優(yōu)缺點(diǎn)；
1.迭代搜索效率高，并且可以并行化；
2.算法中有一定概率接受比當(dāng)前解較差的解，因此一定程度上可以跳出局部最優(yōu)；
3.算法求得的解與初始解狀態(tài)S無(wú)關(guān)，因此有一定的魯棒性；
4.具有漸近收斂性，已在理論上被證明是一種以概率l收斂于全局最優(yōu)解的全局優(yōu)化算法。
參考
https://blog.csdn.net/daaikuaichuan/article/details/81381875
https://blog.csdn.net/wordsin/article/details/79915328
https://blog.csdn.net/chenlnehc/article/details/51933802
總結(jié)
以上是生活随笔為你收集整理的模拟退火算法（TSP问题）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。
- 上一篇： 6-搭建一个私有registry
- 下一篇： Angular6_服务端渲染SSR