%《MATLAB及在電子信息課程中的應(yīng)用》第七章例7.24程序q724
% 用各種窗函數(shù)設(shè)計FIR數(shù)字濾波器q724
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear;close all
N=21; wc=pi/4;						    % 理想低通濾波器參數(shù)
n=0:N-1; r=(N-1)/2;
hdn=sin(wc*(n-r))/pi./(n-r); 	        % 計算理想低通單位脈沖響應(yīng)hd(n)
if rem(N,2)~=0   hdn(r+1)=wc/pi; end	% N為奇數(shù)時，處理n=r點的0/0型
wn1=boxcar(N); 		                    % 矩形窗
hn1=hdn.*wn1'; 		                    % 加窗
% 以上兩條語句可代以fir函數(shù)：hn1=fir1(N-1,wc/pi,boxcar(N));
wn2=hamming(N);		                    % hamming 窗
hn2=hdn.*wn2'; 		                    % 加窗
% 以上兩條語句可代以fir函數(shù)：hn2=fir1(N-1,wc/pi,hamming(N));k=3;
k=1;                                    % 繪圖函數(shù)p751f的位置參數(shù)
hnwplot(hn1,k);title('矩形窗設(shè)計的h(n)')
set(gcf,'color','w')                    % 設(shè)置圖形背景色為白色
%hamming 窗
wn2=hamming(N);
hn2=hdn.*wn2';
k=3;
hnwplot(hn2,k);title('hamming設(shè)計的h(n)')