元胞自动机也是建模中应用潜力很大的一种仿真方法。
即通过设置元胞的变化规则构成整体变化。
这里是一个森林火灾的仿真,可调参数有雷击(自燃)概率和重新生长的概率,以及仿真规模。
特别之处在于设置了对角引燃的概率
,可以使得燃烧范围呈现圆形扩散而不是菱形。
演示图像:
MATLAB脚本代码:
%%
%元胞自动机,燃烧森林模型
% 2016/10/8
% 规则:
% (1)正在燃烧的树变成空格位;
% (2)如果绿树格位的最近邻居中有一个树在燃烧,则它变成正在燃烧的树
% (对角线以sqrt(1/2)-0.5的概率引燃);
% (3)在空格位,树以概率p生长;
% (4)在最近的邻居中没有正在燃烧的树的情况下树在每一时步以概率f(闪
% 电)变为正在燃烧的树。
% 参考文献:
% 祝玉学,赵学龙译,<<物理系统的元胞自动机模拟>>, p23
%%
close all;clc;clear;
figure;
s=300; %树林大小
p=0.02; % 概率p
f=1e-5;