基于模糊马尔可夫过程的路网交通状态预测方法与流程

本发明涉及交通领域，具体涉及一种路网交通状态预测方法。

背景技术：

城市路网交通运行状态对于城市交通控制起着至关重要的作用，同时又为交通管理及决策提供直观信息。路网交通状态判别是一个复杂的过程，表现为多尺度、多变量、随机性和时变性。其判别的关键是建立合理的时空分布模型。

目前的交通状态判别方法主要可以分为两类：第一类是从时间角度出发，充分考虑历史及动态信息，运用数理统计、模式识别、聚类分析等方法，从路网检测器中提取表征路网交通状态的交通参数进行处理与分析；第二类是从空间角度出发，考虑构成路网物理结构的静态拓扑特性，采用评价指标如饱和度、通行能力等对路网交叉口进行评价。

在上述交通状态的判别研究方法中，大多避开其建模问题，一是因为存在合理性问题，二是很难对路网的实时运行状态进行判别。而研究路网交通状态判别的建模过程有利于达到时间与空间的结合，增强学科间的联系。

技术实现要素：

本发明提出了一种路网交通状态判别模型，考虑其动态信息和历史状态，并结合其时间信息和空间信息获得其整体状态，从而对其交通状态进行合理地分析和评价。

本发明所述的基于模糊马尔可夫过程的路网交通状态预测方法包括以下步骤：

1)在选定待预测路网交通状态的路网区域以后，以固定时间间隔采集此区域内各个交叉口的历史交通流量数据；

2)为便于量化计算，根据道路服务水平及饱和度之间的关系，建立饱和度值和量化值的映射，由此，可根据饱和度值界定交通状态；

3)根据步骤1)中所得交叉口历史交通流量数据，利用马尔可夫过程特性计算得到交叉口状态转移矩阵如式(1)所示：

其中，p_ij代表当前时刻交叉口状态为j(j＝1,2,3,4)，下一时刻状态为i(i＝1,2,3,4)的概率；

4)由式(1)并假设已知交叉口当前时刻状态向量为A_t，则可以计算得出交叉口下一时刻状态向量为A_t+1＝A_tP，从而可得交叉口下一时刻交通状态；

5)对所述路网区域内的各个交叉口进行分类，并确定各类型的数目b₁，b₂，b₃，b₄，b₅，b₆，及权重向量如式(2)所示：

B＝(B₁ B₂ B₃ B₄ B₅ B₆) (2)

其中，B_i＝b_i/(b₁+b₂+b₃+b₄+b₅+b₆)；

6)考虑步骤2)中所述四种交通状态分别在预测以及历史下一时刻中所占的数目情况，可以计算出模糊评价矩阵如式(3)所示：

其中，矩阵因子m_qi表示路网中第q(q＝1,2,3,4,5,6)类交叉口的状态为i(i＝1,2,3,4)的概率；

7)根据式(2)交叉口类型权重向量以及式(3)模糊评价矩阵，可以计算得到路网运行状态向量A＝(A₁ A₂ A₃ A₄)，从而可得所述路网区域的交通状态。

优选地，步骤1)中所述采集历史交通流量数据的固定时间间隔是指5分钟或15分钟或30分钟或60分钟，最好15分钟；

步骤2)中所述映射关系具体是指：饱和度小于等于0.3为畅通状态，饱和度大于0.3小于等于0.6为轻微拥挤状态，饱和度大于0.6小于等于0.9为拥挤状态，饱和度大于0.9为严重拥挤状态；

步骤5)中所述交叉口分类方法按照如下原则进行：主干路与主干路相交的交叉口作为第一类交叉口，主干路与次干路相交的交叉口作为第二类交叉口，主干路与支路相交的交叉口作为第三类交叉口，次干路与次干路相交的交叉口作为第四类交叉口，次干路与支路相交的交叉口作为第五类交叉口，支路与支路相交的交叉口作为第六类交叉口；

步骤6)中所述模糊评价矩阵的因子由式(4)计算得到：

