基于遗传算法的旅行商求解方案.doc

基于遗传算法的旅行商问题求解方案 摘要 采用MATLAB，对TSP问题进行基于遗传算法的求解。TSP 问题是典型的NP完全问题，通过MATLAB进行遗传算法编程，从而有效提出一个较好的TSP解，实现对问题的解答。进而讨论遗传算法的特点，以及对本问题的可行性。 关键词： TSP问题 遗传算法 一．问题重述 假设有一个旅行商人要拜访n个城市，他必须选择所要走的路径，路径的限制是每个城市只能拜访一次，而且最后要回到原来出发的城市。路径的选择目标是要求得的路径路程为所有路径之中的最小值。TSP问题是一个组合优化问题。该问题可以被证明具有NPC计算复杂性。因此，任何能使该问题的求解得以简化的方法，都将受到高度的评价和关注。 遗传算法（Genetic Algorithm）是模拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型，是一种通过模拟自然进化过程搜索最优解的方法，它最初由美国Michigan大学J.Holland教授于1975年首先提出来的，并出版了颇有影响的专著《Adaptation in Natural and Artificial Systems》，GA这个名称才逐渐为人所知，J.Holland教授所提出的GA通常为简单遗传算法（SGA）。 election）、交叉（Crossover）和变异（Mutation）是遗传算法的3个主要操作算子，它们构成了所谓的遗传操作（genetic operation），使遗传算法具有了其它传统方法所没有的特性。遗传算子包含如下6个基本因素： 参数编码：由于遗传算法不能直接处理解空间的解数据，因此必须通过编码将它们表示成遗传空间的基因型串结构数据。 生成初始群体：由于遗传算法的群体型操作需要，所以必须为遗传操作准备一个由若干初始解组成的初始群体。初始群体的每个个体都是通过随机方法产生。 适应度评估检测：遗传算法在搜索进化过程中一般不需要其他外部信息，仅用适应度（fitness）值来评估个体或解的优劣，并作为以后遗传操作的依据。 选择（selection）：选择或复制操作是为了从当前群体中选出优良的个体， 使它们有机会作为父代为下一代繁殖子孙。个体适应度越高，其被选择的机会就越多。此处采用与适用度成比例的概率方法进行选择。具体地说，就是首先计算群体中所有个体适应度的总和（），再计算每个个体的适应度所占的比例（），并以此作为相应的选择概率。 交叉操作：交叉操作是遗传算法中最主要的遗传操作。简单的交叉（即一点交叉）可分两步进行：首先对种群中个体进行随机配对;其次，在配对个体中随机设定交叉处，配对个体彼此交换部分信息。 变异：变异操作是按位（bit）进行的，即把某一位的内容进行变异。变异操作同样也是随机进行的。一般而言，变异概率都取得较小。变异操作是十分微妙的遗传操作，它需要和交叉操作配合使用，目的是挖掘群体中个体的多样性，克服有可能限于局部解的弊病。 图1 遗传算法的基本流程 遗传算法解题的基本步骤如下： Step1：参数设置及种群初始化； Step2：对不可行解进行贪婪修复； Step3：适应度评价； Step4：轮盘赌选择； Step5：交叉； Step6：变异； Step7：对不可行解进行贪婪修复； Step8：适应度评价； Step9：终止条件判断，若未达到终止条件，则转到Step4； Step10：输出结果。 图2 遗传算法具体步骤 四．程序源代码 %遗传算法求解旅行商问题 %初始化 a=[1304 2312;3639 1315;4177 2244;3712 1399;3488 1535;3326 1556;... 3238 1229;4196 1044;4312 790;2864 570;1927 1970;2562 1756;... 2788 1491;2381 1676;1332 695;3715 1678;3918 2179;4061 2370;... 3780 2212;3676 2578;1537 2838;2745 2931;3429 1908;3507 2376;... ];%a：假定的24个城市的坐标 n=100;%n：种群个数 C=200;%C：停止代数 m=2;%m：适配值淘汰加速指数,不宜太大 Pc=0.9;%Pc：交叉概率 Pm=0.2;%Pm： 变异概率 D=distance(a);%生成距离矩阵 [R,Rlength]=GeneTSP(D,a,n,C,m,Pc,Pm); %返回值：最优路径R % 总距离Rlength %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %D：距离矩阵 %n