q725a.m

来自「<matlab在电子信息中的应用>一书中个章节的配套源代码。」· M 代码 · 共 16 行

%《MATLAB在电子信息课程中的应用》(第三版)第七章例7.25程序q725a
% 用各种窗函数设计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)');