hc529.m

来自「这是数字信号处理——matlab释义与实现一书的源代码」· M 代码 · 共 13 行

%《数字信号处理教程——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           % 极坐标形式