【优化调度】基于matlab遗传算法求解公交车调度优化问题【含Matlab源码 040期】

✅博主简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,Matlab项目合作可私信。
🍎个人主页:海神之光
🏆代码获取方式:
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭:行百里者,半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理(进阶版)
路径规划(Matlab)
神经网络预测与分类(Matlab)
优化求解(Matlab)
语音处理(Matlab)
信号处理(Matlab)
车间调度(Matlab)

⛄一、遗传算法简介

1 引言
在这里插入图片描述
在这里插入图片描述
2 遗传算法理论
2.1 遗传算法的生物学基础
在这里插入图片描述
在这里插入图片描述
2.2 遗传算法的理论基础
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
2.3 遗传算法的基本概念
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
2.4 标准的遗传算法
在这里插入图片描述
在这里插入图片描述
2.5 遗传算法的特点
在这里插入图片描述
在这里插入图片描述
2.6 遗传算法的改进方向
在这里插入图片描述
3 遗传算法流程
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
4 关键参数说明
在这里插入图片描述

⛄二、部分源代码

clear all;
close all;
clc;

%变量初始化
populationSize=80; %群体大小
generationSize=800; %进化代数
codeLength=35; %编码长度
children=zeros(populationSize,7); %存储子代分量及适应度
pr=0.15; %复制概率
pm=0.1; %变异概率
pc=0.60; %交叉概率
E=round(rand(populationSize,codeLength)); %随机产生二进制串

%主函数
for k=1:1:generationSize

for s=1:1:populationSize
x=E(s,:);
temp=zeros(1,7);%临时保存子代各分量的十进制值
%编码
for i=1:35
if i<=6
temp(1)=temp(1)+x(i)2^(i-1);
elseif i>=7 && i<=13
temp(2)=temp(2)+x(i)2^(i-6-1);
elseif i>=14 && i<=19
temp(3)=temp(3)+x(i)2^(i-13-1);
elseif i>=20 && i<=25
temp(4)=temp(4)+x(i)2^(i-19-1);
elseif i>=26 && i<=30
temp(5)=temp(5)+x(i)2^(i-25-1);
elseif i>=31 && i<=35
temp(6)=temp(6)+x(i)2^(i-30-1);
end
end
children(s,1)=round(60
temp(1)/(2^6-1));
children(s,2)=round(70
temp(2)/(2^7-1));
children(s,3)=round(60
temp(3)/(2^6-1));
children(s,4)=round(50
temp(4)/(2^6-1));
children(s,5)=round(20
temp(5)/(2^5-1));
children(s,6)=round(30
temp(6)/(2^5-1));
children(s,7)=round(children(s,1)+children(s,2)+children(s,3)+children(s,4)+children(s,5)+children(s,6));

while children(s,1)+children(s,6)<60||children(s,1)+children(s,2)<70||children(s,2)+children(s,3)<60||…
children(s,3)+children(s,4)<50||children(s,4)+children(s,5)<20||children(s,5)+children(s,6)<30
E(s,:)=round(rand(1,codeLength));
x=E(s,:);
temp=zeros(1,7);
for i=1:35
if i<=6
temp(1)=temp(1)+x(i)2^(i-1);
elseif i>=7 && i<=13
temp(2)=temp(2)+x(i)2^(i-6-1);
elseif i>=14 && i<=19
temp(3)=temp(3)+x(i)2^(i-13-1);
elseif i>=20 && i<=25
temp(4)=temp(4)+x(i)2^(i-19-1);
elseif i>=26 && i<=30
temp(5)=temp(5)+x(i)2^(i-25-1);
elseif i>=31 && i<=35
temp(6)=temp(6)+x(i)2^(i-30-1);
end
end
children(s,1)=round(60
temp(1)/(2^6-1));
children(s,2)=round(70
temp(2)/(2^7-1));
children(s,3)=round(60
temp(3)/(2^6-1));
children(s,4)=round(50
temp(4)/(2^6-1));
children(s,5)=round(20
temp(5)/(2^5-1));
children(s,6)=round(30
temp(6)/(2^5-1));
children(s,7)=round(children(s,1)+children(s,2)+children(s,3)+children(s,4)+children(s,5)+children(s,6));
end
end

