ms4_11.m

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

a=[1, 0, 0.34319, 0, 0.60439, 0, 0.20407];
b=[1, 0, -3, 0, 3, 0, - 1];
n=30;
hn=impz(b,a,n);                      %求时域冲激响应
gn=dstep(b,a,n);					      %求时域阶跃响应
subplot(1,2,1),stem(hn,'k');             %显示冲激响应曲线
title('系统的单位冲激响应');
ylabel('h(n)');xlabel('n');		      %显示阶跃响应曲线
subplot(1,2,2),stem(gn,'k');
title('系统的单位阶跃响应');
ylabel('g(n)');xlabel('n');