example2_1.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 12 行

b=[0,0,0,0.25,-0.5,0.0625];
a=[1,-1,0.75,-0.25,0.0625];
[delta,n]=impseq(0,0,7);
x=filter(b,a,delta);
subplot(2,1,1);
stem(n,x);
disp(x);
xa=(n-2).*(0.5).^(n-2).*cos(pi*(n-2)/3).*stepseq(2,0,7);
subplot(2,1,2);
stem(n,xa);
disp(xa);