example2_28.m

来自「原名：数字信号处理基础及MATLAB实现,解压缩密码：lfjk2002@163.」· M 代码 · 共 15 行

%例2-27程序，example2_27
b=[1,0];a=[1,-0.8];
[H,w]=freqz(b,a,100);
magH=abs(H);phaH=angle(H);
subplot(1,1,1)
subplot(2,2,1);plot(w/pi,magH);grid
xlabel('频率（以pi为单位）');ylabel('幅度')
title('幅度响应')
subplot(2,2,2);plot(w/pi,phaH/pi);grid
xlabel('频率（以pi为单位）');ylabel('相位 单位：pi')
title('相位响应')
subplot(2,2,3);zplane(b,a)
title('零极点图')
subplot(2,2,4);impz(b,a,40)
title('单位抽样响应')