《MATLAB的 FIR數(shù)字濾波器》由會(huì)員分享，可在線閱讀，更多相關(guān)《MATLAB的 FIR數(shù)字濾波器（27頁(yè)珍藏版）》請(qǐng)?jiān)谘b配圖網(wǎng)上搜索。

1、畢業(yè)論文 基于 MATLAB的 FIR數(shù)字濾波器設(shè)計(jì) [摘要]數(shù)字濾波器是由數(shù)字乘法器、加法器和延時(shí)單元組成的一種裝置。其功能是對(duì)輸入的數(shù)字信號(hào)進(jìn)行運(yùn)算處理，以達(dá)到改變信號(hào)頻譜的目的。隨著電子計(jì)算機(jī)技術(shù)和大規(guī)模集成電路的發(fā)展，數(shù)字濾波器既可用計(jì)算機(jī)軟件實(shí)現(xiàn)，也可用大規(guī)模集成數(shù)字硬件實(shí)時(shí)實(shí)現(xiàn)。目前常用的 FIR 濾波器的設(shè)計(jì)方法主要有三種：窗函數(shù)法、頻率采樣法和最優(yōu)化設(shè)計(jì)法。本文針對(duì)三種設(shè)計(jì)方法基于 MATLAB進(jìn)行 FIR各種形式濾波器的設(shè)計(jì)與仿真，并比 較了三種方法的特點(diǎn)。結(jié)果表明，在同樣的設(shè)計(jì)指標(biāo)下，利用等波紋切比雪夫逼近法則的設(shè) 計(jì)可以獲得最佳的頻率特性和衰耗特性，具有通帶和阻帶

2、平坦，過渡帶窄等優(yōu)點(diǎn)。 [關(guān)鍵詞] FIR 濾波器 MATLAB 窗函數(shù)法 頻率采樣法 等波紋切比雪夫逼近法 Design of FIR Digital Filter Based on MATLAB [Abstract] Digital filter is a kind of installation that consists of digital multiplier, adder and the unit of delay time. It has a funct

3、ion of operating the inputting digital signal in order to change signal frequency spectra. With the development of the computer technology and LSI, digital filter can be achieved by computer software, and can be achieved by digital hardware real time of large scale integration too. At present, there

4、’re three methods to design FIR digital filter: window method, frequency sampling method and optimal design method. In this paper we use three methods to design and simulate various FIR digital filter based on MATLAB, and compares the characteristics of them. The result reveals that Chebyshev approx

5、imation method has a great deal of advantages including optimal frequency response and attenuate response, smooth pass band、forbidden band and narrow transition band under the same design order. [Key words] FIR digital filter MATLAB Window method requency sampling method Chebyshev

6、 approximation method 引言 隨著信息時(shí)代和數(shù)字世界的到來，數(shù)字信號(hào)處理己成為當(dāng)今極其重要的學(xué)科和技術(shù)領(lǐng)域。數(shù)字信號(hào)處理在通信、語(yǔ)音、圖像、自動(dòng)控制雷達(dá)、軍事、航空航天、醫(yī)療和家用電器等眾多領(lǐng)域得到了廣泛的應(yīng)用。數(shù)字濾波器是數(shù)字信號(hào)處理的重要基礎(chǔ)，在對(duì)信號(hào)的濾波、檢測(cè)及參數(shù)的估計(jì)等信號(hào)應(yīng)用中，數(shù)字濾波器是使用最為廣泛的一種線性系統(tǒng)。 在許多數(shù)字信號(hào)處理系統(tǒng)中，F(xiàn)IR濾波器是最常用的組件之一，它完成信號(hào)預(yù)調(diào)、頻帶選擇和濾波等功能。FIR濾波器在截止頻率的邊沿陡峭性能雖然不及II

7、R濾波器，但是，考慮到FIR濾波器嚴(yán)格的線性相位特性和不像IIR濾波器存在穩(wěn)定性的問題，F(xiàn)IR濾波器能夠在數(shù)字信號(hào)處理領(lǐng)域得到廣泛的應(yīng)用。 本文主要研究了FIR數(shù)字濾波器設(shè)計(jì)方法，從數(shù)字濾波器的概念及基本原理的分析，得出數(shù)字濾波器設(shè)計(jì)的基本步驟，并探究了常用的三種設(shè)計(jì)方法:窗函數(shù)法、頻率抽樣法、切比雪夫逼近法,找出設(shè)計(jì)方法存在的問題。并提出了一些解決方案。數(shù)字濾波器按照單位取樣響應(yīng)h(n)的時(shí)域特性可以分為無限脈沖響應(yīng)（IIR）系統(tǒng)和有限脈沖響應(yīng)（FIR）系統(tǒng)。FIR數(shù)字濾波器的優(yōu)點(diǎn)在于它可以做成具有嚴(yán)格線性相位,而同時(shí)可以具有任意的幅度特性;它的傳遞函數(shù)沒有極點(diǎn);這保證了設(shè)計(jì)出的F

8、IR數(shù)字濾波器是穩(wěn)定的。 目 錄 第一章 濾波器的簡(jiǎn)介 1 第二章 國(guó)內(nèi)外研究現(xiàn)狀和發(fā)展趨勢(shì) 3 2.1研究意義 3 2.2研究現(xiàn)狀 3 2.3發(fā)展趨勢(shì) 4 第三章 數(shù)字濾波器設(shè)計(jì)方法 5 3.1 數(shù)字濾波器設(shè)計(jì)的基本步驟 5 3.2 濾波器的MATLAB設(shè)計(jì) 6 第四章 FIR數(shù)字濾波器的MATLAB設(shè)計(jì) 7 4.1 窗函數(shù)設(shè)計(jì)FIR濾波器 7 4.2 頻率取樣法設(shè)計(jì)FIR濾波器 14 4.3 切比雪夫逼近法設(shè)計(jì)FIR濾波器 18 4.4 結(jié)論 20 結(jié)束語(yǔ) 21 致謝語(yǔ) 22 參考文獻(xiàn) 23

