%% load and process data
[Y, FS] = audioread('alemi.mp3');

[S,F,T,P] = spectrogram(Y(:,1),4096,4000,4096,FS,'yaxis');

%% make a plot
maxF = 1e3;
myF = F(F<maxF);
figure
C  = log10(P(F<maxF,:));
imagesc(T,myF,C)
set(gca,'YDir','Normal')
caxis([min(min(C)) max(max(C))])

set(gca,'FontSize',16)

axis tight
xlabel 'Time (s)', ylabel 'Frequency (Hz)'
colormap gray
ylim([0 maxF])

print -dpng boing_alemi.png