TSP问题遗传算法通用Matlab程序.docVIP

下载本文档

9
0
约 5页
2017-08-10 发布于广东
举报
版权申诉

TSP问题遗传算法通用Matlab程序.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

TSP问题遗传算法通用Matlab程序（附图）(2007-03-10 09:56:14) 程序一：主程序%TSP问题（又名：旅行商问题，货郎担问题）遗传算法通用matlab程序%D是距离矩阵，n为种群个数%参数a是中国31个城市的坐标%C为停止代数，遗传到第 C代时程序停止,C的具体取值视问题的规模和耗费的时间而定%m为适应值归一化淘汰加速指数,最好取为1,2,3,4,不宜太大%alpha为淘汰保护指数,可取为0~1之间任意小数,取1时关闭保护功能,建议取0.8~1.0之间的值%R为最短路径,Rlength为路径长度function [R,Rlength]=geneticTSP(D,a,n,C,m,alpha)[N,NN]=size(D);farm=zeros(n,N);%用于存储种群for i=1:n??? farm(i,:)=randperm(N);%随机生成初始种群endR=farm(1,:);subplot(1,3,1)scatter(a(:,1),a(:,2),x)pause(1)subplot(1,3,2)plotaiwa(a,R)pause(1)farm(1,:)=R;len=zeros(n,1);%存储路径长度fitness=zeros(n,1);%存储归一化适应值counter=0;while counter??? for i=1:n??????? len(i,1)=myLength(D,farm(i,:));%计算路径长度??? end??? maxlen=max(len);??? minlen=min(len);??? fitness=fit(len,m,maxlen,minlen);%计算归一化适应值??? rr=find(len==minlen);??? R=farm(rr(1,1),:);%更新最短路径?????? FARM=farm;%优胜劣汰，nn记录了复制的个数??? nn=0;??? for i=1:n??????? if fitness(i,1)=alpha*rand??????????? nn=nn+1;??????????? FARM(nn,:)=farm(i,:);??????? end??? end??? FARM=FARM(1:nn,:);?????? [aa,bb]=size(FARM);%交叉和变异??? while aa??????? if nn=2??????????? nnper=randperm(2);??????? else??????????? nnper=randperm(nn);??????? end??????? A=FARM(nnper(1),:);??????? B=FARM(nnper(2),:);??????? [A,B]=intercross(A,B);??????? FARM=[FARM;A;B];??????? [aa,bb]=size(FARM);??? end??? if aan??????? FARM=FARM(1:n,:);%保持种群规模为n??? end?????? farm=FARM;??? clear FARM??? counter=counter+1???endRlength=myLength(D,R);subplot(1,3,3)plotaiwa(a,R) 程序二：计算邻接矩阵%输入参数a是中国31个城市的坐标%输出参数D是无向图的赋权邻接矩阵function D=ff01(a)[c,d]=size(a);D=zeros(c,c);for i=1:c??? for j=i:c??????? bb=(a(i,1)-a(j,1)).^2+(a(i,2)-a(j,2)).^2;??????? D(i,j)=bb^(0.5);??????? D(j,i)=D(i,j);??? endend 程序三：计算归一化适应值%计算归一化适应值的子程序function fitness=fit(len,m,maxlen,minlen)fitness=len;for i=1:length(len)??? fitness(i,1)=(1-((len(i,1)-minlen)/(maxlen-minlen+0.0001))).^m;end 程序四：交叉和变异的子程序%交叉算法采用的是由Goldberg和Lingle于1985年提出的PMX(部分匹配交叉)function [a,b]=intercross(a,b)L=length(a);if L=10%确定交叉宽度??? W=9