基于MATLAB的FIR數(shù)字濾波器設(shè)計論文

1、摘要摘 要FIR數(shù)字濾波器是數(shù)字信號處理的一個重要組成部分，由于FIR數(shù)字濾波器具有嚴(yán)格的線性相位，因此在信息的采集和處理過程中得到了廣泛的應(yīng)用。本文介紹了FIR數(shù)字濾波器的概念和線性相位的條件，分析了窗函數(shù)法、頻率采樣法和等波紋逼近法設(shè)計FIR濾波器的思路和流程。在分析三種設(shè)計方法原理的基礎(chǔ)上，借助Matlab仿真軟件工具箱中的fir1、fir2和remez子函數(shù)分別實現(xiàn)窗函數(shù)法、頻率采樣法和等波紋逼近法設(shè)計FIR濾波器。然后檢驗濾波器的濾波效果，采用一段音頻進行加噪聲然后用濾波器濾，對比三段音頻效果進而對濾波器的濾波效果進行檢驗。仿真結(jié)果表明，在相頻特性上，三種方法設(shè)計的FIR濾波器在通帶

2、內(nèi)都具有線性相位；在幅頻特性上，相比窗函數(shù)法和頻率采樣法，等波紋逼近法設(shè)計FIR濾波器的邊界頻率精確，通帶和阻帶衰減控制。關(guān)鍵詞：FIR數(shù)字濾波器；窗函數(shù)法；頻率采樣法；等波紋逼近法；加噪聲；MatlabABSTRACTIn 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 digit

3、al filter 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 compare

4、ding the 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 requireme

5、nts It is 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

6、filtered signal are contrasted and analied.Key Words:FIR digital filters，Window function，F(xiàn)requency sampling，Equiripple approximation method，Matlab目 錄第1章緒論11.1數(shù)字濾波器的研究背景與意義1,511.2數(shù)字濾波器的應(yīng)用現(xiàn)狀與發(fā)展趨勢1,511.3數(shù)字濾波器的實現(xiàn)方法分析1，51第2章 數(shù)字濾波器的簡介12.1FIR數(shù)字濾波器的基本概念1(2.1 2.3) 2(2.1 2.3)12.2FIR數(shù)字濾波器的特點512.3FIR數(shù)字濾波器種類512.

7、4FIR數(shù)字濾波器線性相位條件31第3章FIR數(shù)字濾波器的設(shè)計33.1 窗函數(shù)法2,3，533.2 頻率采樣法2,3，533.3 最優(yōu)化設(shè)計3,533.4 三種設(shè)計法之間的對比1第4章基于Matlab的FIR數(shù)字濾波器設(shè)計44.1Matlab簡介2,514.1.1Matlab簡介2,514.1.2 本文用到函數(shù)的說明14.2窗函數(shù)法的Matlab實現(xiàn)2，314.2.1各種窗函數(shù)之間的對比414.3頻率采樣法的Matlab實現(xiàn)2,314.4等波紋逼近法的Matlab實現(xiàn)2,31第5章濾波器使用的實例說明55.1語音信號的采集155.2語音信號加噪與頻譜分析155.3設(shè)計巴特沃斯低通濾波器165.

8、4用濾波器對加噪語音濾波165.5比較濾波前后語音信號波形及頻譜19第6章結(jié)論4參考文獻10致謝12附錄13附錄一：窗函數(shù)程序13附錄二：抽樣函數(shù)程序13附錄三：最優(yōu)化設(shè)計程序13附錄四：實例說明程序13外文資料原文14譯文15第1章 緒論1.1數(shù)字濾波器的研究背景和意義 當(dāng)今，數(shù)字信號處理技術(shù)正飛速發(fā)展，它不但自成一門學(xué)科，更是以不同形式影響和滲透到其他學(xué)科；它與國民經(jīng)濟息息相關(guān)，與國防建設(shè)緊密相連；它影響或改變著我們的生產(chǎn)、生活方式，因此受到人們普遍的關(guān)注。 數(shù)字化、智能化和網(wǎng)絡(luò)化是當(dāng)代信息技術(shù)發(fā)展的大趨勢，而數(shù)字化是智能化和網(wǎng)絡(luò)化的基礎(chǔ)，實際生活中遇到的信號多種多樣，例如廣播信號、電視信

9、號、雷達信號、通信信號、導(dǎo)航信號等等。上述這些信號大部分是模擬信號，也有小部分是數(shù)字信號。模擬信號是自變量的連續(xù)函數(shù)，自變量可以是一維的，也可以是二維或多維的。大多數(shù)情況下一維模擬信號的自變量是時間，經(jīng)過時間上的離散化（采樣）和幅度上的離散化（量化），這類模擬信號便成為一維數(shù)字信號。因此，數(shù)字信號實際上是數(shù)字序列表示的信號，語音信號經(jīng)采樣和量化后，得到的數(shù)字信號是一個一維離散時間序列；而圖像信號的經(jīng)采樣和量化后，得到的是數(shù)字信號是一個二維離散空間序列。數(shù)字濾波技術(shù)是數(shù)字信號分析、處理技術(shù)的重要分支。無論是信號的獲取、傳輸，還是信號的處理和交換都離不開濾波技術(shù)，它對信號安全可靠和有效靈活地傳輸時

10、至關(guān)重要的。在所有的電子系統(tǒng)中，使用最多技術(shù)最復(fù)雜的要算數(shù)字濾波器了。數(shù)字濾波器的優(yōu)劣直接決定產(chǎn)品的優(yōu)劣。1.2數(shù)字濾波器的應(yīng)用現(xiàn)狀與發(fā)展趨勢數(shù)字濾波器精確度高、使用靈活、可靠性高，具有模擬設(shè)備所沒有的許多優(yōu)點，已廣泛應(yīng)用于各個學(xué)科技術(shù)領(lǐng)域，例如數(shù)字電視、語音、通信、雷達、聲納、遙感、圖像、生物醫(yī)學(xué)以及許多工程應(yīng)用領(lǐng)域。隨著信息時代數(shù)字時代的到來，數(shù)字濾波技術(shù)已經(jīng)成為一門極其重要的學(xué)科和技術(shù)領(lǐng)域。以往的濾波器大多采用模擬電路技術(shù)，但是，模擬電路技術(shù)存在很多難以解決的問題，例如，模擬電路元件對溫度的敏感性，等等。而采用數(shù)字技術(shù)則避免很多類似的難題，當(dāng)然數(shù)字濾波器在其他方面也有很多突出的優(yōu)點，所以

