c9ex1.m

为了更好的便于增加交流

n1 = [1 -1];                    % Define numerator factor
num = conv(n1,conv(n1,n1));     % Calculate numerator
d1 = [1 -1/2];                  % Define denominator factor
d2 = [1 -1/8];                  % Define second den. factor
den = conv(d1,d2);              % Calculate denominator
[r,p,k] = residuez(num,den);    % Calculate residues
r                               % Display residues
p                               % Display poles
k                               % Display k vector