Intercorrelation. Zero padding. Obtenir une fenetre de plot plus grande.

Version imprimable

Bonjour tout le monde,

Comme decrit dans le titre je suis en train de faire de l'intercorrélation et tout marche plutot bien.:mrgreen:

J ai en réalité deux petits soucis : :aie:

1. Je pense qu'il y a moyen de reduire le code mais je ne vois pas trop comment etant donné les 'if'.

2. Lorsque je 'plot' mon intercorrélation finale ('Z' dans le code), il me met (logiquement) ma courbe sur 1024 points or je voudrais la voir sur 4096 points avec des zeros autour...8O
Ce qui implique une intercorrelation par segments...('1-->1024' ; '1025 --> 2048' ; '2049 --> 3072' ; '3073 --> 4096'). :?

Voila, j aimerais savoir si quelqu'un a une idée pour avoir ma fenetre complete sur 1 plot avec mes 4096 points...et s'il y a chevauchement (exemple : de 1000 points à 1256 points) si cela pourrait se gerer facilement...:lol::lol:

Merci d avance pour vos reponses.

J espere avoir ete assez clair

Voici mon code :
Code:

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 51 52 53 54 55 function balayageIntercorrelation t = [0:4096]; % de 0 à 1000 points sig = sin(t/(2*pi)); % Original signal, a sine wave % Signal servant à l'intercorrelation (sur 1024 points obligatoirement) sig2 = zeros(1, 1024); sig2(1,1:256)= 0.000075*sig(1,1:256); % Signal à intercorreler (taille sans zero padding : 256 points obligatoirement) sig3 = zeros(1,4096); sig3(1,1100:1356)=sig(1,1100:1356); % Boucle pour le premier segment (1 --> 1024) % for N=1:1024:4096; for i=1:1024; if sig3(1,i) ~= 0 X = fft(sig2);% fft du signal de reference Y = fft(sig3(1:1024));% fft du signal à intercorreler end end % Boucle pour le deuxieme segment for i=1025:2048; if sig3(1,i) ~= 0 X = fft(sig2);% fft du signal de reference Y = fft(sig3(1025:2048));% fft du signal à intercorreler end end % Boucle pour le 3e segment for i=2049:3072; if sig3(1,i) ~= 0 X = fft(sig2);% fft du signal de reference Y = fft(sig3(2049:3072));% fft du signal à intercorreler end end % Boucle pour le 4e segment for i=3073:4096; if sig3(1,i) ~= 0 X = fft(sig2); % fft du signal de reference Y = fft(sig3(3073:4096)); % fft du signal à intercorreler end end % Calcul de l'intercorrelation Z = real(ifft(conj(X).*Y)); subplot(2,4,1), plot(sig), title('sig') subplot(2,4,2), plot(sig2), title('sig2') subplot(2,4,3), plot(sig3), title('sig3') subplot(2,4,7), plot(Z), title('fft')

Pour répondre au point 1, il y a effectivement plus simple puisque tu ne dois faire ton calcul que si tu trouves au moins une valeur de sig3 différente de zéro sur le segment :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
% Boucle pour le premier segment (1 --> 1024)
if any(sig3(1:1024)~=0)
    X = fft(sig2);% fft du signal de reference
    Y = fft(sig3(1:1024));% fft du signal à intercorreler
end
 
% Boucle pour le deuxieme segment
if any(sig3(1025:2048)~=0)
    X = fft(sig2);% fft du signal de reference
    Y = fft(sig3(1025:2048));% fft du signal à intercorreler
end
 
% Boucle pour le 3e segment
if any(sig3(2049:3072)~=0)
    X = fft(sig2);% fft du signal de reference
    Y = fft(sig3(2049:3072));% fft du signal à intercorreler
end
 
% Boucle pour le 4e segment
if any(sig3(3073:4096)~=0)
    X = fft(sig2); % fft du signal de reference
    Y = fft(sig3(3073:4096)); % fft du signal à intercorreler
end

