基于Matlab的FIR低通濾波器設(shè)計(jì)畢業(yè)論文

1、浙江萬里學(xué)院 本科畢業(yè)設(shè)計(jì)本科畢業(yè)設(shè)計(jì)( (論文論文) ) (2012 屆) 論文題目論文題目 基于 matlab的fir低通濾波器設(shè)計(jì) (英文) design of fir low-pass digital filter based on matlab 畢業(yè)設(shè)計(jì)（論文）原創(chuàng)性聲明和使用授權(quán)說明畢業(yè)設(shè)計(jì)（論文）原創(chuàng)性聲明和使用授權(quán)說明 原創(chuàng)性聲明原創(chuàng)性聲明 本人鄭重承諾：所呈交的畢業(yè)設(shè)計(jì)（論文） ，是我個(gè)人在指導(dǎo)教師 的指導(dǎo)下進(jìn)行的研究工作及取得的成果。盡我所知，除文中特別加以標(biāo) 注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果， 也不包含我為獲得 及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使

2、用過的 材料。對(duì)本研究提供過幫助和做出過貢獻(xiàn)的個(gè)人或集體，均已在文中作 了明確的說明并表示了謝意。 作 者 簽 名： 日 期： 指導(dǎo)教師簽名： 日期： 使用授權(quán)說明使用授權(quán)說明 本人完全了解 大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計(jì)（論文） 的規(guī)定，即：按照學(xué)校要求提交畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版本； 學(xué)校有權(quán)保存畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版，并提供目錄檢索與 閱覽服務(wù)；學(xué)?？梢圆捎糜坝?、縮印、數(shù)字化或其它復(fù)制手段保存論文； 在不以贏利為目的前提下，學(xué)校可以公布論文的部分或全部內(nèi)容。 作者簽名： 日 期： 學(xué)位論文原創(chuàng)性聲明學(xué)位論文原創(chuàng)性聲明 本人鄭重聲明：所呈交的論文是本人在導(dǎo)師的指導(dǎo)下獨(dú)立

3、進(jìn)行研究 所取得的研究成果。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包 含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫的成果作品。對(duì)本文的研究做出 重要貢獻(xiàn)的個(gè)人和集體，均已在文中以明確方式標(biāo)明。本人完全意識(shí)到 本聲明的法律后果由本人承擔(dān)。 作者簽名： 日期： 年 月 日 學(xué)位論文版權(quán)使用授權(quán)書學(xué)位論文版權(quán)使用授權(quán)書 本學(xué)位論文作者完全了解學(xué)校有關(guān)保留、使用學(xué)位論文的規(guī)定，同 意學(xué)校保留并向國家有關(guān)部門或機(jī)構(gòu)送交論文的復(fù)印件和電子版，允許 論文被查閱和借閱。本人授權(quán) 大學(xué)可以將本學(xué)位論文的 全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進(jìn)行檢索，可以采用影印、縮印或掃描 等復(fù)制手段保存和匯編本學(xué)位論文。 涉密論文按學(xué)校規(guī)

4、定處理。 作者簽名：日期： 年 月 日 導(dǎo)師簽名： 日期： 年 月 日 摘摘 要要 fir數(shù)字濾波器是數(shù)字信號(hào)處理的一個(gè)重要組成部分，由于fir數(shù)字濾波器具有 嚴(yán)格的線性相位，因此在信息的采集和處理過程中得到了廣泛的應(yīng)用。本文介紹了 fir數(shù)字濾波器的概念和線性相位的條件，分析了窗函數(shù)法、頻率采樣法和等波紋逼 近法設(shè)計(jì)fir濾波器的思路和流程。在分析三種設(shè)計(jì)方法原理的基礎(chǔ)上，借助matlab 仿真軟件工具箱中的fir1、fir2和remez子函數(shù)分別實(shí)現(xiàn)窗函數(shù)法、頻率采樣法和等波 紋逼近法設(shè)計(jì)fir濾波器。仿真結(jié)果表明，在相頻特性上，三種方法設(shè)計(jì)的fir濾波 器在通帶內(nèi)都具有線性相位；在幅頻特

5、性上，相比窗函數(shù)法和頻率采樣法，等波紋 逼近法設(shè)計(jì)fir濾波器的邊界頻率精確，通帶和阻帶衰減控制。 關(guān)鍵詞：關(guān)鍵詞：fir數(shù)字濾波器；窗函數(shù)法；頻率采樣法；等波紋逼近法；matlab abstract in a digital signal, there will be some interference signal which system does not needed often. so we need to use digital filtert to filter the signal to extract the useful signal. fir digital filter

6、 is an important component of digital filter.as the fir digital filters linear phase response is required, it have been widely used in information collection and processing. this paper introduced the concept of fir digital filters and the general idea of the design, introducting and compareding the

7、advantages and disadvantages of three kinds of the fir filter design method, and respectively design filter using the window function method, the frequency sampling method and wave approximation method. especially in digital signal processing, matlab is widly used which can design requirements it is

8、 intuitive and simple so that greatly reduce the workload. this article is designed by matlab function and simulated by matlab software to make the tedious step in the calculation simplify.it also drawed the amplitude-frequency characteristic curve relatly. finally, the original signal and filtered

9、signal are contrasted and analied. key words:fir digital filters，window function，frequency sampling，equiripple approximation method，matlab 目目 錄錄 1引言.1 2fir 數(shù)字濾波器線性相位條件.2 2.1 fir 數(shù)字濾波器概述.2 2.2 fir 數(shù)字濾波器線性相位定義.3 2.3 fir 數(shù)字濾波器線性相位時(shí)域約束條件.3 3fir 數(shù)字濾波器設(shè)計(jì)方法.5 3.1 fir 數(shù)字濾波器的窗函數(shù)設(shè)計(jì)法.5 3.1.1 窗函數(shù)設(shè)計(jì)法的設(shè)計(jì)思路.5 3.1.

10、2 吉布斯效應(yīng).7 3.1.3 常見窗函數(shù)介紹.7 3.2 fir 數(shù)字濾波器的頻率采樣設(shè)計(jì)法.10 3.2.1 頻率采樣法的基本思路.10 3.2.2 頻率采樣法的設(shè)計(jì)步驟.11 3.3 fir 數(shù)字濾波器的等波紋逼近設(shè)計(jì)法.13 3.4 不同設(shè)計(jì)方法的比較 .16 4 基于 matlab 的 fir 數(shù)字濾波器設(shè)計(jì).17 4.1 matlab 簡介.17 4.2 窗函數(shù)法的 matlab 實(shí)現(xiàn).18 4.2.1 fir1 函數(shù)介紹.18 4.2.2 基于 fir1 函數(shù)的窗函數(shù)法 fir 濾波器設(shè)計(jì).18 4.3 頻率采樣法的 matlab 實(shí)現(xiàn).21 4.4 等波紋逼近法的 matlab

