[发明专利]一种基于耦合隐马尔可夫模型的城市道路异常交通流检测方法

权利要求书

1.一种基于耦合隐马尔可夫模型的城市道路异常交通流检测方法，所述交通流是指交通管理部门的中心计算机中的道路交通信息；其特征在于：依据车辆GPS轨迹数据进行异常交通流建模和在线实时检测两部分内容。

2.根据权利要求1所述的基于耦合隐马尔可夫模型的城市道路异常交通流检测方法，其特征在于异常交通流建模的步骤有：

步骤一：时间段划分；

从提取交通流量的开始时间至结束时间记为总时间T，当使用长度为τ的单位时间对所述的T进行等长划分称为时间段，即T＝{t₁,t₂,…,t_j-1,t_j,t_j+1,…,t_l}，j,l表示时间段的标识号，t₁表示第一个时间段，t₂表示第二个时间段，t_j-1表示时间段t_j的前一个时间段，t_j表示当前一个时间段，也称为任意一个时间段，t_j+1表示时间段t_j的后一个时间段，t_l表示最后一个时间段；

步骤二：统计各个时间段下各个道路的交通流量；

对于任意一时间段t_j，任意一道路Road_M拥有一个状态变量和一个观察变量所述观察变量是指道路Road_M在时间段t_j里统计得到的交通流量，而交通流量是行驶于道路Road_M上的车辆数目；

对于任意一道路Road_M在总时间T的所有观察变量记为AOTRoadM={Ot1RoadM,Ot2RoadM,...,OtjRoadM,...,OtlRoadM};]]>

对于所有道路在总时间T的所有观察变量记为BO={AOTRoad1,AOTRoad2,...,AOTRoadM};]]>

步骤三：进行道路交通流状态和交通流量之间概率关系的获取；

在步骤三中，分别记先后两次迭代计算得到的状态-流量似然值为lh_bm_old和lh_bm_new；

在步骤三中，设置状态-流量阈值为ε-BM，且ε-BM＞0；所述的ε-BM是用来衡量lh_bm_old和lh_bm_new之间的差别是否为足够小；

步骤301：初始化状态-平均流量

通过任意道路Road_M在总时间T内的车辆GPS轨迹数据中统计得出的最大流量最小流量然后采用μStateNRoadM=(StateN-12)×OmaxRoadM-OminRoadMN-1+OminRoadM]]>得到初始化的状态-平均流量执行步骤302；

步骤302：初始化状态先验概率

任意道路Road_M处于状态State_N下的初始化状态先验概率为执行步骤303；

步骤303：依据状态-平均流量与观察变量计算，得到状态-观测概率！表示阶乘运算，e为自然底数；执行步骤304；

步骤304：依据状态先验概率和状态-观测概率计算，得到道路交通流状态后验概率分布的充分统计量执行步骤305；

步骤305：依据道路交通流状态后验概率分布的充分统计量和道路的流量计算，更新状态-平均流量为μStateNRoasM=Σtj=1TγStateNRoadMOtjRoadMT;]]>执行步骤306；

步骤306：依据道路交通流状态后验概率分布的充分统计量计算，更新状态先验概率为swStateNRoadM=Σtj=1TγStateN,tjRoadMT;]]>执行步骤307；

步骤307：依据状态先验概率和状态-观测概率计算，得到本次迭代的状态-流量概率分布似然值lh_bmnew=Σtj=1TΣRoadM=1Mln(ΣStateN=1NswStateNRoadMgStateN,tjRoadM);]]>

若不是第一次执行步骤307，则执行步骤308；

若是第一次执行步骤307，则执行步骤309；

步骤308：若lh_bm_new-lh_bm_old＜ε-BM，则结束步骤三；若lh_bm_new-lh_bm_old≥ε-BM，执行步骤309；

步骤309：依据本次迭代的状态-流量概率似然值lh_bm_new更新上次迭代的状态-流量概率似然值lh_bm_old；返回执行步骤303；

步骤四：进行相邻道路交通流状态转移概率关系的获取；

在步骤四中，分别记前后两次迭代计算得到的道路模型似然值为lh_cl_old和lh_cl_new；

在步骤四中，设置道路模型阈值为ε-CL，满足ε-CL＞0；所述道路模型阈值ε-CL用来衡量lh_cl_old和lh_cl_new之间的差别是否为足够小；

步骤401：初始化状态转移概率；

任意道路RoadM的任意一条相邻道路处于状态下并使得道路Road_M处于状态state_N的初始化状态转移概率为执行步骤402；

步骤402：初始化线性耦合系数

任意道路Road_M的所有相邻道路中，任意一相邻道路对道路Road_M的影响占有的初始化权重为执行步骤403；

步骤403：依据状态-平均流量与观察变量计算，得到状态-观测概率其中！表示阶乘运算，e为自然底数；执行步骤404；

