MATLAB常用信号分析处理函数演示示例

 编写一个M文件，分别用filter和impz函数获取制定离散系统（b=[0 0.5 0.3 0.2],a=[1 0.2 0.4 -0.8]）的单位抽样响应；用freqz函数获取系统的频率响应；分别用root及zplane函数获取系统的零/极点图及增益。程序的源代码如下：

% SignalProcess 文件源码 
L = 128;                            %单位抽样序列的长度
Fs = 1000;                          %采样频率为1khz
b = [0.8 0.5 0.6];                  %系统函数的分子系数向量
a = [1 0.2 0.4 -0.8];               %系统函数的分母系数向量
delta = [1 zeros(1,L-1)];           %生成长度为L的单位抽样序列

FilterOut = filter(b,a,delta);      %filter函数获取单位抽样响应
ImpzOut = impz(b,a,L);              %impz函数获取单位抽样响应
[h,f] = freqz(b,a,L,Fs);            %freqz函数求频率响应
mag = 20*log(abs(h)/log(10));       %幅度转换成dB单位
ph = angle(h)*180/pi;               %相位值单位转换
zr = roots(b);                      %求系统的零点，并显示在命令窗口
pk = roots(a);                      %求系统的极点，并显示在命令窗口
g = b(1)/a(1);                      %求系统的增益，并显示在命令窗口

%绘图
figure(1);
subplot(221);stem(FilterOut);
subplot(222);stem(ImpzOut);
subplot(223);plot(f,mag);
subplot(224);plot(f,ph);
figure(2);
freqz(b,a);                         %用feqz函数绘制系统频率响应
figure(3);
zplane(b,a);                        %用zplane函数绘制系统零极点图

程序运行结果如下图，同时可以得到系统的零极点以及增益值。

由上图可知，filter()函数和impz()函数获取的系统单位抽样响应序列完全相同。freqz()函数两种用法的频率响应完全相同，其中函数自动绘制的幅频响应横坐标为归一化的数值。