%《數字信號處理教程——MATLAB釋義與實現》第三章例5.2.9程序hc529
% z反變換求差分方程全解（復極點）
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
b=[1,1,1]/3; a=[1,-0.95,0.9025];     % 分子分母系數向量
Y=[-2,-3]; X=[1,1];                 % 初始條件
xic=filtic(b,a,Y,X)                 % 求初始條件引起的xic
bxplus=[1,-0.5]; axplus=[1,-1,1];   % X(z)的變換系數
ayplus=conv(a,axplus)           % Yplus(z)的分母
byplus=conv(b,bxplus)+conv(xic,axplus)     % Yplus(z)的分子
[R,p,C]=residuez(byplus,ayplus)     % 求極點留數
Mp=abs(p), Ap=angle(p)/pi           % 極坐標形式