| 查看: 406 | 回复: 0 | ||
[求助]
在Matlab上模拟DLA的吸附，请大神帮忙修改下代码，谢谢
|
|
% clc;clear;close all; rand('state',0); set(gcf,'DoubleBuffer','on'); N=256; % 生成状态矩阵大小的控制参数 S=ones(N); % 生成状态矩阵S % S(N/2,N/4:N*3/4)=0; % 设置状态矩阵的初始值 S(N/2:N*3/4,N/2)=0; % 设置状态矩阵的初始值 Ii=imshow(S); % 显示状态矩阵 ti=title(['time = 0',', N=0'],'Fontsize',14,'Fontname','Times New Roman'); % 显示时间 T=0; % 记录时间的参数 for k=1:30000; % 循环计算 xt=4; % 产生粒子的位置 yt=N/2; % 产生粒子的位置 Ss=0; % 控制下面循环是否终止的参数 while Ss<1; % 计算当前粒子的吸附过程 [dx,dy]=choose_D(0.25,0.25,0.25,0.25); % 计算粒子位移量 xt=xt+dx; % 计算下一时刻粒子的位置 yt=yt+dy; % 计算下一时刻粒子的位置 if xt<2.5|xt>N-1.5||yt<2.5|yt>N-1.5; % 判断粒子是否逸出边界 Ss=2; % 把Ss赋值为2,从而可以退出循环 end if [S(xt+1,yt)+S(xt-1,yt)+S(xt,yt+1)+S(xt,yt-1)]<3.5; % 判断近邻是否有粒子 S(xt,yt)=0; % 固定当前运动的粒子 Ss=2; % 把Ss赋值为2,从而可以退出循环 end end set(Ii,'CData',S); % 显示当前状态矩阵 T=T+1; % 累加时间参数 set(ti,'String',['time = ',num2str(T),', N=',num2str(sum(sum(1-S))-[N/4+1])]); % 更新时间参数 pause(0.01); % 暂停一下，显示动画效果 end |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有144人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
点击这里搜索更多相关资源