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