example2_17.m

来自「数字信号处理的比较全面的程序代码」· M 代码 · 共 8 行

b=1;
a=[1,-3/2,1/2];
Y=[4,10];
xic=filtic(b,a,Y)
n=[0:7];
x=(1/4).^n;
y1=filter(b,a,x,xic)
y2=(1/3)*(1/4).^n+(1/2).^n+(2/3)*ones(1,8)