ms23_1.m

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

h=[3,-1,-5,4,6,4,-5,-1,3];
M=length(h);n=0:M-1;
[A,w,type,tao]=amplres(h);type
subplot(2,1,1),stem(n,h,'k');
ylabel('h(n)');xlabel('n');
subplot(2,1,2),plot(w/pi,A,'k');
ylabel('A');xlabel('\pi');
figure(2);
rz=roots(h)
zplane(h,1);
figure(3);
plot(n,angle(h)/pi,'k');