元胞自动机：森林火灾模型代码及讲解.docVIP

PAGE PAGE 1 元胞自动机：森林火灾模型% 规则：%? (1)正在燃烧的树变成空格位；%? (2)如果绿树格位的最近邻居中有一个树在燃烧，则它变成正在燃烧的树；%? (3)在空格位，数以概率p生长；%? (4)在最近的邻居中没有正在燃烧的树的情况下树在每一时步以概率f(闪%???? 电)变为正在燃烧的树。% 参考文献:% 祝玉学，赵学龙译，物理系统的元胞自动机模拟, p23close all;clc;clear;figure;p=0.3;?? % 概率pf=6e-5;? % 概率faxes;rand(state,0);set(gcf,DoubleBuffer,on);% S=round((rand(300)/2+0.5)*2);S=round(rand(300)*2);% \copyright: zjliu% Authors email: hepenghpen@163.com Sk=zeros(302);Sk(2:301,2:301)=S;% 红色表示正在燃烧(S中等于2的位置)% 绿色表示绿树(S中等于1的位置)% 黑色表示空格位(S中等于0的位置)C=zeros(302,302,3);R=zeros(300);G=zeros(300);R(S==2)=1;G(S==1)=1;C(2:301,2:301,