11、 實(shí)現(xiàn).22 5 結(jié)論.24 致謝.25 參考文獻(xiàn).26 1引言引言 隨著信息科學(xué)和計(jì)算機(jī)技術(shù)的不斷發(fā)展，數(shù)字信號(hào)處理(dsp，digital signal processing)的理論和技術(shù)也得到了飛速的發(fā)展，并逐漸成為一門重要的學(xué)科，它 的重要性在日常通信、圖像處理、遙感、聲納、生物醫(yī)學(xué)、地震、消費(fèi)電子、國 防軍事、醫(yī)療方面等顯得尤為突出。在我們面臨的信息革命中，數(shù)字信號(hào)處理幾 乎涉及了所有的工程技術(shù)領(lǐng)域1。 數(shù)字信號(hào)處理是一種將信號(hào)以數(shù)字形式進(jìn)行處理的一種理論和技術(shù)，它的目 的是將真實(shí)世界中的一些信號(hào)進(jìn)行分析并濾波，最后得出其中的有用的信號(hào)。數(shù) 字濾波器是數(shù)字信號(hào)處理的一種，一般根據(jù)單位

12、脈沖響應(yīng)h(n)分為無限脈沖響應(yīng) (iir)和有限脈沖響應(yīng)(fir)系統(tǒng)。iir數(shù)字濾波器的設(shè)計(jì)方法簡單，特別是采用雙線 性變換法來設(shè)計(jì)的數(shù)字濾波器不存在頻域混疊的現(xiàn)象，但是iir濾波器存在一個(gè)較 為明顯的缺憾，就是它的相位響應(yīng)一般都是非線性的，而在傳輸頻帶內(nèi)的相位響 應(yīng)如果不是線性的，就會(huì)造成有用信號(hào)的傳輸失真，而fir數(shù)字濾波器不僅可以設(shè) 計(jì)成任意的幅度響應(yīng)，而且可以設(shè)計(jì)成在通頻帶內(nèi)具有良好的線性相位響應(yīng)。fir 數(shù)字濾波器的單位脈沖響應(yīng)h(n)有限長，所以fir數(shù)字濾波器是穩(wěn)定的，不存在穩(wěn) 定性的問題，且可以通過快速傅里葉變換(fft)的算法來實(shí)現(xiàn)信號(hào)濾波，大大的提 高的運(yùn)算效率。因此，

13、fir數(shù)字濾波器日益引起了人們的關(guān)注。 fir數(shù)字濾波器的設(shè)計(jì)方法有很多，比較常用的有窗函數(shù)設(shè)計(jì)法、頻率采樣設(shè) 計(jì)法、等波紋逼近法等。本課題通過運(yùn)用窗函數(shù)設(shè)計(jì)法、頻率采樣設(shè)計(jì)法和等波 紋逼近法來設(shè)計(jì)fir數(shù)字低通濾波器，并實(shí)現(xiàn)對(duì)給定的信號(hào)進(jìn)行濾波比較。其中窗 函數(shù)設(shè)計(jì)法是最基本的數(shù)字濾波方法，是利用傅里葉反變換(idtft)計(jì)算給定的頻 響的理想單位脈沖響應(yīng)，再加以窗函數(shù)進(jìn)行截?cái)嗪推交?。等波紋逼近法又稱切比 雪夫逼近法，是一種最小化最大誤差的頻域逼近方法2。matlab軟件的信號(hào)處理 工具箱提供了fir數(shù)字濾波器設(shè)計(jì)的子函數(shù)，運(yùn)用matlab軟件設(shè)計(jì)可以避免繁雜的 數(shù)學(xué)運(yùn)算，而且具有豐富的繪圖

14、功能，可以方便地查看所設(shè)計(jì)的數(shù)字濾波器的幅 度響應(yīng)和相位響應(yīng)是否滿足設(shè)計(jì)要求。因此，本課題在理論分析各種fir數(shù)字濾波 器設(shè)計(jì)方法的基礎(chǔ)上，運(yùn)用matlab軟件進(jìn)行仿真分析。 2fir 數(shù)字濾波器線性相位條件數(shù)字濾波器線性相位條件 2.1 fir 數(shù)字濾波器概述 一般來說一個(gè)經(jīng)典的數(shù)字濾波器是一個(gè)線性時(shí)不變系統(tǒng)，其數(shù)學(xué)模型可以用 z域系統(tǒng)函數(shù)來表示：( )h z （2-1） -1 0 1 ( ) 1 n r r r m k k k b z h z a z 其中均為濾波器參數(shù)。, kr a b n m 在（2-1）中，當(dāng)值不全為零值時(shí)，z域系統(tǒng)函數(shù)的必定含有一個(gè)或一 k a( )h z 個(gè)以上的

15、極值點(diǎn)，此時(shí)單位脈沖響應(yīng)為無限長，對(duì)于一個(gè)穩(wěn)定的數(shù)字濾波器來 說，z域系統(tǒng)函數(shù)必須在單位圓內(nèi)，因而把含有極值點(diǎn)的z域系統(tǒng)函數(shù)( )h z 的數(shù)字濾波器稱為無限脈沖響應(yīng)數(shù)字濾波器(infinite impulse response)，即( )h z iir數(shù)字濾波器。 而當(dāng)值全為零時(shí)，z域的系統(tǒng)函數(shù)只有一個(gè)零點(diǎn)， （2-1）表示的系統(tǒng) k a( )h z 函數(shù)可以寫成：( )h z （2-2） -1 0 ( ) n r r r h zb z 公式（2-2）表明，fir濾波器的系統(tǒng)函數(shù)是的階多項(xiàng)式，在有限 1 z(1)n 平面上有個(gè)零點(diǎn)，而在 平面原點(diǎn)處有階極點(diǎn)。z(0)n (1)n zz=0(1

16、)n （2-2）式表示的系統(tǒng)，其單位脈沖響應(yīng)可以表示為： （2-3） -1 0 ( )( )() n r r h ny nbnr 在（2-3）中，只有當(dāng)，才有非零值，所以數(shù)字濾波器的脈01nn( )h n 沖響應(yīng)是有限長的，因此在數(shù)字信號(hào)處理中把這種數(shù)字濾波器稱為有限脈沖響應(yīng) 數(shù)字濾波器(finite impulese response)，即fir數(shù)字濾波器。 fir數(shù)字濾波器最突出的兩個(gè)優(yōu)點(diǎn)是： （1）只要對(duì)附加一定的條件，就很容易獲得嚴(yán)格的線性相位。( )h n （2）由于的極點(diǎn)位于原點(diǎn)處，所以fir數(shù)字濾波器不存在穩(wěn)定性問( )h zz=0 題。 2.2 fir 數(shù)字濾波器線性相位定義

