【MDVRP】遗传算法求解多车辆路径规划问题【含Matlab源码 704期】

本文链接：https://blog.csdn.net/weixin_63266434/article/details/135500750

在这里插入图片描述

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【路径规划】基于matlab遗传算法多车辆路径规划【含Matlab源码 704期】

备注：
点击上面蓝色字体付费专栏Matlab路径规划（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab路径规划（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、遗传算法简介

1 引言
在这里插入图片描述

2 遗传算法理论
2.1 遗传算法的生物学基础

2.2 遗传算法的理论基础

2.3 遗传算法的基本概念

2.4 标准的遗传算法

2.5 遗传算法的特点

2.6 遗传算法的改进方向

3 遗传算法流程
在这里插入图片描述

4 关键参数说明

⛄三、部分源代码

clear
clc
close all
dmax=40; %单车的最大行驶距离
qmax=30; %单车的最大货物携带量
c0=10; %单车的出发成本
c1=1; %单车的行驶成本
X=[18.70,15.29
16.47,8.45
20.07,10.14
19.39,13.37
25.27,14.24
22.00,10.04
25.47,17.02
15.79,15.10
16.60,12.38
14.05,18.12
17.53,17.38
23.52,13.45
19.41,18.13
22.11,12.51
11.25,11.04
14.17,9.76
24.00,19.89
12.21,14.50];
Q=[0 3.0 2.5 5.5 3.0 1.5 4.0 2.5 3.0 2.0 2.5 3.5 3.0 5.0 4.5 2.0 3.5 4.0];
NIND=100; %种群大小
MAXGEN=200;
Pc=0.9; %交叉概率
Pm=0.05; %变异概率
GGAP=0.9; %代沟
D=Distance(X); %生成距离矩阵
N=size(D,1); %客户点数
K=10; %初始的车辆数
%生成初始种群
Chrom=InitPop(NIND,N,K);
%优化
gen=1;
figure(1);
hold on;
box on;
xlim([0,MAXGEN])
title(‘优化过程’)
xlabel(‘代数’)
ylabel(‘最优值’)
ObjV = PathCost(Chrom,Q,D,dmax,qmax,c1,c0,K); %计算总花费
[preObjV,BestIndex] = min(ObjV); %找出最小的花费
BestChrom = Chrom(BestIndex,:);

while gen<MAXGEN
%计算适应度
ObjV=PathCost(Chrom,Q,D,dmax,qmax,c1,c0,K);
line([gen-1,gen],[preObjV,min(ObjV)]);pause(0.001)
[preObjV,BestIndex]=min(ObjV);
BestObjV(gen)=preObjV;
AveObjV(gen)=sum(ObjV)/NIND;
BestChrom(gen,:) = Chrom(BestIndex,:);

FitnV = Fitness(ObjV);
%选择
SelCh1 = Select(Chrom,FitnV,GGAP);
%交叉
SelCh2 = Recombin(SelCh1,Pc);
 %变异
SelCh3 = Mutate(SelCh2,Pm);
%逆转操作
SelCh4 = Reverse(SelCh3,D,Q,dmax,qmax,c1,c0,K);
%重新插入新的种群
Chrom =Reins(Chrom,SelCh4,ObjV);
gen = gen+1;

end
%画出最优解的路线图
ObjV=PathCost(Chrom,Q,D,dmax,qmax,c1,c0,K);
[minObjV,minInd]=min(ObjV);
DrawPath(Chrom(minInd(1)😅,X);
%输出最优解
disp(‘最优服务顺序:’)
p=OutputPath(Chrom(minInd(1)😅);
disp([‘总花费：’,num2str(minObjV)]);
s=0;
R=Chrom(minInd(1)😅;
for i=1:size(R,2)-1
s=s+D(R(i),R(i+1));
end

function DrawPath(Chrom,X)
%%画路线图函数
%输入
%Chrom 待画路线
%X 各服务点的坐标位置

R=Chrom;
figure;
hold on
plot(X(:,1),X(:,2),‘o’,‘color’,[0.5,0.5,0.5])
plot(X(Chrom(1,1),1),X(Chrom(1,1),2),‘rv’,‘MarkerSize’,20)
for i=1:size(X,1)
text(X(i,1)+0.05,X(i,2)+0.05,num2str(i),‘color’,[1,0,0]);
end
A=X(R,:);
row=size(A,1);
for i=2:row
[arrowx,arrowy]=dsxy2figxy(gca,A(i-1:i,1),A(i-1:i,2));
annotation(‘textarrow’,arrowx,arrowy,‘HeadWidth’,8,‘color’,[0,0,1]);
end
function varargout = dsxy2figxy(varargin)
if length(varargin{1}) == 1 && ishandle(varargin{1}) …
&& strcmp(get(varargin{1},‘type’),‘axes’)
hAx = varargin{1};
varargin = varargin(2:end);
else
hAx = gca;
end;
if length(varargin) == 1
pos = varargin{1};
else
[x,y] = deal(varargin{:});
end
axun = get(hAx,‘Units’);
set(hAx,‘Units’,‘normalized’);
axpos = get(hAx,‘Position’);
axlim = axis(hAx);
axwidth = diff(axlim(1:2));
axheight = diff(axlim(3:4));
if exist(‘x’,‘var’)
varargout{1} = (x - axlim(1)) * axpos(3) / axwidth + axpos(1);
varargout{2} = (y - axlim(3)) * axpos(4) / axheight + axpos(2);
else
pos(1) = (pos(1) - axlim(1)) / axwidth * axpos(3) + axpos(1);
pos(2) = (pos(2) - axlim(3)) / axheight * axpos(4) + axpos(2);
pos(3) = pos(3) * axpos(3) / axwidth;
pos(4) = pos(4) * axpos(4 )/ axheight;
varargout{1} = pos;
end