功能描述：
联系：highspeedlogic
QQ ：1224848052
微信：HuangL1121
邮箱：1224848052@qq.com
网站：http://www.mat7lab.com/
网站：http://www.hslogic.com/
微信扫一扫：
车辆路径问题（VRP）是运筹学里重要的研究问题之一。VRP关注有一个供货商与K个销售点的路径规划的情况，可以简述为：给定一个或多个中心（中心车库）一个车辆集合和一个顾客集合，车辆和顾客各有自己的属性，每辆车都有容量，所载的货物不能超过它的容量。
%D是距离矩阵，n为种群个数
%C为停止代数，遗传到第 C代时程序停止,C的具体取值视问题的规模和耗费的时间而定
%m为适值淘汰加速指数,最好取为1,2,3,4,不宜太大
%交叉概率Pc,变异概率Pm
%R为最短路径,Rlength为路径长度
function [R,Rlength]=geneticVRP(D,demand,n,C,m,Pc,Pm)
[N,NN]=size(D);%(31*31)
farm=zeros(n,N);%用于存储种群
for i=1:n
% flag=0;
%while flag==0
%tem=randperm(N); %随机生成初始种群
%if validate(tem,demand,N)==1
% flag=1;
%end
%end
farm(i,:)=randperm(N);
end
R=farm(1,:);%一个随机解(个体)
%farm(1,:)=R;
len=zeros(n,1);%存储路径长度
fitness=zeros(n,1);%存储适配值
counter=0;
while counter<C
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);%返回的是在len中路径最短的路径坐标(i,1)
R=farm(rr(1,1),:);%更新最短路径
%disp('~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
FARM=farm;%优胜劣汰，nn记录了复制的个数
%选择,
K=30;
[aa,bb]=size(FARM);
FARM2=FARM;
len2=len;
[len]=sort(len);
for i=1:aa
tt= find(len2==len(i,1));
FARM(i,:)=FARM2(tt(1,1),:);
end
for i=1:K
j=aa+1-i;
FARM(j,:)=FARM(i,:);
end
% 交叉操作
[aa,bb]=size(FARM);
FARM2=FARM;
for i=1:2:aa
if Pc>rand&&i<aa %交叉概率Pc
A=FARM(i,:);
B=FARM(i+1,:);
[A,B]=intercross(A,B);
FARM(i,:)=A;
FARM(i+1,:)=B;
end
end
%交叉检验 (可省去)
for i=1:aa
if myLength(D,FARM(i,:))>myLength(D,FARM2(i,:))
FARM(i,:)=FARM2(i,:);
end
end
clear FARM2
[aa,bb]=size(FARM); %aa=nn2
% 变异
FARM2=FARM;
for i=1:aa
if Pm>=rand
FARM(i,:)=mutate(FARM(i,:));
end
end
%变异检验(可省略)
for i=1:aa
if myLength(D,FARM(i,:))>myLength(D,FARM2(i,:))
FARM(i,:)=FARM2(i,:);
end
end
clear FARM2
%群体的更新
%FARM2=zeros(n-aa+1,N);
%if n-aa>=1
% for i=1:n-aa
% FARM2(i,:)=randperm(N);%随机生成n-aa种群
% end
% end
FARM=[R;FARM];%将随机产生的n-aa个体加入从后面种群,将上次迭代的最优解从前面加入种群
[aa,bb]=size(FARM);
%disp('~~~~~~~~~~~~~~~~4~~~~~~~验证zong~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
%保持种群规模为n
if aa>n
FARM=FARM(1:n,:);
end
%disp('~~~~~~~~~~~~~~~~~~~5~~~~验证zong~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~')
%更新farm
farm=FARM;
clear FARM
%更新迭代次数
counter=counter+1 ;
end
%结果输出
Rlength=myLength(D,R)
R
Rlength=myLength(D,R)%结果输出
