元胞自动机程序

! 元胞自动机程序（生命永不停止） USE MSFLIB PARAMETER IR=400,JR=400,NMAX=10000 !NMAX-随机产生的生命种子 INTEGER IS(0:1001,0:1001),IS1(0:1001,0:1001),ISN(1:8), TMAX, NUM !IS-基体的二维数组 w rite (*,*) 'PLEASE INPUT LOOP(TMAX)' READ (*,*) TMAX ISEED=RTC() IS=15 !"死"的状态，基体为白色 !赋予生命的种子，"活"的状态 1 DO I=1, NMAX IX0=IR*RAN(ISEED)+1 JY0=JR*RAN(ISEED)+1 IS(IX0,JY0)=1 END DO IS1=IS !EXECUTE THE RULE DO T=1,TMAX !边界条件 IS(0,0:JR+1)=IS(R, 0:JR+1) IS(IMAX+1, 0:JR+1)=IS(1, 0:JR+1) IS(0:IR+1,0)=IS(0:IR+1,JMAX) IS(0:IR+1,JR+1)=IS(0:IR+1,1) !搜索生命存在的位置 DO IX=1,IR DO JY=1,JR !判断邻居状态 ISN=(/IS(IX-1,JY-1),IS(IX-1,JY),IS(IX-1,JY+1),IS(IX,JY-1) & ,IS(IX,JY+1),IS(IX+1,JY-1),IS(IX+1,JY),IS(IX+1,JY+1)/) NUM=COUNT(ISN.EQ.1) !赋予生存的条件 IF((IS(IX,JY)==15.AND.NUM==3).OR.(IS(IX,JY)==1.AND & .(NUM==3.OR.NUM==2))) THEN IS1(IX,JY)=1 ELSE IS1(IX,JY)=15 END IF !画图 ISRE=SETCOLOR(IS1(IX,JY)) ISRE=SETPIXEL(IX,JY) END DO END DO IS=IS1 END DO END 元胞自动机—产生和发展. 四个阶段: 1940s 诞生：Von Neumann 自我复制机. 1960-70s 起步：JH.Conway 生命游戏. 1980s 理论研究：S.Wolfram CA 分类. 1980-90s 应用：HPP-FHP 格子气自动机. C.Langton N.Packard 人工生命 元胞自动机（Cellu lar Au tomata，简称CA，也有人称其为细胞自动机、点格自动机、分子自动机或单元自动机）是一种建立在离散的时间和空间上的动力学系统。散布在规则格网（Lattice Grid）中的每一元胞（Cell）取有限的离散状态，遵循同样的作用规则，依据确定的局部规则作同步更新。大量元胞通过简单的相互作用而构成动态系统的演化。 与一般的动力学模型不同，元胞自动机不是由严格定义的物理方程或函数确定，而是由一系列模型构造的规则构成，凡是满足这些规则的模型都可以算做是元胞自动机模型。因此，元胞自动机是一类模型的总称，或者是一个方法框架。其特点是时间、空间、状态都是离散的，每个变量只取有限多个状态，且其状态改变的规则在时间和空间上都是局部的。 元胞自动机是一种对具有局域连通性的格点，应用局部（有时为中等范围）确定性或概率性的转换规则来描述...