发明名称 基于FPGA及改进Xue-Dai模型的在线交通瓶颈预测控制方法
摘要 本发明公开了一种基于FPGA及改进Xue-Dai模型的在线交通瓶颈控制方法,用于解决现有方法难以在实际的高速公路或封闭道路中对交通瓶颈进行在线预测调控的技术问题,该方法对Xue-Dai模型进行改进,把可变信息显示牌融入到Xue-Dai模型中,基于FPGA平台通过改进的Xue-Dai模型对高速公路或封闭道路整体进行预测分析,根据定义的状态变量找到道路瓶颈,进而给出匝口控制和可变信息显示牌的控制方案,并把这些控制方案按优先级带入预测模型,找到合理的控制方案,对交通瓶颈进行在线控制,使得高速公路或封闭道路中的交通瓶颈可以得到有效控制。
申请公布号 CN102938205A 申请公布日期 2013.02.20
申请号 CN201210470885.2 申请日期 2012.11.19
申请人 西安费斯达自动化工程有限公司 发明人 史忠科;刘通
分类号 G08G1/08(2006.01)I;G06Q10/04(2012.01)I;G06Q50/30(2012.01)I 主分类号 G08G1/08(2006.01)I
代理机构 代理人
主权项 1.一种基于FPGA及改进Xue-Dai模型的在线交通瓶颈预测控制方法,其特征在于包括以下步骤:步骤一、根据Xue-Dai模型:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>&rho;</mi></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>&PartialD;</mo><mrow><mo>(</mo><mi>&rho;v</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>=</mo><mi>&pi;</mi><mo>[</mo><mi>r</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>v</mi></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><mi>v</mi><mo>-</mo><mi>c</mi><mo>)</mo></mrow><mfrac><mrow><mo>&PartialD;</mo><mi>v</mi></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>&rho;</mi><mo>)</mo></mrow><mo>-</mo><mi>v</mi></mrow><mrow><mi>T</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mi>c</mi><mo>=</mo><mo>-</mo><mi>&rho;</mi><mfrac><msub><mi>&tau;</mi><mi>r</mi></msub><mrow><mi>T</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>)</mo></mrow></mrow></mfrac><mfrac><mrow><mi>d</mi><msub><mi>V</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>&rho;</mi><mo>)</mo></mrow></mrow><mi>d&rho;</mi></mfrac><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>T</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&tau;</mi><mi>r</mi></msub><mo>[</mo><mn>1</mn><mo>+</mo><mfrac><mi>E</mi><mrow><mn>1</mn><mo>+</mo><mrow><mo>(</mo><mi>&rho;</mi><mo>/</mo><msub><mi>&rho;</mi><mi>m</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>]</mo></mtd></mtr></mtable></mfenced></math>]]></maths>式中,ρ为交通流密度、v为车辆平均速度、t为时间,x为与仿真道路起始点的距离,π[r(x,t),s(x,t)]为由于匝口进入或驶出的车流量引起的密度变化率函数,r(x,t)=r<sub>0</sub>(x,t)-r<sub>q</sub>(x,t)为t时刻、x路段由匝口进入的车流量,s(x,t)=s<sub>0</sub>(x,t)+s<sub>q</sub>(x,t)为t时刻、x路段由匝口驶出的车流量,r<sub>0</sub>(x,t)、s<sub>0</sub>(x,t)为由匝口驶入驶出的正常车流量,r<sub>q</sub>(x,t)为匝口控制禁止驶入高速路造成的流量降低量,s<sub>q</sub>(x,t)为匝口控制强制驶出车辆造成的流量增量,V<sub>e</sub>(ρ)为等价速度且与自由流速度v<sub>f</sub>和交通流密度ρ相关,c为交通流的音速,T(ρ)为松弛时间,τ<sub>r</sub>、ρ<sub>m</sub>、E均为常数,其中τ<sub>r</sub>是驾驶员反映时间,ρ<sub>m</sub>是交通出现阻塞时的交通流密度,E&gt;0代表在拥挤和非拥挤的交通情况反应时间的不同,全申请书符号定义相同;把可变显示牌显示速度融入Xue-Dai模型,用可变显示牌显示速度v<sub>ind</sub>代替等价速度中的自由流速度v<sub>f</sub>,得到改进的Xue-Dai模型如下:<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>&rho;</mi></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mfrac><mrow><mo>&PartialD;</mo><mrow><mo>(</mo><mi>&rho;v</mi><mo>)</mo></mrow></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>=</mo><mi>&pi;</mi><mo>[</mo><mi>r</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>,</mo><mi>s</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>v</mi></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mrow><mo>(</mo><mi>v</mi><mo>-</mo><mi>c</mi><mo>)</mo></mrow><mfrac><mrow><mo>&PartialD;</mo><mi>v</mi></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><msub><mi>v</mi><mi>ind</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>v</mi></mrow><mrow><mi>T</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><mi>c</mi><mo>=</mo><mo>-</mo><mi>&rho;</mi><mfrac><msub><mi>&tau;</mi><mi>r</mi></msub><mrow><mi>T</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>)</mo></mrow></mrow></mfrac><mfrac><mrow><mi>d</mi><msub><mi>V</mi><mi>e</mi></msub><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><msub><mi>v</mi><mi>ind</mi></msub><mo>)</mo></mrow></mrow><mi>d&rho;</mi></mfrac><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><mi>T</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&tau;</mi><mi>r</mi></msub><mo>[</mo><mn>1</mn><mo>+</mo><mfrac><mi>E</mi><mrow><mn>1</mn><mo>+</mo><mrow><mo>(</mo><mi>&rho;</mi><mo>/</mo><msub><mi>&rho;</mi><mi>m</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>]</mo></mtd></mtr></mtable></mfenced></math>]]></maths>步骤二、定义两个新的状态变量η(x,t)、σ(x,t),当状态变量<img file="FDA00002429203400022.GIF" wi="478" he="112" />趋于无穷时,代表交通密度趋于饱和交通密度,产生交通拥堵,当状态变量<img file="FDA00002429203400023.GIF" wi="300" he="114" />趋于无穷时,代表车辆平均速度趋于零,产生交通拥堵;式中,ρ(x,t)表示t时刻x处平均交通流密度,v(x,t)表示t时刻x处车辆平均速度;步骤三、a.根据步骤一得到的改进的Xue-Dai模型,用差分格式表示微分项并略去高阶项,得到:<maths num="0003"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>&rho;</mi></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>&rho;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>+</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>-</mo><mi>&rho;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>&xi;</mi></mfrac><mo>+</mo><mi>o</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>&rho;</mi><mi>i</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup></mrow><mi>&xi;</mi></mfrac></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>&rho;</mi></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>&rho;</mi><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mi>h</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>&rho;</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>h</mi></mfrac><mo>+</mo><mi>o</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>&rho;</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup></mrow><mi>h</mi></mfrac></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>v</mi></mrow><mrow><mo>&PartialD;</mo><mi>t</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>+</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>-</mo><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>&xi;</mi></mfrac><mo>+</mo><mi>o</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>v</mi><mi>i</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>-</mo><msubsup><mi>v</mi><mi>i</mi><mi>n</mi></msubsup></mrow><mi>&xi;</mi></mfrac></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>v</mi></mrow><mrow><mo>&PartialD;</mo><mi>x</mi></mrow></mfrac><mo>=</mo><mfrac><mrow><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>+</mo><mi>h</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow></mrow><mi>h</mi></mfrac><mo>+</mo><mi>o</mi><mrow><mo>(</mo><mi>h</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>v</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>v</mi><mi>i</mi><mi>n</mi></msubsup></mrow><mi>h</mi></mfrac></mrow></math>]]></maths>式中:ξ为t的微分,h为x的微分,o(ξ)为ξ的高阶无穷小,o(h)为h的高阶无穷小,把道路分成多个路段,每个路段长度为h,采样周期为ξ,<img file="FDA00002429203400028.GIF" wi="45" he="47" />为第i个路段在[nξ,(n+1)ξ]内车辆的平均密度,<img file="FDA00002429203400029.GIF" wi="46" he="60" />为第i个路段在[nξ,(n+1)ξ]车辆的平均速度;得到改进的Xue-Dai模型的差分形式为:<maths num="0007"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>&rho;</mi><mi>i</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><mi>&xi;&pi;</mi><mrow><mo>(</mo><msubsup><mi>r</mi><mi>i</mi><mi>n</mi></msubsup><mo>,</mo><msubsup><mi>s</mi><mi>i</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mfrac><mi>&xi;</mi><mi>h</mi></mfrac><mo>[</mo><msubsup><mi>v</mi><mi>i</mi><mi>n</mi></msubsup><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>+</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mrow><mo>(</mo><msubsup><mi>v</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>v</mi><mi>i</mi><mi>n</mi></msubsup><mo>)</mo></mrow><mo>]</mo><mo>+</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup></mtd></mtr><mtr><mtd><msubsup><mi>v</mi><mi>i</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mi>v</mi><mi>i</mi><mi>n</mi></msubsup><mo>+</mo><mfrac><mrow><mi>&xi;</mi><mo>[</mo><msub><mi>V</mi><mi>e</mi></msub><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mo>,</mo><msub><mi>v</mi><mi>ind</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><msubsup><mi>v</mi><mi>i</mi><mi>n</mi></msubsup><mo>]</mo></mrow><mrow><mi>T</mi><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mrow><mi>&xi;</mi><mo>{</mo><msubsup><mi>v</mi><mi>i</mi><mi>n</mi></msubsup><mo>+</mo><mfrac><mrow><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><msub><mi>&tau;</mi><mi>r</mi></msub></mrow><mrow><mi>T</mi><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mo>)</mo></mrow></mrow></mfrac><mo>&CenterDot;</mo><mfrac><mrow><mi>d</mi><mo>[</mo><msub><mi>V</mi><mi>e</mi></msub><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mo>,</mo><msub><mi>v</mi><mi>ind</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>]</mo></mrow><mi>d&rho;</mi></mfrac><mo>}</mo><mrow><mo>(</mo><msubsup><mi>v</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mo>-</mo><msubsup><mi>v</mi><mi>i</mi><mi>n</mi></msubsup><mo>)</mo></mrow></mrow><mi>h</mi></mfrac></mtd></mtr></mtable></mfenced></math>]]></maths>式中:<img file="FDA00002429203400034.GIF" wi="45" he="60" />表示第i个路段在[nξ,(n+1)ξ]由匝口进入的车流量,<img file="FDA00002429203400032.GIF" wi="48" he="51" />表示第i个路段在[nξ,(n+1)ξ]由匝口驶出的车流量,v<sub>ind</sub>(i,n)表示第i个路段在[nξ,(n+1)ξ]内可变显示牌显示速度;b.建立等价速度模型:<maths num="0008"><![CDATA[<math><mrow><msub><mi>V</mi><mi>e</mi></msub><mo>[</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mo>,</mo><msub><mi>v</mi><mi>ind</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>]</mo><mo>=</mo><mfrac><mrow><msub><mi>v</mi><mi>ind</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mo>/</mo><msub><mi>&rho;</mi><mi>jam</mi></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>+</mo><mi>G</mi><msup><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>i</mi><mi>n</mi></msubsup><mo>/</mo><msub><mi>&rho;</mi><mi>jam</mi></msub><mo>)</mo></mrow><mn>4</mn></msup></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>式中G为常数;c.在FPGA中编写基于改进Xue-Dai模型的计算模块,包括数据接收模块、控制方案选择及数据分配模块、计算模块1-计算模块N、同步模块、数据输出模块,把道路分成N个路段,每个路段对应一个计算模块,计算模块1-计算模块N为按照前述偏微分方程组的差分解法使用浮点数运算器组合而成的交通流预测计算模块,预测控制模块的数据流向为:数据接收模块接收上位机传来的各个路段的交通流数据(交通流密度、车辆平均速度),然后传给控制方案选择及数据分配模块,控制方案选择及数据分配模块根据这些数据确定交通瓶颈,并制定调控方案,然后将使能信号、控制方案和交通流数据传给各个计算模块,各个计算模块接收到使能信号后同时对交通流密度和车辆平均速度进行预测并把结果存入寄存器,各个模块计算结束后把各自的计算结束信号传给同步模块,同步模块在所有计算模块完成计算后发送信号通知控制方案选择及数据分配模块接收交通流数据的预测结果,继续进行预测,在预测时间T<sub>c</sub>内,如果交通瓶颈解除,则采用该方案对实际交通进行调控,如果不能解除,控制方案选择及数据分配模块根据交通流数据和上次调控方案制定新的调控方案,并将交通流数据和调控方案传给各个计算模块,重新进行预测,在多次预测和调整调控方案后选择一个合适的调控方案输出对交通瓶颈进行调控,且已调控的路段在时间T<sub>c</sub>内不再进行调控,然后继续对交通进行预测,寻找新的交通瓶颈,并进行控制;所述步骤三中确定交通瓶颈并对其进行控制的方法为:求解||η(x,t)||<sub>∞</sub>=η<sub>m</sub>(x<sub>m</sub>,t<sub>m</sub>),当η<sub>m</sub>大于给定阈值η<sub>M</sub>时,说明路段x<sub>m</sub>在t<sub>m</sub>时刻将成为交通瓶颈,则在t<sub>m</sub>-T<sub>0</sub>时刻对车辆行驶方向的x<sub>m</sub>前后方入、出匝口及可变信息显示牌进行限速(瓶颈路段前方路段速度降低,后方路段速度提高)、限制进入瓶颈路段甚至强制驶出控制瓶颈路段;或求解||σ(x,t)||<sub>∞</sub>=σ<sub>m</sub>(x<sub>m</sub>,t<sub>m</sub>),当σ<sub>m</sub>大于给定阈值σ<sub>M</sub>时,说明路段x<sub>m</sub>在t<sub>m</sub>时刻将成为交通瓶颈,则在t<sub>m</sub>-T<sub>1</sub>时刻对车辆行驶方向的x<sub>m</sub>前后方出、入匝口及可变信息显示牌进行限速(瓶颈路段前方路段速度降低,后方路段速度提高)、限制进入瓶颈路段甚至强制驶出瓶颈路段;式中T<sub>0</sub>、T<sub>1</sub>为提前施加控制的时间使得||η(x,t)||<sub>∞</sub>=η<sub>m</sub>(x<sub>m</sub>,t<sub>m</sub>)≤η<sub>M</sub>、||σ(x,t)||<sub>∞</sub>=σ<sub>m</sub>(x<sub>m</sub>,t<sub>m</sub>)≤σ<sub>M</sub>,η<sub>M</sub>、σ<sub>M</sub>分别为根据道路密度最大饱和度、最小速度限制得到的正数;控制的优先级原则为:①首先通过可变信息显示牌调整路段速度,使进入瓶颈路段的车辆速度降低,驶出瓶颈路段的车辆速度提高,②仅仅通过可变信息显示牌调整路段速度不能达到控制指标时,则通过匝口限制进入瓶颈路段流量并与可变信息显示牌调整路段速度同时进行控制,③当通过匝口限制进入瓶颈路段流量及可变信息显示牌调整路段速度同时控制也不能达到控制要求时,通过匝口控制在断续时间强制部分路段车辆驶出道路、同时对匝口限制进入瓶颈路段车流量及可变信息显示牌调整路段速度以达到控制指标要求。
地址 710075 陕西省西安市高新区科技路金桥国际广场12101号