11、采用數(shù)字濾波器對信號進行處理是目前的發(fā)展方向。在信號處理過程中，所處理的信號往往混有噪音，從接收到的信號中消除或減弱噪音是信號傳輸和處理中十分重要的問題。根據(jù)有用信號和噪音的不同特性，提取有用信號的過程稱為濾波，實現(xiàn)濾波功能的系統(tǒng)稱為濾波器。在近代電信設(shè)備和各類控制系統(tǒng)中，數(shù)字濾波器應(yīng)用極為廣泛。數(shù)字濾波器的實現(xiàn)方法分析數(shù)字濾波器的實現(xiàn)，大體上有如下幾種方法：(1) 在通用的微型機上用軟件來實現(xiàn)。軟件可以由使用者自己編寫或使用現(xiàn)成的。自IEEE DSP Comm.于1979年推出第一個信號處理軟件包以來，國外的研究機構(gòu)、公司也陸續(xù)推出不同語言不同用途的信號處理軟件包。這種實現(xiàn)方法速度較慢，多用

12、于教學(xué)與科研。(2) 用單片機來實現(xiàn)。目前單片機的發(fā)展速度很快，功能也很強依靠單片機的硬件環(huán)境和信號處理軟件可用于工程實際，如數(shù)字控制、醫(yī)療儀器等。(3) 利用專門用于信號處理的DSP片來實現(xiàn)。DSP芯片較之單片機有著更為突出的優(yōu)點，如內(nèi)部帶有乘法器、累加器，采用流水線工作方式及并行結(jié)構(gòu)，多總線，速度快，配有適于信號處理的指令等，DSP芯片的問世及飛速發(fā)展，為信號處理技術(shù)應(yīng)用于工程實際提供了可能。第2章 數(shù)字濾波器的簡介2.1 FIR數(shù)字濾波器的基本概念FIR濾波器：有限長單位沖激響應(yīng)濾波器，是數(shù)字信號處理系統(tǒng)中最基本的元件，它可以在保證任意幅頻特性的同時具有嚴(yán)格的線性相頻特性，同時其單位抽樣

13、響應(yīng)是有限長的，因而濾波器是穩(wěn)定的系統(tǒng)。數(shù)字濾波器是一個離散的系統(tǒng)。它可以對輸入的離散信號進行一系列運算處理，從輸入的信號中獲得所需要的信息。數(shù)字濾波器的系統(tǒng)函數(shù)通常表示為 數(shù)字濾波器分為有限沖激響應(yīng)數(shù)字濾波器，即FIR數(shù)字濾波器和無限沖激響應(yīng)，即IIR數(shù)字濾波器。從公式的角度來看，F(xiàn)IR數(shù)字濾波器的 始終為零；IIR數(shù)字濾波器至少有一個非零。實現(xiàn)數(shù)字濾波器的方法一般有兩種：一種方法是吧濾波器所要完成的運算編成程序并讓計算機執(zhí)行，也就是采用計算機軟件來實現(xiàn)；另一種方法是設(shè)計專用的數(shù)字硬件、專用的數(shù)字信號處理器或采用通用的數(shù)字信號處理器來實現(xiàn)。實現(xiàn)一個數(shù)字濾波器一般需要三個基本的運算單元：加法器

14、、單位延時器和常數(shù)乘法器。這些基本的單元可以有兩種表示方法方框圖法和信號流程圖法。不論軟件實現(xiàn)還是硬件實現(xiàn)，在濾波器設(shè)計過程中，由同一系統(tǒng)函數(shù)可以構(gòu)成很多不同的運算結(jié)構(gòu)。對于無限精度的系數(shù)和變量，不同結(jié)構(gòu)可能是等效的，與其輸入和輸出特性無關(guān)；但是在系數(shù)和變量精度是有限的情況下，不同運算結(jié)構(gòu)的性能就有很大的差異。因此，有必要對離散時間系統(tǒng)的結(jié)構(gòu)有一基本認(rèn)識。FIR濾波器是指在有限范圍內(nèi)系統(tǒng)的單位脈沖響應(yīng)hk僅有非零值的濾波器。M階FIR濾波器的系統(tǒng)函數(shù)H(z)為 其中H(z)是的M階多項式，在有限的z平面內(nèi)H(z)有M個零點，在z平面原點z=0有M個極點. FIR濾波器的頻率響應(yīng)為 它的另外一種

15、表示方法為 2.2 FIR濾波器的特點有限長單位沖激響應(yīng)（FIR）濾波器有以下特點： （1） 系統(tǒng)的單位沖激響應(yīng)在有限個n值處不為零； （2） 系統(tǒng)函數(shù)在處收斂，極點全部在z = 0處（因果系統(tǒng)）； （3） 結(jié)構(gòu)上主要是非遞歸結(jié)構(gòu)，沒有輸出到輸入的反饋，但有些結(jié)構(gòu)中（例如頻率抽樣結(jié)構(gòu)）也包含有反饋的遞歸部分。 設(shè)FIR濾波器的單位沖激響應(yīng)為一個N點序列，則濾波器的系統(tǒng)函數(shù)為 就是說，它有（N1）階極點在z = 0處，有（N1）個零點位于有限z平面的任何位置。優(yōu)點 ：（1）很容易獲得嚴(yán)格的線性相位，避免被處理的信號產(chǎn)生相位失真，這一特點在寬頻帶信號處理、陣列信號處理、數(shù)據(jù)傳輸?shù)认到y(tǒng)中非常重要；（

16、2）可得到多帶幅頻特性；（3）極點全部在原點（永遠(yuǎn)穩(wěn)定），無穩(wěn)定性問題；（4）任何一個非因果的有限長序列，總可以通過一定的延時，轉(zhuǎn)變?yōu)橐蚬蛄?，所以因果性總是滿足；（5）無反饋運算，運算誤差小。缺點：（1）因為無極點，要獲得好的過渡帶特性，需以較高的階數(shù)為代價；（2）無法利用模擬濾波器的設(shè)計結(jié)果，一般無解析設(shè)計公式，要借助計算機輔助設(shè)計程序完成。2.3 FIR濾波器的種類（1）數(shù)字集成電路FIR濾波器一種是使用單片通用數(shù)字濾波器集成電路，這種電路使用簡單，但是由于字長和階數(shù)的規(guī)格較少，不易完全滿足實際需要。雖然可采用多片擴展來滿足要求，但會增加體積和功耗，因而在實際應(yīng)用中受到限制。設(shè)計數(shù)字濾波

