交通流理论元胞自动机模型概要课件

1、元胞自动机模型主讲人： 李新刚办公地点：8710（51684936）Email: lixingang第1页，共77页。主要内容2 元胞自动机的定义和构成3 184号规则4 NS模型简介 5 BML模型简介 元胞自动机交通流模型6 双车道模型简介 1 绪论第2页，共77页。1 绪论Stephen Wolfram. A New Kind of Science. Wolfram Media, 2002.第3页，共77页。1 绪论“三个世纪以前，人们发现建立在数学方程基础上的规律能够用于对自然界的描述，伴随着这种新观念，科学发生了转变。 在此书中我的目的是将要用简单的电脑程序来表达更为一般类型的规律，

2、并在此种规律基础上建立一种新的科学，从而启 动另一场科学变革。”著名的物理学家、数学家和计算机科学家S. Wolfram以这样的惊世之言开始了他的宏篇巨著一种新科学。Wolfram认为传统科学未能建立起解释宇宙复杂性的理论，靠数学方程做不到这一点。所以他要发动一场新的科学革命，革命的内容就是要用简单的电脑程序取代数学方程。Wolfram所钟情的这种简单电脑程序的核心基础就是我们将要介绍的元胞自动机。第4页，共77页。 元胞自动机（Cellular Automata，简称CA）实质上是定义在一个由具有离散、有限状态的元胞组成的元胞空间上，并按照一定的局部规则，在离散的时间维度上演化的动力学系统。

3、1 绪论第5页，共77页。在元胞自动机中，空间被一定形式的规则网格分割为许多单元。这些规则网格中的每一个单元都称为元胞(cell)，并且它只能在有限的离散状态集中取值。所有的 元胞遵循同样的作用规则，依据确定的局部规则进行更新。大量的元胞通过简单的相互作用而构成动态系统的演化。1 绪论第6页，共77页。1 绪论元胞自动机发展历程20世纪50年代，John von Neumann 最早提出； (von Neumann,J.1963，collected works, edited by A.H.Taub)1970年，John Conway 提出生命游戏 (Conway, J. (1970). In

4、 M. Gardner, (Ed.), Scientific American, 223(4), pp. 120-123.)1983年，Stephen Wolfram 初等元胞自动机 (Stephen Wolfram. Reviews of Modern Physics,1983,Vol.55. Stephen Wolfram. Nature,1984,Vol.311)1986年至今，理论及应用 第7页，共77页。1 绪论元胞自动机自产生以来，被广泛地应用到社会、经济、军事和科学研究的各个领域。到目前为止，其应用领域涉及生物学、 生态学、物理学、化学、交通科学、计算机科学、信息科学、地理、环境

5、、社会学、军事学以及复杂性科学等。 第8页，共77页。1 绪论元胞自动机应用生物学领域：因为元胞自动机的设计思想本身就来源于生物学自繁殖的现象，所以它在生物学上的应用更为自然而广泛。 例如元胞自动机用于肿瘤细胞的增长机理和过程模拟、人类大脑的机理探索、爱滋病病毒HIV的感染过程、自组织、自繁殖等生命现象的研究以及最新流行的克隆 (clone)技术的研究等。另外， 元胞自动机还可以用来模拟植物的生长过程以及贝壳上的色素沉积图案。第9页，共77页。1 绪论元胞自动机应用生态学领域：元胞自动机被用于兔子-草，鲨鱼-小鱼等生态系统动态变化过程的模拟，展示出令人满意的动态效果;元胞自动机还成功地应用于蚂

6、蚁的行走路径，大雁、鱼类洄游等动物的群体行为的模拟；另外，基于元胞自动机模型的生物群落的扩散模拟也是当前的一个应用热点。第10页，共77页。1 绪论元胞自动机应用物理学领域：在元胞自动机基础之上发展出来的格子气自动机(LGA)和格子-波尔兹曼方法(LBM)在计算流体领域获得了 巨大的成功。不仅能够解决传统流体力学计算方法所能解决的绝大多数问题，并且在多孔介质、 多相流、微小尺度方面具有其独特的优越性。格子-波尔兹曼方法还被成功地应用于磁场、电场、热扩散和热传导的模拟。另外，元胞自动机还被用来模拟雪花等枝晶的形成、液态金属材料的凝固结晶过程以及颗粒材料的垮塌现象等。第11页，共77页。1 绪论元

