ms9r31.m

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

z=[1,-3]';
p=[2,-4]';
k=5;
[b,a]=zp2tf(z,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)');