基于模拟退火算法的TSP算法.docx

《基于模拟退火算法的TSP算法.docx》由会员分享，可在线阅读，更多相关《基于模拟退火算法的TSP算法.docx（10页珍藏版）》请在第一文库网上搜索。

1、专业综合设计报告课程名称：电子专业综合设计设计名称：基于模拟退火算法的TSP算法姓名：学号：班级：电子0903指导教师：朱正为起止日期：专业综合设计任务书学生班级：电子0903学生姓名：学号：20235830设计名称：基于模拟退火算法的TSP算法起止日期：指导教师设计要求：旅行商问题，即TSP问题fTrave11ingSa1esmanProb1em）又译为旅行推销员问题、货郎担问题，是数学领域中著名问题之一。假设有一个旅行商人要拜访n个城市，他必须选择所要走的路径，路径的限制是每个城市只能拜访一次，而且最后要回到原来出发的城市。路径的选择目标是要求得的路径路程为所有路径之中的最小值。此设计是用

2、模拟退火算法来实现TSP问题的寻求最优解。专业综合设计学生日志时间设计内容2023.11.9初步了解模拟退火算法的TSP算法2023.11.12设计算法流程、确定解题思路2023.11.20讨论算法流程及解题思路的可行性，为仿真做准备2023.12.2运用MAT1AB软件进行实验仿真，分析仿真结果整理实验报告辩论专业综合设计考勤表周星期一星期二星期三星期四星期五专业综合设计评语表指导教师评语:成绩：指导教师：-设计目的和意义4二设计原理错误!未定义书签。2.1 模拟退火算法的根本原理52.2 TSP问题介绍6三详细设计步骤错误!未定义书签。3.1 .算法流程83.2 模拟退火算法实现步骤错误!

3、未定义书签。四设计结果及分析94.1MAT1AB程序实现及主函数9计算距离矩阵94.1.2 初始解104.1.3 生成新解104.1.4 Metropo1is准那么函数104.1.5 画路线轨迹图114.1.6 输出路径函数124.1.7 可行解路线长度函数124.1.8 模拟退火算法的主函数134.2.仿真结果15五体会9六参考文献10基于模拟退火算法的TSP算法一、设计目的和意义旅行商问题是组合优化领域里的一个典型的、易于描述却难以处理的NP难题，其可能的路径数目与城市数目是呈指数型增长的，求解非常困难。首先介绍了旅行商问题,给出了其数学描述以及实际应用，进而给出解决TSP的一种比拟精确的

4、算法一一模拟退火算法。然后阐述了模拟退火算法的根本原理，重点说明了其根本思想及关键技术。最后运用MAT1AB语言实现了该算法，并将其运用到解决旅行商问题的优化之中。数值仿真的结果说明了该方法能够对数据进行全局寻优，有效克服了基于导数的优化算法容易陷入局部最优的问题。了解模拟退火算法的TSP算法的根本思路及原理，并应用MAT1AB实现仿真，熟练掌握MAT1AB的操作方式及应用，能正确书写报告。二、设计原理2.1 模拟退火算法的根本原理模拟退火算法足20世纪80年代初提出的一种基于蒙特卡罗(MenteCarh)迭代求解策略的启发式随机优化算法。它通过MetrOPO1iS接受准那么概率接受劣化解并以

5、此跳出局部最优，通过温度更新函数的退温过程进行趋化式搜索并最终进入全局最优解集。其出发点是基于物理中固体物质的退火过程与一搬的组合优化问题之间的相似性。模拟退火法是一种通用的优化算法，其物理退火过程由以下三局部组成。(1)加温过程。其目的是增强粒子的热运动，使其偏离平衡位置。当温度足够高时，固体将熔为液体，从而消除系统原先存在的非均匀状态。(2)等温过程。对于与周围环境交换热量而温度不变的密封系统，系统状态的自发变化总是朝自由能减少的方向进行的，当自由能到达最小时，系统到达平衡状态。(3)冷却过程。使粒子热运动减弱，系统能量下降，得到晶体结构。其中，加热过程对应算法的设定初温，等温过程对应算法

6、的Metropo1is抽样过程，冷却过程对应控制参数的下降。这里能量的变化就是目标函数，要得到的最优解就是能量最低态。MetroPOIiS准那么是SA算法收敛于全局最优解的关键所在，Metropo1is准那么以一定的概率接受恶化解，这样就使算法跳离局部最优的陷阱。模拟退火算法为求解传统方法难处理的TSP问题提供了一个有效的途径和通用框架，并逐渐开展成一种迭代自适应启发式概率性搜索算法。模拟退火算法可以用以求解不同的非线性问题，对不可微甚至不连续的函数优化，能以较大的概率求的全局有化解,该算法还具有较强的鲁棒性、全局收敛性、隐含并行性及广泛的适应性，并且能处理不同类型的优化设计变量离散的、连续的

7、和混合型的)，不需要任何的辅助信息，对目标函数和约束函数没有任何要求。利用Metropo1is算法并适当的控制温度下降过程，在优化问题中具有很强的竞争力，此设计即为基于模拟退火算法的TSP算法。SA算法实现过程如下以最小化问题为例)：初始化：取初始温度TO足够大，令T=To任取初始解S1确定每个T时的迭代次数，即Metropo1is链长1。,1重复步骤3)。对当前温度T和k=1,2,3）对当前S1随机扰动产生一个新解S2。计算S2的增量df=f-f）其中f为S1的代价函数。15）假设dfrand也接受S2作为新的当前解，S1=S2;否那么保存当前解S1。16）如果满足终止条件Stop,那么输出

