CN109325625B - 一种基于二元高斯非齐次泊松过程的单车数量预测方法 - Google Patents

一种基于二元高斯非齐次泊松过程的单车数量预测方法 Download PDF

Info

Publication number
CN109325625B
CN109325625B CN201811141977.XA CN201811141977A CN109325625B CN 109325625 B CN109325625 B CN 109325625B CN 201811141977 A CN201811141977 A CN 201811141977A CN 109325625 B CN109325625 B CN 109325625B
Authority
CN
China
Prior art keywords
bicycle
user
time
riding
returning
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.)
Active
Application number
CN201811141977.XA
Other languages
English (en)
Other versions
CN109325625A (zh
Inventor
乔少杰
黄飞虎
韩楠
彭京
魏军林
温敏
冉先进
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chengdu University of Information Technology
Original Assignee
Chengdu University of Information Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chengdu University of Information Technology filed Critical Chengdu University of Information Technology
Priority to CN201811141977.XA priority Critical patent/CN109325625B/zh
Publication of CN109325625A publication Critical patent/CN109325625A/zh
Application granted granted Critical
Publication of CN109325625B publication Critical patent/CN109325625B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Data Mining & Analysis (AREA)
  • Strategic Management (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Operations Research (AREA)
  • Human Resources & Organizations (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Economics (AREA)
  • Mathematical Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Evolutionary Biology (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Game Theory and Decision Science (AREA)
  • Algebra (AREA)
  • Development Economics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Traffic Control Systems (AREA)

Abstract

本发明公开一种基于二元高斯非齐次泊松过程的单车数量预测方法，属于数据挖掘预测技术领域，根据单车站点历史数据，利用泊松理论，建立非齐次泊松模型，并考虑环境参数的影响，建立最终的预测模型，对时间序列进行模拟并对单车数量进行修正，只需要单车站点历史的单车使用数据就能进行预测，不仅能够在相应数据获取受限的情况下，实现站点单车数量的预测，且可减小天气等不实时更新数据对预测结果的影响，显著提高站点单车在未来的一段时间内的数量预测的准确性和及时性。

Description

一种基于二元高斯非齐次泊松过程的单车数量预测方法
技术领域
本发明涉及数据挖掘预测技术领域，特别是涉及一种基于二元高斯非齐次泊松过程的单车数量预测方法。
背景技术
共享单车作为一种新型的交通工具，近年得到了快速的发展。在使用的过程中，单车经常出现分布不均匀的现象，为了方便单车系统管理员的管理，则需要动态地对单车进行调度。而对单车进行调度，则需要提前对站点的单车数量进行预测。由于单车使用受到很多不确定因素的影响，预测单车数量也是一个具有挑战性的任务。
现有的技术中，常用的方法是挖掘影响单车使用的因素，然后结合一些数据挖掘方法(比如回归模型，决策树，神经网络，支持向量机等)，训练预测模型，从而达到预测目的。由于影响单车使用的因素较多，且相关数据的获取受到各种限制，通常研究者能够得到的数据有限，导致并不能很全面地研究影响单车使用的因素。同时，天气数据也不是实时地更新，这对于实施性要求高的单车调度任务来说，会产生延时。导致对站点单车的数量预测准确性低，参考价值不高，对单车的调度和用户对单车的使用提供的帮助有限。
发明内容
为了解决上述问题，本发明提供一种基于二元高斯非齐次泊松过程的单车数量预测方法，只需要单车站点历史的单车使用数据就能进行预测，不仅能够在相应数据获取受限的情况下，实现站点单车数量的预测，且可减小天气等不实时更新数据对预测结果的影响，从而显著提高单车数量预测的准确性。为此，本发明采用的技术方案是：
提供一种基于二元高斯非齐次泊松过程的单车数量预测方法，该方法包括以下步骤：
S1、从单车管理系统中提取单车站点预测当日之前的历史数据，所述历史数据包括过去10-15天的用户骑走单车的时间序列和用户归还单车的时间序列；
预测当日的当前时刻用t表示，未来待预测时刻用t+s表示，t取值范围为00:00:00点到23:59:59点；将用户到达并归还单车的过程和用户到达并骑走单车的过程看作泊松过程，对其进行强度函数为μi(x)和λi(x)的非齐次泊松过程建模，并利用该模型对时间变量x求导，得到用户骑走单车和归还单车的概率密度函数；其中，i表示单车站点编号；
S2、根据步骤S1中的历史数据，用二元高斯函数拟合得到μi(x)和λi(x)的表达式，并结合概率密度函数，采用迭代法学习得到μi(x)和λi(x)的参数值，得到参数赋值后的强度函数，同时将所述参数值代入步骤S1得到的概率密度函数；
S3、根据所述单车站点在当前时刻之前的1-3小时内用户骑走单车的时间序列和用户归还单车的时间序列数据，计算环境参数值，并用环境参数值对步骤S2得到的参数赋值后的强度函数进行修正，同时实现对概率密度函数的修正，得到修正后的概率密度函数；
S4、利用S3得到修正后的概率密度函数，分别模拟所述单车站点中用户骑走单车的时间序列和用户归还单车的时间序列，所述模拟的次数为5-10次；
S5、在[t,t+s]时间段，分别计算步骤S4中模拟的用户骑走单车的时间序列和用户归还单车的时间序列中事件的平均次数，所述事件是指用户骑走单车和用户归还单车事件，0<s≤2.5小时；
S6、计算单车数量的预测值，具体计算公式为：
Xbike(t+s)＝Xbike(t)+Pbike-Pdock
其中，Xbike(t+s)表示在t+s时刻的单车数量，Xbike(t)表示在t时刻的单车数量，Pbike表示步骤S5用户骑走单车事件的平均次数，Pdock表示步骤S5用户归还单车事件的平均次数；
S7、对步骤S6所得的预测值进行修正，其修正方法为：当所述预测值大于所述单车站点的单车容量，则令所述预测值等于所述单车站点的单车容量；当所述预测值小于零，则令所述预测值等于零；本领域技术人员很容易理解，由于前述过程中，根据泊松理论进行模拟时，没有设定事件发生次数的上下限，因此，步骤S6中计算得到的单车数量就可能会大于该站点单车容量，或小于零，因此，根据实际情况，用上述修正方法进行修正后，可保证预测方法的可行性和准确性。
S8、根据步骤S6和S7，计算得到所述单车站点在t+s时刻的单车数量。
进一步的是，步骤S1中，用户骑走单车的概率密度函数表达式为：
其中，表示单车站点i中，当用户在时刻骑走一辆单车，那么下一刻在时刻骑走一辆单车的概率；μi(xn)和μi(x)为强度函数，分别表示xn时刻和x时刻用户骑走单车的数量；
用户归还单车的概率密度函数表达式为：
其中，表示单车站点i中，当用户在时刻骑走一辆单车，那么下一刻在时刻骑走一辆单车的概率；λi(xn)和λi(x)为强度函数，分别表示xn时刻和x时刻用户归还单车的数量。
进一步的是，所述步骤S2中，
其中，c1,c2,a1,a2,b1,b2为强度函数表达式中的参数。
进一步的是，在步骤S3中，所述环境参数值的计算公式为：
其中的计算公式为：
表示在时间范围[t-△t,t]内单车的使用数量，△t取值范围为1～3h；
则修正后的用户骑走单车的概率密度函数公式为：
则修正后的用户归还单车的概率密度函数公式为：
进一步的是，在步骤S4中，模拟用户骑走单车的时间序列的过程包括以下步骤：
S41、设初始时间t0为t，初始化用户骑走单车的时间序列变量Q；
S42、利用拒绝采样方法生成间隔时间τ，模拟用户骑走单车采用的公式为：
S43、令t0＝t0+τ，然后将t0存入Q；
S44、判断t0是否大于t+s，如果t0大于t+s，则结束模拟，否则，返回步骤S42；
模拟用户归还单车的时间序列的过程包括以下步骤：
S51、设初始时间t0为t，初始化用户骑走单车的时间序列变量Q′；
S52、利用拒绝采样方法生成间隔时间τ，模拟用户归还单车采用的公式为：
S53、令t0＝t0+τ，然后将t0存入Q′；
S54、判断t0是否大于t+s，如果t0大于t+s，则结束模拟，否则，返回步骤S42。
进一步的是，步骤S2中，所述迭代法为梯度下降法。
本发明方法的理论过程如下：
目前的共享单车系统，在每个单车站点，都有一定的停车位，也有一定的单车数量。每个用户会随机地到达该站点，要么骑走一辆单车，要么归还一辆单车。单车系统会实时地记录每个站点单车被骑走的时间和单车归还的时间。
通过数据分析，我们发现对每个单车站点而言，用户到达并骑走单车的时间序列和用户到达并归还单车的时间序列这两个序列均满足泊松过程。xn和xm分别表示时间序列的第n个点和第m个点。根据泊松过程，我们需要定义一个强度函数。对于用户到达并骑走单车的过程，其强度函数用μi(x)表示。对于用户到达并归还单车的过程，其强度函数用λi(x)表示。μi(x)和λi(x)均为时间x的函数。强度函数表示的是x时刻单车骑走或者归还的数量。在我们的方案中，x的取值范围是从00:00:00点到23:59:59点。因为通过研究发现，每个站点每天的单车骑走或者归还数量的变化可以用二元高斯函数进行拟合：
为了简化表达，对于用户骑走单车的过程和用户归还单车的两个过程，我们只选择用户骑走单车的过程进行介绍。因为这两个过程使用发方法是一样的。
根据泊松过程理论，在时间范围[v,v’)内，每个站点用户到达满足以下分布：
其中表示单车站点i，从0时刻到xn时刻为止，累计骑走的单车数量。表示站点i，从0时刻到xm时刻为止，累计骑走的单车数量。表示[xm,xn]时间内，发生了k次单车骑走事件。公式(2)中就表示[xm,xn]时间内，发生k次事件的概率。
对于事件:单车站点在时刻有一辆单车被骑走，紧接着在时刻又有一辆单车被骑走。
可以理解为，在和时刻之间没有事件发生，即k＝0，所以有：
对x求导得到其概率密度函数为：
用户到达并骑走单车的时间序列是由很多这样的事件序列组成，因此，对于每一天来说，事件的联合概率为：
如果用该单车站点W天的历史数据进行训练，则对于W天的时间序列来说，似然函数为:
其中L为每天事件发生的个数，通过迭代方法学习强度函数的参数，即c,a1,a2,b1,b2。
此外，由于强度函数估算的是第W+1天的单车使用量变化趋势，不能反应短时间内现实因素(环境因素)对单车使用量的影响，因此我们寻求结合环境参数对强度函数进行修正，以达到提高预测准确率的效果。
根据图2所示的计算方法示意图，可以利用概率密度函数和[t，t-△t]时间范围内单车估计全天单车使用量，用这个估计值作为的取值：
对式(7)进行线性拟合修正，得到环境参数：
对于式(8)的说明：
为根据当天[t-△t,t]时间段内用车数量得到的关于当天总共用车量的估计，在数值模拟中，理想的情况就是，我们希望模拟能产生个随机数，即当天用车量为然而，通过模拟发现，模拟程序产生的当天的事件个数为num_i.我们对num_i和 即要产生num_i个事件，对应的系数为所以，我们希望产生个事件，则系数需要等于：
利用式(8)得到的值对强度函数进行修正，得到用户骑走单车的概率密度函数：
再利用式(9)，进行算法模拟，得到用户骑走单车事件发生的次数。
同理可得，用户归还单车事件发生的次数。
再利用步骤S6-S7，则可得出单车站点在t+s时刻的单车数量。
采用本技术方案的有益效果：
第一，相对于现有技术，本发明方法利用非齐次泊松过程建模单车使用过程，利用学习训练方法进行参数训练，并利用环境参数进行修正，得到的预测模型显著提高站点单车在未来的一定时间(2小时)内的数量预测的准确性和及时性，且只需要单车站点历史的单车使用数据就能进行预测，即能够在相应数据获取受限的情况下，实现站点单车数量的准确预测，能满足适用范围更广。
第二，本发明方法考虑了环境参数的影响，并利用环境参数进行修正，显著提高了单车数量预测的准确性。
第三，本发明利用二元高斯函数建模非齐次泊松过程的强度函数，准确的反映了单车系统的实际使用情况。
附图说明
图1是本发明方法的流程图；
图2是本发明方法的环境参数推导原理图；
图3是在Bay Area单车系统上，本发明方法(BGIP)与对比方法的预测效果对比图；
图4是在Hubway单车系统上，本发明方法(BGIP)与对比方法的预测效果对比图；
图5是在Bay Area单车系统上，本发明方法(BGIP)与本发明方法(BGIP)未考虑环境参数时的预测效果对比图；
图6是在Hubway单车系统上，本发明方法(BGIP)与本发明方法(BGIP)未考虑环境参数时的预测效果对比图。
图
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚，下面结合附图和具体实施例对本发明作进一步阐述。
在本实施例中，如图1-6所示，一种基于二元高斯非齐次泊松过程的单车数量预测方法，该方法包括以下步骤：
S1、从单车管理系统中提取单车站点预测当日之前的历史数据，所述历史数据包括过去10-15天的用户骑走单车的时间序列和用户归还单车的时间序列；
预测当日的当前时刻用t表示，未来待预测时刻用t+s表示，t取值范围为00:00:00点到23:59:59点；将用户到达并归还单车的过程和用户到达并骑走单车的过程看作泊松过程，对其进行强度函数为μi(x)和λi(x)的非齐次泊松过程建模，并利用该模型对时间变量x求导，得到用户骑走单车和归还单车的概率密度函数；其中，i表示单车站点编号；
S2、根据步骤S1中的历史数据，用二元高斯函数拟合得到μi(x)和λi(x)的表达式，并结合概率密度函数，采用迭代法学习得到μi(x)和λi(x)的参数值，得到参数赋值后的强度函数，同时将所述参数值代入步骤S1得到的概率密度函数；
S3、根据所述单车站点在当前时刻之前的1-3小时内用户骑走单车的时间序列和用户归还单车的时间序列数据，计算环境参数值，并用环境参数值对步骤S2得到的参数赋值后的强度函数进行修正，同时实现对概率密度函数的修正，得到修正后的概率密度函数；
S4、利用S3得到修正后的概率密度函数，分别模拟所述单车站点中用户骑走单车的时间序列和用户归还单车的时间序列，所述模拟的次数为5-10次；
S5、在[t,t+s]时间段，分别计算步骤S4中模拟的用户骑走单车的时间序列和用户归还单车的时间序列中事件的平均次数，所述事件是指用户骑走单车事件和用户归还单车事件，0<s≤2.5小时；
S6、计算单车数量的预测值，具体计算公式为：
Xbike(t+s)＝Xbike(t)+Pbike-Pdock
其中，Xbike(t+s)表示在t+s时刻的单车数量，Xbike(t)表示在t时刻的单车数量，Pbike表示步骤S5用户骑走单车事件的平均次数，Pdock表示步骤S5用户归还单车事件的平均次数；
S7、对步骤S6所得的预测值进行修正，其修正方法为：当所述预测值大于所述单车站点的单车容量，则令所述预测值等于所述单车站点的单车容量；当所述预测值小于零，则令所述预测值等于零；
S8、根据步骤S6和S7，计算得到所述单车站点在t+s时刻的单车数量。
作为一种优化方案，步骤S1中，用户骑走单车的概率密度函数表达式为：
其中，表示单车站点i中，当用户在时刻骑走一辆单车，那么下一刻在时刻骑走一辆单车的概率；μi(xn)和μi(x)为强度函数，分别表示xn时刻和x时刻用户骑走单车的数量；
用户归还单车的概率密度函数表达式为：
其中，表示单车站点i中，当用户在时刻骑走一辆单车，那么下一刻在时刻骑走一辆单车的概率；λi(xn)和λi(x)为强度函数，分别表示xn时刻和x时刻用户归还单车的数量。
作为一种优化方案，所述步骤S2中，
其中，c1,c2,a1,a2,b1,b2为强度函数表达式中的参数。
作为一种优化方案，在步骤S3中，所述环境参数值的计算公式为：
其中的计算公式为：
表示在时间范围[t-△t,t]内单车的使用数量，△t取值范围为1～3h；
则修正后的用户骑走单车的概率密度函数公式为：
则修正后的用户归还单车的概率密度函数公式为：
作为一种优化方案，在步骤S4中，模拟用户骑走单车的时间序列的过程包括以下步骤：
S41、设初始时间t0为t，初始化用户骑走单车的时间序列变量Q；
S42、利用拒绝采样方法生成间隔时间τ，模拟用户骑走单车采用的公式为：
S43、令t0＝t0+τ，然后将t0存入Q；
S44、判断t0是否大于t+s，如果t0大于t+s，则结束模拟，否则，返回步骤S42；
模拟用户归还单车的时间序列的过程包括以下步骤：
S51、设初始时间t0为t，初始化用户骑走单车的时间序列变量Q′；
S52、利用拒绝采样方法生成间隔时间τ，模拟用户归还单车采用的公式为：
S53、令t0＝t0+τ，然后将t0存入Q′；
S54、判断t0是否大于t+s，如果t0大于t+s，则结束模拟，否则，返回步骤S42。
作为一种优化方案，步骤S2中，所述迭代法为梯度下降法。
下面以具体的实施例来做进一步的说明。
实施例1
本实施例的运行环境为：Matlab2016b,Java(JDK 1.7),Windows 10,CPU(Core i7(7500U)，8GB RAM。
在本实施例中，对于bay area单车系统，单车站点A(i＝10)，假设当前的日期为2013年8月23号，上午10点,单车数量为10，即Xbike(10点)＝10，预测单车站点A在未来的12点30分时刻的单车数量。
具体步骤包括：
S1、根据单车站点A的历史数据，从bay area单车管理系统中提取提取过去10-15天(即13号～22号，不包括23号)的单车骑行时间序列数据(用户骑走单车的时间序列和用户到达并归还单车的时间序列)，将用户到达并归还单车的过程和用户到达并骑走单车的过程看作泊松过程，对其进行强度函数为μi(x)和λi(x)的非齐次泊松过程建模，并利用该模型对时间变量x求导，得到用户骑走单车和归还单车的概率密度函数；
S2、根据公式(4)-(6)并结合时间序列和采用梯度下降法分别学习得到强度函数μi(x)和λi(x)的参数c1,c2,a1,a2,b1,b2，取值分别为c1＝0.1214,a1＝18,b1＝3.253,c2＝0.092,a2＝8.2,b2＝1.838；
S3、根据公式(7)和(8)计算环境参数ei(x)的值，采用的数据是当日(23号)时间范围7点至10点时间段内的单车数据，并用环境参数值对步骤S2得到的强度函数进行修正，进而实现对概率密度函数的修正；
S4、利用S3得到的概率密度函数，分别模拟所述单车站点A中用户骑走单车的时间序列和用户归还单车的时间序列，所述模拟的次数为5次；
S5、在10点至12点30分时间段，分别计算步骤S4中模拟的用户骑走单车的时间序列和用户归还单车的时间序列中事件的平均次数，所述事件是指用户骑走单车和用户归还单车事件，其具体次数为16次和12次；
S6、计算单车数量的预测值，具体计算公式为：
Xbike(12点30分)＝Xbike(10点)+Pbike-Pdock
S7、对步骤S6所得的预测值进行修正，其修正方法为：当所述预测值大于所述单车站点A的单车容量，则令所述预测值等于单车站点A的单车容量；当所述预测值小于零，则令所述预测值等于零；
S8、根据步骤S6和S7，计算得到单车站点A在12点30分时刻的单车数量为14。
实施例2
本实施例的运行环境为：Matlab2016b,Java(JDK 1.7),Windows 10,CPU(Core i7(7500U)，8GB RAM。
在本实施例中，对于hubway单车系统，单车站点A(i＝20)，假设当前的日期为2011年8月23号，上午9点，单车数量为8，即Xbike(9点)＝8，预测单车站点A在未来的上午10点时刻的单车数量。
具体步骤包括：
S1、根据单车站点A的历史数据，从hubway单车管理系统中提取提取过去10天(即1号～10号，不包括11号)的单车骑行时间序列数据(用户骑走单车的时间序列和用户到达并归还单车的时间序列)，将用户到达并归还单车的过程和用户到达并骑走单车的过程看作泊松过程，对其进行强度函数为μi(x)和λi(x)的非齐次泊松过程建模，并利用该模型对时间变量x求导，得到用户骑走单车和归还单车的概率密度函数；
S2、根据公式(4)-(6)并结合时间序列和采用梯度下降法分别学习得到强度函数μi(x)和λi(x)的参数c1,c2,a1,a2,b1,b2，取值分别为c1＝0.0694,a1＝9,b1＝3.253,c2＝0.07254,a2＝16,b2＝4.67；
S3、根据公式(7)和(8)计算环境参数ei(x)的值，采用的数据是当日(11号)时间范围8点至9点时间段内的单车数据，并用环境参数值对步骤S2得到的强度函数进行修正，进而实现对概率密度函数的修正；
S4、利用S3得到的概率密度函数，在9点至10点时间段，分别模拟所述单车站点A中用户骑走单车的时间序列和用户归还单车的时间序列，所述模拟的次数为5次；
S5、分别计算步骤S4中模拟的用户骑走单车的时间序列和用户归还单车的时间序列中事件的平均次数，所述事件是指用户骑走单车和用户归还单车事件，其具体次数为8次和7次；
S6、计算单车数量的预测值，具体计算公式为：
Xbike(10)＝Xbike(9)+Pbike-Pdock
S7、对步骤S6所得的预测值进行修正，其修正方法为：当所述预测值大于所述单车站点A的单车容量，则令所述预测值等于单车站点A的单车容量；当所述预测值小于零，则令所述预测值等于零；
S8、根据步骤S6和S7，计算得到单车站点A在10点时刻的单车数量为9。
实施例3
如图3和4所示，在实施例1和2两个单车系统中，对当前时刻(实施例1为2013年8月23号10点，实施例2为2011年8月23号上午9点)之后的时刻进行连续性的单车数量预测，验证了本方法的可用性以及其达到的效果，验证指标为均方根误差(root-mean-square-error,RMSE)，其计算公式为：
其中，K表示预测时间集合，为实际使用量，X(t+s)为预测使用量。
RMSE值越小，表示准确率越高。
单车预测效果如图3和图4所示，在Bay Area单车系统和Hubway单车系统上，对比方法为HA、ARMA和QMP，其中：
HA(Historical Average)，记载于Z.Yang,J.Hu,Y.Shu,P.Cheng,J.Chen,andT.Moscibroda,“Mobilitymodeling and prediction in bike-sharing systems,”inInternationalConference on Mobile Systems,Applications,and Services,2016,pp.165–178，文献名称中文翻译为：自行车共享系统的移动建模和预测；
ARMA(Auto-Regressive and Moving Average)，记载于P.Vogel andD.C.Mattfeld,“Strategic and operational planning ofbike-sharing systems bydata mining:a case study,”in InternationalConference on ComputationalLogistics,2011,pp.127–141，文献名称中文翻译为：以数据挖掘为目的的自行车共享系统的战略和业务规划：案例研究。
QMP(Queuing model predictor)，记载于N.Gast,G.Massonnet,D.Reijsbergen,and M.Tribastone,“Probabilis-tic forecasts of bike-sharing systems forjourney planning,”2015,pp.703–712，文献名称中文翻译为：关于旅行计划的自行车共享系统的概率预测。
从图3和图4的对比结果来看，在短期的预测时间内(至少在2小时内)，本发明方法的RMSE值均小于对比方法，也就是说，其预测的效果均好于对比方法，预测准确性更高。2小时对于用户或者系统管理员来说，已经可以满足要求。因为，如果用户查询一个站点是否有车，那么他应该是在很短的时间之后会去该站点骑车；对系统管理员来说，他们对单车的调度，必须尽快地完成，以满足用户的使用需求，所以需要的预测的时间也不会太久。因此，本方法具有实际的较高实际使用意义。
实施例4
如图5和图6所示，在实施例1和2的基础上，分别在考虑和不考虑环境参数对预测结果的影响的情况下，对当前时刻(实施例1为2013年8月23号10点，实施例2为2011年8月23号上午9点)之后的时刻进行连续性的单车数量预测则，同样用均方根误差(root-mean-square-error,RMSE)进行计算，结果如图5和图6所示，图中BGIP表示本发明的方法考虑环境参数的预测结果，BGIP without ei(x)表示本发明的方法不考虑环境参数的预测结果。
结果显示，在任何时间，BGIP都比BGIP without ei(x)的均方根误差小，因此，BGIP比BGIP without ei(x)的预测准确率更高，说明本发明方法考虑的环境参数对发明的技术效果贡献较大，使本发明在对站点单车进行预测时，预测准确率进一步的提高。
需要说明的是，在图3-6中，横坐标的时间初始值0，与本申请中的当前时刻t对应，在不同的实施例中，t不同，因此，初始值0对应的时刻也不同。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解，本发明不受上述实施例的限制，上述实施例和说明书中描述的只是说明本发明的原理，在不脱离本发明精神和范围的前提下，本发明还会有各种变化和改进，这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。

Claims (6)

1.一种基于二元高斯非齐次泊松过程的单车数量预测方法，其特征在于，包括以下步骤：
S1、从单车管理系统中提取单车站点预测当日之前的历史数据，所述历史数据包括过去10-15天的用户骑走单车的时间序列和用户归还单车的时间序列；
预测当日的当前时刻用t表示，未来待预测时刻用t+s表示，t取值范围为00:00:00点到23:59:59点；将用户到达并归还单车的过程和用户到达并骑走单车的过程看作泊松过程，对其进行强度函数为μi(x)和λi(x)的非齐次泊松过程建模，并利用该模型对时间变量x求导，得到用户骑走单车和归还单车的概率密度函数；其中，i表示单车站点编号；
S2、根据步骤S1中的历史数据，用二元高斯函数拟合得到μi(x)和λi(x)的表达式，并结合概率密度函数，采用迭代法学习得到μi(x)和λi(x)的参数值，得到参数赋值后的强度函数，同时将所述参数值代入步骤S1得到的概率密度函数；
S3、根据所述单车站点在当前时刻之前的1-3小时内用户骑走单车的时间序列和用户归还单车的时间序列数据，计算环境参数值，并用环境参数值对步骤S2得到的参数赋值后的强度函数进行修正，同时实现对概率密度函数的修正，得到修正后的概率密度函数；
S4、利用S3得到修正后的概率密度函数，分别模拟所述单车站点中用户骑走单车的时间序列和用户归还单车的时间序列，所述模拟的次数为5-10次；
S5、在[t,t+s]时间段，分别计算步骤S4中模拟的用户骑走单车的时间序列和用户归还单车的时间序列中事件的平均次数，所述事件是指用户骑走单车事件和用户归还单车事件，0<s≤2.5小时；
S6、计算单车数量的预测值，具体计算公式为：
Xbike(t+s)＝Xbike(t)-Pbike+Pdock
其中，Xbike(t+s)表示在t+s时刻的单车数量，Xbike(t)表示在t时刻的单车数量，Pbike表示步骤S5用户骑走单车事件的平均次数，Pdock表示步骤S5用户归还单车事件的平均次数；
S7、对步骤S6所得的预测值进行修正，其修正方法为：当所述预测值大于所述单车站点的单车容量，则令所述预测值等于所述单车站点的单车容量；当所述预测值小于零，则令所述预测值等于零；
S8、根据步骤S6和S7，计算得到所述单车站点在t+s时刻的单车数量。
2.根据权利要求1所述的一种基于二元高斯非齐次泊松过程的单车数量预测方法，其特征在于，步骤S1中，用户骑走单车的概率密度函数表达式为：
其中，表示单车站点i中，当用户在时刻骑走一辆单车，那么下一刻在时刻骑走一辆单车的概率；μi(xn)和μi(x)为强度函数，分别表示xn时刻和x时刻用户骑走单车的数量；
用户归还单车的概率密度函数表达式为：
其中，表示单车站点i中，当用户在时刻骑走一辆单车，那么下一刻在时刻骑走一辆单车的概率；λi(xn)和λi(x)为强度函数，分别表示xn时刻和x时刻用户归还单车的数量。
3.根据权利要求2所述的一种基于二元高斯非齐次泊松过程的单车数量预测方法，其特征在于，所述步骤S2中，
其中，c1,c2,a1,a2,b1,b2为强度函数表达式中的参数。
4.根据权利要求2所述的一种基于二元高斯非齐次泊松过程的单车数量预测方法，其特征在于，在步骤S3中，所述环境参数值的计算公式为：
其中的计算公式为：
表示在时间范围[t-△t,t]内单车的使用数量，△t取值范围为1～3h；
则修正后的用户骑走单车的概率密度函数公式为：
则修正后的用户归还单车的概率密度函数公式为：
5.根据权利要求4所述的一种基于二元高斯非齐次泊松过程的单车数量预测方法，其特征在于，步骤S4中，模拟用户骑走单车的时间序列的过程包括以下步骤：
S41、设初始时间t0为t，初始化用户骑走单车的时间序列变量Q；
S42、利用拒绝采样方法生成间隔时间τ，模拟用户骑走单车采用的公式为：
S43、令t0＝t0+τ，然后将t0存入Q；
S44、判断t0是否大于t+s，如果t0大于t+s，则结束模拟，否则，返回步骤S42；
模拟用户归还单车的时间序列的过程包括以下步骤：
S51、设初始时间t0为t，初始化用户骑走单车的时间序列变量Q′；
S52、利用拒绝采样方法生成间隔时间τ，模拟用户归还单车采用的公式为：
S53、令t0＝t0+τ，然后将t0存入Q′；
S54、判断t0是否大于t+s，如果t0大于t+s，则结束模拟，否则，返回步骤S42。
6.根据权利要求1所述的一种基于二元高斯非齐次泊松过程的单车数量预测方法，其特征在于，步骤S2中，所述迭代法为梯度下降法。
CN201811141977.XA 2018-09-28 2018-09-28 一种基于二元高斯非齐次泊松过程的单车数量预测方法 Active CN109325625B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811141977.XA CN109325625B (zh) 2018-09-28 2018-09-28 一种基于二元高斯非齐次泊松过程的单车数量预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811141977.XA CN109325625B (zh) 2018-09-28 2018-09-28 一种基于二元高斯非齐次泊松过程的单车数量预测方法

Publications (2)

Publication Number Publication Date
CN109325625A CN109325625A (zh) 2019-02-12
CN109325625B true CN109325625B (zh) 2019-12-17

Family

ID=65265100

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811141977.XA Active CN109325625B (zh) 2018-09-28 2018-09-28 一种基于二元高斯非齐次泊松过程的单车数量预测方法

Country Status (1)

Country Link
CN (1) CN109325625B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110191184B (zh) * 2019-05-31 2022-02-01 国信优易数据股份有限公司 一种时长信息推送方法、装置、设备和介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104835318A (zh) * 2015-03-23 2015-08-12 中南大学 一种基于贝叶斯联合模型量化径路安全水平的方法
CN105718490A (zh) * 2014-12-04 2016-06-29 阿里巴巴集团控股有限公司 一种用于更新分类模型的方法及装置
US9519912B2 (en) * 2013-09-20 2016-12-13 Xerox Corporation Hybrid system for demand prediction
CN106326998A (zh) * 2015-06-29 2017-01-11 株式会社日立制作所 路径规划方法和路径规划装置
CN106503840A (zh) * 2016-10-17 2017-03-15 中国科学院深圳先进技术研究院 停车场可用车位预测方法及系统
CN108460537A (zh) * 2018-03-23 2018-08-28 华南理工大学 一种基于大数据与用户心理预测共享自行车需求量的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9519912B2 (en) * 2013-09-20 2016-12-13 Xerox Corporation Hybrid system for demand prediction
CN105718490A (zh) * 2014-12-04 2016-06-29 阿里巴巴集团控股有限公司 一种用于更新分类模型的方法及装置
CN104835318A (zh) * 2015-03-23 2015-08-12 中南大学 一种基于贝叶斯联合模型量化径路安全水平的方法
CN106326998A (zh) * 2015-06-29 2017-01-11 株式会社日立制作所 路径规划方法和路径规划装置
CN106503840A (zh) * 2016-10-17 2017-03-15 中国科学院深圳先进技术研究院 停车场可用车位预测方法及系统
CN108460537A (zh) * 2018-03-23 2018-08-28 华南理工大学 一种基于大数据与用户心理预测共享自行车需求量的方法

Also Published As

Publication number Publication date
CN109325625A (zh) 2019-02-12

Similar Documents

Publication Publication Date Title
CN108197739B (zh) 一种城市轨道交通乘客流量预测方法
CN110458336B (zh) 一种基于深度学习的网约车供需预测方法
CN105023437B (zh) 一种公交od矩阵的构建方法及系统
CN108629503B (zh) 一种基于深度学习的出租车上车需求量的预测方法
CN109920248B (zh) 一种基于gru神经网络的公交到站时间预测方法
CN107194491A (zh) 一种基于公交客流和站间行程时间预测的动态调度方法
CN110555990A (zh) 一种基于lstm神经网络的有效停车时空资源预测方法
CN104778508A (zh) 一种基于多源数据融合的公共自行车租借预测方法
CN111160622A (zh) 基于混合神经网络模型的景区客流预测方法和装置
CN112766600B (zh) 一种城市区域人群流量预测方法及系统
CN113780665B (zh) 一种基于增强循环神经网络的私家车停留位置预测方法和系统
CN113672846A (zh) 网约车调度方法、装置、电子设备及存储介质
CN112417753A (zh) 一种基于城市公共交通资源联合调度方法
CN110889530A (zh) 基于递归神经网络的目的地预测方法及服务器
CN109858681A (zh) 一种基于ic卡的交通短时客流预测方法和相关装置
CN110443422B (zh) 基于od吸引度的城市轨道交通od客流预测方法
CN109325625B (zh) 一种基于二元高斯非齐次泊松过程的单车数量预测方法
CN109118412B (zh) 城市轨道交通网络客流在线控制系统
CN112613630B (zh) 一种融合多尺度时空统计信息的短时交通需求预测方法
CN116662815B (zh) 时间预测模型的训练方法以及相关设备
Li et al. Cycle-based signal timing with traffic flow prediction for dynamic environment
Xiangdong et al. Prediction of short-term available parking space using LSTM model
JP6926978B2 (ja) パラメータ推定装置、トリップ予測装置、方法、及びプログラム
CN113326972B (zh) 基于实时公交车速度统计数据的公交专用道短时速度预测方法
CN112667394B (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