CN107170233B - 一种基于矩阵分解的典型日交通需求od矩阵获取方法 - Google Patents
一种基于矩阵分解的典型日交通需求od矩阵获取方法 Download PDFInfo
- Publication number
- CN107170233B CN107170233B CN201710262143.3A CN201710262143A CN107170233B CN 107170233 B CN107170233 B CN 107170233B CN 201710262143 A CN201710262143 A CN 201710262143A CN 107170233 B CN107170233 B CN 107170233B
- Authority
- CN
- China
- Prior art keywords
- matrix
- distribution
- traffic demand
- time
- class
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0125—Traffic data processing
- G08G1/0129—Traffic data processing for creating historical data or processing based on historical data
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于矩阵分解的典型日交通需求OD矩阵获取方法,包括以下步骤:1)将日期为i的OD矩阵为Di展开为行向量di,并将多天的OD矩阵展开得到行向量按时间堆叠成为时间×OD对的矩阵M,日期为i的OD行向量di对应矩阵M的第i行;2)采用奇异值分解方法将矩阵M分解为三个子矩阵的乘积;3)根据规律性和爆发性指标对OD分布进行分类;4)将各类别的OD分布进行重组,得到典型日交通需求OD矩阵。与现有技术相比,本发明具有可以识别交通需求OD矩阵的结构特征和提取典型日OD矩阵、为交通需求分析和预测以及突发事件影响分析提供依据等优点。
Description
技术领域
本发明涉及交通需求分析领域,尤其是涉及一种基于矩阵分解的典型日交通需求OD矩阵获取方法。
背景技术
在交通需求分析中,通常采用出行OD矩阵来表示城市居民的出行需求的空间分布情况。出行OD矩阵的元素dij,表示第i个交通小区(Traffic Analysis Zone,TAZ)到第j个交通小区之间的出行量。
在传统交通规划和管理中,通常采用居民出行调查方法获取出行OD矩阵,该方法的成本较高,其采样代表性引起了不少学者的质疑,也难以胜任交通需求时变规律分析的要求。一方面,用于交通需求建模的出行OD矩阵通常通过5~10年一次的城市居民出行调查得到,这种城市居民出行调查的成本很高,往往只能获取少部分居民(1~5%)一天的出行信息,这种方法的得到的OD矩阵是否有代表性,是交通需求分析需要回答的问题,在实际工程中,工程师们常用“核查线”来检验OD调查的合理性。另一方面,由于日期类型(节假日)、大型活动、天气(雨雪)等,每天的出行需求在出行总量、空间分布上都可能存在差异,但具体是怎样的差异,受到各种因素的影响程度如何,是否有规律,这些问题受制于数据限制过去都无法被很好地回答。
近年新兴的数据源如公交卡、浮动车、手机信令等,可以得到各种时间粒度下连续多天、1个月、甚至1年的出行需求信息,为解决上述问题提供了可能。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种计算简单、准确、为交通需求分析、预测及突发事件影响分析提供依据的基于矩阵分解的典型日交通需求OD矩阵获取方法。
本发明的目的可以通过以下技术方案来实现:
一种基于矩阵分解的典型日交通需求OD矩阵获取方法,包括以下步骤:
1)将日期为i的OD矩阵为Di展开为行向量di,并将多天的OD矩阵展开得到行向量按时间堆叠成为时间×OD对的矩阵M,日期为i的OD行向量di对应矩阵M的第i行;
2)采用奇异值分解方法将矩阵M分解为三个子矩阵的乘积;
3)根据规律性和爆发性指标对OD分布进行分类;
4)将各类别的OD分布进行重组,得到典型日交通需求OD矩阵。
所述的步骤2)中,矩阵M分解为三个子矩阵的乘积的表达式为:
其中,r为矩阵M的秩,S为对角阵,δi为对角阵S对角线上的第i个元素,ui为矩阵U的第i列,vi为矩阵V的第i列。
所述的步骤3)具体包括以下步骤:
31)定义矩阵V的第i列vi为第i个交通需求分布模式,矩阵U的第i列ui为时间变化模式,且交通需求分布模式的值与时间变化模式的值一一对应;
32)当采用快速傅里叶变换判断时间变化模式的值具有周期性时,则该OD分布标为第一类,当时间变化模式的值中存在偏离均值超过三倍标准差的值时,则该OD分布标为第二类,其余的标记为第三类。
所述的步骤32)中,当对于持续时间小于7周的数据,则通过判断周中和周末时间变化模式的值是否存在差异,若存在,则标为第一类。
与现有技术相比,本发明具有以下优点:
本发明采用SVD方法识别OD矩阵的结构特征,通过从多天的OD矩阵提取几种典型的OD分布,然后对各种分布进行分类,提取当中特定的部分作为典型的日交通OD需求矩阵,可以为交通需求分析、预测及突发事件影响分析提供依据。
附图说明
图1为上海2011年9月的地铁站点OD矩阵的奇异值分布图。
图2为上海2011年9月的典型OD分布实例,其中,图(2a)为第一类的第一种典型OD分布实例,图(2b)为第一类的第二种典型OD分布实例,图(2c)为第二类的第一种典型OD分布实例,图(2d)为第二类的第二种典型OD分布实例,图(2e)为第三类的第一种典型OD分布实例,图(2f)为第三类的第二种典型OD分布实例。
图3为上海2011年9月的典型OD分布类型判断结果。
图4为上海2011年9月28日上海市地铁增加的典型OD矩阵期望线图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。
实施例:
本方法从连续多天的交通需求OD数据提取了几种典型的OD分布,并分析各种分布的含义,最后根据实际应用的需求将各个成分重组为典型日交通需求OD矩阵。其特征在于:使用矩阵分解方法将多天的OD矩阵分为多个成分;通过分析各个成分随时间的变化,对各个成分进行分类;并根据实际问题将各成分重组为典型日交通需求OD矩阵。
具体包括以下步骤:
1)构建“时间×OD对”矩阵M
设日期i的OD矩阵为Di,将其展开为行向量di,将多天的OD矩阵得到行向量按时间堆叠成为“时间×OD对”矩阵M,形如表1所示。
表1时间×OD对矩阵示例
2)奇异值分解,分解结果如附图1所示,图1为矩阵M的奇异值分布,图中出现明显的“陡坡”,说明少量维度可以解释多天OD矩阵中的大部分变化。
通过SVD方法,将m×n的“时间×OD对”矩阵M分解为三个矩阵的乘积:
其中矩阵下标为矩阵的维度,VT为V的转置,r为矩阵M的秩,S为对角阵,对角线上的第i个元素为δi,ui为矩阵U的第i列,vi为矩阵V的第i列。
定义vi为第i个DemandPattern(交通需求分布模式),定义ui为第i个TemporalFlow(时间变化模式)。其中DemandPattern和TemporalFlow一一对应,Demand Pattern代表了某一种特定的OD分布,Temporal Flow是一个时间序列,表征了某种Demand Pattern在各天OD中的贡献在分析时段内的变化。
3)典型OD分布分析
使用规律性、爆发性指标对典型OD分布进行分类。具体判断流程如下:
i.使用快速傅里叶变换(FFT)判断TemporalFlow是否具有周期性,即FFT结果存在峰值,若有,该OD分布标为第一类。(对于持续时间小于7周的数据,该步骤替换为使用ks检验判断周中和周末Temporal Flow的值是否有显著差异,若有,标为第一类);
ii.判断Temporal Flow中是否有偏离均值超过三倍标准差的值,若有,该OD分布标为第二类;
iii.余下的OD分布标为第三类。
不同类别的OD分布如附图2所示,附图2中第一列(图2a、2c和2e):与日期类型有关(工作日、周末),工作日较大,休息日明显降低,呈现周期性变化;第二列(图2b、2的、和2f):与某些事件有关,在某一天或两天出现爆发性需求,图(2c)和图(2d)对应的事件分别为国庆前一天以及2011年9月27日上海市地铁十号线追尾;第三列:变化无明显规律。所有OD分布的类别判断结果如附图3所示,第一类OD分布主要分布在前几位;第二类和第三类的分布无明显规律。
4)OD分布重组为典型OD矩阵
根据实际需求将OD分布重组为典型OD矩阵,例如:将第一类OD分布重组,从而滤除原始OD矩阵中爆发及随机波动的成分,可以作为预测模型的输入;将第二类OD分布重组可以得到特殊事件下OD变化的典型OD矩阵,可用于交通状况变化的建模及分析,为交通需求分析、预测及突发事件影响分析提供依据,如附图4所示,附图4为典型OD矩阵增量的期望线图,从中可以发现在公交接驳的起终点(伊犁路站和海伦路站)增加了大量的OD,OD的另一头主要是位于两端未停运的10号线站点;其次,市中心的换乘站OD大量增长,包括10号线上换乘站及当时作为站外换乘站的上海火车站。
本发明使用矩阵分解方法(SVD)从多天的OD矩阵分为多个成分,每一个成分代表一种典型的OD分布,通过分析各成分随时间的变化确定其含义,主要分为呈现规律性变化的成分、突然爆发的成分及随机发生的成分,最后根据实际应用的需求将各个成分重组为典型日交通需求OD矩阵,可以为交通需求分析、预测及突发事件影响分析及提供依据,具有良好应有价值。
Claims (1)
1.一种基于矩阵分解的典型日交通需求OD矩阵获取方法,其特征在于,包括以下步骤:
1)将日期为i的OD矩阵为Di展开为行向量di,并将多天的OD矩阵展开得到行向量按时间堆叠成为时间×OD对的矩阵M,日期为i的OD行向量di对应矩阵M的第i行;
2)采用奇异值分解方法将矩阵M分解为三个子矩阵的乘积,矩阵M分解为三个子矩阵的乘积的表达式为:
其中,r为矩阵M的秩,S为对角阵,δi为对角阵S对角线上的第i个元素,ui为矩阵U的第i列,vi为矩阵V的第i列;
3)根据规律性和爆发性指标对OD分布进行分类,具体包括以下步骤:
31)定义矩阵V的第i列vi为第i个交通需求分布模式,矩阵U的第i列ui为时间变化模式,且交通需求分布模式的值与时间变化模式的值一一对应;
32)当采用快速傅里叶变换判断时间变化模式的值具有周期性时,则该OD分布标为第一类,当时间变化模式的值中存在偏离均值超过三倍标准差的值时,则该OD分布标为第二类,其余的标记为第三类,所述的步骤32)中,当对于持续时间小于7周的数据,则通过判断周中和周末时间变化模式的值是否存在差异,若存在,则标为第一类;
4)将各类别的OD分布进行重组,得到典型日交通需求OD矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710262143.3A CN107170233B (zh) | 2017-04-20 | 2017-04-20 | 一种基于矩阵分解的典型日交通需求od矩阵获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710262143.3A CN107170233B (zh) | 2017-04-20 | 2017-04-20 | 一种基于矩阵分解的典型日交通需求od矩阵获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107170233A CN107170233A (zh) | 2017-09-15 |
CN107170233B true CN107170233B (zh) | 2020-08-18 |
Family
ID=59813296
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710262143.3A Active CN107170233B (zh) | 2017-04-20 | 2017-04-20 | 一种基于矩阵分解的典型日交通需求od矩阵获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107170233B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108615360B (zh) * | 2018-05-08 | 2022-02-11 | 东南大学 | 基于神经网络的交通需求逐日演变预测方法 |
CN110009175A (zh) * | 2018-12-25 | 2019-07-12 | 阿里巴巴集团控股有限公司 | Od需求分析算法的性能评估方法和装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279950A (zh) * | 2011-08-02 | 2011-12-14 | 中国铁道科学研究院电子计算技术研究所 | 基于数据挖掘和神经网络的轨道交通清分方法 |
CN102724078B (zh) * | 2012-06-29 | 2014-12-10 | 蒋定德 | 一种动态网络下基于压缩感知的端到端网络流量重构方法 |
CN103279534B (zh) * | 2013-05-31 | 2016-06-01 | 西安建筑科技大学 | 基于智能公交***数据的公交卡乘客通勤od分布估计方法 |
CN104200286B (zh) * | 2014-09-10 | 2017-06-06 | 东南大学 | 一种城市轨道交通时刻表优化技术应用框架 |
CN105488581B (zh) * | 2015-11-13 | 2019-09-27 | 清华大学 | 一种基于模拟退火算法的交通需求量估计方法 |
CN105489056B (zh) * | 2015-12-28 | 2018-01-26 | 中兴软创科技股份有限公司 | 一种基于od矩阵的停车需求预测方法 |
CN106301950B (zh) * | 2016-09-07 | 2019-08-09 | 中国联合网络通信集团有限公司 | 一种od流量的分析方法及分析装置 |
-
2017
- 2017-04-20 CN CN201710262143.3A patent/CN107170233B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107170233A (zh) | 2017-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111540198A (zh) | 基于有向图卷积神经网络的城市交通态势识别方法 | |
Ugille et al. | Multilevel meta-analysis of single-subject experimental designs: A simulation study | |
Xydas et al. | Electric vehicle load forecasting using data mining methods | |
CN105493109A (zh) | 使用多个数据源的空气质量推断 | |
CN109462853B (zh) | 一种基于神经网络模型的网络容量预测方法 | |
CN107170233B (zh) | 一种基于矩阵分解的典型日交通需求od矩阵获取方法 | |
CN102306297A (zh) | 一种空中交通管制员工作负荷测量方法 | |
CN103955596A (zh) | 一种基于交通事故采集技术的事故热点综合判定方法 | |
WO2010086325A1 (de) | Verfahren und sensornetz zur merkmalsauswahl für eine ereigniserkennung | |
CN105825269A (zh) | 一种基于并行自动编码机的特征学习方法及*** | |
CN104112366B (zh) | 基于隐语义模型的交通信号优化方法 | |
CN107766983A (zh) | 一种城市轨道交通车站应急救援驻车点的设置方法 | |
Schorr et al. | Safety propensity index for signalized and unsignalized intersections: Exploration and assessment | |
CN109543922B (zh) | 用于有桩共享单车站点群的分时段借还量预测方法 | |
Bouman et al. | Detecting activity patterns from smart card data | |
Hamidun et al. | Assessing pedestrian crossing risk at signalised intersection | |
CN104217263A (zh) | 一种配电开关设备关键状态量的优化方法 | |
Wang et al. | Collective behavior in the North Rhine-Westphalia motorway network | |
EP1758065B1 (de) | Verfahren zur Prognose eines Verkehrszustandes in einem Straßennetz und Verkehrsmanagementzentrale | |
Hewson | Quantile regression provides a fuller analysis of speed data | |
Jiang et al. | Heterogeneous platoon flow dispersion model based on truncated mixed simplified phase‐type distribution of travel speed | |
Kan et al. | Comparing designing across different domains: An exploratory case study | |
CN109725121B (zh) | 一种生成采样任务的方法和*** | |
Tan et al. | Statistical analysis and prediction of regional bus passenger flows | |
CN116245365A (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 |