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