17、設(shè)fir數(shù)字濾波器脈沖響應(yīng)的長度為n，則其頻率響應(yīng)可以表示為： （2-4） 1 0 ()( ) n jj n n h eh n e 上式通過歐拉恒等式展開可得到的相位特性，有兩種線性相位特() j h e ( ) 性，通常稱為第一類線性相位和第二類線性相位。 第一類線性相位特性： 是一個(gè)與無關(guān)的常數(shù)( )= 第二類線性相位特性： 是起始相位 0 ( )= 0 嚴(yán)格地說第二種情況時(shí)的是不具有線性相位特性的，但上述兩種情況都( ) 滿足群延遲是一個(gè)常數(shù)，仍可以視為具有線性相位的，在第二類線性相位中 是常用的一種情況。 0 / 2 2.3 fir 數(shù)字濾波器線性相位時(shí)域約束條件 對(duì)于第一類線性相位，

18、即，通過一系列的運(yùn)算整理之后可得到一( )= 個(gè)三角函數(shù)求和公式： （2-5） 1 0 ( )sin()0 n n h nn 式中正弦函數(shù)為奇對(duì)稱，當(dāng)時(shí)，對(duì)稱中心為( )sin()h nn=(1)/ 2n ，需要滿足關(guān)于偶對(duì)稱，即要求：n=(1)/ 2n ( )h n(1)/ 2n ， （2-6）( )(1)h nh nn 10nn 對(duì)于第二類線性相位，即時(shí)，通過運(yùn)算得到公式：( )=- /2 （2-7） 1 0 ( )cos()0 n n h nn 函數(shù)為偶對(duì)稱，當(dāng)時(shí)，對(duì)稱中心也為( )cos()h nn=(1)/ 2n 。若要使上式成立，則要使關(guān)于奇對(duì)稱，即要求：(1)/ 2nn( )h

19、n(1)/ 2n ， （2-8）( )(1)h nh nn 10nn 從上述分析看來，線性相位fir數(shù)字濾波器的時(shí)域約束條件是指滿足線性相位 時(shí)對(duì)的約束條件，對(duì)于第一類線性相位，沖激響應(yīng)h(n)滿足（2-6）式；對(duì)于( )h n 第二類線性相位，沖激響應(yīng)h(n)滿足（2-8）式3-5。 3fir 數(shù)字濾波器設(shè)計(jì)方法數(shù)字濾波器設(shè)計(jì)方法 fir 數(shù)字濾波器的設(shè)計(jì)方法主要有窗函數(shù)設(shè)計(jì)法、頻率采樣設(shè)計(jì)法以及等波 紋逼近設(shè)計(jì)法三種，其中窗函數(shù)設(shè)計(jì)法是最常用的，其次是頻率采樣法，但這兩 種方法在設(shè)計(jì)中還會(huì)存在一些不足之處，所以需要優(yōu)化的設(shè)計(jì)方法，而等波紋逼 近法很好的彌補(bǔ)了窗函數(shù)法和頻率采樣法的不足6-7

20、。 3.1 fir 數(shù)字濾波器的窗函數(shù)設(shè)計(jì)法 3.1.1 窗函數(shù)設(shè)計(jì)法的設(shè)計(jì)思路 窗函數(shù)設(shè)計(jì)法是fir數(shù)字濾波器里最簡單的一種設(shè)計(jì)法，又叫傅里葉級(jí)數(shù)法， 為了設(shè)計(jì)簡單方便，通常選擇所希望逼近的濾波器的頻率響應(yīng)函數(shù)為具() j d he 有片段常數(shù)特性的理想濾波器，尋找一組，確定其頻率響應(yīng)( )h n ，然后用來逼近 1。窗函數(shù)法設(shè)計(jì)fir濾波器 1 0 ()( ) n jj n n h eh n e () j h e () j d he 是在時(shí)域中進(jìn)行的，那么可以通過傅里葉反變換得到得到頻率響應(yīng)，() j d he 即： （3-1） 1 ( )() 2 jjn dd hnheed 在實(shí)際中，一

21、般是處于逐段恒定的，在邊界頻率處有不連續(xù)點(diǎn)，因() j d he 而單位脈沖響應(yīng)是無限長的非因果序列，不能直接作為fir數(shù)字濾波器的單( ) d h n 位脈沖響應(yīng)，因此需要對(duì)進(jìn)行階段，轉(zhuǎn)換為有限長的一段因果序列，也就是( ) d h n 用一個(gè)有限長度的窗函數(shù)序列來截取，即，并將非因( )n( ) d h n( )= ( )*( ) d h nnh n 果序列轉(zhuǎn)變?yōu)橐粋€(gè)因果序列。截取的長度和加權(quán)窗函數(shù)的形狀都直接影響到逼近 精度。 窗函數(shù)法設(shè)計(jì)fir濾波器過程如圖3-1所示： 加窗： 計(jì)算 開始 正確逼近理想濾波器頻率響應(yīng)函數(shù) () j d he 求理想濾波器的單位脈沖響應(yīng)( ) d h n

22、 () j h e 選擇窗函數(shù)w(n)和窗長度n ( )( )*( ) d h nh nn 是否滿足要求 () j h e 計(jì)算h(n)或h(z) 結(jié)束 是 否 圖3-1 窗函數(shù)法設(shè)計(jì)fir濾波器流程 以截止頻率為，相位為零的理想低通濾波器為例，其頻率特性為： c （3-2） 1 () 0 cj d c he 通過傅里葉反變換得到對(duì)應(yīng)的為：( ) d h n （3-3） sin() ( ) c d n h n n 此時(shí)的是一個(gè)無限長的非因果序列，我們需要對(duì)其進(jìn)行截?cái)?，變成一個(gè)( ) d h n 有限長的因果序列?？梢韵劝严蛴移揭苽€(gè)點(diǎn)，得到為：( ) d h n(1)/ 2an( ) d hn

