CN115953104B - 一种基于蜣螂优化算法的混合车队调度方法 - Google Patents

Abstract

本发明公开了一种基于蜣螂优化算法的混合车队调度方法包括：获取用户的订单数据；根据历史出行数据得到拟合函数，计算不同温度下的电车的电量消耗；建立目标函数，设置约束条件；使用改进的蜣螂优化算法求解混合车队的调度问题模型，得到算法迭代完成后输出的目标函数值和配送路线方案。能够提升电动汽车出行的可靠性，减少驾驶员产生里程焦虑的几率，提高电动汽车这种清洁能源交通工具路网渗透率，算法在局部搜索和全局搜索之间保持平衡，针对算法容易陷入局部最优的情况，本发明引入了具有跳跃性质的Levy飞行策略，在进行可行域搜索时，可以搜索到更全面的范围，有效保持了种群的多样性，防止算法过早收敛。

Description

一种基于蜣螂优化算法的混合车队调度方法

技术领域

本发明涉及车辆路径技术领域，具体为一种基于蜣螂优化算法的混合车队调度方法。

背景技术

近年来，城市物流业快速发展，城市配送作为物流行业重要的表现形式，贯穿于物流环节的各个阶段，逐渐成为物流行业关注的重点。但在城市物流配送发展同时，仍存在着物流运营成本高，运力资源浪费、配送服务水平低等问题。此外，配送组织模式的落后也引发了严重的城市道路拥堵和大气污染等其他问题。据国际能源署统计，交通占全世界能源消耗28％，其中交通能源40％集中在城市。“十四五”规划更是提出在“双碳”目标下，降低城市能耗、加快推动绿色低碳发展。根据上述政策可知，创建低碳环保、科学合理的配送方式已成为城市物流发展所急需解决的关键问题，其中配送中的规划方法是解决问题的关键。

当前多数物流企业面临着严重的运营成本压力，配送成本几乎占到物流总成本的一半以上，因此，配送成本的大小严重影响着综合物流成本的高低。现阶段，随着电动物流车的日渐普及，越来越多物流企业为减少配送成本而开始选择传统燃油车和电动物流车共同完成配送任务。

目前，车辆路径规划问题主要的研究对象是传统燃油车辆，对纯电动车的研究并不深入。现有电车路径问题的研究内容主要包括：考虑能量消耗和充电需求、考虑碳排放、考虑客户时间窗、求解方法研究等。电量消耗模型大多采用基于传统汽车能耗模型，包含牵引力克服四项阻力做功：加速阻力；空气阻力；爬坡阻力；滚动阻力，作为电动汽车实际消耗的电能。

但是，基于动力学原理能耗估计模型仅考虑了传统汽车行驶过程中动力学共性，实际出行中电动汽车能耗特性在模型中并未得以体现，如空调、暖气等电器消耗的电量；不同温度影响下，能耗波动导致的额外消耗等情况。这些情况导致了估计结果与实际出行电能消耗偏差。

发明内容

本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊，而这种简化或省略不能用于限制本发明的范围。

鉴于上述存在的问题，提出了本发明。

因此，本发明解决的技术问题是：现有的混合车队调度方法存在能耗波动导致的额外消耗，以及如何解决估计结果与实际出行电能消耗偏差的优化问题。

为解决上述技术问题，本发明提供如下技术方案：一种基于蜣螂优化算法的混合车队调度方法，包括：

获取用户的订单数据；

根据历史出行数据得到拟合函数，计算不同温度下的电车的电量消耗；

建立目标函数，设置约束条件；

使用改进的蜣螂优化算法求解混合车队的调度问题模型，得到算法迭代完成后输出的目标函数值和配送路线方案；

其中E_b为动力学耗电速率，公式如下：

E_t为温度影响下的额外能耗，公式如下：

E_t＝η₀+η₁T+η₂T²+η₃T³

E_c为空调使用的功耗，公式如下：

电车的电量消耗公式如下：

其中，T为当日气温；E为电量消耗系数；v是车速；m为车重；g为重力常量；ρ为空气密度；Y为迎风面积；C_d为空气动力阻力系数；C_r为滚动阻力系数；ε为车辆传动系效率；m_ijk是车辆k从i到j处的货物重量；代表使用时间，由两地距离y_ij和速度v计算而得；AC是空调的功率；η₀、η₁、η₂、η₃是公式参数，其数据值由历史出行数据拟合而得；

