CN106503869A - 一种基于站点短时需求预测的公共自行车动态调度方法 - Google Patents
一种基于站点短时需求预测的公共自行车动态调度方法 Download PDFInfo
- Publication number
- CN106503869A CN106503869A CN201611023798.7A CN201611023798A CN106503869A CN 106503869 A CN106503869 A CN 106503869A CN 201611023798 A CN201611023798 A CN 201611023798A CN 106503869 A CN106503869 A CN 106503869A
- Authority
- CN
- China
- Prior art keywords
- website
- public bicycles
- time
- dispatching
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000003062 neural network model Methods 0.000 claims abstract description 5
- 239000010410 layer Substances 0.000 claims description 62
- 238000012549 training Methods 0.000 claims description 19
- 230000006870 function Effects 0.000 claims description 13
- 238000013528 artificial neural network Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 7
- 230000007935 neutral effect Effects 0.000 claims description 4
- 230000004913 activation Effects 0.000 claims description 3
- 238000005553 drilling Methods 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 239000011229 interlayer Substances 0.000 claims description 3
- 230000009897 systematic effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001537 neural effect Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- YBMRDBCBODYGJE-UHFFFAOYSA-N germanium dioxide Chemical compound O=[Ge]=O YBMRDBCBODYGJE-UHFFFAOYSA-N 0.000 description 2
- 210000005036 nerve Anatomy 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 230000001680 brushing effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000004089 microcirculation Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000002407 reforming Methods 0.000 description 1
- 230000004044 response 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
- G06Q30/00—Commerce
- G06Q30/06—Buying, selling or leasing transactions
- G06Q30/0645—Rental transactions; Leasing transactions
-
- 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)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Accounting & Taxation (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Finance (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (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:采集某一公共自行车站点一天内各时间段内公共自行车租赁的数量Rt和归还的数量Ht,t=1,2,…,T,t为时间段编号,T为划分的时间段个数,则各时间段末公共自行车站点空闲桩位数Kt为Kt=Kt-1-Ht+Rt,从而得到公共自行车站点空闲桩位数的时间序列K=﹛K1,K2,…,Kt…,KT﹜;采集公共自行车站点在上个星期相同时间段的空闲桩位数Klt,则上个星期站点空闲桩位数的时间序列为Kl=﹛Kl1,Kl2,…,Klt…,KlT﹜;
步骤2:计算选取的第t个时间段末站点空闲桩位数Kt与Kt-4、Kt-3、Kt-2、Kt-1、Kl(t-2)、Kl(t-1)、Klt、Kl(t+1)、Kl(t+2)之间的斯皮尔曼相关系数,取与Kt间的相关系数大于要求阈值的所有时间段作为输入向量,从而确定输入向量的维数M,所述要求阈值在0.8~1之间;
步骤3:建立BP神经网络模型,设定层与层间的权值、初始阈值、最大学习次数,并提供训练样本和测试样板,所述训练样本和测试样板均包括从步骤2中挑选出的M个空闲桩位数作为输入向量X=(X1,X2,…,Xt,…,XM)和输出向量Kt;
步骤4:通过BP神经网络的梯度训练法对步骤3提供的训练样本进行训练,然后用测试样本测试,直至BP神经网络的总误差小于期望误差值,网络训练完成,最后对未来某天的第t时间段末公共自行车站点空闲桩位数Kt(future)进行预测;
步骤5:根据步骤4预测的结果计算各站点调度时间窗和调度需求量范围,所述调度时间窗为[hmin,hmax],其中hmin为某一站点i在未来时间段t末时刻的可供租赁的公共自行车数量qi(t)与该站点库存能力Qi之比的最小值,hmax为站点i在未来时间段t末时刻的可供租赁的公共自行车数量qi(t)与该站点库存能力Qi之比的最大值,其中qi(t)=Qi-Kt(future);未来时间段t末时刻站点i的调度需求量di的范围则为qi(t)-hmax·Qi≤di≤qi(t)-hmin·Qi;
步骤6:根据各站点的调度需求,以调度成本最低和用户满意度最大为目标,计算调度方案,调度车辆执行调度计划。
进一步的,本发明方法中,步骤4中的通过BP神经网络梯度训练法对训练样本进行训练,包括如下步骤:
步骤41:确定输入层为M层,输出层为1层,取隐含层层数为N层,根据公式On=f(∑ωmn×Xm-θm)和Yt=f(∑ωnm×On-θm)逐层计算样本从输入层向隐含层的输出、隐含层向输出层的输出,其中On为输入层向隐含层的输出,n为隐含层的任意一层,n=1,2,…,N,Yt为隐含层向输出层的输出,ωmn为连接结点m到下一层结点n的权重,m为输入层的任意一层,m=1,2,…,M,Xm为网络初始化选取的任意一输入向量,θm为神经单元阙值,f(·)为激活函数
步骤42:根据公式和δ=Yt(1-Yt)(Yt-Kt)Om分别计算网络的总误差和误差信号δ;
步骤43:根据公式ω(n+1)=ω(n)+(-ηδ),按梯度对输出层到隐含层、隐含层到输入层的权重和阈值进行调整,其中η为学习速率;
步骤44:训练完全部的输入样本后,判断网络的总误差是否满足精度要求,如果已达到精度要求或已进行的学习次数到达预先设定的最大学习次数,则训练结束,否则开始新一轮学习。
进一步的,本发明方法中,所述步骤43中,学习速率η的取值范围在0.01~0.8之间,η通过选取多个不同值对神经网络进行训练后,根据观察系统误差的下降速率来确定。
进一步的,本发明方法中,步骤6中通过采用遗传算法,在约束条件下求解调度目标为运输成本最低和顾客满意度最大的目标函数,计算得到调度方案的方法;
所述目标函数为:其中为调度成本,k为需要调度的次数,cij表示调度车辆从站点i向站点j调度一次的总成本,cij=c0+Tij,其中c0为调度车辆启动所花费的固定成本,Tij为调度车辆从站点i向站点j调度的运输成本,xij是个二值函数,当调度车辆在站点i完成调度任务并驶向站点j时,xij取1,否则取0;
Si(t)i是用户满意度函数,表示为:
[Bi,Ci]是站点最佳调度时间范围,[Ai,Di]为可容忍调度时间范围;
λ为变换因子,λ=a*k*c0,a为[0,1]之间的常数,按照调度的侧重点进行设置;
所述约束条件如下:
Ⅰ.调度车辆上公共自行车的载运量始终位于零和车辆的最大载运能力之间,即:
式中di为公共自行车站点i的调度需求,该数值随时间不断变化,当其为正数时,表示站点i需要调出公共自行车;当其为负数时,表示站点i需要调入公共自行车;yi是个二值函数,当站点i完成调度任务时,yi取1,否则取0;Q为调度车辆的最大载运能力;
Ⅱ.调度开始时,调度车辆应从停车场出发,并最终返回初始停车场,即:
Ⅲ.调度车辆对任一站点进行调度后,必定会从该站点离开,即:
Ⅳ.到达站点j的时刻tj应等于到达上一个站点i的时刻加上在该站点的服务时间以及两个站点间调度车辆的行程时间,即:tj=ti+tui+tij,ti为车辆到达站点i并开始服务的时间,tui为车辆在站点i的服务时间,tij为调度车辆从站点i行驶到站点j的时间。
进一步的,本发明方法中,所述步骤2中的阈值为0.9。
本发明方法将公共自行车站点短时需求预测与动态调度相结合,通过信息更新和用户满意度的量化,实时调整调度计划,保障公共自行车调度的高效与准确性,提高公共自行车服务水平,增强公共自行车的吸引力。
有益效果:本发明与现有技术相比,具有以下优点:
传统方法一般通过尽可能多的罗列影响因素来对公共自行车需求进行预测,然而由于影响因素众多,难以罗列齐全且数据收集困难,使得预测的精度较低。如步骤1所示,本发明以最直接的站点空闲桩位数的时间序列作为因变量,数据仅从刷卡数据中就能获得,无需进行大规模的实地调研,并建立神经网络模型对未来的站点空闲桩位数进行预测,运算速度快,预测精确高。
本发明所提出的基于神经网络的公共自行车站点短时需求预测方法,能够预测未来短时间站点空闲桩位数,运算速度快,预测精确高。如步骤1中取T=144,则将1天24小时划分为144个长度为10分钟的时间段来作为站点空闲桩位数的时间序列,再通过步骤4的预测手段就能够预测未来10分钟内站点需求变化情况。可以真正反映站点实时的需求变化。相比传统方法只能得到每天的需求变化,本发明能够得到一天之内各个时段的需求变化,为站点的实时智能调度提供了依据。
传统的调度方法仅基于调配前各站点的调入调出需求量来进行车辆调度,缺乏用户满意度对系统影响的考虑,不能准确反应站点需求的变化。本发明在步骤5中得到的调度时间窗和调度需求量的基础上,还充分考虑步骤6中引入的调度成本最低和用户满意度最大的约束函数来制定车辆执行调度计划,能够同时满足公共自行车运营管理者和用户使用者的双重需求。
同时本发明本发明通过设定更新周期,在调度过程中实时调整调度计划,提高调度的效率与准确度,并能充分考虑调度延误对公共自行车系统的影响,相比传统方法更能满足市场的需求。
附图说明
图1为本发明方法的流程图;
图2为本发明方法中站点短时需求预测的预测值与实际值效果对比图。
具体实施方式
结合附图和实施例,对本发明技术方案详细说明如下:
本发明是一种基于短时需求预测的公共自行车动态调度方法,具体预测过程如下:
1)获取站点空闲桩位时间序列
以宁波市某一公共自行车站点为例,以10min为间隔,采集2014年5月12日—5月23日两周工作日的站点公共自行车租赁的数量Rt和归还的数量Ht,t=1,2,…,T,T=720,则各时间段末公共自行车站点空闲桩位数Kt可表示为Kt=Kt-1-Ht+Rt,从而得到5月19日—5月23日公共自行车站点空闲桩位数的时间序列K=﹛K1,K2,…,K720﹜;用Klt表示公共自行车站点在上个星期5月12日-5月16日相同时间段的空闲桩位数,则该星期站点空闲桩位数的时间序列可表示为Kl=﹛Kl1,Kl2,…,Kl720﹜。
2)计算各时间段spearman相关系数
计算Kt时间段与Kt-4、Kt-3、Kt-2、Kt-1以及Kl(t-2)、Kl(t-1)、Klt、Kl(t+1)、Kl(t+2)之间的斯皮尔曼相关系数如表1所示,Kt间的相关系数大于要求阈值的所有时间段作为输入向量,本发明取该阈值为0.9,斯皮尔曼相关系数大于0.9,分别为Kt-1、Kt-2、Kt-3、Kl(t-1)、Klt、Kl(t+1)和Kl(t+2),从而确定输入向量的维数M(本例中M=7)。
表1相邻公交停靠站之间交叉口数量
3)网络初始化
建立BP神经网络模型,随机分配层与层间的初始权值和阈值,并设定最大学习次数、期望误差等参数;并提供训练样本,每个样本需要提供输入向量X=(X1,X2,…,XM),和输出向量Kt;
4)网络训练和预测
根据公式On=f(∑ωmn×Xm-θm)和Yt=f(∑ωnm×On-θm)逐层计算样本从输入层向隐含层的输出、隐含层向输出层的输出,其中On为输入层向隐含层的输出,n为隐含层的任意一层,n=1,2,…,N,Yt为隐含层向输出层的输出,ωmn为连接结点m到下一层结点n的权重,m为输入层的任意一层,m=1,2,…,M,Xm为网络初始化选取的任意一输入向量,θm为神经单元阙值,f(·)为激活函数
基于计算的输出值和实际值之间的误差,计算系统总误差,公式为:
计算误差信号,公式为δ=Yt(1-Yt)(Yt-Kt)Om;
为修正误差值,使总误差E按误差信号δ梯度下降,应按梯度对输出层到隐含层、隐含层到输入层的权重和阈值进行调整,假设神经元j到神经元k的修正值用Δjk表示,则Δjk=-ηδ,其中η称为学习速率取值范围在0.01~0.8之间,选取多个不同的学习速率,对神经网络进行训练,若网络的总误差E和误差信号δ呈现出下降趋势,则说明选取的学习速率合适,若误差出现震荡现象,则说明选取的学习速率过大,需重新选取。在选取多个不同的值对神经网络进行训练后,本发明中取0.01;
设定相关参数,输入样本进行训练,判断网络的总误差是否满足精度要求,如果已达到精度要求或已进行的学习次数到达预先设定的最大学习次数,训练结束,否则开始新一轮学习。
以周二到周四之间的所有数据为样本用于训练BP神经网络,BP神经网络训练完成后,为了预测周五第t时间段的站点空闲桩位数,重新选取周五第t时间段的输入向量X(future),输入训练好的BP神经网络中,此时输出的就是周五第t时间段末公共自行车站点空闲桩位数Kt(future),得到预测结果如图2所示。对预测结果进行误差分析可知,90.2%的预测误差在5%以内,表明本发明所建立的预测方法可进行公共自行车站点短时需求预测,有良好的精度。
5)计算调度时间窗和调度需求量
在本例中选取的调度区域内共有17个公共自行车站点编号分别1-17,因此根据需求预测调度区域内17个站点周五各个时段的空闲桩位数结果,结果显示在早高峰时段,2号、9号和10号站点首先出现调度需求,计算需要调度站点的调度时间窗和调度需求量如表1所示:
表1站点调度需求情况
编号 | 库存能力(辆) | 最佳调配时间窗 | 可容忍服务时间窗 | 调度需求(辆) |
0 | — | — | — | — |
2 | 20 | 7:53-7:55 | 7:53-8:01 | 10 |
9 | 25 | 7:46-8:00 | 7:46-8:03 | 5 |
10 | 25 | 7:52-7:57 | 7:48-8:00 | 15 |
6)求解调度方案
根据站点的调度需求,以调度成本最低和用户满意度最大为目标。
所述目标函数为:其中为调度成本,cij表示调度车辆从站点i向站点j调度一次的总成本,cij=c0+Tij,其中c0为调度车辆启动所花费的固定成本,在本例中宁波调度车辆启动花费时间的固定成本为40元。Tij为调度车辆从站点i向站点j调度的运输成本,xij是个二值函数,当调度车辆在站点i完成调度任务并驶向站点j时,xij取1,否则取0;
Si(t)i是用户满意度函数,表示为:
[Bi,Ci]是站点最佳调度时间范围,[Ai,Di]为可容忍调度时间范围;
λ为变换因子,λ=a*k*c0,a为[0,1]之间的常数,按照调度的侧重点进行设置,当较为看重运输成本时,a取较小值;当更为看重满意度时,a取较大值;本例中为保证顾客满意度为100%,因此取a=1.
所述约束条件如下:
Ⅰ.调度车辆上公共自行车的载运量始终位于零和车辆的最大载运能力之间,即:
式中di为公共自行车站点i的调度需求,该数值随时间不断变化,当其为正数时,表示站点i需要调出公共自行车;当其为负数时,表示站点i需要调入公共自行车;yi是个二值函数,当站点i完成调度任务时,yi取1,否则取0;Q为调度车辆的最大载运能力,本例中Q取30辆;
Ⅱ.调度开始时,调度车辆应从停车场出发,并最终返回初始停车场,即:
Ⅲ.调度车辆对任一站点进行调度后,必定会从该站点离开,即:
Ⅳ.到达站点j的时刻tj应等于到达上一个站点i的时刻加上在该站点的服务时间以及两个站点间调度车辆的行程时间,即:tj=ti+tui+tij,ti为车辆到达站点i并开始服务的时间,tui为车辆在站点i的服务时间,tij为调度车辆从站点i行驶到站点j的时间。
通过遗传算法求解,得到调度方案为:0-9-10-2,调度车辆执行调度计划。
7)信息更新
对所有信息进行第一次更新后,得到站点需求如下表2所示,此时车辆已完成站点9和站点10的调度,按照同样的方法计算得到新的调度方案为:9-10-2-12。
表2第一次信息更新后站点调度需求情况
8)检查调度进程
在下一个周期继续进行信息更新,车辆已完成站点2的调度,正在前往站点12的途中,站点需求产生变化如表3所示,计算得到新的调度方案为:9-10-2-12-1-3,调度车辆执行新的调度计划。
表3第二次信息更新后站点调度需求情况
编号 | 库存能力(辆) | 最佳调配时间窗 | 可容忍服务时间窗 | 调度需求(辆) |
0 | — | — | — | — |
1 | 20 | 8:03-8:10 | 8:03-8:16 | 15 |
2 | 20 | 7:53-7:55 | 7:53-8:01 | 10 |
3 | 25 | 8:05-8:13 | 8:05-8:16 | -20 |
9 | 25 | 7:46-8:00 | 7:46-8:03 | 5 |
10 | 25 | 7:52-7:57 | 7:48-8:00 | 15 |
12 | 25 | 7:58-8:05 | 7:58-8:08 | -20 |
按照相同的周期进行信息更新,得到站点调度需求如表4所示,可得最终调度方案为9-10-2-12-1-7-3-6-4-15,当调度车辆在站点15完成调度,即可认为调度结束。
表4最终站点调度需求情况
编号 | 库存能力(辆) | 最佳调配时间窗 | 可容忍服务时间窗 | 调度需求(辆) |
0 | — | — | — | — |
1 | 20 | 8:03-8:10 | 8:03-8:16 | 15 |
2 | 20 | 7:53-7:55 | 7:53-8:01 | 10 |
3 | 25 | 8:05-8:13 | 8:05-8:16 | -20 |
4 | 30 | 8:15-8:21 | 8:15-8:23 | -15 |
6 | 50 | 8:13-8:20 | 8:13-8:20 | 20 |
7 | 25 | 8:06-8:10 | 8:06-8:13 | 5 |
9 | 25 | 7:46-8:00 | 7:46-8:03 | 5 |
10 | 25 | 7:52-7:57 | 7:48-8:00 | 15 |
12 | 25 | 7:58-8:05 | 7:58-8:08 | -20 |
15 | 25 | 8:15-8:26 | 8:15-8:29 | 8 |
上述实施例仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和等同替换,这些对本发明权利要求进行改进和等同替换后的技术方案,均落入本发明的保护范围。
Claims (5)
1.一种基于站点短时需求预测的公共自行车动态调度方法,其特征在于,该方法包括如下步骤:
步骤1:采集某一公共自行车站点一天内各时间段内公共自行车租赁的数量Rt和归还的数量Ht,t=1,2,…,T,t为时间段编号,T为划分的时间段个数,则各时间段末公共自行车站点空闲桩位数Kt为Kt=Kt-1-Ht+Rt,从而得到公共自行车站点空闲桩位数的时间序列K=﹛K1,K2,…,Kt…,KT﹜;采集公共自行车站点在上个星期相同时间段的空闲桩位数Klt,则上个星期站点空闲桩位数的时间序列为Kl=﹛Kl1,Kl2,…,Klt…,KlT﹜;
步骤2:计算选取的第t个时间段末站点空闲桩位数Kt与Kt-4、Kt-3、Kt-2、Kt-1、Kl(t-2)、Kl(t-1)、Klt、Kl(t+1)、Kl(t+2)之间的斯皮尔曼相关系数,取与Kt间的相关系数大于要求阈值的所有时间段作为输入向量,从而确定输入向量的维数M,所述要求阈值在0.8~1之间;
步骤3:建立BP神经网络模型,设定层与层间的权值、初始阈值、最大学习次数,并提供训练样本和测试样板,所述训练样本和测试样板均包括从步骤2中挑选出的M个空闲桩位数作为输入向量X=(X1,X2,…,Xt,…,XM)和输出向量Kt;
步骤4:通过BP神经网络的梯度训练法对步骤3提供的训练样本进行训练,然后用测试样本测试,直至BP神经网络的总误差小于期望误差值,网络训练完成,最后对未来某天的第t时间段末公共自行车站点空闲桩位数Kt(future)进行预测;
步骤5:根据步骤4预测的结果计算各站点调度时间窗和调度需求量范围,所述调度时间窗为[hmin,hmax],其中hmin为某一站点i在未来时间段t末时刻的可供租赁的公共自行车数量qi(t)与该站点库存能力Qi之比的最小值,hmax为站点i在未来时间段t末时刻的可供租赁的公共自行车数量qi(t)与该站点库存能力Qi之比的最大值,其中qi(t)=Qi-Kt(future);未来时间段t末时刻站点i的调度需求量di的范围则为qi(t)-hmax·Qi≤di≤qi(t)-hmin·Qi;
步骤6:根据各站点的调度需求,以调度成本最低和用户满意度最大为目标,计算调度方案,调度车辆执行调度计划。
2.根据权利要求1所述的一种基于站点短时需求预测的公共自行车动态调度方法,其特征在于,所述步骤4中的通过BP神经网络梯度训练法对训练样本进行训练,包括如下步骤:
步骤41:确定输入层为M层,输出层为1层,取隐含层层数为N层,根据公式On=f(∑ωmn×Xm-θm)和Yt=f(∑ωnm×On-θm)逐层计算样本从输入层向隐含层的输出、隐含层向输出层的输出,其中On为输入层向隐含层的输出,n为隐含层的任意一层,n=1,2,…,N,Yt为隐含层向输出层的输出,ωmn为连接结点m到下一层结点n的权重,m为输入层的任意一层,m=1,2,…,M,Xm为网络初始化选取的任意一输入向量,θm为神经单元阙值,f(·)为激活函数
步骤42:根据公式和δ=Yt(1-Yt)(Yt-Kt)Om分别计算网络的总误差和误差信号δ;
步骤43:根据公式ω(n+1)=ω(n)+(-ηδ),按梯度对输出层到隐含层、隐含层到输入层的权重和阈值进行调整,其中η为学习速率;
步骤44:训练完全部的输入样本后,判断网络的总误差是否满足精度要求,如果已达到精度要求或已进行的学习次数到达预先设定的最大学习次数,则训练结束,否则开始新一轮学习。
3.根据权利要求2所述的一种基于站点短时需求预测的公共自行车动态调度方法,其特征在于,所述步骤43中,学习速率η的取值范围在0.01~0.8之间,η通过选取多个不同值对神经网络进行训练后,根据观察系统误差的下降速率来确定。
4.根据权利要求1、2或3所述的一种基于站点短时需求预测的公共自行车动态调度方法,其特征在于,所述步骤6中通过采用遗传算法,在约束条件下求解调度目标为运输成本最低和顾客满意度最大的目标函数,计算得到调度方案的方法;
所述目标函数为:其中为调度成本,k为需要调度的次数,cij表示调度车辆从站点i向站点j调度一次的总成本,cij=c0+Tij,其中c0为调度车辆启动所花费的固定成本,Tij为调度车辆从站点i向站点j调度的运输成本,xij是个二值函数,当调度车辆在站点i完成调度任务并驶向站点j时,xij取1,否则取0;
Si(t)i是用户满意度函数,表示为:
[Bi,Ci]是站点最佳调度时间范围,[Ai,Di]为可容忍调度时间范围;
λ为变换因子,λ=a*k*c0,a为[0,1]之间的常数,按照调度的侧重点进行设置;
所述约束条件如下:
Ⅰ.调度车辆上公共自行车的载运量始终位于零和车辆的最大载运能力之间,即:
式中di为公共自行车站点i的调度需求,该数值随时间不断变化,当其为正数时,表示站点i需要调出公共自行车;当其为负数时,表示站点i需要调入公共自行车;yi是个二值函数,当站点i完成调度任务时,yi取1,否则取0;Q为调度车辆的最大载运能力;
Ⅱ.调度开始时,调度车辆应从停车场出发,并最终返回初始停车场,即:
Ⅲ.调度车辆对任一站点进行调度后,必定会从该站点离开,即:
Ⅳ.到达站点j的时刻tj应等于到达上一个站点i的时刻加上在该站点的服务时间以及两个站点间调度车辆的行程时间,即:tj=ti+tui+tij,ti为车辆到达站点i并开始服务的时间,tui为车辆在站点i的服务时间,tij为调度车辆从站点i行驶到站点j的时间。
5.根据权利要求1、2或3所述的一种基于站点短时需求预测的公共自行车动态调度方法,其特征在于,所述步骤2中的阈值为0.9。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611023798.7A CN106503869A (zh) | 2016-11-14 | 2016-11-14 | 一种基于站点短时需求预测的公共自行车动态调度方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611023798.7A CN106503869A (zh) | 2016-11-14 | 2016-11-14 | 一种基于站点短时需求预测的公共自行车动态调度方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106503869A true CN106503869A (zh) | 2017-03-15 |
Family
ID=58327270
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611023798.7A Pending CN106503869A (zh) | 2016-11-14 | 2016-11-14 | 一种基于站点短时需求预测的公共自行车动态调度方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106503869A (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107103388A (zh) * | 2017-04-26 | 2017-08-29 | 杨毅 | 基于需求预测的机器人调度系统与方法 |
CN107145714A (zh) * | 2017-04-07 | 2017-09-08 | 浙江大学城市学院 | 基于多因素的公共自行车使用量预测方法 |
CN107239883A (zh) * | 2017-05-12 | 2017-10-10 | 同济大学 | 一种汽车共享系统车辆的调度方法 |
CN107609698A (zh) * | 2017-09-06 | 2018-01-19 | 上海享骑电动车服务有限公司 | 虚拟桩 |
CN107657353A (zh) * | 2017-11-09 | 2018-02-02 | 东峡大通(北京)管理咨询有限公司 | 租赁物的调度方法及系统 |
CN107766994A (zh) * | 2017-12-04 | 2018-03-06 | 长沙理工大学 | 一种共享自行车调度方法与调度系统 |
CN107886159A (zh) * | 2017-08-28 | 2018-04-06 | 上海享骑电动车服务有限公司 | 虚拟站点 |
CN108764606A (zh) * | 2018-04-08 | 2018-11-06 | 浙江广播电视大学 | 基于动态调度时间域的共享自行车系统调度方法 |
CN109146129A (zh) * | 2018-07-03 | 2019-01-04 | 蔚来汽车有限公司 | 充电桩群推荐方法、可用桩预测、结果获取方法及控制器 |
CN109829649A (zh) * | 2019-01-31 | 2019-05-31 | 北京首汽智行科技有限公司 | 一种车辆调度方法 |
CN110059839A (zh) * | 2019-04-25 | 2019-07-26 | 北京首汽智行科技有限公司 | 一种基于需求的提取未满足需求的方法 |
CN110444008A (zh) * | 2019-08-22 | 2019-11-12 | 欧拉信息服务有限公司 | 一种车辆调度方法及装置 |
CN110866591A (zh) * | 2019-10-28 | 2020-03-06 | 浙江大学 | 基于需求预测进行前瞻式云制造服务租赁配置的方法 |
CN110930687A (zh) * | 2017-12-18 | 2020-03-27 | 北京物资学院 | 高峰用车时间段车辆调度的方法及系统 |
CN111325425A (zh) * | 2018-12-14 | 2020-06-23 | 顺丰科技有限公司 | 调拨方案计算方法和装置 |
CN107045673B (zh) * | 2017-03-31 | 2020-09-29 | 杭州电子科技大学 | 基于堆模型融合的公共自行车流量变化量预测方法 |
CN111815095A (zh) * | 2019-06-20 | 2020-10-23 | 北京嘀嘀无限科技发展有限公司 | 共享汽车调度方法、电子设备及存储介质 |
CN111915209A (zh) * | 2020-08-11 | 2020-11-10 | 上海钧正网络科技有限公司 | 共享交通工具的调度方法、装置、计算机设备和存储介质 |
CN113610455A (zh) * | 2021-07-06 | 2021-11-05 | 中科云谷科技有限公司 | 用于调拨工程机械配件的系统及方法 |
-
2016
- 2016-11-14 CN CN201611023798.7A patent/CN106503869A/zh active Pending
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107045673B (zh) * | 2017-03-31 | 2020-09-29 | 杭州电子科技大学 | 基于堆模型融合的公共自行车流量变化量预测方法 |
CN107145714A (zh) * | 2017-04-07 | 2017-09-08 | 浙江大学城市学院 | 基于多因素的公共自行车使用量预测方法 |
CN107145714B (zh) * | 2017-04-07 | 2020-05-22 | 浙江大学城市学院 | 基于多因素的公共自行车使用量预测方法 |
CN107103388A (zh) * | 2017-04-26 | 2017-08-29 | 杨毅 | 基于需求预测的机器人调度系统与方法 |
CN107239883A (zh) * | 2017-05-12 | 2017-10-10 | 同济大学 | 一种汽车共享系统车辆的调度方法 |
CN107886159A (zh) * | 2017-08-28 | 2018-04-06 | 上海享骑电动车服务有限公司 | 虚拟站点 |
CN107609698A (zh) * | 2017-09-06 | 2018-01-19 | 上海享骑电动车服务有限公司 | 虚拟桩 |
CN107657353A (zh) * | 2017-11-09 | 2018-02-02 | 东峡大通(北京)管理咨询有限公司 | 租赁物的调度方法及系统 |
CN107766994A (zh) * | 2017-12-04 | 2018-03-06 | 长沙理工大学 | 一种共享自行车调度方法与调度系统 |
CN110930687B (zh) * | 2017-12-18 | 2020-12-25 | 北京物资学院 | 高峰用车时间段车辆调度的方法及系统 |
CN110930687A (zh) * | 2017-12-18 | 2020-03-27 | 北京物资学院 | 高峰用车时间段车辆调度的方法及系统 |
CN108764606A (zh) * | 2018-04-08 | 2018-11-06 | 浙江广播电视大学 | 基于动态调度时间域的共享自行车系统调度方法 |
CN108764606B (zh) * | 2018-04-08 | 2020-07-17 | 浙江广播电视大学 | 基于动态调度时间域的共享自行车系统调度方法 |
CN109146129A (zh) * | 2018-07-03 | 2019-01-04 | 蔚来汽车有限公司 | 充电桩群推荐方法、可用桩预测、结果获取方法及控制器 |
CN111325425A (zh) * | 2018-12-14 | 2020-06-23 | 顺丰科技有限公司 | 调拨方案计算方法和装置 |
CN109829649A (zh) * | 2019-01-31 | 2019-05-31 | 北京首汽智行科技有限公司 | 一种车辆调度方法 |
CN110059839A (zh) * | 2019-04-25 | 2019-07-26 | 北京首汽智行科技有限公司 | 一种基于需求的提取未满足需求的方法 |
CN111815095A (zh) * | 2019-06-20 | 2020-10-23 | 北京嘀嘀无限科技发展有限公司 | 共享汽车调度方法、电子设备及存储介质 |
CN110444008A (zh) * | 2019-08-22 | 2019-11-12 | 欧拉信息服务有限公司 | 一种车辆调度方法及装置 |
CN110444008B (zh) * | 2019-08-22 | 2021-06-04 | 欧拉信息服务有限公司 | 一种车辆调度方法及装置 |
CN110866591A (zh) * | 2019-10-28 | 2020-03-06 | 浙江大学 | 基于需求预测进行前瞻式云制造服务租赁配置的方法 |
CN111915209A (zh) * | 2020-08-11 | 2020-11-10 | 上海钧正网络科技有限公司 | 共享交通工具的调度方法、装置、计算机设备和存储介质 |
CN113610455A (zh) * | 2021-07-06 | 2021-11-05 | 中科云谷科技有限公司 | 用于调拨工程机械配件的系统及方法 |
CN113610455B (zh) * | 2021-07-06 | 2024-05-07 | 中科云谷科技有限公司 | 用于调拨工程机械配件的系统及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106503869A (zh) | 一种基于站点短时需求预测的公共自行车动态调度方法 | |
CN103632212B (zh) | 一种时变用户均衡动态网络演化客流预测系统和方法 | |
CN109034566B (zh) | 一种基于公交站点上下客流的智能调度方法及装置 | |
CN102044149B (zh) | 一种基于时变客流的城市公交运营协调方法与装置 | |
CN110555990B (zh) | 一种基于lstm神经网络的有效停车时空资源预测方法 | |
CN104636828B (zh) | 基于马尔科夫链的公共自行车站点供需预测方法 | |
Huang | Pricing and logit-based mode choice models of a transit and highway system with elastic demand | |
CN109886491A (zh) | 一种海量运行数据驱动的地面公交两阶段优化调度方法 | |
Chen et al. | Uncertainty in urban mobility: Predicting waiting times for shared bicycles and parking lots | |
CN102324128A (zh) | 基于ic卡记录的公交站间od客流预测方法及装置 | |
CN106127357A (zh) | 一种基于预约数据的定制公交自动布线系统及方法 | |
US10635916B2 (en) | Determining vehicle crowdedness using real-time location data | |
CN110991607A (zh) | 一种地铁客流预测方法、装置、电子设备及存储介质 | |
Tsiaras et al. | parkitsmart: Minimization of cruising for parking | |
CN110570656A (zh) | 定制公共交通线路方法及装置 | |
CN114418606B (zh) | 基于时空卷积网络的网约车订单需求预测方法 | |
Kamel et al. | A modelling platform for optimizing time-dependent transit fares in large-scale multimodal networks | |
CN113821547B (zh) | 快速高效的停车场占有率短时预测方法、系统及存储介质 | |
CN107316096A (zh) | 一种轨道交通一票通乘客进站量预测方法 | |
CN106355541A (zh) | 一种基于gps数据的出租车打车难易度识别方法 | |
CN109447352B (zh) | 一种公交出行od矩阵反推的组合方法 | |
CN106204404A (zh) | 一种公交运力的调整方法和系统 | |
CN107331156A (zh) | 城市交通综合调查云平台与系统 | |
Sang et al. | Research on optimization of customized bus routes based on uncertainty theory | |
CN114925911A (zh) | 一种基于无人驾驶公交车精准客流预测的自适应动态排班方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170315 |
|
RJ01 | Rejection of invention patent application after publication |