[发明专利]基于FPGA及改进离散宏观P模型的在线交通瓶颈预测控制方法有效
申请号: | 201210470946.5 | 申请日: | 2012-11-19 |
公开(公告)号: | CN102945611A | 公开(公告)日: | 2013-02-27 |
发明(设计)人: | 史忠科;刘通 | 申请(专利权)人: | 西安费斯达自动化工程有限公司 |
主分类号: | G08G1/08 | 分类号: | G08G1/08;G06Q10/04;G06Q50/30 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 710075 陕西省西安*** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明公开了一种基于FPGA及改进离散宏观P模型的在线交通瓶颈控制方法,用于解决现有方法难以在实际的高速公路或封闭道路中对交通瓶颈进行在线预测调控的技术问题,该方法对离散宏观P模型进行改进,把可变信息显示牌融入到离散宏观P模型中,基于FPGA平台通过改进的离散宏观P模型对高速公路或封闭道路整体进行预测分析,根据定义的状态变量找到道路瓶颈,进而给出匝口控制和可变信息显示牌的控制方案,并把这些控制方案按优先级带入预测模型,找到合理的控制方案,对交通瓶颈进行在线控制,使得高速公路或封闭道路中的交通瓶颈可以得到有效控制。 | ||
搜索关键词: | 基于 fpga 改进 离散 宏观 模型 在线 交通 瓶颈 预测 控制 方法 | ||
【主权项】:
1.一种基于FPGA及改进离散宏观P模型的在线交通瓶颈预测控制方法,其特征在于包括以下步骤:步骤一、根据离散宏观交通流P模型:x(n+1)=x(n)+f[x(n)]+B(n)u(n)式中:x(n+1)=[k1(n+1) v1(n+1) k2(n+1) v2(n+1) … kN(n+1) vN(n+1)]Tu ( n ) = u 1 ( n ) u 2 ( n ) . . . u N T ( n ) T , ]]> u1(n)=ak0(n)v0(n)+r1(n)-s1(n),ui(n)=ri(n)-si(n),(i=2,…,N-1),u N T ( n ) = [ r N ( n ) - s N ( n ) - ( 1 - a ) k out ( n ) v out ( n ) ] , ]]>f [ x ( n ) ] = - T L 1 [ a k 1 ( n ) v 1 ( n ) + ( 1 - a ) k 2 ( n ) v 2 ( n ) ] T τ [ v e ( k 1 ( n ) ) - v 1 ( n ) ] + T L 1 k 0 ( n ) k 1 ( n ) + λ v 1 ( n ) [ v 0 ( n ) v 1 ( n ) - v 1 ( n ) ] - μ ( n ) T τ L 1 ω 1 ( n ) T L 2 [ a k 1 ( n ) v 1 ( n ) + ( 1 - 2 a ) k 2 ( n ) v 2 ( n ) - ( 1 - a ) k 3 ( n ) v 3 ( n ) ] T τ [ v e ( k 2 ( n ) ) - v 2 ( n ) ] + T L 2 k 1 ( n ) k 2 ( n ) + λ v 2 ( n ) [ v 1 ( n ) v 2 ( n ) - v 2 ( n ) ] - μ ( u ) T τ L 2 ω 2 ( n ) · · · T L N - 1 [ a k N - 2 ( n ) v N - 2 ( n ) + ( 1 - 2 a ) k N - 1 ( n ) v N - 1 ( n ) - ( 1 - a ) k N ( n ) v N ( n ) ] T τ [ v e ( k N - 1 ( n ) ) - v N - 1 ( n ) ] + T L N - 1 k N - 2 ( n ) k N - 1 ( n ) + λ v N - 1 ( n ) [ v N - 2 ( n ) v N - 1 ( n ) - v N - 1 ( n ) ] - μ ( u ) T τ L N - 1 ω N - 1 ( n ) T L N [ a k N - 1 ( n ) v N - 1 ( n ) + ( 1 - 2 a ) k N ( n ) v N ( n ) ] T τ [ v e ( k N ( n ) ) - v N ( n ) ] + T L N k N - 1 ( n ) k N ( n ) + λ v N ( n ) [ v N - 1 ( n ) v N ( n ) - v N ( n ) ] - μ ( u ) T τ L N ω N ( n ) ]]>B ( n ) = B 1 T B 2 T · · · B N - 1 T B N T T , ]]>B i T = T L i 0 0 ]]> (i=1,…,N)ω i ( n ) = k i + 1 ( n ) - k i ( n ) k i ( n ) + λ , ]]> i=1,2,…N;n=0,1,2,…式中,T为采样周期,Li表示第i个路段,ki(n)表示第i个路段在[nT,(n+1)T]内的平均交通流密度,vi(n)表示第i个路段在[nT,(n+1)T]内车辆的平均速度,ri(n)=ri0(n)-riq(n)为第i个路段在[nT,(n+1)T]内由匝口进入的车流量,si(n)=si0(n)+siq(n)是第i个路段在[nT,(n+1)T]内由匝口驶出的车流量,ri0(n)、si0(n)为由匝口驶入驶出的正常车流量,riq(n)为匝口控制禁止驶入高速路造成的流量降低量,siq(n)为匝口控制强制驶出车辆造成的流量增量,ve(k)是等价速度,τ,ξ,λ,a,Γ是常数,全申请书符号定义相同;把可变显示牌显示速度融入离散宏观P模型,用可变显示牌显示速度vind代替自由流速度vf,得到改进的离散宏观P模型如下:x(n+1)=x(n)+f[x(n)]+B(n)u(n)式中:x(n+1)=[k1(n+1) v1(n+1) k2(n+1) v2(n+1) … kN(n+1) vN(n+1)]Tu ( n ) = u 1 ( n ) u 2 ( n ) . . . u N T ( n ) T , ]]> u1(n)=ak0(n)v0(n)+r1(n)-s1(n),ui(n)=ri(n)-si(n),(i=2,…,N-1),u N T ( n ) = [ r N ( n ) - s N ( n ) - ( 1 - a ) k out ( n ) v out ( n ) ] , ]]>f [ x ( n ) ] = - T L 1 [ a k 1 ( n ) v 1 ( n ) + ( 1 - a ) k 2 ( n ) v 2 ( n ) ] T τ [ v e ( k 1 ( n ) , v ind ( 1 , n ) ) - v 1 ( n ) ] + T L 1 k 0 ( n ) k 1 ( n ) + λ v 1 ( n ) [ v 0 ( n ) v 1 ( n ) - v 1 ( n ) ] - μ ( n ) T τ L 1 ω 1 ( n ) T L 2 [ a k 1 ( n ) v 1 ( n ) + ( 1 - 2 a ) k 2 ( n ) v 2 ( n ) - ( 1 - a ) k 3 ( n ) v 3 ( n ) ] T τ [ v e ( k 2 ( n ) , v ind ( 2 , n ) ) - v 2 ( n ) ] + T L 2 k 1 ( n ) k 2 ( n ) + λ v 2 ( n ) [ v 1 ( n ) v 2 ( n ) - v 2 ( n ) ] - μ ( u ) T τ L 2 ω 2 ( n ) · · · T L N - 1 [ a k N - 2 ( n ) v N - 2 ( n ) + ( 1 - 2 a ) k N - 1 ( n ) v N - 1 ( n ) - ( 1 - a ) k N ( n ) v N ( n ) ] T τ [ v e ( k N - 1 ( n ) , v ind ( N - 1 , n ) ) - v N - 1 ( n ) ] + T L N - 1 k N - 2 ( n ) k N - 1 ( n ) + λ v N - 1 ( n ) [ v N - 2 ( n ) v N - 1 ( n ) - v N - 1 ( n ) ] - μ ( u ) T τ L N - 1 ω N - 1 ( n ) T L N [ a k N - 1 ( n ) v N - 1 ( n ) + ( 1 - 2 a ) k N ( n ) v N ( n ) ] T τ [ v e ( k N ( n ) , v ind ( N , n ) ) - v N ( n ) ] + T L N k N - 1 ( n ) k N ( n ) + λ v N ( n ) [ v N - 1 ( n ) v N ( n ) - v N ( n ) ] - μ ( u ) T τ L N ω N ( n ) ]]>B ( n ) = B 1 T B 2 T · · · B N - 1 T B N T T , ]]>B i T = T L i 0 0 ]]> (i=1,…,N)ω i ( n ) = k i + 1 ( n ) - k i ( n ) k i ( n ) + λ , ]]> i=1,2,…N;n=0,1,2,…式中,vind(i,n)表示第i个路段在[nT,(n+1)T]内可变信息显示牌显示速度;步骤二、定义两个新的状态变量η(i,n)、σ(i,n),当状态变量
趋于无穷时,代表交通密度趋于饱和交通密度,产生交通拥堵,当状态变量
趋于无穷时,代表车辆平均速度趋于零,产生交通拥堵;式中,kjam是交通出现阻塞时的交通流密度;步骤三、a.建立等价速度模型:v e ( k i ( n ) , v ind ( i , n ) ) = v ind ( i , n ) { 1 - [ k i ( n ) k jam ] l } m , ]]> 式中,l,m为常数;b.在FPGA中编写基于改进离散宏观P模型的预测控制模块,包括数据接收模块、控制方案选择及数据分配模块、计算模块1-计算模块N、同步模块、数据输出模块,把道路分成N个路段,每个路段对应一个计算模块,计算模块1-计算模块N为按照前述偏微分方程组的差分解法使用浮点数运算器组合而成的交通流预测计算模块,预测控制模块的数据流向为:数据接收模块接收上位机传来的各个路段的交通流数据(交通流密度、车辆平均速度),然后传给控制方案选择及数据分配模块,控制方案选择及数据分配模块根据这些数据确定交通瓶颈,并制定调控方案,然后将使能信号、控制方案和交通流数据传给各个计算模块,各个计算模块接收到使能信号后同时对交通流密度和车辆平均速度进行预测并把结果存入寄存器,各个模块计算结束后把各自的计算结束信号传给同步模块,同步模块在所有计算模块完成计算后发送信号通知控制方案选择及数据分配模块接收交通流数据的预测结果,继续进行预测,在预测时间Tc内,如果交通瓶颈解除,则采用该方案对实际交通进行调控,如果不能解除,控制方案选择及数据分配模块根据交通流数据和上次调控方案制定新的调控方案,并将交通流数据和调控方案传给各个计算模块,重新进行预测,在多次预测和调整调控方案后选择一个合适的调控方案输出对交通瓶颈进行调控,且已调控的路段在时间Tc内不再进行调控,然后继续对交通进行预测,寻找新的交通瓶颈,并进行控制;所述步骤三中确定交通瓶颈并对其进行控制的方法为:求解||η(i,n)||∞=ηm(im,nm),当ηm(im,nm)大于给定阈值ηM时,说明路段im在nmT时将成为交通瓶颈,则在nmT-T0时刻对车辆行驶方向的路段im前后方入、出匝口及可变信息显示牌进行限速(瓶颈路段前方路段速度降低,后方路段速度提高)、限制进入瓶颈路段甚至强制驶出瓶颈路段;或求解||σ(i,n)||∞=σm(im,nm),当σm(im,nm)大于给定阈值σM时,说明路段im在nmT时刻将成为交通瓶颈,则在nmT-T1时刻对车辆行驶方向的路段im前后方出、入匝口及可变信息显示牌进行限速(瓶颈路段前方路段速度降低,后方路段速度提高)、限制进入瓶颈路段甚至强制驶出瓶颈路段;式中T0、T1为提前施加控制的时间使得||η(i,n)||∞=ηm(im,nm)≤ηM、||σ(i,n)||∞=σm(im,nm)≤σM,ηM、σM分别为根据道路密度最大饱和度、最小速度限制得到的正数;控制的优先级原则为:①首先通过可变信息显示牌调整路段速度,使进入瓶颈路段的车辆速度降低,驶出瓶颈路段的车辆速度提高,②仅仅通过可变信息显示牌调整路段速度不能达到控制指标时,则通过匝口限制进入瓶颈路段流量并与可变信息显示牌调整路段速度同时进行控制,③当通过匝口限制进入瓶颈路段流量及可变信息显示牌调整路段速度同时控制也不能达到控制要求时,通过匝口控制在断续时间强制部分路段车辆驶出道路、同时对匝口限制进入瓶颈路段车流量及可变信息显示牌调整路段速度以达到控制指标要求。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西安费斯达自动化工程有限公司,未经西安费斯达自动化工程有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201210470946.5/,转载请声明来源钻瓜专利网。