ms4r321.m

《数字信号处理实验(MATLAB版)》这

N=32;
a=[1,0,0,0,0,0];
b=[1,0.5,-0.5,-1,-0.5,1];
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('系统单位阶跃响应');