1.2t.m

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

b=[1,1,1];a=[1,0.5,-0.25];
subplot(3,1,1);zplane(b,a);title('因果系统零极图');
n=0:50;
x=3*cos(pi*n/3);
y=filter(b,a,x);
subplot(3,1,2);stem(n,x,'.');title('输入x的波形');
subplot(3,1,3);stem(n,y,'.');title('输出y的波形');