9、 第一章 濾波器的簡(jiǎn)介 數(shù)字濾波器是一種用來過濾時(shí)間離散信號(hào)的數(shù)字系統(tǒng)，通過對(duì)抽樣數(shù)據(jù)進(jìn)行數(shù)學(xué)處理來達(dá)到頻域?yàn)V波的目的?？梢栽O(shè)計(jì)系統(tǒng)的頻率響應(yīng)，讓它滿足一定的要求，從而對(duì)通過該系統(tǒng)的信號(hào)的某些特定的頻率成分進(jìn)行過濾，這就是濾波器的基本原理。如果系統(tǒng)是一個(gè)連續(xù)系統(tǒng)，則濾波器稱為模擬濾波器。如果系統(tǒng)是一個(gè)離散系統(tǒng)，則濾波器稱為數(shù)字濾波器。 數(shù)字濾波器可以理解為是一個(gè)計(jì)算程序或算法 ,將代表輸入信號(hào)的數(shù)字時(shí)間序列轉(zhuǎn)化為代表輸出信號(hào)的數(shù)字時(shí)間序列 ,并在轉(zhuǎn)化過程中 ,使信號(hào)按預(yù)定的形式變化。 數(shù)字濾波器按單位脈沖響應(yīng)h(n)的時(shí)域特性可分為

10、無限長(zhǎng)脈沖響IIR（Infinite Impulse Response）濾波器和有限長(zhǎng)脈沖響應(yīng)FIR（Finite Impulse Response）濾波器。 IIR濾波器一般采用遞歸型的實(shí)現(xiàn)結(jié)構(gòu)。其N階遞歸型數(shù)字濾波器的差分方程為 (1-1) 式中的系數(shù)至少有一項(xiàng)不為零。≠0 說明必須將延時(shí)的輸出序列反饋回來，也即遞歸系統(tǒng)必須有反饋環(huán)路。相應(yīng)的IIR濾波器的系統(tǒng)函數(shù)為 ：

11、 (1-2) FIR濾波器在保證幅度特性滿足技術(shù)要求的同時(shí), 比IIR 濾波器更容易做到有嚴(yán)格的線性相位特性。設(shè)FIR濾波器單位脈沖響應(yīng)h(n)長(zhǎng)度為N,其系統(tǒng)函數(shù)H(z)為: (1-3) FIR濾波器的設(shè)計(jì)任務(wù)是選擇有限長(zhǎng)度的h(n),使函數(shù)H(z)滿足技術(shù)要求。 FIR數(shù)字濾波器的特點(diǎn)(與IIR數(shù)字濾波器比較)： （1）很容易獲得嚴(yán)格的線性相位，避免被處理的信號(hào)產(chǎn)生相位失真，這一特點(diǎn)在寬頻帶信號(hào)處理、陣列信號(hào)處理、數(shù)據(jù)傳輸?shù)认到y(tǒng)中非常重

12、要； （2）可得到多帶幅頻特性； 信號(hào)通過線性系統(tǒng)后，其輸出就是輸入信號(hào)和系統(tǒng)沖激響應(yīng)的卷積。從頻域分析來看，信號(hào)通過線性系統(tǒng)后，輸出信號(hào)的頻譜將是輸入信號(hào)的頻譜與系統(tǒng)傳遞函數(shù)的乘積。線性系統(tǒng)對(duì)信號(hào)頻譜的這種加工，可以有意識(shí)地用來按照人們的需要改變信號(hào)的頻譜結(jié)構(gòu)。這就是所謂頻率濾波，是信號(hào)處理的一種方法。此時(shí)，人們也往往稱此系統(tǒng)為濾波器。由于計(jì)算機(jī)和大規(guī)模集成電路技術(shù)的進(jìn)步，依靠傳統(tǒng)的模擬電路來實(shí)現(xiàn)的電子系統(tǒng)已不適應(yīng)。現(xiàn)在都在開始采用數(shù)字化技術(shù)，傳統(tǒng)的模擬濾波器，正在被數(shù)字濾波器所代替，數(shù)字濾波器的輸入是個(gè)數(shù)字序列。從本質(zhì)上說它只是一個(gè)序列的運(yùn)算加工過程.但另方面因?yàn)樗?/p>

13、是一個(gè)離散系統(tǒng)，而一個(gè)離散系統(tǒng)具有一定的頻率響應(yīng)特性，適當(dāng)?shù)乜刂齐x散系統(tǒng)結(jié)構(gòu)使其頻率特性滿足一定的要求，就可以起到和模擬濾波器同樣的作用。但數(shù)字濾波器卻具有精度高，可靠性強(qiáng)，靈活性大，適應(yīng)范圍廣(在甚低頻范圍)，快速等優(yōu)點(diǎn)。而且可以分時(shí)復(fù)用，同時(shí)處理若干不同信號(hào)，因此已得到越來越廣泛的應(yīng)用。 第二章 國(guó)內(nèi)外研究現(xiàn)狀和發(fā)展趨勢(shì) 2.1研究意義 　　從事電子通信業(yè)而不能熟練操作使用MATLAB電子線路設(shè)計(jì)軟件，在工作和學(xué)習(xí)中將是寸步難行的。在數(shù)學(xué)、電子、金融等行業(yè)，使用MATLAB等計(jì)算機(jī)軟件對(duì)產(chǎn)品進(jìn)行設(shè)計(jì)、仿

14、真在很早以前就已經(jīng)成為了一種趨勢(shì)，這類軟件的問世也極大地提高了設(shè)計(jì)人員在通信、電子等行業(yè)的產(chǎn)品設(shè)計(jì)質(zhì)量與效率。眾所周知，實(shí)際過程中信號(hào)傳輸都要經(jīng)過調(diào)制與解調(diào)這一過程，由于消息傳過來的原始信號(hào)即調(diào)制信號(hào)具有頻譜較低的頻譜分量，這種信號(hào)在許多信道中不宜傳輸。因而，在通信系統(tǒng)的發(fā)送端通常需要有調(diào)制過程，反之在接收端則需要有解調(diào)過程。利用matlab的強(qiáng)大運(yùn)算功能，基于matlab信號(hào)處理工具箱（signal processing toolbox）的數(shù)字濾波器設(shè)計(jì)法可以快速有效的設(shè)計(jì)由軟件組成的常規(guī)數(shù)字濾波器，設(shè)計(jì)方便、快捷，極大的減輕了工作量。在設(shè)計(jì)過程中可以對(duì)比濾波器特性，隨時(shí)更改參數(shù)，以達(dá)到濾波

