39t.m

来自「包含50个数字信号处理方面的matlab程序」· M 代码 · 共 9 行

h=[-4,1,-1,-2,5,6,6,5,-2,-1,1,-4];
[Hr,w,b,L]=Hr_type2(h);
w=2*pi*n/500;
plot(w/pi,Hr,'b');xlabel('w/pi');hold on;%振幅响应
axis([0,2,-70,40]);n=0:1:500;
w=2*pi*n/500;
x=freqz(h,1,w);
X=abs(x);
plot(w/pi,x,'r');xlabel('w/pi');ylabel('dB');%频域响应