matlab布朗运动完整答案

模拟布朗运动

在二维空间中产生一个花粉（用球代替），设定恒定速度，给定初始位置。

每个时间步随机给定一个方向，计算下一个时间点的位置。

判断是否结束条件（如边界、最大步数）

t=1;x=100;y=100;vx=17;vy=0;

for k=1:40

p=2*pi*rand(1,1)

vx=17*cos(p);

vy=17*sin(p);

x=x+vx*t;

y=y+vy*t;

line('ydata',y,'xdata',x,'Color' ,[1 0 0], 'Marker' ,'.' , 'MarkerSize' ,12, 'EraseMode' , 'non');

plot(x,y);

axis([0 200 0 200])

if(x>200||y>200||x<0||y<0)

break

end

hold on

pause(0.4);

end