【正文】 1 湖北理工學院 本科生畢業(yè)論文（設計） 題 目 ： 基于 MATLAB的 FIR 帶通數(shù)字濾波器設 計 姓 名： 鄭潔 學 院： 機電工程學院 專 業(yè)： 機械電子工程 班 級： 2020 學 號： 202040140133 指導教師： 孫川 完成時間 ： 2020年 5 月 25 畢業(yè)論文（ 設計） 2 目 錄 摘要 ABSTRACT…………………………………………………………………………………… 4 第 1章 緒論 ……………………………………………………………………………… … … 6 課題背景 … …………………………………………………… ………… … 6 課題研究的目的和意義 …… ………………………………………… …… 7 MATLAB 概述 …………………………………………………………… 8 國內(nèi)外研究現(xiàn)狀和發(fā)展趨勢 …… …………………………………… 9 論文的主要研究內(nèi)容 …… ………………………………………… 10 預期達到的目標 …… ……………………………………………… 11 第 2 章 濾波器的特性 數(shù)字濾波器的定義、分類及特點 …… ……………………………………… 12 FIR 濾波器的結(jié)構(gòu)特點 …… ……………………………………………… 13 FIR 線性相位濾波器的特點 …… ………………………………………… 14 實際濾波器的設計指標 …… ……………………………………………… 16 第 3 章 FIR 數(shù)字濾波器 主要設計方法的 比較 窗函數(shù)法 …… ……………………………………………………… 17 概述 …… ……………………………………………………… 17 優(yōu)缺點 …… ……………………………………………………… 18 頻率抽樣法 概述 …… ……………………………………………………… 19 優(yōu)缺點 …… …………………………………………………… 20 其它設計方法 …… ……………………………………………………… 22 設計方法的選定 …… …………………………………………………… 22 FIR 數(shù)字濾波器的最優(yōu)化設計 …… …………………………………… 23 第 4 章 FIR 濾波器的 MATLAB 設計 設計的基本思路 …… …………………………………………………… 25 設計的總流程圖 …… …………………………………… 26 窗函數(shù)的種類及理論特點 …… ………………………………………… 27 畢業(yè)論文（ 設計） 3 設計函數(shù) fir1 和 fir2…… …………………………………………… 28 FIR 數(shù)字濾波器的具體設計步驟 …… ………………………………… 29 線性相位 FIR 低通濾波器 的 設計 …… ………………………………… 31 帶通濾波器 的 設計 …… ………………………………………………… 37 概述 …… ……………………………………………………… 37 帶通濾波器的獲取 …… …………………………………… 40 結(jié)論分析 …… ……………………………………………………… 43 第 5 章 總結(jié)與展望 …… ……………………………………………………… 44 參考文獻 ……………………………………………………………………………………… 45 致謝 ……………………………………………………………………………………………… 46 畢業(yè)論文（ 設計） 4 摘 要 傳統(tǒng)的數(shù)字濾波器的設計過程復雜計算工作量大，濾波特性調(diào)整困難，影響了它的應用。它是一個離散時間系統(tǒng)，其功能 是對輸入離散信號的數(shù)字代碼進行運算處理，以達到改變信號頻譜的目的。數(shù)字濾波器按單位脈沖響應的性質(zhì)可分為無限長單位脈沖響應濾波器 IIR 和有限長單位脈沖響應 濾波器 (FIR)兩種?？煽啃愿?，具有模擬設備所沒有的許多優(yōu)點。它將數(shù)值分析、 矩陣計算 、科學數(shù)據(jù)可視化以及非 線性 動態(tài)系統(tǒng)的 建模 和仿真等諸多強大功能集成在一個易于使用的視窗環(huán)境中，為科學研究、工程設計以及必須進行有效數(shù)值計算的眾多科學 領域 提供了一種全面的解決方案，并在很大程度上擺脫了傳統(tǒng)非交互式程序設計語言（如 C、Fortran）的編輯模式，代表了當今國際科學計算軟件的先進水平。這里介紹了一種利用 MATLAB信號處理工具箱（ Signal Processing Toolbox）快速有效的設計由軟件組成的常規(guī)數(shù)字濾波器的設計方法。 4. 語法限制不嚴格，程序設計 自由度大 5. 程序的可移植性很好，基本上不做修改就可以在各種型號的計算機和操作系統(tǒng)上運行。功能性工具箱能用于多種學科。 因此，運用計算機對數(shù)字信號進行處理就顯得尤為重要， MATLAB 是一種面向科學和工程計 算的高級語言，現(xiàn)已成為國際公認的最優(yōu)秀的科技界應用軟件，在世界范圍內(nèi)廣為流行和使用。在軍事上廣泛用于導航、制導、電子對抗、戰(zhàn)場偵察；在電力系統(tǒng)中被應用于能源分布和自動檢測；在環(huán)境保護中被用于對空氣污染和噪聲干擾的自動監(jiān)控，在經(jīng)濟領域中被用于股票市場預測和經(jīng)濟效益分析等等。 (2) 明確研究的意義、課題來源及設計目標。 預期達到的目標 本研究課題主要研究目標是，利用 MATLAB 實現(xiàn)低通濾波器的設計，再通過頻率變換將低通濾波器轉(zhuǎn)換成希望類型的濾波器。 數(shù)字濾波器分類 按照不同的分類方法，數(shù)字濾波器有許多種類，但是總起來可以分為兩大類：經(jīng)典濾波器和現(xiàn)代濾波器。選擇 FIR濾波器還是IIR 濾波器取決于不同類型濾波器的優(yōu)點在設計中的重要性。在很多情況下， FIR數(shù)字濾波器的線性相位與它的高階數(shù)帶來的額外成本相比是非常值得的。 (5)從設計手段上看， IIR 濾波器可以利用模擬濾波器設計的成果，一般有大量有理函數(shù)的設計公式、曲線、圖表等可供計算、查找，設計簡單工作量較小。因此 FIR 濾波器可設計出理想正交變換器、理想微分器、線性調(diào)頻器等各種網(wǎng)絡，適應性很廣。有限長單位沖激響應（ FIR）濾波器有以下幾個特點 (1)系統(tǒng)的單位沖激響應 h(n)在有限個 數(shù) 值處不為零； (2)系統(tǒng)函數(shù) H(z)在 0z? 處收斂，在 0z? 處只有零點，有限 z平面只有零點，而全部極點都在 z=0 處（因果系統(tǒng)）； (3)結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu)，沒有輸出到輸入的反饋，但有些結(jié)構(gòu)中（例如頻率抽樣結(jié)構(gòu)）也包含有反饋的遞歸部分。 FIR 線性相位濾波器的特點 設 h(n) (0≤n≤M 1)是長度為 M的線性相位濾波器的脈沖響應，它的系統(tǒng)函數(shù)為 : H(z)= ?? ????????? ? 101)1(10 )()(MnnMMMnn znhzznh 它的頻率響應函數(shù)為 : ????? ???? ???? ,)()( 10 njMnj enheH 它具有一個線性相位約束 : ??????? ?????? ,)( jeH 其中， β =0或 ? π /2， α 為一個常數(shù)。 （ 3）線性相位（ linear phase） ， 造成 h(n)是偶對稱（ even）或奇對稱（ odd）且有限長。參數(shù) 定義了通帶波紋，及濾波器通帶內(nèi)偏移單位增益的最大值。 為了減少頻譜能量泄漏，可采用不同的截取函數(shù)對信號進行截短，截短函數(shù)稱為窗函數(shù)，簡稱為窗。窗函 數(shù)法由于簡單、物理意義清晰，因而得到了較為廣泛的應用 不同的窗函數(shù)對信號頻譜的影響是不一樣的，這主要是因為不同的窗函數(shù)，產(chǎn)生泄漏的大小不一樣，頻率分辨能力也不一樣。頻率取樣法是指在脈沖響應 h(n)為有限長度的條件下，根據(jù)頻域取樣定理，對所要求的頻率響應進行取樣，從樣點中恢復原來的頻率特性，達到設計 濾波器的目的 頻 率抽樣法從頻域出發(fā) ， 把給定的理想頻率響應 加 以等間隔抽 樣得到 Hd(k) 頻率取樣法先對理想頻響抽樣，得到樣值 H(k)[6]。前面已經(jīng)敘述過前兩種，現(xiàn)在說說切比雪夫等波紋逼近法 。它 采用 “ 最大誤差最小化 ” 優(yōu)化準則，即 min（ max|E（ ω） |） ，其中權函數(shù)誤差 E（ ω） =W（ ω） [Hd（ ω） H（ ω） ]，W（ ω） 為加權函數(shù) ， Hd（ ω） 為期望頻率響應， H（ ω） 為實際頻率響應。 FIR 數(shù)字濾波器的 最優(yōu)化設計 本節(jié)將介紹優(yōu)化設計的三種最優(yōu)化設計方法：均方誤差最小化準則、最大誤差最小化準則、切比雪夫最佳一致逼近。 這一方法注重的是整個 ??— 頻率區(qū)間內(nèi)總誤差的全局最小，但不能保證局部頻率點的性能，有些頻點可能會有較大的誤差。顯然，第二種濾波器的階數(shù)一定比第一種高。窗函數(shù)法的設計過程如下 ()( ) ( ) ( ) ( )I D TF T n D TF TjjddH e h n h n H e?? ? ????? ???? ???? 其中 )( ?jd eH 是理想濾波器的頻率響應， H( je? )是設計出的實際濾波器的頻響。1 ， ? ， 177。 2） fir2函數(shù)，其調(diào)用格式為： ? b=fir2(n,f,m) ? b=fir2(n,f,m,window) ? b=fir2(n,f,m,npt) ? b=fir2(n,f,m,npt,window) ? b=fir2(n,f,m,npt,lap) ? b=fir2(n,f,m,npt,lap,window) 格式中， n為濾波器階數(shù)； f和 m分別為濾波器期望幅頻響應的頻率向量和幅值向量，取值在 0～ 1間， m和 f長度相同； window為窗函數(shù)，列向量，長度必須是（ n+1）；默認 時自動取 hamming； npt為對頻率響應進行內(nèi)插點數(shù)，默認時為 512；lap定義一個區(qū)域尺寸，默認為 25； b為 FIR濾波器系數(shù)向量，長度為 n+1。 （ 3） 求濾波器的單位脈沖響應 h(n) h(n)= hd(n)?W(n) h(n)即為所設計的 FIR濾波器系數(shù)向量 h(n)式中 W(n)是上面選擇好的窗函數(shù)。 31 ( ) 2 4 1 0 ( / s e c )2c p s t r a d?? ? ? ? ? ? ? ? 其對應的數(shù)字頻率為 2 0 .3 7 5ccc ssf? ? ???? ? ?? 由此可得 ? ?11() 22 ccjnj j ndh n e e d e d?? ??? ? ? ???? ???????? ? ?1 sin ( ) ,(){,ccnnnn? ? ???? ??????? 其中， ? 為線性相位所必需的移位，我們已經(jīng)知道其滿足 12N? ?? （ 3） 由阻帶衰減 As來確定窗形狀，由過渡帶寬確定 N 畢業(yè)論文（ 設計） 33 由 As=40dB 可知，漢寧窗、海明窗、布萊克曼窗均滿足要求。 %定義窗函數(shù)的長度 N=M1。)。時間 /秒 39。 subplot(212) Fs=fft(s,512)。)。 畢業(yè)論文（ 設計） 35 ylabel(39。 %濾波后的信號頻域圖 AFsf=abs(Fsf)。 title(39。但是實際上，沒有 真正意義的理想帶通濾波器。 %采樣頻率 fcuts=[1000 1300 2210 2410]。 %濾波器的幅頻特性圖 [H,W]=freqz(b,a,M) %當 N 是整數(shù)時函數(shù)返回 N點的頻率向量和幅頻響應向量 figure(1) plot(f*fs/(2*pi),20*log10(abs(h)))。濾波器的增益響應 3