example2_18.m

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

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)