17、器的任務(wù)就是尋求一個因果穩(wěn)定的線性時不變系統(tǒng)，使其系統(tǒng)函數(shù)H(z)具有指定的頻率特性。（2）DSP芯片F(xiàn)IR濾波器另一種是使用DSP芯片。DSP芯片有專用的數(shù)字信號處理函數(shù)可調(diào)用，實現(xiàn)FIR濾波器相對簡單，但是由于程序順序執(zhí)行，速度受到限制。而且，就是同一公司的不同系統(tǒng)的DSP芯片，其編程指令也會有所不同，開發(fā)周期較長。 .y(n)(a)y(n)(b)圖1 DSP FIR濾波器（3）可編程FIR濾波器還有一種是使用可編程邏輯器件，F(xiàn)PGACPLD。FPGA有著規(guī)整的內(nèi)部邏輯塊整列和豐富的連線資源，特別適合用于細(xì)粒度和高并行度結(jié)構(gòu)的FIR濾波器的實現(xiàn)，相對于串行運算主導(dǎo)的通用DSP芯片來說，并行

18、性和可擴展性都更好。2.4 FIR數(shù)字濾波器線性相位條件設(shè)FIR數(shù)字濾波器脈沖響應(yīng)的長度為N，則其頻率響應(yīng)可以表示為： （2-4）上式通過歐拉恒等式展開可得到的相位特性，有兩種線性相位特性，通常稱為第一類線性相位和第二類線性相位。第一類線性相位特性： 是一個與無關(guān)的常數(shù)第二類線性相位特性： 是起始相位嚴(yán)格地說第二種情況時的是不具有線性相位特性的，但上述兩種情況都滿足群延遲是一個常數(shù)，仍可以視為具有線性相位的，在第二類線性相位中是常用的一種情況。對于第一類線性相位，即，通過一系列的運算整理之后可得到一個三角函數(shù)求和公式： （2-5）式中正弦函數(shù)為奇對稱，當(dāng)時，對稱中心為，需要滿足關(guān)于偶對稱，即要

19、求：， （2-6）對于第二類線性相位，即時，通過運算得到公式： （2-7）函數(shù)為偶對稱，當(dāng)時，對稱中心也為。若要使上式成立，則要使關(guān)于奇對稱，即要求：， （2-8）從上述分析看來，線性相位FIR數(shù)字濾波器的時域約束條件是指滿足線性相位時對的約束條件，對于第一類線性相位，沖激響應(yīng)h(n)滿足（2-6）式；對于第二類線性相位，沖激響應(yīng)h(n)滿足（2-8）式。第3章FIR數(shù)字濾波器的設(shè)計3.1 FIR數(shù)字濾波器的窗函數(shù)設(shè)計法窗函數(shù)設(shè)計法是FIR數(shù)字濾波器里最簡單的一種設(shè)計法，又叫傅里葉級數(shù)法，為了設(shè)計簡單方便，通常選擇所希望逼近的濾波器的頻率響應(yīng)函數(shù)為具有片段常數(shù)特性的理想濾波器，尋找一組，確定其

20、頻率響應(yīng)，然后用來逼近1。窗函數(shù)法設(shè)計FIR濾波器是在時域中進行的，那么可以通過傅里葉反變換得到得到頻率響應(yīng)，即： （3-1）在實際中，一般是處于逐段恒定的，在邊界頻率處有不連續(xù)點，因而單位脈沖響應(yīng)是無限長的非因果序列，不能直接作為FIR數(shù)字濾波器的單位脈沖響應(yīng)，因此需要對進行階段，轉(zhuǎn)換為有限長的一段因果序列，也就是用一個有限長度的窗函數(shù)序列來截取，即，并將非因果序列轉(zhuǎn)變?yōu)橐粋€因果序列。截取的長度和加權(quán)窗函數(shù)的形狀都直接影響到逼近精度。窗函數(shù)法設(shè)計FIR濾波器過程如圖3-1所示： 圖3-1 窗函數(shù)法設(shè)計FIR濾波器流程以截止頻率為，相位為零的理想低通濾波器為例，其頻率特性為： （3-2）通過傅

21、里葉反變換得到對應(yīng)的為： （3-3）此時的是一個無限長的非因果序列，我們需要對其進行截斷，變成一個有限長的因果序列?？梢韵劝严蛴移揭苽€點，得到為： （3-4）相應(yīng)的傳輸函數(shù)為： （3-5）然后對截取從0到的N個點，N為窗函數(shù)的長度，所得的結(jié)果表示為： （3-6）表示窗函數(shù)，一般用下標(biāo)來表示窗函數(shù)的類型，矩形窗記為。3.1.2 吉布斯效應(yīng)用窗函數(shù)對進行直接截斷，得到有限長序列，并以代替，肯定會引起誤差，表現(xiàn)在頻域就是通常所說的吉布斯(Gibbs)效應(yīng)。對于一個在有限區(qū)間分布的信號，其連續(xù)頻譜在頻域上分布往往是無限的，而在實際信號處理時，我們通常只能在有限的區(qū)間內(nèi)做傅里葉分析，也就是說，我們只能用

22、有限區(qū)間來代替理論分析中的無限區(qū)間，多數(shù)情況下，我們總是選擇信號的低頻部分，而舍棄高頻部分。而信號的高頻部分往往是反應(yīng)信號快速變化特征，如果信號本身是連續(xù)的，這樣做一般不會引起信號的顯著變化，但實際中的信號往往是比較豐富的，特別是信號本身存在劇烈變化的，這樣做必定會引起一些誤差。該誤差引起過渡帶加寬以及通帶和阻帶內(nèi)的波動8。為了減小吉布斯效應(yīng)帶來的影響，需要調(diào)整窗口的長度來控制過渡帶的寬度，但要減小帶內(nèi)的波動以及增大阻帶衰減，還需要從窗函數(shù)的形狀上尋找解決方法。為了減少序列因截斷而產(chǎn)生的Gibbs效應(yīng)，窗函數(shù)在設(shè)計時需要注意：(1)頻率特性的主瓣要盡可能的窄，并且盡量把能量都集在主瓣內(nèi)。(2)

