CN109344529A - 一种基于两阶段启发式算法的定制公交线路设计方法 - Google Patents
一种基于两阶段启发式算法的定制公交线路设计方法 Download PDFInfo
- Publication number
- CN109344529A CN109344529A CN201811227668.4A CN201811227668A CN109344529A CN 109344529 A CN109344529 A CN 109344529A CN 201811227668 A CN201811227668 A CN 201811227668A CN 109344529 A CN109344529 A CN 109344529A
- Authority
- CN
- China
- Prior art keywords
- vehicle
- website
- passenger
- time
- customization
- 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
- 238000013461 design Methods 0.000 title claims abstract description 38
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000032258 transport Effects 0.000 claims description 29
- 238000005457 optimization Methods 0.000 claims description 18
- 238000003780 insertion Methods 0.000 claims description 11
- 230000037431 insertion Effects 0.000 claims description 11
- 238000013459 approach Methods 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 3
- 238000009825 accumulation Methods 0.000 claims description 2
- 238000003032 molecular docking Methods 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 7
- 238000007726 management method Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明公开了一种基于两阶段启发式算法的定制公交线路设计方法,包括定制公交需求获取、带时间窗的定制公交线路设计模型的构建、基于两阶段启发式算法的设计以及求解模型得到定制公交线路方案等步骤。本方法考虑了多区域定制公交线路设计问题,在兼顾乘客出行成本和定制公交运营成本的基础上构建定制公交线路设计模型并输出线路方案,为定制公交的实际规划与运营提供参考。与现有定制公交线路设计问题的研究相比,本发明考虑了乘客在时间与空间需求,并设计了一种两阶段启发式算法,算法的输出包括线路途经的站点,车辆在每个站点到达时间与运营时间等信息。本发明所得结果可为实际定制公交线路的规划,以及时刻表编制提供理论方法。
Description
技术领域
本发明涉及公共交通信息处理技术领域,具体地说是一种基于两阶段启发式算法的定制公交线路设计方法。
背景技术
近年来,在城市居民出行需求多样化和常规公交运营服务模式单一化的矛盾下,鼓励构建多元化的公共交通服务模式已经成为城市交通发展的必然趋势。在此背景下,定制公交作为以缓解地面交通拥堵以及满足乘客多元化出行需求为目标的一种需求响应式公共交通服务模式应运而生。
定制公交作为多元化公共交通服务模式的组成部分,是指通过集中整合个体的交通出行需求,为出行起终点、出行时间、服务需求相同或相似的人群提供专门定制的公共交通服务方式。定制公交是城市公共交通体系中的一种辅助客运公交服务系统,是在多元化公交服务模式的前提下,作为个性化和精细化出行需求市场的一种补充,是与城市主要公共交通方式相辅相成、相互配合的一种形式。
定制公交作为一种新型的公共交通服务模式,由于在国内的发展时间较短,国外类似需求响应的交通方式无法直接套用在国内现状中,而国内的经营方式也多依赖于常规公交的理论与实践经验,因此从目前的情况来看还存在一些问题。在线路设计层面,多数研究集中在单一起终点之间的路径规划,直达专线的设计,对多起终点之间的线路设计问题分析不够。且模型的求解多直接采用传统的启发式方法,如遗传算法,蚁群算法等。
发明内容
针对目前技术存在的不足,本发明的目的在于提供一种基于两阶段启发式算法的定制公交线路设计方法,即在定制公交乘客出行需求提取的基础上,在满足乘客出行需求的前提下,兼顾乘客的出行成本与运营成本,构建带时间窗的定制公交线路设计模型,并设计一种两阶段启发式算法求解模型,输出定制公交线路方案。本发明的实现包括以下步骤:
第一步,获取定制公交出行需求。该需求可以包括以下几个方面。第一,由出行调查获取的出行需求。第二,由定制公交相关运营机构通过网页或手机软件获取的出行需求数据。第三,由地图导航软件所提取的用户出行需求数据,一般包括驾车导航数据与出行规划数据两个方面。
第二步,将所获得的出行需求进行整理,提取出定制公交线路设计所需要的出行信息。具体包括乘客的出发地、目的地、期望出发时间及期望到达时间。
第三步,构建带时间窗的定制公交线路设计模型,该模型在满足所有乘客需求的前提下,以乘客的出行时间成本最小和定制公交运营成本最小为优化目标。
第四步,采用一种两阶段启发式算法求解该模型,算法第一阶段采用一种基于概率的站点选择策略优化乘客出行成本。算法第二阶段通过不同算子优化运营成本。
最后,通过求解线路设计模型,输出定制公交线路集合,每条线路包括途经站点,每个站点到达时间与运营时间等信息。
本发明的优点在于:
1.本发明综合考虑了定制公交乘客的出行需求,不仅考虑了考虑乘客在空间上的出行需求,同时模型还考虑了乘客在时间上的出行需求,即线路的设计要遵循乘客的出行时空窗。2.不同于多数研究集中在单区域或单起终点的定制公交线路设计问题,本发明所构建的模型考虑多区域(多起终点)线路设计问题。
3.综合考虑乘客的出行成本与定制公交运营成本构建定制公交多目标线路设计模型。不同于大部分研究多采用传统启发式算法直接求解模型,本发明设计了一种两阶段启发式方法生成定制公交线路集合,可为定制公交的实际运营提供支撑。
附图说明
图1为本发明的方法流程图;
图2为两阶段启发式算法中插入算子示意图
图3为实施例中定制公交线路设计模型示意图
图4为实施例中出行需求示意图
图5为定制公交线路类型示意图。
具体实施方式
下面结合附图和实施例对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
第一步,获取定制公交出行需求。该需求可以包括以下几个方面。第一,由出行调查获取的出行需求。第二,由定制公交相关运营机构通过网页或手机软件获取的出行需求数据。第三,由地图导航软件所提取的用户出行需求数据,一般包括驾车导航数据与出行规划数据两个方面。
第二步,将所获得的出行需求进行整理,提取出定制公交线路设计所需要的出行信息。对于每位乘客p,提取乘客出发地o(p)、目的地d(p)、期望出发时间期望到达时间因此每个乘客p的出行需求信息可由下式表示:
因此,每位乘客p都有一个时间窗。考虑实际情况,在线路设计时给增加一个15到30分钟的缓冲时间w,即车辆要在时间段内到达乘客上车站点,且要在之前到达目的地d(p)。
第三步,构建带时间窗的定制公交线路设计模型,该模型在满足所有乘客需求的前提下,以乘客的出行时间成本最小和定制公交运营成本最小为优化目标。
模型中的变量如表1和表2所示。本发明主要针对早高峰通勤乘客需求所构建模型。
表1模型中各参数含义
参数 | 含义 |
N | 定制公交站点集合 |
P | 乘客集合 |
p | 乘客标识 |
K | 运营车辆集合 |
k | 车辆标识 |
H | 居住地区域定制公交候选站点 |
W | 工作地区域定制公交候选站点 |
Q(p) | 乘客p出行需求 |
C(k) | 车辆k座位容量 |
b(k) | 车辆k出发站点(场站) |
e(k) | 车辆k到达站点(场站) |
T(i,j) | 站点i和j之间的直达旅行时间 |
SE<sub>p</sub>(i) | 乘客p在站点i的服务时间 |
v(k) | 车辆k的平均行驶速度 |
r<sub>p</sub>(i) | 站点i的上车客流量 |
L<sub>k</sub>(i) | 车辆k在访问过站点i后的累积客流 |
D(i,j) | 站点i和j之间的距离 |
V<sub>k</sub> | 由车辆k所服务的线路 |
S | 由V<sub>k</sub>组成的线路集合 |
tt(i,j) | 站点i和j之间的旅行时间 |
表2模型中决策变量含义
(1)优化乘客出行成本
第一个目标函数为优化乘客出行成本,即最小化乘客的损失时间,该目标函数定义如下:
目标函数(2)表示最小化选择定制公交出行的乘客的人均损失时间。损失时间f(p,k)可由下式来表示:
其中表示车辆k在站点i的出发时间,可由车辆k到达站点i的时间与服务时间SEp(i)计算:
为乘客p到达目的地j的时间,T(i,j)为站点i与站点j之间不经过中间站点而直达所需要的旅行时间,也可以视作私家车或者出租车所用的时间。这样,式(3)定义的乘客损失时间可由两部分组成,一部分是提前到达目的地所损失的时间(需要等待工作时间开始);第二部分为与出租车的旅行时间成本相比较所损失的时间。|p|表示具有相同出行需求特征Q(p)的人数;为决策变量。因此,目标函数(2)表示优化乘客的平均损失时间。
(2)优化运营成本
第二个目标函数为优化运营成本,本发明中,运营成本目标为最小化运营车辆数:
其中b(k)为车辆k的出发站点,为决策变量,(5)式表示最小化线路的运营车辆数。
(3)模型包含以下约束条件:
约束(6)表示每位乘客只能由一辆车辆服务。
约束(7)为时间约束,表示每位乘客p到达目的地时间不能晚于期望到达时间
约束(8)限定车辆在运营时应先访问乘客p的出发地o(p),再访问目的地d(p)。
约束(9)为时间约束,若i和j为一条线路两个站点,若车辆k先访问i再访问j,则时间上应满足约束条件(9),其中G为一个值很大的常数。约束(9)表示车辆到达j的时间应大于到达站点i的时间、服务时间和i到j的旅行时间之和。
约束(10)和(11)与约束(9)类似,为车辆出发站点b(k)与到达站点e(k)时间约束。
D(si+1,s1)>D(si,s1), (12)
D(si+1,sn)<D(si,sn), (13)
约束(12)与(13)为站点间的方向约束,此约束是为了保证生成的线路非直线系数较低,以降低站点之间的旅行时间。约束(12)表示在线路的设计中,新插入的站点si+1离出发地s1相较于si要远。约束(13)表示新插入的站点si+1离目的地sn相较于si要近。
约束(14)限定车辆的负载。不同于常规公交,定制公交需保证“一人一座”以提升乘客的舒适性。因此,当车辆访问某个站点后,车辆上的乘客数量不能超过车辆的容量限制。
约束(15)为流平衡,也表示为车辆k访问站点i后需要离开i。
式(16-18)为变量约束。
第四步,本发明采用一种两阶段启发式算法求解该模型,算法为迭代算法,每次迭代包括两个阶段。算法第一阶段为优化乘客出行成本。算法第二阶段通过不同算子优化运营成本(车辆数)。
算法的第一阶段为优化出行成本,主要原则是将同一个时间段到达同一个目的地(工作地)的乘客一起运送。具体包括以下步骤:
(1)对于同一个时间段到达同一个目的地d的乘客,令其出发地集合为so={s1,s2,…si…,sn}。在车辆从起点b(k)出发之后,从so随机选择一个乘客的出行起点si(k)。这样,由车辆k服务的线路可以初始化定义为Vk={a1(k),d},a1(k)=si(k)。
(2)采用一种基于概率的选择算法选择后续车辆k要访问的定制公交站点。站点sj(k)被选择的概率为:
其中概率P(sj)主要包含两个部分,第一部分为前后两个候选站点sj-1(k)和sj(k)的旅行时间,即T(sj-1(k),sj(k))。第二个部分是候选站点sj(k)到目的地d的旅行时间。因此,式(19)中分子表示这两部分旅行时间之和。分母表示对每个待选择的候选站点计算两部分旅行时间并加和。对于式(19),若sj(k)到sj-1(k)之间旅行时间越小,同时sj(k)到目的地d旅行时间越小,则sj(k)被选择的概率越大。
(3)在选择sj(k)之后,要判断新插入的站点sj(k)是否满足模型的约束条件(6)-(15),包括时间约束、车辆容量约束以及方向约束。如果满足,将站点放入车辆k所服务的线路Vk之中,并更新线路。更新之后将sj(k)从so中删除。否则,返回步骤2直到选择一个可行的站点插入。
(4)重复步骤(2)和(3)直到没有站点可以添加到Vk之中,这样该线路的目标函数(2)的值可以被确定。由车辆k运送的线路Vk也被确定。
(5)重复步骤(1)-(4)直到每个目的地d的需求都被覆盖,这样就得到一个初始的线路集合S。
算法的第二阶段为优化运营成本,即最小化运营车辆数量。在第一阶段生成初始的线路集合基础上,第二阶段对于车辆数的优化主要分为以下两种情形。
第一种情形是对两个线路进行合并操作,主要包括以下步骤。
在定制公交的实际运营中,对于运营的车辆k,在服务完一个车次后,可以选择回到场站(定制公交可以灵活选择场站,常站可以设置在居住区),也可以继续服务下个车次。因此,若某车辆k在服务完一个车次后空驶到达下个车次或线路的起点如果满足以下条件,可以对两个车次或线路进行合并由同一车辆服务:
Tk1(dk1)+se(dk1)+T(dk1,ok2)≤Ok2, (20)
其中Tk1(dk1)表示车辆k1到达第一个车次目的地dk1的时间,se(dk1)表示车辆k1在目的地站点dk1的服务时间,T(dk1,ok2)为第一个车次目的地dk1到第二个车次出发地ok2的旅行时间,Ok2是第二个车次的运营车辆k2的出发时间。该不等式(20)表示如果第二个车次与第一个车次相互间没有任何影响,那么两个车次可以进行合并,由一辆车进行服务。因此车辆的数量由2变为1。
第二种情形,采用一种插入算子将一个车次中的站点分配到其他的车次中,以减少车辆数,具体步骤如下。
令TR={TR1,TR2,…,TRn}为第一阶段产生的车次集合。
(1)从TR中随机选择一个车次TRi,TRi中的站点将被分配到其他车次中。
(2)从TR\{TRi}中再选择一个车次TR*,检测TRi中的目的地是否与TR*中的目的地一样。如果不一致,算子将尝试在TR*中搜索可以插入的位置。
(3)在TR*中搜索可以插入的位置,将TRi中的站点分配到TR*中。注意到增加车次或线路中的中间站点数量会影响出行成本。因此,为平衡乘客出行成本与运营成本,可对插入站点的数量进行限制。如果TR*中没有可以插入的位置,那么算子返回步骤(2)继续搜索可以插入的车次。
(4)重复步骤(2)和(3)直到TRi中的站点全部被分配到其余站点。如果TRi中的站点全部被分配到其余站点,从TR中删除TRi和TR*;否则,TR与TRi保持不变。
(5)重复步骤(1)到(4)直到TR中没有多余的车次可以被分配到其他的车次中,这样得到新的线路集合。算子的步骤可由图2表示。
另外,插入站点时需要考虑模型的约束条件;且当新的站点插入时,会影响目标函数(2)的值。因此,在每个站点插入后,需更新每条车次或线路的信息。
实施例
一种基于两阶段启发式算法的定制公交线路设计方法,具体如下:
本实施例选择在北京国贸地区工作,居住地分布在通州区、北三县、廊坊地区的出行需求。本发明中,模型可以被看作多区域对单区域的问题,如图3所示。每个居住地区域有若干个定制公交乘客上车站点,每个工作地区域有若干个定制公交乘客下车站点,居住地与工作地之间可以由不同的线路进行运送,即多区域(多起终点)定制公交线路设计问题。实施例所包含的出行需求如图4所示。经过对需求的整理,在本例中,三个居住地区域共有63个上车站点,一个工作地区域有5个主要的下车站点。此外,对于车型的设置本例为单一车型,每辆车限定容量为50人。站点之间的旅行时间按下式计算:
D(i,j)表示两个站点间的直线距离,v(k)代表车辆的平均速度,这里假设车辆匀速行驶。私家车的行驶速度为40公里/小时,定制公交的行驶速度为30公里/小时。上车乘客的服务时间按下式计算:
se1=19+2.6·ni, (22)
其中se1以秒为单位,ni为在站点i的上车乘客数量。下车乘客的服务时间按下式计算:
se2=29+1.9·ni, (23)
其中se2以秒为单位,ni为在站点i的下车乘客数量。
接下来应用两阶段启发式算法对模型进行求解,迭代次数为5000次。由于两个目标函数相互冲突,因此模型求解的目标为非支配解。经过计算,最优乘客出行成本和最优运营成本的结果如表3所示。
表3定制公交线路设计模型测试结果
参数 | 方案A | 方案B |
出行成本 | 15.83min | 17.12min |
运营成本 | 30辆 | 28辆 |
如表3所示,方案A为乘客出行成本最优的方案,方案B为运营成本最低的方案。可以看出,根据不同的优化目标,算法可以输出不同的线路方案。在实际运营中,规划人员可根据不同情况在算法生成的线路集中选择方案并进行调整以适应需要。本发明所采用的算法所产生的线路主要包括三种类型。图5对这三种类型进行示例。第一种包含3个出发地和1个目的地,线路走向为A-B-C-W1(图5(a))。第二种对应合并两个车次(线路)的情况,包括5个出发地和2个目的地,第一个车次为A-B-C-D-W1,之后车辆到达站点E执行第二个车次E-W2(图5(b))。第三种对应插入站点的情况,包括3个出发地与2个目的地,线路走向为A-B-C-W1-W2,其中C与W2为插入站点(图5(c))。三种情况的线路基本信息如表4-6所示。
以上详细描述了本发明的实施案例,但是本发明并不局限于上述实施案例的具体细节,在本发明的整体结构范围内,可以对本发明的部分步骤进行多种变换并重新组合,本发明对各种可能的组合方式不再列举,这些变换组合均属于本发明的保护范围。
表4线路基本信息(对应图5(a))
表5线路基本信息(对应图5(b))
表6线路基本信息(对应图5(c))
Claims (4)
1.一种基于两阶段启发式算法的定制公交线路设计方法,其特征在于,包括以下步骤:
第一步,获取定制公交出行需求
该需求可以包括:第一,由出行调查获取的出行需求;第二,由定制公交相关运营机构通过网页或手机软件获取的出行需求数据;第三,由地图导航软件所提取的用户出行需求数据,包括驾车导航数据与出行规划数据;
第二步,将所获得的出行需求进行整理,提取出定制公交线路设计所需要的出行信息
所述信息包括乘客的出发地、目的地、期望出发时间及期望到达时间;
第三步,构建带时间窗的定制公交线路设计模型,该模型在满足所有乘客需求的前提下,以乘客的出行时间成本最小和定制公交运营成本最小为优化目标;
第四步,采用一种两阶段启发式算法求解该模型
算法第一阶段采用一种基于概率的站点选择策略优化乘客出行成本;算法第二阶段通过不同算子优化运营成本;
第五步,通过求解线路设计模型,输出定制公交线路集合,每条线路包括途经站点,每个站点到达时间与运营时间等信息。
2.根据权利要求1所述的一种基于两阶段启发式算法的定制公交线路设计方法,其特征在于,第一步中,需求数据的来源包括出行调查数据、定制公交相关运营机构通过网页或手机软件获取的出行需求数据、地图导航软件所提取的用户出行需求数据,也可以包含其他相关数据。
3.根据权利要求1所述的一种基于两阶段启发式算法的定制公交线路设计方法,其特征在于,第三步中,构建的定制公交线路设计模型如下所述。
第一个目标函数为优化乘客出行成本:
第二个目标函数为优化定制公交运营成本:
其中,P为乘客集合;p为乘客标识;K为运营车辆集合;k为车辆标识;f(p,k)为损失时间;若乘客p由车辆k从站点i运送到站点j,则为1,否则为0;b(k)为车辆k出发站点为如果i和j分别为由车辆k服务的前后两个相邻站点,则为1;否则为0;
此外,该模型包含以下约束条件:
约束(1)表示每位乘客只能由一辆车辆服务
约束(2)为时间约束,表示每位乘客p到达目的地时间不能晚于期望到达时间
约束(3)限定车辆在运营时应先访问乘客p的出发地o(p),再访问目的地d(p)
约束(4)为时间约束,若i和j为一条线路两个站点,若车辆k先访问i再访问j,则时间上应满足约束条件(4),其中G为一个值很大的常数。约束(4)表示车辆到达j的时间应大于到达站点i的时间、服务时间和i到j的旅行时间之和。
约束(5)和(6)与约束(4)类似,为车辆出发站点b(k)与到达站点e(k)时间约束
约束(7)与(8)为站点间的方向约束,此约束是为了保证生成的线路非直线系数较低,以降低站点之间的旅行时间。约束(7)表示在线路的设计中,新插入的站点si+1离出发地s1相较于si要远。约束(8)表示新插入的站点si+1离目的地sn相较于si要近。
D(si+1,s1)>D(si,s1), (7)
D(si+1,sn)<D(si,sn), (8)
约束(9)限定车辆的负载。不同于常规公交,定制公交需保证“一人一座”以提升乘客的舒适性。因此,当车辆访问某个站点后,车辆上的乘客数量不能超过车辆的容量限制。
约束(10)为流平衡,也表示为车辆k访问站点i后需要离开i:
约束(11)(12)(13)为变量约束:
式(1)-(13)中,N为定制公交站点集合;P为乘客集合;p为乘客标识;K为运营车辆集合;k为车辆标识;H为居住地区域定制公交候选站点;W为工作地区域定制公交候选站点;Q(p)为乘客p出行需求;C(k)为车辆k座位容量;b(k)为车辆k出发站点;e(k)为车辆k到达站点;T(i,j)为站点i和j之间的直达旅行时间;SEp(i)为乘客p在站点i的服务时间;v(k)为车辆k的平均行驶速度;rp(i)为站点i的上车客流量;Lk(i)为车辆k在访问过站点i后的累积客流;D(i,j)为站点i和j之间的距离;Vk为由车辆k所服务的线路;S为由Vk组成的线路集合;tt(i,j)为站点i和j之间的旅行时间;若乘客p由车辆k从站点i运送到站点j,则为1,否则为0;b(k)为车辆k出发站点为如果i和j分别为由车辆k服务的前后两个相邻站点,则为1;否则为0,Ttk(i)为线路t中车辆k到达站点i的时间。
4.根据权利要求3所述的一种基于两阶段启发式算法的定制公交线路设计方法,其特征在于,第四步中,基于两阶段的启发式算法,步骤如下:
算法为迭代算法,每次迭代包括两个阶段;算法第一阶段为优化乘客出行成本;算法第二阶段通过不同算子优化运营成本(车辆数);
算法的第一阶段为优化出行成本,是将同一个时间段到达同一个目的地(工作地)的乘客一起运送;具体包括以下步骤:
S101对于同一个时间段到达同一个目的地d的乘客,令其出发地集合为so={s1,s2,…si…,sn};在车辆从起点b(k)出发之后,从so随机选择一个乘客的出行起点si(k);这样,由车辆k服务的线路可以初始化定义为Vk={a1(k),d},a1(k)=si(k)
S102采用一种基于概率的选择算法选择后续车辆k要访问的定制公交站点;站点sj(k)被选择的概率为:
其中概率P(sj)主要包含两个部分,第一部分为前后两个候选站点sj-1(k)和sj(k)的旅行时间,即T(sj-1(k),sj(k));第二个部分是候选站点sj(k)到目的地d的旅行时间;因此,式中分子表示这两部分旅行时间之和;分母表示对每个待选择的候选站点计算两部分旅行时间并加和;若sj(k)到sj-1(k)之间旅行时间越小,同时sj(k)到目的地d旅行时间越小,则sj(k)被选择的概率越大;
S103在选择sj(k)之后,要判断新插入的站点sj(k)是否满足模型的约束条件(1)-(13),包括时间约束、车辆容量约束以及方向约束;如果满足,将站点放入车辆k所服务的线路Vk之中,并更新线路;更新之后将sj(k)从so中删除;否则,返回步骤2直到选择一个可行的站点插入;
S104重复步骤S102和S103直到没有站点可以添加到Vk之中,这样该线路的出行成本目标函数的值可以被确定;由车辆k运送的线路Vk也被确定;
S105重复步骤S101和S104直到每个目的地d的需求都被覆盖,这样就得到一个初始的线路集合S;
算法的第二阶段为优化运营成本,即最小化运营车辆数量;在第一阶段生成初始的线路集合基础上,第二阶段对于车辆数的优化主要分为以下两种情形;
第一种情形是对两个线路进行合并操作,包括以下步骤:
在定制公交的实际运营中,对于运营的车辆k,在服务完一个车次后,可以选择回到场站(定制公交可以灵活选择场站,常站可以设置在居住区),也可以继续服务下个车次;因此,若某车辆k在服务完一个车次后空驶到达下个车次或线路的起点如果满足以下条件,可以对两个车次或线路进行合并由同一车辆服务:
Tk1(dk1)+se(dk1)+T(dk1,ok2)≤Ok2
其中Tk1(dk1)表示车辆k1到达第一个车次目的地dk1的时间,se(dk1)表示车辆k1在目的地站点dk1的服务时间,T(dk1,ok2)为第一个车次目的地dk1到第二个车次出发地ok2的旅行时间,Ok2是第二个车次的运营车辆k2的出发时间,本不等式表示如果第二个车次与第一个车次相互间没有任何影响,那么两个车次可以进行合并,由一辆车进行服务;因此车辆的数量由2变为1;
第二种情形,采用一种插入算子将一个车次中的站点分配到其他的车次中,以减少车辆数,具体步骤如下:
令TR={TR1,TR2,…,TRn}为第一阶段产生的车次集合;
S201从TR中随机选择一个车次TRi,TRi中的站点将被分配到其他车次中;
S202从TR\{TRi}中再选择一个车次TR*,检测TRi中的目的地是否与TR*中的目的地一样;如果不一致,算子将尝试在TR*中搜索可以插入的位置;
S203在TR*中搜索可以插入的位置,将TRi中的站点分配到TR*中;注意到增加车次或线路中的中间站点数量会影响出行成本;因此,为平衡乘客出行成本与运营成本,可对插入站点的数量进行限制;如果TR*中没有可以插入的位置,那么算子返回步骤(2)继续搜索可以插入的车次;
S204重复步骤(2)和(3)直到TRi中的站点全部被分配到其余站点;如果TRi中的站点全部被分配到其余站点,从TR中删除TRi和TR*;否则,TR与TRi保持不变;
S205重复步骤(1)到(4)直到TR中没有多余的车次可以被分配到其他的车次中,这样得到新的线路集合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811227668.4A CN109344529B (zh) | 2018-10-22 | 2018-10-22 | 一种基于两阶段启发式算法的定制公交线路设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811227668.4A CN109344529B (zh) | 2018-10-22 | 2018-10-22 | 一种基于两阶段启发式算法的定制公交线路设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109344529A true CN109344529A (zh) | 2019-02-15 |
CN109344529B CN109344529B (zh) | 2023-11-03 |
Family
ID=65310623
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811227668.4A Active CN109344529B (zh) | 2018-10-22 | 2018-10-22 | 一种基于两阶段启发式算法的定制公交线路设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109344529B (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109934401A (zh) * | 2019-03-09 | 2019-06-25 | 吉林大学 | 一种基于移动互联的定制公交方法 |
CN110598971A (zh) * | 2019-07-25 | 2019-12-20 | 中山大学 | 一种基于蚁群算法的响应式公交服务规划方法 |
CN110793535A (zh) * | 2019-08-22 | 2020-02-14 | 腾讯科技(深圳)有限公司 | 一种通勤线路规划方法、装置 |
CN111105078A (zh) * | 2019-11-29 | 2020-05-05 | 同济大学 | 一种定制公交网络优化方法 |
CN111144618A (zh) * | 2019-12-04 | 2020-05-12 | 东南大学 | 一种基于两阶段优化模型的需求响应型定制公交网络规划方法 |
CN111832882A (zh) * | 2020-04-30 | 2020-10-27 | 北京嘀嘀无限科技发展有限公司 | 一种交通运输的控制方法、装置、存储介质和电子设备 |
CN112085249A (zh) * | 2020-07-27 | 2020-12-15 | 北京工业大学 | 一种基于强化学习的定制公交线路规划方法 |
CN112085349A (zh) * | 2020-08-19 | 2020-12-15 | 大连海事大学 | 一种基于乘客出行时间窗约束的需求响应公交调度方法 |
CN112419128A (zh) * | 2020-12-16 | 2021-02-26 | 武汉元光科技有限公司 | 一种线路规划方法以及相关设备 |
CN112419704A (zh) * | 2020-11-06 | 2021-02-26 | 杭州图软科技有限公司 | 一种基于大数据的公共交通路线规划方法及其系统 |
CN112466122A (zh) * | 2021-01-28 | 2021-03-09 | 深圳市城市交通规划设计研究中心股份有限公司 | 一种公交线网的备选线路集生成、线路规划方法及装置 |
CN112700850A (zh) * | 2020-12-30 | 2021-04-23 | 杭州电子科技大学 | 一种基于三阶段混合启发式算法的家庭医疗护理调度随机优化方法 |
CN112734188A (zh) * | 2020-12-30 | 2021-04-30 | 杭州电子科技大学 | 一种基于两阶段混合元启发式算法的家庭医疗护理调度优化方法 |
CN112750063A (zh) * | 2021-01-04 | 2021-05-04 | 李璐 | 基于随机规划的公交车队设施选址-路径规划-调度方法 |
CN113096430A (zh) * | 2021-03-09 | 2021-07-09 | 东南大学 | 基于欧氏距离收费的弹性公交即时需求动态响应方法 |
CN113469451A (zh) * | 2021-07-19 | 2021-10-01 | 杭州数知梦科技有限公司 | 基于启发式算法的定制公交线路生成方法 |
CN113538886A (zh) * | 2021-06-10 | 2021-10-22 | 广东工业大学 | 一种实时响应式定制公交递阶调度方法 |
CN114518763A (zh) * | 2022-04-20 | 2022-05-20 | 清华大学深圳国际研究生院 | 一种协同公交车的路径规划方法 |
CN114897445A (zh) * | 2022-07-12 | 2022-08-12 | 深圳市城市交通规划设计研究中心股份有限公司 | 公交线网的停靠点调整优化方法、装置及可读存储介质 |
CN115641704A (zh) * | 2022-12-26 | 2023-01-24 | 东风悦享科技有限公司 | 一种智能公交调度方法及系统 |
CN116090785A (zh) * | 2023-02-21 | 2023-05-09 | 哈尔滨工业大学 | 针对大型活动散场场景两阶段的定制公交规划方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106448233A (zh) * | 2016-08-19 | 2017-02-22 | 大连理工大学 | 基于大数据的公交线路时刻表协同优化方法 |
CN107194128A (zh) * | 2017-06-30 | 2017-09-22 | 东南大学 | 基于中心辐射型网络框架的多模式公共交通网络设计方法 |
CN107330559A (zh) * | 2017-07-03 | 2017-11-07 | 华南理工大学 | 一种多起终点多车型混载的定制公交线路规划方法 |
CN108053062A (zh) * | 2017-12-11 | 2018-05-18 | 北京航空航天大学 | 一种基于多源数据的定制公交线路生成方法 |
-
2018
- 2018-10-22 CN CN201811227668.4A patent/CN109344529B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106448233A (zh) * | 2016-08-19 | 2017-02-22 | 大连理工大学 | 基于大数据的公交线路时刻表协同优化方法 |
CN107194128A (zh) * | 2017-06-30 | 2017-09-22 | 东南大学 | 基于中心辐射型网络框架的多模式公共交通网络设计方法 |
CN107330559A (zh) * | 2017-07-03 | 2017-11-07 | 华南理工大学 | 一种多起终点多车型混载的定制公交线路规划方法 |
CN108053062A (zh) * | 2017-12-11 | 2018-05-18 | 北京航空航天大学 | 一种基于多源数据的定制公交线路生成方法 |
Non-Patent Citations (2)
Title |
---|
XUESONG ZHOU 等: "Finding optimal solutions for vehicle routing problem with pickup and delivery services with time windows: A dynamic programming approach based on state–space–time network representations", 《HTTPS://ARXIV.ORG/PDF/1507.02731.PDF》 * |
王健等: "考虑出行时间窗的定制公交线路车辆调度方法", 《中国公路学报》 * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109934401A (zh) * | 2019-03-09 | 2019-06-25 | 吉林大学 | 一种基于移动互联的定制公交方法 |
CN110598971A (zh) * | 2019-07-25 | 2019-12-20 | 中山大学 | 一种基于蚁群算法的响应式公交服务规划方法 |
CN110793535A (zh) * | 2019-08-22 | 2020-02-14 | 腾讯科技(深圳)有限公司 | 一种通勤线路规划方法、装置 |
CN111105078B (zh) * | 2019-11-29 | 2023-09-29 | 同济大学 | 一种定制公交网络优化方法 |
CN111105078A (zh) * | 2019-11-29 | 2020-05-05 | 同济大学 | 一种定制公交网络优化方法 |
CN111144618A (zh) * | 2019-12-04 | 2020-05-12 | 东南大学 | 一种基于两阶段优化模型的需求响应型定制公交网络规划方法 |
CN111832882A (zh) * | 2020-04-30 | 2020-10-27 | 北京嘀嘀无限科技发展有限公司 | 一种交通运输的控制方法、装置、存储介质和电子设备 |
CN112085249A (zh) * | 2020-07-27 | 2020-12-15 | 北京工业大学 | 一种基于强化学习的定制公交线路规划方法 |
CN112085249B (zh) * | 2020-07-27 | 2022-11-04 | 北京工业大学 | 一种基于强化学习的定制公交线路规划方法 |
CN112085349A (zh) * | 2020-08-19 | 2020-12-15 | 大连海事大学 | 一种基于乘客出行时间窗约束的需求响应公交调度方法 |
CN112419704A (zh) * | 2020-11-06 | 2021-02-26 | 杭州图软科技有限公司 | 一种基于大数据的公共交通路线规划方法及其系统 |
CN112419128A (zh) * | 2020-12-16 | 2021-02-26 | 武汉元光科技有限公司 | 一种线路规划方法以及相关设备 |
CN112419128B (zh) * | 2020-12-16 | 2024-03-05 | 武汉元光科技有限公司 | 一种线路规划方法以及相关设备 |
CN112734188A (zh) * | 2020-12-30 | 2021-04-30 | 杭州电子科技大学 | 一种基于两阶段混合元启发式算法的家庭医疗护理调度优化方法 |
CN112700850B (zh) * | 2020-12-30 | 2024-04-02 | 杭州电子科技大学 | 一种基于三阶段混合启发式算法的家庭医疗护理调度随机优化方法 |
CN112700850A (zh) * | 2020-12-30 | 2021-04-23 | 杭州电子科技大学 | 一种基于三阶段混合启发式算法的家庭医疗护理调度随机优化方法 |
CN112734188B (zh) * | 2020-12-30 | 2023-11-24 | 杭州电子科技大学 | 一种基于两阶段混合元启发式算法的家庭医疗护理调度优化方法 |
CN112750063A (zh) * | 2021-01-04 | 2021-05-04 | 李璐 | 基于随机规划的公交车队设施选址-路径规划-调度方法 |
CN112750063B (zh) * | 2021-01-04 | 2023-12-05 | 李璐 | 基于随机规划的公交车队设施选址-路径规划-调度方法 |
CN112466122B (zh) * | 2021-01-28 | 2021-06-15 | 深圳市城市交通规划设计研究中心股份有限公司 | 一种公交线网的备选线路集生成、线路规划方法及装置 |
CN112466122A (zh) * | 2021-01-28 | 2021-03-09 | 深圳市城市交通规划设计研究中心股份有限公司 | 一种公交线网的备选线路集生成、线路规划方法及装置 |
CN113096430A (zh) * | 2021-03-09 | 2021-07-09 | 东南大学 | 基于欧氏距离收费的弹性公交即时需求动态响应方法 |
CN113096430B (zh) * | 2021-03-09 | 2022-04-01 | 东南大学 | 基于欧氏距离收费的弹性公交即时需求动态响应方法 |
CN113538886A (zh) * | 2021-06-10 | 2021-10-22 | 广东工业大学 | 一种实时响应式定制公交递阶调度方法 |
CN113469451A (zh) * | 2021-07-19 | 2021-10-01 | 杭州数知梦科技有限公司 | 基于启发式算法的定制公交线路生成方法 |
CN113469451B (zh) * | 2021-07-19 | 2024-01-05 | 杭州数知梦科技有限公司 | 基于启发式算法的定制公交线路生成方法 |
CN114518763B (zh) * | 2022-04-20 | 2022-07-12 | 清华大学深圳国际研究生院 | 一种协同公交车的路径规划方法 |
CN114518763A (zh) * | 2022-04-20 | 2022-05-20 | 清华大学深圳国际研究生院 | 一种协同公交车的路径规划方法 |
CN114897445A (zh) * | 2022-07-12 | 2022-08-12 | 深圳市城市交通规划设计研究中心股份有限公司 | 公交线网的停靠点调整优化方法、装置及可读存储介质 |
CN115641704A (zh) * | 2022-12-26 | 2023-01-24 | 东风悦享科技有限公司 | 一种智能公交调度方法及系统 |
CN116090785A (zh) * | 2023-02-21 | 2023-05-09 | 哈尔滨工业大学 | 针对大型活动散场场景两阶段的定制公交规划方法 |
CN116090785B (zh) * | 2023-02-21 | 2024-03-08 | 哈尔滨工业大学 | 针对大型活动散场场景两阶段的定制公交规划方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109344529B (zh) | 2023-11-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109344529A (zh) | 一种基于两阶段启发式算法的定制公交线路设计方法 | |
CN109409599B (zh) | 基于实时需求的定制公交线路开通优化方法 | |
CN103198104B (zh) | 一种基于城市智能公交系统的公交站点od获取方法 | |
CN108053062B (zh) | 一种基于多源数据的定制公交线路生成方法 | |
CN104620296B (zh) | 按需车辆运行管理装置、按需车辆运行管理方法以及按需车辆运行管理系统 | |
CN101216913B (zh) | 合乘动态匹配多级筛选方法 | |
CN1969170B (zh) | 导航系统、路径搜索服务器及路径搜索方法 | |
TWI735292B (zh) | 一種為從源頭位置到目標位置的旅程提供路線的方法 | |
CN103440780B (zh) | 一种基于定位标签的城市复杂交通环境路径诱导系统和方法 | |
US11514796B2 (en) | System and method for determining and recommending vehicle pick-up location | |
CN111144618B (zh) | 一种两阶段需求响应型定制公交网络规划方法 | |
CN101187996B (zh) | 一种基于自适应骨干网的快速公交选线方法 | |
CN107944611B (zh) | 面向个性化出行需求的跨模式组合出行规划方法 | |
CN105070044A (zh) | 一种基于乘客预约的定制公交合乘车辆动态调度方法 | |
CN110222277B (zh) | 一种基于大数据分析的旅游信息推荐方法及装置 | |
CN106203725A (zh) | 基于启发式搜索的门到门出行路径方案个性化推荐方法 | |
CN107092986B (zh) | 基于站点及共线运行的公交乘客出行路径选择方法 | |
CN103530694A (zh) | 基于时空网络构建的城市地铁动态客流分配方法 | |
US9366541B2 (en) | AMI visual routing | |
Garcia et al. | Hybrid approach for the public transportation time dependent orienteering problem with time windows | |
CN107202587B (zh) | 规划出行路线的方法和装置 | |
JP7227862B2 (ja) | ポート開設場所提案装置 | |
CN115660157A (zh) | 智慧城市充电桩建造规划方法和物联网系统、装置及介质 | |
Madleňák et al. | The optimization approach of postal transportation network based on uncapacitated fixed charge location model in conditions of Slovak Republic | |
Dalkılıç et al. | A gradual approach for multimodel journey planning: A case study in Izmir, Turkey |
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 |