DSP課程設(shè)計(jì)-FIR高通濾波器設(shè)計(jì)

1、fir高通濾波器設(shè)計(jì)南京師范大學(xué)物科院摘要從實(shí)現(xiàn)方法方面考慮，將濾波器分為兩種，一種是iir濾波器，另一種是fir濾波器。firdf的最大優(yōu)點(diǎn)是可以實(shí)現(xiàn)線性相位濾波。而iirdf主要對幅頻特性進(jìn)行逼近，相頻特性會存在不同程度的非線性。我們知道，無失真?zhèn)鬏斉c濾波處理的條件是，在信號的有效頻譜范圍內(nèi)系統(tǒng)幅頻響應(yīng)應(yīng)為常數(shù)，相頻響應(yīng)為頻率的線性函數(shù)。另外，fir是全零點(diǎn)濾波器，硬件和軟件實(shí)現(xiàn)結(jié)構(gòu)簡單，不用考慮穩(wěn)定性問題。所以，firdf是一種很重要的濾波器，在數(shù)字信號處理領(lǐng)域得到廣泛應(yīng)用。firdf設(shè)計(jì)方法主要分為兩類：第一類是基于逼近理想濾波器特性的方法，包括窗函數(shù)法、頻率采樣法和等波紋最佳逼近法；

2、第二類是最優(yōu)設(shè)計(jì)法。其中窗函數(shù)計(jì)法的基本思想是用firdf逼近希望的濾波特性。本次設(shè)計(jì)主要采用窗函數(shù)設(shè)計(jì)法，對理想濾波器進(jìn)行逼近，從而實(shí)現(xiàn)高通濾波器的設(shè)計(jì)。在matlab軟件中，有一系列函數(shù)用于設(shè)計(jì)濾波器，應(yīng)用時十分方便。因此，在本次設(shè)計(jì)中，濾波器的設(shè)計(jì)主要采用matlab軟件，編寫適當(dāng)?shù)某绦?，得到濾波器的單位脈沖響應(yīng)。本設(shè)計(jì)對濾波器的硬件仿真主要使用ccs軟件，通過對濾波器的硬件仿真，可以較為真實(shí)的看出濾波器的濾波效果。關(guān)鍵字：高通、firdf、線性相位、hanning窗、matlab、ccs正文1. 設(shè)計(jì)目標(biāo)產(chǎn)生一個多頻信號，設(shè)計(jì)一個高通濾波器消除其中的低頻成分，通過ccs的graph v

3、iew波形和頻譜顯示，并和matlab計(jì)算結(jié)果比較2. 設(shè)計(jì)原理122.1 數(shù)字濾波器數(shù)字濾波器(digital filter)是由數(shù)字乘法器、加法器和延時單元組成的一種裝置。其功能是對輸入離散信號的數(shù)字代碼進(jìn)行運(yùn)算處理，以達(dá)到改變信號頻譜的目的。由于電子計(jì)算機(jī)技術(shù)和大規(guī)模集成電路的發(fā)展，數(shù)字濾波器已可用計(jì)算機(jī)軟件實(shí)現(xiàn)，也可用大規(guī)模集成數(shù)字硬件實(shí)時實(shí)現(xiàn)。數(shù)字濾波器廣泛用于數(shù)字信號處理中，如電視、vcd、音響等。按照濾波電路的工作頻帶為其命名：設(shè)截止頻率為fp，頻率低于fp的信號可以通過，高于fp的信號被衰減的電路稱為低通濾波器，頻率高于fp的信號可以通過，低于fp的信號被衰減的電路稱為高通濾波

4、器；而帶通嗎，就是頻率介于低頻段截止頻率和高頻段截止頻率的信號可以通過的電路。2.2 高通濾波器高通濾波器是容許高頻信號通過、但減弱（或減少）頻率低于截止頻率信號通過的濾波器。對于不同濾波器而言,每個頻率的信號的減弱程度不同。它有時被稱為低頻剪切濾波器；在音頻應(yīng)用中也使用低音消除濾波器或者噪聲濾波器。高通濾波器與低通濾波器特性恰恰相反。這樣的濾波器能夠把高頻率的聲音引導(dǎo)至專用高音喇叭(tweeter),并阻止可能干擾或者損害喇叭的低音信號。使用線圈而不是電容的低通濾波器也可以同時把低頻信號引導(dǎo)至低音喇叭(woofer)。高通和低通濾波器也用于數(shù)字圖像處理中在頻域中進(jìn)行變換。2.3 高通濾波器的

