Mouvement brownien avec plusieurs molécules/particules [Matlab]

Question

Je veux modéliser le mouvement brownien avec plusieurs molécules/particules et les animer. Jusqu'à présent, j'ai trouvé cette solution de code, mais je suis incapable d'obtenir le résultat souhaité. Ci-dessous est mon code

N = 500; % number of samples 
tau = .1; % time interval in seconds 
D = 10; % diffusion coefficient 
NumMolecules = 500; % number of moelcues to be released 
k = sqrt(2*D*tau); % scaling factor 
dx = k * randn(N,NumMolecules); 
dy = k * randn(N,NumMolecules); 
x = cumsum(dx); 
y = cumsum(dy); 
h = gscatter(x(1,:),y(1,:),1:NumMolecules,'b'); 
for k = 2:N 
    for p = 1:NumMolecules 
     h(p).XData = x(k,p); 
     h(p).YData = x(k,p); 
    end 
    drawnow limitrate 
end 
drawnow 

Quelqu'un peut-il me dire pourquoi je suis incapable d'obtenir les résultats souhaités?

Modifier

Résultat souhaité: Pour animer le mouvement des molécules en utilisant « drawnow » fonction

Answer 1

vous avez h(p).YData = x(k,p); changement à h(p).YData = y(k,p);

Il complote les mêmes numéros les uns contre les autres droits maintenant, c'est pourquoi vous obtenez une ligne droite.