Matlab - Espectro de Frecuencias

[x,fs]=wavread('pcg.wav');%lee señal de audio, x=señal, fs=frecuencia de muestreo

espectro=fft(x);%transformada rapida de fourier de x

espectro=espectro(1:end/2);

freq=linspace(0,fs/2,length(espectro));

espreal=real(espectro);

espimg=imag(espectro);

magnitud=abs(espectro);

fase=angle(espectro);

t=0:1/fs:(length(x)-1)/fs;

figure(1)

subplot(3,1,1);

plot(t,x);

xlabel('tiempo de muestreo')

subplot(3,1,2);

plot(freq,espreal);

title('espectro real');

subplot(3,1,3);

plot(freq,espimg);

title('espectro imaginario');

figure(2)

subplot(2,1,1);

plot(freq,magnitud);

title('espectro de magnitud');

subplot(2,1,2);

plot(freq,fase);

title('espectro de fase');