matlab:＜三＞freqz函数

【转载:博客】

一 函数调用介绍

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

(1)[H,w]=freqz(B,A,N);  %N默认值为512

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

(3)[H,w]=freqz(B,A,[自定义区间]);

上式中B和A分别对应离散系统的系统函数H(z)的分子、分母多项式的系数向量：

       B=[b1,b2,...];    A=[a1,a2,...];

返回量H则包含了离散系统对应区间内N（N为正整数）个频率等分点的频率响应，w为N个频率等分点的值。

(1)式调用默认区间 0：pi

(2)式调用主值区间 -pi：pi

(3)式自定义区间，如[0:2*pi/n:2*pi]

最后利用abs()和angle()函数及plot()函数，即可绘制出系统在频率区间的频响曲线。

二 示例

绘制如下系统的频响曲线： H(z)=(1-0.5z^-1)

MATLAB命令如下：

B=[1 -0.5]; 

A =[1];

[H,w]=freqz(B,A,[0:2*pi/400:2*pi]);

Hf=abs(H);  %取幅度值实部

Hx=angle(H);  %取相位值对应相位角

clf

figure(1)

plot(w,20*log(Hf))  %幅值变换为分贝单位

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

figure(2)

plot(w,Hx)

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