ms4r311.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 15 行

N=32;
a=[1,-1,1];
b=[1,-0.5,0];
n=0:N-1;
hn=impz(b,a,n);                       %求时域单位冲激响应
gn=dstep(b,a,n);					        %求时域单位阶跃响应
subplot(1,2,1),stem(n,hn);
axis([0,N,1.1*min(hn),1.1*max(hn)]);
title('系统单位冲激响应');
subplot(1,2,2),stem(n,gn);
axis([0,N,1.1*min(gn),1.1*max(gn)]);
title('系统单位阶跃响应');