实验六遗传算法求解TSP问题实验Word下载.docx

1、适应度评分函数,为给定假设赋予一个评估分数 Fitness_threshold:指定终止判据的阈值 p:群体中包含的假设数量 r:每一步中通过交叉取代群体成员的比例 m:变异率初始化群体:P随机产生的p个假设评估:对于P中的每一个h,计算Fitness（h）当maxFitness（h）Fitness_threshold,做产生新的一代Ps: （1）选择:用概率方法选择P的（1-r）p个成员加入Ps.从P中选择假设hi的概率用下面公式计算: （2）交叉:根据上面给出的,从P中按概率选择r（p/2）对假设.对于每对假设,应用交叉算子产生两个后代.把所有的后代加入Ps （3）变异:使用均匀的概率从P

2、s中选择m%的成员.对于选出的每个成员,在它表示中随机选择一个为取反 （4）更新:PPs （5）评估:对于P中的每个h计算Fitness（h）从P中返回适应度最高的假设3. TSP问题的遗传算法设计与实现对于n个城市的问题，每个个体即每个解的长度为n，用s行, t列的pop矩阵，表示初始群体，s表示初始群体的个数，t为n+1，矩阵的每一行的前n个元素表示城市编码，最后一个元素表示这一路径的长度。城市的位置可以手动输入，使用一个NN矩阵D存储，D（i，j）代表城市i和城市j之间的距离。 D（i，j）=sqrt（Xi-Xj）.2+（Yi-Yj）.2）。在TSP的求解中，可以直接用距离总和作为适应度

3、函数。个体的路径长度越小，所得个体优越，距离的总和越大，适应度越小，进而探讨求解结果是否最优。选择就是从群体中选择优胜个体、淘汰劣质个体的操作，它是建立在群体中个体适应度评估基础上。这里采用方法是最优保存方法。本实例中交叉采用部分匹配交叉策略，先随机选取两个交叉点，然后将两交叉点中间的基因段互换，将互换的基因段以外的部分中与互换后基因段中元素冲突的用另一父代的相应位置代替，直到没有冲突。变异操作是以变异概率Pm对群体中个体串某些基因位上的基因值作变动，若变异后子代的适应度值更加优异，则保留子代染色体，否则，仍保留父代染色体。这有助于增加种群的多样性，避免早熟收敛（非全局最优）。判断结束准则是固

4、定指定了迭代的次数当算法达到迭代次数时，算法结束，输出当前的最优解。在根据适配值计算并选择的时候，记录下来的当前最优值，在变异后加入跟新的群体，保证新的迭代循环中TSP解越来越好（不会变差）。 在选择的一种操作是拿最优的K个替换最差的K个个体，本例是按适配值选择，并使群体数目变少，当每次变异操作后，产生随机路径补充群体是群体数目不变，再次循环，一定程度上防止因初始群体的选择问题而陷入局部最优。4. TSP问题的遗传算法的具体步骤解最短路径的遗传算法如下： Generatep（n）;表示程序开始时要首先产生一个群体，群体个数为nEvaluatep（h）;表示计算每个个体适应度，h是种群中的一个个

5、体Repeat roof Generations times;重复下面的操作，直到满足条件为止Select p（h） from p（n-1）;表示从前一代群体中选择一对双亲，用于交叉、变异 操作，P（n）代表第n代群体Crossover and mutation p（n）;进行交叉和变异操作Learningp（n）;自学习过程计算新生成的种群中每个个体的适应度End;具体流程图如下所示：流程图5.遗传算法求解不同规模的TSP问题的算法性能（1） 遗传算法执行方式说明： 适应度值计算方法：当前路线的路径长度 个体选择概率分配方法：适应度比例方法 选择个体方法：轮盘赌选择 交叉类型：PMX交叉 变

6、异类型: 两点互换变异（2）实验模拟结果：城市个数时间（ms）51692510166301518833202259625241593030289353523940386084540032504375755477466058143655994270643617571417图1-1（3）分析 由图1-1可知，遗传算法执行时间随着TSP问题规模的增大而增大，并且大致为线性增长。五、不同参数下的计算结果对比（1）种群规模对算法结果的影响实验次数：最大迭代步数:100交叉概率：0.85变异概率：0.15表1-1种群规模适应度值最优路径25.2644-5-8-7-6-3-1-0-9-226.34282-9

7、-1-0-3-6-7-5-8-425.16521-3-6-7-5-8-4-2-9-00-1-3-6-7-5-8-4-2-9809-0-1-3-6-7-5-8-4-21-0-9-2-4-8-5-7-6-31505-8-4-2-9-0-1-3-6-72002503-1-0-9-2-4-8-5-7-6300如表1-1所示，显然最短路径为25.1652m,最优路径为1-0-9-1-3-6-7-5-8-4-2或3-1-0-9-2-4-8-5-7-6，注意到这是一圈，顺时针或者逆时针都可以。当种群规模为10，20时，并没有找到最优解。（2）交叉概率对算法结果的影响种群规模：实验结果：表1-2交叉概率最好适

