北大讲义空间分析与机助制图元胞自动机.ppt

第五章 基于元胞自动机的林火蔓延模拟 引线 中央电视台关于美国和国内森林大火的报道 美国相关机构进行林火蔓延模拟研究的报道 提纲 1、元胞自动机的概念 2、元胞自动机的构成 3、元胞自动机的定义 4、元胞自动机的特征 5、 Rothermel林火模型 6、基于元胞自动机的林火蔓延模拟 元胞自动机的定义 元胞自动机是物理学、数学、计算机科学和生物学相交叉的结晶 物理学：一个离散的无穷维的动力学系统 数学：一个时空离散的数学模型，描述连续对象的偏微分方程的的对立体 计算机科学：新兴的人工智能、人工生命的分支 生物学：生命现象的一种描述 物理学的定义 具有离散、有限状态的元胞组成的元胞空间上，并按照一定局部规则，在离散的时间维上演化的动力学系统 基于集合论的定义 设d代表空间维数，k代表元胞的状态，并在一个有限的集合S中取值，r代表元胞邻居半径，Z是整数集，表示一维空间，t代表时间 在一维空间上考虑元胞自动机，即假定d=1。那么整个元胞空间就是在一维空间，将整数集Z上的状态集S的分布计为Sz。元胞自动机的动态演化就是在时间上状态组合的变化，计为： F：StZ St+1Z 这个动态演化又由各个元胞的局部演化规则f所决定，这个局部函数f通常又称为局部规则。对于一维空间，元胞及其邻居可计为S2r+1，局部函数则可以计为：f: St 2r+1 St+1 对于局部规则f 来讲，函数的输入、输出集均为有限集合,实际上，它是一个。 对于元胞空间内的元胞，独立施加上述局部函数，则可得到全局的演化：F(Cit+1)=f(Cti-r,…Cit,…，Cti+r)， Cit+1表示位置i处的元胞。上述即为一元胞自动机。 基于拓扑学的定义 假定维数d为1,设S为k个符号的有限集，Z为整数全体的集合，称Z到S的映射的全体SZ为构形空间。显然，SZ就是用S中的符号组成的双侧无限的符号序列的全体，即一维元胞自动机的所有构形的集合。称a=(…a-1,a0,a1,…)为构形空间中的点。 在中引进任意两点x和y之间的距离d(x,y)=Σβ(xi,yi)2-|i| 其中，当xi=yi时β(xi,yi)＝0，当xiyi时，β(xi,yi)1。然后，在SZ中可以建立起开、闭、紧等拓扑概念。 在SZ中定义移位算子β为β(xi)=xi+1,i€Z。若连续映射F： SZ—SZ与β可交换，即F β＝ βF，或对于任意x € SZ,有F(β(x))= β(F (x))，则称F为元胞自动机。 任意维空间： Sz计为SZd,S2r+1计为S2r+1d 元胞自动机的构成 元胞、元胞空间、邻居和规则四部分。元胞自动机可视为一个元胞空间和定义于该空间的变化函数所组成。 元胞自动机的构成 1. 元胞 元胞自动机的最基本组成部分。它分布在离散的一维、二维或多维欧几里德空间的晶格点上。 2. 状态 状态可以是{0,1}的二进制形式，或是{s0,s1,s2,…，si}整数形式的离散集。严格意义上，元胞只能有一个状态变量，但在实际应用中，往往将其进行扩展。 元胞自动机的构成（续） 3. 元胞空间 元胞所分布在的空间网点集合就是元胞空间 A、元胞空间的几何划分 任意维数的欧几里德空间规则划分。对于一维元胞自动机，元胞空间划分只有一种。而高维的元胞自动机，元胞空间的划分则可能有多种形式。对于常见的二维自动机，元胞空间通常可按三角形、四边形或六边形三种网格排列。 元胞自动机的构成（续） 元胞自动机的构成（续） 三角网格拥有较少的邻居数目，这在某些时候很有用。缺点是计算机的表达与显示不方便。 四边形网格直观简单，特别适合于计算机环境下进行表达显示。 六边形网格能较好的模拟各向同性的现象，因此，模型能更加自然而真实。其缺点同正三角网格一样，在表达显示上较为困难和复杂。 元胞自动机的构成（续） B、边界条件 理论上，元胞空间在各个维向上是无限延展的。实际应用过程中，无法在计算机上实现这一理想条件。 C、构形 在某个时刻，在元胞空间上所有元胞状态的空间分布组合。在数学上，它通常可以表示为一个多维的整数矩阵。 元胞自动机的构成（续） 4. 邻居 元胞和元胞空间只表示了系统的静态成分，为了将动态引入系统，必须加入演化规则。这些规则是定义在局部空间范围内的，即一个元胞下一时刻的状态决定于本身的状态和它的邻居元胞的状态。因此，在指定规则之前，必须定义一定的邻居规则，明确哪些元胞属于该元胞的邻居。在一维元胞自动机中，通常以半径r来确定邻居，距离一个元胞r内的所有元胞都属于该元胞的邻居。二维元