7、胞自动机应用交通科学领域：1986年，M. Cremer和J. Ludwig初次将元胞自动机运用到车辆交通的研究中。随后，元胞自动机在车辆 交通中的应用主要沿着两条主线展开：对城市道路交通流的研究，以Nagel-Schreckenberg模型为代表;对城市交通网络 的研究，以BML模型为代表。另外，80年代以来，计算机水平日新月异的发展为元胞自动机的 应用提供了强有力的支持。因此，在进入上个世纪90年代后，元胞自动机在交通流理论研究领域中得到了广泛的应用。第12页，共77页。1 绪论元胞自动机应用计算机科学与信息学领域：元胞自动机的逻辑思维方法为并行机的发展提供了另一个理论框架。20世纪80年

8、代，T. Toffoli和N.H. Margolus 制造出第一台通用元胞自动机计算机CAM6，其性能可与当时的巨型计算机相比拟，并且其图形显示功能明显优于其他类型的计算机。元胞自动机还被用来研究信息的保存、传递、扩散的过程。除此之外，元胞自动机在图像处理和模式识别中也体现出了其独到的优势 。第13页，共77页。 元胞自动机（Cellular Automata，简称CA）实质上是定义在一个由具有离散、有限状态的元胞组成的元胞空间上，并按照一定的局部规则，在离散的时间维度上演化的动力学系统。元胞自动机的定义：2 元胞自动机的定义和构成第14页，共77页。2 元胞自动机的定义和构成元胞自动机的构成

9、：元胞自动机最基本的组成：元胞、元胞空间、邻居及规则四部分。另外，还应包含状态和时间。 可以视为由一个元胞空间和定义于该空间的变换函数所组成。 第15页，共77页。元胞自动机的构成示意图2 元胞自动机的定义和构成第16页，共77页。元胞 元胞又可称为单元、细胞或基元，是元胞自动机的最基本的组成部分。元胞分布在离散的一维、二维或多维欧几里德空间的晶格点上。 2 元胞自动机的定义和构成第17页，共77页。元胞状态元胞的状态可以是二进制形式，如：（0，1），（生，死），（黑、白）等 ；也可以在一个有限整数集内S内取值：如交通领域的CA模型中，有时元胞状态可在-(Vmax+1)Vmax+1)之间取值。

10、状态参量：严格意义上的CA只能有一个状态参量；但是，在实际应用中，可以具有多个状态参量。2 元胞自动机的定义和构成第18页，共77页。元胞空间元胞在空间中分布的空间格点的集合就是元胞空间。 元胞空间的几何划分元胞空间的边界条件2 元胞自动机的定义和构成第19页，共77页。元胞空间的几何划分理论上，它可以是任意维数的欧几里德空间规则划分。常用的元胞自动机一般是一维和二维的。 二维元胞自动机通常有三种划分方式三角形正方形正六边形 一维元胞自动机的元胞空间只有一种划分 2 元胞自动机的定义和构成第20页，共77页。二维元胞自动机的三种网格划分2 元胞自动机的定义和构成第21页，共77页。网格类型优点

11、缺点三角形拥有相对较少的邻居数目，易于处理复杂边界在计算机的表达与显示不方便，需要转换为四方网格。正方形直观而简单，而且特别适合于在现有计算机环境下进行表达显示不能较好地模拟各向同性的现象正六边形能较好地模拟各向同性的现象，因此，模型能更加自然而真实在表达显示上较为困难、复杂三类网格划分的优缺点对比2 元胞自动机的定义和构成第22页，共77页。元胞空间边界条件理论上，元胞空间是无限的；实际应用中无法达到这一理想条件。常用的边界条件如下：周期型定值型绝热型反射型2 元胞自动机的定义和构成第23页，共77页。元胞空间边界条件周期型边界条件（periodic boundary）定义：周期型是指相对边