8、当前解S1为最优解，结束程序。终止条件Stop通常为：在连续假设干个Metropo1is链中新解s2都没有被接受时终止算法，或是设定结束温度。否那么按衰减函数衰减T后返回步骤2）。以上步骤成为Metropo1is过程。逐渐降低控制温度，重复Metropo1is过程，直至满足结束准那么Stop,求出最优解。2.2 TSP问题介绍旅行商问题Trave1ingSa1esmanProb1em,简称TSP）又名货郎担问题，是威廉哈密尔顿爵士和英国数学家克克曼（TRKirkman）于19世纪初提出的一个数学问题，也是著名的组合优化问题。问题是这样描述的：一名商人要到假设干城市去推销商品，城市个数和各城市间

9、的路程或旅费），要求找到一条从城市1出发，经过所有城市且每个城市只能访问一次，最后回到城市1的路线，使总的路程或旅费）最小。TSP刚提出时,不少人认为这个问题很简单。后来人们才逐步意识到这个问题只是表述简单，易于为人们所理解，而其计算复杂性却是问题的输入规模的指数函数，属于相当难解的问题。这个问题数学描述为：假设有n个城市，并分别编号，给定一个完全无向图G=V,E）,V=12,n,n1o其每一边（i,j）E有一非负整数消耗C,j（即上的权记为C中i,jV）oG的一条巡回路线是经过V中的每个顶点恰好一次的回路。一条巡回路线的消耗是这条路线上所有边的权值之和。TSP问题就是要找出G的最小消耗回路。

10、人们在考虑解决这个问题时,一般首先想到的最原始的一种方法就是:列出每一条可供选择的路线（即对给定的城市进行排列组合），计算出每条路线的总里程，最后从中选出一条最短的路线。假设现在给定的4个城市分别为A、B、C和D,各城市之间的消耗为己知数，如图1所示。我们可以通过一个组合的状态空间图来表示所有的组合，如图图1顶点带权图图2TSP问题的解空间树从图中不难看出，可供选择的路线共有6条，从中很快可以选出一条总消耗最短的路线:顶点序列为fA,C,B,D,A1由此推算，假设设城市数目为n时，那么组合路径数那么为（n-1）!o很显然，当城市数目不多时要找到最短距离的路线并不难，但随着城市数目的不断增大，组

11、合路线数将呈指数级数规律急剧增长，以至到达无法计算的地步，这就是所谓的“组合爆炸问题。假设现在城市的数目增为20个，组合路径数那么为(2-i)!.2i6o17,如此庞大的组合数目，假设计算机以每秒检索Iooo万条路线的速度计算，也需要花上386年的时间。三、详细设计步骤3.1 算法流程模拟退火算法求解流程框图如图1所示。图3模拟退火算法求解流程框图3.2 模拟退火算法实现步骤如下：11)控制参数的设置需要设置的主要控制参数有降温速率q、初始温度To、结束温度Tend以及链长1o初始解对于n个城市TSP问题，得到的解就是对1n的一个排序，其中每个数字为对应城市的编号，如对10个城市的TSP问题1

12、,2,3,4,5,6,7,8,9,10,那么11024568793就是一个合法的解，采用产生随机排列的方法产生一个初始解So13)解变换生成新解通过对当前解S1进行变换，产生新的路径数组即新解，这里采用的变换是产生随机数的方法来产生将要交换的两个城市，用二邻域变换法产生新的路径，即新的可行解S2。例如n=10时，产生两个1,10范围内的随机整数门和迨，确定两个位置，将其对换位置，如=4,2=795163871042得到的新解为95173861042Metropo1is准那么假设路径长度函数为/),新解的路径为/62),路径差为d/=/fS2S。，那么Metropo1is准那么为Idf0P-df

13、1 exp(4)W0如果dfO,那么以概率1接受新路线，否那么以概率exp-dfT)接受新路线。降温利用降温速率q进行降温，即TFT,假设T小于结束温度，那么停止迭代输出当前状态，否那么继续迭代。四、设计结果及分析4.1MAT1AB程序实现及主函数4.1.1计算距离矩阵利用给出的N个城市的坐标，算出N个城市的两两之间的距离，得到距离矩阵Nxn）。计算函数为DiStance,得到初始群种。程序如下functionD=Distanse(a)4.%计算两两城市之间的距离%输入a各城市的位置坐标%输出D两两城市之间的距离4二row=size(a,1);S1=randperm(N)；%随机产生一个初始路

14、线4.-forj=i+1:row解TX打仁”a(tn-a凡n)A2+(a(i.2)-a伍22人2)八06:functionS2=NewAnswer(S1)%输入%S1:当前解4.%输出%S2：新解N=1ength(S1);fnncf；CnqR1=N/Tafrcnci;/q1qz7TrR2=Path1ength(D,S2);%计算路线长度dC=R2-R1;%计算能力之差ifdC=rand%以exp(-dC/T)概率接受新路线S=S2;4R=R2;e1se%不接受新路线S=S1;R=R1;4endtgure;ho1donP1Ot(X(:,1),X(:,2),o;Cok)r,0.5,0.5,0.5)P1ot(X(ChrOm(1,1),1),X(Chrom(1,1),2),Tv?MarkerSiZe,20)fori=1:size(X,1)text(X(i,1)+0.05,X(i,2)+0.05,num2str(i),co1or,1,0,0);functionIen=Path1ength(D9Chrom)4