ms4_8.m

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

a=[1,0,1/3,0];
b=[1/6,1/2,1/2,1/6];
n=32;
hn=dimpulse(b,a,n);                       %求时域单位冲激响应
gn=dstep(b,a,n);					        %求时域单位阶跃响应
subplot(1,2,1),stem(hn,'k');          %显示冲激响应曲线
title('系统的单位冲激响应');
ylabel('h(n)');xlabel('n');		     %显示阶跃响应曲线
axis([0,n,1.1*min(hn),1.1*max(hn)]);
subplot(1,2,2),stem(gn,'k');
title('系统的单位阶跃响应');
ylabel('g(n)');xlabel('n');
axis([0,n,1.1*min(gn),1.1*max(gn)]);