23、窗函數(shù)頻率特性的旁瓣趨于的過程中，其能量迅速減小為零。3.1.3 常見窗函數(shù)介紹常見的窗函數(shù)有矩形窗（Rectangle Window）、三角形窗（Bartlerr Window）、漢寧（Hanning）窗升余弦窗、哈明（Hamming）窗改進的升余弦窗、布萊克曼（Blackman）窗、凱塞貝塞爾窗（Kaiser-Basel Window）。矩形窗的窗函數(shù)為： （3-7）其頻譜的幅度函數(shù)為 （3-8）矩形窗的主瓣寬度為，用矩形窗設(shè)計的FIR數(shù)字濾波器的過渡帶寬度近似為。三角形窗的窗函數(shù)為： （3-9）其頻譜的幅度函數(shù)為 （3-10）三角窗的主瓣寬度為。漢寧窗窗函數(shù)為 （3-11）漢寧窗的頻譜的

24、幅度函數(shù)為 （3-12）漢寧窗的主瓣寬度為，漢寧窗在其兩個端點都為零，實際中這兩個端點的數(shù)據(jù)是不可用的。哈明窗的窗函數(shù)為 （3-13）其幅度函數(shù)為 （3-14）哈明窗是一種改進的余弦窗，能量更加集中在主瓣，是一種高效的窗函數(shù)，主瓣寬度與漢寧窗的相同。布萊克曼窗窗函數(shù)為 （3-15）其頻譜的幅度函數(shù)為 （3-16）該窗函數(shù)位移不同，幅度函數(shù)也不同，會使旁瓣進一步抵消，主瓣寬度為。凱塞窗是一種最優(yōu)窗函數(shù)，不同于前面五種窗函數(shù)，凱塞窗是一種參數(shù)可調(diào)的窗函數(shù)，其函數(shù)形式如下： （3-17）其中 （3-18）一般取15-25項可以滿足精度要求。參數(shù)可以控制窗的形狀。一般越大，主瓣越寬，而旁瓣幅度會隨之減

25、小，典型的數(shù)據(jù)在4到9之間。各種窗函數(shù)的性能比較如表3-1所示：表3-1 不同窗函數(shù)性能比較窗函數(shù)旁瓣峰值/dB近似過渡帶寬精確過渡帶寬阻帶最小衰減/dB矩形窗-134/N1.8/N21三角形窗-258/N6.1/N25漢寧窗-318/N6.2/N44哈明窗-418/N6.6/N53布萊克曼窗-5712/N11/N74凱澤窗（）-5710/N80雖然窗函數(shù)設(shè)計法設(shè)計思路簡單，但是它的邊界頻率不容易控制，而且窗函數(shù)還有吉布斯效應(yīng)，需要選擇不同的窗函數(shù)來減小吉布斯效應(yīng)對結(jié)果的影響，但無論哪種窗函數(shù)，都無法很好的解決這一問題，所以我們需要通過其他的設(shè)計方法來進行濾波，便于滿足實際工程中的不同要求。3

26、.2 FIR數(shù)字濾波器的頻率采樣設(shè)計法窗函數(shù)設(shè)計法是從時域出發(fā)來設(shè)計FIR數(shù)字濾波器的，而頻率采樣法是從頻域出發(fā)設(shè)計FIR數(shù)字濾波器的。和窗函數(shù)設(shè)計法相同，頻率采樣法也需要預(yù)先構(gòu)造一個希望逼近的濾波器頻率響應(yīng)函數(shù)，對其加以等間隔采樣后，作為FIR數(shù)字濾波器的頻率響應(yīng)。對在到之間等間隔采樣點，得到頻率采樣值： （3-19）再對進行點IDFT，得到： （3-20）將作為所涉及的FIR數(shù)字濾波器的單位脈沖響應(yīng)，其系統(tǒng)函數(shù)為為 （3-21）由于濾波器頻率響應(yīng)是理想的，即有間斷點，那么其單位沖激響應(yīng)是無限長的。這樣，由于時域混疊，引起所設(shè)計的h(n)和有偏差。因此，采樣點處與相等，逼近誤差為0，而在采樣

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

28、帶采樣點的個數(shù)m與濾波器阻帶最小衰減的經(jīng)驗數(shù)據(jù)列于表3-2中，我們可以根據(jù)給定的阻帶最小衰減，選擇過渡帶采樣點的個數(shù)m。表3-2 過渡帶采樣點的個數(shù)m與濾波器阻帶最小衰減的經(jīng)驗數(shù)據(jù)m1234454dB6575dB8595dB3.2.2 頻率采樣法的設(shè)計步驟首先根據(jù)阻帶最小衰減按照表3-2選擇過渡帶采樣點的個數(shù)，再確定過渡帶寬度，估算頻域采樣點數(shù)，如果增加個過渡帶采樣點，則過渡帶寬度近似變成。當(dāng)確定時，過渡帶會隨著的增大而變寬。如果給定的過渡帶寬度為，則要求，濾波器的長度必須滿足以下公式： （3-22）接著，構(gòu)造一個希望逼近的頻率響應(yīng)函數(shù)： （3-23）設(shè)計標(biāo)準(zhǔn)型片段常數(shù)特性的FIR數(shù)字濾波器時

29、，一般構(gòu)造幅度特性函數(shù)為相應(yīng)的理想頻響特性，且滿足的對稱情況。對（3-23）進行頻域采樣： （3-24） （3-25）并加入過渡帶采樣。過渡帶采樣值一般為經(jīng)驗值，或者用累試法確定，也可以采用優(yōu)化算法估算。對進行點IDFT，得到第一類線性相位FIR數(shù)字濾波器的單位脈沖響應(yīng)： （3-26） 圖3-2 頻率采樣法設(shè)計FIR濾波器流程最后檢驗設(shè)計結(jié)果，如果阻帶最小衰減未達到指標(biāo)要求，則要改變過渡帶的采樣值，直到滿足指標(biāo)要求為止。如果濾波器的邊界頻率未達到指標(biāo)要求，則需要微調(diào)的邊界頻率。頻率采樣法設(shè)計FIR濾波器流程如圖3-2所示。頻率采樣法最大的優(yōu)點就是直接從頻率進行設(shè)計，比較直觀，也適合于設(shè)計具有任

30、意幅度特性的濾波器。但是頻率采樣法在邊界頻率不容易控制，如果采樣點數(shù)N增加，對確定邊界頻率有好處，但同樣會增加了濾波器的成本，因此只適合窄帶濾波，且這種設(shè)計方法理解起來比較困難。3.3 FIR數(shù)字濾波器的等波紋逼近設(shè)計法窗函數(shù)設(shè)計法和頻率采樣設(shè)計法雖然設(shè)計方法簡單，但都存在濾波器邊緣頻率不易精確控制缺點，且這兩種設(shè)計方法設(shè)計出來的濾波器的通帶和阻帶的波動幅度都是相等的，兩種設(shè)計方法都不能分別控制通帶和阻帶的波動幅度，而現(xiàn)實工程中往往對二者都有不同的要求，需要分別進行控制。等波紋逼近法是一種優(yōu)化設(shè)計方法，它克服了窗函數(shù)設(shè)計法和頻率采樣法的缺陷，是最大誤差最小化設(shè)計方法，并在整個逼近頻段上均勻分布

