CN108805335A - 一种公共自行车调度方法 - Google Patents
一种公共自行车调度方法 Download PDFInfo
- Publication number
- CN108805335A CN108805335A CN201810475790.7A CN201810475790A CN108805335A CN 108805335 A CN108805335 A CN 108805335A CN 201810475790 A CN201810475790 A CN 201810475790A CN 108805335 A CN108805335 A CN 108805335A
- Authority
- CN
- China
- Prior art keywords
- website
- public bicycles
- indicate
- soil
- scheduling
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000002689 soil Substances 0.000 claims abstract description 69
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 56
- 238000002922 simulated annealing Methods 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 210000003813 thumb Anatomy 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000012360 testing method Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 210000004247 hand Anatomy 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/40—Business processes related to the transportation industry
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开一种公共自行车调度方法,包括,获取子区域内各公共自行车站点的位置信息以及需求量信息;根据所述公共自行车站点的位置信息以及需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度,得到最优路径;获取本次迭代计算得到的最优解,更新所述最优解所经过的路径上的泥土量;判断所述迭代次数是否大于预先设定值,得到第二判断结果,若第二判断结果表示所述迭代次数大于所述预先设定值,则执行解的优化;直到迭代次数达到预先设定值,当迭代次数达到预先设定值后输出所有迭代中最优的调度方案,进而实现对自行车的精确调度。
Description
技术领域
本发明涉及城市智能公共交通系统技术领域,特别是涉及一种公共自行车调度方法。
背景技术
随着城市的快速发展,城市人口不断增加,机动车的数量也随着大量增长,使得交通拥堵和环境污染问题日益严峻。充分发挥城市公共自行车的作用,能够有效地缓解这些问题。但是目前公共自行车在运营过程中出现了一些问题影响着运营效率,“租车难、换车难”就是在自行车使用过程中用户反馈最强烈的问题,即某些自行车租赁站点在一些时间段内的自行车数量不够多,使得用户不能租到自行车。某些自行车站点的在一些时间段内没有停车位,导致用户不能归还自行车。而解决这一问题的关键就是对公共自行车进行智能调度,合理的自行车调度可以提高用户对公共自行车使用的满意度,并且节约调度成本,提高运营效率,对促进公民绿色出行、缓解交通拥堵具有重大的意义。
目前,国内外在此领域的研究并不多,但也取得了一定的研究成果,然而现有技术中的的求解效率不够高效,对问题建模不够合理,对调度策略的设计不够完善,导致调度方法的整体求解效率不佳。
发明内容
本发明的目的是提供一种公共自行车调度方法,来实现对公共自行车的调度,解决生活中公共自行车站点的自行车数量过饱和以及自行车数量过少不够用的问题。
为实现上述目的,本发明提供了如下方案:
一种公共自行车调度方法,所述方法包括:
获取子区域内各公共自行车站点的位置信息以及需求量信息;所述子区域为对各所述公共自行车站点进行划分后形成的子区域;每个所述子区域包含一个调度中心;
根据所述公共自行车站点的位置信息以及需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度,得到最优路径,具体步骤如下:
步骤1、将智能水滴算法中的参数进行初始化,并初始化迭代次数a=1;
步骤2、构造调度方案,具体包括:
步骤2-1、水滴根据概率矩阵,优先选择概率较大的公共自行车站点作为要服务站点,所述概率矩阵公式如下:
其中,τij为启发式算子,dij表示从站点i到站点j之间的距离;ε为大于零的自然数,δ为大零的自然数,FV为需要服务的站点集合,uij为节约算子,表示路径的节约量,uij=di0+d0j-dij,0∈depot,i,j∈(1,2,…N),di0表示站点i到所述调度中心的距离,d0j表示站点j到所述调度中心的距离;其中r=0.01,sw(i,j)为从站点i到站点j的路径arc(i,j)中的泥土量;
判断所述服务站点是否满足约束条件,得到第一判断结果,当所述第一判断结果表示满足约束条件时,则对所述服务站点进行调度,当所述第一判断结果表示不满足约束条件时,则按照水滴到其他所述公共自行车站点的概率值从大到小对站点进行访问,直到找到满足所述约束条件的站点;
所述约束条件为:
0≤dj+wijk≤Q,(i,j∈N∪M),k∈(1,2,…,K)
ej≤ti+wi+tij-k(1-xijk)≤fj,
其中,Q表示每量运输车的最大承载量,dj表示运输车到达公共自行车站点的j距离,tij表示运输车从公共自行车站点i到公共自行车站点j的时间,ti表示运输车到达站点i的时间,wijk表示运输车从公共自行车站点i到公共自行车站点j的自行车数量,wi表示运输车在站点i的停留时间;M表示调度中心的集合,M={SN+1,SN+2,…,SN+M},N表示公共自行车站点的集合,N={v1,v2,…vN};ei、fj表示调度时刻;
如果所有可以访问的公共自行车站点都不满足所述约束条件,则水滴返回调度中心,由调度中心派遣其他运输车辆,当被访问的公共自行车站点满足所述约束条件时,确定当前被访问的公共自行车站点为被服务的下一个自行车站点;
步骤2-2.更新水滴的流速vely、水滴中的泥土量soiliwd y、路径arc(i,j)中的泥土量sw(i,j),具体公式如下:
其中vely为当前水滴的流速,av、bv、cv为速度更新参数,vely-1为上一次迭代中水滴的流速;
soiliwd y=soiliwd y-1+Δsoil(i,j),其中soiliwd y为当前水滴中的泥土量,soiliwd y-1为上一次迭代水滴中的泥土量, as、bs、cs为泥土更新参数,D(i,j)表示从站点i到站点j的路程,ε为大于零的自然数;
sw(i,j)y=sw(i,j)y-1-α*Δsoil(i,j),soilmin≤sw(i,j)y≤soilmax;其中,sw(i,j)y为当前路径中的泥土量,sw(i,j)y-1为上一次迭代路径中的泥土量,ρ为大于零的自然数;TTB表示全局最优解;α的为大于零的自然数;as、bs、cs为泥土更新参数,n为水滴到达目的地的个数;
重复所述步骤2-1和所述步骤2-2直到所有需要调度服务的站点都被服务;
步骤3、获取本次迭代计算得到的最优解,迭代次数a=a+1;
步骤4、更新所述最优解所经过的路径上的泥土量;
步骤5、判断所述迭代次数是否大于预先设定值,得到第二判断结果,若第二判断结果表示所述迭代次数大于所述预先设定值,则执行解的优化;
步骤6、重复所述步骤2-所述步骤5直到迭代次数达到预先设定值,当迭代次数达到预先设定值后输出所有迭代中最优的调度方案。
可选的,所述对各所述公共自行车站点进行划分具体包括:
按照最邻近原则将距离每个调度中心前K个最近的公共自行车站点划分为一个子区域,具体公式如下:
其中N表示公共自行车站点的数量,M表示子区域的数量。
可选的,所述公共自行车站点的位置信息为坐标信息,所述需求量信息为各个公共自行车站点的调度量di、所述调度量di对应的时域[ei,fi]、所述调度量di对应的最佳时域[ai,bi]。
可选的,所述方法还包括:在所述根据所述公共自行车站点的位置信息以及需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度,得到最优路径之后判断对所述公共自行车的调度是否满足调度条件,得到第三判断结果,当第三判断结果表示对所述公共自行车的调度满足所述调度条件时,则根据所述公共自行车站点的位置信息和需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度:
当在最佳调度时域[ai,bi],di≥0时,通过所述调度中心调入di辆自行车至所述公共自行车站点i,当在最佳调度时域[ai,bi],当di<0时,通过所述调度中心从所述公共自行车站点i调出di辆自行车;
当所述第三判断结果表示对所述公共自行车的调度不满足所述调度条件时,则对不满足调度条件的情况进行惩罚:
当对公共自行车的调度时刻在[ai,bi]外且在[ei,fi]内时,此时为不满足调度条件的情况,对所述不满足调度条件的情况进行惩罚,惩罚函数如下:
其中Cw表示每秒给工作人员的薪资,表示单位里程的油价,ai、bi表示调度时刻,ti表示对站点i的服务时间;
可选的,所述执行解的优化具体包括:
步骤4-1、获取每次迭代中所有最优解,将所述所有最优解中代价值最小的最优解作为全局最优解;所述代价值用cost(Xnew)表示,
其中Cw表示每秒给工作人员的薪资,Cb表示单位里程的油价,ai、bi表示调度时刻,dij表示站点i到站点j的距离;ti表示对站点i的服务时间;
步骤4-2、将所述全局最优解作为模拟退火算法的初始解,所述初始解用X表示,初始化模拟退火算法中的所有参数;
步骤4-3、对所述初始解执行3-opt变换,得到新解,用Xnew表示;
步骤4-4、计算所述Xnew的代价值,所述代价值用cost(Xnew)表示,其中Cw表示每秒给工作人员的薪资,Cb表示单位里程的油价,ai、bi表示调度时刻,dij表示站点i到站点j的距离;ti表示对站点i的服务时间;
步骤4-4、按照Metropolis准则接收新解。
重复步骤(2)-(4)直到迭代次数满足预先设定值。
可选的,所述更新所述最优解所经过的路径上的泥土量具体包括:
其中TIB表示最优解,sw(i,j)b-1为上一次迭代路径中的泥土量,表示最优路径上的水滴携带的泥土量,NIB表示当前迭代得到的最优解中的站点的个数。
可选的,所述更新全局最优解,具体包括:
当cost(TIB)<cost(TTB)时,TTB=TIB,其中TIB表示最优解,TTB表示全局最优解。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明中的该方法能够根据各个站点的自行车需求信息,制定合理的调度方案,能够满足站点的需求。
本发明将智能水滴算法和模拟退火算法融合一起,提出了新的混合智能水滴算法,改善了原算法容易过早收敛的缺点,提高了算法的全局搜索能力以及对目标函数优化能力。并且,引入启发式算子和节约算子,提高了整体方法的运行速率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例公共自行车调度方法流程图;
图2为本发明实施例对公共自行车站点划分的区域图;
图3为本发明实施例惩罚函数图;
图4为本发明实施例表格2的算例数据求解结果图;
图5为本发明实施例每代路径最小代价对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种公共自行车调度方法,来实现对公共自行车的调度,解决生活中公共自行车站点的自行车数量过饱和以及自行车数量过少不够用的问题。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例公共自行车调度方法流程图,如图1所示,所述方法包括:
步骤101:获取子区域内各公共自行车站点的位置信息以及需求量信息;所述子区域为对各所述公共自行车站点进行划分后形成的子区域;每个所述子区域包含一个调度中心;
步骤102:根据所述公共自行车站点的位置信息以及需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度,得到最优路径,具体步骤如下:
步骤1011、将智能水滴算法中的参数进行初始化,并初始化迭代次数a=1;
步骤1012、构造调度方案,具体包括:
步骤10121、水滴根据概率矩阵,优先选择概率较大的公共自行车站点作为要服务站点,所述概率矩阵公式如下:
其中,τij为启发式算子,dij表示从站点i到站点j之间的距离;ε为大于零的自然数,δ为大零的自然数,FV为需要服务的站点集合,uij为节约算子,表示路径的节约量,uij=di0+d0j-dij,0∈depot,i,j∈(1,2,…N),di0表示站点i到所述调度中心的距离,d0j表示站点j到所述调度中心的距离;其中r=0.01,sw(i,j)为从站点i到站点j的路径arc(i,j)中的泥土量;
判断所述服务站点是否满足约束条件,得到第一判断结果,当所述第一判断结果表示满足约束条件时,则对所述服务站点进行调度,当所述第一判断结果表示不满足约束条件时,则按照水滴到其他所述公共自行车站点的概率值从大到小对站点进行访问,直到找到满足所述约束条件的站点;
所述约束条件为:
0≤dj+wijk≤Q,(i,j∈N∪M),k∈(1,2,…,K)
ej≤ti+wi+tij-k(1-xijk)≤fj,
其中,Q表示每量运输车的最大承载量,dj表示运输车到达公共自行车站点的j距离,tij表示运输车从公共自行车站点i到公共自行车站点j的时间,ti表示运输车到达站点i的时间,wijk表示运输车从公共自行车站点i到公共自行车站点j的自行车数量,wi表示运输车在站点i的停留时间;M表示调度中心的集合,M={SN+1,SN+2,…,SN+M},N表示公共自行车站点的集合,N={v1,v2,…vN};ei、fj表示调度时刻;
如果所有可以访问的公共自行车站点都不满足所述约束条件,则水滴返回调度中心,由调度中心派遣其他运输车辆,当被访问的公共自行车站点满足所述约束条件时,确定当前被访问的公共自行车站点为被服务的下一个自行车站点;
步骤10122.更新水滴的流速vely、水滴中的泥土量soiliwd y、路径arc(i,j)中的泥土量sw(i,j),具体公式如下:
其中vely为当前水滴的流速,av、bv、cv为速度更新参数,vely-1为上一次迭代中水滴的流速;
soiliwd y=soiliwd y-1+Δsoil(i,j),其中soiliwd y为当前水滴中的泥土量,soiliwd y-1为上一次迭代水滴中的泥土量, as、bs、cs为泥土更新参数,D(i,j)表示从站点i到站点j的路程,ε为大于零的自然数;
sw(i,j)y=sw(i,j)y-1-α*Δsoil(i,j),soilmin≤sw(i,j)y≤soilmax;其中,sw(i,j)y为当前路径中的泥土量,sw(i,j)y-1为上一次迭代路径中的泥土量,ρ为大于零的自然数;TTB表示全局最优解;α的为大于零的自然数;as、bs、cs为泥土更新参数,n为水滴到达目的地的个数;
重复所述步骤10121和所述步骤10122直到所有需要调度服务的站点都被服务;
步骤103、获取本次迭代计算得到的最优解,迭代次数a=a+1;
步骤104、更新所述最优解所经过的路径上的泥土量;
步骤105、判断所述迭代次数是否大于预先设定值,得到第二判断结果,若第二判断结果表示所述迭代次数大于所述预先设定值,则执行解的优化;
步骤106、重复所述步骤102至所述步骤105直到迭代次数达到预先设定值,当迭代次数达到预先设定值后输出所有迭代中最优的调度方案。
具体的,所述方法还包括:在所述根据所述公共自行车站点的位置信息以及需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度,得到最优路径之后判断对所述公共自行车的调度是否满足调度条件,得到第三判断结果,当第三判断结果表示对所述公共自行车的调度满足所述调度条件时,则根据所述公共自行车站点的位置信息和需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度:
当在最佳调度时域[ai,bi],di≥0时,通过所述调度中心调入di辆自行车至所述公共自行车站点i,当在最佳调度时域[ai,bi],当di<0时,通过所述调度中心从所述公共自行车站点i调出di辆自行车;
当所述第三判断结果表示对所述公共自行车的调度不满足所述调度条件时,则对不满足调度条件的情况进行惩罚:
当对公共自行车的调度时刻在[ai,bi]外且在[ei,fi]内时,此时为不满足调度条件的情况,对所述不满足调度条件的情况进行惩罚,惩罚函数如下:
其中Cw表示每秒给工作人员的薪资,表示单位里程的油价,ai、bi表示调度时刻,ti表示对站点i的服务时间;
具体的,所述执行解的优化具体包括:
步骤4-1、获取每次迭代中所有最优解,将所述所有最优解中代价值最小的最优解作为全局最优解;所述代价值用cost(Xnew)表示,
其中Cw表示每秒给工作人员的薪资,Cb表示单位里程的油价,ai、bi表示调度时刻,dij表示站点i到站点j的距离;ti表示对站点i的服务时间;
步骤4-2、将所述全局最优解作为模拟退火算法的初始解,所述初始解用X表示,初始化模拟退火算法中的所有参数;
步骤4-3、对所述初始解执行3-opt变换,得到新解,用Xnew表示;
步骤4-4、计算所述Xnew的代价值,所述代价值用cost(Xnew)表示,其中Cw表示每秒给工作人员的薪资,Cb表示单位里程的油价,ai、bi表示调度时刻,dij表示站点i到站点j的距离;ti表示对站点i的服务时间;
步骤4-4、按照Metropolis准则接收新解,具体包括:在降温过程中根据以下公式接受新的解:
重复步骤(2)-(4)直到迭代次数满足预先设定值。
具体的,所述更新所述最优解所经过的路径上的泥土量具体包括:
其中TIB表示最优解,sw(i,j)b-1为上一次迭代路径中的泥土量,表示最优路径上的水滴携带的泥土量,NIB表示当前迭代得到的最优解中的站点的个数。
具体的,更新全局最优解,具体包括:
当cost(TIB)<cost(TTB)时,TTB=TIB,其中TIB表示最优解,TTB表示全局最优解。
如图2所示,图2为本发明实施例对公共自行车站点划分的区域图。
按照最邻近原则将距离每个调度中心前K个最近的公共自行车站点划分为一个子区域,具体公式如下:
其中N表示公共自行车站点的数量,M表示子区域的数量,通过采用上述方法对自行车站点进行规划,将整个调度区域分为多个子区域,规划好每个调度中心对应管理的子区域。
具体的,所述公共自行车站点的位置信息为坐标信息,所述需求量信息为各个公共自行车站点的调度量di、所述调度量di对应的时域[ei,fi]、所述调度量di对应的最佳时域[ai,bi]。
当在最佳调度时域[ai,bi],di≥0时,通过所述调度中心调入di辆自行车至所述公共自行车站点i,当在最佳调度时域[ai,bi],当di<0时,通过所述调度中心从所述公共自行车站点i调出di辆自行车;
当对公共自行车的调度时刻在[ai,bi]外且在[ei,fi]内时,此时为不满足调度条件的情况,并对所述不满足调度条件的情况进行惩罚,惩罚函数如下:
如图3所示,图3为本发明实施例惩罚函数图。
为了测试本发明中改进的智能水滴算法的有效性以及算法的高效性,本发明设计10个样本数据用来测试混合智能水滴算法。采用Matlab进行算法编程和对问题求解结果模拟仿真。测试所用的10个样本是通过程序随机生成的。所有测试样本的数据格式如表3.4所示。表3.4展示了样本p25的数据。样本数据的类型、对算法的评价指标以及模型的参数设置与混合蚁群算法实验中的一致不再赘述。对于混合智能水滴算法里的相关参数设置如下:
as=1,bs=0.1,cs=1,av=1,bv=0.1,cp=0.6,cw=0.3,α=1,β=1,
表1实验数据格式说明
对于模型和算法的重要参数做如下设置,Tco=21:00,cp=0.6(RMB/KM),tij=150*dij,ei=ai-120,fi=bi+120,Tac=20(minute),cw=0.3。时间片的长度设置为15分钟,运输车的容量可以装载30辆自行车。根据时间限定规则:fi≤Tac+now。
表2算例说明
本发明中采用10个测试样例来对改进的智能水滴算法进行测试,将其与传统智能水滴算法进行比较其求解结果下表所示。
如图4所示,图4为本发明实施例表格2的算例数据求解结果图。
其中IHIWD是本发明最终的算法,IWD-1是本发明提出增加了启发式因子和节约算子以及最大最小调制机制后的改进智能水滴算法。IWD是修改的传统智能水滴算法,对这个算法修改是为了使其可以用于求解本文中的模型,可以用来对比。通过对比每个算法的对每个具体样例的求解时间和求解结果的代价值,可以检测到算法的求解效率。算法的计算时间越短,得到的解决方案的代价值越低,则算法的求解效率越高。解决方案的代价值根据模型的目标函数来计算。
图5为本发明实施例每代路径最小代价对比图,如图5所示,可以看到IWD-1优于IWD,如果不对水滴中的泥土量加以限制,就会出现某条路线泥土量或多或者过少,算法会过早陷入停滞。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (7)
1.一种公共自行车调度方法,其特征在于,所述方法包括:
获取子区域内各公共自行车站点的位置信息以及需求量信息;所述子区域为对各所述公共自行车站点进行划分后形成的子区域;每个所述子区域包含一个调度中心;
根据所述公共自行车站点的位置信息以及需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度,得到最优路径,具体步骤如下:
步骤1、将智能水滴算法中的参数进行初始化,并初始化迭代次数a=1;
步骤2、构造调度方案,具体包括:
步骤2-1、水滴根据概率矩阵,优先选择概率较大的公共自行车站点作为要服务站点,所述概率矩阵公式如下:
其中,τij为启发式算子,dij表示从站点i到站点j之间的距离;ε为大于零的自然数,δ为大零的自然数,FV为需要服务的站点集合,uij为节约算子,表示路径的节约量,uij=di0+d0j-dij,0∈depot,i,j∈(1,2,…N),di0表示站点i到所述调度中心的距离,d0j表示站点j到所述调度中心的距离;其中r=0.01,sw(i,j)为从站点i到站点j的路径arc(i,j)中的泥土量;
判断所述服务站点是否满足约束条件,得到第一判断结果,当所述第一判断结果表示满足约束条件时,则对所述服务站点进行调度,当所述第一判断结果表示不满足约束条件时,则按照水滴到其他所述公共自行车站点的概率值从大到小对站点进行访问,直到找到满足所述约束条件的站点;
所述约束条件为:
0≤dj+wijk≤Q,(i,j∈N∪M),k∈(1,2,…,K)
其中,Q表示每量运输车的最大承载量,dj表示运输车到达公共自行车站点的j距离,tij表示运输车从公共自行车站点i到公共自行车站点j的时间,ti表示运输车到达站点i的时间,wijk表示运输车从公共自行车站点i到公共自行车站点j的自行车数量,wi表示运输车在站点i的停留时间;M表示调度中心的集合,M={SN+1,SN+2,…,SN+M},N表示公共自行车站点的集合,N={v1,v2,…vN};ei、fj表示调度时刻;
如果所有可以访问的公共自行车站点都不满足所述约束条件,则水滴返回调度中心,由调度中心派遣其他运输车辆,当被访问的公共自行车站点满足所述约束条件时,确定当前被访问的公共自行车站点为被服务的下一个自行车站点;
步骤2-2.更新水滴的流速vely、水滴中的泥土量soiliwd y、路径arc(i,j)中的泥土量sw(i,j),具体公式如下:
其中vely为当前水滴的流速,av、bv、cv为速度更新参数,vely-1为上一次迭代中水滴的流速;
soiliwd y=soiliwd y-1+Δsoil(i,j),其中soiliwd y为当前水滴中的泥土量,soiliwd y-1为上一次迭代水滴中的泥土量, as、bs、cs为泥土更新参数,D(i,j)表示从站点i到站点j的路程,ε为大于零的自然数;
sw(i,j)y=sw(i,j)y-1-α*Δsoil(i,j),soilmin≤sw(i,j)y≤soilmax;其中,sw(i,j)y为当前路径中的泥土量,sw(i,j)y-1为上一次迭代路径中的泥土量,ρ为大于零的自然数;TTB表示全局最优解;α的为大于零的自然数;as、bs、cs为泥土更新参数,n为水滴到达目的地的个数;
重复所述步骤2-1和所述步骤2-2直到所有需要调度服务的站点都被服务;
步骤3、获取本次迭代计算得到的最优解,迭代次数a=a+1;
步骤4、更新所述最优解所经过的路径上的泥土量;
步骤5、判断所述迭代次数是否大于预先设定值,得到第二判断结果,若第二判断结果表示所述迭代次数大于所述预先设定值,则执行解的优化;
步骤6、重复所述步骤2-所述步骤5直到迭代次数达到预先设定值,当迭代次数达到预先设定值后输出所有迭代中最优的调度方案。
2.根据权利要求1所述的公共自行车调度方法,其特征在于,所述对各所述公共自行车站点进行划分具体包括:
按照最邻近原则将距离每个调度中心前K个最近的公共自行车站点划分为一个子区域,具体公式如下:
其中N表示公共自行车站点的数量,M表示子区域的数量。
3.根据权利要求1所述的公共自行车调度方法,其特征在于,所述公共自行车站点的位置信息为坐标信息,所述需求量信息为各个公共自行车站点的调度量di、所述调度量di对应的时域[ei,fi]、所述调度量di对应的最佳时域[ai,bi]。
4.根据权利要求1所述的公共自行车调度方法,其特征在于,所述方法还包括:在所述根据所述公共自行车站点的位置信息以及需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度,得到最优路径之后判断对所述公共自行车的调度是否满足调度条件,得到第三判断结果,当第三判断结果表示对所述公共自行车的调度满足所述调度条件时,则根据所述公共自行车站点的位置信息和需求量信息采用智能水滴算法对同一子区域内的公共自行车进行调度:
当在最佳调度时域[ai,bi],di≥0时,通过所述调度中心调入di辆自行车至所述公共自行车站点i,当在最佳调度时域[ai,bi],当di<0时,通过所述调度中心从所述公共自行车站点i调出di辆自行车;
当所述第三判断结果表示对所述公共自行车的调度不满足所述调度条件时,则对不满足调度条件的情况进行惩罚:
当对公共自行车的调度时刻在[ai,bi]外且在[ei,fi]内时,此时为不满足调度条件的情况,对所述不满足调度条件的情况进行惩罚,惩罚函数如下:
其中Cw表示每秒给工作人员的薪资,表示单位里程的油价,ai、bi表示调度时刻,ti表示对站点i的服务时间。
5.根据权利要求1所述的公共自行车调度方法,其特征在于,所述执行解的优化具体包括:
步骤4-1、获取每次迭代中所有最优解,将所述所有最优解中代价值最小的最优解作为全局最优解;所述代价值用cost(Xnew)表示,
其中Cw表示每秒给工作人员的薪资,Cb表示单位里程的油价,ai、bi表示调度时刻,dij表示站点i到站点j的距离;ti表示对站点i的服务时间;
步骤4-2、将所述全局最优解作为模拟退火算法的初始解,所述初始解用X表示,初始化模拟退火算法中的所有参数;
步骤4-3、对所述初始解执行3-opt变换,得到新解,用Xnew表示;
步骤4-4、计算所述Xnew的代价值,所述代价值用cost(Xnew)表示,其中Cw表示每秒给工作人员的薪资,Cb表示单位里程的油价,ai、bi表示调度时刻,dij表示站点i到站点j的距离;ti表示对站点i的服务时间;
步骤4-4、按照Metropolis准则接收新解。
重复步骤(2)-(4)直到迭代次数满足预先设定值。
6.根据权利要求1所述的公共自行车调度方法,其特征在于,所述更新所述最优解所经过的路径上的泥土量具体包括:
其中TIB表示最优解,sw(i,j)b-1为上一次迭代路径中的泥土量,表示最优路径上的水滴携带的泥土量,NIB表示当前迭代得到的最优解中的站点的个数。
7.根据权利要求1所述的公共自行车调度方法,其特征在于,更新全局最优解,具体包括:
当cost(TIB)<cost(TTB)时,TTB=TIB,其中TIB表示最优解,TTB表示全局最优解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810475790.7A CN108805335B (zh) | 2018-05-17 | 2018-05-17 | 一种公共自行车调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810475790.7A CN108805335B (zh) | 2018-05-17 | 2018-05-17 | 一种公共自行车调度方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108805335A true CN108805335A (zh) | 2018-11-13 |
CN108805335B CN108805335B (zh) | 2021-10-22 |
Family
ID=64092521
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810475790.7A Expired - Fee Related CN108805335B (zh) | 2018-05-17 | 2018-05-17 | 一种公共自行车调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108805335B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111080006A (zh) * | 2019-12-13 | 2020-04-28 | 河海大学常州校区 | 一种社区居家养老服务调度模型与求解方法 |
CN111091286A (zh) * | 2019-12-13 | 2020-05-01 | 河海大学常州校区 | 一种公共自行车调度模型与求解方法 |
CN111090935A (zh) * | 2019-12-13 | 2020-05-01 | 河海大学常州校区 | 一种公共自行车预约调度与路径规划方法 |
CN111598343A (zh) * | 2020-05-18 | 2020-08-28 | 武汉轻工大学 | 配送路径优化方法、设备及可读存储介质 |
CN112085637A (zh) * | 2020-08-17 | 2020-12-15 | 杭州金通科技集团股份有限公司 | 一种公共自行车驿站一键停运方法 |
CN114202911A (zh) * | 2022-02-18 | 2022-03-18 | 北京阿帕科蓝科技有限公司 | 一种车辆调度方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105426997A (zh) * | 2015-11-18 | 2016-03-23 | 西南交通大学 | 城市公共自行车智能调度与运用管理系统及其实现方法 |
US20170270448A1 (en) * | 2016-03-18 | 2017-09-21 | Jusda International Logistics (TAIWAN) CO.,LTD | Vehicle scheduling device and method for transportation systems |
CN107766994A (zh) * | 2017-12-04 | 2018-03-06 | 长沙理工大学 | 一种共享自行车调度方法与调度系统 |
CN108038578A (zh) * | 2017-12-28 | 2018-05-15 | 东南大学 | 基于需求预测和中心辐射网络的公共自行车静态调度方法 |
-
2018
- 2018-05-17 CN CN201810475790.7A patent/CN108805335B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105426997A (zh) * | 2015-11-18 | 2016-03-23 | 西南交通大学 | 城市公共自行车智能调度与运用管理系统及其实现方法 |
US20170270448A1 (en) * | 2016-03-18 | 2017-09-21 | Jusda International Logistics (TAIWAN) CO.,LTD | Vehicle scheduling device and method for transportation systems |
CN107766994A (zh) * | 2017-12-04 | 2018-03-06 | 长沙理工大学 | 一种共享自行车调度方法与调度系统 |
CN108038578A (zh) * | 2017-12-28 | 2018-05-15 | 东南大学 | 基于需求预测和中心辐射网络的公共自行车静态调度方法 |
Non-Patent Citations (1)
Title |
---|
胡云清: "改进智能水滴算法在车辆调度问题中的应用", 《包装工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111080006A (zh) * | 2019-12-13 | 2020-04-28 | 河海大学常州校区 | 一种社区居家养老服务调度模型与求解方法 |
CN111091286A (zh) * | 2019-12-13 | 2020-05-01 | 河海大学常州校区 | 一种公共自行车调度模型与求解方法 |
CN111090935A (zh) * | 2019-12-13 | 2020-05-01 | 河海大学常州校区 | 一种公共自行车预约调度与路径规划方法 |
CN111080006B (zh) * | 2019-12-13 | 2022-08-16 | 河海大学常州校区 | 一种社区居家养老服务调度方法 |
CN111091286B (zh) * | 2019-12-13 | 2022-08-16 | 河海大学常州校区 | 一种公共自行车调度方法 |
CN111090935B (zh) * | 2019-12-13 | 2022-09-23 | 河海大学常州校区 | 一种公共自行车预约调度与路径规划方法 |
CN111598343A (zh) * | 2020-05-18 | 2020-08-28 | 武汉轻工大学 | 配送路径优化方法、设备及可读存储介质 |
CN112085637A (zh) * | 2020-08-17 | 2020-12-15 | 杭州金通科技集团股份有限公司 | 一种公共自行车驿站一键停运方法 |
CN112085637B (zh) * | 2020-08-17 | 2024-02-02 | 杭州金通科技集团股份有限公司 | 一种公共自行车驿站一键停运方法 |
CN114202911A (zh) * | 2022-02-18 | 2022-03-18 | 北京阿帕科蓝科技有限公司 | 一种车辆调度方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108805335B (zh) | 2021-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108805335A (zh) | 一种公共自行车调度方法 | |
Sayarshad et al. | Non-myopic dynamic routing of electric taxis with battery swapping stations | |
CN109711630A (zh) | 一种基于出行概率矩阵的电动汽车快充站选址定容方法 | |
CN109034465A (zh) | 考虑充电站选址与出行路径耦合的充电站两层规划方法 | |
CN107490386A (zh) | 一种电动汽车最优路径和驾驶方式的规划方法及系统 | |
CN107187442A (zh) | 基于工况预测的插电式混合动力电动汽车能量管理系统 | |
CN109558988A (zh) | 一种基于大数据融合的电动汽车能耗预测方法及系统 | |
CN106530680B (zh) | 一种基于大站快车的公交线路组合服务方法 | |
CN108765948A (zh) | 一种弹性的公车调度方法和系统 | |
CN106228841A (zh) | 一种城市级停车诱导方法、装置及电子设备 | |
CN105678425B (zh) | 基于多节拍组合的城际铁路列车开行的优化方法 | |
CN110222906A (zh) | 电动车能耗预测方法、计算机可读存储介质和电子设备 | |
CN104778832B (zh) | 一种多方式的城市道路横断面资源协调配置方法 | |
CN113393690A (zh) | 考虑自动驾驶专用车道的混合交通均衡分配方法 | |
CN115689083A (zh) | 一种预见性能量管理方法、装置、电子设备及存储介质 | |
CN110598971A (zh) | 一种基于蚁群算法的响应式公交服务规划方法 | |
CN114842641B (zh) | 面向省域的多模式链式交通分配方法 | |
CN113450564A (zh) | 一种基于narx神经网络和c-v2x技术的交叉路口的通行方法 | |
Annema | Policy implications of the potential carbon dioxide (CO2) emission and energy impacts of highly automated vehicles | |
Cai et al. | Predicting the carbon emission reduction potential of shared electric bicycle travel | |
Zhao et al. | Optimization of intensive land use in blocks of Xi’an from the perspective of bicycle travel | |
Pan et al. | Research on multi-lane energy-saving driving strategy of connected electric vehicle based on vehicle speed prediction | |
CN111127076A (zh) | 一种出租车共享模式下的定价方法 | |
CN112632374B (zh) | 一种考虑定制公交的居民出行方式选择分析方法 | |
CN107909201A (zh) | 基于广义出行费用的交通方式优势出行距离的量化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211022 |