matlab下三维dla模型模拟,Matlab下三维DLA模型模拟.doc

模型模拟,Matlab下三维DLA模型模拟.doc"/>

matlab下三维dla模型模拟,Matlab下三维DLA模型模拟.doc

Matlab下三维DLA模型模拟　2007-01-11 19:18

分类：science

字号： 大 中 小小

function dla3dv5(Nsum,Wstep) %定义dla函数，Nsum为所生成絮体包含的颗粒数，Wstep为计算过程中所采取的步长

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%本程序内变量的定义 %%

%%radius为颗粒半径，release为起始释放半径 %%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

radius=0.5; %%颗粒半径

release=2; %%起始释放半径

L=200;

Xhalf=floor(L/2);

Yhalf=Xhalf;

Zhalf=Xhalf;

n=1; %粒子计数

N(1)=1;

N(Xhalf)=0;

p(1,:)=[Xhalf Yhalf Zhalf];

szpoints=zeros(L,L,L); %%网格点阵

szpoints(Xhalf,Yhalf,Zhalf)=1; %%种子位置标志

%%释放初始粒子

theta=2*pi*rand;

gama=pi*rand;

M=p(1,:)+Wstep*[cos(theta) sin(theta) cos(gama)];

while n(release+15)^2 %%判断是否逃逸

theta=2*pi*ra