% Step 1 : 确定适应度函数
fitness=1./children(:,7);
[orderFitness,indexFitness]=sort(fitness); %orderFitness按行从小到大排序的新矩阵,indexFitness元素位置构成的新矩阵
bestFitness=orderFitness(populationSize); %bestFitness=max(fitness)
bestS=E(indexFitness(populationSize)😅; %bestS 精英保留

%Step 2 : 复制操作
sumFitness=sum(fitness);
fitnessSize=(orderFitness/sumFitness)*populationSize;
fitnessSelect=floor(fitnessSize); %朝负无穷大方向取整
kk=1;
for i=1:1:populationSize
for j=1:1:fitnessSelect(i) %Select and Reproduce
TempE(kk,:)=E(indexFitness(i)😅;
kk=kk+1; %kk is used to reproduce
end
end

%Step 3 : 交叉操作
n=ceil(35*rand); %朝正无穷大方向取整
for i=1:2:(populationSize-1)
temp=rand;
if pc>temp %Crossover Condition
for j=n:1:35
TempE(i,j)=E(i+1,j);
TempE(i+1,j)=E(i,j);
end
end
end
TempE(populationSize,:)=bestS;
E=TempE;

%Step 4: 变异操作
for i=1:1:populationSize
%Mutation Condition
if TempE(i,j)==0
TempE(i,j)=1;
else
TempE(i,j)=0;
end
end
end
end
TempE(populationSize,:)=bestS;
E=TempE;
disp([‘第’,num2str(k),‘代:’,num2str(children(populationSize,:))]);
end
disp([‘二进制编码:’,num2str(bestS)]);
disp([‘最优选择:’,num2str(children(populationSize,:))]);

⛄三、运行结果

在这里插入图片描述

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 包子阳,余继周,杨杉.智能优化算法及其MATLAB实例(第2版)[M].电子工业出版社,2016.
[2]张岩,吴水根.MATLAB优化算法源代码[M].清华大学出版社,2017.

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 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合

  • 6
    点赞
  • 36
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 2
    评论
摘 要 公共交通是城市交通的重要组成部分,作好公交车调度对于完善城市交通环境、改进市民出行状况、提高公交公司的经济和社会效益,都具有重要意义。 本题要求我们得出一个最优的公交车调度方案,通过分析题中所给数据与要求,得出此题为多目标线性规划问题。由于涉及到多目标,为了便于求解要将多目标转化为单目标规划模型,由于题中所给目标具有抽象性,因此需要先根据所给资料对目标进行定量化分析,即选择合适的参数表示目标函数,本文采用等车时间与乘车时拥挤程度表示乘客的利益,平均载客率表示公交公司的利益,显然这两者的关系是互斥的,即公交车数量越少,虽公交公司利益最大化,但等车时间越长,拥挤程度更大,造成乘客抱怨度更高。由此,本文建立了两个模型,关于模型Ⅰ,将公交公司的利益作为目标函数,为了照顾乘客的满意度,考虑在早高峰时间将等车时间尽量缩短至少于5分钟,其他时间段都控制在5到10分钟之内,建立单目标的最小车次线性规划模型并运用LINGO软件求得公交车调度车次。最终得到最优方案中需要62辆公交车,具体调度方案见附录中表5-1。关于模型Ⅱ,我们主要利用模糊数学及层次分析法设置加权因子将其转化为单目标优化问题,并判断在模型Ⅰ的情况下该放宽还是加紧约束条件才可得到最优解,在其他时段则优先考虑公交车公司的满意度,并用MATLAB拟合公交公司与乘客的满意度。最终得到最优方案中需要64辆公交车,具体调度方案见附录中表5-2。 最后,我们对两个模型进行了误差分析与稳定性分析,通过误差分析得到本文模型仍存在误差,误差与提供的数据本身以及模型的假设均有关。因此为了减少误差,使结果更加准确,我们需要更加准确的各时间段客流量数据,且时间间隔越短可以使结果更加精确,更有实用价值。通过稳定性分析发现在对最大满载率及乘客在一般时间内的等待时间做微小变动时,模型Ⅰ、Ⅱ 结果都没有太大变化,由此可知两个模型的稳定性都比较好。 本文中的模型通俗易懂,易于用软件得到确切结果,具有实际意义,但考虑因素还不够全面,应当结合实际多方面的情况,设计出更加优化调度方案,模型较稳定,具有推广意义。
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

海神之光

有机会获得赠送范围1份代码

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值