matlab levy 飞行程序

matlab levy 飞行程序，可以有广泛的应用

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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

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

----------- ---------  ---------- -----  ----

                 5374                    1