下载本文档
版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领
文档简介
1、关于dla生长模型分类:默认栏目2007.4.9 19:14 作者：lionking0083 | 评论：0 | 阅读：420编程是一件漫长而痛苦的事情,尤其是不会用vb和vc的时候，有时候我真的感到很无助，要 是没有matlab,一切真是不敢想象.dla凝聚模型DLA模型采用了在晶格原点上置一个初始粒子作为种子，以该点为圆心(原点)，r为半径 作一个大圆.在该圆附近随机产生一个微粒在圆上作近似于Brown运动的随机行走，即以 1/ 4的概率向上下左右方向行走.若微粒与种子颗粒相碰，则令其附着于种子微粒之上， 与之结合形成凝聚集团；若微粒行走到圆的边界或离开这个圆，令其被边界吸收而消失， 接着再
2、随机地产生第二个微粒并重复以上步骤，直至附着于凝聚集团或被边界吸收.如此 不断地进行，当种子集团长大到一定程度后，可以发现，它就形成了以种子微粒为中心的 无规分叉图形.死伤n个脑细胞后总算搞定.function DLAmodel%Fractal growth modeling by DLAL=150;axis ( 0 L 0 L )set(gcf,color,w)%将图形背景设置为白色 xhalf=floor(L/2)+1;yhalf=xhalf;hold onX=sparse(L,L);X(xhalf,yhalf) = 1; % 设定“种子”位置i,j=find(X);plothdl = p
3、lot(i,j,b.);axis offaxis 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 kxL-1&ky1&ky1&flag =0,N=sum(sum(X(kx+p,ky+p); %统计与随机粒子相邻的“种子团簇”粒子的个数 if N0X(kx,ky) = 1;flag = 1;elseran = rand;%布朗运动控制if ran2/3kx=kx+1;elseif ran2/3ky=ky+1;elseif ran1/3ky=ky-1;endset(plot2,xdata,kx,
温馨提示
- 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
评论
0/150
提交评论