23、（3-4） 1 sin() 2 ( )() 1 () 2 c dd n n hnh na n n 相應(yīng)的傳輸函數(shù)為： () j d he （3-5）()()( ) jjj aj a ddd heh eehe 然后對(duì)截取從0到的n個(gè)點(diǎn)，n為窗函數(shù)的長度，所得的結(jié)果( ) d hn1n 表示為：( )h n （3-6）( )( )*( ) d h nnhn 表示窗函數(shù)，一般用下標(biāo)來表示窗函數(shù)的類型，矩形窗記為。( )n( ) r n 3.1.2 吉布斯效應(yīng) 用窗函數(shù)對(duì)進(jìn)行直接截?cái)?，得到有限長序列，并以代( ) d hn( )h n( )h n 替，肯定會(huì)引起誤差，表現(xiàn)在頻域就是通常所說的吉布斯(g

24、ibbs)效應(yīng)。對(duì)于( ) d h n 一個(gè)在有限區(qū)間分布的信號(hào)，其連續(xù)頻譜在頻域上分布往往是無限的，而在實(shí)際 信號(hào)處理時(shí)，我們通常只能在有限的區(qū)間內(nèi)做傅里葉分析，也就是說，我們只能 用有限區(qū)間來代替理論分析中的無限區(qū)間，多數(shù)情況下，我們總是選擇信號(hào)的低 頻部分，而舍棄高頻部分。而信號(hào)的高頻部分往往是反應(yīng)信號(hào)快速變化特征，如 果信號(hào)本身是連續(xù)的，這樣做一般不會(huì)引起信號(hào)的顯著變化，但實(shí)際中的信號(hào)往 往是比較豐富的，特別是信號(hào)本身存在劇烈變化的，這樣做必定會(huì)引起一些誤 差。該誤差引起過渡帶加寬以及通帶和阻帶內(nèi)的波動(dòng)8。為了減小吉布斯效應(yīng)帶 來的影響，需要調(diào)整窗口的長度來控制過渡帶的寬度，但要減小帶

25、內(nèi)的波動(dòng)以及 增大阻帶衰減，還需要從窗函數(shù)的形狀上尋找解決方法。 為了減少序列因截?cái)喽a(chǎn)生的gibbs效應(yīng)，窗函數(shù)在設(shè)計(jì)時(shí)需要注意： (1)頻率特性的主瓣要盡可能的窄，并且盡量把能量都集在主瓣內(nèi)。 (2)窗函數(shù)頻率特性的旁瓣趨于的過程中，其能量迅速減小為零。 3.1.3 常見窗函數(shù)介紹 常見的窗函數(shù)有矩形窗（rectangle window）、三角形窗（bartlerr window）、漢寧（hanning）窗升余弦窗、哈明（hamming）窗改進(jìn)的 升余弦窗、布萊克曼（blackman）窗、凱塞貝塞爾窗（kaiser-basel window）。 矩形窗的窗函數(shù)為： （3-7） 1 1 (

26、)2 0 otherwise r n n n 其頻譜的幅度函數(shù)為 （3-8） sin(/ 2) ( ) sin(/ 2) g r n w 矩形窗的主瓣寬度為，用矩形窗設(shè)計(jì)的fir數(shù)字濾波器的過渡帶寬度近4/ n 似為。1.8/ n 三角形窗的窗函數(shù)為： （3-9） 21 0n(1) 12 ( ) 21 2 (1)1 12 b n n n n n nnn n 其頻譜的幅度函數(shù)為 （3-10） 2 2sin(/ 4) ( ) sin(/ 2) g b n w n 三角窗的主瓣寬度為。8/ n 漢寧窗窗函數(shù)為 （3-11） 2 0.50.5cos 0 ( )1 0 h n nn nn otherwi

27、se 漢寧窗的頻譜的幅度函數(shù)為 （3-12） 22 ( )0.5( )0.25()() 11 hrrr wwww nn 漢寧窗的主瓣寬度為，漢寧窗在其兩個(gè)端點(diǎn)都為零，實(shí)際中這兩個(gè)端8/ n 點(diǎn)的數(shù)據(jù)是不可用的。 哈明窗的窗函數(shù)為 （3-13） 2 0.540.46cos 0 ( )1 0 h n nn nn otherwise 其幅度函數(shù)為 （3-14） 22 ( )0.54( )0.23()() 11 hrrr wwww nn 哈明窗是一種改進(jìn)的余弦窗，能量更加集中在主瓣，是一種高效的窗函數(shù)， 主瓣寬度與漢寧窗的相同。 布萊克曼窗窗函數(shù)為 （3-15） 241 (0.420.5cos0.8c

28、os) ( )112 0 b nnn n nnn otherwise 其頻譜的幅度函數(shù)為 （3-16） 22 ( )0.42( )0.25()() 11 44 0.04()() 11 brrr rr wwww nn ww nn 該窗函數(shù)位移不同，幅度函數(shù)也不同，會(huì)使旁瓣進(jìn)一步抵消，主瓣寬度為 。12/ n 凱塞窗是一種最優(yōu)窗函數(shù)，不同于前面五種窗函數(shù)，凱塞窗是一種參數(shù)可調(diào) 的窗函數(shù)，其函數(shù)形式如下： （3-17） ( ) ( ) 01 ( ) o k o i nnn i 其中 （3-18） 2 2 1 (1) 1 n n 2 1 1 ( )1() ! 2 k o k i k 一般取15-25項(xiàng)

29、可以滿足精度要求。參數(shù)可以控制窗的形狀。一般( ) o i 越大，主瓣越寬，而旁瓣幅度會(huì)隨之減小，典型的數(shù)據(jù)在4到9之間。 各種窗函數(shù)的性能比較如表3-1所示： 表3-1 不同窗函數(shù)性能比較 窗函數(shù)類型旁瓣峰值(db)阻帶最小衰減(db)過渡帶寬度(p/n) 矩形窗13214 三角窗25258 漢寧窗31448 續(xù)表3-1 哈明窗41538 布萊克曼窗577412 凱塞窗(5.6)51607.442 雖然窗函數(shù)設(shè)計(jì)法設(shè)計(jì)思路簡單，但是它的邊界頻率不容易控制，而且窗函 數(shù)還有吉布斯效應(yīng)，需要選擇不同的窗函數(shù)來減小吉布斯效應(yīng)對(duì)結(jié)果的影響，但 無論哪種窗函數(shù)，都無法很好的解決這一問題，所以我們需要通

30、過其他的設(shè)計(jì)方 法來進(jìn)行濾波，便于滿足實(shí)際工程中的不同要求。 3.2 fir 數(shù)字濾波器的頻率采樣設(shè)計(jì)法 3.2.1 頻率采樣法的基本思路 窗函數(shù)設(shè)計(jì)法是從時(shí)域出發(fā)來設(shè)計(jì)fir數(shù)字濾波器的，而頻率采樣法是從頻域 出發(fā)設(shè)計(jì)fir數(shù)字濾波器的。和窗函數(shù)設(shè)計(jì)法相同，頻率采樣法也需要預(yù)先構(gòu)造一 個(gè)希望逼近的濾波器頻率響應(yīng)函數(shù)，對(duì)其加以等間隔采樣后，作為fir數(shù)() j d he 字濾波器的頻率響應(yīng)。對(duì)在到之間等間隔采樣點(diǎn)，得到頻率采() j d he =02n 樣值：( ) d hk （3-19） 2 ( )() 0,1,2,.,1 j dd k n hkhekn 再對(duì)進(jìn)行點(diǎn)idft，得到：( ) d

