| 查看: 1578 | 回复: 10 | |||
[求助]
在matlab上运行，大神帮忙调试下，是搞分形模拟的DLA
|
|||
|
function DLAmodel %Fractal growth modeling by DLA L = 151 ; axis ( [ 0 L 0 L ]) set (gcf ,‘color’ ,‘w’ ) %将图形背景设置为白色 Xhalf = floor (L/ 2) + 1 ; Yhalf = Xhalf ; hold on X = sparse (L , L) ; X(xhalf , yhalf ) = 1 ; %设定 “种子” 位置 [ i , j ] = find (X) ; plothdl = plot (i , j ,‘b’ ) ; axis off axis square p = - 1 : 1 ; Times = 15000 ; %设定随机产生粒子的个数 for count = 1 ∶Times sita = rand*2*pi ; %生成一个随机点 R = floor (L/2) ; kx = floor (Xhalf + R*cos( sita) ) ; ky = floor ( Yhalf + R*sin( sita) ) ; plot2 = plot (kx , ky ,‘r’ ) ; set (plot2 ,‘marker’ ,‘d’,‘markersize’, 8) drawnow %动态演示随机产生粒子的运动模式 flag = 0 ; [ i , j ] = find (X) ; lr = lengt h (i) ; while kx < L - 1 &ky < L - 1 &kx > 1 &ky > 1 &flag = = 0 , N = sum( sum(X(kx + p , ky + p) ) ) ; %统计与随机粒子相邻的 “种子团簇” 粒子的个数 if N > 0 X(kx , ky) = 1 ; flag = 1 ; else ran = rand ; %布朗运动控制 if ran > 2/3 kx = kx + 1 ; elseif ran < 1/3 kx = kx - 1 ; end ran = rand ; if ran > 2/ 3 ky = ky + 1 ; el seif ran < 1/ 3 ky = ky - 1 ; end set (plot2 ,‘xdata’ , kx ,‘ydata’, ky) ; end end cla ; [ i , j ] = find (X) ; plot (i , j ,‘b.’ ) drawnow end |
回复此楼» 本帖已获得的红花（最新10朵）
» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有135人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab程序在学校的曙光服务器上运行，我用ssh进行操作，每次都是证书问题？
已经有4人回复
- 有关matlab运行的一个问题
已经有1人回复
- linux上运行comsol调用matlab程序的错误
已经有1人回复
- 投稿要求图片为黑白的，如何把matlab运行的程序改成黑白的？
已经有8人回复
- 如何将matlab限定在一个CPU核上运行
已经有12人回复
- 关于matlab GUI打包生产exe文件，并使其在没有安装matlab的计算机上运行的问题。。
已经有4人回复
- matlab 数值积分速度上不去
已经有4人回复
- [视频教程][matlab混合编程视频教学][WMV]
已经有92人回复
- 求教：Matlab在集群上如何实现并行计算？
已经有5人回复
- 请教matlab在linux下的并行问题
已经有18人回复
- MATLAB可以在Android手机上运行了 Android版MATLAB下载使用
已经有33人回复
- asp.net调用matlab的制作网站，在vs2010上正常运行，但在服务器IIS无法显示计算的结果
已经有3人回复
- MATLAB中的混沌信号仿真程序
已经有52人回复
- 【教程】应用MATLAB建模与仿真
已经有227人回复
- 同样的MATLAB程序，为何在别人的电脑上运行不了
已经有21人回复
zhangzhiguang
至尊木虫 (文坛精英)
- 应助: 121 (高中生)
- 金币: 53088.8
- 红花: 24
- 沙发: 222
- 帖子: 26519
- 在线: 881.4小时
- 虫号: 492270
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
★
csgt0: 金币+1, 谢谢 2013-06-28 11:30:36
csgt0: 金币+1, 谢谢 2013-06-28 11:30:36
| 这是文件,防止复制时字体变化 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : DLAmodel.m
2013-06-27 12:36:53, 1.34 K
9楼2013-06-27 12:37:23
 |
2楼2013-05-28 09:16:39
zhangzhiguang
至尊木虫 (文坛精英)
- 应助: 121 (高中生)
- 金币: 53088.8
- 红花: 24
- 沙发: 222
- 帖子: 26519
- 在线: 881.4小时
- 虫号: 492270
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
3楼2013-05-28 13:43:44
4楼2013-06-04 08:44:17
5楼2013-06-26 20:37:06
zhangzhiguang
至尊木虫 (文坛精英)
- 应助: 121 (高中生)
- 金币: 53088.8
- 红花: 24
- 沙发: 222
- 帖子: 26519
- 在线: 881.4小时
- 虫号: 492270
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
6楼2013-06-27 09:39:22
7楼2013-06-27 11:13:37
zhangzhiguang
至尊木虫 (文坛精英)
- 应助: 121 (高中生)
- 金币: 53088.8
- 红花: 24
- 沙发: 222
- 帖子: 26519
- 在线: 881.4小时
- 虫号: 492270
- 注册: 2008-01-08
- 性别: GG
- 专业: 基础物理学
【答案】应助回帖
★ ★
csgt0: 金币+2, 应助指数+1, 谢谢 2013-06-28 11:30:30
csgt0: 金币+2, 应助指数+1, 谢谢 2013-06-28 11:30:30
|
这是修改后的代码,改动不大,还有运行后的图 function DLAmodel %Fractal growth modeling by DLA L = 151 ; axis ([ 0 L 0 L ]) set (gcf,'color','w' ) %将图形背景设置为白色 Xhalf = floor (L/2) + 1 ; Yhalf = Xhalf ; hold on X = sparse (L,L) ; X(Xhalf,Yhalf ) = 1 ; %设定 “种子” 位置 [i,j] = find (X) ; plot(i,j,'b' ) ; axis off axis square p = - 1 : 1 ; Times = 15000 ; %设定随机产生粒子的个数 for count = 1 :Times sita = rand*2*pi ; %生成一个随机点 R = floor (L/2) ; kx = floor (Xhalf + R*cos( sita) ) ; ky = floor (Yhalf + R*sin( sita) ) ; plot2 = plot (kx , ky ,'r') ; set (plot2 ,'marker' ,'d','markersize', 8) drawnow%动态演示随机产生粒子的运动模式 flag = 0 ; [ i,j] = find (X) ; lr = length (i) ; while kx < L - 1 &&ky < L - 1 &&kx > 1 &&ky > 1 &&flag == 0 , N = sum( sum(X(kx+p ,ky + p))) ; %统计与随机粒子相邻的 “种子团簇” 粒子的个数 if N > 0 X(kx , ky) = 1 ; flag = 1 ; else ran = rand ; %布朗运动控制 if ran > 2/3 kx = kx + 1 ; elseif ran < 1/3 kx = kx - 1 ; end ran = rand ; if ran > 2/3 ky = ky + 1 ; elseif ran < 1/ 3 ky = ky - 1 ; end set (plot2 ,'xdata' , kx ,'ydata', ky) ; end end cla ; [ i,j ] = find (X) ; plot (i,j,'b.' ) drawnow end  untitled.jpg |
8楼2013-06-27 12:30:12
送红花一朵
 |
10楼2013-09-11 11:23:11
