Normalisation de la fonction spectrogramme

Bonjour à tous et à toutes.

Mon tuteur m'a demander de normaliser un spectrogramme provenant d'un signal temporel.

J'ai un peu scripté, j'obtiens un spectrogramme mais je n'arrive pas à le normaliser correctement.

Code:

---

1 2 3 4 5 6 7 8 9 10

x=signal; nlap=4000 nfft=13000 longueur=26000 figure spectrogram(x,longueur,nlap,nfft,Fs,'yaxis'); [B,f,t]=spectrogram(x,longueur,nlap,nfft,Fs,'yaxis'); L=numel(B(:,1)); M = median(B,2);

---

Le fait de le normaliser comme on peux normaliser une fft classique permet de retrouver une valeur RMS identique aux signaux temporelles et fréquentiel.

Si vous avez une idée me permettant d'avancer, je suis preneur !

Nb: ce n'est pas une normalisation vis à vis de l'amplitude max, mais sur le domaine fréquentiel ou temporelle ( vis à vis du nombre d'élément dans mon signal)



Merci,


Thomasgeo

Salut, j'ai finalement trouvé

Il suffit de normaliser par la fenêtre de hamming :) et par la fréquence d’échantillonnage.

Code:

---

1 2 3 4

[B,f,t]=spectrogram(signal,longueur,[],[],Fs,'yaxis'); % Calcul du spectrogramme U=mean(window(@hamming,longueur).^2); % Calcul d'un élément de normalisation vis à vis de la fenetre de Hamming B=2*abs(B).^2/(longueur*Fs*U);% Normalisation du spectrogramme

---

Résolut !