网站设计师,网站上传的工具,安徽建站管理系统价格,展览设计制作公司clc
clear
n 300; % 定义表示森林的矩阵大小
Plight 5e-6; Pgrowth 1e-2; % 定义闪电和生长的概率
UL [n,1:n-1]; DR [2:n,1]; % 定义上左#xff0c;下右邻居
vegzeros(n,n); % 初始化表示森林的矩阵
imh ima…
clc
clear
n 300; % 定义表示森林的矩阵大小
Plight 5e-6; Pgrowth 1e-2; % 定义闪电和生长的概率
UL [n,1:n-1]; DR [2:n,1]; % 定义上左下右邻居
vegzeros(n,n); % 初始化表示森林的矩阵
imh image(cat(3,veg,veg,veg)); % 可视化表示森林的矩阵
% veg 空地为0 着火为1 树木为2
for i1:3000 sum (veg(UL,:)1) (veg(:,UL)1)(veg(:,DR)1) (veg(DR,:)1); % 计算出所有格子有几个邻居是着火的
% 根据规则更新森林矩阵是否树是否树-是否着火的树是否新生的树0-1运算veg 2*(veg2) - ( (veg2) (sum0 |(rand(n,n)Plight)) ) 2*((veg0) rand(n,n)Pgrowth) ; set(imh, cdata, cat(3,(veg1),(veg2),zeros(n)) )drawnow % 可视化表示森林的矩阵
endclc
clear all;
n200;
Sezeros(n);
Zzeros(n)
Se(n/2-2:n/22,n/2-2:n/22)1;
Chimagesc(cat(3,Se,Z,Z));
axis square
Sdzeros(n2);
while(1)Sd(2:n1,2:n1)Se;sumSd(1:n,2:n1)Sd(3:n2,2:n1)Sd(2:n1,1:n)Sd(2:n1,3:n2);Semod(sum,2);set(Ch,cdata,cat(3,Se,Z,Z))pause(0.05)
endfigureclc
clear;
n200;
p0.4;
zzeros(n)
Serand(n)p;
Sdzeros(n2);%矩阵初始化
Phimage(cat(3,Se,z,z));%初始可视化
while(1)Sd(2:n1,2:n1)Se;%SumSd(1:n,2:n1)Sd(3:n2,2:n1)Sd(2:n1,1:n)Sd(2:n1,3:n2)Sd(1:n,1:n)Sd(3:n2,1:n)Sd(1:n,3:n2)Sd(3:n2,3:n2);
%邻居之和邻居中生的元胞的数目for i1:nfor j1:nif Sum(i,j)3||(Sum(i,j)2Se(i,j)1)%生的条件Se(i,j)1;else Se(i,j)0;endend
end
set(Ph,cdata,cat(3,Se,z,z))%更新可视化drawnow
end