12、界连接起来的元胞空间 对一维空间，首尾相接形成一个圆环对二维空间，上下相接，左右相接，而形成一个拓扑圆环面，形似车胎或甜点圈 周期型空间与无限空间最为接近，因而在理论探讨时，常以此类空间作为试验。 2 元胞自动机的定义和构成第24页，共77页。元胞空间边界条件定值型边界条件(Constant Boundary) 定义：所有边界外元胞均取某一固定常量 绝热型边界条件(Adiabatic Boundary) 定义：在指边界外邻居元胞的状态始终和边界元胞的状态保持一致，即具有状态的零梯度。 2 元胞自动机的定义和构成第25页，共77页。定义：在边界外邻居的元胞状态是以边界元胞为轴的镜面反射。 元胞空

13、间边界条件反射型边界条件(Constant Boundary) 2 元胞自动机的定义和构成第26页，共77页。构形(Configuration)定义：构形是在某个时刻，在元胞空间上所有元胞状态的空间分布组合。在数学上，它通常可以表示为一个多维的整数矩阵。2 元胞自动机的定义和构成第27页，共77页。邻居(Neighbor) 冯-诺依曼(Von. Neumann)型定义如下：分别表示邻居元胞的行坐标和列坐标：分别表示中心元胞的行坐标和列坐标：邻居的数目2 元胞自动机的定义和构成第28页，共77页。邻居(Neighbor) 摩尔(Moore)型 定义如下：邻居的数目2 元胞自动机的定义和构成第29

14、页，共77页。邻居(Neighbor) 扩展的摩尔(Moore)型 定义如下：邻居的数目2 元胞自动机的定义和构成第30页，共77页。邻居(Neighbor) 马哥勒斯 (Margolus)型 与前几种邻居的本质区别：以22的元胞块为单元进行处理，而不是向前面几种，对每个元胞分别处理。主要应用领域：格子气流体，颗粒流等Margolus邻居的表现形式和几个演化规则2 元胞自动机的定义和构成第31页，共77页。规则(Rule) 根据元胞当前状态及其邻居状况确定下一时刻该元胞状态的动力学函数，简单讲，就是一个状态转移函数。 称f为元胞自动机的局部映射或局部规则 2 元胞自动机的定义和构成第32页，共

15、77页。时间 元胞自动机中的时间是离散的，是一系列的整数值，是一个无量纲的整数。 若时间步长为dt=1,t=0为初始时刻，则t+1就为下一个时刻。2 元胞自动机的定义和构成第33页，共77页。根据上面对元胞自动机的组成分析，我们可以更加深入地理解元胞自动机的概念。 可以将元胞自动机概括为一个用数学符号来表示的四元组。A:代表一个元胞自动机系统；Ld：代表元胞空间；d：为空间维数；S：是元胞有限的离散的状态集合；N：表示邻域内所有元胞的组合（包括中心元胞在内）；f：是局部转换函数，也就是规则。2 元胞自动机的定义和构成第34页，共77页。1986年，Cremer和Ludwig初次将元胞自动机运用

16、到车辆交通的研究中。交通流的元胞自动机模型大致可分为两大类：研究高速公路交通的模型（以NS模型为代表）；研究城市网络交通的模型（以BML模型为代表）这两类模型是以Wolfram命名的184号模型为基础发展而来的。3 184号规则第35页，共77页。184号模型道路被划分为等距格子，每个格点表示一个元胞；某个时刻，元胞或者是空的，或者被一辆车占据；所有车辆的行进方向都是一致的（如向右）；在每一个时间步内：若第n辆车的前方元胞是空的，则该车可以向前行驶一步；若前面的元胞被另一辆车n+1所占据，即使第n+1辆车在本时间步内离开此元胞，第n辆车也停在原地不动；整个系统采用周期性边界条件以确保车辆数守恒

17、。3 184号规则第36页，共77页。10001110111 110 101 100 011 010 001 000 tt+1Rule 184：3 184号规则第37页，共77页。tt+1t+223+24+25+27=184 3 184号规则第38页，共77页。 作为对184号规则的推广，Nagel和Schreckberg在1992年提出了一个模拟车辆交通的元胞自动机模型，即NS模型（也有人称它为NaSch模型）。时间、空间和车辆速度都被整数离散化道路被划分为等距离的离散的格子，即元胞每个元胞或者是空的，或者被一辆车所占据车辆的速度可以在（0Vmax）之间取值4 NS模型第39页，共77页。在

