【CVRP】遗传算法求解带容量和体积的车辆路径规划问题【含Matlab源码 3275期】

Matlab武动乾坤

已于 2024-11-12 11:58:01 修改

阅读量935

点赞数 12

分类专栏： Matlab路径规划（进阶版）文章标签： matlab

于 2024-02-22 23:28:36 首次发布

本文链接：https://blog.csdn.net/KeepingMatlab/article/details/136245140

版权

Matlab路径规划（进阶版）专栏收录该内容

998 篇文章 40 订阅

订阅专栏

💥💥💥💥💥💥💥💥💞💞💞💞💞💞💞💞💞Matlab武动乾坤博客之家💞💞💞💞💞💞💞💞💞💥💥💥💥💥💥💥💥
🚀🚀🚀🚀🚀🚀🚀🚀🚀🚀🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚀🚀🚀🚀🚀🚀🚀🚀🚀🚀
在这里插入图片描述
🔊博主简介：985研究生，Matlab领域科研开发者；

🚅座右铭：行百里者，半于九十。

🏆代码获取方式：
CSDN Matlab武动乾坤—代码获取方式

更多Matlab路径规划仿真内容点击👇
①Matlab路径规划（进阶版）

⛳️关注CSDN Matlab武动乾坤，更多资源等你来！！

⛄一、VRP简介

1 VRP基本原理
车辆路径规划问题(Vehicle Routing Problem，VRP)是运筹学里重要的研究问题之一。VRP关注有一个供货商与K个销售点的路径规划的情况，可以简述为：对一系列发货点和收货点，组织调用一定的车辆，安排适当的行车路线，使车辆有序地通过它们，在满足指定的约束条件下（例如：货物的需求量与发货量，交发货时间，车辆容量限制，行驶里程限制，行驶时间限制等），力争实现一定的目标（如车辆空驶总里程最短，运输总费用最低，车辆按一定时间到达，使用的车辆数最小等）。
VRP的图例如下所示：
在这里插入图片描述
2 问题属性与常见问题
车辆路径问题的特性比较复杂，总的来说包含四个方面的属性：
（1）地址特性包括：车场数目、需求类型、作业要求。
（2）车辆特性包括：车辆数量、载重量约束、可运载品种约束、运行路线约束、工作时间约束。
（3）问题的其他特性。
（4）目标函数可能是总成本极小化，或者极小化最大作业成本，或者最大化准时作业。

3 常见问题有以下几类：
（1）旅行商问题
（2）带容量约束的车辆路线问题(CVRP)
在这里插入图片描述

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

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

模型2(参考2017 A generalized formulation for vehicle routing problems)：
该模型为2维决策变量

（4）收集和分发问题
（5）多车场车辆路线问题
参考(2005 lim，多车场车辆路径问题的遗传算法_邹彤, 1996 renaud)

由于车辆是同质的，这里的建模在变量中没有加入车辆的维度。
在这里插入图片描述

（6）优先约束车辆路线问题
（7）相容性约束车辆路线问题
（8）随机需求车辆路线问题

4 解决方案
（1）数学解析法
（2）人机交互法
（3）先分组再排路线法
（4）先排路线再分组法
（5）节省或插入法
（6）改善或交换法
（7）数学规划近似法
（8）启发式算法

5 VRP与VRPTW对比
在这里插入图片描述

⛄二、遗传算法简介

1 遗传算法
遗传算法是一种基于生物进化论模型的优化算法，通过模拟生物进化的过程，通过复制、交叉、突变等操作产生下一代的解，并逐步淘汰掉适应度函数值低的解，增加适应度函数值高的解。遗传算法可以用于解决各种优化问题，如函数优化、组合优化、机器学习等。在遗传算法中，个体的适应度函数值越高，就越有可能被选择为下一代的父代，从而进化出更优秀的解。遗传算法的优点是可以在大规模搜索空间中找到全局最优解，但是也存在一些缺点，如收敛速度慢、参数设置困难等。

2 遗传算法步骤
遗传算法是一种模拟自然进化过程的优化算法，其步骤如下：
（1）初始化种群：随机生成一定数量的个体，每个个体都是由若干个基因组成的染色体。
（2）评估适应度：对于每个个体，通过一个适应度函数来评估其适应度，即该个体在解决问题中的表现好坏。
（3）选择操作：根据适应度函数的值，选择一部分个体作为父代，用于产生下一代个体。
（4）交叉操作：对于选出的父代个体，进行交叉操作，生成新的个体。
（5）变异操作：对于新生成的个体，进行变异操作，引入新的基因，增加种群的多样性。
（6）重复步骤2-5，直到达到预设的终止条件（如达到最大迭代次数或找到满足要求的解）。

3 遗传算法求解带容量和体积的车辆路径规划问题
遗传算法是一种基于生物进化原理的优化算法，可以用于求解带容量和体积的车辆路径规划问题。具体实现过程如下：
（1）定义染色体编码方式，将车辆路径规划问题转化为染色体的编码问题。
（2）初始化种群，随机生成一定数量的染色体。
（3）评估适应度，根据染色体编码计算每个染色体的适应度。
（4）选择操作，根据适应度选择一定数量的染色体作为下一代的父代。
（5）交叉操作，对父代染色体进行交叉操作，生成新的染色体。
（6）变异操作，对新生成的染色体进行变异操作，增加种群的多样性。
（7）评估适应度，计算新生成的染色体的适应度。
（8）选择操作，根据适应度选择一定数量的染色体作为下一代的种群。
（9）重复步骤4-8，直到满足停止条件。
在带容量和体积的车辆路径规划问题中，染色体编码可以采用基因串编码，其中每个基因表示一个客户点，基因串表示车辆的路径。同时，需要引入容量和体积约束条件，确保每个车辆的容量和体积不超过限制。在评估适应度时，可以考虑车辆的行驶距离和满足约束条件的程度。