其中，x_i(i＝1,2,3,4)为马尔可夫过程预测的交叉口交通状态界定值，y_i(i＝1,2,3,4)为历史数据中下一时刻交叉口交通状态界定值，u为调整参数，一般取0.75；

步骤7)中所述路网运行状态向量的各分量由式(5)计算得到：

优选地，步骤2)中所述量化值具体是指将四类交通状态：畅通，轻微拥挤，拥挤，严重拥堵分别对应为数值1,2,3,4；

步骤4)中所述交叉口当前或下一时刻状态向量的四个分量分别代表步骤2)中所述四种交通状态的概率值，其中最大的分量代表当前或下一时刻交叉口的状态；

步骤7)中所述路网运行状态向量的四个分量分别代表步骤2)中所述四种交通状态的概率值，其中最大的分量代表路网交通状态。

本发明采用马尔可夫过程对每个交叉口的交通状态进行预测，能够给出每个交叉口节点的交通状态变化情况。本发明在马尔可夫过程预测的交叉口交通状态的基础上，结合历史数据，利用模糊综合评价法分析整体路网的交通运行状态，能够给出整个路网的交通变化趋势。本发明所述路网交通状态预测方法适合于在线实时地预测交通状态并为交通管理部门提供管控决策指导。

本发明的优点是：计算过程简单、耗时短、预测精度高。

附图说明

图1是本发明中模糊马尔可夫过程预测路网交通状态流程图；

图2是本发明的实施路网区域示意图；

图3是本发明所采集的历史数据中某交叉口状态示意图；

图4是本发明对实验路网中各交叉口状态真实值预测值对比示意图；

图5是本发明中预测的路网状态描述示意图。

具体实施方式

下面结合附图进一步说明本发明。

本发明所述的基于模糊马尔可夫过程的路网交通状态预测方法包括以下步骤，见图1所示：

1)在选定待预测路网交通状态的路网区域以后，以固定时间间隔采集此区域内各个交叉口的历史交通流量数据；

2)为便于量化计算，根据道路服务水平及饱和度之间的关系，建立饱和度值和量化值的映射，由此，可根据饱和度值界定交通状态；

3)根据步骤1)中所得交叉口历史交通流量数据，利用马尔可夫过程特性计算得到交叉口状态转移矩阵如式(1)所示：

其中，p_ij代表当前时刻交叉口状态为j(j＝1,2,3,4)，下一时刻状态为i(i＝1,2,3,4)的概率；

4)由式(1)并假设交叉口当前时刻状态向量为A_t，则可以计算得出交叉口下一时刻状态向量为A_t+1＝A_tP，从而可得交叉口下一时刻交通状态；

5)对所述路网区域内的各个交叉口进行分类，并确定各类型的数目b₁，b₂，b₃，b₄，b₅，b₆，及其权重向量如式(2)所示：

B＝(B₁ B₂ B₃ B₄ B₅ B₆) (2)

其中，B_i＝b_i/(b₁+b₂+b₃+b₄+b₅+b₆)；

6)考虑步骤2)中所述四种交通状态分别在预测以及历史下一时刻中所占的数目情况，可以计算出模糊评价矩阵如式(3)所示：

其中，矩阵因子m_qi表示路网中第q(q＝1,2,3,4,5,6)类交叉口的状态为i(i＝1,2,3,4)的概率；

7)根据式(2)交叉口类型权重向量以及式(3)模糊评价矩阵，可以计算得到路网运行状态向量A＝(A₁ A₂ A₃ A₄)，从而可得所述路网区域的交通状态。

步骤1)中所述采集历史交通流量数据的固定时间间隔是指5分钟或15分钟或30分钟或60分钟，最好15分钟；

步骤2)中所述映射关系具体是指：饱和度小于等于0.3为畅通状态，饱和度大于0.3小于等于0.6为轻微拥挤状态，饱和度大于0.6小于等于0.9为拥挤状态，饱和度大于0.9为严重拥挤状态；