15、器設(shè)計(jì)的最優(yōu)化。利用matlab設(shè)計(jì)數(shù)字濾波器在電力系統(tǒng)二次信號(hào)處理軟件和微機(jī)保護(hù)中，有著廣泛的應(yīng)用前景。 2.2研究現(xiàn)狀 MATLAB是由MATH WORKS公司于1984年推出的一種面向科學(xué)與工程的計(jì)算軟件，通過MATLAB和相關(guān)工具箱，工程師、科研人員、數(shù)學(xué)家和教育工作者可以在統(tǒng)一的平臺(tái)下完成相應(yīng)的科學(xué)計(jì)算工作。MATLAB系統(tǒng)的強(qiáng)大功能是由其核心內(nèi)容（語(yǔ)言系統(tǒng)、開發(fā)環(huán)境、圖形系統(tǒng)、數(shù)學(xué)函數(shù)庫(kù)、應(yīng)用程序接口等）和輔助工具箱（符號(hào)計(jì)算、圖象處理、優(yōu)化、統(tǒng)計(jì)和控制等工具箱）兩大部分構(gòu)成。 MATLAB 本身包含了 600 余個(gè)用于數(shù)學(xué)計(jì)算、統(tǒng)計(jì)和工程處理的函數(shù)，這樣，就可以迅速完成科學(xué)

16、計(jì)算任務(wù)而不必進(jìn)行額外的開發(fā)。業(yè)內(nèi)領(lǐng)先的工具箱算法極大的擴(kuò)展了 MATLAB 的應(yīng)用領(lǐng)域，所以MATLAB自推出以來就受到廣泛的關(guān)注，信號(hào)處理工具箱就是其中之一，在信號(hào)處理工具箱中，MATLAB提供了濾波器分析、濾波器實(shí)現(xiàn)、FIR濾波器實(shí)現(xiàn)、IIR數(shù)字濾波器設(shè)計(jì)、IIR數(shù)字濾波器階次估計(jì)等方面的函數(shù)命令。 隨著信息時(shí)代和數(shù)字世界的到來，數(shù)字信號(hào)處理在通信、語(yǔ)音、圖像、自動(dòng)控制等眾多領(lǐng)域得到了廣泛的應(yīng)用。在數(shù)字信號(hào)處理應(yīng)用中，數(shù)字濾波器十分重要并已獲得廣泛應(yīng)用，并將得到進(jìn)一步的發(fā)展。 數(shù)字濾波器的實(shí)現(xiàn)方法一般有以下幾種: ①采用加法器、乘法器、延時(shí)器設(shè)計(jì)專用的濾波電路 ②用通用的可編程D

17、SP芯片實(shí)現(xiàn) ③采用計(jì)算機(jī)上的專用軟件進(jìn)行設(shè)計(jì)。 2.3發(fā)展趨勢(shì) 由于我們所面對(duì)的工程問題越來越復(fù)雜，過去所依賴分析的技術(shù)已逐漸不敷使用。利用電腦來分析及解決工程問題已是當(dāng)今工程師的必要工具。使用MATLAB軟件進(jìn)行科學(xué)計(jì)算，能夠極大加快科研人員進(jìn)行研究開發(fā)的進(jìn)度，減少在編寫程序和開發(fā)算法方面所消耗的時(shí)間和有限的經(jīng)費(fèi)，從而獲得最大的效能。隨著科學(xué)技術(shù)的發(fā)展，新技術(shù)應(yīng)用于濾波器上成為必然趨勢(shì)。利用流水線技術(shù)（PI技術(shù)），在基于多項(xiàng)式并行表示的并行結(jié)構(gòu)基礎(chǔ)上，構(gòu)造了新的FIR濾波器并行結(jié)構(gòu)。新結(jié)構(gòu)具有與原結(jié)構(gòu)相當(dāng)?shù)氖諗克俣?。但新結(jié)構(gòu)的復(fù)雜程度遠(yuǎn)遠(yuǎn)低于原結(jié)構(gòu)，與串行處理（單入單出）的濾波結(jié)構(gòu)相

18、當(dāng)。在保證較高的數(shù)據(jù)通過率和較低的系統(tǒng)功耗的同時(shí)，因?yàn)榻Y(jié)構(gòu)簡(jiǎn)單而更易于實(shí)現(xiàn),特別是在對(duì)現(xiàn)代信號(hào)處理要求高度集成的今天，對(duì)實(shí)踐具有一定的理論指導(dǎo)意義。 第三章 數(shù)字濾波器設(shè)計(jì)方法 所謂數(shù)字濾波器設(shè)計(jì)，簡(jiǎn)單地說，就是要找到一組能滿足特定濾波要求的系數(shù)向量 a 和 b。而濾波器設(shè)計(jì)完成后還需要進(jìn)一步考慮如何將其實(shí)現(xiàn)，即選擇什么樣的濾波器結(jié)構(gòu)來完成濾波運(yùn)算。FIR數(shù)字濾波器的設(shè)計(jì)方法很多，其中較為常用的是窗函數(shù)設(shè)計(jì)法、頻率采樣設(shè)計(jì)法和最優(yōu)化設(shè)計(jì)法。本文討論利用窗函數(shù)法、頻率采樣法和等波紋切比雪夫逼近法（調(diào)用 remez函數(shù)）來分別實(shí)現(xiàn)各種 FIR濾波器的設(shè)計(jì)。

19、FIR濾波器設(shè)計(jì)的任務(wù)是選擇有限長(zhǎng)度的h(n)，使傳輸函數(shù)H(e)滿足一定的幅度特性和線性相位要求。由于FIR濾波器很容易實(shí)現(xiàn)嚴(yán)格的線性相位，所以FIR數(shù)字濾波器設(shè)計(jì)的核心思 想是求出有限的脈沖響應(yīng)來逼近給定的頻率響應(yīng)。 目前在 MATLAB信號(hào)處理工具箱中，一共有 10種 FIR數(shù)字濾波器設(shè)計(jì)函數(shù)，對(duì)應(yīng) 5種設(shè)計(jì)方法，如表 3-1。 表3-1 FIR 濾波器設(shè)計(jì)函數(shù) 函數(shù) 設(shè)計(jì)方法 fir1,fir2,kaiserord 窗函數(shù)法 firls,remez,remezord 等波紋最小平方誤差設(shè)計(jì) fircls,fircls1 最小二乘約束設(shè)計(jì) cremez 任意頻響“