5、分析2.3.1 高通濾波器的時域分析在時域,信號經(jīng)過系統(tǒng)的響應(yīng)y (n)體現(xiàn)為激勵x(n)跟系統(tǒng)單位抽樣響應(yīng)h(n)的卷積和y(n)=(n)h(n)=n1m=0h(m)x(n-m)223 。對于長度為n 的fir系統(tǒng), h(n)可以看成一個長度為n 點(diǎn)的固定窗口,而x(n)則看成一個隊(duì)列以齊步走的方式穿過h(n)窗口,每走一步,位于窗口中的x(n)部分的點(diǎn)跟h(n)的對應(yīng)點(diǎn)的值相乘(即加權(quán))再求和,所得結(jié)果構(gòu)成此時系統(tǒng)的響應(yīng)值y(n), x(n)隊(duì)列每走一步就得到一個響應(yīng)值y(n),即y(n)是h(n)對位于其窗口中的x(n)的加權(quán)求和。高通濾波要求h(n)窗口具有波形銳化作用,即利用h(n)

6、窗口加權(quán)和使得變化快的(即高頻)正弦分量保留(理想高通)或衰減幅度小(實(shí)際高通) ,而變化緩慢(即低頻)的正弦分量正負(fù)抵消(理想高通)或衰減幅度大(實(shí)際高通) 。設(shè)其中n 必須取奇數(shù),其中: fs = 11khz, n為整數(shù),即x ( n)由100hz的x1(n)和1khz的x2(n)兩種頻率的信號組成。高通濾波的目的就是要盡可能地去掉x(n)中的低頻分量x1(n) ,同時盡可能地保留x(n)中的高頻分量x2(n)。2.3.2 高通濾波器的頻域分析在頻域,信號經(jīng)過系統(tǒng)的響應(yīng)y(n)的頻譜y(ej)體現(xiàn)為激勵x(n)的頻譜x (ej)跟系統(tǒng)單位抽樣響應(yīng)h(n)的頻譜h(ej)(即系統(tǒng)的頻譜)的乘

7、積y(ej)=|h(ej)|x(ej)|ejh+x,即響應(yīng)的頻譜y(ej)的幅值由系統(tǒng)頻譜h(ej)的幅值對激勵頻譜x(ej)的幅值相乘(加權(quán))得到,響應(yīng)的頻譜y(ej)的幅角由系統(tǒng)的頻譜h(ej)的幅角跟激勵頻譜x(ej)的幅角相加(移相)得到122 。高通濾波要求系統(tǒng)幅度函數(shù)|h(jf)|對需要保留的高頻信號頻譜加權(quán)權(quán)重較大(理想時為1) ,對需要濾除的低頻信號頻譜加權(quán)權(quán)重較小(理想時為0) 。其中l(wèi)為x(n)的長度(l=100), 0kl- 1,n為h(n)的有值長度,m、k均為整數(shù),跟數(shù)字頻率k相對應(yīng)的模擬頻率為f=fsk/l (hz)。h(n)在不同n值時的頻譜如圖2所示(橫軸單位為

8、khz),當(dāng)n=11時,在f=m(khz)即f=1khz、2khz、3khz等處為1,而在f=(2m-0.5)(khz)即f=1.5khz、3.5khz等處幅度最大。而在f=(2m+0.5)(khz)即f=2.5khz、4.5khz等處幅度最小,如圖2(b);當(dāng)n=5時, 在f=mfs/5處為1,即f=2.2khz、3.3khz等處為1,在f=(4m-1)fs/10即f=3.3khz處幅度最大。而在f=(4m+1)fs/10即f=5.5khz處幅度最小,如圖2(c);當(dāng)n=21時,在f=mfs/21即f=524hz、1047hz等處為1, 在f=(4m-1)fs/42處幅度最大,而在f=(4m

9、+1)fs/42處幅度最小,如圖1圖1用不同寬度的h(n)對x(n)的濾波在頻域上表現(xiàn)如圖2.2.1所示,圖2.2.1(a)為x(n)的頻譜|x(jf)| ,從圖2中可以看出x(n)中含有100hz和1khz兩種頻率的信號,圖1(b)(d)為不同長度的h(n)對同一x(n)的濾波情況,這跟在時域中分析的結(jié)論是一致的。2.4 fir濾波器從實(shí)現(xiàn)方法方面考慮，將濾波器分為兩種，一種是iir濾波器，另一種是fir濾波器。iir濾波器保留了模擬濾波器較好的幅度特性，設(shè)計(jì)簡單有效。但這些特性是以犧牲相位特性為代價而獲得的，然而現(xiàn)在許多數(shù)據(jù)傳輸，圖像處理系統(tǒng)都越來越多的要求系統(tǒng)具有線性相位特性。firdf

