試驗用窗函數(shù)法設計FIR濾波器分析解析

1、實驗六用窗函數(shù)法設計FIR濾波器一、實驗目的(1)掌握用窗函數(shù)法設計FIR數(shù)字濾波器的原理和方法。(2)熟悉線性相位FIR數(shù)字濾波器特性。(3)了解各種窗函數(shù)對濾波特性的影響。二、實驗原理濾波器的理想頻率響應函數(shù)為H(ej3),則其對應的單位脈沖響應為:兀hd(n)=Hd(ej)e1nd-ji窗函數(shù)設計法的基本原理是用有限長單位脈沖響應序列h(n)逼hd(n)。由于hd(n)往往是無限長序列，且是非因果的，所以用窗函數(shù)。w(n)將hd(n)截斷，并進行加權處理：h(n)=hd(n)w(n)h(n)就作為實際設計的FIR數(shù)字濾波器的單位脈沖響應序列，其頻率響應函數(shù)H(ej")為:NJH

2、(ejw)=2h(n)en=0如果要求線性相位特性，則h(n)還必須滿足：h(n)=h(N-1-n)可根據(jù)具體情況選擇h(n)的長度及對稱性。用窗函數(shù)法設計的濾波器性能取決于窗函數(shù)w(n)的類型及窗口長度N的取值。設計過程中，要根據(jù)對阻帶最小衰減和過渡帶寬度的要求選擇合適的窗函數(shù)類型和窗口長度No三、實驗步驟1,寫出理想低通濾波器的傳輸函數(shù)和單位脈沖響應。2.寫出用四種窗函數(shù)設計的濾波器的單位脈沖響應。3,用窗函數(shù)法設計一個線性相位FIR低通濾波器，用理想低通濾波器作為逼近濾波器，截止頻率Wc=U/4rad,選擇窗函數(shù)的長度N=15,33兩種情況。要求在兩種窗口長度下，分別求出h(n),打印出

3、相應的幅頻特性和相頻特性曲線，觀察3dB帶寬和阻帶衰減；4用其它窗函數(shù)(漢寧窗(升余弦窗)、哈明窗(改進的升余弦窗)、布萊克曼窗)設計該濾波器，要求同1;比較四種窗函數(shù)對濾波器特性的影響。四、實驗用MATLAIB數(shù)可以調(diào)用MATLAB:具箱函數(shù)fir1實現(xiàn)本實驗所要求的線性相位FIR-DF的設計，調(diào)用一維快速傅立葉變換函數(shù)fft來計算濾波器的頻率響應函數(shù)。firl是用窗函數(shù)法設計線性相位FIRDF的上其箱函數(shù)，調(diào)用格式如卜:hn=fir1(N,wc,'ftype',window)firl實現(xiàn)線性相位FIR濾波器的標準窗函數(shù)法設計。hn=fir1(N,wc)可得到6dB截止頻率為

4、wc的N階(單位脈沖響應h(n)長度為N+1)FIR低通濾波器，默認(缺省參數(shù)windows)選用hammiing窗。其單位脈沖響應h(n)滿足線性相位條件：h(n)=h(N-1-n)其中wc為對冗歸一化的數(shù)字頻率，0&wc01。當wc=wc1,wc21時，彳J到的姑.帶通淤波器hn=fir1(N,wc,'ftype')U設計而通和帶阻濾波器。當ftype=high時，設計高通FIR濾波器：當ftype=stop時，設計帶阻FIR濾波器：應當注意，在設計高通和帶阻濾波器時，階數(shù)N只能取偶數(shù)(h(n)長度N+1為奇數(shù))。不過，當用戶將N設置為奇數(shù)時，fir1會自動對N加1

5、。hn=fir1(N,wc,window)可以指定窗函數(shù)向量window。如果缺省window參數(shù)，則fir1默認為hamming窗?？捎玫钠渌昂瘮?shù)有Boxcar,Hanning,Bartlett,Blackman,Kaiser和Chebwin窗。例如：hn=fir1(N,wc,bartlett(N+1)使用Bartlett窗設計；hn=fir1(N,wc,chebwin(N+1,R)使用Chebyshev窗設計。ftype和window參數(shù)(含hn=fir1(N,wc,'ftype',window)通過選擇wo義同上)，可以設計各種加窗濾波器。五、實驗內(nèi)容<1>

6、理想低通濾波器傳輸函數(shù):H(ej)e7mls|<®cQ,%<|®H冗a=(N-1)/2單位脈沖響應:隨尸六二Hd(ej)ej»=/=%小二(n一二)<2>四種窗函數(shù)的設計濾波器的單位脈沖響應矩形窗：R(n)=RN(n)h(n)=hd(n)R4n)2-n一(2)漢寧窗：WHn(n)=0.51-cos()R4n)N-12-nh(n)=hd(n)-0.51-cos(-n-)RN(n)N-1(3)哈明窗：wHn(n)=0.54-0.46cos(-n)RN(n)N-12-n一h(n)=hd(n)0.54-0.46cos(-n)R«n)N-1

