Processing math: 10%

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

快速路系统反馈控制器设计及算法

谢劲松 韩印 范炳全 干宏程

谢劲松, 韩印, 范炳全, 干宏程. 快速路系统反馈控制器设计及算法[J]. 交通运输工程学报, 2008, 8(4): 99-103.
引用本文: 谢劲松, 韩印, 范炳全, 干宏程. 快速路系统反馈控制器设计及算法[J]. 交通运输工程学报, 2008, 8(4): 99-103.
XIE Jing-song, HAN Yin, FAN Bing-quan, GAN Hong-cheng. Feedback controller design and algorithm of freeway system[J]. Journal of Traffic and Transportation Engineering, 2008, 8(4): 99-103.
Citation: XIE Jing-song, HAN Yin, FAN Bing-quan, GAN Hong-cheng. Feedback controller design and algorithm of freeway system[J]. Journal of Traffic and Transportation Engineering, 2008, 8(4): 99-103.

快速路系统反馈控制器设计及算法

基金项目: 

上海市重点学科建设项目 T0502

上海市科委国际合作研究项目 062107043

上海理工大学博士启动基金项目 X693

详细信息
    作者简介:

    谢劲松(1974-), 男, 湖南新化人, 上海理工大学管理学博士研究生, 从事智能交通系统规划与管理研究

    范炳全(1942-), 男, 江苏丹阳人, 上海理工大学教授

  • 中图分类号: U491.1

Feedback controller design and algorithm of freeway system

