%《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)');