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

Matlab下三维DLA模型模拟 2007-01-11 19:18 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;