【正文】 寬度/(2π/N)過渡帶寬度Δω/(2π/N)阻帶最小衰減/dB矩形窗132221三角窗254425漢寧窗314444漢明窗414453布萊克曼窗576674凱澤窗（β=）575580最小阻帶衰減只由窗形狀決定，不受N的影響，而過渡帶的寬度則隨窗函數(shù)的增加而減小。利用Simulink工具進行數(shù)字濾波系統(tǒng)的仿真。一個簡單的FIR數(shù)字濾波系統(tǒng)如圖2所示。經(jīng)過所設計的FIR低通數(shù)字濾波器濾濾波。仿真結果分別如圖4所示。該設計方法簡單易行，對于了解數(shù)字濾波器的原理和設計方法有一定的參考價值。參數(shù)sflag是用來控制窗函數(shù)首尾的兩個元素值：取值為symmetric或periodic，默認值為symmetric。而調用函數(shù) w = hanning(n,‘periodic’)則回包括為零點的窗函數(shù)的首尾兩個元素。 %清除工作空間clc。 %濾波去階數(shù)a=(N1)/2。 %截止頻率n=0:N1。 %避免被零除，eps為MATLAB的系統(tǒng)精度hd=sin(wc*m)./(pi*m)。 [H10,w10]=freqz(wd1)。[H1,w1]=freqz(h1)。 [H20,w20]=freqz(wd2)。[H2,w2]=freqz(h2)。 [H30,w30]=freqz(wd3)。[H3,w3]=freqz(h3)。 [H40,w40]=freqz(wd4)。[H4,w4]=freqz(h4)。 [H50,w50]=freqz(wd5)。[H5,w5]=freqz(h5)。 [H60,w60]=freqz(wd6)h6=hd.*wd6。%凱澤窗%%%繪制窗函數(shù)的沖激響應、幅頻響應subplot(3,2,1)。.39。 %矩形窗沖激響應title(39。)。n39。ylabel(39。)。plot(w10/pi,20*log10(abs(H10)/max(abs(H10))))。 %幅頻響應title(39。)。歸一化頻率39。ylabel(39。)。plot(w10/pi,180/pi*unwrap(angle(H10)))。 %相頻響應title(39。)。歸一化頻率39。ylabel(39。39。subplot(3,2,2)。.39。title(39。)。n39。ylabel(39。)。plot(w20/pi,20*log10(abs(H20)/max(abs(H20))))。title(39。)。歸一化頻率39。ylabel(39。)。plot(w20/pi,180/pi*unwrap(angle(H20)))。title(39。)。歸一化頻率39。ylabel(39。39。subplot(3,2,1)。.39。title(39。)。n39。ylabel(39。)。plot(w30/pi,20*log10(abs(H30)/max(abs(H30))))。title(39。)。歸一化頻率39。ylabel(39。)。plot(w30/pi,180/pi*unwrap(angle(H30)))。title(39。)。歸一化頻率39。ylabel(39。39。subplot(3,2,2)。.39。title(39。)。n39。ylabel(39。)。plot(w40/pi,20*log10(abs(H40)/max(abs(H40))))。title(39。)。歸一化頻率39。ylabel(39。)。plot(w40/pi,180/pi*unwrap(angle(H40)))。title(39。)。歸一化頻率39。ylabel(39。39。subplot(3,2,1)。.39。title(39。)。n39。ylabel(39。)。plot(w50/pi,20*log10(abs(H50)/max(abs(H50))))。title(39。)。歸一化頻率39。ylabel(39。)。plot(w50/pi,180/pi*unwrap(angle(H50)))。title(39。)。歸一化頻率39。ylabel(39。39。subplot(3,2,2)。.39。title(39。)。n39。ylabel(39。)。plot(w60/pi,20*log10(abs(H60)/max(abs(H60))))。title(39。)。歸一化頻率39。ylabel(39。)。plot(w60/pi,180/pi*unwrap(angle(H60)))。title(39。)。歸一化頻率39。ylabel(39。39。%%%繪制各種窗函數(shù)的幅頻響應subplot(2,2,1)。.39。axis([0 30 0 1])。矩形窗波形39。xlabel(39。)。h(n)39。subplot(2,2,2)。.39。title(39。)。n39。ylabel(39。)。subplot(2,2,3)。.39。title(39。)。n39。ylabel(39。)。subplot(2,2,4)。.39。title(39。)。n39。ylabel(39。)。subplot(2,2,1)。.39。title(39。)。n39。ylabel(39。)。subplot(2,2,2)。.39。title(39。)。n39。ylabel(39。)。%對輸入信號進行采樣濾波NN=1000。%數(shù)據(jù)總數(shù)和采樣頻率nn=[0:NN1]。f1=50。x=sin(2*pi*f1*t)+cos(2*pi*f2*t)。t1=t(n1)。%與采樣點對應的輸入信號plot(t1,x1)。帶濾波信號39。xlabel(39。)。幅度39。%%%把輸入信號經(jīng)濾波器處理y1=fftfilt(h1,x,256)。 y3=fftfilt(h3,x,256)。 y5=fftfilt(h5,x,256)。%%%處理輸出信號n2=n1(N1)/2。 %輸出信號扣除了相位延遲y11=y1(n2)。 y13=y3(n2)。 y15=y5(n2)。%%%繪制窗函數(shù)設計濾波器的沖激響應、幅頻特性、相頻特性、輸出信號%%矩形窗設計濾波器subplot(2,2,1)。.39。 %濾波器沖激響應title(39。)。n39。ylabel(39。)。plot(w1,20*log10(abs(H1)))。 %設定坐標軸范圍，以方便比較title(39。)。頻率 w39。ylabel(39。)。 %濾波器相頻響應axis([0 pi 2000 0])。矩形窗濾波器相頻響應39。xlabel(39。)。相位/176。)。plot(t2,y11)。矩形窗濾波器輸出信號39。xlabel(39。)。幅度39。%三角窗設計濾波器subplot(2,2,1)。.39。title(39。)。n39。ylabel(39。)subplot(2,2,2)。axis([0 pi 150 0])。三角窗濾波器幅頻響應20lg|H(ejw)|曲線39。xlabel(39。)。幅頻響應/dB39。subplot(2,2,3),plot(w2,180/pi*unwrap(angle(H2)))。title(39。)。頻率 w39。ylabel(39。39。subplot(2,2,4)。title(39。)。時間/t39。ylabel(39。)。stem(n,h3,39。)。漢寧窗設計濾波器h(n)39。xlabel(39。)。h(n)39。plot(w3,20*log10(abs(H3)))。title(39。)。頻率39。ylabel(39。)。axis([0 pi 2000