Et même en fait, pour gagner encore du temps, puisque si une valeur est différente dans le dernier segment ce sont ces valeurs-là qui sont importantes :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
if any(sig3(3073:4096)~=0)
    X = fft(sig2); % fft du signal de reference
    Y = fft(sig3(3073:4096)); % fft du signal à intercorreler
elseif any(sig3(2049:3072)~=0)
    X = fft(sig2);% fft du signal de reference
    Y = fft(sig3(2049:3072));% fft du signal à intercorreler
elseif any(sig3(1025:2048)~=0)
    X = fft(sig2);% fft du signal de reference
    Y = fft(sig3(1025:2048));% fft du signal à intercorreler
elseif any(sig3(1:1024)~=0)
    X = fft(sig2);% fft du signal de reference
    Y = fft(sig3(1:1024));% fft du signal à intercorreler
end

En espérant avoir bien compris ton code.
Pour le 2 aucune idée...

14/01/2008, 17h01
DOC doc doc

Tres bien merci beaucoup pour cette simplification...:yaisse2:

Tu as bien compris ce que je voulais faire.

Maintenant il ne me reste plus qu'à essayer de trouver le 2e point...:aie::aie:

Encore un petit pas en avant pour la rapidité d'exécution :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
%Calcul de X
X = fft(sig2);
%plus grand indice pour lequel sig3 est différent de zéro
index=max(find(sig3));
if index>3072
    %on a un sig3 différent de zéro dans le segment 3073:4096
    Y = fft(sig3(3073:4096));
elseif index>2048
    Y = fft(sig3(2049:3072));
elseif index>1024
    Y = fft(sig3(1025:2048));
else
    Y = fft(sig3(1:1024));
end
%Note : que se passe-t-il si sig3 est toujours égal à zéro ?

14/01/2008, 17h24
DOC doc doc

en fait sig3 ne sera jamais egal à '0'...:yaisse2: donc pas de prise de tete...:mrgreen::mrgreen:

Merci encore pour ce gain de rapidité d'execution. Tout marche bien...:king:
14/01/2008, 17h25
Jerome Briot

Citation:

Code:

index=max(find(sig3));

Mouais... je verrais plutôt :

Code:

[index,index]=max(sig3);

Ou avec les dernières versions :

Code:

index=find(sig3,1,'last');

A vérifier...
14/01/2008, 17h29
DOC doc doc

Oki merci beaucoup a vous deux... :)

Mais Dut aurais tu une idée pour le 2e petit probleme?? :roll::roll:

Merci d avance

Comme ceci ?

Code:

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
function balayageIntercorrelation
 
t = [0:4096]; % de 0 à 1000 points
sig = sin(t/(2*pi)); % Original signal, a sine wave
 
% Signal servant à l'intercorrelation (sur 1024 points obligatoirement)
sig2 = zeros(1, 1024);
sig2(1,1:256)= 0.000075*sig(1,1:256);
 
% Signal à intercorreler (taille sans zero padding : 256 points obligatoirement)
sig3 = zeros(1,4096);
sig3(1,1100:1356)=sig(1,1100:1356);
 
%Calcul de X
X = fft(sig2);
%plus grand indice pour lequel sig3 est différent de zéro
index=max(find(sig3));
% ou
% [index]=find(sig3,1,'last');
 
if index>3072
    idx = 3073:4096;
elseif index>2048
    idx = 2049:3072;
elseif index>1024
    idx = 1025:2048;
else
    idx = 1:1024;
end
 
Y = fft(sig3(idx));
 
% Calcul de l'intercorrelation
 
Z=zeros(size(sig));
Z(idx) = real(ifft(conj(X).*Y)); 
 
subplot(2,4,1), plot(sig), title('sig')
subplot(2,4,2), plot(sig2), title('sig2')
subplot(2,4,3), plot(sig3), title('sig3')
subplot(2,4,7), plot(Z), title('fft')

14/01/2008, 22h05
DOC doc doc

Ouiiiiiiiiiii....merci beaucoup à vous deux!!!

Maintenant je saurai reduire mon code et en plus je gagne de la rapidité d execution !!! :mrgreen::mrgreen:

Merci encore...Ca gere!! :king::king: