第7章 元胞自动机.pptVIP

参考书：《分形算法与程序设计》 * 从逻辑运算谈起 7.1 逻辑异或 本 行：0 0 0 1 1 0 1 1 下一行： 0 1 1 0 7.2 一维元胞自动机 元胞按等间隔方式分布在一条向两侧无限延伸的直线中，称为一维元胞自动机。 本 行：001 100 其他 下一行： 1 1 0 7.2 一维元胞自动机 变量：i ，j（循环变量） y1[ ]（本行元胞数组） y2[ ]（下一行元胞数组） 标题： 一维元胞自动机 参数： n （水平方向最大值） m（垂直方向最大值） 函数：SPost(x, y, color) （画点） 算法：cell_one BEGIN FOR i=1 TO m FOR j=1 TO n IF y1[j-1]=0 AND y1[j]=0 AND y1[j+1]=1 THEN y2[j]=1 ELSE IF y1[j-1]=1 AND y1[j]=0 AND y1[j+1]=0 THEN y2[j]=1 ELSE y2[j]=0 ENDIF IF y2[j]=1 THEN SPost(x, y, color) ENDIF ENDFOR ENDFOR END 7.3 二维元胞自动机 在一个二维网格中，如果抛下一粒种子（元胞着色），然后考察一下种子身边的格子中的元胞状态会发生什么事情。给一个规则，即每一个格子的状态，由其周围的八个格子的状态（0或1）来决定，如果它周围八个格子中的状态值相加为奇数时，则此格子下一个状态为1；如果它周围八个格子中的状态值相加为偶数时，则此格子下一个状态为0。就这样一步一步演化下去，会看到图案。 7.3 二维元胞自动机 变量：i ，j（循环变量） a[ ][ ]（本步元胞数组） b[ ][ ]（下一步元胞数组） 标题： 二维元胞自动机 参数： n （元胞的范围） 函数：SPost(x, y, color) （画点） 算法：cell_two BEGIN FOR i=1 TO n FOR j=1 TO n k=a[i-1][j-1]+a[i-1][j]+a[i-1][j+1] +a[i][j-1]+a[i][j+1] +a[i+1][j1] +a[i+1][j]+a[i+1][j+1] IF k mod 2=1 THEN b[i][j]=1 ELSE b[i][j]= 0 ENDIF ENDFOR ENDFOR END 7.4 分形演化的DLA模型 自然界中有很多种这样的生长集团 7.4 分形演化的DLA模型 变量：r0（凝聚半径） r（释放粒子的半径） rmax（粒子逃逸半径） cyt,cyb,cyl,cyr （游走粒子周围4点） 标题： DLA模型 参数： BAX ，BAY（基点坐标） m（垂直方向最大值） 函数：SPost(x, y, color) （画点） getpixel(x,y)（获取点） read( )（随机函数） sgn( ) （符号函数） 算法：DLA BEGIN r0=5 WHILE(n1000) r=r0+5 rmax=r0+5 rx=(2*r+1)*read()-r rv=r-abs(rx) ry=rv*sgn(read()-0.5) x=rx y=ry FOR( ; ;) xb=x yb=y dist=abs(x)