20、復(fù)濾波器”設(shè)計(jì) firrcos 升余弦設(shè)計(jì) 3.1 數(shù)字濾波器設(shè)計(jì)的基本步驟 數(shù)字濾波器的設(shè)計(jì)要經(jīng)過四步完成: 1. 確定指標(biāo) 在設(shè)計(jì)一個(gè)濾波器前，要根據(jù)應(yīng)用要求確定一些指標(biāo)，如以幅度指標(biāo)給定，一般有絕對(duì)指標(biāo)和相對(duì)指標(biāo)兩種方式。絕對(duì)指標(biāo)提供對(duì)幅度響應(yīng)││的要求，一般用于FIR數(shù)字濾波器。相對(duì)指標(biāo)，以分貝(dB)的形式提出要求，一般用于IIR 數(shù)字濾波器。 2 .逼近 確定技術(shù)指標(biāo)后，利用已有的數(shù)學(xué)和基本原理，提出一個(gè)濾波器的模型來逼近給定的指標(biāo)。 3. 實(shí)現(xiàn) 上面兩步得到的結(jié)果是系統(tǒng)函數(shù)H(z)或者單位沖激響應(yīng)h(n)描述的濾波器，根據(jù)此描述在計(jì)算機(jī)上用軟件實(shí)現(xiàn)。 4 .

21、軟件傳真 應(yīng)用MATLAB的傳真功能進(jìn)行算法級(jí)傳真，驗(yàn)證設(shè)計(jì)是否達(dá)到要求，如達(dá)到要求。則可以進(jìn)行后繼設(shè)計(jì); 否則要重新設(shè)計(jì)直到達(dá)到要求為止。 3.2 濾波器的MATLAB設(shè)計(jì) 數(shù)字濾波器是具有一定傳輸特性的數(shù)字信號(hào)處理裝置，其輸入與輸出都是離散的數(shù)字信號(hào)，借助于數(shù)字器件或一定的數(shù)值計(jì)算方法，對(duì)輸入信號(hào)進(jìn)行處理，改變輸入信號(hào)的波形或頻譜，達(dá)到保留信號(hào)中有用成分去除無用成分的目的。實(shí)際應(yīng)用中，數(shù)字濾波器往往作為模擬濾波器的一個(gè)組成部分，如圖所示。 A/D 數(shù)字濾波器 D/A x(n) y(n) Y(t) X(t) 圖3-1濾波系統(tǒng)框圖 輸入x(n)是一個(gè)時(shí)間序列，系

22、統(tǒng)函數(shù)為H(z)，其沖激響序列h(n)，則: y(n)=x(n)×h(n) 兩邊作傅立葉變換得頻域表達(dá)式:Y(jω)=H(j ω)X(jω)。H(jω)為數(shù)字濾波器的頻率特性，x(jω)和y(jω)分別為x(n)和y(n)的頻譜，數(shù)字濾波器的H(jω)改變了輸入x(n)的頻率特性。y(n)保留信號(hào)x(n)中有用頻率成分，去除無用的頻率成分，完成了數(shù)字濾波作用。 第四章 FIR數(shù)字濾波器的MATLAB設(shè)計(jì) 窗函數(shù)法設(shè)計(jì)的基本思想是把給定的頻率響應(yīng)通過IDTF逆離散傅里葉變換（Inverse Discrete Time Fourier Transform），求得脈

23、沖響應(yīng)，然后利用加窗函數(shù)對(duì)它進(jìn)行截?cái)嗪推交?，以?shí)現(xiàn)一個(gè)物理可實(shí)現(xiàn)且具有線性相位的FIR數(shù)字濾波器的設(shè)計(jì)目的。其核心是從給定的頻率特性，通過加窗確定有限長(zhǎng)單位取樣響應(yīng)h(n)；頻率采樣法設(shè)計(jì)的基本思想是把給出的理想頻率響應(yīng)進(jìn)行取樣，通過IDFT從頻譜樣點(diǎn)直接求得有限脈沖響應(yīng)；等波紋切比雪夫逼近法，利用MATLAB提供的remez函數(shù)實(shí)現(xiàn)Remez算法，設(shè)計(jì)濾波器逼近理想頻率響應(yīng)。 4.1 窗函數(shù)設(shè)計(jì)FIR濾波器 DTFT *ω(n) IDTFT 窗函數(shù)法設(shè)計(jì)FIR數(shù)字濾波器的最簡(jiǎn)單的方法是窗函數(shù)法，通常也稱之為傅立葉級(jí)數(shù)法。FIR數(shù)字濾波器的設(shè)計(jì)首先給出要求的理想濾波器的頻率響應(yīng),設(shè)

24、計(jì)一個(gè)FIR 數(shù)字濾波器頻率響應(yīng)，去逼近理想的濾波響應(yīng)。然而，窗函數(shù)法設(shè)計(jì)FIR數(shù)字濾波器是在時(shí)域進(jìn)行的，因而必須由理想的頻率響應(yīng)推導(dǎo)出對(duì)應(yīng)的單位取樣響應(yīng)hd(n)，再設(shè)計(jì)一個(gè)FIR數(shù)字濾波器的單位取樣響應(yīng)h(n)去逼近hd(n)。設(shè)計(jì)過程如下： hd(n) h(n) 加窗的作用是通過把理想濾波器的無限 長(zhǎng)脈沖響應(yīng)hd(n)乘以窗函數(shù)w(n)來產(chǎn)生一個(gè)被截?cái)嗟拿}沖響應(yīng)，即h(n)= hd(n)w(n)，并且對(duì)頻率響應(yīng)進(jìn)行平滑。MATLAB工具箱 提供的窗 函數(shù)有： 矩形窗 (Rectangular window)、三角窗(Tria

25、ngular window)、布拉克曼窗(Blackman window)、漢寧窗(Hanning window)、海明窗(Hamming window)、凱塞窗 (Kaiser window) 、巴特里特窗 (Bartlett window)、切比雪夫窗(Chebyshev window), 性能對(duì)比如表4-1. 表4-1 幾種窗函數(shù)的性能比較 窗函數(shù) 旁瓣峰值幅度（dB） 過度帶寬 （π/N） 阻帶最小衰減 （dB） 矩形窗 13 4 21 三角窗 25 8 25 漢寧窗 31 8 44 海明窗 41 8 53 凱塞窗（β=7.856）

