ms4_2.m

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

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