31、。設(shè)為希望逼近的幅度特性函數(shù)，且要求設(shè)計線性相位的FIR數(shù)字濾波器時，必須滿足線性相位約束條件。用表示實際設(shè)計的幅度特性函數(shù)，定義加權(quán)誤差函數(shù)為 （3-27）式中，被稱為誤差加權(quán)函數(shù)，是由設(shè)計者定義的，用來控制不同頻段的逼近精度。經(jīng)過推導(dǎo)可把統(tǒng)一標(biāo)示為：式中，是系數(shù)不同的余弦組合式，記；是不同的常數(shù)，在設(shè)計FIR濾波器時存在四種線性相位，當(dāng)且奇對稱時，N為奇數(shù)，為1；N為偶數(shù)時，為；而當(dāng)偶對稱時，不管N為奇數(shù)還是偶數(shù)，都取。圖3-3 remez算法搜索最佳頻率點流程等波紋逼近的問題就是選擇N個系數(shù)a(n)，使加權(quán)誤差的最大值最小。設(shè)定誤差加權(quán)函數(shù)、通帶波紋峰值及N+1個初始頻率i,i=0,1,

32、2,N，計算各個頻率點上的加權(quán)誤差函數(shù)值，如果這些加權(quán)誤差函數(shù)值小于等于給定的通帶波紋峰值，則為最佳頻率點；如果某個頻率點上的加權(quán)誤差函數(shù)值大于給定的通帶波紋峰值，就用Remez算法更新迭代查找最佳頻率點。找到最佳頻率點后，計算系數(shù)及。Remez算法搜索最佳頻率點流程如圖3-3所示。借助Remez算法來設(shè)計等波紋的FIR數(shù)字濾波器的步驟如下：給出所需的頻率響應(yīng)、誤差加權(quán)函數(shù)和單位脈沖響應(yīng)的長度N，根據(jù)具體情況，進行統(tǒng)一表達式的轉(zhuǎn)換；設(shè)定初始頻率點和通帶波紋值，采用Remez算法計算最佳頻率點組，由最佳頻率點組計算及系數(shù)，根據(jù)求出的的表達式再利用傅里葉逆變換求得濾波器的單位脈沖響應(yīng)。用等波紋逼近

33、法設(shè)計線性相位FIR濾波器的流程如圖3-4所示：圖3-4 等波紋逼近法設(shè)計FIR濾波器流程利用等波紋逼近法設(shè)計FIR濾波器，其誤差均勻分布在頻帶中，可以得優(yōu)良的濾波特性，它在同樣過渡帶較窄的情況下，通帶最穩(wěn)定，阻帶有最大化的最小衰減。3.4 三種設(shè)計法之間的對比窗函數(shù)設(shè)計法是這三種設(shè)計法里操作最簡單的，但是窗函數(shù)的設(shè)計階數(shù)相對其他兩種更多，如果降低了設(shè)計的階數(shù)，則最后得到的設(shè)計結(jié)果往往不能達到實際設(shè)計的要求和指標(biāo)，除矩形窗之外，其他窗可以滿足一般的設(shè)計指標(biāo)，但在要求精度較高的時候，窗函數(shù)設(shè)計往往不能很好的達到要求。頻率采樣法設(shè)計濾波器相對于窗函數(shù)，其階數(shù)稍微比窗函數(shù)要少，但要通過優(yōu)化過渡帶的采

34、樣點取值來控制阻帶的波紋幅度，所以并不能很好的滿足工程中的要求。而等波紋逼近法是這三種方法中階數(shù)最少的，且通帶最大衰減是最小的，阻帶衰減最大。第4章 基于Matlab的FIR數(shù)字濾波器設(shè)計1.1 Matlab簡介1.1.1 MATLAB的發(fā)展MATLAB是英文MATrix LABoratory（矩陣實驗室）的縮寫。早期的MATLAB是用FORTRAN語言編寫的，盡管功能十分簡單，但作為免費軟件，還是吸引了大批使用者。經(jīng)過幾年的校際流傳，在John Little。Cleve Moler和Steve Banger合作，于1984年成立MathWorks公司，并正式推出MATLAB第一版版。從這時起

35、，MATLAB的核心采用C語言編寫，功能越來越強大，除原有的數(shù)值計算功能外，還新增了圖形處理功能。Matlab主要面對科學(xué)計算、可視化以及交互式程序設(shè)計的高科技計算環(huán)境。它將數(shù)值分析、矩陣計算、科學(xué)數(shù)據(jù)可視化以及非線性動態(tài)系統(tǒng)的建模和仿真等諸多強大功能集成在一個易于使用的視窗環(huán)境中，為科學(xué)研究、工程設(shè)計以及必須進行有效數(shù)值計算的眾多科學(xué)領(lǐng)域提供了一種全面的解決方案，代表了當(dāng)今國際科學(xué)計算軟件的先進水平。MATLAB語言被稱為第四代計算機語言，其利用豐富的函數(shù)資源，使程序員從繁瑣的程序代碼中解放出來，其最突出的特點就是簡潔。MATLAB用更直觀的、符合人們思維習(xí)慣的代碼，代替了C和FORTRAN

36、語言的冗長代碼，給用戶帶來最直觀、最簡潔的程序開發(fā)環(huán)境，下面簡單介紹一下MATLAB的主要特點。語言簡潔緊湊，使用方便，庫函數(shù)十分豐富。MATLAB程序書寫的形式自由，利用豐富的庫函數(shù)避開了繁瑣的子程序編程任務(wù)，由于庫函數(shù)都是由本領(lǐng)域的專家編寫，所以不必?fù)?dān)心函數(shù)的可靠性。高效方便的矩陣和數(shù)組運算，MATLAB語言不需要定義數(shù)組的維數(shù)，并給出了矩陣函數(shù)、特殊矩陣函數(shù)、特殊矩陣專門的庫函數(shù)，使得在求解信號處理、建模、系統(tǒng)識別、優(yōu)化和控制等領(lǐng)域的問題時，顯得大為簡潔、方便、高效，這是其他高級語言所不能的。MATLAB既具有結(jié)構(gòu)化的控制語句，又具有面向?qū)ο缶幊痰奶匦?。MATLAB語法限制不嚴(yán)格，程序設(shè)

