Question

Matlab中给定滤波器传函，如何绘制幅频响应 ？

问题是这样的，如果已知一个低通IIR数字滤波器的传函为：

(z+0.2)*(z^2-0.5z+0.8)
Hlp（z）=0.35* ----------------------
(z-0.3)*(z^2-1.5z+0.7)

如何绘制该低通IIR数字滤波器的幅频响应特性？各位大哥帮帮吧，matlab语言！

如何得到它的幅频响应特性如下：

就是说，那几条语句怎么写的，嗨！

Answer 1

一般来说，分析离散系统频率响应特性就要绘制频率响应曲线，而这是相当麻烦的。

虽然可以通过几何矢量法来定性画出频率响应特性曲线，但一般来说这也是很麻烦的。

值得庆幸的是,MATLAB提供了专门用于求解离散系统频率响应的函数freqz() ，其调用格式如下：

                 

其中，B和A分别是表示待分析的离散系统的系统函数的分子，分母多项式的向量，如公式所示：

N为正整数，返回向量H则包含了离散系统频率响应函数的N个频率等分点的值。

默认情况下N=512。

由于调用freqz()函数只能求出离散系统频率响应的数值，不能直接绘制曲线图，因此，可以先用freqz()函数求出系统频率响应的值，然后再利用MATLAB的abs()和angle()函数以及plot()命令，即可绘制出系统在或范围内的幅频特性和相频特性曲线。

Answer 2

[H,W]=freqz(b,a,N,'whole',Fs)
b=[0.35 0.245 0.245 0.056];
a=[1 1.8 0.7 1.15 0.21];
[H,W]=freqz(b,a,256,'whole',1);
hr=abs(H)

plot（Hr）