% Program 11_5
% Development of an AR model of an FIR filter
%

colordef black
b = remez(20, [0 0.5 0.6 1], [1 1 0 0]);
[h,w] = freqz(b,1,512);
[d,p0] = lpc(b,30);
[h1,w] = freqz(p0,d,512);
plot(w/pi,abs(h),w/pi,abs(h1),'--');
xlabel('\omega/\pi');ylabel('Magnitude');