最新基于一维交通流元胞自动机计算机编程模拟.docx

基于一维交通流元胞自动机计算机编程模拟【摘　要】主要介绍了元胞自动机的特点及定义,元胞自动机在道路交通流模拟中的应用,设计了交通元胞自动机的结构,建立了运用元胞自动机的交通流模拟模型, 运用C语言编程实现模拟,并利用Matlab得到流率与交通密度的关系;探讨了元胞自动机在道路交通流模拟中应用前景。【关键词】一维交通流；元胞自动机; C语言编程;模拟一、引言1992年,Nagel-Schreckenberg采用元胞自动机模拟公路交通流,从此以后,相继有许多研究机构开始了交通元胞自动机的研究,并已有许多研究成果及相关模拟软件产生。元胞自动机(Cellular Automata,简称CA)是一种时间和空间都离散的动力系统模型。散布在栅格(Lattice Grid)中的每一元胞(Cell)取有限的离散状态,遵循同样的作用规则,依据确定的局部规则作同步更新,大量元胞通过简单的相互作用而构成动态系统的演化。不同于一般的动力学模型,元胞自动机不是由严格定义的物理方程或函数确定,而是用一系列模型构造的规则构成。凡是满足这些规则的模型都可以算作是元胞自动机。因此,元胞自动机是一类模型的总称,或者说是一个方法框架。其特点是时间、空间、状态都离散,每个变量只取有限多个状态,且其状态改变的规则在时间和空间上都是局部的。二、元胞自动机介绍１、元胞自动机在元胞自动机的发展过程中,许多学者构造了不同的模型。本文主要介绍交通流模拟常用的初等元胞自动机,又称Wolfram元胞自动机。初等元胞自动机(Elementary Cellular Automata,简称E-CA)是状态集S只有两个元素{s1,s2},即状态个数k=2,邻居半径r=1的一维元胞自动机。2、元胞自动机的构成元胞自动机最基本的组成是元胞、元胞空间、邻居及规则四部分。即元胞自动机是由一个元胞空间和定义于该空间的变换函数所组成的。（１）元胞元胞又可称为单元或基元,是元胞自动机的最基本的组成部分。元胞分布在离散的一维或二维欧几里德空间的晶格点上。（２）状态状态可以是{0,1}的二进制形式。或是{s0,s1,…,sk}整数形式的离散集。对于一维元胞自动机,定义元胞的两个状态,即{0,1}。（３）元胞空间元胞分布的空间网点集合就是所谓的元胞空间。元胞空间的几何划分:理论上,它可以是任意维数的欧几里德空间规则划分。目前研究多集中在一维和二维元胞自动机上。对于一维元胞自动机,元胞空间的划分只有一种。边界条件:在理论上,元胞空间在各维上是无限的。但是在实际应用过程中,无法在计算机上实现这一理想条件,因此,需要定义不同的边界条件。构形:构形是指某个时刻,元胞空间上所有元胞状态的空间分布组合。通常,在数学上,它可以表示为一个多维的整数矩阵。（４）邻居以上的元胞及元胞空间只表示了系统的静态成分,为将动态引入系统,必须加入演化规则。在元胞自动机中,这些规则是定义在局部空间范围内的,即一个元胞下一时刻的状态取决于其本身状态和它的邻居元胞的状态。因而,在指定规则之前,必须定义一定的邻居规则,明确哪些元胞属于该元胞的邻居。在一维元胞自动机中,通常以半径来确定邻居,距离一个元胞内的所有元胞均被认为是该元胞的邻居。3、规则根据元胞当前状态及其邻居状况确定下一时刻该元胞状态的动力学函数,即一个状态转移函数。将一个元胞的所有可能状态连同负责该元胞的状态变换的规则一起称为一个变换函数。这个函数构造了一种简单的、离散的空间/时间的局部物理成分。４、时间元胞自动机是一个动态系统,它在时间维上的变化是离散的,即时间t是一个整数值,而且连续等间距。假设时间间距dt=1,若t=0为初始时刻,t=1为其下一时刻。在上述转换函数中,一个元胞在t+1时刻的状态只决定于t时刻的该元胞及其邻居元胞的状态。三、一维交通元胞自动机的设计１、元胞空间的离散一维交通元胞自动机适宜于单车道交通流的模拟,为此以车辆行驶时所占道路长度(包括安全间距)l为单位,将长度为L的道路均匀划分成连续的栅格,构成元胞空间,形成一列元胞。２、模拟时间的离散假定模拟总时间为T,在该时间空间上,以时间间距dt推进模拟时间,当dt=1时,从t=0初始时间开始,依次顺序加1,推进模拟时间,直至模拟结束。3、状态定义元胞状态为{0,1},其中1表示该元胞为有车状态,0表示该元胞为无车状态。且车辆的运行方向是相同的。4、邻居对于单车道交通流的模拟,可假设元胞t+1的状态仅受到其最近的邻居元胞的状态及该元胞t时刻状态的影响。因此其邻居半径取为r=1。5、动态进化规则一维交通元胞自动机状态更新规则实际上是Wolfram规则的扩展,统称为规则184。(1)Wolfram规则根据实际的交通流变化,对于实际的单车道交通流模拟而言,可设计初等元胞自动机的一个规则如右表所示 (黑色方块代表1,白色方块代表