所述目标函数包括：

其中，x_ijk表示车辆k是否经过路段ij，y_ij为i到j的距离，P_f、P_e是燃油车和电车每公里行驶成本，ω_k为车辆k碳排放系数，C₀为每公里碳排放成本；

所述约束条件包括：

其中，A为用户下单点集合，O为配送起、终点{0，n+1}，M为充电点集合，V_f为燃油车集合，V_e为电动车集合，V为所有车辆集合，S₂为点集合A∪O，S₁为点集合A∪O∪M，W_e为电动车最大载重，W_f为燃油车最大载重，C_f燃油车最大容积，E_k为电车每公里耗电量，A_ik表示车辆k到达i点的剩余电量，L_ik表示车辆k离开i点的剩余电量，E_ek表示车辆k在充电站e的充电量，B表示电动车电池容量，W_i表示用户i的货物重量；

采用反向学习策略，解决随机生成初始种群多样性减少、种群质量不高、影响算法的收敛速度的问题：

设X是n维空间一个点，X＝|x₁,x₂...,x_n|，x_i∈[a,b]，i＝1,2,…,n，其反向点为x'_i；

根据随机策略产生最初种群，并生成其对应的反向种群；

分别计算最初种群和反向种群的个体适应度值；

比较两类种群中所对应个体的适应度值，保留适应度值较高的个体，进入最终初始种群；

其中，蜣螂优化算法中蜣螂的行为包括：

滚球：进行位置更新模型：

x_i(t+1)＝x_i(t)+α×k×x_i(t-1)+b×Δx

Δx＝|x_i(t)-X^W|

跳舞：使用正切函数进行模拟：

x_i(t+1)＝x_i(t)+tanθ|x_i(t)-x_i(t-1)|

产卵：通过边界选择策略模拟产卵区域：

Bt^*＝min(X^*×(1-R)，Bt)

Bl^*＝max(X^*×(1+R)，Bl)

边界范围随着R值的变化而变化，所以下一代的窝球位置也会变化，增强了局部搜索能力，其中下一代的位置信息计算公式为：

B_i(t+1)＝X^*+b₁×(B_i(t)-Bl^*)+b₂×(B_i(t)-Bt^*)

觅食：觅食行为的种群称为小蜣螂，其觅食区域的边界如下：

Bt′＝min(X′×(1-R)，Bt)

Bl′＝max(X′×(1+R)，Bl)

X′表示全局最佳位置，Bt′和Bl′分别代表最佳觅食区的上下限；对于小蜣螂种群，为提高全局搜索能力，引入了Levy飞行策略，由此可知，小蜣螂的位置更新公式如下：

其中，Levy(β)为服从参数β的Levy分布，0＜β＜2，μ服从N(0，σ²)分布，V服从N(0，1)分布，标准差σ表示为；

偷窃：围绕X′展开的位置变化公式如下：

x_i(t+1)＝X′+S×g×(|x_i(t)-X^*|+|x_i(t)-X′|)

其中，t代表当前迭代次数；x_i(t)代表第t次迭代时蜣螂i的位置信息；k为一个常量，表示偏转系数；b为常数，b的范围为(0，1)；α表示是否偏离方向，用1或-1表示；X^w表示全局最差位置；Δx是光照强度的变化；θ表示偏转角；|x_i(t)-x_i(t-1)|表示蜣螂i在第t次迭代时，其位置与第t-1次时的位置差；X^*表示当前局部最佳位置；Bt^*和Bl^*分别代表产卵区的上下限；Bt和Bl分别代表优化问题的上下限；T_max表示最大迭代次数；B_i(t)是窝球i第t次迭代时的位置；b₁和b₂是大小为1×D的随机向量；D是优化问题维度；x_i(t)是小蜣螂i在第t次迭代时的位置信息；C₁是[0，1]之间的随机数；表示元素相乘；Γ代表伽马分布函数；β的值为1.5；x_i(t)代表小偷i在第t次迭代时的位置信息；S是常量；g是大小为1×D的服从正态分布的随机向量；

判断各种群是否超出边界，当没有超出边界时，更新当前最优解及适应度值；

若满足终止条件，则输出最优解；若不满足终止条件，则转到对个体进行适应度评价步骤。

一种计算机设备，包括：存储器和处理器；所述存储器存储有计算机程序，其特征在于：所述处理器执行所述计算机程序时实现本发明中任一项所述的方法的步骤。

