dan@0: fs = 441;;
dan@0: N = 2^10-1;
dan@0: 
dan@0: ply = [1000, -65, -1.2];
dan@0: 
dan@0: t = [-(N-1)/2:(N-1)/2];
dan@0: 
dan@0: p_val = polyval(ply,t);
dan@0: f = log(p_val);
dan@0: f_fit = polyfit(t,f,20);
dan@0: real(f_fit(end - 3:end))
dan@0: 
dan@0: ply(end-1)/ply(end)
dan@0: log(ply(end))
dan@0: 
dan@0: figure(1);clf;
dan@0: subplot(311);
dan@0: plot(p_val);
dan@0: subplot(312);hold on
dan@0: plot(real(f)); 
dan@0: plot(real(polyval(f_fit,t)),'r--');
dan@0: subplot(313);
dan@0: plot(imag(f)/pi*-2+1);