26、57 10 80 布拉克曼窗 57 12 74 窗函數(shù)主要用來減少序列因截?cái)喽a(chǎn)生的Gibbs效應(yīng)。但當(dāng)這個(gè)窗函數(shù)為矩形時(shí)，得到的FIR濾波器幅頻響應(yīng)會(huì)有明顯的Gibbs效應(yīng)，并且任意增加窗函數(shù)的長(zhǎng)度（即 FIR濾波器的抽頭數(shù)）Gibbs效應(yīng)也不能得到改善。為了克服這種現(xiàn)象,窗函數(shù)應(yīng)該使設(shè)計(jì)的濾波器：（1）頻率特性的主瓣寬度應(yīng)盡量窄,且盡可能將能量集中在主瓣內(nèi)；(2)窗函數(shù)頻率特性的旁瓣ω趨于π的過程中，其能量迅速減小為零。 例：根據(jù)下列技術(shù)指標(biāo)，設(shè)計(jì)一個(gè)FIR低通濾波器。  通帶截止頻率ωp=0.2π，通帶允許波動(dòng)Ap=0.25dB；  阻帶截止頻率ω

27、s=0.3π，阻帶衰減As=50dB。  解：根據(jù)表4-1可知，海明窗和布拉克曼窗均可提供大于50dB的衰減。但海明窗具有較小的過渡帶從而具有較小的長(zhǎng)度N。 根據(jù)題意，所要設(shè)計(jì)的濾波器的過渡帶為 (4-1) 由表4-1可知，利用海明窗設(shè)計(jì)的濾波器的過渡帶寬Δω=8π/N，所以低通濾波器單位脈沖響應(yīng)的長(zhǎng)度為 (4-2) 3dB通帶截止頻率為

28、 (4-3) 理想低通濾波器的單位脈沖響應(yīng)為 （4-4） 海明窗為 （4-5） 則所設(shè)計(jì)的濾波器的單位脈沖響應(yīng)為

29、 （4-6） 所設(shè)計(jì)的濾波器的頻率響應(yīng)為 （4-7） FREQZ_M. M function [db, mag, pha, grd, w ] = freqz_m ( b, a) % db為相對(duì)振幅 (dB) ; mag為絕對(duì)振幅 ;pha為相位響應(yīng) ;grd為群延時(shí) ; W為頻率樣本點(diǎn)向量。 [ H , w ] = freqz( b, a, 1024 ,‘w hole’) ; % b

30、 和a 分別為H(z) 分子和分母多項(xiàng)式系數(shù) (對(duì)FIR 而言 , b = a, a = 1 ) . H=(H( 1: 501 )) ’; w=(w( 1: 501 )) ’; mag=abs(H) ; db= 203log10((mag+ep s)/max(mag)) pha=angle( H ) ; grd=grp delay(b,a,w) ; DEAL_L P.M Function hd = ideal_ lp(,M) %理想低通濾波器；hd為0到M-1之間的理想脈沖響應(yīng) ; 為截止頻率 ;M為理想濾波器長(zhǎng)度 alpha =(M-1)/2; n=0: M - 1;

31、 m=n-alpha+ eps; hd=sin(*m )./( pi*m ); ACTUAL_LP.m =0.10*pi; =0.20*pi; tr_width = - ; M=ceil(6.63pi/tr_width)+ 1; Disp([濾波器的長(zhǎng)度為：num2str(M)]) n=0:M-1; =(+)/2; =DEAL_LP(.M); w_han = ( hamming(M) ); h=hd * w_han; [db, mag, pha, grd, w ] = FREQZ_M (h,[ 1 ]) ; delta_w = 2*pi/1000; =-(min

32、(db(1:1:/delta_w+1))); Disp([最小組帶衰減為：num2str(As)]); subplot(221); stem ( n, ) ; title(‘理想沖激響應(yīng)’)； axis ( [ 0 M - 1 - 0. 1 0. 3 ] ) ; ylabel (‘ hd ( n) ’)； subplot ( 222 ) stem( n, w _han)； title（‘海明窗’） axis ( [ 0 M - 1 0 1. 1 ] ) ; ylabe l(‘w(n)’) subp lo t ( 223 ) stem ( n, h) ; title

33、(‘實(shí)際沖擊響應(yīng)’) axis ( [ 0 M - 1 - 0. 1 0. 3 ] ) ; xlabel (‘h’) ylabel(‘h(n)’) subplot ( 224 ) plot (w/pi, db) ; title(‘幅度響應(yīng)（db）’) axis ( [ 0 1 - 100 10 ] ) grid; xlabel(‘以pi為單位的頻率’) ylabel(‘分貝數(shù)’) 利用計(jì)算機(jī)編程實(shí)現(xiàn)，結(jié)果如圖所示。 圖4-1 理想低通濾波器的特性曲線 圖（a）是理想低通濾波器的單位脈沖響應(yīng)(n); 圖

34、(b)是海明窗函數(shù); 圖(c) 是實(shí)際低通濾波器的單位脈沖響應(yīng)h(n); 圖(d)是實(shí)際低通濾波器的幅頻特性|H()|，以dB為單位。濾波器長(zhǎng)N=80，實(shí)際阻帶衰減為As=53dB，通帶波動(dòng)為Ap=0.0316 dB，均滿足設(shè)計(jì)要求。 窗函數(shù)法還可以用 MATLAB信號(hào)處理工具箱中的 FIR濾波器設(shè)計(jì)函數(shù) fir1、fir2和kaiserord來實(shí)現(xiàn)。函數(shù) fir1 用來設(shè)計(jì)線性相位的低通、高通、帶通、帶阻 FIR濾波器（默認(rèn)為低通），它使用標(biāo)準(zhǔn)的窗函數(shù)法進(jìn)行設(shè)計(jì)（默認(rèn)窗函數(shù)為 Hamming），濾波器的階數(shù)由參數(shù) n指定（fir1在設(shè)計(jì)高通和帶阻濾波器時(shí)得到的濾波器的階數(shù)總是偶階），截止