一种计算机可读存储介质，其上存储有计算机程序，其特征在于：所述计算机程序被处理器执行时实现本发明中任一项所述的方法的步骤。

本发明的有益效果：本发明提供的基于蜣螂优化算法的混合车队调度方法能够提升电动汽车出行的可靠性，减少驾驶员产生里程焦虑的几率，提高电动汽车这种清洁能源交通工具路网渗透率，算法在局部搜索和全局搜索之间保持平衡，针对算法容易陷入局部最优的情况，本发明引入了具有跳跃性质的Levy飞行策略，在进行可行域搜索时，可以搜索到更全面的范围，有效保持了种群的多样性，防止算法过早收敛。

附图说明

为了更清楚地说明本发明实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其它的附图。其中：

图1为本发明第一个实施例提供的一种基于蜣螂优化算法的混合车队调度方法的整体流程图；

图2为本发明第二个实施例提供的一种基于蜣螂优化算法的混合车队调度方法在不同温度下电车的电量消耗折线图；

图3为本发明第二个实施例提供的一种基于蜣螂优化算法的混合车队调度方法的模型运算的收敛图。

具体实施方式

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合说明书附图对本发明的具体实施方式做详细的说明，显然所描述的实施例是本发明的一部分实施例，而不是全部实施例。基于本发明中的实施例，本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明的保护的范围。

在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以采用其他不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本发明内涵的情况下做类似推广，因此本发明不受下面公开的具体实施例的限制。

其次，此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例，也不是单独的或选择性的与其他实施例互相排斥的实施例。

本发明结合示意图进行详细描述，在详述本发明实施例时，为便于说明，表示器件结构的剖面图会不依一般比例作局部放大，而且所述示意图只是示例，其在此不应限制本发明保护的范围。此外，在实际制作中应包含长度、宽度及深度的三维空间尺寸。

同时在本发明的描述中，需要说明的是，术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一、第二或第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。

本发明中除非另有明确的规定和限定，术语“安装、相连、连接”应做广义理解，例如：可以是固定连接、可拆卸连接或一体式连接；同样可以是机械连接、电连接或直接连接，也可以通过中间媒介间接相连，也可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。

实施例1

参照图1，为本发明的一个实施例，提供了一种基于蜣螂优化算法的混合车队调度方法，包括：

S1：获取用户的订单数据。

更进一步的，订单数据包括货物的重量、体积、位置信息等，及当日气温，用于判断空调是否开启；根据历史出行数据，得到拟合函数，用于计算不同温度下的电车的电量消耗。

应说明的是，假设物流企业的配送系统拥有一个配送中心和一支电动车与燃油车并存的混合车队,有n个已知位置和需求量的客户需要取货服务。要求每辆车在完成配送后都要返回起点且每辆车只能为一条路径服务,每个客户也只能被一辆车服务。由于电动车续航里程较短,可在给定的充电站候选集合中选择适合的位置充电，需进行混合车队配送路径规划,使得车队配送成本最低。

为使问题更加简洁明了，预先进行以下假设；只有一个配送中心，所有车辆从配送中心出发最后返回原配送中心；客户的数量、客户的位置及每个客户的需求量均已知，且车辆的总载重不得超过车辆的额定载重；配送车辆从配送中心出发，电池为满电状态，且在货物装卸过程中不存在电量的消耗；货物流向为单向，只进行货物收集；车辆匀速行驶，不考虑天气，限行，道路状况等；配送过程中客户的需求一次满足。

总配送成本考虑以下三部分：

1、车辆行驶费用。行驶费用与行驶距离相关，燃油和电动车辆的单位行驶成本为P_f和P_e。

2、充电费用P_ch。充电费用的多少与电池剩余电量相关，充电时的单位充电成本为商业电价每千瓦时p元，即P_ch＝p.E_ek。

3、碳排放成本。借助碳税量化碳排放成本C₀。两种车型都设置一个碳排放系数ω_k，虽然电车没有碳排放产生，但考虑到市面上火力发电依旧会产生碳排放，因此给电车也设置一个碳排放系数，但远低于燃油车的数值。每公里成本可以表示为ω_kC₀y_ij。

S2：根据历史出行数据得到拟合函数，计算不同温度下的电车的电量消耗。

