exa1_44.m

来自「matlab应用实例。如：FIR、IIR数字滤波器设计等。」· M 代码 · 共 16 行

% exa1-44_filter_impz_unit_pulse.m for example 1-44
% filter and impz function realize unit pulse response

clear all;
pulse=[1 zeros(1,63)]; %generate unit sample sequence
b=[0.3 0.2]; %sytem parameter
a=[1 -0.4 -0.7]; %sytem parameter
h1=filter(b,a,pulse);
h2=impz(b,a,64);
subplot(211)
stem(h1)
title('filter function')
subplot(212)
stem(h2)
title('impz function')