q724a.m

来自「dsp相关的matlab程序」· M 代码 · 共 16 行

%《MATLAB在电子信息课程中的应用》第七章例7.24程序q724a
% 用各种窗函数设计FIR数字滤波器
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

clear;close all
f=[0,0.6,0.6,1]; m=[0,0,1,1];   % 设定预期幅频响应
b=fir2(30,f,m); n=0:30;         % 设计FIR数字滤波器系数
subplot(3,2,1),stem(n,b,'.')
xlabel('n'); ylabel('h(n)');
axis([0,30,-0.4,0.5]),line([0,30],[0,0])
[h,w]=freqz(b,1,256);
subplot(3,2,2),plot(w/pi,20*log10(abs(h)));grid
axis([0,1,-80,0]),
xlabel('w/pi'); ylabel('幅度(dB)');