%《MATLAB在電子信息課程中的應用》(第三版)第七章例7.25程序q725
% 用各種窗函數設計FIR數字濾波器q724
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

clear;close all
N=21; wc=pi/4;						% 理想低通濾波器參數
n=0:N-1; r=(N-1)/2;
hdn=sin(wc*(n-r))/pi./(n-r); 	% 計算理想低通單位脈沖響應hd(n)
if rem(N,2)~=0   hdn(r+1)=wc/pi; end		% N為奇數時，處理n=r點的0/0型
wn1=boxcar(N); 		% 矩形窗
hn1=hdn.*wn1'; 		% 加窗
% 以上兩條語句可代以fir函數：hn1=fir1(N-1,wc/pi,boxcar(N));
wn2=hamming(N);		% hamming 窗
hn2=hdn.*wn2'; 		% 加窗
% 以上兩條語句可代以fir函數：hn2=fir1(N-1,wc/pi,hamming(N));k=3;
k=1;    %繪圖函數p751f的位置參數
hnwplot(hn1,k);title('矩形窗設計的h(n)')
%hamming 窗
wn2=hamming(N);
hn2=hdn.*wn2';
k=3;
hnwplot(hn2,k);title('hamming設計的h(n)')