의 디지털 공명기에서 CEMTool을 사용하여 극점의 위치 변화를 구해보자.

ex6_25.cem

x = [-0.9, 0.81];
[y,L,B] = QCoeff(x,3) 
%Unquantized parameters
r =0.9; theta = pi/3; a1 = -2*r*cos(theta); a2 = r*r;
p1 = r*exp(j*theta); p2 = p1';
%Quantized parameters : N = 3;
[ahat, L, B] = QCoeff([a1,a2],3); rhat = sqrt(ahat(2));
thetahat = acos(-ahat(1)/(2*rhat)); p1hat = rhat*exp(j*thetahat);
p2 = p1';
%changes in pole locations
Dp1 = abs(p1-p1hat)