蚁群算法 C语言程序(已运行).docVIP

下载本文档

56
0
约3.8千字
约 5页
2017-12-25 发布于江西
举报
版权申诉

蚁群算法 C语言程序(已运行).doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

蚁群算法 C语言程序(已运行)

//Basic Ant Colony Algorithm for TSP #include iostream.h #include fstream.h #include math.h #include time.h #include conio.h #include stdlib.h #include iomanip.h #define N 31 //city size #define M 31 //ant number double inittao=1; double tao[N][N]; double detatao[N][N]; double distance[N][N]; double yita[N][N]; int tabu[M][N]; int route[M][N]; double solution[M]; int BestRoute[N]; double BestSolution=10000000000; double alfa,beta,rou,Q; int NcMax; void initparameter(void); // initialize the parameters of basic ACA double EvalueSolution(int *a); // evaluate the solution of TSP, and calculate the length of path void InCityXY( double x[], double y[], char *infile ); // input the nodes coordinates of TSP void initparameter(void) { alfa=1; beta=5; rou=0.9; Q=100; NcMax=200; } void main(void) { int NC=0; initparameter(); double x[N]; double y[N]; InCityXY( x, y, city31.tsp ); for(int i=0;iN;i++) for(int j=i+1;jN;j++) { distance[j][i]=sqrt((x[i]-x[j])*(x[i]-x[j])+(y[i]-y[j])*(y[i]-y[j])); distance[i][j]=distance[j][i]; } // calculate the heuristic parameters for(i=0;iN;i++) for(int j=0;jN;j++) { tao[i][j]=inittao; if(j!=i) yita[i][j]=100/distance[i][j]; } for(int k=0;kM;k++) for(i=0;iN;i++) route[k][i]=-1; srand(time(NULL)); for(k=0;kM;k++) { route[k][0]=k%N; tabu[k][route[k][0]]=1; } //each ant try to find the optiamal path do { int s=1; double partsum; double pper; double drand; //ant choose one whole path while(sN) { for(k=0;kM;k++) { int jrand=rand()%3000; drand=jrand/3001.; partsum=0; pper=0; for(int j=0;jN;j++) { if(tabu[k][j]==0) partsum+=pow(tao[route[k][s-1]][j],alfa)*pow(yita[route[k][s-1]][j],beta); } for(j=0;jN;j++) { if(tabu[k][j]==0) pper+=pow(tao[route[k][s-1]][j],alfa)*pow(yita[route[k][s-1]][j],beta)/partsum; if(pperdrand) break; } tabu[k][j]=1; route[k][s]=j; } s++; } // the pheromone is updated for(i=0;iN;i++) for(int j=0;jN;j++) detatao[i][j]=0; for(k=0;kM;k+