La Fonction spectrogramm

Bonjour,

j'ai un problème sur Matlab, je n'ai pas bien compris la fonction spectrogramm, on a:

[S,F,T] = SPECTROGRAM(...) returns a vector of frequencies F and a
vector of times T at which the spectrogram is computed. F has length
equal to the number of rows of S. T has length k (defined above) and
its value corresponds to the center of each segment

Donc F ce sont les fréquences d'un signal et T les portions du temps
Si je fait:

Code:

---

1 2 3

t=0:0.1:10; x=cos(2*pi*0.2*t); [S,F,T] = SPECTROGRAM(x)

---

F devrait contenir une seule valeur, d'une seule fréquence, mais il me donne beaucoup de valeurs dans une colonne

et T contient 8 valeurs dans une ligne

je n'ai pas tres bien compris cette fonction

il ne donne pas la fréquence à un temps donné. il effectue une suite de sft sur une certaine période avec un certain nombre de points.
le nombre de points de la fft donne la résolution fréquentielle : la longuer du vecteur F
les valeurs qu'ils te donnent poue toutes ces fréquences à un temps t donné te permettront de savoir si la fréquence est présente ou non. je ne sais pas si je suis clair mais si tu fais une simple fft de ton signal, il ne te donnera pas qu'une valeur, il te donnera une valeur par fréquence.

Il existe un moyen d'extraire la fréquence que j'ai mis?

Il y a déjà une longue discussion sur ce sujet :
http://www.developpez.net/forums/d87...s-fichier-wav/

Vous devriez en parler entre vous :roll:

avec ton exemple j'ai des résultats étranges.
peut-être qu'il faut régler des paramètres de fft et de résoilution temporelle du fait de ta fréquence d'échantillonnage basse.
cad : pour certin tn on retrouve bien 0.2, pour d'autres on a 0 Hz
rq : 1<=tn<=length(T)
sinon pour avoir UNE fréquence prédominante à UN instant T(tn) je ferai

Code:

---

1 2 3 4 5 6 7 8 9 10 11

t=0:0.1:10; x=cos(2*pi*0.2*t); [s,f,t,p] = spectrogram(x); freq=f/pi; tn=1; %si tu n'utilises pas p écrire p=s.*conj(s); [m,i]=max(p(1:size(p)/2,tn)); freq=f(i); %pour tracer ler spectre au temps tn plot(f(1:length(f)/2),p(1:size(p,1)/2,tn))

---

oui j'ai lu l'autre post, mais j'ai pas trouvé de solution

avec ce dernier code, si j'augment la frep, si je met 0.8 au lieu de 0.2, ca me donne pas le résultat exacte, ca me donne 0.4 dans freq:?

edit: autant pour moi ta fréquence d'échantillonage est bien de 10Hz ce qui devrait modifier la définition de f (multiplication par un facteur 10 à priori). je ne comprends pas trop car il me semblait avoir vérifié, aussi avec une simle fft (et ça n'aurait pas du me donner le bon résultat puisque je me trompais de fréquence d'échantillonage)...
pour être sûr relie la doc de la fonction et teste différents cas pour voir

tu peux consulter des cours de traitement de signal sinon, par exemple ici partie traitement du signal