ms23_54.m

来自「《数字信号处理实验(MATLAB版)》这」· M 代码 · 共 10 行

h=[3,-1,-5,4,-4,5,1,-3];
%h=[3,-1,2,-3,-3,2,-1,3];
M=length(h);n=0:M-1;
rz=roots(h)
for i=1:7
   r1(i)=1/rz(i);
end
r'
zplane(h,1)