发明名称 基于隐马尔科夫模型的短时交通流状态预测方法
摘要 本发明涉及智能交通系统领域,尤其是涉及利用路段的参数值序列对短时交通流状态进行预测。一种基于隐马尔科夫模型短时交通流状态预测方法,包括以下步骤:对采集的数据进行处理和统计,通过设定预测窗口,对预测窗口起始时刻测得值以及预测窗口内参数平均值和序列对比度离散化,构成隐马尔科夫模型的隐状态和观察状态集合,然后利用Baum-Welch算法结合训练数据对模型参数进行学。最后,对于一定的预测窗口,在已知观察状态序列的基础上,利用Viterbi算法求得最优的隐状态序列,则最优隐状态序列的最后的状态即为预测状态。本发明方法可以对未来短时的交通状态进行预测,是一种有效的短时交通状态预测方法。
申请公布号 CN103413443A 申请公布日期 2013.11.27
申请号 CN201310276581.7 申请日期 2013.07.03
申请人 太原理工大学 发明人 谢刚;阎高伟;续欣莹;陈泽华;窦寿军;杨江波
分类号 G08G1/065(2006.01)I 主分类号 G08G1/065(2006.01)I
代理机构 太原科卫专利事务所(普通合伙) 14100 代理人 朱源
主权项 1.一种基于隐马尔科夫模型的短时交通流状态预测方法,其特征在于:包括如下步骤:(1)、确定隐马尔科夫模型的隐状态集合和观察状态集合,具体如下:I、以采集周期δ对通过某路段横截面的某一交通流状态参数进行采集,得到对应于该参数在已检测时段内的以采集周期δ为间隔的数据序列;Ⅱ、设定固定的时段长度作为预测窗口Φ,即短时预测时长,所述预测窗口Φ是采集周期δ的整数倍,因此预测窗口Φ内含有Φ/δ个某一交通流状态参数值组成的数据序列;设定转移窗口Δ,表示预测窗口Φ以转移窗口Δ为单位在时间轴上依次向后滑动转移,所述转移窗口Δ是采集周期δ的整数倍,范围为δ≤Δ≤Φ;据此确定在已检测时段内预测窗口Φ的数量;利用灰度联合共生矩阵C,确定每个预测窗口Φ内数据序列的对比度CON;如下:灰度联合共生矩阵C中的元素c<sub>ij</sub>表示数据点的强度值为i及其相邻的数据点的强度值为j这样的数据组合出现的频率,即<img file="FDA00003456458900011.GIF" wi="751" he="158" />则<maths num="0001"><![CDATA[<math><mrow><mi>CON</mi><mo>=</mo><msub><mi>&Sigma;</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>|</mo><mi>j</mi><mo>-</mo><mi>i</mi><mo>|</mo><mrow><mo>(</mo><mi>j</mi><mo>-</mo><mi>i</mi><mo>)</mo></mrow><msub><mi>c</mi><mi>ij</mi></msub><mo>;</mo></mrow></math>]]></maths>每个预测窗口Φ均对应有一个参数平均值<img file="FDA00003456458900013.GIF" wi="77" he="92" />每个预测窗口Φ内的起始时刻参数值θ<sub>t</sub>作为观察值,则所有观察值构成观察值序列O={O<sub>1</sub>,O<sub>2</sub>,...,O<sub>T</sub>};Ⅲ、统计所有观察值的变化范围,根据统计结果,将观察值的变化范围进行离散化为M个区间,同时得到对应于区间的等级,即将所有观察值离散化为M级,设定等级即为观察状态o<sub>i</sub>(i=1,2,...,M),得到观察状态集合O={o<sub>1</sub>,o<sub>2</sub>,...,o<sub>M</sub>};同理,统计预测窗口Φ内的参数平均值<img file="FDA00003456458900021.GIF" wi="45" he="75" />和对比度CON的变化范围,根据统计结果,将参数平均值<img file="FDA00003456458900022.GIF" wi="42" he="75" />和对比度CON分别进行离散化为m个和n个区间,同时得到对应于区间的等级,即将参数平均值<img file="FDA00003456458900023.GIF" wi="59" he="78" />离散化为m级、对比度CON离散化为n级,则利用参数平均值<img file="FDA00003456458900024.GIF" wi="42" he="72" />和对比度CON的等级的二维全因子联合描述隐状态就是m×n个,得到隐状态集合S={s<sub>1</sub>,s<sub>2</sub>,...,s<sub>N</sub>},N=m×n;(2)、在确定了隐马尔科夫模型的隐状态集合和观察状态集合之后,对隐马尔科夫模型进行训练,得到适于交通流的隐马尔科夫模型<maths num="0002"><![CDATA[<math><mrow><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mo>=</mo><mrow><mo>(</mo><mover><mi>&Pi;</mi><mo>&OverBar;</mo></mover><mo>,</mo><mover><mi>A</mi><mo>&OverBar;</mo></mover><mo>,</mo><mover><mi>B</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>具体如下:首先利用随机赋值对隐马尔科夫模型参数进行初始化,得到了隐马尔科夫初始化模型λ<sub>initial</sub>=(Π,A,B),根据λ<sub>initial</sub>=(Π,A,B)和已知的观察值序列O={O<sub>1</sub>,O<sub>2</sub>,...,O<sub>T</sub>},利用隐马尔科夫重估公式迭代得到新的隐马尔科夫模型<img file="FDA00003456458900026.GIF" wi="313" he="78" />可以证明<img file="FDA00003456458900027.GIF" wi="431" he="78" />对重估过程继续迭代直到<img file="FDA00003456458900028.GIF" wi="171" he="85" />收敛,此时的<img file="FDA00003456458900029.GIF" wi="279" he="84" />即为所求的适于交通流的隐马尔科夫模型<maths num="0003"><![CDATA[<math><mrow><mover><mi>&lambda;</mi><mo>&OverBar;</mo></mover><mo>=</mo><mrow><mo>(</mo><mover><mi>&Pi;</mi><mo>&OverBar;</mo></mover><mo>,</mo><mover><mi>A</mi><mo>&OverBar;</mo></mover><mo>,</mo><mover><mi>B</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>(3)、在给定的适用于交通流的隐马尔科夫模型和观察值序列的基础上,利用Viterbi算法求得与观察值序列对应的最优隐状态序列,则隐状态序列中的最后状态即为已检测时段后的所预测的交通流状态,并依次推移进行短时交通流状态预测。
地址 030024 山西省太原市迎泽西大街79号
您可能感兴趣的专利