%4、離散系統的頻率響應分析和零、極點分布
k=512;
num=[0.0528 0.0797 0.1295 0.1295 0.0797 0.0528];
den=[1 -1.8107 2.4947 -1.8801 0.9537 -0.2336];
w=0:pi/k:pi;
h=freqz(num,den,w);
subplot(2,1,1);plot(w/pi,abs(h));
title('幅度譜')
[z,p,k]=tf2zp(num,den);
disp('零點');disp(z);
disp('極點');disp(p);
subplot(2,1,2);zplane(num,den)
title('零極點圖');