(完整版)元胞自动机-森林火灾模型MATLAB代码.pdf
版权申诉
5星 · 超过95%的资源 12 浏览量
2022-07-13
23:33:51
上传
评论 1
收藏 64KB PDF 举报 
% 元胞自动机：森林火灾模型
% 规则：
% (1)正在燃烧的树变成空格位；
% (2)如果绿树格位的最近邻居中有一个树在燃烧，则它变成正在燃烧的树；
% (3)在空格位，树以概率 p 生长；
% (4)在最近的邻居中没有正在燃烧的树的情况下树在每一时步以概率f(闪
%? 电)变为正在燃烧的树。
% 参考文献:
% 祝玉学，赵学龙译，<<物理系统的元胞自动机模拟>>, p23
close all;
clc;
clear;
figure;
p=0.3;% 概率 p
f=6e-5;% 概率 f
axes;
rand('state',0);
set(gcf,'DoubleBuffer','on');
% S=round((rand(300)/2+0.5)*2);
S=round(rand(300)*2);
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,1)=R;
C(2:301,2:301,2)=G;
Ci=imshow(C);
ti=0;
tp=title(['T = ',num2str(ti)]);%%时间记录
while 1;
ti=ti+1;
St=Sk; %%St 表示 t 时刻的森林情况
St(Sk==2)=0; % for rule (1)
Su=zeros(302);
Sf=Sk;%%Sf 表示模拟着火的过程
Sf(Sf<1.5)=0;%%只留下着火点
- 1
- 2
- 3
- 4
前往页