Discussion :

[sky-mars]

Bonjour
J'ai un léger problème de compréhension sur Matlab.
Je souhaite tracer l'intercorrelation entre deux signaux mais je sais pas si je dois utiliser "crosscorr" ou bien "xcorr"
J'ai consulté la documentation matlab mais bon j'y comprends pas grand chose.
En utilisant les deux, j'obtiens des résultats complétements différents.
Lequel dois je choisir si je veux mon intercorrelation ?
Par contre pour l'autocorrelation j'ai des schémas tout moches, est-ce normal ?

merci de votre aide

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
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
clear all
n=0:8;
 
%Function
x=n
h=n.*(n<=3)
%Plot of x and h
subplot(1,2,1)
stem(n,x);
xlabel('x(n)')
grid on
subplot(1,2,2)
stem(n,h)
xlabel('h(n)')
grid on
figure(2)
 
rhh=autocorr(h)
rxh=crosscorr(x,h)
rxx=autocorr(x)
rhx=crosscorr(h,x)
figure(2)
subplot(2,2,1)
plot(rxx)
xlabel('auto-Correlation Rxx')
grid on
subplot(2,2,2)
plot(rhh)
xlabel('auto-Correlation Rhh')
grid on
subplot(2,2,3)
plot(rxh)
xlabel('cross-Correlation Rxh')
grid on
subplot(2,2,4)
plot(rhx)
xlabel('cross-Correlation Rhx')
grid on
%Part 2
%h(n-4)=h(n') for 4<n'<7
h1=n.*(n>=4 & n<=7)
figure(3);
stem(n,h1);
xlabel('h(n-3)')
grid on
rhh1=crosscorr(h,h1)
figure(4);
plot(rhh1);
grid on
xlabel('Intercorrelation between h(n) and h(n-3)')

[Jean Dumoncel]

Bonjour,

pour l'inter-correaltion et l'auto-correlation, il faut utiliser xcorr.

crosscorr fait partie de la toolbox GARCH. Je ne crois pas que cette toolbox ait un rapport avec ce que tu cherches à faire. Tu fais du traitement du signal, utilises les fonctions de la toolbox traitement du signal (comme xcorr)

[sky-mars]

Salut

Quand j'utilise xcorr entre x et x pourquoi mon maximum n'est pas 0 ?
On est censé avoir une auto correlation égale à 1 puisque c'est le même signal et sa devrait tendre vers 0 quand Rxx(h) quand h tend vers l'infini non ?

[Jean Dumoncel]

Pourquoi le maximum devrait être égal à 1? Il n'y a pas de raison si tu regardes la formule théorique (dans le help de matlab ou sur le net).

Pour l'histoire de tendre vers 0, ça c'est la théorie, le nombre d'élément dans ton vecteur de sortie est de 2N-1 (si ton vecteur d'entrée est de longueur N). Mais normalement tu devrait tout de même avoir une courbe croissante puis décroissante.

[sky-mars]

Sa devrait être égale à 1 car il est nécessaire de normalisée l'auto/inter correlation si par exemple on a x(t) et 10000x(t) il est nécessaire de normalisé non ?
en faisant

Code :

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

xcorr(x,x,'coeff')

[Jean Dumoncel]

Je dirais que c'est plus utile pour l'intercorrelation, mais oui, tu peux normaliser en utilisant 'coeff'.
Petite remarque :

Code :

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

xcorr(x,x,'coeff');

peut directement s'écrire :

Code :

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

xcorr(x,'coeff');