在MATLAB中，用filter函数直接实现IIR的直接形式。

例5.1 用直接型实现系统函数为http://s15/mw690/001DLjb7gy6PtFSVDMG9e&690的IIR数字滤波器，求单位脉冲响应和单位阶跃信号的输出。

 图1 直接型单位脉冲响应和输出信号

http://s14/mw690/001DLjb7gy6PtEmwQEJcd&690

解：程序清单如下：

% 2015-1-26 余成波-IIR-滤波器
% p157 例5.1
% filename:  ycb_p157_5_1

%-----------------------------------

b=[1,-3,11,27,18];a=[16,12,2,-4,-1];
N=25;
h=impz(b,a,N);  %直接型单位脉冲响应
x=[ones(1,5),zeros(1,N-5)]; %单位阶跃信号
y=filter(b,a,x);    %直接型输出信号
subplot(1,2,1);stem(h);title('直接型h(n)');
subplot(1,2,2);stem(y);title('直接型y(n)');

参考资料：余成波等，数字信号处理及MATLAB实现（第二版）[M].北京：清华大学出版社，2008.1