7、(4)布萊克曼窗：Bi(n)=0.42-0.5cos(-l)+0.08cos(&1)RN(n)N-1N-12二n4二nh(n)=hd(n)0.42-0.5cos()+0.08cos(-)RN(n)N-1N-1<3>MATALAB序設計實驗程序1:N=1533,即N變化時，對濾波器特性的影響:MATLA程序：clearall;N1=15;N2=33;Wc=pi/4;wc=Wc/pi;%頻率歸一化b1=fir1(N1-1,wc,boxcar(N1);b2=fir1(N2-1,wc,boxcar(N2);H1,W=freqz(b1,1);H1_db=20*log10(abs(H1

8、);magH1=abs(H1);phaH1=angle(H1);H2,W=freqz(b2,1);H2_db=20*log10(abs(H2);magH2=abs(H2);時，實際低通濾波器的h(n),)副頻衰減特性,)相頻特性,)副頻特性,)時，實際低通濾波器的h(n),)副頻衰減特性,)相頻特性,)副頻特性,)phaH2=angle(H2);subplot(2,4,1);stem(b1);title('N=15subplot(2,4,2);plot(W/pi,H1_db);title('subplot(2,4,3);plot(W/pi,phaH1);title('s

9、ubplot(2,4,4);plot(W/pi,magH1);title('subplot(2,4,5);stem(b2);title('N=33subplot(2,4,6);plot(W/pi,H2_db);title('subplot(2,4,7);plot(W/pi,phaH2);title('subplot(2,4,8);plot(W/pi,magH2);title('實驗波形截圖:N=15時、實際低通濾波器的MR副頻亮牖特性實驗程序2:N=1§即N固定不變時，不同類型的窗函數(shù)對濾波器特性的影響:MATLA程序:clearall;N=1

10、5;Wc=pi/4;wc=Wc/pi;%頻率歸一化b1=fir1(N-1,wc,boxcar(N);b2=fir1(N-1,wc,hanning(N);b3=fir1(N-1,wc,bartlett(N);b4=fir1(N-1,wc,blackman(N);H1,W=freqz(b1,1);H1_db=20*log10(abs(H1);magH1=abs(H1);phaH1=angle(H1);H2,W=freqz(b2,1);H2_db=20*log10(abs(H2);magH2=abs(H2);phaH2=angle(H2);H3,W=freqz(b3,1);H3_db=20*log1

11、0(abs(H3);magH3=abs(H3);phaH3=angle(H3);H4,W=freqz(b4,1);H4_db=20*log10(abs(H4);magH4=abs(H4);phaH4=angle(H4);subplot(4,4,1);stem(b1);title('矩形窗:subplot(4,4,2);plot(W/pi,H1_db);title('subplot(4,4,3);plot(W/pi,phaH1);title('subplot(4,4,4);plot(W/pi,magH1);title('實際低通濾波器的副頻衰減特性')相頻

12、特性')副頻特性')h(n)1)subplot(4,4,5);stem(b2);title('漢寧窗:subplot(4,4,6);plot(W/pi,H2_db);title('subplot(4,4,7);plot(W/pi,phaH2);title('subplot(4,4,8);plot(W/pi,magH2);title('實際低通濾波器的副頻衰減特性')相頻特性')副頻特性')h(n)1)subplot(4,4,9);stem(b1);title('三角窗:subplot(4,4,10);plot(W/

13、pi,H1_db);title('subplot(4,4,11);plot(W/pi,phaH1);title('subplot(4,4,12);plot(W/pi,magH1);title('實際低通濾波器的副頻衰減特性')相頻特性')副頻特性')h(n)1)subplot(4,4,13);stem(b1);title('subplot(4,4,14);plot(W/pi,H1_db);title('subplot(4,4,15);plot(W/pi,phaH1);title('subplot(4,4,16);plot(

14、W/pi,magH1);title('布拉克曼窗：實際低通濾波器的h(n)')副頻衰減特性')相頻特性')副頻特性')實驗波形截圖:副期裝展后性相先特性-fl50510060.6助相賽屬特性0.6副演蓑展特性0510漢寧茵；實際低通謔波器的h間祖陽窗：實際低通涯,波器的hw5/心口uJq15三角茵：實際低通遮波器的h0515布拉克同密：其際低通淞波器朋六、實驗結(jié)論1 .N的大小決定了窗譜的主瓣寬度，N越大，窗譜的主瓣寬度越大。2 .最小阻帶衰減只有窗行決定，不受N的影響，過渡帶寬度與N和窗形都有關，N越大，過渡帶寬越小。3 .由實驗可知濾波特性：布拉克曼窗漢寧窗三角窗矩形窗。七、思考題(1)如果給定通帶截止頻率和阻帶截止頻率以及阻帶最小衰減，如何用窗函數(shù)法設計線性相位低通濾波器？寫出設計步驟。答：根據(jù)小缶=%一8P計算出所要設計的濾波器的過渡帶查表找到滿足的阻帶最小衰減As所需的窗函數(shù)利用已選的窗函數(shù)的過渡帶寬除以Acc得到低通濾波器單位脈沖響應的長度N根據(jù)hn)=hd(n0(n)求出h(n)N再根據(jù)H(ejt0)=£hn)e"j求出響應函數(shù)H(ejtt)n6(2)如果要求用窗函數(shù)法設計帶通濾波器，且給定上、下邊帶截止頻率為i和2,試求理想帶通的單位脈沖響應