q724a.m

来自「matlab在电子信息中的应用和实例中的源代码」· M 代码 · 共 18 行

%《MATLAB及在电子信息课程中的应用》第七章例7.24程序q724a
% 用各种窗函数设计FIR数字滤波器
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月初版，2003年7月第二版
%
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,'.')
set(gcf,'color','w')                    % 设置图形背景色为白色
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)');