4t.m

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

b=[1 1];
a=[1 -0.9 0.81];
[H,w]=freqz(b,a,400,'whole');
magH=abs(H);
phaH=angle(H);
subplot(4,1,1);
plot(w/pi,magH);
title('系统的幅频响应');
subplot(4,1,2);
plot(w/pi,phaH/pi);
title('系统的相频响应');
n=0:200;
x=sin(pi*n/3)+5*cos(pi*n);
y=filter(b,a,x);
subplot(4,1,3);
plot(n,x);
title('输入信号X');
subplot(4,1,4);
plot(n,y);
title('输出信号Y');
grid;