example2_18.m

西交大刘树棠写的《数字信号处理》那本书的所有源代码

b=[1,1,1]/3;
a=[1,-0.95,0.9025];
Y=[-2,-3];
X=[1,1];
xic=filtic(b,a,Y,X)
bx=[1,-0.5];
ax=[1,-1,1];
ay=conv(a,ax);
by=conv(b,bx)+conv(xic,ax);
[R,P,C]=residuez(by,ay)
Mp=abs(P),Ap=angle(P)/pi
n=[0:7];x=cos(pi*n/3);
y=filter(b,a,x,xic)