[x, fs] = auload('res_short.wav');
[X, freqs] = freqz(x(1:512), 1, 512, fs);
plot(freqs, 20*log10(abs(X)), '.;input;'); hold on;
A1 = lpc(x, 8);
[H1, freqs] = freqz(1, A1, 512, fs);
plot(freqs, 20*log10(abs(H1)), ';LPC: 8;'); 
A2 = lpc(x, 16);
[H2, freqs] = freqz(1, A2, 512, fs);
plot(freqs, 20*log10(abs(H2)), ';LPC: 16;');  
A3 = lpc(x, 32);
[H3, freqs] = freqz(1, A3, 512, fs);
plot(freqs, 20*log10(abs(H3)), ';LPC: 32;');  hold off;
xlabel('frequency [Hz]');
ylabel('[dB]');
axis([0, 20000, -20, 25]);
replot;
gset term postscript
gset output "lpc_analysis.eps"
replot;
gset term x11