一、VRP简介

1 VRP基本原理

车辆路径规划问题(Vehicle Routing Problem,VRP)是运筹学里重要的研究问题之一。VRP关注有一个供货商与K个销售点的路径规划的情况,可以简述为:对一系列发货点和收货点,组织调用一定的车辆,安排适当的行车路线,使车辆有序地通过它们,在满足指定的约束条件下(例如:货物的需求量与发货量,交发货时间,车辆容量限制,行驶里程限制,行驶时间限制等),力争实现一定的目标(如车辆空驶总里程最短,运输总费用最低,车辆按一定时间到达,使用的车辆数最小等)。

VRP的图例如下所示:

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_优化算法

2 问题属性与常见问题

车辆路径问题的特性比较复杂,总的来说包含四个方面的属性:

(1)地址特性包括:车场数目、需求类型、作业要求。

(2)车辆特性包括:车辆数量、载重量约束、可运载品种约束、运行路线约束、工作时间约束。

(3)问题的其他特性。

(4)目标函数可能是总成本极小化,或者极小化最大作业成本,或者最大化准时作业。

3 常见问题有以下几类:

(1)旅行商问题

(2)带容量约束的车辆路线问题(CVRP)

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_02

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_matlab_03

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_matlab_04

该模型很难拓展到VRP的其他场景,并且不知道具体车辆的执行路径,因此对其模型继续改进。

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_05

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_06

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_07

(3)带时间窗的车辆路线问题

由于VRP问题的持续发展,考虑需求点对于车辆到达的时间有所要求之下,在车辆途程问题之中加入时窗的限制,便成为带时间窗车辆路径问题(VRP with Time Windows, VRPTW)。带时间窗车辆路径问题(VRPTW)是在VRP上加上了客户的被访问的时间窗约束。在VRPTW问题中,除了行驶成本之外, 成本函数还要包括由于早到某个客户而引起的等待时间和客户需要的服务时间。在VRPTW中,车辆除了要满足VRP问题的限制之外,还必须要满足需求点的时窗限制,而需求点的时窗限制可以分为两种,一种是硬时窗(Hard Time Window),硬时窗要求车辆必须要在时窗内到达,早到必须等待,而迟到则拒收;另一种是软时窗(Soft Time Window),不一定要在时窗内到达,但是在时窗之外到达必须要处罚,以处罚替代等待与拒收是软时窗与硬时窗最大的不同。

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_08

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_matlab_09

模型2(参考2017 A generalized formulation for vehicle routing problems):

该模型为2维决策变量

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_10

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_11

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_12

(4)收集和分发问题

(5)多车场车辆路线问题

参考(2005 lim,多车场车辆路径问题的遗传算法_邹彤, 1996 renaud)

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_优化算法_13

由于车辆是同质的,这里的建模在变量中没有加入车辆的维度。

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_14

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_遗传算法_15

(6)优先约束车辆路线问题

(7)相容性约束车辆路线问题

(8)随机需求车辆路线问题

4 解决方案

(1)数学解析法

(2)人机交互法

(3)先分组再排路线法

(4)先排路线再分组法

(5)节省或插入法

(6)改善或交换法

(7)数学规划近似法

(8)启发式算法

5 VRP与VRPTW对比

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_遗传算法_16

二、遗传算法简介

1 引言

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_遗传算法_17

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_18

2 遗传算法理论

2.1 遗传算法的生物学基础

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_优化算法_19

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_优化算法_20

2.2 遗传算法的理论基础

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_21

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_22

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_23

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_遗传算法_24

2.3 遗传算法的基本概念

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_25

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_matlab_26

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_matlab_27

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_28

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_29

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_遗传算法_30

2.4 标准的遗传算法

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_优化算法_31

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_遗传算法_32

2.5 遗传算法的特点

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_33

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_matlab_34

2.6 遗传算法的改进方向

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_算法_35

3 遗传算法流程

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_36

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_37

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_matlab_38

4 关键参数说明

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_遗传算法_39

三、部分源代码

function [ globalMin, optPop1, optPop2] =dtsp_ga_basic( popSize, numIter, xy, range, speed )

% 如果为空参数
showprogress=false;
if nargin < 5 % 如果输入变量小于 6,则使用默认值
clear; clc;
nStops=6; numIter=500; speed=2; range=2;
xy = 1*rand([nStops, 2]);


popSize=length(xy)*5;
showprogress=true;
end

% 将变量初始化为整数/标准大小,创建距离矩阵
nPoints = size(xy,1);
popSize = 5*floor(popSize/5);
numIter = max(1,round(real(numIter(1))));
meshg = meshgrid(1:nPoints);
dmat = reshape(sqrt(sum((xy(meshg,:)-xy(meshg',:)).^2,2)),nPoints,nPoints);

%初始化种群
[n, ~]= size(xy);
pop = zeros(popSize,n);
pop2 = zeros(popSize,n);
pop(1,:) = 1:n;
for k = 2:popSize
pop(k,:) = randperm(n);
end


% 运行GA
globalMin = Inf;
totalDist = zeros(1,popSize);
distHistory = zeros(1,numIter);
tmpPop = zeros(5,n);
newPop = zeros(popSize,n);

四、运行结果

【配送路径规划】基于matlab遗传算法求解无人机和货车配送路径规划问题【含Matlab源码 361期】_自动驾驶_40

五、matlab版本及参考文献

1 matlab版本

2014a

2 参考文献

[1] 包子阳,余继周,杨杉.智能优化算法及其MATLAB实例(第2版)[M].电子工业出版社,2016.

[2]张岩,吴水根.MATLAB优化算法源代码[M].清华大学出版社,2017.