%《數字信號處理教程——MATLAB釋義與實現》第三章例5.2.3程序hc523
% 極點留數法求z反變換（實極點）及長除法
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
b=1;a=poly([0.9,0.9,-0.7]);         % 分子分母系數向量
[r,p,k]=residuez(b,a)               % 求極點留數
n=0:7;                              % 給出時間序列
x=(r(1)+r(2))*p(1).^n+r(2)*n.*p(2).^n+r(3)*p(3).^n  % 求脈沖響應
x1=deconv([b,zeros(1,10)],a)        % 長除法求脈沖響應