10、的最大優(yōu)點(diǎn)是可以實(shí)現(xiàn)線性相位濾波。而iirdf主要對幅頻特性進(jìn)行逼近，相頻特性會存在不同程度的非線性。我們知道，無失真?zhèn)鬏斉c濾波處理的條件是，在信號的有效頻譜范圍內(nèi)系統(tǒng)幅頻響應(yīng)應(yīng)為常數(shù)，相頻響應(yīng)為頻率的線性函數(shù)。另外，fir是全零點(diǎn)濾波器，硬件和軟件實(shí)現(xiàn)結(jié)構(gòu)簡單，不用考慮穩(wěn)定性問題。所以，firdf是一種很重要的濾波器，在數(shù)字信號處理領(lǐng)域得到廣泛應(yīng)用。firdf設(shè)計(jì)方法主要分為兩類：第一類是基于逼近理想濾波器特性的方法，包括窗函數(shù)法、頻率采樣法和等波紋最佳逼近法；第二類是最優(yōu)設(shè)計(jì)法。其中窗函數(shù)計(jì)法的基本思想是用firdf逼近希望的濾波特性。本次設(shè)計(jì)主要采用窗函數(shù)設(shè)計(jì)法，對理想濾波器進(jìn)行逼近，從

11、而實(shí)現(xiàn)高通濾波器的設(shè)計(jì)。2.5 各種窗函數(shù)參數(shù)各種窗函數(shù)的基本參數(shù)窗函數(shù)類型旁瓣峰值(db)過渡帶寬度阻帶最小衰減 (db)近似值精確值矩形窗-134/n1.8/n-21三角窗-258/n6.1/n-25漢寧窗-318/n6.2/n-44哈明窗-418/n6.6/n-74布萊克曼窗-5712/n11/n-80hanning（漢寧）窗的窗函數(shù)為：圖2 漢寧窗3. 設(shè)計(jì)流程1231233.1 應(yīng)用matlab設(shè)計(jì)濾波器：為了滿足高通濾波的效果，設(shè)置線性相位高通fir的參數(shù)為：通帶截止頻率，阻帶截止頻率，通帶最大衰減，阻帶最小衰減。由以上參數(shù)，通過查表可得，hanning窗和haming窗都滿足條件

12、，我們選擇hanning窗。過渡帶寬，hanning窗的精確過渡帶寬為=6.2/n，所以要求，解之得n=25。輸入程序：從而得到滿足條件的高通firdf：h(n)=-0.0004,-0.0006,0.0028,0.0071,-0.0000,-0.0185,-0.0210,0.0165,0.0624, 0.0355,-0.1061,-0.2898,0.6249,-0.2898,-0.1061,0.0355,0.0624,0.0165,-0.0210,0.0185,-0.0000,0.0071,0.0028,-0.0006,-0.0004其波形如圖3所示，頻譜如圖4，其幅頻特性和相頻特性如圖5。圖

13、3 圖4圖53.2 用ccs對高通濾波器進(jìn)行仿真設(shè)置輸入信號fin=sin(2n/30)/2+cos(2n*1.6),使用matlab運(yùn)算得到的高通fir濾波器進(jìn)行濾波，濾去低頻分量。選擇“debug”菜單的“run”項(xiàng)，或按f12 鍵運(yùn)行程序。觀察“fin”、“fout”窗口中時域圖形；觀察濾波效果。鼠標(biāo)右鍵單擊“input”和“output”窗口，選擇“properties”項(xiàng)，設(shè)置“display type”為“fft magitude”，再單擊“ok”按鈕結(jié)束設(shè)置。觀察“input”、“output”窗口中頻域圖形；理解濾波效果。 4. 仿真波形12344.1 matlab對高通fir