31、 hkn( )h n （3-20） 1 0 1 ( )( ) 0,1,2,.,1 n kn dn k h nhk wnn n 將作為所涉及的fir數(shù)字濾波器的單位脈沖響應(yīng)，其系統(tǒng)函數(shù)為為( )h n( )h z （3-21） 1 0 ( )=( ) n n n h zh n z 由于濾波器頻率響應(yīng)是理想的，即有間斷點(diǎn)，那么其單位沖() j d he () j d he 激響應(yīng)是無限長的。這樣，由于時(shí)域混疊，引起所設(shè)計(jì)的h(n)和有偏)(nhd)(nhd 差。 因此，采樣點(diǎn)處與相等，逼近誤差為0，而在采樣()(2/) k j k h ek n ( )h k 點(diǎn)之間，由有限項(xiàng)的之和形成。其誤差和特

32、性() j h e ( )2/h kk n （）() j d he 的平滑程度有關(guān)，特性愈平滑誤差愈?。惶匦郧€間斷點(diǎn)處，誤差越大。誤差表 現(xiàn)形式為間斷點(diǎn)用傾斜線取代，且間斷點(diǎn)附近形成振蕩特性，使阻帶衰減減小， 往往不能滿足實(shí)際工程中的技術(shù)要求。當(dāng)然，增大n值，可以減小逼近誤差，但 間斷點(diǎn)附近誤差仍然最大，且n太大會(huì)增加濾波器級(jí)數(shù)與成本。提高阻帶衰減最 有效的方法是在頻響間斷點(diǎn)附近區(qū)間內(nèi)插一個(gè)或幾個(gè)過渡采樣點(diǎn)，使不連續(xù)點(diǎn)變 成緩慢過渡。 過渡帶采樣點(diǎn)個(gè)數(shù)與阻帶最小衰減的關(guān)系以及使阻帶最小衰減最大化的 s s 每個(gè)過渡帶采樣值求解都要用優(yōu)化算法解決。其基本思路是將過渡帶采樣值設(shè)為 一個(gè)自由量，用

33、一種優(yōu)化算法改變它們，最終使阻帶最小衰減最大。將過渡帶 s 采樣點(diǎn)的個(gè)數(shù)m與濾波器阻帶最小衰減的經(jīng)驗(yàn)數(shù)據(jù)列于表3-2中，我們可以根據(jù) s 給定的阻帶最小衰減，選擇過渡帶采樣點(diǎn)的個(gè)數(shù)m。 s 表3-2 過渡帶采樣點(diǎn)的個(gè)數(shù)m與濾波器阻帶最小衰減的經(jīng)驗(yàn)數(shù)據(jù) s m123 s 4454db6575db8595db 3.2.2 頻率采樣法的設(shè)計(jì)步驟 首先根據(jù)阻帶最小衰減按照表3-2選擇過渡帶采樣點(diǎn)的個(gè)數(shù)，再確定過 s m 渡帶寬度，估算頻域采樣點(diǎn)數(shù)，如果增加個(gè)過渡帶采樣點(diǎn)，則過渡帶寬度 t bnm 近似變成。當(dāng)確定時(shí)，過渡帶會(huì)隨著的增大而變寬。如果給定(1)2/mnnm 的過渡帶寬度為，則要求，濾波器的

34、長度必須滿足以下公 t b(1)2/ t mnbn 式： （3-22） 2 (1) t nm b 接著，構(gòu)造一個(gè)希望逼近的頻率響應(yīng)函數(shù)： （3-23） (1)/2 ()( ) jjn ddg hehe 設(shè)計(jì)標(biāo)準(zhǔn)型片段常數(shù)特性的fir數(shù)字濾波器時(shí)，一般構(gòu)造幅度特性函數(shù) 為相應(yīng)的理想頻響特性，且滿足的對(duì)稱情況。( ) dg h( )h n 對(duì)（3-23）進(jìn)行頻域采樣： （3-24） 1 2 ( )()( ) 0,1,2,.,1 n jk j n dg k n h khehk ekn （3-25） 2 ( )() k=0,1,2,.,n-1 gdg hkhk n 并加入過渡帶采樣。過渡帶采樣值一般為

35、經(jīng)驗(yàn)值，或者用累試法確定，也可 以采用優(yōu)化算法估算。 對(duì)進(jìn)行點(diǎn)idft，得到第一類線性相位fir數(shù)字濾波器的單位脈沖響( )h kn 應(yīng)： （3-26） 1 0 1 ( )( )( ) n=0,1,2,.,n-1 n kn n k h nidft h kh k w n 求出： 開始 確定所逼近的 ，估算采樣點(diǎn)數(shù)n與過渡帶采樣點(diǎn)數(shù)m() j d he 采樣： 2/ ( )(),1,2,.,1 jk n dd hkhekn 插入過渡帶采樣點(diǎn)數(shù) ( )( ),0,1,2,.1 d h nidft hknn 求出：() ( ) j h efft h n 是否滿足要求() j h e 輸出設(shè)計(jì)結(jié)果 或

36、( )h n( )h z 結(jié)束 增加過渡帶采樣點(diǎn)數(shù) 改變n的值 是 否 圖3-2 頻率采樣法設(shè)計(jì)fir濾波器流程 最后檢驗(yàn)設(shè)計(jì)結(jié)果，如果阻帶最小衰減未達(dá)到指標(biāo)要求，則要改變過渡帶的 采樣值，直到滿足指標(biāo)要求為止。如果濾波器的邊界頻率未達(dá)到指標(biāo)要求，則需 要微調(diào)的邊界頻率。頻率采樣法設(shè)計(jì)fir濾波器流程如圖3-2所示。( ) dg h 頻率采樣法最大的優(yōu)點(diǎn)就是直接從頻率進(jìn)行設(shè)計(jì)，比較直觀，也適合于設(shè)計(jì) 具有任意幅度特性的濾波器。但是頻率采樣法在邊界頻率不容易控制，如果采樣 點(diǎn)數(shù)n增加，對(duì)確定邊界頻率有好處，但同樣會(huì)增加了濾波器的成本，因此只適 合窄帶濾波，且這種設(shè)計(jì)方法理解起來比較困難。 3.3

37、 fir 數(shù)字濾波器的等波紋逼近設(shè)計(jì)法 窗函數(shù)設(shè)計(jì)法和頻率采樣設(shè)計(jì)法雖然設(shè)計(jì)方法簡單，但都存在濾波器邊緣頻 率不易精確控制缺點(diǎn)，且這兩種設(shè)計(jì)方法設(shè)計(jì)出來的濾波器的通帶和阻帶的波動(dòng) 幅度都是相等的，兩種設(shè)計(jì)方法都不能分別控制通帶和阻帶的波動(dòng)幅度，而現(xiàn)實(shí) 工程中往往對(duì)二者都有不同的要求，需要分別進(jìn)行控制。 等波紋逼近法是一種優(yōu)化設(shè)計(jì)方法，它克服了窗函數(shù)設(shè)計(jì)法和頻率采樣法的 缺陷，是最大誤差最小化設(shè)計(jì)方法，并在整個(gè)逼近頻段上均勻分布。 設(shè)為希望逼近的幅度特性函數(shù)，且要求設(shè)計(jì)線性相位的fir數(shù)字濾波器( ) d h 時(shí)，必須滿足線性相位約束條件。用表示實(shí)際設(shè)計(jì)的幅度特性函( ) d h( ) g h

38、數(shù)，定義加權(quán)誤差函數(shù)為( )e （3-27）( )( )( )( ) dg ewhh 式中，被稱為誤差加權(quán)函數(shù)，是由設(shè)計(jì)者定義的，用來控制不同( )0w 頻段的逼近精度。經(jīng)過推導(dǎo)可把統(tǒng)一標(biāo)示為：( ) g h ( )=q( )p( ) g h 式中，是系數(shù)不同的余弦組合式，記；p( ) n nnap)cos()()( 是不同的常數(shù)，在設(shè)計(jì)fir濾波器時(shí)存在四種線性相位，當(dāng)q( ) 且奇對(duì)稱時(shí)，n為奇數(shù)，為1；n為偶數(shù)時(shí)，為( )(1)h nh nnq( )q( ) ；而當(dāng)偶對(duì)稱時(shí)，不管n為奇數(shù)還是偶數(shù)，都取cos(/ 2)( )(1)h nh nn q( ) 。sin 開始 設(shè)定加權(quán)函數(shù)及通帶

39、波紋峰值 ( )e( )e ( )p 設(shè)定初始頻率 計(jì)算 ，并求出 的新極值點(diǎn) 多于 個(gè)極值 選 個(gè) 最大極值點(diǎn) 得到最優(yōu)化的 個(gè)極值 計(jì)算 的系數(shù) 結(jié)束 是 否 否 是 1n 1n 1n 極值點(diǎn)位置是否改變 圖3-3 remez算法搜索最佳頻率點(diǎn)流程 等波紋逼近的問題就是選擇n個(gè)系數(shù)a(n)，使加權(quán)誤差的最大值最小。( )e 設(shè)定誤差加權(quán)函數(shù)、通帶波紋峰值及n+1個(gè)初始頻率i,i=0,1,2,n，計(jì)算各 個(gè)頻率點(diǎn)上的加權(quán)誤差函數(shù)值，如果這些加權(quán)誤差函數(shù)值小于等于給定的通帶波 紋峰值，則為最佳頻率點(diǎn)；如果某個(gè)頻率點(diǎn)上的加權(quán)誤差函數(shù)值大于給定的通帶 波紋峰值，就用remez算法更新迭代查找最佳頻

40、率點(diǎn)。找到最佳頻率點(diǎn)后，計(jì)算 系數(shù)及。remez算法搜索最佳頻率點(diǎn)流程如圖3-3所示。p( )p( ) 借助remez算法來設(shè)計(jì)等波紋的fir數(shù)字濾波器的步驟如下： 給出所需的頻率響應(yīng)、誤差加權(quán)函數(shù)和單位脈沖響應(yīng)的長( ) d h( )w( )h n 度n，根據(jù)具體情況，進(jìn)行統(tǒng)一表達(dá)式的轉(zhuǎn)換；設(shè)定初始頻率點(diǎn)和通帶波( ) g h 紋值，采用remez算法計(jì)算最佳頻率點(diǎn)組，由最佳頻率點(diǎn)組計(jì)算及系( )p( )p 數(shù)，根據(jù)求出的的表達(dá)式再利用傅里葉逆變換求得濾波器的單位脈沖響應(yīng)( )p 。用等波紋逼近法設(shè)計(jì)線性相位fir濾波器的流程如圖3-4所示：( )h n 輸入濾波器技術(shù)要 求： 設(shè)定誤差加權(quán)

41、函數(shù) 及通帶波紋峰值 ,( ) d n h 給出n+1個(gè)交錯(cuò)點(diǎn)組頻率初 始值：,0,1,2,. i in 調(diào)用remez算法程序求解 最佳極值頻率，計(jì)算 ( )p 計(jì)算單位脈沖響應(yīng)( )h n 輸出最佳誤差和 ( )h n 圖3-4 等波紋逼近法設(shè)計(jì)fir濾波器流程 利用等波紋逼近法設(shè)計(jì)fir濾波器，其誤差均勻分布在頻帶中，可以得優(yōu)良的 濾波特性，它在同樣過渡帶較窄的情況下，通帶最穩(wěn)定，阻帶有最大化的最小衰 減。 3.4 不同設(shè)計(jì)方法的比較 窗函數(shù)設(shè)計(jì)法是這三種設(shè)計(jì)法里操作最簡單的，但是窗函數(shù)的設(shè)計(jì)階數(shù)相對(duì) 其他兩種更多，如果降低了設(shè)計(jì)的階數(shù)，則最后得到的設(shè)計(jì)結(jié)果往往不能達(dá)到實(shí) 際設(shè)計(jì)的要求和

42、指標(biāo)，除矩形窗之外，其他窗可以滿足一般的設(shè)計(jì)指標(biāo)，但在要 求精度較高的時(shí)候，窗函數(shù)設(shè)計(jì)往往不能很好的達(dá)到要求。 頻率采樣法設(shè)計(jì)濾波器相對(duì)于窗函數(shù)，其階數(shù)稍微比窗函數(shù)要少，但要通過 優(yōu)化過渡帶的采樣點(diǎn)取值來控制阻帶的波紋幅度，所以并不能很好的滿足工程中 的要求。 而等波紋逼近法是這三種方法中階數(shù)最少的，且通帶最大衰減是最小的，阻 帶衰減最大。 4 基于基于 matlab 的的 fir 數(shù)字濾波器設(shè)計(jì)數(shù)字濾波器設(shè)計(jì) 4.1 matlab 簡介 matlab 是 matrix laboratory（“矩陣實(shí)驗(yàn)室”）的縮寫，是由美國 mathworks 公司開發(fā)的集數(shù)值計(jì)算、符號(hào)計(jì)算和圖形可視化三大基

43、本功能于一體 的，功能強(qiáng)大、操作簡單的語言。matlab 是國際公認(rèn)的優(yōu)秀數(shù)學(xué)應(yīng)用軟件之一。 它集中了日常數(shù)學(xué)處理中的各種功能，包括高校的數(shù)值計(jì)算、矩陣運(yùn)算、信號(hào)處 理和圖像生成等功能。在 matlab 環(huán)境下，用戶可以進(jìn)行程序設(shè)計(jì)、數(shù)值計(jì)算、圖 形繪制、輸入輸出、文件管理等各項(xiàng)操作。除此之外，matlab 易于擴(kuò)充。除內(nèi)部 函數(shù)外，所有 matlab 的核心文件和工具箱文件都是可讀可改的源文件，用戶可修 改源文件和加入自己的文件，它們可以與庫函數(shù)一樣被調(diào)用9。 matlab是一種矩形運(yùn)算為基礎(chǔ)的交互式程序語言，著重針對(duì)科學(xué)計(jì)算、工程 計(jì)算和繪畫的需求，與其他機(jī)器語言相比，其特點(diǎn)是簡單和智能化

44、，適應(yīng)科技專 業(yè)人員的思維方式和書寫習(xí)慣，使得編程和調(diào)試效率大大提高。 matlab由一系列工具組成。這些工具方便用戶使用 matlab的函數(shù)和文 件，其中許多工具采用的是圖形用戶界面。包括 matlab桌面、命令窗口、歷 史命令窗口、編輯器和調(diào)試器、路徑搜索和用于用戶瀏覽幫助、工作空間、 文件的瀏覽器。隨著 matlab的商業(yè)化以及軟件本身的不斷升級(jí)， matlab的用 戶界面也越來越精致，更加接近 windows的標(biāo)準(zhǔn)界面，人機(jī)交互性更強(qiáng)，操作 更簡單。matlab自產(chǎn)生之日起就具有方便的數(shù)據(jù)可視化功能，將向量和矩陣用 圖形表現(xiàn)出來，并且可以對(duì)圖形進(jìn)行標(biāo)注和打印。高層次的作圖包括二維和 三

45、維的可視化、圖象處理、動(dòng)畫和表達(dá)式作圖，可用于科學(xué)計(jì)算和工程繪圖。 matlab 具有功能強(qiáng)大的工具箱，工具箱可分為兩類：功能性工具箱和學(xué)科性 工具箱。功能性工具箱主要用來擴(kuò)充其符號(hào)計(jì)算功能、圖示建模仿真功能、文字 處理功能以及與硬件實(shí)時(shí)交互的功能。而學(xué)科性工具箱是專業(yè)性比較強(qiáng)的，如優(yōu) 化工具箱、統(tǒng)計(jì)工具箱、控制工具箱、小波工具箱、圖象處理工具箱、通信工具 箱等。 4.2 窗函數(shù)法的 matlab 實(shí)現(xiàn) matlab 信號(hào)處理工具箱提供了基于加窗的線性相位 fir 濾波器設(shè)計(jì)函數(shù) fir1。 4.2.1 fir1 函數(shù)介紹 fir1 調(diào)用格式為：b=fir1(n, wc,ftype,windo

46、w) 函數(shù)參數(shù)說明如下： 1.n 表示濾波器的階數(shù) 2. wc 為濾波器的歸一化截止頻率，它是一個(gè)大于 0 小于 1 的一個(gè)數(shù) 3.ftype表示所設(shè)計(jì)濾波器的類型，如果 ftype=high，則表示高通濾波器，如 果 ftype=stop，則表示帶阻濾波器，如果此時(shí)沒有參數(shù)，就表示低通濾波器 4.window 表示的是指定的窗函數(shù)，如矩形窗為 rectwin(n)，三角窗為 bartlett(n)， 如果缺省 window 參數(shù)，則 fir1 默認(rèn)為是哈明窗 hamming(n)。 4.2.2 基于 fir1 函數(shù)的窗函數(shù)法 fir 濾波器設(shè)計(jì) 下面給出利用 fir1 函數(shù)的窗函數(shù)法數(shù)字低通

47、濾波器。利用 fir1 函數(shù)進(jìn)行設(shè)計(jì)， 這種設(shè)計(jì)方法只需要給出濾波器的階數(shù)、截止頻率、窗函數(shù)等參數(shù)，matlab 即可 自行完成設(shè)計(jì)，并可通過 freqz 函數(shù)查看濾波器的幅頻響應(yīng)和相頻響應(yīng)，以驗(yàn)證濾 波器是否滿足設(shè)計(jì)要求。設(shè)計(jì)程序如下： n=27; %設(shè)濾器的階數(shù)為 27 wn=0.2; %歸一化截止頻率為 0.2 window=rectwin(28); %矩形窗 b=fir1(n,wn,window); freqz(b); 運(yùn)行結(jié)果如圖 4-1 所示： 從幅頻響應(yīng)上看，通帶波動(dòng)約為 1db 且?guī)?nèi)有過沖現(xiàn)象，0.2 截止頻率處的 衰減為 6.5db，即截止頻率發(fā)生了漂移，3db 截止頻率為

