
zplane([1 1.655 1.655 1],[1 -1.57 1.264 -0.4]);
 
pause
figure
[h,w]=freqz([1 1.655 1.655 1],[1 -1.57 1.264 -0.4],40);

plot(w, abs(h));
pause
plot(w, (-20*log(abs(h))));
title('Magnitude in db');
pause
figure
plot(w,angle(h));
title('Phase in radians');

pause
figure
grpdelay([1 1.655 1.655 1],[1 -1.57 1.264 -0.4]);
