%-----------------------------------------------------------------
% 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;