实验六遗传算法求解TSP问题实验Word文档下载推荐.doc
- 文档编号:13167313
- 上传时间:2022-10-07
- 格式:DOC
- 页数:24
- 大小:331KB
实验六遗传算法求解TSP问题实验Word文档下载推荐.doc
《实验六遗传算法求解TSP问题实验Word文档下载推荐.doc》由会员分享,可在线阅读,更多相关《实验六遗传算法求解TSP问题实验Word文档下载推荐.doc(24页珍藏版)》请在冰豆网上搜索。
通过模拟自然选择和自然遗传过程中发生的繁殖、交叉和基因突变现象,在每次迭代中都保留一组候选解,并按某种指标从解群中选取较优的个体,利用遗传算子(选择、交叉和变异)对这些个体进行组合,产生新一代的候选解群,重复此过程,直到满足某种收敛指标为止。
遗传算法在本质上是一种不依赖具体问题的直接搜索方法,是一种求解问题的高效并行全局搜索方法。
其假设常描述为二进制位串,位串的含义依赖于具体应用。
搜索合适的假设从若干初始假设的群体集合开始。
当前种群成员通过模仿生物进化的方式来产生下一代群体,如随机变异和交叉。
每一步,根据给定的适应度评估当前群体的假设,而后使用概率方法选出适应度最高的假设作为产生下一代的种子。
下面介绍几个遗传算法的几个基本概念:
(1)染色体(Chromosome):
在使用遗传算法时,需要把问题的解编成一个适合的码子。
这种具有固定结构的符号串既是染色体,符号串的每一位代表一个基因。
符号串的总位数成为染色体的长度,一个染色体就代表问题的一个解,每个染色体也被称为一个个体。
(2)群体(Population):
每代所产生的染色体总数成为群体,一个群体包含了该问题在这一代的一些解的集合。
(3)适应度(Fitness):
对群体中每个染色体进行编码后,每个个体对应一个具体问题的解,而每个解对应于一个函数值。
该函数值即适应函数,就是衡量染色体对环境适应度的指标,也是反映实际问题的目标函数
在前一代群体的基础上产生新一代群体的工作成为遗传操作,基本的遗传操作有:
(1)选择(Select):
按一定的概率从上代群体中选择M对个体作为双亲,直接拷贝到下一代,染色体不发生变化。
(2)交叉(Crossover):
对于选中进行繁殖的两个染色体X,Y,以X,Y为双亲作交叉操作,从而产生两个后代X1,Y1.
(3)变异(Mutation):
对于选中的群体中的个体(染色体),随机选取某一位进行取反运算,即将该染色体码翻转。
用遗传算法求解的过程是根据待解决问题的参数集进行编码,随机产生一个种群,计算适应函数和选择率,进行选择、交叉、变异操作。
如果满足收敛条件,此种群为最好个体,否则,对产生的新一代群体重新进行选择、交叉、变异操作,循环往复直到满足条件。
遗传算法原型:
GA(Fitness,Fitness_threshold,p,r,m)
Fitness:
适应度评分函数,为给定假设赋予一个评估分数
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)对假设.对于每对假设<
h1,h2>
应用交叉算子产生两个后代.把所有的后代加入Ps
(3)变异:
使用均匀的概率从Ps中选择m%的成员.对于选出的每个成员,在它表示中随机选择一个为取反
(4)更新:
P←Ps
(5)评估:
对于P中的每个h计算Fitness(h)
从P中返回适应度最高的假设
3.TSP问题的遗传算法设计与实现
对于n个城市的问题,每个个体即每个解的长度为n,用s行,t列的pop矩阵,表示初始群体,s表示初始群体的个数,t为n+1,矩阵的每一行的前n个元素表示城市编码,最后一个元素表示这一路径的长度。
城市的位置可以手动输入,使用一个N×
N矩阵D存储,D(i,j)代表城市i和城市j之间的距离。
D(i,j)=sqrt((Xi-Xj).^2+(Yi-Yj).^2)。
在TSP的求解中,可以直接用距离总和作为适应度函数。
个体的路径长度越小,所得个体优越,距离的总和越大,适应度越小,进而探讨求解结果是否最优。
选择就是从群体中选择优胜个体、淘汰劣质个体的操作,它是建立在群体中个体适应度评估基础上。
这里采用方法是最优保存方法。
本实例中交叉采用部分匹配交叉策略,先随机选取两个交叉点,然后将两交叉点中间的基因段互换,将互换的基因段以外的部分中与互换后基因段中元素冲突的用另一父代的相应位置代替,直到没有冲突。
变异操作是以变异概率Pm对群体中个体串某些基因位上的基因值作变动,若变异后子代的适应度值更加优异,则保留子代染色体,否则,仍保留父代染色体。
这有助于增加种群的多样性,避免早熟收敛(非全局最优)。
判断结束准则是固定指定了迭代的次数当算法达到迭代次数时,算法结束,输出当前的最优解。
在根据适配值计算并选择的时候,记录下来的当前最优值,在变异后加入跟新的群体,保证新的迭代循环中TSP解越来越好(不会变差)。
在选择的一种操作是拿最优的K个替换最差的K个个体,本例是按适配值选择,并使群体数目变少,当每次变异操作后,产生随机路径补充群体是群体数目不变,再次循环,一定程度上防止因初始群体的选择问题而陷入局部最优。
4.TSP问题的遗传算法的具体步骤
解最短路径的遗传算法如下:
Generate[p(n)];
表示程序开始时要首先产生一个群体,群体个数为n
Evaluate[p(h)];
表示计算每个个体适应度,h是种群中的一个个体
RepeatroofGenerationstimes;
重复下面的操作,直到满足条件为止
Selectp(h)fromp(n-1);
表示从前一代群体中选择一对双亲,用于交叉、变异操作,P(n)代表第n代群体
Crossoverandmutationp(n);
进行交叉和变异操作
Learning[p(n)];
自学习过程
计算新生成的种群中每个个体的适应度
End;
具体流程图如下所示:
流程图
5.遗传算法求解不同规模的TSP问题的算法性能
(1)遗传算法执行方式说明:
l适应度值计算方法:
当前路线的路径长度
l个体选择概率分配方法:
适应度比例方法
l选择个体方法:
轮盘赌选择
l交叉类型:
PMX交叉
l变异类型:
两点互换变异
(2)实验模拟结果:
城市个数
时间(ms)
5
16925
10
16630
15
18833
20
22596
25
24159
30
30289
35
35239
40
38608
45
40032
50
43757
55
47746
60
58143
65
59942
70
64361
75
71417
图1-1
(3)分析
由图1-1可知,遗传算法执行时间随着TSP问题规模的增大而增大,并且大致为线性增长。
五、不同参数下的计算结果对比
(1)种群规模对算法结果的影响
实验次数:
最大迭代步数:
100
交叉概率:
0.85
变异概率:
0.15
表1-1
种群规模
适应度值
最优路径
25.264
4-5-8-7-6-3-1-0-9-2
26.3428
2-9-1-0-3-6-7-5-8-4
25.1652
1-3-6-7-5-8-4-2-9-0
0-1-3-6-7-5-8-4-2-9
80
9-0-1-3-6-7-5-8-4-2
1-0-9-2-4-8-5-7-6-3
150
5-8-4-2-9-0-1-3-6-7
200
250
3-1-0-9-2-4-8-5-7-6
300
如表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
交叉概率
最好适应度
最差适应度
平均适应度
最优解
运行时间
0.001
28.0447
36.6567
32.6002
9-2-6-0-5-4-8-7-3-1
310
0.01
27.0935
34.9943
32.1495
7-8-3-1-9-2-6-0-5-4
260
0.1
35.3033
31.9372
7-3-1-9-2-6-0-5-4-8
34.1175
31.2183
0-5-4-8-7-3-1-9-2-6
270
0.2
28.7108
33.9512
30.9035
3-1-9-2-6-5-0-4-7-8
280
0.25
35.1623
30.7456
1-3-7-8-4-5-0-6-2-9
0.3
31.9941
29.9428
8-3-1-9-2-6-0-5-4-7
290
0.35
32.8085
30.9945
9-1-3-8-7-4-5-0-6-2
0.4
32.5313
30.1534
1-3-8-7-4-5-0-6-2-9
279
0.45
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 遗传 算法 求解 TSP 问题
![提示](https://static.bdocx.com/images/bang_tan.gif)