14、df的設(shè)計(jì)4.2 ccs對高通firdf的應(yīng)用仿真5. 設(shè)計(jì)結(jié)果由matlab和ccs的仿真結(jié)果可以看出，此高通firdf可以實(shí)現(xiàn)高通濾波的功能，在輸入為fin=sin(2n/30)/2+cos(2n*1.6)時，輸出波形中已不含w= 1/30的低頻分量，濾波效果良好，輸出頻譜中低頻部分基本沒有信號。6. 心得通過這次課程設(shè)計(jì)，使我掌握了課程設(shè)計(jì)的基本思路和方法，掌握了課程設(shè)計(jì)說明書的基本撰寫方法，了解了基本的設(shè)計(jì)思想和設(shè)計(jì)方法。同時此次課程設(shè)計(jì)使我課程的基礎(chǔ)知識和基本理論有了更深的理解和掌握，鍛煉了我們綜合運(yùn)用所學(xué)知識的能力，并在理論分析設(shè)計(jì)、計(jì)算制圖運(yùn)用標(biāo)準(zhǔn)和規(guī)范查閱設(shè)計(jì)手冊與資料以及計(jì)算

15、機(jī)應(yīng)用能了等方面得到了初步的訓(xùn)練和提高，培養(yǎng)了我們嚴(yán)謹(jǐn)求實(shí)的科學(xué)態(tài)度。這次課程設(shè)計(jì)使我進(jìn)一步熟悉了計(jì)算工具軟件-matlab.并進(jìn)一步掌握了matlab的使用方法.對matlab語言的發(fā)展和特點(diǎn)有了更深的了解，熟悉其工作環(huán)境,在其應(yīng)用方面也有了更深的了解，了解了matlab程序書寫的一般思路及一般步驟，學(xué)會了用matlab解決復(fù)雜信號處理等問題的方法。在課程設(shè)計(jì)的這段時間里，我認(rèn)為收獲還是很多的，不但進(jìn)一步掌握了數(shù)字信號處理的基礎(chǔ)知識及一門專業(yè)仿真軟件的基本操作，還提高了自己的設(shè)計(jì)能力及動手能力，同時對于模擬濾波器來了個系統(tǒng)的總結(jié)。更多的是讓我看清了自己，明白了凡事需要耐心，實(shí)踐是檢驗(yàn)學(xué)習(xí)的唯

16、一標(biāo)準(zhǔn)。理論知識的不足在這次課設(shè)中表現(xiàn)的很明顯。這將有助于我今后的學(xué)習(xí)，端正自己的學(xué)習(xí)態(tài)度，從而更加努力的學(xué)習(xí)。參考文獻(xiàn)1劉泉 闕大順主編. 數(shù)字信號處理原理與實(shí)現(xiàn). 北京：電子工業(yè)出版社，20052程佩青. 數(shù)字信號處理教程（第二版）. 北京：清華大學(xué)出版社，20013 高西全 數(shù)字信號處理-原理、實(shí)現(xiàn)及應(yīng)用. 北京：電子工業(yè)出版社，2010附錄matlab程序-生成高通firdfwp=pi/2;ws=pi/4;db=wp-ws;n0=ceil(6.2*pi/db);n=n0+mod(n0+1,2);m=1:n;wc=(wp+ws)/2/pi;hn=fir1(n-1,wc,high,hann

17、ing(n);fw=abs(fft(hn);figure(1)subplot(1,2,1)plot(m,hn);grid;subplot(1,2,2)plot(m,fw);grid;figure(2)freqz(hn) xlabel(omegapi)ylabel(|h(ejomega|db)subplot(2,1,2)xlabel(omegapi)ccs程序#include dsp281x_device.h / dsp281x headerfile include file#include dsp281x_examples.h / dsp281x examples include file#i

18、nclude f2812a.h#includemath.h#define firnumber 25#define signal1f 1000#define signal2f 4500#define samplef 10000#define pi 3.1415926float inputwave();float fir();float fhnfirnumber= -0.0004,-0.0006,0.0028,0.0071,-0.0000,-0.0185,-0.0210,0.0165,0.0624, 0.0355,-0.1061,-0.2898,0.6249,-0.2898,-0.1061,0.0355,0.0624,0.0165, -0.0210,0.0185,-0.0000,0.0071,0.0028,-0.0006,-0.0004 ;float fxnfirnumber= 0.0 ;float finput,foutput;float fsignal1,fsignal2;float fstepsignal