⛄三、部分源代码

clc
close all
clear all
Maxgen=500;
%每次就改改上面这个迭代次数就好了！
G=[2000 2800 3200 3700 4000 5000 5800 7000];%每辆车的额定载重量
V=[12 19 25 32 30 38 42 50];%每辆车的额定容积
g=[251 306 212 158 166 251 183 288 216 158 121 255 188 114 161 358 232 272 69 116.5 181 196 219 248 158];%每个顾客需求总重量
v=[2.44 2.83 1.82 1.46 1.47 2.27 1.61 2.65 2.19 1.52 1.07 2.42 1.83 1.05 1.4 3.39 2.22 2.42 0.68 1.15 1.72 1.78 2.08 2.35 1.52];%每个顾客需求总体积
vend_location=[76 125.47 119.47 117.47 13.47 73 -143 -62 -17 121 70 54 -23 -6 1 -118 -118 125 83 -148 121 116 -54 -136 -77;
100.47 -142.53 -138.53 -137.53 -19.53 138 -131 104 -125 74 113 53 -109 32 39 121 -40 56 -147 36 -5 -108 -16 -95 129];%客户点坐标
%上面这些也可以改，下面这些和其它文件就别碰了，一碰全是错误
center=[0;0];
coordinate=[center vend_location];
N=25;
d=zeros(N+1,N+1);
for i=1:N+1
for j=1:N+1
if ij
d(i,j)=0;
else
d(i,j)=sqrt(sum(((coordinate(:,i)-coordinate(:,j)).^2)));
end
end
end
K=length(G);
gen=0;
GGAP=1;
Popsize=1200;
lada1=0.7;
lada2=0.9;
lada3=0.005;
lada4=0.01;
FieldDR=[ones(1,N);(K)*ones(1,N)];
Chorm=inti_pop(Popsize,G,V,g,v,K)+rand(Popsize,N);
[Objv, F1, F2, F3, F4]=Objective_compute(Chorm,G,V,g,v,d,K);
trace=zeros(5,Maxgen);
[c111, index]=max(Objv);
enlist_chorm=Chorm(index,:);
h=waitbar(0,‘等去吧’);
while (gen<Maxgen)
FitV=ranking(-Objv);
Selch=select(‘rws’,Chorm,FitV,GGAP);
Selch=recombin1(Selch,lada1,lada2,FitV);
Selch=mutbga1(Selch,FieldDR,lada3,lada4,FitV);
[ObjVsel, F1sel, F2sel, F3sel, F4sel,car_route]=Objective_compute(Selch,G,V,g,v,d,K);
Ind=find(ObjVsel100000);
NN=length(Ind);
tem=inti_pop(NN,G,V,g,v,K)+rand(NN,N);
Selch(Ind,:)=tem;
[ObjVsel(Ind), F1sel(Ind), F2sel(Ind), F3sel(Ind) ,F4sel(Ind)]=Objective_compute(tem,G,V,g,v,d,K);
clear tem;
[Chorm, Objv]=reins(Chorm,Selch,1,1,Objv,ObjVsel);
[Objv1,Idx]=sort(Objv,‘ascend’);
Chorm(Idx(1),:)=enlist_chorm;
new_Chorm=inti_pop(round(Popsize0.1),G,V,g,v,K)+rand(round(Popsize0.1),N);
ObjV_new=Objective_compute(new_Chorm,G,V,g,v,d,K);
gen=gen+1;
[trace(1,gen), index]=max(Objv);
enlist_chorm=Chorm(index,:);
[a, trace(2,gen), trace(3,gen), trace(4,gen), trace(5,gen)]=Objective_compute(Chorm(index,:),G,V,g,v,d,K);
waitbar(gen/Maxgen);
end
close(h);
figure(1)
plot(trace(1,:))
xlabel(‘迭代次数’)
ylabel(‘遗传算法适应度’)
figure(2)
plot(trace(2,:))
xlabel(‘迭代次数’)
ylabel(‘运输距离’)
figure(4)
plot(trace(4,:))
xlabel(‘迭代次数’)
ylabel(‘重量利用率’)
figure(5)
plot(trace(5,:))
xlabel(‘迭代次数’)
ylabel(‘容积利用率’)
str={‘-r’,‘-g’,‘-k’,‘-b’,‘-y’,‘-m’,‘-c’,‘–b’,‘-.b’};
figure(6)
plot(vend_location(1,:),vend_location(2,:),‘k*’);
for i=1:25
text(vend_location(1,i)+1,vend_location(2,i),num2str(i),‘FontSize’,14)
end
hold on
str2={‘节点’};
kk=1;
for i=1:K
if ~isempty(car_route{i})
tem=[center vend_location(:,car_route{i}) center];
plot(tem(1,:),tem(2,:),str{kk})
kk=kk+1;
str2=[str2,strcat(‘车辆’,num2str(i))];
end
end
legend(str2)
for i=1:K
if ~isempty(car_route{i})
disp([‘派出第’,num2str(i),‘车运送的客户为：’])
disp(num2str(car_route{i}));
end
end

⛄四、运行结果

在这里插入图片描述

⛄五、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]黄戈文,蔡延光,戚远航,陈厚仁,王世豪.自适应遗传灰狼优化算法求解带容量约束的车辆路径问题[J].电子学报. 2019,47(12)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流人群疏散病毒扩散晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合