ex4p15.m

来自「关于数字信号处理的一些matlab例程」· M 代码 · 共 22 行

b=[1,1,1]/3;a=[1,-0.95,0.9025];
Y=[-2,-3];X=[1,1];
xic=filtic(b,a,Y,X)

bxplus=[1,-0.5];axplus=[1,-1,1];
ayplus=conv(a,axplus)

byplus=conv(b,bxplus)+conv(xic,axplus)

[R,P,C]=residuez(byplus,ayplus)

Mp=abs(p),Ap=angle(p)/pi

n=[0:7];x=cos(pi*n/3)
y=filter(b,a,x,xic)


A=real(2*R(1));B=imag(2*R(1));C=real(2*R(4));D=imag(2*R(3));

y=A*cos(pi*n/3)+B*sin(pi*n/3)+((0.95).^n).*(C*cos(pi*n/3)+D*sin(pi*n/3))