48、 0.182；阻帶波紋較大， 因而阻帶特性較不理想，阻帶最小衰減約為 21db；相頻響應(yīng)曲線在通帶內(nèi)為直線， 效果較好，信號(hào)失真小。 01 -600 -400 -200 0 normalized frequency ( rad/sample) phase (degrees) 01 -100 -50 0 50 normalized frequency ( rad/sample) magnitude (db) 圖 4-1 矩形窗 fir 濾波器幅頻和相頻特性曲線 利用三角窗進(jìn)行設(shè)計(jì)時(shí)，原理與

49、矩形窗基本相同，只不過生成窗函數(shù)時(shí)采用 triang（n）函數(shù)生成三角窗，程序運(yùn)行結(jié)果如圖 4-2 所示： 01 -2500 -2000 -1500 -1000 -500 0 normalized frequency ( rad/sample) phase (degrees) 01 -100 -50 0 50 normalized frequency ( rad/sample) magnitude (db) 圖 4-2 三角窗設(shè)計(jì)的 fir 濾波器幅頻和相頻特性曲線 從幅頻響應(yīng)上看，通帶

50、平坦沒有波動(dòng)，阻帶最小衰減約為 25db；相頻響 應(yīng)曲線在整個(gè)頻率范圍內(nèi)為直線，滿足線性相位要求。 利用哈明窗進(jìn)行設(shè)計(jì)時(shí)，原理與矩形窗相同，只不過生成窗函數(shù)時(shí)采用 hamming（n）函數(shù)生成三角窗，程序運(yùn)行結(jié)果如下： 01 -800 -600 -400 -200 0 normalized frequency ( rad/sample) phase (degrees) 01 -150 -100 -50 0 50 normalized frequency ( rad/sample) mag

