Matlab动画模拟分子布朗运动的示例

Matlab动画模拟分子布朗运动的示例

%Brownian motion 
clf; 
n=20; 
s=0.02; 
x = rand(n,1)-0.5; 
y = rand(n,1)-0.5; 
h = plot(x,y,'.'); 
axis([-1 1 -1 1]) 
axis square 
grid off 
set(h,'EraseMode','xor','MarkerSize',18) 
grid on; 
title('Press Ctl-C to stop'); 
while 1 
 drawnow 
 x = x + s*randn(n,1); 
 y = y + s*randn(n,1); 
 set(h,'XData',x,'YData',y) 
end

来源：http://www.bime.ntu.edu.tw/~dsfon/Mechanism2/%E5%8B%95%E4%BD%9C%E7%9A%84%E5%90%88%E6%88%90.PDF