37、計自由度大，通過建立M后綴名文件的形式，與用戶已經(jīng)編好的FORTRAN、C語言成語混合編程，方便地調(diào)用有關(guān)的FORTRAN、C語言的子程序?？梢浦残院芎?，基本上不做修改就可以在各種型號的計算機和操作系統(tǒng)上面運行。MATLAB的圖形功能強大。在C和FORTRAN語言里，繪圖都很不容易，但在MATLAB里，數(shù)據(jù)的可視化非常簡單。此外，MATLAB還具有較強的編輯圖形界面的能力。MATLAB擁有功能強大的工具箱，主要用來擴充其符號計算功能、圖示建模仿真功能、文字處理功能以及與硬件實施交互功能。源程序的開放性強。除內(nèi)部函數(shù)以外，所有MATLAB的核心文件和工具箱文件都是可讀可改變的源文件，用戶可通過對

38、源文件的修改以及加入自己的文件構(gòu)成新的工具箱。MATLAB軟件自1984年推向市場以來，歷經(jīng)十幾年的發(fā)展和競爭，現(xiàn)已成為國際公認(rèn)的最優(yōu)秀的科技應(yīng)用軟件。它功能強大、界面友好、語言自然、開放性強，很快成為應(yīng)用學(xué)科計算機輔助分析、設(shè)計、仿真、教學(xué)乃至科技文字吹不可缺少的基礎(chǔ)軟件。1.1.1 窗函數(shù)法的Matlab實現(xiàn) 實際設(shè)計時，一般采用MATLAB工具箱函數(shù)實現(xiàn)。步驟（1）由設(shè)計者根據(jù)設(shè)計指標(biāo)完成，步驟（2）-（4）的解題過程可調(diào)用設(shè)計函數(shù)fir1來實現(xiàn)。Fir1時用窗函數(shù)法設(shè)計線性相位FIRDF的工具箱函數(shù)，調(diào)用格式如下：hn=fir1(N,wc,ftype,window) fir實現(xiàn)線性相位

39、FIR濾波器的標(biāo)準(zhǔn)窗函數(shù)法設(shè)計?！皹?biāo)準(zhǔn)”時指再設(shè)計低通、高通、帶通、帶阻FIR濾波器時，分別取相應(yīng)的理想低通、高通、帶通、帶阻濾波器，故而設(shè)計的濾波器的頻率響應(yīng)稱為標(biāo)注頻率響應(yīng)。hn=fir1(N,wc)可得到6dB截至頻率為wc的N階（單位脈沖響應(yīng)h(n)長度為N+1）FIR低通濾波器，默認(rèn)（缺省參數(shù)windows）選用hammiing窗。其單位脈沖響應(yīng)h(n)為：h(n)=h(n+1),n=0,1,2,N而且滿足線性相位條件：h(n)=h(N-1-n)其中wc為對歸一化的數(shù)字頻率，0wc1。hn=fir1(N,wc,ftype)可設(shè)計高通和帶阻濾波器。·當(dāng)ftype=high時，

40、設(shè)計高通FIR濾波器；·當(dāng)ftype=stop時，設(shè)計帶阻FIR濾波器。hn=fir1(N,wc,ftype,window)通過選擇wc、ftype和window參數(shù)，可以設(shè)計各種加窗濾波器。Fir2可以指定任何形狀的Hd（），用help命令查閱其調(diào)用格式。下面給出利用fir1函數(shù)的窗函數(shù)法數(shù)字低通濾波器。利用fir1函數(shù)進行設(shè)計，這種設(shè)計方法只需要給出濾波器的階數(shù)、截止頻率、窗函數(shù)等參數(shù)，Matlab即可自行完成設(shè)計，并可通過freqz函數(shù)查看濾波器的幅頻響應(yīng)和相頻響應(yīng)，以驗證濾波器是否滿足設(shè)計要求。設(shè)計程序如下：n=27; %設(shè)濾器的階數(shù)為27Wn=0.2; %歸一化截止頻率為0

41、.2window=rectwin(28); %矩形窗b=fir1(n,Wn,window);freqz(b);運行結(jié)果如圖4-1所示：從幅頻響應(yīng)上看，通帶波動約為1dB且?guī)?nèi)有過沖現(xiàn)象，0.2截止頻率處的衰減為6.5dB，即截止頻率發(fā)生了漂移，3dB截止頻率為0.182；阻帶波紋較大，因而阻帶特性較不理想，阻帶最小衰減約為21dB；相頻響應(yīng)曲線在通帶內(nèi)為直線，效果較好，信號失真小。 圖4-1矩形窗FIR濾波器幅頻和相頻特性曲線利用三角窗進行設(shè)計時，原理與矩形窗基本相同，只不過生成窗函數(shù)時采用triang（n）函數(shù)生成三角窗，程序運行結(jié)果如圖4-2所示：圖4-2三角窗設(shè)計的FIR濾波器幅頻和相頻

42、特性曲線從幅頻響應(yīng)上看，通帶平坦沒有波動，阻帶最小衰減約為25dB；相頻響應(yīng)曲線在整個頻率范圍內(nèi)為直線，滿足線性相位要求。利用哈明窗進行設(shè)計時，原理與矩形窗相同，只不過生成窗函數(shù)時采用hamming（n）函數(shù)生成三角窗，程序運行結(jié)果如下：圖4-3哈明窗設(shè)計的FIR濾波器幅頻和相頻特性曲線利用布萊克曼窗進行設(shè)計時，原理與矩形窗基本相同，只不過生成窗函數(shù)時采用blackman（n）函數(shù)生成布萊克曼窗，程序運行結(jié)果如下：圖4-4布萊克曼窗設(shè)計的FIR濾波器幅頻和相頻特性曲線利用凱塞窗進行設(shè)計時，原理與矩形窗基本相同，只不過生成窗函數(shù)時采用kaiser（n,beta）函數(shù)生成凱塞窗，程序運行結(jié)果：圖4

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