18、时刻t到时刻t+1的过程中按照下面的规则进行更新：4 NS模型第40页，共77页。在时刻t到时刻t+1的过程中按照下面的规则进行更新：4 NS模型第41页，共77页。更新过程图示：4 NS模型第42页，共77页。更新过程图示：4 NS模型第43页，共77页。更新过程图示：4 NS模型第44页，共77页。更新过程图示：4 NS模型第45页，共77页。边界条件周期性边界条件 在每次更新结束后，我们要监测道路上头车的位置X_lead，如果 X_leadL_road，那么这两车将从道路的另一端进入系统，变为道路上的尾车，并且X_lead = X_lead-L_road, V_last= V_lead。

19、开口边界条件假设道路最左边的元胞对应于X=1，并且道路的入口端包含V_max个元胞，也就是说，车辆可以从元胞 (1, 2, V_max)进入到道路中。在t t+1时刻，当道路上的车辆更新完成后，监测道路上的头车和尾车的位置X_lead和X_last。如果X_lastV_max，则一辆速度为V_max的车将以概率a进入元胞minX_last-V_max, V_max。在道路的出口处，如果X_leadL_road，那么道路上的头车以概率b驶出路段，而紧跟其后的第二辆车成为新的头车。4 NS模型第46页，共77页。与184号模型相比，NS模型的主要改进之处是引入了慢化概率和最大车速不再是1。NS模型

20、虽然具有十分简单的形式，但却可以描述一些实际交通现象。比如NS模型可以模拟出自发产生的堵塞现象以及拥挤交通情况下的时走时停波等。4 NS模型第47页，共77页。堵塞的形成4 NS模型第48页，共77页。时走时停波4 NS模型第49页，共77页。模拟结果：d=0.05 ,p=0.34 NS模型第50页，共77页。模拟结果：d=0.1 ,p=0.34 NS模型第51页，共77页。模拟结果：d=0.2 ,p=0.34 NS模型第52页，共77页。模拟结果：d=0.4 ,p=0.34 NS模型第53页，共77页。模拟结果：d=0.8,p=0.34 NS模型第54页，共77页。模拟结果：d=0.2，p=

21、04 NS模型第55页，共77页。NS模型的衍生和发展以NS模型为基础，研究者主要做了以下三个方面的工作理论方面为了能够模拟出交通实测中的各种现象，他们提出了各种各样的改进模型。如通过改进慢化规则，提出慢启动模型，巡航控制极限(Cruise Control Limit)模型，密度相关的慢启动概率模型；通过改进加速规则，提出Fukui-Ishibashi(FI)模型；以及考虑前车速度效应的模型，三相交通状态模型，舒适驾驶模型等等 。4 NS模型第56页，共77页。NS模型的衍生和发展理论方面他们对元胞自动机模型(主要是NS模型及FI模型)做了相关的解析分析(主要是平均场分析)，加深了人们对元胞自

22、动机模型的理解。引入换道规则用来模拟多车道可以超车的情况，考察了双向交通情况.4 NS模型第57页，共77页。NS模型的衍生和发展应用方面 利用NS模型对大型交通系统的模拟近年来取得了可喜的进展。到目前，NS模型已经应用于美国城市智能交通项目TRANSIMS，杜伊斯堡的内城交通，达拉斯/福斯-华斯地区的交通规划，以及北莱茵-魏斯特伐利亚地区的交通公路网中。4 NS模型第58页，共77页。NS模型及其衍生和发展介于理论和应用之间 利用NS模型操作灵活简单易行的特点，考察了各种交通瓶颈对交通的影响。一方面发现了一些新的物理现象，另一方面又为交通工程建设及交通管理规则的制定提供了合理建议和科学依据。4 NS模型第59页，共77页。/applets/RoadApplet/RoadApplet.html/mm/traffic/ca.htmlhttp:/bolay.de/kai/RoadApplet/http:/smtca.dyns.cx/4 NS模型第60页，共77页。BML模型 1992年，Biham, Middleton, Levine提出了第一个二维交通流元胞自动机模型。 模型定义于一个NN的方形格点的网络上；每