%《數字信號處理教程——MATLAB釋義與實現》第七章例7.3.5程序hc735
% 用哈明窗設計FIR低通濾波器
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
wp = 0.2*pi; ws = 0.3*pi;  deltaw= ws - wp;		% 過渡帶寬Δω的計算
N0 = ceil(6.6*pi/ deltaw) + 1;	% 按哈明窗(表7.3.2)計算所需的濾波器階數N0
N=N0+mod(N0+1,2);		% 為了實現第一類偶對稱濾波器，應使其長度N為奇數
wdham = (hamming(N))';	% 求窗函數
wc = (ws+wp)/2	;		% 截止頻率取為兩邊界頻率的平均值
hd = ideallp(wc,N);		% 求理想脈沖響應	
h = hd .* wdham;		% 設計的脈沖響應(即系數)為理想脈沖響應與窗函數乘積
[db,mag,pha,grd,w] = myfreqz(h,[1]);	% 對設計結果進行檢驗
dw = 2*pi/1000;
Rp = -(min(db(1:wp/dw+1)))        	% 檢驗通帶波動
As = -round(max(db(ws/dw+1:501))) 	% 檢驗最小阻帶衰減
%繪圖
n=0:N-1;
subplot(2,2,1); stem(n,hd,'.'); title('理想脈沖響應')
axis([0 N-1 -0.1 0.3]); ylabel('hd(n)');text(N+1,-0.1,'n')
subplot(2,2,2); stem(n,wdham,'.');title('哈明窗')
axis([0 N-1 0 1.1]); ylabel('w(n)');text(N+1,0,'n')
subplot(2,2,3); stem(n,h,'.');title('實際脈沖響應')
axis([0 N-1 -0.1 0.3]); xlabel('n'); ylabel('h(n)')
subplot(2,2,4); plot(w/pi,db);title('幅度響應（單位： dB）');grid
axis([0 1 -100 10]); xlabel('頻率（單位：pi）'); ylabel('分貝')
set(gca,'XTickMode','manual','XTick',[0,0.2,0.3,1])
set(gca,'YTickMode','manual','YTick',[-50,0])
set(gca,'YTickLabelMode','manual','YTickLabels',['50';' 0'])
set(gcf,'color','w');                % 置圖形背景色為白