35、頻率由參數(shù)Wn（歸一化截止頻率，對(duì)應(yīng)于頻響曲線上的-6dB點(diǎn)）定義；函數(shù) fir2用來設(shè)計(jì)多通帶任意響應(yīng) FIR 濾波器，該濾波器的幅頻特性由向量對(duì) f 和 m 確定，f 為歸一化頻率向量，m 為對(duì)應(yīng)頻率點(diǎn)上的幅度。 當(dāng)設(shè)計(jì)的濾波器在頻率為π的幅度響應(yīng)不是 0時(shí)，濾波器的階數(shù) n為偶數(shù)；當(dāng)函數(shù)kaiserord中各參數(shù)得到后可以利用返回值 n和beta設(shè)計(jì)凱塞窗函數(shù)，然后利用返回值和 ftype傳輸給 fir1進(jìn)行濾波器的設(shè)計(jì)。 三種函數(shù)比較完整的語(yǔ)法形式如下： b= fir1(n, , ‘ftype’, window)； b= fir2(n, f, m, npt, lap, wi

36、ndow)； [n, , beta, ftype] = kaiserord (f,a,dev, fs)； 下面分別用這三種函數(shù)設(shè)計(jì)： （1） 用 fir1 設(shè)計(jì)一個(gè) 27階 FIR低通濾波器， 其中通帶截止頻率為 0.2π，用漢寧窗函數(shù)。 n=27; =0.2; window=hanning(28); b=fir1(n, ,window); freqz(b); 幅度 相位 圖 4-2 fir1 設(shè)計(jì)的FIR低通濾波器頻率響應(yīng) （2） 用 fir2 設(shè)計(jì)一個(gè) 60階的 FIR濾波器，要求濾波器 0到π/4的幅度響應(yīng)為 0，π/4 到π/2的幅度響應(yīng)為 1/4，π/2到 3

37、π/4的幅度響應(yīng)為 0，3π/4 到 1的幅度響應(yīng)為 1。 n=60; f=[0 0.25 0.25 0.50 0.50 0.75 0.75 1]; m=[0 0 1/4 1/4 0 0 1 1]; %對(duì)幅頻響應(yīng)插值時(shí)插值點(diǎn)的個(gè)數(shù) npt=1024; %插值時(shí)不連續(xù)點(diǎn)轉(zhuǎn)變成連續(xù)時(shí)的點(diǎn)數(shù) lap=50; %衰減為 30dB的切比雪夫窗函數(shù) window=chebwin(61,30); b=fir2(n,f,m,npt,lap,window); freqz(b) 相位 幅度 圖 4-3 fir2 設(shè)計(jì)的FIR多通帶濾波器頻率響應(yīng) （3） 用 kaiserord 函

38、數(shù)設(shè)計(jì)一個(gè)帶通 FIR濾波器，通帶范圍是1500HZ到2500HZ，通帶的波紋最大為 0.03，阻帶范圍是0HZ到1000HZ和 3500HZ到 5000HZ，阻帶的波紋最大為 0.01，信號(hào)的采樣頻率為 10KHZ. f=[1000 1500 2500 3500]; a=[0 1 0]; fs=10000; dev=[0.01 0.03 0.01]; [n, ,beta,ftype]=kaiserord(f,a,dev, fs); hh=fir1(n, ,ftype,kaiser(n+1,beta),'noscale') freqz(hh); 查看 Workspace得到濾

39、波器的階數(shù)為45，截止頻率為 0.25π和 0.6π，凱塞窗函數(shù)的β值為 3.3953。 相位 幅度 圖 4-4 kaiserord 設(shè)計(jì)的 FIR 帶通濾波器頻率響應(yīng) 窗口法設(shè)計(jì)的主要優(yōu)點(diǎn)是簡(jiǎn)單，使用方便。窗口函數(shù)大多有封閉的公式可循，性能、參數(shù)都已有表格、資料可供參考，計(jì)算程序簡(jiǎn)便，所以很實(shí)用。缺點(diǎn)是通帶和阻帶的截止頻率不易控制。 4.2 頻率取樣法設(shè)計(jì)FIR濾波器 頻率采樣法是從頻域出發(fā)，把給定的理想頻率響應(yīng)（e）以等間隔采樣 （4-8） 以此Hd(

40、k)作為實(shí)際FIR數(shù)字濾波器的頻率響應(yīng)的采樣值H(k)，|=2πk/N k=0, 1, 2, …, N-1 （4-9） 知道H(k)后, 由IDFT定義, 可以用這N個(gè)采樣值H(k)來惟一確定有限長(zhǎng)序列h(n)，即 n=0, 1, 2, …, N-1 （4-10） 式中，h(n)為待設(shè)計(jì)的濾波器的單位脈沖響應(yīng)。其系統(tǒng)函數(shù)H(z)為 （4-1

41、1） 以上就是頻率采樣法設(shè)計(jì)濾波器的基本原理。此外，由頻域內(nèi)插公式知道，利用這N個(gè)頻域采樣值H(k)同樣可求得FIR濾波器的系統(tǒng)函數(shù)H(z) （4-12） (4-13) 其中，是一個(gè)內(nèi)插函數(shù)： （4-14） 從以上公式可以看出，在每個(gè)采樣頻率 點(diǎn)w k=2πk/N處，濾波器的實(shí)際頻率響應(yīng)是嚴(yán)格地和理想頻率響應(yīng)數(shù)值相等,即： （4-15） 而在各采樣點(diǎn)間的頻率響應(yīng)則是其的加權(quán)內(nèi)插函數(shù)

