ex2p9.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 22 行

b=[1];
a=[1,-1,0.9];
x=impseq(0,-20,120);
n=[-20:120];
h=filter(b,a,x);
subplot(2,1,1);stem(n,h);
title('Impulse response');
xlabel('n');
ylabel('h(n)')

x=stepseq(0,-20,120);
s=filter(b,a,x);
subplot(2,1,2);stem(n,s);
title('step response');
xlabel('n');
ylabel('s(n)')

sum(abs(h))
z=roots(a);
magz=abs(z)