ms9r33.m

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

r=[2,-1,1]';
p=[1,0.5,-0.5]';
k=0;
[b,a]=residuez(r,p,k);
rz=roots(b)
rp=roots(a)
subplot(1,2,1),zplane(b,a);
title('系统的零极点分布图');
subplot(1,2,2),impz(b,a,20);
title('系统的冲激响应');
xlabel('n');ylabel('h(n)');