42、延伸疊加的結(jié)果。但對(duì)于一個(gè) 無限長(zhǎng)的序列，用頻率采樣法必然有一定的逼近誤差,誤差的大小取決于理想頻響曲線的形狀,理想頻響特性變換越平緩,則內(nèi)插函數(shù)值越接近理想值,誤差越小。為了提高逼近的質(zhì)量，可以通過在頻率相應(yīng)的過渡帶內(nèi) 插入比較連續(xù)的采樣點(diǎn)，擴(kuò)展過渡帶使其比較連續(xù)，從而使得通帶和阻帶之間變換比較緩慢，以達(dá)到減少逼近誤差的目的。 選取w∈[0,2π]內(nèi)N個(gè)采樣點(diǎn)的約束條件為： (4-16) 下面用頻率采樣法設(shè)計(jì)一個(gè)FIR通帶截止頻率為 0.2π，通帶波紋最大為 0.04，阻帶截止頻率為 0.3π，阻帶波紋最大為 0.02，濾波器的階數(shù)通過remezord函數(shù)估算，程序如下： rem

43、ezord函數(shù)用于估算FIR數(shù)字濾波器的等波紋最佳一致逼近設(shè)計(jì)的最低階數(shù)N，從而使 濾波器在滿足指標(biāo)的前提下造價(jià)最低。 f=[0.2 0.3]; a=[1 0]; dev=[0.04 0.02]; %給出濾波器的參數(shù) [n f0 a0 w]=remezord(f,a,dev); N=n; alpha=(N-1)/2; k=0:N-1; =0.2*pi; =0.3*pi; %計(jì)算理想低通濾波器的截止頻率 =( +)/2; m=fix(*N/(2*pi)+1); %在兩邊過渡帶取值為 0.5 的采樣點(diǎn) T = 0.5; Hrs=[ones(1,m),T,zeros(1,N

44、-2*m-1),T,ones(1,m-1)] k1 = 0:floor(alpha); k2 = floor(alpha+1):N-1; phai=[-alpha*(2*pi)/N*k1,alpha*(2*pi)/N*(N-k2)]; H =Hrs.*exp(j*phai); %計(jì)算單位沖激響應(yīng) h =ifft(H,N); [h1,w1] = freqz(h,1,256,1); hr = abs( h1); h1 = 20* log10(hr); %畫出 FIR DF的單位取樣響應(yīng) figure(1); k=0:N-1; stem(k,h,'k.') axis([0

45、,N-1,1.1*min(real(h)),1.1*max(real(h))]); xlabel('n'); ylabel('h(n)'); grid on; %畫出 FIR DF 的低通衰減幅頻特性 figure(2); plot(w1,h1); xlabel('Normalized Frequency(×πrad/sample)'); ylabel('Magnitude(dB)'); grid on; 由 Workspace 查出濾波器參數(shù)為 n=27,f0=[0 0.2 0.31], a0=[1 1 0 0 ],w=[1 2 ]。在理想濾波器的設(shè)計(jì)中，若不增加過渡點(diǎn)，阻帶和

46、通帶之間的衰減約為21dB。如果在通帶和阻帶之間增加一個(gè)采樣點(diǎn)，阻帶的最小衰減可以達(dá)到65dB（衰減程度由選擇的采樣點(diǎn)決定），本程序在過渡帶取值為 0.5 的采樣點(diǎn)，阻帶衰減約為 40dB. 圖 4-5 FIR 的單位取樣響應(yīng) 圖 4-6 FIR 的低通衰減幅頻特性 頻率采樣法的優(yōu)點(diǎn)是可以在頻域直接設(shè)計(jì)，并且適合最優(yōu)化設(shè)計(jì); 缺點(diǎn)是采樣頻率只能等于2π/N的整數(shù)倍，因而不能確保截止頻率ωc的自由取值，要想實(shí)現(xiàn)自由地選擇截止頻率，必須增加采樣點(diǎn)數(shù)N，但這又使計(jì)算量加大。 4.3 切比雪夫逼近法設(shè)計(jì)FIR濾波器 盡管窗函數(shù)法與頻率采樣法在FIR數(shù)字濾波器的設(shè)計(jì)中有著廣泛

47、的應(yīng)用, 但兩者都不是最優(yōu)化的設(shè)計(jì)。通常線性相位濾波器在不同的頻帶內(nèi)逼近的最大容許誤差要求不同。等波紋切比雪夫逼近準(zhǔn)則就是通過對(duì)通帶和阻帶使用不同的加權(quán)函數(shù),實(shí)現(xiàn)在不同頻段(通常指的是通帶和阻帶)的加權(quán)誤差最大值相同,從而實(shí)現(xiàn)其最大誤差在滿足性能指標(biāo)的條件下達(dá)到最小值，即使得和之間的最大絕對(duì)誤差最小。 等波紋切比雪夫逼近是采用加權(quán)逼近誤差 E它可以表示為： E=W（-H） (4-17) 其中，W為逼近誤差加權(quán)函數(shù)，在誤差要求高的頻段上，可以取較大的加權(quán)值，否則，應(yīng)當(dāng)取較小的加權(quán)值。 盡管按照 FIR數(shù)字濾波器單位取樣響應(yīng) h(n)的對(duì)稱性和

48、 N的奇、偶性，F(xiàn)IR 數(shù)字濾波器可以分為 4種類型，但濾波器的頻率響應(yīng)可以寫成統(tǒng)一的形式： （4-18） 其中，k∈{0 ,1} , H (ω)為幅度函數(shù)， 且是一個(gè)純實(shí)數(shù)，表達(dá)式也可以寫成統(tǒng)一的 形式： H=Q（ω）P（ω） (4-19) 其中，Q（ω）為ω的固定函數(shù)，P（ω）為M個(gè)余弦函數(shù)的線性組合。若令：W(ω)= WQ（ω），Hd(ω)= / Q（ω），因此，由式（4-19）、（4-20）將E改寫成： E=W(ω)[ Hd(ω)- P（ω）] （4-20）

49、 故等波紋切比雪夫逼近法設(shè)計(jì)FIR數(shù)字濾波器的步驟是： ①給出所需的頻率響應(yīng)，加權(quán)函數(shù)W和濾波器的單位取樣響應(yīng)h(n)的長(zhǎng)度N。 ②由①中給定的參數(shù)來形成所需的W(ω)、Hd(ω)和P（ω）的表達(dá)式。 ③根據(jù)Remez 算法,求解逼近問題。 ④利用傅立葉逆變換計(jì)算出單位取樣響應(yīng)h( n) 。 Remez 算法是由Parks 和McClellan 等人在1972年推導(dǎo)出來的。它是將FIR數(shù)字濾波器中的五個(gè)參數(shù)( N ,σ1 , σ2 ,ωp ,ωs) 中的N ,ωp , ωs 和σ1/σ2 固定,而視σ1(或σ2)為變量的一種迭代方法。在MATLAB工具箱中可以直接調(diào)用re