应说明的是，目前由于动力电池能量密度和成本的限制，电动汽车的续驶里程较传统车仍有较大差距。而且实际使用工况复杂、车辆能耗变化大。现有车辆路径问题的研究对电量的估计精度不高，里程估计值与实际行驶距离相差很大，导致乘客担心现有的电量不能保证车辆到达目的地，产生里程焦虑，降低电动汽车的使用信心。因此，提高电动汽车续驶里程估计的精度是提高电动汽车普及率的重要因素。对企业选择电车配送也起到了关键作用。

影响电动车的电能损耗的因素有很多，包括电动车行驶的速度和距离、路况以及自身重量等。其中，车速的动态变化受到空气阻力和道路拥堵等的影响，进而影响电量的消耗；电动车的重量包括空车重量和载重，载重随着对客户点的取送而不断变化。

动力学耗电速率如下所示：

其中，E为电量消耗系数；v是车速；m为车重；g为重力常量；ρ为空气密度；Y为迎风面积；C_d为空气动力阻力系数；C_r为滚动阻力系数；ε为车辆传动系效率；m_ijk是车辆k从i到j处的货物重量。

在实际运行条件下，环境温度会影响电动汽车的续航里程。主要由于空调功率与环境温度有关，使得在其他影响因素相同情况下，电动汽车开启空调的耗电量比关闭空调的耗电量变化速率更大。且不同温度下的电池实际最大载电量发生变化，低温段最大载电量随温度下降而减小。可知，环境温度、空调运行功率影响电动汽车的耗电量。

因此，本发明在电车基于动力学的电量消耗模型，做了进一步的优化，增加了温度和空调使用的参数。对历史出行数据做函数拟合，经观察比对拟合结果，综合考虑曲线贴合程度和后续计算难度，最终进行3次函数拟合，得到二者的曲线关系为：

E_t＝η₀+η₁T+η₂T²+η₃T³

η为影响因子，根据数据可以拟合出具体数值；η₀、η₁、η₂、η₃是公式参数，其数据值由历史出行数据拟合而得。本文中其具体数值为η₀＝0.226，η₁＝7.7×10^-3，η₂＝-8.4×10^-3，η₃＝2.95×10^-6

空调使用的功耗公式：

代表使用时间，由两地距离y_ij和速度v计算而得，AC是空调的功率；

综上，电车的电量消耗公式如下：

规定目标函数为总成本最小，根据问题建立数学模型如下：

约束条件如下：

以上公式中，(1)为目标函数，包括三部分，行驶费用，碳排放费用，电车的充电费用；(2)表示每辆车从起点出发，且只执行一条路径；(3)表示每个用户点只被访问一次；(4)为流平衡约束，表示进入某个节点的车辆必须离开该节点，同时保证离开配送中心的车辆,最终回到配送中心；(5)-(8)表示燃油和电动车辆的载重和容积约束；(9)表示电车到达任意点剩余电量大于等于0；(10)电车离开配送中心或充电站应为满电；(11)电车在充电站的充电量；(12)电车装货时不损耗电量；(13)电车在行驶过程中电量损耗关系式，(14)0-1变量约束；其中，A为用户下单点集合，O为配送起、终点{0，n+1}，M为充电点集合，V_f为燃油车集合，V_e为电动车集合，V为所有车辆集合，S₂为点集合A∪O，S₁为点集合点集合A∪O∪M，W_e为电动车最大载重，W_f为燃油车最大载重，S为车辆速度，C_f燃油车最大容积，E_k为电车每公里耗电量，A_ik表示车辆k到达i点的剩余电量，L_ik表示车辆k离开i点的剩余电量，E_ek表示车辆k在充电站e的充电量，B表示电动车电池容量，W_i表示用户i的货物重量。

进行初始化操作，结合反向学习策略，优化种群，并设置算法的相关参数。初始化种群操作中，通常是随机生成初始种群，生成的种群分布不均匀，会导致种群多样性减少，种群质量不高，影响算法的收敛速度，因此采用反向学习策略来解决这一问题。

空间反向点的生成方式说明：设X是n维空间一个点，X＝|x₁,x₂...,x_n|，x_i∈[a,b]，i＝1,2,…,n，其反向点为x'_i：

生成初始种群的基本步骤如下：1、采用随机策略生成最初种群,并生成最初种群的反向种群。2、分别计算最初种群和反向种群的个体适应度值。3、将最初种群和反向种群中所对应个体的适应度值进行对比,选择适应度值更高的个体进入最终初始种群。

