发明名称 考虑多车的OVDM交通流跟驰模型稳定性建模方法
摘要 本发明公开了一种考虑多车的OVDM交通流跟驰模型稳定性建模方法,用于解决现有的OVDM微观交通流跟驰模型对ITS环境适应性差的技术问题。技术方案是在OVDM微观交通流跟驰模型的基础上,通过建立新的交通模型,得到交通拥堵问题和系统稳定性的关系,采用局部稳定性角度从微观上判定交通是否会出现拥堵或其他异常现象。为交通控制、决策提供基本依据,可以在ITS环境中直接应用处理交通拥堵问题。
申请公布号 CN102800193A 申请公布日期 2012.11.28
申请号 CN201210308469.2 申请日期 2012.08.27
申请人 西北工业大学 发明人 史忠科;周杰
分类号 G08G1/00(2006.01)I;G06F19/00(2011.01)I 主分类号 G08G1/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种考虑多车的OVDM交通流跟驰模型稳定性建模方法,其特征在于包括以下步骤:步骤一、建立微观交通流模型<maths num="0001"><![CDATA[<math><mrow><mfrac><mrow><msup><mi>d</mi><mn>2</mn></msup><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msup><mi>dt</mi><mn>2</mn></msup></mfrac><mo>=</mo><mi>a</mi><mo>[</mo><mi>V</mi><mrow><mo>(</mo><mi>&Delta;</mi><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><msub><mi>dx</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>]</mo><mo>+</mo><mi>&lambda;</mi><mfrac><mrow><mi>d&Delta;</mi><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><mo>+</mo><mi>&gamma;</mi><mo>[</mo><mi>V</mi><mrow><mo>(</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mi>V</mi><mrow><mo>(</mo><mi>&Delta;</mi><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>式中,x<sub>j</sub>(t)是时刻t第j辆车所在位置,<img file="FDA00002062276800013.GIF" wi="128" he="115" />是时刻t第j辆车的速度,Δx<sub>j</sub>(t)是连续的两辆车之间的车头间距,V(·)是优化速度函数,a是驾驶员的敏感系数,λ是相对速度差的反应参数,γ是优化速度差的反应参数;步骤二、若给定优化速度函数<maths num="0003"><![CDATA[<math><mrow><mi>V</mi><mrow><mo>(</mo><mi>&Delta;</mi><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mi>v</mi><mi>max</mi></msub><mn>2</mn></mfrac><mo>[</mo><mi>tanh</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&alpha;</mi><mi>l</mi></msub><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>l</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>h</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>tanh</mi><mrow><mo>(</mo><msub><mi>h</mi><mi>c</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>式中,v<sub>max</sub>为最大速度,h<sub>c</sub>为安全间距;α<sub>l</sub>满足以下条件(1)α<sub>l</sub>随l递增而递减,即α<sub>l</sub>&lt;α<sub>l-1</sub>;(2)<img file="FDA00002062276800015.GIF" wi="196" he="119" />当n=1时,α<sub>l</sub>=1;令<img file="FDA00002062276800016.GIF" wi="189" he="106" />l≠n-1;<img file="FDA00002062276800017.GIF" wi="194" he="106" />l=n-1;为方便行文,记<maths num="0004"><![CDATA[<math><mrow><mi>V</mi><mrow><mo>(</mo><mi>&Delta;</mi><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>n</mi><mo>-</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>&equiv;</mo><mi>V</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&alpha;</mi><mi>l</mi></msub><msub><mi>&Delta;x</mi><mrow><mi>j</mi><mo>+</mo><mi>l</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>步骤三、将上式改写为<maths num="0005"><![CDATA[<math><mrow><mfrac><mrow><msup><mi>d</mi><mn>2</mn></msup><mi>&Delta;</mi><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><msup><mi>dt</mi><mn>2</mn></msup></mfrac><mo>=</mo><mi>a</mi><mo>[</mo><mi>V</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&alpha;</mi><mi>l</mi></msub><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>l</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mi>d&Delta;</mi><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>]</mo><mo>+</mo><mi>&lambda;</mi><mo>[</mo><mfrac><mrow><mi>d</mi><msub><mi>&Delta;x</mi><mrow><mi>j</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>-</mo><mfrac><mrow><mi>d&Delta;</mi><msub><mi>x</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>dt</mi></mfrac><mo>]</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mo>+</mo><mi>&gamma;</mi><mo>[</mo><mi>V</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&alpha;</mi><mi>l</mi></msub><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>l</mi><mo>+</mo><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mn>2</mn><mi>V</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&alpha;</mi><mi>l</mi></msub><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>l</mi><mo>+</mo><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>+</mo><mi>V</mi><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>&alpha;</mi><mi>l</mi></msub><mi>&Delta;</mi><msub><mi>x</mi><mrow><mi>j</mi><mo>+</mo><mi>l</mi></mrow></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo></mrow></math>]]></maths>式中,Δx<sub>j</sub>(t)为状态变量;根据所建立的微观交通流模型,画出(Δx,a)的二维相平面图,根据线性稳定性条件以及扭结-反扭结波的结果,分别得到共存曲线及中性稳定曲线,将二维相平面分成稳定区域、亚稳定区域及不稳定区域;求出临界点(h<sub>c</sub>,a<sub>c</sub>);在不同的参数、状态、输入条件下做相平面及临界点;当状态变量Δx趋于h<sub>c</sub>时,说明交通流趋于稳定;当状态变量Δx在h<sub>c</sub>附近来回波动时,说明交通流处于不稳定状态。
地址 710072 陕西省西安市友谊西路127号