Discussion :

[Thomasgeo]

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 :

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

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

[Thomasgeo]

Salut, j'ai finalement trouvé

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

Code :

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

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 !