3t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 21 行

b=[2];
a=[1 -0.8 -0.5];
subplot(4,1,1);
zplane(b,a);
title('系统的零极图');
[H,w]=freqz(b,a,100,'whole');
magH=abs(H);
phaH=angle(H);
subplot(4,1,2);
plot(w/pi,magH);
title('系统的幅频响应');
subplot(4,1,3);
plot(w/pi,phaH/pi);
title('系统的相频响应');
n=0:100;
x=impseq(0,0,100);
h=filter(b,a,x);
subplot(4,1,4);
stem(n,h,'.');
title('系统的冲激响应');