44、特性曲線。下面用頻率采樣法設(shè)計一個通帶截止頻率為0.2的FIR數(shù)字低通濾波器，設(shè)濾波器階數(shù)n=27，編寫程序：n=27; %設(shè)濾器的階數(shù)為33F=0 0.2 0.3 1; %理想濾波器的頻率響應(yīng)A=1 0. 707 0 0; b=fir2(n,F,A);freqz(b);經(jīng)過運算得FIR濾波器的幅頻特性如圖4-6所示：圖4-6 頻率采樣法設(shè)計的FIR低通濾波器的幅頻和相頻特性曲線從4-6幅頻特性曲線可知，0.2截止頻率處的衰減為4.7dB，即截止頻率發(fā)生了漂移，3dB截止頻率為0.165；通帶內(nèi)的相頻特性為直線，滿足線性相位要求。4.4 等波紋逼近法的Matlab實現(xiàn)頻率采樣法設(shè)計FIR濾波器

45、流程如圖3-4所示，因此利用Matlab軟件設(shè)計FIR濾波器時，可以根據(jù)3-4的流程編寫程序?qū)崿F(xiàn)。但Matlab軟件工具箱提供了remez子函數(shù)，可以利用remez子函數(shù)快速方便地設(shè)計FIR濾波器，避免了繁雜的數(shù)學(xué)運算。remez函數(shù)有幾種調(diào)用格式，常見的一種調(diào)用格式為：b= remez (n, f, a, w); 函數(shù)參數(shù)說明如下：1.n表示濾波器的階數(shù)2. f為濾波器的歸一化頻率點，在0,1范圍內(nèi)取值，1表示實際頻率為fs/23.a為頻率響應(yīng)幅度參數(shù)f和a共同描述了預(yù)先設(shè)定的理想濾波器的頻率響應(yīng)。4.w表示通帶和阻帶波紋控制。下面用remez函數(shù)設(shè)計一個27階的FIR低通濾波器，其通帶截止

46、頻率為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); 程序運行結(jié)果如圖4-7所示：圖4-7 等波紋逼近法設(shè)計的FIR低通濾波器幅頻和相頻特性曲線同樣是設(shè)計一個27階的FIR低通數(shù)字濾波器，比較fir1窗函數(shù)法、fir2頻率采樣法和remez等波逼近法的設(shè)計結(jié)果，可以看出： （1）窗函數(shù)法和頻率采樣法設(shè)計濾波器時，濾波器邊緣頻率不易精確控制，通帶和阻帶的波動幅度 也不易控制。（2）利用等波紋逼近法可以較好的控制

47、濾波器的邊緣頻率及通帶和阻帶的波動幅度。 第5章濾波器使用的實例說明4.1語音信號的采集基于聲卡進行數(shù)字信號的采集。將話筒插入計算機的語音輸入插口上,啟動錄音機。按下錄音按鈕，對話筒說話,說完后停止錄音。要保存文件時,利用了計算機上的A/D轉(zhuǎn)換器,把模擬的聲音信號變成了離散的量化了的數(shù)字信號,放音時,它又通過D/A轉(zhuǎn)換器,把保存的數(shù)字?jǐn)?shù)據(jù)恢復(fù)為原來的模擬的聲音信號。在 Matlab軟件平臺下可以利用函數(shù)wavread對語音信號進行采樣,得到了聲音數(shù)據(jù)變量x1,同時把x1的采樣頻率fs=22050Hz和數(shù)據(jù)位Nbits=16Bit放進了MATALB的工作空間。圖4.1為原始語音信號的時域圖形。

48、圖4.1 原始語音信號從圖中可以看出在時域環(huán)境下,信號呈現(xiàn)出不規(guī)則的信號峰值。通過freqz函數(shù)繪制原始語音信號的頻率響應(yīng)圖4.2。圖4.2 語音信號頻率響應(yīng)圖然后對語音信號進行頻譜分析,在Matlab中可以利用函數(shù)fft對信號行快速傅里葉變換,得到信號的頻譜圖4.3， 圖4.3 原始語音信號FFT與信號頻譜從圖中可以看出對各個頻點上的隨機信號在頻域進行抽樣 ,抽樣頻率為 22050Hz。4.2語音信號加噪與頻譜分析 圖4.4 加噪后的信號時域圖在Matlab中人為設(shè)計一個固定頻率5500Hz的噪聲干擾信號。噪聲信號通常為隨機序列,在本設(shè)計中用正弦序列代替,干擾信號構(gòu)建命令函數(shù)為d=Au*si

49、n(2*pi*5500*t),給出的干擾信號為一個正弦信號，針對上面的語音信號 ,采集了其中一段。再對噪音信號進行頻譜變換得到其頻譜圖,從圖中可以看出干擾信號，在4000Hz和 6000Hz頻點處有一高峰 ,其中 5500Hz 正是本設(shè)計所要利用的。4.3設(shè)計巴特沃斯低通濾波器 圖4.5 原始信號與加噪后信號頻譜對比由模擬濾波器變換為數(shù)字濾波器時,采用的是雙線性變換法,它保留的是從模擬到數(shù)字域的系統(tǒng)函數(shù)表示。在濾波器的選取時, 由于設(shè)計方法的側(cè)重點不同, 作出比較是困難的。如果FIR濾波器情況下,最優(yōu)的設(shè)計將是橢圓濾波器。用雙線性變換法設(shè)計低通濾波器。的本設(shè)計是用雙線性變換法設(shè)計BW帶阻型濾波

50、器。在 MATLAB中 ,可以利用函數(shù) butterworth,設(shè)計FIR濾波器,利用MATLAB 中的函數(shù)freqz 畫出各濾波器的頻率響應(yīng)。用設(shè)計好的帶阻濾波器對含噪語音信號進行濾波,在Matlab用FIR濾波器利用函數(shù)fftfilter對信號進行濾波,FIR濾波器利用函數(shù)filter對信號進行濾波。在一個窗口同時畫出濾波前后的波形及頻譜。從圖中可以看出 ,5500Hz看到的高峰消失了 ,語音信號與開始的一樣 ,濾波器成功的濾除了干擾信號。利用MATLAB對語音信號進行分析和處理，采集語音信號后，利用MATLAB軟件平臺進行頻譜分析；并對所采集的語音信號加入干擾噪聲，對加入噪聲的信號進行頻

51、譜分析，設(shè)計合適的濾波器濾除噪聲，恢復(fù)原信號4。4.4用濾波器對加噪語音濾波 圖4.6 低通濾波器頻率響應(yīng)曲線 給原始的語音信號加上一個高頻余弦噪聲，頻率為5500Hz。對加噪后的語音進行分析，并畫出其信號時域和頻譜圖。與原始信號對比，區(qū)別：先原始信號沒加噪音之前0到2000有幅值，在4000到6000之間沒有幅值，但是在加了噪音之后4000到6000之間出現(xiàn)最大幅值12，超出正常值。4.5比較濾波前后語音信號波形及頻譜 圖4.7 濾波前后時域波形圖 通過對比分析可知，濾波后的輸出波形和原始語音加噪聲信號的圖形發(fā)生了一些變化，在3000-6000Hz之間的信號消失，出現(xiàn)0-1000和7000-

