| 查看: 278 | 回复: 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 |
回复此楼» 猜你喜欢
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有218人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收10人
已经有0人回复
-
天津城建大学理学院化学、环境工程、材料物理等方向接收调剂，10人
已经有0人回复
点击这里搜索更多相关资源
5