8、应度最差适应度平均适应度最优解运行时间0.00128.044736.656732.60029-2-6-0-5-4-8-7-3-13100.0127.093534.994332.14957-8-3-1-9-2-6-0-5-42600.135.303331.93727-3-1-9-2-6-0-5-4-834.117531.21830-5-4-8-7-3-1-9-2-62700.228.710833.951230.90353-1-9-2-6-5-0-4-7-82800.2535.162330.74561-3-7-8-4-5-0-6-2-90.331.994129.94288-3-1-9-2-6-0-

9、5-4-72900.3532.808530.99459-1-3-8-7-4-5-0-6-20.432.531330.15341-3-8-7-4-5-0-6-2-92790.4533.201430.17574560.528.093433.630730.90265-0-2-6-9-1-3-8-7-46630.5533.523329.13041-9-2-6-0-5-4-7-8-35200.633.251230.78363-1-9-2-6-0-5-4-7-85460.6533.700330.93715-4-8-7-3-1-9-2-6-05960.732.092729.95025710.7532.448

10、830.36995590.832.155129.93827-4-5-0-6-2-9-1-3-835834.539930.35945-0-6-2-9-1-3-8-7-43600.932.627330.696-0-5-4-7-8-3-1-9-23750.9532.467229.9196-2-9-1-3-8-7-4-5-0476（注:红色表示非最优解）在该情况下，交叉概率过低将使搜索陷入迟钝状态，得不到最优解。（3）变异概率对算法结果的影响表1-3变异概率29.471734.73232.49110-6-2-1-9-3-8-7-4-524529.044634.659132.37148-4-5-0-2-

11、6-9-1-3-727434.01130.941732.09330.256824632.234930.31448-7-4-5-0-6-2-9-1-328232.71830.15724-5-0-6-2-9-1-3-8-730.28540-5-4-7-8-3-1-9-2-625233.316730.774826634.370531.30412-0-5-4-8-7-3-1-9-636231.37429.68162-6-0-5-4-7-8-3-1-943832.375230.22112-9-1-3-8-7-4-5-0-643133.381930.662349230.361-3-8-7-4-5-0-2-

12、6-941732.749130.020143432.423830.7851-3-8-7-4-0-5-6-2-943231.892830.245147531.613530.34719-1-3-8-7-4-5-0-2-632729.66233.239231.15852-9-1-3-7-8-4-0-5-631432.038730.415239631.303630.00679-1-3-7-8-4-5-0-6-2436又表1-3可知，当变异概率过大或过低都将导致无法得到最优解。注：（2）（3）的实验数据与（1）的实验数据不同，详见附录。六、不同变异策略和个体选择概率分配策略对算法结果的影响（1）两点互换

13、变异与插入变异的比较： 试验次数（CASNUM）：10 城市数（POINTCNT）: 种群规模（POPSIZE）： 最大迭代步数（GENERATIONS）: 交叉概率（PC）： 变异概率（PM）:a. 变异类型:表1-4两点互换变异程序结果序号130.422929.08916-2-0-5-4-7-8-3-1-91199231.141728.98411678330.422829.06041940430.370328.8787175631.061929.07551885631.158929.39421936729.76486-2-9-1-3-7-8-4-5-01772831.347529.8415

14、4-5-0-2-6-9-1-3-7-81980930.614329.0590-6-2-9-1-3-8-7-4-530.558529.08119-2-6-0-5-4-7-8-3-118721131.017129.426415171229.241415411332.025529.07891431.51628.8906134529.022613771630.408128.908118531729.331615221830.020328.524316011931.140429.5672-9-1-3-7-8-4-5-0-6160929.53591311平均值27.336130.878229.187716

15、57b. 变异类型: 插入变异表1-5插入变异程序结果31.475328.8453138828.9168135529.663128.902163730.524129.51194-5-0-6-2-9-1-3-7-8116431.057529.4682124528.5546122230.820529.7483-1-9-2-6-0-5-4-8-7114829.3907174230.42328.6878206428.72151829.3282124030.52328.5544120429.0508173431.117729.5905153229.1904148328.8061128231.763929.4591148529.161428.5974150728.8036123427.1886230.646529.06431439分析： 两点互换变异20次模拟中，4次得到非最优解；而插入变异只有2次；插入变异的最好适应度平均值比两点互换变异小0.14755，最差适应度平均值和总的适应度平均值都比两点互换下，并且在Release下，运行时间前者比后者快218.3ms。可见在该条件下（交