52、8000之間的信號。濾波后的輸出波形明顯在對應(yīng)時間幅度比原語音加噪聲信號的要小，而且濾波的效果也與濾波器的選擇有關(guān)，可以看出濾波器的性能差異以及參數(shù)的調(diào)節(jié)會對濾波器產(chǎn)生一定的影響，并且通過回放可以發(fā)現(xiàn)濾波前后的聲音有變化.低通濾波后,已很接近原來的聲音，人耳幾乎辨別不出 。從頻譜圖中我們還可以看出聲音的能量信號主要集中在低頻部分，說明高頻語音信號被濾出，濾波器達到既定要求。信噪比為0.62。第6章結(jié)論論文正文主要簡單介紹了Matlab、數(shù)字濾波器及利用matlab實現(xiàn)FIR濾波器的多種技術(shù)設(shè)計。Matlab語言簡潔緊湊，使用方便，庫函數(shù)十分豐富。MATLAB程序書寫的形式自由，利用豐富的庫函數(shù)

53、避開了繁瑣的子程序編程任務(wù)。在數(shù)字信號處理中 ,由于信號中經(jīng)?；煊懈鞣N復(fù)雜成分，所以很多信號分析都是基于濾波器而進行的， 因此數(shù)字濾波器占有極其重要的地位 。數(shù)字濾波器分為有限沖激響應(yīng)數(shù)字濾波器，即FIR數(shù)字濾波器和無限沖激響應(yīng)，即IIR數(shù)字濾波器。我們主要介紹了FIR數(shù)字濾波器。目前FIR濾波器的設(shè)計方法主要有三種：窗函數(shù)法、頻率抽樣法和優(yōu)化設(shè)計方法。我們主要介紹前兩種方法。涉及FIR濾波器的多種技術(shù)設(shè)計。各種方法都有其優(yōu)點和缺點,需根據(jù)不同的濾波器類型選擇不同的方法。窗函數(shù)法在設(shè)計標(biāo)準(zhǔn)濾波器，例如低通、高通、帶通，是很有用的。另一方面, 頻率抽樣法的優(yōu)點是可以在頻域直接設(shè)計，并且適合于最優(yōu)

54、化設(shè)計；缺點是抽樣頻率只能等于2pi/M的整數(shù)倍或等于2pi/M的整數(shù)倍上加上pi/M，因而不能確保截止頻率Wc的自由取值。要想實現(xiàn)自由選擇頻率，則必須增加抽樣點數(shù)M，但這種計算量加大。本文實現(xiàn)了基于Matlab的數(shù)字低通和數(shù)字高通濾波器。高通濾波器是容許高頻信號通過、但減弱（或減少）頻率低于截止頻率信號通過的濾波器。對于不同濾波器而言,每個頻率的信號的減弱程度不同。它有時被稱為低頻剪切濾波器；在音頻應(yīng)用中也使用低音消除濾波器或者噪聲濾波器。低通濾波器與高通濾波器特性恰恰相反。此次設(shè)計當(dāng)中有很多問題困擾我，通過查閱資料、同學(xué)和邵霞老師的幫助逐步解決了問題，在此艱難的過程中讓我懂得了很多。萬事開

55、頭難，不要畏懼，做好了開頭也就成功了一半；其發(fā)現(xiàn)問題要立即解決問題；自己鉆研不出來的，要敢于問問題；做事認(rèn)真仔細(xì)，不要怕麻煩，否則只會更麻煩。致 謝 首先，在邵霞老師開的第一次會議上，邵老師把她自己的一些構(gòu)想給我說了一下，這些構(gòu)想使我對完成FIR濾波器的設(shè)計有了很大的信心，在此特別感謝邵老師對我的悉心指導(dǎo)。在與邵老師接觸的過程中，她給我們的感覺就是我們就像朋友間的交流，毫無距離感。另外，邵老師思維縝密、知識淵博、生活上平易近人，使我從邵老師身上學(xué)到了很多專業(yè)的知識和做人的道理，這都是我以后做人做事的榜樣。其次，感謝通信教研室老師和領(lǐng)導(dǎo)，他們對我們在設(shè)計的過程中出現(xiàn)的問題同樣給與了一些的建議，我

56、也從他們身上學(xué)到了許多東西，這對我以后的生活與工作都有巨大的作用。再則，我要感謝和我一組的同學(xué)，由于我們一組的幾個同學(xué)用到的知識都有重疊，所以我們也在互幫互助中共同成長，在此過程中我深深體會到了團隊的重要性，雖然我們不是一個團隊，一個和諧團結(jié)的團隊是一只非?？膳碌年犖?，眾志成城，其利斷金。感謝和班里和系里的一些同學(xué)，他們也為我提供了無私的幫助，而且很重要，我們相互交流，相互學(xué)習(xí)，這個過程中我們一起解決了很多問題，同時也學(xué)到了很多新的知識。 最后我要感謝我的父母家人，他們對我的影響是最大的。謝謝，謝謝你們！參考文獻1 劉令普.數(shù)字信號處理M.哈爾濱：哈爾濱工業(yè)大學(xué)出版社,2002。2 羅軍輝、羅

57、勇江等. MATLAB7.0在數(shù)字信號處理中的應(yīng)用M.北京:機械工業(yè)出版社，20053 楊永昌,李晨輝,王凱. FIR數(shù)字濾波器的設(shè)計方法J. 桂林航天工業(yè)高等專科院校學(xué)報, 2006, 11, 43(3): 22-31. 附 錄clc;clear;close all;fs=22050; %ÓïÒôÐźŲÉÑùƵÂÊΪ22050x1=wavread('D:start_e.wav'); %

58、82;ÁÈ¡ÓïÒôÐźŵÄÊý¾Ý£¬¸³¸ø±äÁ¿x1£¬ sound(x1,22050); %²¥·ÅÓïÒôÐźÅy1=fft(x1,1024); %¶Ô

59、ÐźÅ×ö1024µãFFT±ä»»f=fs*(0:511)/1024; %½«0µ½511£¬²½³¤Îª1µÄÐòÁеÄÖµÓëfsÏà³Ë²¢³ýÒÔ1024µÄÖµ£¬¸³Öµ¸