应该知道的是，这样既增大了种群的多样性又保留了优秀个体，避免算法陷入局部最优的同时加快了算法的运算速度。种群数为N，其中分为四类种群，分别为滚球蜣螂、窝球、小蜣螂和小偷，四类种群的位置更新方式不同。

S3：使用改进的蜣螂优化算法求解混合车队的调度问题模型，得到算法迭代完成后输出的目标函数值和配送路线方案。

更进一步的，对个体进行适应度评价后，判断各种群是否超出边界；更新当前最优解及其适应度值；最后判断是否满足终止条件，如若满足，就停止迭代，输出最优解，否则转到对个体进行适应度评价步骤。

适应度值应该反映实际对应路径的所有信息(包括油耗、总配送距离，约束处理情况)，所以本方法采用算术相加的方法来定义适应度值，计算公式如下：

式中，z为式(1)的目标函数，f(i)即某个种群所对应的适应度值，f(i)越大，被选中的几率就越大；

应该知道的是，更新蜣螂位置，通过模拟蜣螂的行为，一共有五种行为方式，如下：

方式1滚球，蜣螂在滚动过程中通过外界信息来导航，以保持粪球在直线上滚动。为了模拟滚动球的行为，需要蜣螂在整个搜索空间中沿着给定的方向移动。按下列公式进行位置更新：

x_i(t+1)＝x_i(t)+α×k×x_i(t-1)+b×Δx (17)

Δx＝|x_i(t)-X^w| (18)

其中，t代表当前迭代次数，x_i(t)代表第t次迭代蜣螂i的位置信息，k偏转系数是一个常量，b也是常数范围在(0，1)，α是自然系数值为1或-1，X^w表示全局最差位置，Δx是光照强度的变化。

α的存在相当于变异操作，可以扩大问题空间的探索范围，使得局部搜索能力更好，降低陷入局部最优的可能性。通过概率法模拟是否遇到障碍物，设置一个值λ∈[0，1]，依此为标准，随机生成[0，1]之间的数∈，α的变化过程如下：

当α＝1时，表示无偏差。α＝-1时，会偏离方向。

方式2跳舞，蜣螂遇到障碍物无法前进时，需要通过跳舞行为来改变方向，以获得新的路线。使用正切函数进行模拟，如下所示：

x_i(t+1)＝x_i(t)+tanθ|x_i(t)-x_i(t-1)| (19)

θ表示偏转角，|x_i(t)-x_i(t-1)|表示蜣螂i在第t次迭代时，其位置与第t-1次时的位置差，以起到参考历史信息的作用。

当θ值为0，π/2，π时，正切函数无值，此时表示不更新位置。

滚球蜣螂基于方式1和方式2进行位置更新，通过随机数的方式，判断使用何种行为方式进行位置更新，生成随机数δ∈[0，1]

方式3产卵，蜣螂会滚动粪球到安全位置，繁衍后代，每次产卵产生一个后代。通过边界选择策略模拟产卵区域，公式如下：

Bt^*＝min(X^*×(1-R)，Bt) (20)

Bl^*＝max(X^*×(1+R)，B1) (21)

X^*表示当前局部最佳位置，Bt^*和Bl^*分别代表产卵区的上下限，Bt和Bl分别代表优化问题上下限。T_max表示最大迭代次数。

边界范围随着R值的变化而变化，所以下一代的窝球位置也会变化，增强了局部搜索能力，下一代的位置信息计算公式如下：

B_i(t+1)＝X^*+b₁×(B_i(t)-Bl^*)+b₂×(B_i(t)-Bt^*) (23)

式中，B_i(t)是窝球i第t次迭代的位置，b₁和b₂是大小为1×D的随机向量，D是优化问题维度。窝球的位置需要限制在产卵区，如有超过，通过以下方式进行转换：

方式4觅食，觅食行为的种群称为小蜣螂，其觅食区域的边界如下：

Bt′＝min(X′×(1-R)，Bt) (24)

Bl′＝max(X′×(1+R)，Bl) (25)

X′表示全局最佳位置，Bt′和Bl′分别代表最佳觅食区的上下限。对于小蜣螂种群，为提高全局搜索能力，引入了Levy飞行策略，由此可知，小蜣螂的位置更新公式如下：

式中，x_i(t)是小蜣螂i在第t次迭代的位置信息，C₁是[0，1]之间的随机数，是元素相乘。