50、mez函數(shù)（采用Remez算法），來進(jìn)行FIR數(shù)字濾波器的設(shè)計(jì)。其具體算法有幾種，常見的一種算法格式為： b= remez (n, f, a, w,‘ftype’); 下面用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); 相位 幅度 圖4-7 remez設(shè)計(jì)的FIR低通濾波器頻率響應(yīng)

51、4.4 結(jié)論 比較窗函數(shù)法的程序I、頻率采樣法和等波紋切比雪夫逼近法中的程序，同樣是設(shè)計(jì)一個(gè)27階的FIR低通數(shù)字濾波器，可以看出： （1）窗函數(shù)法在階數(shù)較低時(shí)，阻帶特性不滿足設(shè)計(jì)要求，只有當(dāng)濾波器階數(shù)較高時(shí)，使用海明窗和凱塞窗基本可以達(dá)到阻帶衰耗要求； （2）頻率采樣法偏離設(shè)計(jì)指標(biāo)最明顯，阻帶衰減最小，而且設(shè)計(jì)比采用窗函數(shù)法復(fù)雜。只有適當(dāng)選取過渡帶樣點(diǎn)值，才會(huì)取得較好的衰耗特性； （3）利用等波紋切比雪夫逼近法則的設(shè)計(jì)可以獲得最佳的頻率特性和衰耗特性，具有通帶和阻帶平坦，過渡帶窄等優(yōu)點(diǎn)。 綜上所述，F(xiàn)IR濾波器很容易實(shí)現(xiàn)具有嚴(yán)格線性相位的系統(tǒng), 使信號(hào)經(jīng)過處理后不產(chǎn)生相位失

52、真,舍入誤差小,而且穩(wěn)定，因此越來越受到廣泛的重視。 結(jié)束語(yǔ) 與 C語(yǔ)言的程序設(shè)計(jì)相比較 , 利用 MATLAB可節(jié)省大量的編程時(shí)間。將其用于數(shù)字濾波器的編程 ,提高了編程效率。隨著版本的不斷提高 ,相信 MATLAB在數(shù)字濾波技術(shù)中會(huì)發(fā)揮更大的作用。FIR數(shù)字濾波器在數(shù)字處理領(lǐng)域有非常重要的地位, 因?yàn)樵谠S多信號(hào)處理系統(tǒng), 如圖像處理、數(shù)字音頻等系統(tǒng)要求具有線性相位特性。在這方面FIR 濾波器有獨(dú)特的優(yōu)點(diǎn)。此外, FIR濾波器的沖激響應(yīng)是有限長(zhǎng)序列, 其系統(tǒng)函數(shù)為一個(gè)多項(xiàng)式,它所含的極點(diǎn)多為原點(diǎn),

53、 所以FIR濾波器永遠(yuǎn)是穩(wěn)定的。在應(yīng)用Matlab語(yǔ)言進(jìn)行FIR濾波器的設(shè)計(jì)時(shí), 可以隨時(shí)對(duì)比設(shè)計(jì)要求和濾波器特性, 并可通過不斷調(diào)整設(shè)計(jì)參數(shù), 獲得較合適的沖激響應(yīng)和幅度響應(yīng), 以使濾波器達(dá)到最優(yōu)化。 致謝語(yǔ) 本設(shè)計(jì)及畢業(yè)論文是在老師的親切關(guān)懷和悉心指導(dǎo)下完成的。他嚴(yán)肅的科學(xué)態(tài)度，嚴(yán)謹(jǐn)?shù)闹螌W(xué)精神，精益求精的工作作風(fēng)，深深地感染和激勵(lì)著我。從課題的選擇到項(xiàng)目的最終完成，老師都始終給予我細(xì)心的指導(dǎo)和不懈的支持。在此謹(jǐn)向老師致以誠(chéng)摯的謝意和崇高的敬意。 在此，我還要感謝在一起愉快的度過本科兩年的宿舍的各位姐妹們，

54、正是由于你們的幫助和支持，我才能克服一個(gè)一個(gè)的困難和疑惑，直至本文的順利完成此外，從開始進(jìn)入課題到論文的順利完成，有多少可敬的師長(zhǎng)、同學(xué)、朋友給了我無言的幫助，在這里請(qǐng)接受我誠(chéng)摯的謝意！最后我還要感謝培養(yǎng)我長(zhǎng)大含辛茹苦的父母，謝謝你們！ 參考文獻(xiàn) [1] 張亞妮. 基于MATLAB的數(shù)字濾波器設(shè)計(jì)[J]. 遼寧工程技術(shù)大學(xué)學(xué)報(bào)， 2005，(05):5～15 [2] 饒志強(qiáng)，葉念渝. FIR和IIR數(shù)字濾波器的探討與實(shí)現(xiàn)[J].計(jì)算機(jī)與數(shù)字工程，2005，(06):1～6 [3] 孫強(qiáng). 運(yùn)用MATLAB實(shí)現(xiàn)數(shù)

55、字濾波器的設(shè)計(jì)[J].電腦學(xué)習(xí)，2005，(02):2～22 [4] 何小河. 基于MATLAB的FIR數(shù)字濾波器設(shè)計(jì)[J]. 科技資訊， 2008，(03):15～19 [5] 林彥杰. 基于MATLAB的FIR數(shù)字濾波器設(shè)計(jì)與仿真[J]. 寧夏工程技術(shù)，2004，(01):7～10 [6] 張登奇， 李望移. 基于MATLAB的FIR數(shù)字濾波器設(shè)計(jì)[J]. 計(jì)算機(jī)時(shí)代，2007(11):9～21 [7] 趙紅怡.數(shù)字信號(hào)處理及其MATLAB實(shí)現(xiàn)[M].北京化學(xué)工業(yè)出版社，工業(yè)裝備與信息工程出版中心，2007 [8] 程佩青. 數(shù)字信號(hào)處理教程[M].北京清華大學(xué)出版社，2003 26