Sur l'utilisation de spectrogramme de la fonction dans matlab

Considérons l'exemple suivant:

Fs = 40;% sampling frequency
x = 0:(1/Fs):4;% time domain
y = [sin(2 * pi * 5 * x(x <= 2)), sin(2 * pi * 10 * x(x > 2))];% signal
N = length(x);                     % Length of signal
NFFT = 2^nextpow2(N); % Next power of 2 from length of y
Y = fft(y,NFFT)/N;
f = Fs/2*linspace(0,1,NFFT/2+1);
% Generate the plot, title and labels.
fh = figure(1);
% set(fh,'color','white','visible','off');
subplot(311);
plot(x,y,'k');
xlabel('Time (s)','FontName','Times New Roman','fontsize',10);
ylabel('Amplitude','FontName','Times New Roman','fontsize',10);
set(gca,'FontName','Times New Roman','fontsize',10);
% # Frequency domain plots
subplot(312);
plot(f,2*abs(Y(1:NFFT/2+1))) 
xlabel('Frequency (cycles/second)','FontName','Times New Roman','fontsize',10);
ylabel('Amplitude','FontName','Times New Roman','fontsize',10);
set(gca,'FontName','Times New Roman','fontsize',10);
subplot(313);
window = x(1:10:end);
[S,F,T] = spectrogram(y,window);
pcolor(T,F,abs(S));shading interp;
xlabel('Time (s)');
ylabel('Frequency (cycles/second)');

La TFCT (de courte durée à transformée de Fourier) ici n'est pas la preuve de ce que j'attendais. Je me serais attendu à l'axe des y être le même que le xaxis dans l'intrigue secondaire(312) et la xaxis être le même que le xaxis dans l'intrigue secondaire(311).

InformationsquelleAutor KatyB | 2013-10-03
  1. 3

    Le problème avec l'axe des x est parce que vous êtes de ne pas spécifier la fréquence d'échantillonnage lors de l'appel de spectrogram. Essayez

    [S,F,T] = spectrogram(y,window,[],[],Fs);

    Aussi, je pense que vous êtes une mauvaise interprétation de la window argument. Il se réfère à la fenêtre de la forme (la façon dont vous le définir, peut-être que vous pensez à elle définit les positions de la fenêtre?). Je voudrais donc utiliser quelque chose comme

    window = ones(1,10);

    (ou d'essayer de la fenêtre de forme autre que rectangulaire).

    Comme pour l'axe y, en 313, l'axe y est (instantanée) de la fréquence, alors qu'en 312, il est spectral amplitude. Ils ne sont donc pas comparables. En principe, vous pouvez comparer les couleur (axe z) de 313 avec l'axe y de 312, car les deux sont spectrale d'amplitude; cependant, leur normalisation est probablement différente.

    InformationsquelleAutor Luis Mendo