ms4_12.m

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

x01=0;y01=0;N=32;
a=[1,0,1/3,0];
b=[1/6,1/2,1/2,1/6];
xi=filtic(b,a,0);
n=0:N-1;
x1=[n==0];
hn=filter(b,a,x1,xi);
x2=[n>=0];
gn=filter(b,a,x2,xi);
subplot(1,2,1),stem(n,hn,'k');
axis([0,N,1.1*min(hn),1.1*max(hn)]);
title('系统单位冲激响应');
subplot(1,2,2),stem(n,gn,'k');
axis([0,N,1.1*min(gn),1.8*max(gn)]);
title('系统单位阶跃响应');