MATLAB提供了专门用于求离散系统频响特性的函数freqz()，调用freqz()的格式有以下两种：

(1)[H,w]=freqz(B,A,N)

(2)[H,w]=freqz(B,A,N,’whole’)

    (1)中B和A分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了离散系统频响在 0~pi范围内N个频率等分点的值（其中N为正整数），w则包含了范围内N个频率等分点。调用默认的N时，其值是512。

(2)中调用格式将计算离散系统在0~pi范内的N个频率等分店的频率响应的值。

因此，可以先调用freqz()函数计算系统的频率响应，然后利用abs()和angle()函数及plot()函数，即可绘制出系统在 或 范围内的频响曲线。

绘制如下系统的频响曲线

H(z)=(z-0.5)/z

 MATLAB命令如下：

    B=[1 -0.5];

    A =[1 0];

    [H,w]=freqz(B,A,400,'whole');

    Hf=abs(H);

    Hx=angle(H);

    clf

    figure(1)

    plot(w,Hf)

title('离散系统幅频特性曲线')

    figure(2)

    plot(w,Hx)

    title('离散系统相频特性曲线')

本文来自CSDN博客，转载请标明出处：http://blog.csdn.net/silenceee/archive/2006/11/01/1361258.aspx