求单位抽样响应.m

来自「数字信号处理工具箱」· M 代码 · 共 11 行

%-----------------------------------------------------------------
% exa020802_impz.m, for example 2.8.2
% to test impz.m and to obtain the impulse response.
%-------------------------------------------------------------
clear;

b=[.001836,.007344,.011016,.007374,.001836];
a=[1,-3.0544,3.8291,-2.2925,.55075];
[h,t]=impz(b,a,40);         % 求单位抽样响应
stem(t,h,'.');grid on;