步骤404：依据任意道路Road_M在t₁时段的状态-观测概率及状态先验概率计算，得到初始状态分布πStateNRoadM=swStateNRoadMgStateN,t1RoadMΣStateN=1NswStateNRoadMgStateN,t1RoadM;]]>执行步骤405；

步骤405：依据任意道路Road_M在t₁时段的状态-观测概率及初始状态分布计算，得到任意道路Road_M在t₁时段的联合概率变量δStateN,t1RoadM=πStateNRoadM×gStateN,t1RoadM,]]>执行步骤406；

步骤406：设置当前的时间段标号j＝1，执行步骤407；

步骤407：依据任意道路Road_M在t_j时段的联合概率变量计算，得到任意道路Road_M在t_j时段的流量概率变量执行步骤408；

步骤408：依据任意道路Road_M在t_j时段的联合概率变量及流量概率变量计算，得到任意道路Road_M在t_j时段的前向概率变量执行步骤409；

步骤409：若当前时间段标号j＜l，则执行步骤410；若当前时间标号j＝l，则执行步骤412；

步骤410：对于任意道路Road_M和其任意一条相邻道路依据线性耦合系数在t_j时段的前向概率变量状态转移概率以及在t_j+1时段的状态-观测概率计算，得到道路Road_M在t_j+1时段的联合概率变量δStateN,tj+1RoadM=gStateN,tj+1RoadMΣm=1||link||(ωlinkmRoadM→RoadMΣStateNlinkmRoadM=1NαStateN,tjRoadMAStateNlinkmRoadM→StateNlinkmRoadM→RoadM),]]>执行步骤411；

步骤411：将当前时间标号j增加1，执行步骤407；

步骤412：对于任意道路Road_M及任意状态State_N，初始化任意道路Road_M在时间段T下的后向概率变量为执行步骤413；

步骤413：设置当前时间段标号为j＝l；执行步骤414；

步骤414：若当前时间段标号j＞1，则将当前时间段标号减一(j-1),并执行步骤415；否则，若当前时间段标号j＝1，执行步骤416；

步骤415：对于任意道路Road_M和其任意一条相邻道路依据线性耦合系数状态转移概率以及在t_j+1时段的状态-观测概率流量概率变量计算，得到道路Road_M在当前时间段的后验概率变量为

βStateN,tjRoadM=1ctj+1RoadMΣm=1||link||(ωlinkmRoadM→RoadMΣStateNlinkmRoadM=1N(AStateNlinkmRoadM→StateNlinkmRoadM→RoadMβStateN,tj+1RoadMgStateN,tj+1RoadM))]]>；执行步骤412；

步骤416：对于任意道路RoadM和其任意一条相邻道路依据状态转移概率线性耦合系数前向概率变量后向概率变量相邻道路的流量概率变量道路的流量概率变量状态-观测概率计算，得到状态转移后验概率ξStateNlinkmRoadM→StateN,tj+1linkmRoadM→RoadM=αStateNlinkmRoadM,tjlinkmRoadMωlinkmRoadM→RoadMAStateNlinkmRoadM→StateNlinkmRoadM→RoadMgStateN,tj+1RoadMβStateN,tj+1RoadMctjlinkmRoadMctj+1RoadM]]>执行步骤417；

步骤417：更新状态转移概率

对于任意道路Road_M和其任意一条相邻道路依据状态转移后验概率计算，得到更新的状态转移概率AStateNlinkmRoadM→StateNlinkmRoadM→RoadM=ξStateNlinkmRoadM→StateN,tjlinkmRoadM→RoadMΣStateN=1NξStateNlinkmRoadM→StateN,tjlinkmRoadM→RoadM,]]>执行步骤418；

步骤418：更新线性耦合系数

对于任意道路Road_M和其任意一条相邻道路依据更新的状态转移概率计算，得到更新的线性耦合系数ωlinkmRoadM→RoadM=ΣStateNlinkmRoadM=1NΣStateN=1NAStateNlinkmRoadM→StateNlinkmRoadM→RoadMΣm=1||link||ΣStateNlinkmRoadM=1NΣStateN=1NAStateNlinkmRoadM→StateNlinkmRoadM→RoadM,]]>执行步骤419；

步骤419：依据任意时段t_j的流量概率变量计算，得到本次迭代的道路模型似然值lh_clnew=Σtj=1TΣRoadM=1MlnctjRoadM;]]>

若不是第一次执行步骤419，则执行步骤420；

若是第一次执行步骤419，则执行步骤421；

步骤420：若lh_cl_new-lh_cl_old＜ε-CL，则结束步骤四；否则，若lh_cl_new-lh_cl_old≥ε-CL，则执行步骤421；

步骤421：依据本次迭代计算得出的道路模型似然值lh_cl_new更新上次迭代计算得出的道路模型似然值lh_cl_old；返回执行步骤404。