51、nitude (db) 圖 4-3 哈明窗設(shè)計(jì)的 fir 濾波器幅頻和相頻特性曲線 利用布萊克曼窗進(jìn)行設(shè)計(jì)時(shí)，原理與矩形窗基本相同，只不過生成窗函數(shù) 時(shí)采用 blackman（n）函數(shù)生成布萊克曼窗，程序運(yùn)行結(jié)果如下： 01 -1500 -1000 -500 0 normalized frequency ( rad/sample) phase (degrees) 01 -200 -150 -100 -50 0 normalized frequency ( rad/sample) magn

52、itude (db) 圖 4-4 布萊克曼窗設(shè)計(jì)的 fir 濾波器幅頻和相頻特性曲線 利用凱塞窗進(jìn)行設(shè)計(jì)時(shí)，原理與矩形窗基本相同，只不過生成窗函數(shù)時(shí)采用 kaiser（n,beta）函數(shù)生成凱塞窗，程序運(yùn)行結(jié)果： 01 -800 -600 -400 -200 0 normalized frequency ( rad/sample) phase (degrees) 01 -150 -100 -50 0 50 normalized frequency ( rad/sample) magnit

53、ude (db) 圖 4-5 凱塞窗設(shè)計(jì)的 fir 濾波器幅頻和相頻特性曲線 4.3 頻率采樣法的 matlab 實(shí)現(xiàn) 頻率采樣法設(shè)計(jì)fir濾波器流程如圖3-2所示，因此利用matlab軟件設(shè)計(jì)fir濾 波器時(shí)，可以根據(jù)3-2的流程編寫程序?qū)崿F(xiàn)。但matlab軟件工具箱提供了fir2子函 數(shù)，可以利用fir2子函數(shù)快速方便地設(shè)計(jì)fir濾波器，避免了繁雜的數(shù)學(xué)運(yùn)算。 fir2 子函數(shù)調(diào)用格式為：b=fir2(n, f, a) 函數(shù)參數(shù)說明如下： 1.n 表示濾波器的階數(shù) 2. f 為濾波器的歸一化頻率點(diǎn)，在0,1范圍內(nèi)取值，1 表示實(shí)際頻率為 fs/2 3.a 為頻率響應(yīng)幅度 參數(shù)f和a共同描述

