ms9r32.m

《数字信号处理实验(MATLAB版)》这

b=[4,-1.6,-1.6,4];a=[1,0.4,0.35,-0.4];
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)');