ms23_52.m

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

h=[3,-1,-5,4,4,-5,-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)