
%% load data
[x,fs]=auload("bruant.wav");
%% or [x,fs]=wavread("bruant.wav");
N=length(x); 

%% define time and frequency axes
t=1:20:N;
Nf=512;
f=(0:Nf/2)/Nf;

%% compute spectrogram
tfr=tfrsp(x,t,Nf);

%% use a logarithmic scale and select positive frequencies
tfr=log10(tfr(1:Nf/2+1,:));

%% make plot
colormap(jet);
imagesc(t/fs,fs*f,flipud(tfr));

%% note: matlab users can use a GUI to adjust the layout
%% and other parameters of the display. this utility is
%% not yet available with GNU octave.