Bonjour,
Voici mon probleme: je cherche a calculer la MSD (mean square displacement) d'un signal a 1 dimension en appliquant la formule ci joint en attachement.
Mais je pense qu'il y a une faute dans mon code, car si je l'applique a des nombres aleatoires avec l'aide de la fonction rand, je ne trouve pas une relation lineraire.
Voici le code:
Merci d'avance,
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23 clc clear all close all d=rand(5000,1); t=1:1:5000; GN=size(t); GN=GN(2); for n=1:GN temp = 0; for i=1:GN-n temp = temp+(d(i+n)-d(i))^2; end msd(n) = 1/(GN-n)*temp; end subplot(2,1,1) plot(t,d) title('Random motion') subplot(2,1,2) plot(t,msd) title('MSD of a random motion')
Naomie
Partager