Discussion :

[marc-alex]

Bonjour,

Je cherche à déterminer le temps de retard entre deux signaux a et b. Pour cela j'ai donc voulu calculé différentes corrélations laggées, afin de trouver celle avec le plus haut score, correspondant logiquement au temps de retard. Par conséquent j'ai employé la fonction xcorr de la manière suivante

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
[c ] = xcorr(a,b,'coeff');
[d,e] = max(c);
timelag = (length(e) - f) * 1

Or j'ai testé ces quelques lignes avec b identique à a, mais translaté de 100 pas de temps bien sur, et pourtant je trouve une valeur timelag éqale à 0.
A partir de là deux possibilités, soit j'utilise la mauvaise fonction pour faire un calcul de retard, soit je l'utilise mal.
Merci d'avance pour toute suggestion.

[Invité]

Bonjour,

As-tu regardé cette discussion?

[marc-alex]

Maintenant oui mais pour le moi décalage mesuré dans ce sujet est juste la différence de taille entre les 2 signaux et qui dans ton exemple, est également le décalage entre les deux signaux. Mais dans le cas présent (cf pièce jointe), cela ne fonctionnera pas puisque ce sont deux choses différentes

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
a=load('temp.txt');
a=calc.tair; b = [a(100:end) a(1:99)]; 
subplot(311);
plot(a);
subplot(312);
plot(b);
[c lags] = xcorr(a,b,'coeff');
[d,f] = max(c);
timelag = (length(a) - f) * 1
subplot(313);
plot(c),hold on
plot(f,c(f),'r*')
N = length(c);
plot(round([N/2 N/2]),ylim,'r:')

On obtient un décalage à 0 pour deux signaux identiques mais translatés.