example1_11.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 21 行

b=[1];
a=[1,-1,0.9];
x=impseq(0,-20,100);
n=[-20:100];
h=filter(b,a,x);
subplot(2,1,1);
stem(n,h);
axis([-20,100,-1.1,1.1]);
title('pulse responese');
xlabel('n');ylabel('h(n)');
x1=stepseq(0,-20,100);
s=filter(b,a,x1);
subplot(2,1,2);
stem(n,s);
axis([-20,100,-0.5,2.5]);
title('step responese');
xlabel('n');ylabel('s(n)');
sumh=sum(abs(h));
disp(sumh);
z=roots(a)
magz=abs(z)