clc
clear?all
%%%%%%%%%%%%%%%%%%%%%%%%
%myobj=VideoWriter（‘particle.avi‘）;
%writerObj.frameRate=10;
%open（myobj）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
[xxxyyyzzz]=sphere（100）;
%break;
R0=1;
rnbr2=（5*R0）^2.0;

k=10.0;
%K=100;
Lx=40.0;
Ly=20.0;
Dr=0.001;
D0=0.001;
v01=0.1;
dt=0.01;
height=10.0;
Lxa=0.25*Lx;
tg1=height/Lxa;
tg2=-height/（Lx-Lxa）;
%phi=0.5;
%file_t=fopen（‘x_y.dat‘‘r‘）;
%n=round（phi*（Lx*Ly-pi*R*R）/（pi*R0*R0））
n=20;
%break;
%n=50;
%n=round（phi*L^2/（pi*d））
%[x0y0]=textscan（file_t‘%f?%f‘）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
tend=20000.0;
m=round（tend/dt）;

for?ii=1:n
d（ii）=R0;
v0（ii）=v01;
end

phi=sum（pi*d.*d）/（Lx*Ly-Lx*height）
for?ii=1:n
?x0（ii）=Lx*rand（1）;
?y0（ii）=Ly*rand（1）;
?theta0（ii）=2*pi*rand（1）;
?ix（ii）=0;
end

x0

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????5374??2015-09-15?11:29??levy.m

-----------?---------??----------?-----??----

?????????????????5374????????????????????1