exa144.m

来自「自己整理的信号与系统中的MATLAB应用： 离散系统模型及其MATLAB实现 」· M 代码 · 共 14 行

%exa144_filter_impz_unit_pulse.m for example 1-44
%filter and impz function realize unit pulse reponse 利用filter函数与impz函数实现单位冲激响应
clear all;
pulse=[1 zeros(1,63)];%generate unit sample sequence 产生一个64点的单位冲激信号
b=[0.3 0.2];%system parameter
a=[1 -0.4 -0.7];%system parameter
h1=filter(b,a,pulse);
h2=impz(b,a,64);
subplot(2,1,1)
stem(h1)
title('filter function')
subplot(2,1,2)
stem(h2)
title('impz function')