Levy飞行是一种随机行走策略，多数情况下，模拟的步长较小，偶尔也会有较大步长的情况，可以帮助跳出局部最优的情况，其公式如下：

式中，Levy(β)为服从参数β的Levy分布，0＜β＜2，μ服从N(0，σ²)分布，V服从N(0，1)分布，σ可以下列公式计算得到；

其中，Γ代表伽马分布函数，β的值为1.5。

方式5偷窃，由方式4可知X′是全局最优值，所以模拟蜣螂的偷窃行为，相当于收敛的过程。此类是算法的重点，因为最终达到模型收敛的过程。围绕X′展开的位置变化公式如下：

x_i(t+1)＝X′+S×g×(|x_i(t)-X^*|+|x_i(t)-X′|) (29)

式中，x_i(t)代表小偷i在第t次迭代的位置信息，S是常量，g是大小为1×D的服从正态分布的随机向量。

更进一步的，判断各种群是否超出边界，当没有超出边界时，更新当前最优解及适应度值；

若满足终止条件，则输出最优解；若不满足终止条件，则转到对个体进行适应度评价步骤本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器、磁变存储器、铁电存储器、相变存储器、石墨烯存储器等。易失性存储器可包括随机存取存储器或外部高速缓冲存储器等。作为说明而非局限，RAM可以是多种形式，比如静态随机存取存储器或动态随机存取存储器等。本申请所提供的各实施例中所涉及的数据库可包括关系型数据库和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数据库等，不限于此。

本申请所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等，不限于此。

实施例2

参照图2-3，为本发明的一个实施例，提供了一种基于蜣螂优化算法的混合车队调度方法，为了验证本发明的有益效果，通过经济效益计算和仿真实验进行科学论证。

在实际运行条件下，环境温度会影响电动汽车的续航里程。主要由于空调功率与环境温度有关，使得在其他影响因素相同情况下，电动汽车开启空调的耗电量比关闭空调的耗电量变化速率更大。且不同温度下的电池实际最大载电量发生变化，低温段最大载电量随温度下降而减小。可知，环境温度、空调运行功率影响电动汽车的耗电量，具体变化情况如图2所示。目标函数模型运算的收敛如图3所示。

选取客户数为30的测试实例，客户需求量、客户坐标和配送中心坐标序号1，如表1所示。

表1

通过以下仿真实验进一步说明本发明的效果:

1.实验条件：在Intel(R)Core(TM)i5-7300HQ CPU@2.50GHz内存16GB、WINDOWS 10系统上使用MATLAB 2016a进行仿真。

2.实验内容：选取客户数为30的测试实例，具体数据值，如表1所示。

3.实验结果：对模型优化仿真30次，优化结果取平均值。使用经典的遗传算法和粒子群算法进行对比，结果如表2所示。

表2

算法名称	总成本/元	行驶距离/km	车辆数量	计算耗时/s
					蜣螂优化算法	1071.99	152.69	7.0	39.6
遗传算法	1098.31	151.68	7.2	51.4
					粒子群算法	1078.52	151.71	7.1	42.4

综合计算结果表明，蜣螂优化算法对模型的求解效果较好。

应说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (3)

1.一种基于蜣螂优化算法的混合车队调度方法，其特征在于，包括：

获取用户的订单数据；

根据历史出行数据得到拟合函数，计算不同温度下的电车的电量消耗；

建立目标函数，设置约束条件；

使用改进的蜣螂优化算法求解混合车队的调度问题模型，得到算法迭代完成后输出的目标函数值和配送路线方案；

电车的电量消耗公式如下：

其中，E_b为动力学耗电速率，公式如下：

E_t为温度影响下的额外能耗，公式如下：

E_t＝η₀+η₁T+η₂T²+η₃T³

E_c为空调使用的功耗，公式如下：

其中，T为当日气温；E为电量消耗系数；v是车速；m为车重；g为重力常量；ρ为空气密度；Y为迎风面积；C_d为空气动力阻力系数；C_r为滚动阻力系数；ε为车辆传动系效率；m_ijk是车辆k从i到j处的货物重量；代表使用时间，由两地距离y_ij和速度v计算而得；AC是空调的功率；η₀、η₁、η₂、η₃是公式参数，其数据值由历史出行数据拟合而得；

所述目标函数包括：

其中，x_ijk表示车辆k是否经过路段ij，y_ij为i到j的距离，P_f、P_e是燃油车和电车每公里行驶成本，ω_k为车辆k碳排放系数，C₀为每公里碳排放成本；

