Discussion :

[naomie]

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:

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')

Merci d'avance,

Naomie