发明名称 基于FPGA及改进LWR模型的在线交通瓶颈预测控制方法
摘要 本发明公开了一种基于FPGA及改进LWR模型的在线交通瓶颈控制方法,用于解决现有方法难以在实际的高速公路或封闭道路中对交通瓶颈进行在线预测调控的技术问题,该方法对LWR模型进行改进,把可变信息显示牌融入到LWR模型中,基于FPGA平台通过改进的LWR模型对高速公路或封闭道路整体进行预测分析,根据定义的状态变量找到道路瓶颈,进而给出匝口控制和可变信息显示牌的控制方案,并把这些控制方案按优先级带入预测模型,找到合理的控制方案,对交通瓶颈进行在线控制,使得高速公路或封闭道路中的交通瓶颈可以得到有效控制。
申请公布号 CN102938207B 申请公布日期 2014.07.02
申请号 CN201210470897.5 申请日期 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及改进LWR模型的在线交通瓶颈预测控制方法,其特征在于包括以下步骤:步骤一、根据LWR模型:<img file="2012104708975100001DEST_PATH_IMAGE002.GIF" wi="301" he="86" />式中, t为时间,x为当前位置与道路起始点之间距离,<img file="2012104708975100001DEST_PATH_IMAGE004.GIF" wi="23" he="21" />是交通流密度且为x、t的函数,<img file="2012104708975100001DEST_PATH_IMAGE006.GIF" wi="74" he="22" />,v为车辆平均速度且为x、t的函数,<img file="2012104708975100001DEST_PATH_IMAGE008.GIF" wi="72" he="24" />,<img file="2012104708975100001DEST_PATH_IMAGE010.GIF" wi="115" he="24" />为由于匝口进入或驶出的车流量引起的密度变化率函数,<img file="2012104708975100001DEST_PATH_IMAGE012.GIF" wi="192" he="32" />为t时刻、x路段由匝口进入的车流量,<img file="2012104708975100001DEST_PATH_IMAGE014.GIF" wi="191" he="32" />为t时刻、x路段由匝口驶出的车流量,<img file="DEST_PATH_IMAGE016.GIF" wi="58" he="30" />、<img file="DEST_PATH_IMAGE018.GIF" wi="41" he="30" />为由匝口驶入驶出的正常车流量,<img file="DEST_PATH_IMAGE020.GIF" wi="56" he="31" />为匝口控制禁止驶入高速路造成的流量降低量,<img file="DEST_PATH_IMAGE022.GIF" wi="44" he="34" />为匝口控制强制驶出车辆造成的流量增量,<img file="2012104708975100001DEST_PATH_IMAGE024.GIF" wi="92" he="26" />为等价速度且与自由流速度<img file="DEST_PATH_IMAGE026.GIF" wi="23" he="32" />和交通流密度<img file="400771DEST_PATH_IMAGE004.GIF" wi="23" he="21" />相关;把可变显示牌显示速度融入LWR模型,用可变显示牌显示速度<img file="DEST_PATH_IMAGE028.GIF" wi="35" he="33" />代替等价速度中的自由流速度<img file="263685DEST_PATH_IMAGE026.GIF" wi="23" he="32" />,得到改进的LWR模型如下:<img file="DEST_PATH_IMAGE030.GIF" wi="247" he="86" />式中,<img file="DEST_PATH_IMAGE032.GIF" wi="34" he="35" />表示可变信息显示牌显示的速度;步骤二、定义两个新的状态变量<img file="DEST_PATH_IMAGE034.GIF" wi="48" he="23" />、<img file="DEST_PATH_IMAGE036.GIF" wi="52" he="24" />,当状态变量<img file="DEST_PATH_IMAGE038.GIF" wi="169" he="46" />趋于无穷时,代表交通密度趋于饱和交通密度,产生交通拥堵,当状态变量<img file="DEST_PATH_IMAGE040.GIF" wi="104" he="45" />趋于无穷时,代表车辆平均速度趋于零,产生交通拥堵;式中,<img file="DEST_PATH_IMAGE042.GIF" wi="38" he="30" />为交通出现阻塞时的交通流密度;步骤三、a. 根据步骤一得到的改进的LWR模型,用差分格式表示微分项并略去高阶项,得到:<img file="DEST_PATH_IMAGE044.GIF" wi="299" he="156" />式中:<img file="DEST_PATH_IMAGE046.GIF" wi="14" he="22" />为对t微分的步长、即采样周期,h为对x的微分步长、即划分的每个路段长度,o(<img file="526476DEST_PATH_IMAGE046.GIF" wi="14" he="22" />)为<img file="774924DEST_PATH_IMAGE046.GIF" wi="14" he="22" />的高阶无穷小,o(h)为h的高阶无穷小,把道路分成多个路段,<img file="DEST_PATH_IMAGE048.GIF" wi="23" he="26" />为第i个路段在[n<img file="655679DEST_PATH_IMAGE046.GIF" wi="14" he="22" />,(n+1)<img file="5889DEST_PATH_IMAGE046.GIF" wi="14" he="22" />]内的交通流密度,<img file="DEST_PATH_IMAGE050.GIF" wi="25" he="33" />为第i个路段在[n<img file="751997DEST_PATH_IMAGE046.GIF" wi="14" he="22" />,(n+1)<img file="340105DEST_PATH_IMAGE046.GIF" wi="14" he="22" />]内车辆的平均速度;得到改进的LWR模型的差分形式为:<img file="DEST_PATH_IMAGE052.GIF" wi="407" he="81" />式中,<img file="DEST_PATH_IMAGE054.GIF" wi="20" he="28" />表示第i个路段在[n<img file="185570DEST_PATH_IMAGE046.GIF" wi="14" he="22" />,(n+1)<img file="272343DEST_PATH_IMAGE046.GIF" wi="14" he="22" />]内由匝口进入的车流量,<img file="DEST_PATH_IMAGE056.GIF" wi="28" he="28" />表示第i个路段在[n<img file="428999DEST_PATH_IMAGE046.GIF" wi="14" he="22" />,(n+1)<img file="668351DEST_PATH_IMAGE046.GIF" wi="14" he="22" />]内由匝口驶出的车流量,<img file="DEST_PATH_IMAGE058.GIF" wi="64" he="26" />表示第i个路段在[n<img file="153559DEST_PATH_IMAGE046.GIF" wi="14" he="22" />,(n+1)<img file="71836DEST_PATH_IMAGE046.GIF" wi="14" he="22" />]内可变显示牌显示速度;b. 建立等价速度模型:<img file="DEST_PATH_IMAGE060.GIF" wi="300" he="61" />,式中E为常数;c. 在FPGA中编写基于改进LWR模型的预测控制模块,包括数据接收模块、控制方案选择及数据分配模块、计算模块1-计算模块N、同步模块、数据输出模块,把道路分成N个路段,计算模块1-计算模块N为按照前述偏微分方程组的差分解法使用浮点数运算器组合而成的交通流预测计算模块,预测控制模块的数据流向为:数据接收模块接收上位机传来的各个路段的交通流数据,其中交通流数据包括交通流密度、车辆平均速度,然后传给控制方案选择及数据分配模块,控制方案选择及数据分配模块根据这些数据确定交通瓶颈,并制定调控方案,然后将使能信号、控制方案和交通流数据传给各个计算模块,各个计算模块接收到使能信号后同时对交通流密度和车辆平均速度进行预测并把结果存入寄存器,各个模块计算结束后把各自的计算结束信号传给同步模块,同步模块在所有计算模块完成计算后发送信号通知控制方案选择及数据分配模块接收交通流数据的预测结果,继续进行预测,在预测时间<img file="DEST_PATH_IMAGE062.GIF" wi="17" he="25" />内,如果交通瓶颈解除,则采用该方案对实际交通进行调控,如果不能解除,控制方案选择及数据分配模块根据交通流数据和上次调控方案制定新的调控方案,并将交通流数据和调控方案传给各个计算模块,重新进行预测,在多次预测和调整调控方案后选择一个合适的调控方案输出对交通瓶颈进行调控,且已调控的路段在时间<img file="753222DEST_PATH_IMAGE062.GIF" wi="17" he="25" />内不再进行调控,然后继续对交通进行预测,寻找新的交通瓶颈,并进行控制;所述步骤三中确定交通瓶颈并对其进行控制的方法为:求解<img file="DEST_PATH_IMAGE064.GIF" wi="177" he="32" />,当<img file="DEST_PATH_IMAGE066.GIF" wi="25" he="28" />大于给定阈值<img file="DEST_PATH_IMAGE068.GIF" wi="29" he="28" />时,说明路段<img file="DEST_PATH_IMAGE070.GIF" wi="25" he="28" />在<img file="DEST_PATH_IMAGE072.GIF" wi="20" he="28" />时刻将成为交通瓶颈,则在<img file="DEST_PATH_IMAGE074.GIF" wi="54" he="28" />时刻对车辆行驶方向的<img file="20649DEST_PATH_IMAGE070.GIF" wi="25" he="28" />前后方入、出匝口及可变信息显示牌进行限速使得瓶颈路段前方路段速度降低、后方路段速度提高,限制进入瓶颈路段甚至强制驶出瓶颈路段;或求解<img file="DEST_PATH_IMAGE076.GIF" wi="182" he="32" />,当<img file="DEST_PATH_IMAGE078.GIF" wi="28" he="28" />大于给定阈值<img file="DEST_PATH_IMAGE080.GIF" wi="30" he="28" />时,说明路段<img file="544602DEST_PATH_IMAGE070.GIF" wi="25" he="28" />在<img file="684596DEST_PATH_IMAGE072.GIF" wi="20" he="28" />时刻将成为交通瓶颈,则在<img file="DEST_PATH_IMAGE082.GIF" wi="53" he="28" />时刻对车辆行驶方向的<img file="451564DEST_PATH_IMAGE070.GIF" wi="25" he="28" />前后方出、入匝口及可变信息显示牌进行限速使得瓶颈路段前方路段速度降低、后方路段速度提高,限制进入瓶颈路段甚至强制驶出瓶颈路段;式中<img file="DEST_PATH_IMAGE084.GIF" wi="21" he="28" />、<img file="DEST_PATH_IMAGE086.GIF" wi="18" he="28" />为提前施加控制的时间使得<img file="DEST_PATH_IMAGE088.GIF" wi="220" he="32" />、<img file="DEST_PATH_IMAGE090.GIF" wi="228" he="32" />,<img file="524562DEST_PATH_IMAGE068.GIF" wi="29" he="28" />、<img file="679469DEST_PATH_IMAGE080.GIF" wi="30" he="28" />分别为根据道路密度最大饱和度、最小速度限制得到的正数;控制的优先级原则为:①首先通过可变信息显示牌调整路段速度,使进入瓶颈路段的车辆速度降低,驶出瓶颈路段的车辆速度提高,② 仅仅通过可变信息显示牌调整路段速度不能达到控制指标时,则通过匝口限制进入瓶颈路段流量并与可变信息显示牌调整路段速度同时进行控制,③ 当通过匝口限制进入瓶颈路段流量及可变信息显示牌调整路段速度同时控制也不能达到控制要求时,通过匝口控制在断续时间强制部分路段车辆驶出道路、同时对匝口限制进入瓶颈路段车流量及可变信息显示牌调整路段速度以达到控制指标要求。
地址 710075 陕西省西安市高新区科技路金桥国际广场12101号