54、了預(yù)先設(shè)定的理想濾波器的頻率響應(yīng)，可以借助plot函數(shù) 畫圖觀察幅頻特性曲線。 下面用頻率采樣法設(shè)計(jì)一個(gè)通帶截止頻率為0.2的fir數(shù)字低通濾波器，設(shè)濾 波器階數(shù)n=27，編寫程序： n=27; %設(shè)濾器的階數(shù)為33 f=0 0.2 0.3 1; %理想濾波器的頻率響應(yīng) a=1 0. 707 0 0; b=fir2(n,f,a); freqz(b); 經(jīng)過運(yùn)算得fir濾波器的幅頻特性如圖4-6所示： 圖4-6 頻率采樣法設(shè)計(jì)的fir低通濾波器的幅頻和相頻特性曲線 從4-6幅頻特性曲線可知，0.2截止頻率處的衰減為4.7db，即截止頻率發(fā)生 了漂移，3db截止頻率為0.165；通帶內(nèi)的相頻特性為

55、直線，滿足線性相位要 求。 4.4 等波紋逼近法的 matlab 實(shí)現(xiàn) 頻率采樣法設(shè)計(jì)fir濾波器流程如圖3-4所示，因此利用matlab軟件設(shè)計(jì)fir濾 波器時(shí)，可以根據(jù)3-4的流程編寫程序?qū)崿F(xiàn)。但matlab軟件工具箱提供了remez子 函數(shù)，可以利用remez子函數(shù)快速方便地設(shè)計(jì)fir濾波器，避免了繁雜的數(shù)學(xué)運(yùn) 算。 remez 函數(shù)有幾種調(diào)用格式，常見的一種調(diào)用格式為：b= remez (n, f, a, w); 函數(shù)參數(shù)說明如下： 1.n 表示濾波器的階數(shù) 2. f 為濾波器的歸一化頻率點(diǎn)，在0,1范圍內(nèi)取值，1 表示實(shí)際頻率為 fs/2 3.a 為頻率響應(yīng)幅度 參數(shù)f和a共同描述了

56、預(yù)先設(shè)定的理想濾波器的頻率響應(yīng)。 4.w表示通帶和阻帶波紋控制。 下面用 remez 函數(shù)設(shè)計(jì)一個(gè) 27 階的 fir 低通濾波器，其通帶截止頻率為 0.2，通帶波紋最大為 0.04，阻帶截止頻率為 0.3，阻帶波紋最大為 0.02，程序 如下： n=27; f=0 0.2 0.3 1; a=1 1 0 0; w=0.04 0.02; b=remez(n,f,a,w); freqz(b); 程序運(yùn)行結(jié)果如圖 4-7 所示： 01 -800 -600 -400 -200 0 normalized frequency ( rad/sample)

57、 phase (degrees) 01 -100 -50 0 50 normalized frequency ( rad/sample) magnitude (db) 圖4-7 等波紋逼近法設(shè)計(jì)的fir低通濾波器幅頻和相頻特性曲線 同樣是設(shè)計(jì)一個(gè) 27 階的 fir 低通數(shù)字濾波器，比較 fir1 窗函數(shù)法、fir2 頻率 采樣法和 remez 等波逼近法的設(shè)計(jì)結(jié)果，可以看出： （1）窗函數(shù)法和頻率采樣法設(shè)計(jì)濾波器時(shí)，濾波器邊緣頻率不易精確控制， 通帶和阻帶的波動(dòng)幅度 也不易控制。 （2）利用等波紋逼近法可以較好的控制濾波器的邊緣頻率及通帶

58、和阻帶的波動(dòng)幅度。 5 結(jié)論結(jié)論 論文主要分析了 fir 數(shù)字濾波器的基本理論，討論了 fir 數(shù)字濾波器的線性 相位種類及其約束條件，分析了窗函數(shù)設(shè)計(jì)法、頻率采樣設(shè)計(jì)法、等波紋逼近法 三種不同的設(shè)計(jì)方法，并借助 matlab 軟件對(duì)這三種設(shè)計(jì)方法進(jìn)行了對(duì)比分析。 窗函數(shù)設(shè)計(jì)法對(duì)信號(hào)加窗之后會(huì)使不連續(xù)點(diǎn)處邊沿加寬形成過渡帶，其寬度 （兩肩峰之間的寬度）等于窗函數(shù)頻率響應(yīng)的主瓣寬度。在處出現(xiàn)2/ c n 肩峰值，兩側(cè)形成起伏振蕩，振蕩的幅度和多少取決于旁瓣的幅度和多少，改變 n 只能改變窗譜的主瓣寬度，但不能改變主瓣與旁瓣的相對(duì)比例，其相對(duì)比例由 窗函數(shù)形狀決定。 窗函數(shù)設(shè)計(jì)法設(shè)計(jì) fir 數(shù)字濾波器是傅里葉變換的典型運(yùn)用，而頻率采樣法 設(shè)計(jì)的指導(dǎo)思想是通過頻域采樣點(diǎn)實(shí)現(xiàn)的，同時(shí)為了保證過渡帶寬的不變，濾波 器的采樣點(diǎn)數(shù)也要相應(yīng)增加，這就要