Discussion :

[yoann_deygas]

Bonjour a tous
un autre problème :
mon signal correspond à un vecteur de 6 000 000 points
Faut-il tous les garder ou faut-il raccourcir le signal ?

Si il faut le raccourcir , comment faire?

Faut-il tester la périodicité, et si c'est periodique on ne garde qu'une seule période ? Si c'est ça, comment calculer la périodicité de mon signal et garder la bonne longueur?

Ou faut-il faire autrement


merci d'avance

[Ivory69]

Salut!
un conseil, il vaut mieux que tu essayes de voir d'abord si ton signal est périodique pour réduire le nombre d'échantillons.
Tu peux essayer la transformée de Fourier et prendre les harmoniques les plus importantes et tester si ton signal est effectivement périodique à cette fréquence là.

Une méthode simple pour repérer la périodicité d'un signal est d'analyser le spectre de puissance de ton signal. Je m'explique:
On suppose que ton signal a été echantillonée à une fréquence fe et il est de longueur L. On note y ton signal dans l'espace des temps (y(t)) et Y ton signal dans l'espace des fréquences (Y(f)).

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
T = 1/fe;                       % Temps d'échantillonnage
L = 6000000;                  % Longueur du signal
t = (0:L-1)*T;                % Vecteur temps
NFFT = 2^nextpow2(L);   % Puissance de 2 qui suit la longueur de ton signal
Y = fft(y,NFFT)/L;
f = fe/2*linspace(0,1,NFFT/2);
% Affichage du spectre de puissance.
plot(f,2*abs(Y(1:NFFT/2))) 
title('Spectre de puissance de y(t)')
xlabel('Fréquence (Hz)')
ylabel('|Y(f)|')

Le pic maximal correspond à la fréquence prépondérante dans ton signal qui doit être la fréquence de périodicité de ton signal. Fais gaffe pcq ton signal peut avoir plusieurs fréquences de périodicité (en principe ce sont des multiples de la fréquence fondamentale).
Une fois que tu as cette fréquence tu peux tester ta périodicité en testant l'intercorrélation entre Y(0:1/fréquencepéridicité) et Y(1/fréquencepéridicité+1:2/fréquencepéridicité).
Voilà ça devrait t'aider.
@ +