步骤5)中所述交叉口分类方法按照如下原则进行：主干路与主干路相交的交叉口作为第一类交叉口，主干路与次干路相交的交叉口作为第二类交叉口，主干路与支路相交的交叉口作为第三类交叉口，次干路与次干路相交的交叉口作为第四类交叉口，次干路与支路相交的交叉口作为第五类交叉口，支路与支路相交的交叉口作为第六类交叉口；

步骤6)中所述模糊评价矩阵的因子由式(4)计算得到：

其中，x_i(i＝1,2,3,4)为马尔可夫过程预测的交叉口交通状态界定值，y_i(i＝1,2,3,4)为历史数据中下一时刻交叉口交通状态界定值，u为调整参数，一般取0.75；

步骤7)中所述路网运行状态向量的各分量由式(5)计算得到：

步骤2)中所述量化值具体是指将四类交通状态：畅通，轻微拥挤，拥挤，严重拥堵分别对应为数值1,2,3,4；

步骤4)中所述交叉口当前或下一时刻状态向量的四个分量分别代表步骤2)中所述四种交通状态的概率值，其中最大的分量代表当前或下一时刻交叉口的状态；

步骤7)中所述路网运行状态向量的四个分量分别代表步骤2)中所述四种交通状态的概率值，其中最大的分量代表路网交通状态。

具体实施时，建立一个交通仿真区域路网如图2所示，该路网由主干路，次干路和支路构成，共产生32个交叉口。按所述原则对32个交叉口进行分类后，六类交叉口对应数目分别为10,8,7,2,4,1。以15分钟为时间间隔，采集并统计此路网中各交叉口近两周内上午和中午共四个时间段：7:30-7:45,7:45-8:00,11:30-11:45,11:45-12:00的交通流量。

首先，将采集到的流量数据根据道路通行能力换算成对应饱和度，再根据所述映射关系界定各交叉口状态等级，图3展示32个交叉口中某一个的状态等级。

重复利用马尔可夫过程计算出各交叉口的状态转移矩阵P_i，其中i＝<1,2,3,…,32>，代表需要计算出32个交叉口各自的状态转移矩阵。

采集路网中各个交叉口于当天7:30-7:45和11:30-11:45两个时间段的交通数据并用所述界定方法将交通数据转化为状态等级。

结合前述各个交叉口的状态转移矩阵P_i，分别由当天7:30-7:45的交叉口状态预测7:45-8:00的交叉口状态，由当天11:30-11:45的交叉口状态预测11:45-12:00的交叉口状态。两次预测的真实值和预测值对比示意图如图4所示。

由各类交叉口的数目确定其权重：

B＝(0.31 0.25 0.22 0.06 0.13 0.03)

统计预测的7:45-8:00时间段中交叉口状态为畅通、轻微拥挤、拥挤、严重拥挤的交叉口数量分别为：11,15,6,0个。

统计预测的11:45-12:00时间段中交叉口状态为畅通、轻微拥挤、拥挤、严重拥挤的交叉口数量分别为：31,1,0,0个。

结合历史数据及所预测数据，将数据输入下式并计算模糊评价矩阵的因子：

取u＝0.75得到两个时间段的模糊评价矩阵分别为：

由权重向量B，模糊评价矩阵M₁，M₂以及式：

计算可得两个预测时间段的路网运行状态向量：

A₁＝(0.18 0.53 0.27 0.02)

A₂＝(0.94 0.06 0.00 0.00)

分别求得最大分量为0.53和0.94，两个分量所对应的位置分别为轻微拥挤和畅通，也即预测的7:45-8:00的路网交通状态为轻微拥挤，预测的11:45-12:00的路网交通状态为畅通。图5以预测的时间段7:45-8:00为例，描述此时间段的路网交通状态。

本说明书实施例所述的内容仅仅是对发明构思的实现形式的列举，本发明的保护范围不应当被视为仅限于实施例所陈述的具体形式，本发明的保护范围也及于本领域技术人员根据本发明构思所能够想到的等同技术手段。