logo资料库

基于元胞自动机的交通流模型研究.pdf

第1页 / 共10页
第2页 / 共10页
第3页 / 共10页
第4页 / 共10页
第5页 / 共10页
第6页 / 共10页
第7页 / 共10页
第8页 / 共10页
资料共10页,剩余部分请下载后查看
基于元胞自动机的交通流模型研究 田园 摘要:本文研究了当前一些主流元胞自动机交通流模型,并对各种不同模型从单、双车 道和更新规则进行了分类与评价。通过计算机数值模拟各种模型,探索模拟结果中出现 的非线性现象,并结合实测交通现象,验证模型的适用性,对各种元胞自动机交通流模 型的特性及各自的优缺点进行了分析和讨论。 关键字:交通流模型;交通仿真;元胞自动机;亚稳态;相分离 Abstract: In this paper, a various kinds of popular traffic cellular automaton models are studied; Classification and evaluation of both the single lane and double lane models are accomplished. Simulation is done to compare the nonlinear phenomenon of the models; Fitness of the models to actual environment is examined. And characteristics of different models are discussed and compared. Key words: Traffic flow model, Traffic simulation, Cellular automaton, Metastable states, Phase separation 1 前言 1.1 研究背景和意义 在我国,城市人口众多、交通工具多样化加剧了交通问题的尖锐性,科学管理和控 制混合交通流是必须解决的迫切问题。而利用交通流模型研究、模拟交通状况,是解决 这一问题的有力手段[1]。近年来元胞自动机交通流模型的研究受到了广泛关注。元胞自 动机是一种时间、空间和变量均离散的数学模型,和以往的模型相比,它不需要具体的 公式,只需要给出对应的演化规则,具有算法简单且灵活可调、计算效率高等特点,更 适合于计算机模拟。通过元胞之间的微观交互作用,可以模拟出宏观的交通流演变规律, 是研究交通流特性的有效工具。 1.2 元胞自动机交通流模型概述 元胞自动机是时间和空间都离散的动力学系统[1]。散布在规则网格中的每一元胞取 有限的离散状态,依据确定的局部规则作同步更新,大量元胞通过简单的相互作用而构 成动态系统的演化。将元胞自动机赋予车辆交通的含义: 在建立元胞自动机交通流模型中,时间、空间以及速度都被整数离散化。道路被划 分为若干个离散的格子(即元胞),一个元胞对应一辆或几辆汽车,或是几个元胞对应 一辆汽车,每个元胞状态或空或是容纳车辆的速度(每辆车的速度可取 0 到 Vmax),每 辆车都同时按照所建立的规则运动。这些规则由车辆运动应遵守的运动规则和交通规则 组成,并且包含驾驶行为、外界干扰等随机变化规则。由于交通元素从本质上来说是离 散的,用元胞自动机理论来研究交通,就避免了离散—连续—离散的近似过程,因此有 其独特的优越性[2]。其主要优点是: (1)模型简单,特别易于在计算机上实现。 (2)能够再现各种复杂的交通现象,反映交通流特性[3]。 2 数值模拟程序流程设计及参数设置
本文用 matlab 对各种元胞自动机交通流模型在周期性边界条件下的演化进行数值 模拟,并对得到的结果进行比较和分析。假设车辆行进拥有理想的交通条件和道路条件: 每部车都在一条道路上连续不断的行驶,不受其他方向车辆的干扰;每条车道的宽度均 大于等于 3.65m,道路的侧向余宽大于等于 1.75m;路面状况良好、视野开阔等等。 在模拟中取单车道由 1000 个格点(元胞)组成,每个元胞对应的实际长度为 7.5m, 则对应的实际道路长度大约为 7.5km。假设道路上由单一类型的车辆组成且最大车速 Vmax=5cell/s,对应的实际车速约为 135km/h。采用周期性边界条件,总演化时步为 1000 步,为了消除初始位形的随机性对结果的影响,取后 500 步的数值结果作时间平均。仿 真程序的流程图设计见图 1: 图 1 仿真程序流程图 3 元胞自动机交通流模型仿真结果对比分析 (1) 单车道交通流模型仿真及分析 ① NaSch 模型的仿真及分析 (a)自由流堵塞现象 (b)行走波现象 图 2 NaSch 模型时空图 开始输入参数随机产生车辆输出车辆位置与速度根据更新规则车辆运动运行步数等于T?结束YN
图 2 显示的是 NaSch 模型的时空图,两图参数中仅初始车辆数不同,图(a)中初始 车辆总数为 140,图(b)为 250。图中均匀部分表示车辆行驶在自由状态,黑色区域表示 车辆处于拥堵状态。由图可知,NaSch 模型可以模拟出自发产生的堵塞现象以及拥挤交 通情况下的时走时停波等。行走波在时空图中是一系列的带状分布[2],如图(b)所示,在 带状之内车辆密度较高,而在带状之间车辆密度较低。反映在实际问题中,对应在路 面常常出现的时走时停的交通现象,即一辆车离开拥挤区之后不久又由于前面的堵塞 而停下来[3]。 图 3 显示了 NaSch 模型在不同初始状态下的基本图。显然,在道路上没有车辆即密 度 k=0 时,q=0;在密度达到最大值即道路上发生致密堵塞时,交通流量也降为 0;交通 流量在中间密度范围内存在一个最大值。密度流量关系图看起来像希腊字母λ 的镜像, 这个反λ 的两个分支分别用来定义自由流和拥挤流[4]。由图 3 可看出:NaSch 模型并不 能模拟出亚稳态和回滞现象,且不同初始状态下对应的车流量相差不大。NaSch 模型模 拟得出的基本图在快要达到通行能力的地方没有呈现出应有的流量间断的现象,交通量 较实际交通量要小,与实际情况有出入。 (a)密度—流量图 (b) 密度—速度图 图 3 NaSch 模型基本图 图 4 显示了在不同最大车速和慢化率影响下的密度—流量关系图,从图(a)中可知, 当最大速度 Vmax 取 1 时,NaSch 模型就退化为 Wolfram 的 184 号模型。图(b)显示,在 NaSch 模型中,随着随机慢化概率的增大,临界最大流量值及最大流量值对应的密度值 都大幅度的减小,可见,随机慢化概率对交通流的影响特别大。 (a) 不同最大车速下密度—流量图 (b) 不同慢化概率下密度—流量图 图 4 NaSch 模型密度—流量图 00.10.20.30.40.50.60.70.80.9100.050.10.150.20.250.30.350.40.45密度—流量图密度流量 曲线一曲线二曲线一初始状态为均匀分布曲线二初始状态为致密堵塞00.10.20.30.40.50.60.70.80.9100.511.522.533.544.55密度—速度图密度速度 曲线1曲线2曲线1初始状态为均匀分布曲线2初始状态为致密静止分布00.10.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.7密度—流量图密度流量 Vmax=1Vmax=2Vmax=3Vmax=4Vmax=500.10.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.7密度—流量图密度流量 P=0.16P=0.25P=0.3P=0.4P=0.5
图 5 显示了在不同最大速度和慢化率影响下的密度—速度关系图,从图中可以看出, 随着车辆的最大速度增大,道路上车辆的平均速度也增大,临界最大流量值也增大,临 界最大流量所对应的密度值减小,这与实际是相符合的。随着随机慢化概率的增大,车 辆的平均速度将大幅度减小,且随着随机慢化概率的增大,车辆的平均速度减小的坡度 将增大。 (a) 不同最大车速下的密度—速度图 (b) 不同慢化概率下的密度—速度图 图 5 NaSch 模型密度—速度图 从模拟得到的基本图和时空图可分析出,NS 模型很好地反映了车流运动的宏观特 性,主要表现在: a.呈现起动一停车波(start-to-stop waves)现象,即在时空图上,可以清晰地看到车流 逐渐由低密度下畅行到高密度下拥挤的运动过程,形象地表现出车流运动象波一样在车 队中传播的景象。 b.模拟得到了连续的流量一密度曲线,与实测结果极为相似,显示出交通流的两种 状态—拥挤状态和非拥挤状态。 综上所述,NaSch 模型利用非常简单的规则再现实际交通中观测到的某些现象,比 如可以模拟出自发产生的堵塞现象等。同时也存在不足之处,比如 NaSch 模型的临界密 度远低于实际交通中的数值、不能重现交通流的重要特性,如相分离和滞后现象等。 ② 基于 NaSch 的衍生模型仿真及分析 A.巡航驾驶极限模型仿真及分析 在巡航驾驶极限模型中,随机慢化只对车速 V
(a)密度—流量图 (b) 密度—速度图 图 6 巡航驾驶极限模型基本图 在巡航驾驶极限模型中,当车流密度增加至一个临界值 k*时,车流突然从自由流状 态转变为拥挤状态。对所有 k
k*会随着慢化概率 P 的减小而增大;而当 P 0 时, 。如图(b)所示, 在给定慢化概率的情况下,k*会随着最大速度 Vmax 的减小而增大。 B.引入慢启动规则模型仿真及分析 慢启动模型只对速度为零的车辆采用慢启动规则,其他车辆均以相同概率减速。慢 启动规则的引入不仅可以模拟出临界点附近亚稳态和回滞现象,在高密度区还可以模拟 出相分离的现象[1]。在 NaSch 模型基础上引入慢启动规则的经典模型有 TT、BJH、VDR 等,由于篇幅限制,只介绍 VDR 模型的仿真及分析: 图8是VDR模型在不同慢化概率下的时空图。如图(b)所示,在慢启动概率较小的情 况下,系统趋近于NaSch模型的演化形态。图(a)显示了VDR模型在慢启动概率较大时, 系统表现出了典型的自由流和堵塞流的分化状态[4]。图中箭头所指的分界处是阻塞消散 区,在这个区域由于阻塞而停止的车辆在慢启动规则下逐渐加速,驶入下游畅行区域, 从而出现了阻塞相和畅行相并存的相位分离现象[2]。同时观察到由于慢启动规则的引入 使消散后的车流密度明显小于原来的初始密度,这是由于模型中当车辆进入阻塞区域而 减速停下后,在阻塞消散时由于慢启动规则的作用而不能及时加速,从而出现了现实交 通流中的“加速一减速不对称”的现象[3],诱使模拟的交通流中出现了滞后现象。图(a) 清晰地展示出:均匀亚稳态具有一定的生存寿命[4],在经过一定的时间演化后,该均匀 态会逐渐衰退并导致相分离现象的发生。这种高密度下相分离状态的微观结构与TT模型 和BJH模型中的情况定性相同,然而却同NaSch模型有着显著的差异。 (a)慢启动概率 P0=0.75 (b) 慢启动概率 P0=0.01 图 8 VDR 模型的典型时空图 以上现象在NaSch模型中是不会出现的,且VDR模型的最大流量和实际交通流流量 相符。可见通过对随机规则的改进,VDR模型在模拟实际交通流现象上明显优于原先的 NaSch模型。仿真结果表明,T—T模型、BJH模型、VDR模型的基本图都具备亚稳态和 回滞特征,在高密度区还可以模拟出相分离的现象,用这些模型所做的数值模拟得到的 交通容量均优于NaSch模型。 C.速度效应模型仿真及分析 图 9 为速度效应模型在确定性条件和非确定性条件下得到的基本图,图(a)能看到亚 稳态和回滞现象,不过在 p 时回滞现象消失。
(a)P=0 (b)P≠0 图 9 VE 模型基本图 VE 模型得到的基本图较之 NaSch 交通量更大,显然与 NaSch 模型相比要更加接近 实测数据。这是因为 NaSch 模型的车辆速度更新规则把前车作为静止的粒子来处理,由 此造成模拟速度小于实际车辆速度,对伴有随机慢化的交通流,所得基本图其流量远小 于实测数据。 D.FI 模型仿真及分析 图 10 显示了 FI 模型的基本图,可看出 FI 模型的密度—流量关系表现出了自由流的 排列特点。在达到临界密度之前,流量保持线性增长,平均速度保持在最大车速附近; 在达到临界密度之后,流量、平均速度均开始逐渐下降,直至为 0。 (a)密度—流量图 (b) 密度—速度图 图 10 FI 模型基本图 将 NaSch 和 FI 模型的基本图进行比较可知,在相同的参数设置条件下,FI 模型仿 真得到的交通流量值大于 NaSch 模型仿真得到的值,相较于 NaSch 模型而言,FI 模型 能得到更符合实际交通流情况的结果。 图 11 显示了 FI 模型在初始分布车流密度为 K=0.3 情况下的时空图。可看出,即 使在较大的全局密度下也没有出现与 NaSch 模型时空图类似的行走波现象。因此可知 FI 模型不能描述真实交通中出现的时停时走现象,但是它可以得到精确解。 00.10.20.30.40.50.60.70.80.9100.20.40.60.811.21.4密度—流量图密度流量 曲线1曲线2曲线1初始状态为均匀分布曲线2初始状态为致密静止分布00.10.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.70.8密度—流量图密度流量 曲线1曲线2曲线1初始状态为均匀分布曲线2初始状态为致密静止分布00.10.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.70.8密度—流量图密度流量 曲线1曲线2曲线1初始状态为均匀分布曲线2初始状态为致密静止分布00.10.20.30.40.50.60.70.80.9100.511.522.533.544.55密度—速度图密度速度 曲线1曲线2曲线1初始状态为均匀分布曲线2初始状态为致密静止分布
图 11 FI 模型的时空图 (2) 双车道交通流模型仿真及分析 图 12 模拟的是双车道道路长度均为 L=1000 个格点(元胞),每个元胞对应的实际 长度为 7.5m,对应的实际道路长度大约为 7.5km,车辆的最大速度 Vmax=5cell/s,对应 的实际车速约为 135km/h,周期性边界条件下,假设交通流由单一均匀车辆组成,模型 演化 1000 步所获得的基本图。 (a) 密度—流量图 (b) 密度—速度图 图 12 均匀系统下双车道模型基本图 由图 12 可见,与 NaSch 模型的基本图比较可得,在相同车流密度的情况下,STNS 模型的车流量要比 NaSch 模型的要大,这也验证了实际交通流中的情况。 图 13 显示的是在不同随机慢化概率下,STNS 的换道频率分布图。如图(a)显示出 了一个有趣的现象,当随机慢化概率 p 取值较小时,换道概率存在着两个峰值,并且在 最大车流量所对应的密度附近存在着一个局部最小值。通过考察换道过程中所必须的空 元胞数来解释这种现象: 在自由流区换道时通常需要在目标道上至少有 2Vmax+1 个空的元胞,如图(b)所示。 这样,在均匀系统的慢化概率比较小时,可以发现在密度 Ks= 附近换道频 率达到一个局部最大值;而对于慢化概率比较大时(P=0.5),没有局部最大值的出现。 对于 P 比较小的情况,随着密度的进一步增加,换道频率出现一个最小值。在 P 的 00.10.20.30.40.50.60.70.80.9100.10.20.30.40.50.60.7密度—流量图密度流量00.10.20.30.40.50.60.70.80.9100.511.522.533.544.55密度—速度图密度速度
分享到:
收藏