所述约束条件包括：

其中，A为用户下单点集合，O为配送起、终点{0，n+1}，M为充电点集合，V_f为燃油车集合，V_e为电动车集合，V为所有车辆集合，S₂代表点集合A∪O，S₁代表点集合A∪O∪M，W_e为电动车最大载重，W_f为燃油车最大载重，C_f为燃油车最大容积，E_k为电车每公里耗电量，A_ik表示车辆k到达i点的剩余电量，L_ik表示车辆k离开i点的剩余电量，E_ek表示车辆k在充电站e的充电量，B表示电动车电池容量，W_i表示用户i的货物重量；

所述混合车队调度方法采用反向学习策略，解决由于随机生成初始种群多样性减少、种群质量不高、影响算法的收敛速度的问题：

设X是n维空间一个点，X＝|x₁,x₂...,x_n|，x_i∈[a,b]，i＝1,2,…,n，其反向点为x'_i；

根据随机策略产生最初种群，并生成其对应的反向种群；

分别计算最初种群和反向种群的个体适应度值；

比较两类种群中所对应个体的适应度值，保留适应度值较高的个体，进入最终初始种群；

所述蜣螂优化算法中蜣螂的行为包括：

滚球：进行位置更新模型：

x_i(t+1)＝x_i(t)+α×k×x_i(t-1)+b×Δx

Δx＝|x_i(t)-X^w|

跳舞：使用正切函数进行模拟：

x_i(t+1)＝x_i(t)+tanθ|x_i(t)-x_i(t-1)|

产卵：通过边界选择策略模拟产卵区域：

Bt^*＝min(X^*×(1-R),Bt)

Bl^*＝max(X^*×(1+R),Bl)

边界范围随着R值的变化而变化，所以下一代的窝球位置也会变化，增强了局部搜索能力，其中下一代的位置信息计算公式为：

B_i(t+1)＝X^*+b₁×(B_i(t)-Bl^*)+b₂×(B_i(t)-Bt^*)

觅食：觅食行为的种群称为小蜣螂，其觅食区域的边界如下：

Bt'＝min(X'×(1-R),Bt)

Bl'＝max(X'×(1+R),Bl)

X'表示全局最佳位置，Bt'和Bl'分别代表最佳觅食区的上下限；对于小蜣螂种群，为提高全局搜索能力，引入了Levy飞行策略，由此可知，小蜣螂的位置更新公式如下：

其中，Levy(β)为服从参数β的Levy分布，0<β<2，μ服从N(0,σ²)分布，v服从N(0,1)分布，标准差σ表示为；

偷窃：围绕X'展开的位置变化公式如下：

x_i(t+1)＝X'+S×g×(|x_i(t)-X^*|+|x_i(t)-X'|)

其中，t代表当前迭代次数；x_i(t)代表第t次迭代时蜣螂i的位置信息；k为一个常量，表示偏转系数；b为常数，b的范围为(0,1)；α表示是否偏离方向，用1或-1表示；X^w表示全局最差位置；Δx是光照强度的变化；θ表示偏转角；|x_i(t)-x_i(t-1)|表示蜣螂i在第t次迭代时，其位置与第t-1次时的位置差；X^*表示当前局部最佳位置；Bt^*和Bl^*分别代表产卵区的上下限；Bt和Bl分别代表优化问题的上下限；T_max表示最大迭代次数；B_i(t)是窝球i第t次迭代时的位置；b₁和b₂是大小为1×D的随机向量；D是优化问题维度；x_i(t)是小蜣螂i在第t次迭代时的位置信息；C₁是[0,1]之间的随机数；⊕表示元素相乘；Γ代表伽马分布函数；β的值为1.5；x_i(t)代表小偷i在第t次迭代时的位置信息；S是常量；g是大小为1×D的服从正态分布的随机向量；

判断各种群是否超出边界，当没有超出边界时，更新当前最优解及适应度值；

若满足终止条件，则输出最优解；若不满足终止条件，则转到对个体进行适应度评价步骤。

2.一种计算机设备，包括：存储器和处理器；所述存储器存储有计算机程序，其特征在于：所述处理器执行所述计算机程序时实现权利要求1中所述的方法的步骤。

3.一种计算机可读存储介质，其上存储有计算机程序，其特征在于：所述计算机程序被处理器执行时实现权利要求1中所述的方法的步骤。