基于matlab的模拟生长卢梦主要内容模拟树生长模拟石墨烯生长树.PPT

基于matlab的模拟生长;主要内容;树生长;function l a=zeros(53,303);a(2,5:80:300)=1; %设置一个零矩阵，并在底部设置凝结核 for i=1:1:50000 %总共产生50000个粒子 t=randperm(301);f=t(1)+1;e=51; %在顶部随机产生粒子 while e52e2f302f2 %限定粒子的运动范围 s=randperm(4);r=rem(s(1),4); %定义随机运动 switch(r) case 0 f=f-1; case 1 f=f+1; case 2 e=e-1; otherwise e=e+1; end if a(e,f+1)+a(e,f-1)+a(e-1,f)+a(e+1,f)0 %凝结判断 a(e,f)=1; %赋值到矩阵 break else continue end end end [x,y]=find(a); %找出矩阵中不为零的元素 scatter(y,x,.);;生长图;模拟石墨烯生长;原理图;function smx L=151; axis([0 L 0 L]); %建立坐标系 x0=floor(L/2)+1; y0=x0; X=sparse(L,L); %建立相对应的矩阵 X(x0,y0)=1; %设定凝结核 R=floor(L/2)-5; for i=1:1:15000 r=rand*2*pi; kx=floor(x0+R*cos(r)); ky=floor(y0+R*sin(r)); %产生随机点 if X(kx,ky)==1 break end while (kx-76)*(kx-76)+(ky-76)*(ky-76)5184 s=randperm(4);t=rem(s(1),4); switch(t) case 0 kx=kx-1; if X(kx,ky)==1 kx=kx+1;continue end case 1 kx=kx+1; if X(kx,ky)==1 kx=kx-1;continue end;;; 结束 谢 谢 ！