FIR高通濾波器設計南京師范大學物科院摘要從實現方法方面考慮，將濾波器分為兩種，一種是IIR濾波器，另一種是FIR濾波器。FIRDF的最大優點是可以實現線性相位濾波。而IIRDF主要對幅頻特性進行逼近，相頻特性會存在不同程度的非線性。我們知道，無失真傳輸與濾波處理的條件是，在信號的有效頻譜范圍內系統幅頻響應應為常數，相頻響應為頻率的線性函數。另外，FIR是全零點濾波器，硬件和軟件實現結構簡單，不用考慮穩定性問題。所以，FIRDF是一種很重要的濾波器，在數字信號處理領域得到廣泛應用。FIRDF設計方法主要分為兩類第一類是基于逼近理想濾波器特性的方法，包括窗函數法、頻率采樣法和等波紋最佳逼近法；第二類是最優設計法。其中窗函數計法的基本思想是用FIRDF逼近希望的濾波特性。本次設計主要采用窗函數設計法，對理想濾波器進行逼近，從而實現高通濾波器的設計。在MATLAB軟件中，有一系列函數用于設計濾波器，應用時十分方便。因此，在本次設計中，濾波器的設計主要采用MATLAB軟件，編寫適當的程序，得到濾波器的單位脈沖響應。本設計對濾波器的硬件仿真主要使用CCS軟件，通過對濾波器的硬件仿真，可以較為真實的看出濾波器的濾波效果。關鍵字高通、FIRDF、線性相位、HANNING窗、MATLAB、CCS正文1設計目標產生一個多頻信號，設計一個高通濾波器消除其中的低頻成分，通過CCS的GRAPHVIEW波形和頻譜顯示，并和MATLAB計算結果比較2設計原理21數字濾波器數字濾波器DIGITALFILTER是由數字乘法器、加法器和延時單元組成的一種裝置。其功能是對輸入離散信號的數字代碼進行運算處理，以達到改變信號頻譜的目的。由于電子計算機技術和大規模集成電路的發展，數字濾波器已可用計算機軟件實現，也可用大規模集成數字硬件實時實現。數字濾波器廣泛用于數字信號處理中，如電視、VCD、音響等。按照濾波電路的工作頻帶為其命名設截止頻率為FP，頻率低于FP的信號可以通過，高于FP的信號被衰減的電路稱為低通濾波器，頻率高于FP的信號可以通過，低于FP的信號被衰減的電路稱為高通濾波器；而帶通嗎，就是頻率介于低頻段截止頻率和高頻段截止頻率的信號可以通過的電路。22高通濾波器高通濾波器是容許高頻信號通過、但減弱（或減少）頻率低于截止頻率信號通過的濾波器。對于不同濾波器而言,每個頻率的信號的減弱程度不同。它有時被稱為低頻剪切濾波器；在音頻應用中也使用低音消除濾波器或者噪聲濾波器。高通濾波器與低通濾波器特性恰恰相反。這樣的濾波器能夠把高頻率的聲音引導至專用高音喇叭TWEETER,并阻止可能干擾或者損害喇叭的低音信號。使用線圈而不是電容的低通濾波器也可以同時把低頻信號引導至低音喇叭WOOFER。高通和低通濾波器也用于數字圖像處理中在頻域中進行變換。23高通濾波器的分析231高通濾波器的時域分析在時域,信號經過系統的響應YN體現為激勵XN跟系統單位抽樣響應HN的卷積和YNNHNN1M0HMXNM223。對于長度為N的FIR系統,HN可以看成一個長度為N點的固定窗口,而XN則看成一個隊列以齊步走的方式穿過HN窗口,每走一步,位于窗口中的XN部分的點跟HN的對應點的值相乘即加權再求和,所得結果構成此時系統的響應值YN,XN隊列每走一步就得到一個響應值YN,即YN是HN對位于其窗口中的XN的加權求和。高通濾波要求HN窗口具有波形銳化作用,即利用HN窗口加權和使得變化快的即高頻正弦分量保留理想高通或衰減幅度小實際高通,而變化緩慢即低頻的正弦分量正負抵消理想高通或衰減幅度大實際高通。設其中N必須取奇數,/102SIN/102SIN821SSFNFXNX其中FS11KHZ,N為整數,即XN由100HZ的X1N和1KHZ的X2N兩種頻率的信號組成。高通濾波的目的就是要盡可能地去掉XN中的低頻分量X1N,同時盡可能地保留XN中的高頻分量X2N。232高通濾波器的頻域分析在頻域,信號經過系統的響應YN的頻譜YEJ體現為激勵XN的頻譜XEJ跟系統單位抽樣響應HN的頻譜HEJ即系統的頻譜的乘積YEJ|HEJ|XEJ|EJHX,即響應的頻譜YEJ的幅值由系統頻譜HEJ的幅值對激勵頻譜XEJ的幅值相乘加權得到,響應的頻譜YEJ的幅角由系統的頻譜HEJ的幅角跟激勵頻譜XEJ的幅角相加移相得到122。高通濾波要求系統幅度函數|HJF|對需要保留的高頻信號頻譜加權權重較大理想時為1,對需要濾除的低頻信號頻譜加權權重較小理想時為0。其中L為XN的長度L100,0KL1,N為HN的有值長度,M、K均為整數,跟數字頻率K相對應的模擬頻率為FFSK/LHZ。HN在不同N值時的頻譜如圖2所示橫軸單位為KHZ,當N11時,在FMKHZ即F1KHZ、2KHZ、3KHZ等處為1,而在F2M05KHZ即F15KHZ、35KHZ等處幅度最大。而在F2M05KHZ即F25KHZ、45KHZ等處幅度最小,如圖2B當N5時,在FMFS/5處為1,即F22KHZ、33KHZ等處為1,在F4M1FS/10即F33KHZ處幅度最大。而在F4M1FS/10即F55KHZ處幅度最小,如圖2C當N21時,在FMFS/21即F524HZ、1047HZ等處為1,在F4M1FS/42處幅度最大,而在F4M1FS/42處幅度最小,如圖1圖1用不同寬度的HN對XN的濾波在頻域上表現如圖221所示,圖221A為XN的頻譜|XJF|,從圖2中可以看出XN中含有100HZ和1KHZ兩種頻率的信號,圖1BD為不同長度的HN對同一XN的濾波情況,這跟在時域中分析的結論是一致的。24FIR濾波器從實現方法方面考慮，將濾波器分為兩種，一種是IIR濾波器，另一種是FIR濾波器。IIR濾波器保留了模擬濾波器較好的幅度特性，設計簡單有效。但這些特性是以犧牲相位特性為代價而獲得的，然而現在許多數據傳輸，圖像處理系統都越來越多的要求系統具有線性相位特性。FIRDF的最大優點是可以實現線性相位濾波。而IIRDF主要對幅頻特性進行逼近，相頻特性會存在不同程度的非線性。我們知道，無失真傳輸與濾波處理的條件是，在信號的有效頻譜范圍內系統幅頻響應應為常數，相頻響應為頻率的線性函數。另外，FIR是全零點濾波器，硬件和軟件實現結構簡單，不用考慮穩定性問題。所以，FIRDF是一種很重要的濾波器，在數字信號處理領域得到廣泛應用。FIRDF設計方法主要分為兩類第一類是基于逼近理想濾波器特性的方法，包括窗函數法、頻率采樣法和等波紋最佳逼近法；第二類是最優設計法。其中窗函數計法的基本思想是用FIRDF逼近希望的濾波特性。本次設計主要采用窗函數設計法，對理想濾波器進行逼近，從而實現高通濾波器的設計。25各種窗函數參數各種窗函數的基本參數過渡帶寬度B窗函數類型旁瓣峰值DBN近似值精確值阻帶最小衰減DBS矩形窗134/N18/N21三角窗258/N61/N25漢寧窗318/N62/N44哈明窗418/N66/N74布萊克曼窗5712/N11/N80HANNING（漢寧）窗的窗函數為12COS5012COS150NNRNNRNNNNNH0510152025303500102030405060708091圖2漢寧窗3設計流程31應用MATLAB設計濾波器為了滿足高通濾波的效果，設置線性相位高通FIR的參數為通帶截止頻率，阻帶截止頻率，通帶最大衰減，阻帶最RADP2/RADS4/DBP1小衰減。BS40由以上參數，通過查表可得，HANNING窗和HAMING窗都滿足條件，我們選擇HANNING窗。過渡帶寬，HANNING窗的精確過渡帶寬為4/PS62/N，所以要求，解之得N25。BN62B輸入程序從而得到滿足條件的高通FIRDFHN00004,00006,00028,00071,00000,00185,00210,00165,00624,00355,01061,02898,06249,02898,01061,00355,00624,00165,00210,00185,00000,00071,00028,00006,00004其波形如圖3所示，頻譜如圖4，其幅頻特性和相頻特性如圖5。0102030030201001020304050607010203000204060811214圖3圖400102030405060708091150010005000500PIPHASEDEGREES0010203040506070809110050050PI|HEJ|DB圖532用CCS對高通濾波器進行仿真設置輸入信號FINSIN2N/30/2COS2N16,使用MATLAB運算得到的高通FIR濾波器進行濾波，濾去低頻分量。選擇“DEBUG”菜單的“RUN”項，或按F12鍵運行程序。觀察“FIN”、“FOUT”窗口中時域圖形；觀察濾波效果。鼠標右鍵單擊“INPUT”和“OUTPUT”窗口，選擇“PROPERTIES”項，設置“DISPLAYTYPE”為“FFTMAGITUDE”，再單擊“OK”按鈕結束設置。觀察“INPUT”、“OUTPUT”窗口中頻域圖形；理解濾波效果。4仿真波形41MATLAB對高通FIRDF的設計05101520250302010010203040506070510152025108060402002040608142CCS對高通FIRDF的應用仿真5設計結果由MATLAB和CCS的仿真結果可以看出，此高通FIRDF可以實現高通濾波的功能，在輸入為FINSIN2N/30/2COS2N16時，輸出波形中已不含W1/30的低頻分量，濾波效果良好，輸出頻譜中低頻部分基本沒有信號。6心得通過這次課程設計，使我掌握了課程設計的基本思路和方法，掌握了課程設計說明書的基本撰寫方法，了解了基本的設計思想和設計方法。同時此次課程設計使我課程的基礎知識和基本理論有了更深的理解和掌握，鍛煉了我們綜合運用所學知識的能力，并在理論分析設計、計算制圖運用標準和規范查閱設計手冊與資料以及計算機應用能了等方面得到了初步的訓練和提高，培養了我們嚴謹求實的科學態度。這次課程設計使我進一步熟悉了計算工具軟件MATLAB并進一步掌握了MATLAB的使用方法對MATLAB語言的發展和特點有了更深的了解，熟悉其工作環境,在其應用方面也有了更深的了解，了解了MATLAB程序書寫的一般思路及一般步驟，學會了用MATLAB解決復雜信號處理等問題的方法。在課程設計的這段時間里，我認為收獲還是很多的，不但進一步掌握了數字信號處理的基礎知識及一門專業仿真軟件的基本操作，還提高了自己的設計能力及動手能力，同時對于模擬濾波器來了個系統的總結。更多的是讓我看清了自己，明白了凡事需要耐心，實踐是檢驗學習的唯一標準。理論知識的不足在這次課設中表現的很明顯。這將有助于我今后的學習，端正自己的學習態度，從而更加努力的學習。參考文獻1劉泉闕大順主編數字信號處理原理與實現北京電子工業出版社，20052程佩青數字信號處理教程（第二版）北京清華大學出版社，20013高西全數字信號處理原理、實現及應用北京電子工業出版社，2010附錄MATLAB程序生成高通FIRDFWPPI/2WSPI/4DBWPWSN0CEIL62PI/DBNN0MODN01,2M1NWCWPWS/2/PIHNFIR1N1,WC,HIGH,HANNINGNFWABSFFTHNFIGURE1SUBPLOT1,2,1PLOTM,HNGRIDSUBPLOT1,2,2PLOTM,FWGRIDFIGURE2FREQZHNXLABELOMEGAPIYLABEL|HEJOMEGA|DBSUBPLOT2,1,2XLABELOMEGAPICCS程序INCLUDE“DSP281X_DEVICEH“/DSP281XHEADERFILEINCLUDEFILEINCLUDE“DSP281X_EXAMPLESH“/DSP281XEXAMPLESINCLUDEFILEINCLUDE“F2812AH“INCLUDE“MATHH“DEFINEFIRNUMBER25DEFINESIGNAL1F1000DEFINESIGNAL2F4500DEFINESAMPLEF10000DEFINEPI31415926FLOATINPUTWAVEFLOATFIRFLOATFHNFIRNUMBER00004,00006,00028,00071,00000,00185,00210,00165,00624,00355,01061,02898,06249,02898,01061,00355,00624,00165,00210,00185,00000,00071,00028,00006,00004FLOATFXNFIRNUMBER00FLOATFINPUT,FOUTPUTFLOATFSIGNAL1,FSIGNAL2FLOATFSTEPSIGNAL1,FSTEPSIGNAL2FLOATF2PIINTIFLOATFIN256,FOUT256INTNIN,NOUTMAINVOIDNIN0NOUT0F2PI2PIFSIGNAL100FSIGNAL2PI01FSTEPSIGNAL12PI/30FSTEPSIGNAL22PI16WHILE1FINPUTINPUTWAVEFINNINFINPUTNINNIN25