实验二蚁群算法(共6页)

1、人工智能课程实验报告之二基于ACO/PSO的TSP问题求解班级:学号：姓名：成绩评定：评阅老师：日 期：实验报告正文一、实验目的 初步掌握ACO对旅行商问题（TSP）的设计与实现,加深对ACO算法的的理解。二、实验内容 根据ACO算法原理，利用带权完全图G=（V，E）其中V = v1,v2,vn 是带有n = |V|个城市的集合，E是完全连接这些点的边的集合，来模拟n个城市的tsp问题，设计并编程实现一个城市的TSP最佳路径求解系统系统。三、实验所用智能算法基本原理与流程ACO算法可以直接应用到TSP中， ACO的构建图为G = (C, L)， 为一个完全图，构建图与描述TSP实例的图形一致，

2、也就是说C = V且L = E；问题的状态集合与所有可能的部分路径集合相对应；约束条件为蚂蚁建立的路径必须是与城市标号的排列相对应的可行路径。AS算法有两个最主要的步骤及蚂蚁构建问题的解和信息素的更新。在AS中m只蚂蚁并行地构建TSP的路径。算法开始时，将m只蚂蚁放在随机选择出来的城市中。在路径构建过程中，每一只蚂蚁都按随机比例规则的概率行为选择规则来随机的选择下一个城市。当所有的蚂蚁都构建好路径之后，开始进行信息素的更新操作。先进行信息素的蒸发，所有边上的信息素按蒸发率 减少四、系统构成与程序设计ACO算法最重要的两个步骤是解的构建和信息素更新，除此之外，在算法中还要加入一些更新统计信息和判

3、断终止条件的方法。另外在解的构建过程中也可以加入局部搜索以提高算法性能，但本次实验中没有加入局部搜索。算法总的框架如下：ACOFORTSP( )1InitializeData( )/ 初始化数据2while (terminate() = false) do/ 判断是否终止3ConstructSolutions( )/ 解的构建4LocalSearch( )/ 局部搜索5UpdatePheromone( )/ 更新信息素6end-while解的构建算法如下：ConstructSolutions( )1for k = 1 to m/ 1 5 行初始化蚂蚁的记忆，把所有蚂蚁访问过的城市设2for I

4、 = 1 to n/ 为空3antk.visitedi false4end-for5end-for6step1/ 6 11行，将所有蚂蚁随机地放到n个城市中7for k = 1 to m8rrandom1,n9antk.tourstep r10antk.visitedr true11end-for12while (step < n)/ 12 17行，并行的构建每只蚂蚁的解13stepstep + 114for k = 1 to m do15ChooseNextCity(k, step)16end-for17end-while18for k = 1 to m/ 18 21 行，把每只蚂蚁最

5、后一步设为起始城市建立回路19antk.tourn + 1 antk.tour120antk.tourlengthComputeTourLength(k)21end-for信息素更新方法如下：UpdatePheromone( )1Evaporate( )/ 信息素蒸发2for k = 1 to m/ 2 4 行， 释放信息素3DepositPheromone( k )4end-for5UpdateGlobalData( )/ 计算一些全局信息其程序流程图为： 五、核心代码TSP中的核心代码为：void Perm(Type list, int k, int m)/递归函数，生成排列，寻找最佳路径

6、 int c=1;int a=m;while(a>1)/定义Mat数组的行数c=c*a*(a-1);a-=2;int d=m+1;/定义Mat数组的列数if ( k = m && list0=0)for(int i = 0; i <= m; i+)Matei=listi;/符合条件时生成排列，存储在数组Mat中e+; if(Matc-1m!=-1 && e=c )/当找出所有可能路径时，计算所有路径长度，并找到最短路径for(int g=0;g<c;g+) for(int h=0;h<d;h+) intx=Matgh;/获取路径坐标 in

7、t y=Matgh+1; Countg=Countg+G.arcsxy;/计算路径长度cout<<"所有路径及其路径长度："<<endl;for(g=0 ;g<c;g+)/输出所有路径及其路径长度for(int h=0;h<=d;h+)int k=Matgh;cout<<G.vexsk<<"->"cout<<" "<<Countg<<endl;for(g=0;g<c;g+)Helpg=Countg;/赋值for(int i=0;i<c;i+)/对辅助数组冒泡排序，找最小值for(int j=1;j<c;j+)if (Helpi>Helpj)int temp=Helpi;Helpi=Helpj;Helpj=temp;int Min_way=Help0;/Min_way为最短路径for (i=0;i<c;i+)if (Counti=Min_way)cout<<"旅行商最佳路径为："for(int j=0;j<=d;j+)int k=Matij;cout<<G.vexsk<<&qu