

[y,fs]=wavread('G.wav');

transf_fourier=fft(y(:,1))/length(y(:,1));
eixo_freq=linspace(0,fs,length(transf_fourier));


figure(1)
plot(eixo_freq,2*abs(transf_fourier))
axis([0 2000 0 max(2*abs(transf_fourier))])


figure(2)
spectrogram(y(:,1),2000,1900,[],fs);