이고 을 만족하는 다음을 풀어보자.

ex4_15.cem

/*   Chapter 4: Example 4.15*/
 
 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]; /* X(z) transform coeff */
 ayplus = convo(a,axplus) /* Denominator of Yplus(z) */
 
 byplus = convo(b,bxplus)+convo(xic,axplus)  /* Numerator of Yplus(z) */
 [R,p,C] = resid_z(byplus,ayplus) Mp = abs(p); Ap = angle(p)/pi /* Polar form */
 
n = [0:7]; x = cos(pi*n/3); y = filtic(b,a,x,xic)/* CEMTool Verifigation*/
A=real(2*R(2)); B=imag(2*R(2)); C=real(2*R(4)); D=imag(2*R(4));
y=A*cos(pi*n/3)+B*sin(pi*n/3)+((0.95).^n).*(C*cos(pi*n/3)+D*sin(pi*n/3))