% lpc_spectrum: plot FFT + superimposed LPC spectrum clear all [x,rate]=wavread('wheel.wav'); x=x(1501:2000); x=x.*hamming(500); p=12; a=lpc(x,p); nfft=512; freq=linspace(0,rate/2,nfft/2+1); fft_logmag = 20 * log10 ( abs( fft( x, nfft ) ) ); fft_logmag = fft_logmag(1:nfft/2+1); plot(freq,fft_logmag,'b'); lpc_logmag = -20 * log10 ( abs( fft( a, nfft ) ) ); lpc_logmag = lpc_logmag(1:nfft/2+1); hold on; plot(freq,lpc_logmag-3,'b'); hold off xlabel('Frequency (Hz)') ylabel('Magnitude (dB)')