More Information
Article Text (Baidu Translation)
  • 摘要: 为了提高快速路系统的利用效率, 设计了一个反馈控制器, 提出了反馈控制器的迭代求解算法, 控制器能实时跟踪车辆的密度, 并与所期望的密度进行比较, 把误差反馈给车辆, 通过给车辆发布恰当的速度命令来控制车辆运行, 使实际车辆密度收敛于期望密度。仿真计算结果表明: 算法取4次迭代时的车辆密度与运行速度最大相对误差为0.32%, 平均耗时为0.32 s, 因此, 设计的控制器使交通流运行趋于平稳, 减少了车辆延误, 从而提高了系统利用效率。

     

  • 拥挤已成为当今快速路系统的普遍现象, 其表现形式之一为不均一的交通流, 包括同一车道的车辆密度、速度、车头时距的不均一及不同车道间的车辆密度、速度、车头时距的不均一。为了消除拥挤, 人们尝试了不同方法, 从供给方面考虑有新建或者改、扩建已有道路的, 但无数事实证明此方法不能从根本上解决问题。从管理与控制方面考虑, 有通过合理管理和控制已有道路, 提高已有道路设施的利用效率来实现的, 如匝道控制、可变限速控制、路线诱导等策略的应用[1-3]; 有通过建立拥挤模型和收费理论来进行[4-5]; 有通过设置路段控制器来实现等[6]。但这些方法中有些是因为模型本身精度不够高, 有些是模型中的算法不能完全满足实时性的要求, 所以在实际应用当中, 还有改进的余地。本文提出的反馈控制器设计方法对算法进行了切合实际的改进与提高, 能够满足适时性的需要, 它通过实时跟踪车辆的密度, 并与期望密度比较, 从而把误差通过恰当的速度命令实时反馈给运行的车辆, 以此控制车辆运行, 达到提高交通流率及降低拥挤的目的。

    快速路宏观交通流模型主要有LWR模型、Payne模型[7]、Papageorgiou模型[1]。由于LWR模型只把密度当成状态变量, 而且使用的是稳态流密关系, 不能很好地描述车辆的瞬时行为和交通流的演进过程[8]; Payne模型和Papageorgiou模型对之做了改进, 尤其是后者在几个地方进行了实际测试[9], 并证明其有效性, 但通过一些仿真试验发现, 此模型本身也存在一些缺点。为了分析问题方便, 本文利用其改进形式[9], 在此只考虑一个车道的情形。快速路网拓扑中(图 1), 快速路被分成N段, 每段长度分别为Δi(i=1, 2, …, N), ρikqikvik分别为采样时刻ki段的密度、流率和空间平均速度, riksik分别为相应的入匝与出匝流率。模型改进形式为

    qik=αρikvik+(1α)ρ(i+1)kv(i+1)k (1)

    ρi(k+1)=ρik+ΤΔi(q(i-1)k-qik+rik-sik)(2)vi(k+1)=vik+Ττ[Vfund(ρik)-vik]+ΤΔivik[v(i-1)k-vik]-uik(3)Vfund(ρik)=vf[1-(ρiρjam)l]m(4)uik=υΤτΔiρ(i+1)k-ρikρik+κ(5)

    图  1  快速路网拓扑
    Figure  1.  Topology of freeway network

    式中: ατ为正的模型参数, 其值由实际的交通数据标定而得; T为采样周期; Vfund(ρik)代表稳态速密关系; vfρjam分别为自由流速度和阻塞密度; lm为模型参数; uik为速度的演进对车辆产生的影响。式(5)为通常的Papageorgiou模型, 式中υκ为模型参数。当uik选取不同的形式, 对车辆就会产生不同的影响, 下面的反馈控制器设计正是基于此而提出的。

    为了跟踪车辆密度, 并与期望密度比较, 以及把误差通过适当方式反馈给行驶的车辆, 需要有反馈控制器来产生恰当的速度命令给车辆, 车辆据此命令调整自身行为, 从而达到期望的密度。为了得到满足要求的反馈控制器, 定义每一路段的跟踪误差为

    ξik=ρik-ρdik(6)

    因为反馈控制器可以通过设计uik来实现, 不同的uik产生不同的速度和密度, 从而实现交通密度ρik指数收敛于期望密度ρdik, 简记为ρikexpρdik,k

    为了简化分析, 可考虑没有出入口的快速路系统。首先给出引理如下。

    引理  考虑离散系统

    x(k+1)=cx(k)+z(k),x(0)=x0(7)

    其中c为常数, 且c < 1, 那么当z(k)exp0,k, 有x(k)exp0,kz(k)为k时刻不同状态变量, x0为初值。

    定理  如果存在满足式(8)的反馈控制器Uk, 则快速路的交通密度指数收敛于期望密度。

    ΡΝkUk=Ek(8)ΡΝk=[αb1(k+1)2α-1c1(k+1)a2(k+1)b2(k+1)c2(k+1)a(Ν-1)(k+1)b(Ν-1)(k+1)c(Ν-1)(k+1)aΝ(k+1)αbΝ(k+1)2α-1]Uk=[u1ku2kuΝk]ΤEk=[e1ke2keΝk]Τaik=ΤΔiαρ(i-1)k,bik=ΤΔi(1-2α)ρikcik=-ΤΔi(1-α)ρ(i+1)kfik=vik+Ττ[Vfund(ρik)-vik]+ΤΔivik[v(i-1)k-vik]

    eik=ai(k+1)f(i-1)k+bi(k+1)fik+ci(k+1)f(i+1)k+di(k+1)-cζζikdik=ρik-ρdi(k+1)-cξξikζik=ΤΔi[q(i-1)k-qik]+ρik-ρdi(k+1)-cξξikθik=-ai(k+1)u(i-1)k-bi(k+1)uik-ci(k+1)u(i+1)k+eik

    证明: 由上述定义得系统

    {ξi(k+1)=cξξik+ζikζi(k+1)=cζζik+θik(9)

    其中cξcζ为常数, 0 < cξ < 0, 0 < cζ < 1。

    由引理知, 如果恰当设计uik使θik为0成立, 则可得ζikexp0,k, 从而ξikexp0kρikexpρdik,k。考虑到边界的交通情况[6]: 进入路段的平均速度等于第一段的平均速度, 退出路段的平均速度、交通密度分别等于最后一段的平均速度、交通密度, 并利用式(1), 得如下的边界条件

    {ρ0k=q0k/v1k-(1-α)ρ1kαv0k=v1kρ(Ν+1)k=ρΝkv(Ν+1)k=vΝk(10)

    由此求得

    {a1(k+1)=Τq0(k+1)Δ1v1(k+1)+1-α2α-1b1(k+1)cΝ(k+1)=1-α2α-1bΝ(k+1)(11)

    于是求解θik=0, 并利用式(10)、(11)即得uik, 从而定理得证。

    由于矩阵PNk可逆, 所以满足要求的uik总存在。求得uik后知k+1时刻给车辆的速度命令为

    vcom=fik-uik

    此时交通密度指数收敛于期望密度。

    由上可知, 为了获得每一段的控制律uik, 须求解

    Uk=ΡΝk-1Ek

    如果N很大, 已有的求解P Νk-1的算法需要花费大量的计算时间。为此下面给出一个改进的迭代算法, 在保证精度的前提下可以大大节省计算时间, 从而提高计算效率。

    分解PNk

    ΡΝk=Lk+Rk

    利用泰劳展式有

    P Νk-1 =L k-1 -L k-1 RkL k-1 +L k-1 RkL k-1 RkL k-1 -…

    Uk=L k-1 Ek-L k-1 RkL-1Ek+L k-1 RkL k-1 RkL k-1 Ek-…

    Lk=[αb1(k+1)2α-1a2(k+1)b2(k+1)a3(k+1)b3(k+1)a(Ν-1)(k+1)b(Ν-1)(k+1)aΝ(k+1)αbΝ(k+1)2α-1]Rk=[0c1(k+1)0c2(k+1)0c3(k+1)0c(Ν-1)(k+1)0]

    Z k0 =L k-1 Ek

    Z kn =L k-1 RkZ kn-1

    Uk=j=0(-1)jΖkj

    Uk可用前n阶泰劳展式近似替代, 即

    Ukn=j=0n-1(-1)jΖkjUk

    于是在采样时刻k, 在路段i处的反馈控制器近似表示为u¯iknu¯iknUkn中的元素, 从而得如下迭代算法。

    步骤1: 求解Z k0

    Z k0 =[Z 1k0 Z 2k0Z Νk0 ]T

    LkZ k0 =Ek

    {α2α-1b1(k+1)Ζ1k0=e1ka2(k+1)Ζ1k0+b2(k+1)Ζ2k0=e2kaΝ(k+1)Ζ(Ν-1)k0+α2α-1bΝ(k+1)ΖΝk0=eΝk

    步骤2:求解Z kn

    Z kn =[Z 1kn Z 2knZ Νkn ]T

    LkZ kn =RkZ kn-1

    {α2α-1b1(k+1)Ζ1kn=c1(k+1)Ζ2kn-1a2(k+1)Ζ1kn+b2(k+1)Ζ2kn=c2(k+1)Ζ3kn-1aΝ(k+1)Ζ(Ν-1)kn+α2α-1bΝ(k+1)ΖΝkn=0

    步骤3:求解Ukn

    从后面的仿真分析知, Ukn虽然只是Uk的近似, 但是当选取适当的模型参数后, 4次迭代就能到达很高的精度。

    本文考虑一个单向单车道的快速路系统, 其由8个路段组成, 每段的长度均为1 km, 初始时刻进入路段1的流率假设为1 500 veh·h-1。速度和密度初始值及相应的模型参数见表 12[10]

    表  1  初始速度和密度
    Table  1.  Initial speeds and densities
    路段 1 2 3 4 5 6 7 8
    速度/(km·h-1) 81 81 29 29 81 81 81 81
    密度/(veh·km-1) 18 18 52 52 18 18 18 18
    下载: 导出CSV 
    | 显示表格
    表  2  模型参数
    Table  2.  Model parameters
    α l m τ κ υ T vf ρjam cξ cζ ρdi
    0.95 1.86 4.05 20 40 60 10 110 110 0.85 0.85 25
    下载: 导出CSV 
    | 显示表格

    显然从初始密度可知, 如果不采取控制措施, 那么路段的拥挤将很快由中间路段蔓延至上游路段, 并最终导致阻塞形成。下面是采纳控制措施的情形, 不同之处在于对上面所提算法分别进行了1~4次迭代, 以此来比较其相应的性能指标, 从而为采取合理可行的措施提供依据。在给定入口处的流率为定值1 500 veh·h-1, 期望密度为25 veh·km-1的情况下, 仿真1 h, 结果见表 3图 2~9

    表  3  迭代后性能指标
    Table  3.  Performance indexes after iterations
    下载: 导出CSV 
    | 显示表格
    图  2  1次迭代后的密度演化曲线
    Figure  2.  Evolution graph of density after one iteration
    图  3  1次迭代后的速度演化曲线
    Figure  3.  Evolution graph of speed after one iteration
    图  4  2次迭代后的密度演化曲线
    Figure  4.  Evolution graph of density after two iterations
    图  5  2次迭代后的速度演化曲线
    Figure  5.  Evolution graph of speed after two iterations
    图  6  3次迭代后的密度演化曲线
    Figure  6.  Evolution graph of density after three iterations
    图  7  3次迭代后的速度演化曲线
    Figure  7.  Evolution graph of speed after three iterations
    图  8  4次迭代后的密度演化曲线
    Figure  8.  Evolution graph of density after four iterations
    图  9  4次迭代后的速度演化曲线
    Figure  9.  Evolution graph of speed after four iterations

    表 3知, 迭代1次虽耗时少, 平均只需0.12 s, 但是误差较大, 其相对误差达到8.48%, 而且在速度演化曲线图上, 出现了不太合理的数值(大于110 km·h-1), 密度曲线离期望密度也有一定距离。虽然迭代次数增加会降低相对误差, 但由此所需的时间会指数增加, 这不利于实际应用。所以取4次迭代既可以满足精确要求(最大相对误差只有0.32%), 又可以满足实际耗时需要(平均耗时只有0.32 s)。同时从密度曲线可知, 在控制器作用下交通拥挤现象很快就被消除, 并且快速路各段的密度都趋于期望密度。从速度曲线可知, 所有车辆以一致的速度行驶。这样快速路上的车辆很快就达到均一的速度和密度, 交通流运行趋于平稳, 从而提高了系统流率, 减少了车辆延误。

    本文基于改进的宏观离散交通流模型, 对快速路系统进行了建模分析, 在实时跟踪车辆密度与期望密度误差的基础上提出了反馈控制器设计方法, 它以给车辆发布恰当速度命令为手段来控制车辆运行, 使系统中车辆达到一致的速度与密度, 然后提出了求解此控制器的有效的迭代算法, 并对一个单车道单方向的快速路系统进行了仿真分析, 结果表明算法在取4次迭代时既能达到很好的精度, 又能满足实时性的需要。仿真分析亦显示快速路系统运行平稳高效, 由此知本文方法有效。

  • 图  1  快速路网拓扑

    Figure  1.  Topology of freeway network

    图  2  1次迭代后的密度演化曲线

    Figure  2.  Evolution graph of density after one iteration

    图  3  1次迭代后的速度演化曲线

    Figure  3.  Evolution graph of speed after one iteration

    图  4  2次迭代后的密度演化曲线

    Figure  4.  Evolution graph of density after two iterations

    图  5  2次迭代后的速度演化曲线

    Figure  5.  Evolution graph of speed after two iterations

    图  6  3次迭代后的密度演化曲线

    Figure  6.  Evolution graph of density after three iterations

    图  7  3次迭代后的速度演化曲线

    Figure  7.  Evolution graph of speed after three iterations

    图  8  4次迭代后的密度演化曲线

    Figure  8.  Evolution graph of density after four iterations

    图  9  4次迭代后的速度演化曲线

    Figure  9.  Evolution graph of speed after four iterations

    表  1  初始速度和密度

    Table  1.   Initial speeds and densities

    路段 1 2 3 4 5 6 7 8
    速度/(km·h-1) 81 81 29 29 81 81 81 81
    密度/(veh·km-1) 18 18 52 52 18 18 18 18
    下载: 导出CSV

    表  2  模型参数

    Table  2.   Model parameters

    α l m τ κ υ T vf ρjam cξ cζ ρdi
    0.95 1.86 4.05 20 40 60 10 110 110 0.85 0.85 25
    下载: 导出CSV

    表  3  迭代后性能指标

    Table  3.   Performance indexes after iterations

    下载: 导出CSV
  • [1] PAPAGEORGIOU M, KOTSI ALOS A. Freeway ramp metering: an overview[J]. IEEE Transactions on Intelligent Transportation Systems, 2002, 3(4): 271-281. doi: 10.1109/TITS.2002.806803
    [2] KOTSI ALOS A, PAPAGEORGIOU M, MANGEAS M, et al. Coordinated and integrated control of motorway net-works via non-linear optimal control[J]. Transportation Research Part C, 2002, 10: 65-84. doi: 10.1016/S0968-090X(01)00005-5
    [3] ZHANG H M, RECKER W W. On optimal freeway rampcontrol policies for congested traffic corridors[J]. Transportation Research Part B, 1999, 33: 417-436. doi: 10.1016/S0191-2615(98)00045-9
    [4] 李志纯, 谷强, 史峰. 弹性需求下拥挤道路收费的模型与算法研究[J]. 交通运输工程学报, 2001, 1(3): 81-85. http://transport.chd.edu.cn/article/id/200103020

    LI Zhi-chun, GU Qiang, SHI Feng. Toll model and algorithm of roadjammed with traffic based on elastic demand[J]. Journalof Traffic and Transportation Engineering, 2001, 1(3): 81-85. (in Chinese) http://transport.chd.edu.cn/article/id/200103020
    [5] 史峰, 李志纯. 拥挤道路使用收费的理论构架[J]. 交通运输工程学报, 2002, 2(2): 78-82. http://transport.chd.edu.cn/article/id/200202019

    SHI Feng, LI Zhi-chun. Theoretic framework and development on research of congested road-use pricing[J]. Journal of Traffic and Transportation Engineering, 2002, 2(2): 78-82. (in Chinese) http://transport.chd.edu.cn/article/id/200202019
    [6] CHIEN C C, ZHANG Y, IOANNOU P A. Traffic density control for automated highway systems[J]. Automatica, 1997, 33(7): 1273-1285.
    [7] PAYNE H J. Models of freeway traffic and control[J]. Mathematical Models of Public Systems, 1971, 1: 51-61.
    [8] KOHAN R R. Robust state estimation and control of highway traffic systems[D]. Toronto: University of Toronto, 2001.
    [9] PAPAGEORGIOU M, BLOSSEVILLE J M, HADJ-SALEM H. Modelling and real time control of traffic flow on the southern part of boulevard periherique in Paris[J]. Transportation Research Part A, 1990, 24: 345-359.
    [10] WANG Y, PAPAGEORGIOU M, MESSMER A. A realtime freeway network traffic surveillance tool[J]. IEEE Transactions on Control Systems Technology, 2006, 14(1): 18-32.
  • 加载中
图(9) / 表(3)
计量
  • 文章访问数:  330
  • HTML全文浏览量:  144
  • PDF下载量:  226
  • 被引次数: 0
出版历程
  • 收稿日期:  2008-01-11
  • 刊出日期:  2008-08-25

目录

/

返回文章
返回