精品【畢業(yè)設(shè)計(jì)】基于MATLAB環(huán)境的數(shù)字濾波器設(shè)計(jì)-定

1、河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）I基于基于 MATLABMATLAB 環(huán)境的數(shù)字濾波器設(shè)計(jì)環(huán)境的數(shù)字濾波器設(shè)計(jì)摘 要傳統(tǒng)的數(shù)字濾波器的設(shè)計(jì)過(guò)程復(fù)雜，計(jì)算工作量大，濾波特性調(diào)整困難，影響了它的應(yīng)用。本文介紹了一種利用 MATLAB 信號(hào)處理工具箱（Signal Processing Toolbox）快速有效的設(shè)計(jì)由軟件組成的常規(guī)數(shù)字濾波器的設(shè)計(jì)方法。利用 MATLAB 設(shè)計(jì)濾波器，可以隨時(shí)對(duì)比設(shè)計(jì)要求和濾波器特性調(diào)整參數(shù)，直觀簡(jiǎn)便，極大的減輕了工作量，有利于濾波器設(shè)計(jì)的最優(yōu)化。在本文中，我們分別研究了在 MATLAB 環(huán)境下 IIR 數(shù)字濾波器的典型設(shè)計(jì)和完全設(shè)計(jì)等方法。典型設(shè)計(jì)是先按一定規(guī)則將給出

2、的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器的技術(shù)指標(biāo)，據(jù)此產(chǎn)生模擬濾波器原型，然后把模擬低通濾波器原型轉(zhuǎn)換成模擬低通、高通、帶通、帶阻濾波器，最后再把模擬濾波器轉(zhuǎn)換成數(shù)字濾波器。完全設(shè)計(jì)方法中我們利用函數(shù)直接設(shè)計(jì)出低通、高通、帶通和帶阻濾波器。對(duì)于 FIR 數(shù)字濾波器的設(shè)計(jì)，可以根據(jù)所給定的頻率特性直接設(shè)計(jì)，文中采用的設(shè)計(jì)方法是窗函數(shù)法。我們給出了一個(gè)用窗函數(shù)法設(shè)計(jì)的帶阻濾波器實(shí)例與MATLAB 程序仿真。本文根據(jù) IIR 濾波器和 FIR 濾波器的特點(diǎn)，在 MATLAB 壞境下分別用雙線性變換法設(shè)計(jì) IIR 和用窗函數(shù)設(shè)計(jì) FIR 數(shù)字濾波器，并對(duì)采集的語(yǔ)音信號(hào)進(jìn)行分析，最后給出了 IIR

3、和 FIR 對(duì)語(yǔ)音濾波的效果。 關(guān)鍵詞：數(shù)字濾波器 ，MATLAB ，IIR， FIR 河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）IIBased on MATLAB environment digital filter designABSTRACTTraditional digital filter design process complex, the computation work load is big, the filter characteristic adjustment difficulty, has affected its application. This article introd

4、uced one kind using the MATLAB signal processing toolbox (Signal Processing Toolbox) the fast effective design the conventional numeral filter design method which is composed by software.Using the MATLAB design filter ,may contrast the design request and the filter characteristic adjustment paramete

5、r as necessary,direct-viewing is simple, enormous reduced the work load, is advantageous optimization which designs to the filter.In this thesis the typical and complete disigns under MATLAB are studied.The technical parameters from digital filters that should be disigned,and then transformed into t

6、he analog parameters of a low-pass analog filter prototype.The prototype is converted into the analog low-pass，high-pass，band-pass and the band-stop filters respeetively，Which are transformed into the digital ones.The complete design uses the given functions and releases the low-Pass，high-pass，band-

7、pass and the band-stop filters direetly. As for the FIR filter, we can design it directly based on the giving frequency, in the paper it uses the design of the window function. we are given a window function method with band-stop filter design examples with MATLAB simulation program. This article ac

8、cording to the IIR filter and the FIR filter characteristic, uses the bilinearity method of transformation under the MATLAB bad boundary to design IIR and to design the FIR numeral filter separately with the window box number, and carries on the analysis to the gathering pronunciation signal, and fi

9、nally gives IIR and FIR to the pronunciation filter effect.Key words: Digital filter，MATLAB ，IIR ，F(xiàn)IR 河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）III目目 錄錄前前言言.1第 1 章 數(shù)字濾波器及 MATLAB 語(yǔ)言概述.21.1 數(shù)字濾波器概述 .21.1.1 數(shù)字濾波器的定義.21.1.2 數(shù)字濾波器的分類.21.1.3 常用濾波器的性能指標(biāo).31.2 MATLAB 概述.6第 2 章 IIR 濾波器的設(shè)計(jì).8 82.1 IIR 濾波器設(shè)計(jì)的常用方法.82.2 雙線性變換法設(shè)計(jì) IIR 數(shù)字濾波器 .9

10、2.3 脈沖響應(yīng)不變法.102.4 完全濾波器的設(shè)計(jì).132.5 語(yǔ)音濾波應(yīng)用.18第 3 章 FIR 數(shù)字濾波器的設(shè)計(jì).243.1 典型的窗函數(shù).243.2 FIR 濾波器對(duì)語(yǔ)音信號(hào)的應(yīng)用.27結(jié)論.32參考文獻(xiàn).33謝辭.35外文資料翻譯.36河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）1前言隨著信息時(shí)代和數(shù)字世界的到來(lái)，數(shù)字信號(hào)處理已成為今一門極其重要的學(xué)科和技術(shù)領(lǐng)域。數(shù)字濾波器是數(shù)字信號(hào)處理的基礎(chǔ)，用來(lái)對(duì)信號(hào)進(jìn)行過(guò)濾、檢測(cè)與參數(shù)估計(jì)等處理，在通信、圖像、語(yǔ)音、雷達(dá)等許多領(lǐng)域都有著十分廣泛的應(yīng)用。尤其在圖像處理、數(shù)據(jù)壓縮等方面取得了令人矚目的進(jìn)展和成就1。鑒于此，數(shù)字濾波器的設(shè)計(jì)就顯得尤為重要。例如，在

11、電力系統(tǒng)微機(jī)保護(hù)和二次控制中，很多信號(hào)的處理與分析都是基于對(duì)正弦基波和某些整次諧波的分析，而在系統(tǒng)電壓電流信號(hào)（尤其是故障舜變過(guò)程）中各種復(fù)雜成分，所以濾波器一直是電力系統(tǒng)二次裝置的關(guān)鍵部件。目前微機(jī)保護(hù)和二次信號(hào)處理軟件主要是采用的數(shù)字濾波器。目前數(shù)字濾波器的設(shè)計(jì)有許多現(xiàn)成的高級(jí)語(yǔ)言設(shè)計(jì)程序，但他們都存在設(shè)計(jì)效率低，不具有可視圖形，不便于參數(shù)修改等缺點(diǎn)，而 MATLAB 的 FDATool 工具為數(shù)字濾波器的研究和應(yīng)用提供了一個(gè)直觀、高效、便捷的利器。它以矩陣運(yùn)算為基礎(chǔ)，把計(jì)算、可視化、仿真設(shè)計(jì)融合到一個(gè)交互式的工作環(huán)境中。尤其是MATLAB 工具箱使各個(gè)領(lǐng)域的研究人員可以直觀方便地進(jìn)行科學(xué)

12、研究和工程應(yīng)用。其中的信號(hào)處理、圖像處理等工具箱為數(shù)字濾波研究的蓬勃發(fā)展提供了有力的工具。本文介紹了一種基于 MATLAB 環(huán)境下對(duì)數(shù)字理濾波器進(jìn)行設(shè)計(jì)與仿真的方法。在基于 MATLAB 環(huán)境的數(shù)字濾波器的設(shè)計(jì)過(guò)程中，MATLAB 工具箱提供了幾種模擬濾波的原型產(chǎn)生函數(shù)貝塞爾低通模擬濾波器原型，巴特沃斯濾波器原型，切比雪夫（型、型）濾波器原型，橢圓濾波器原型等不同的模擬濾波器原型；模擬低通濾波器原型向低通、高通、帶通、帶阻的轉(zhuǎn)變函數(shù)；從模擬濾波器向數(shù)字濾波器轉(zhuǎn)化的雙線性變換和沖激響應(yīng)不變法；模擬數(shù)字 IIR 濾波器階數(shù)選擇函數(shù)以及數(shù)字濾波器直接設(shè)計(jì)函數(shù)等。這一整套設(shè)計(jì)函數(shù)給在 MATLAB 中

13、設(shè)計(jì) IIR 數(shù)字帶來(lái)了極大的方便。利用 MATLAB 設(shè)計(jì)濾波器，可以隨時(shí)對(duì)比設(shè)計(jì)要求和濾波器特性調(diào)整參數(shù)，直觀簡(jiǎn)便，極大地減輕了工作量，有利于濾波器設(shè)計(jì)的最優(yōu)化。河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）2 第 1 章 數(shù)字濾波器及 MATLAB 語(yǔ)言概述1.1 數(shù)字濾波器概述1.1.1 數(shù)字濾波器的定義數(shù)字濾波器是指完成信號(hào)濾波處理功能的，用有限精度算法實(shí)現(xiàn)的離散時(shí)間線性非時(shí)變系統(tǒng)，其輸入是一組數(shù)字量，其輸出是經(jīng)過(guò)變換的另一組數(shù)字量。因此，數(shù)字濾波器本身既可以是用數(shù)字硬件裝配成的一臺(tái)完成給定運(yùn)算的專用的數(shù)字計(jì)算機(jī)，也可以將所需要的運(yùn)算編成程序，讓通用計(jì)算機(jī)來(lái)執(zhí)行。數(shù)字濾波器，輸入輸出均為數(shù)字信號(hào)，通過(guò)

14、一定的運(yùn)算關(guān)系，改變輸入信號(hào)中所含頻率成分的相對(duì)比例，或則濾除某些頻率成分的器件2。對(duì)于數(shù)字濾波器而言，若系統(tǒng)函數(shù)為 H(z),其脈沖響應(yīng)為 h(n),輸入時(shí)間序列為 x(n),則它們?cè)跁r(shí)域內(nèi)的關(guān)系式如下：y(n)=h(n)x(n) （1-1）在 Z 域內(nèi)，輸入和輸出存在如下關(guān)系： Y(z)= H(z)X(z) （1-2）式中，X(z)、Y(z)分別為 x(n)和 y(n)的 Z 變換。在頻域內(nèi)，輸入和輸出則存在如下關(guān)系：Y(j)=H(j)X(j) （1-3）式中，H(j)是數(shù)字濾波器的頻率特性；X(j)Y(j)分別為 x(n)和 y(n)的頻譜，而為數(shù)字角頻率。 1.1.2 數(shù)字濾波器的分類

15、數(shù)字濾波器可以有很多種分類方法，但總體上可分為兩大類。一類稱為經(jīng)典濾波器，即一般的濾波器，其特點(diǎn)是輸入信號(hào)中的有用成分和希望濾除的成分占用不同的頻帶，通過(guò)合適的選頻濾波器可以實(shí)現(xiàn)濾波3。例如，若輸入信號(hào)中有干擾，信號(hào)和干擾的頻帶互不重疊，則可濾出信號(hào)中的干擾得到純信號(hào)。但是，如果輸入信號(hào)中信號(hào)和干擾的頻帶相重疊，則干擾就不能被有效的濾出。另一類稱為現(xiàn)代濾波器，如維納濾波器、卡爾曼濾波器等，其輸入信號(hào)中有用信號(hào)和希望濾除的成分頻帶重疊。對(duì)于經(jīng)典濾波器，從頻域上也可以分為低通、高通、帶河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）3通和帶阻濾波器。從時(shí)域特性上看，數(shù)字濾波器還可以分為有限沖激響應(yīng)數(shù)字濾波器（FIR）

16、和無(wú)限沖激響應(yīng)數(shù)字濾波器（IIR） 。對(duì)于有限沖激響應(yīng)數(shù)字濾波器（FIR） ，其輸出 y(n)只取決于有限個(gè)過(guò)去和現(xiàn)在的輸入，x(n),x(n-1),x(n-m),濾波器的輸入輸出關(guān)系可表示為y(n)=Mrrrnxb0)( （1-4）對(duì)于無(wú)限沖激響應(yīng)數(shù)字濾波器（IIR） ，它的輸出不僅取決于過(guò)去和現(xiàn)在的輸入，而且還取決于過(guò)去的輸出，其差分方程為y(n)+Nkkknya1)(=Mrrrnxb0)( (1-5)該差分方程的單位沖激響應(yīng)是無(wú)限延續(xù)的。 1.1.3 常用濾波器的性能指標(biāo)濾波器性能一般用系統(tǒng)頻率特性)(jeH來(lái)說(shuō)明，常用的性能指標(biāo)主有以下三個(gè)參數(shù)：1. 幅度平方函數(shù)2*()()*()()

17、() ( )( )jjjjjjz eH eH eHeH eH eH z Hz (1-6)該性能指標(biāo)主要用來(lái)說(shuō)明系統(tǒng)的幅頻特性。2. 相位函數(shù)()()Re()Im()()jjjjjjeH eH ejH eH ee (1-7)其中：)(Re)(Im)(jjjeHeHarctge (1-8)該指標(biāo)主要用來(lái)說(shuō)明系統(tǒng)的相位特性。3. 群延時(shí)河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）4dedj)()( （1-9）定義為相位對(duì)角頻率導(dǎo)數(shù)的負(fù)值，說(shuō)明了濾波器對(duì)不同的頻率成分的平均延時(shí)。當(dāng)要求在通帶內(nèi)的群延遲是常數(shù)時(shí)，濾波器相位響應(yīng)特性是線性的。實(shí)際設(shè)計(jì)中所能得到的濾波器的頻率特性與理想濾波器的頻率特性之間存在著一些顯著的差

18、別，現(xiàn)以低通濾波器的頻率特性為例進(jìn)行說(shuō)明。1.理想低通濾波器的特性：設(shè)濾波器輸入信號(hào)為)(tx，信號(hào)中混入噪音)(tu，它們有不同的頻率成分。濾波器的單位脈沖響應(yīng)為)(th。則理想濾波器輸出為：( ) ( )( )( )()y tx tu th tK x t （1-10）即噪音信號(hào)被濾除0)()(thtu，而信號(hào)無(wú)失真只有延時(shí)和線性放大。對(duì)上式作傅里葉變換得：()()()()()()jY jX jH jU jH jKeX j （1-11）假定噪音信號(hào)被濾除，即：()()0U jH j （1-12）整理得：()()()jY jH jKeX j （1-13）圖 1-1 理想低通濾波器頻率特性 圖

19、1-2 理想濾波器的單位脈沖響應(yīng)假定信號(hào)頻率成分為：c ，噪音頻率成分為c 。則完成濾波的理想低通濾波器特性是： |()()()0 |jccK eY jH jX j (1-14)即：()H jcK0h(t)河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）5 |()|0 |ccKH j (1-15)arg()H j (1-16)系統(tǒng)的單位脈沖響應(yīng)為：sin()1( )2()ccjj tcth tKeedKt (1-17)理想低通濾波器的頻率特性如圖 1-1 所示，單位脈沖響應(yīng)的波形如圖 1-2 所示。理想濾波器具有非因果、無(wú)限長(zhǎng)的單位脈沖響應(yīng)和不連續(xù)的頻率特性，要用穩(wěn)定的線性時(shí)不變（LTI）系統(tǒng)來(lái)實(shí)現(xiàn)這樣的特性是不

20、可能的4。工程上是用脈沖響應(yīng)為有限長(zhǎng)的、因果的、穩(wěn)定的線性時(shí)不變系統(tǒng)或具有連續(xù)頻率特性的線性時(shí)不變系統(tǒng)來(lái)逼近理想特性。在滿足一定的誤差要求的情況下來(lái)實(shí)現(xiàn)理想濾波特性。因此實(shí)際的濾波器的頻率特性如圖 1-3 所示。 其中：c截止頻率 s阻帶起始頻率cs過(guò)渡帶寬在通帶內(nèi)幅度響應(yīng)以1的誤差接近于 1，即：111()1jcH e (1-18)111102通 帶阻帶過(guò)度帶渡帶圖 1-3 實(shí)際濾波器頻率特性cs()jH e河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）6s為阻帶起始頻率，在阻帶內(nèi)幅度響應(yīng)以小于2的誤差接近于零，即：2()jsH e (1-19)為了使逼近理想低通濾波器的方法成為可能，還必須提供一帶寬為sc的

21、不為零的過(guò)渡帶。在這個(gè)頻帶內(nèi)，幅度響應(yīng)從通帶平滑的下落到阻帶。1.2 MATLAB 概述MATLAB是美國(guó)MathWorks公司開發(fā)的一種功能極其強(qiáng)大的高技術(shù)計(jì)算語(yǔ)言和內(nèi)容極其豐富的軟件庫(kù)，集數(shù)值計(jì)算、矩陣運(yùn)算和信號(hào)處理與顯示于一身。該軟件最初是由美國(guó)教授Cleve Moler創(chuàng)立的。1980年前后，他在教線性代數(shù)課程時(shí)，發(fā)現(xiàn)用其他高級(jí)語(yǔ)言編程時(shí)極不方便，便構(gòu)思開發(fā)了MATLAB，即矩陣實(shí)驗(yàn)室（Matrix Laboratory）。該軟件利用了當(dāng)時(shí)代表數(shù)值線性代數(shù)領(lǐng)域最高水平的EISPACK和LINPACK兩大軟件包，并且利用Fortran語(yǔ)言編寫了最初的一套交互式軟件系統(tǒng)，MATLAB的最初

22、版本便由此產(chǎn)生了5。 最初的MATLAB由于語(yǔ)言單一，只能進(jìn)行矩陣的運(yùn)算，繪圖也只能用原始的描點(diǎn)法，內(nèi)部函數(shù)只有幾十個(gè)，因此功能十分簡(jiǎn)單。1984年該公司推出了第一個(gè)MATLAB的商業(yè)版，并用C語(yǔ)言作出了全部改寫?，F(xiàn)在的MATLAB程序是MathWorks公司用C語(yǔ)言開發(fā)的，第一版由steve Bangert主持開發(fā)編譯解釋程序，Steve Kleiman完成圖形功能的設(shè)計(jì)，John Little和Cleve Moler主持開發(fā)了各類數(shù)學(xué)分分析的子模塊，撰寫用戶指南和大部分的M文件。接著又添加了豐富的圖形圖像處理、多媒體功能、符號(hào)運(yùn)算和與其它流行軟件的接口功能，使MATLAB的功能越來(lái)越強(qiáng)大。

23、MTALAB系統(tǒng)主要由以下五個(gè)部分組成6：(1)MATALB語(yǔ)言體系。 MATLAB是高層次的矩陣數(shù)組語(yǔ)言，具有條件控制、函數(shù)調(diào)用、數(shù)據(jù)結(jié)構(gòu)、輸入輸出、面向?qū)ο蟮瘸绦蛘Z(yǔ)言特性。利用它既可以進(jìn)行小規(guī)模端程，完成算法設(shè)計(jì)和算法實(shí)驗(yàn)的基本任務(wù)，也可以進(jìn)行大規(guī)模編程，開發(fā)復(fù)雜的應(yīng)用程序。(2)MATLAB工作環(huán)境 。這是對(duì)MATLAB提供給用戶使用的管理功能的總稱。包括管理工作空間中的變量據(jù)輸入輸出的方式和方法，以及開發(fā)、調(diào)試、管理M河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）7文件的各種工具。(3)圖形句相系統(tǒng) 。這是MATLAB圖形系統(tǒng)的基礎(chǔ)，包括完成2D和3D數(shù)據(jù)圖示、圖像處理、動(dòng)畫生成、圖形顯示等功能的高層M

24、ATLAB命令，也包括用戶對(duì)圖形圖像等對(duì)象進(jìn)行特性控制的低層MATLAB命令，以及開發(fā)GUI應(yīng)用程序的各種工具。(4)MATLAB數(shù)學(xué)函數(shù)庫(kù)。這是對(duì)MATLAB使用的各種數(shù)學(xué)算法的總稱。包括各種初等函數(shù)的算法，也包括矩陣運(yùn)算、矩陣分析等高層次數(shù)學(xué)算法。(5)MATLAB應(yīng)用程序接口(API)。這是MATLAB為用戶提供的一個(gè)函數(shù)庫(kù)，使得用戶能夠在MATLAB環(huán)境中使用C程序或FORTRAN程序，包括從MATLAB中調(diào)用于程序(動(dòng)態(tài)鏈接)，讀寫MAT文件的功能7。 除此之外，MATLAB系統(tǒng)還具有如下特點(diǎn)：（1）具有易學(xué)易用的語(yǔ)言體系；（2）具有交互式的工作環(huán)境；（3）具有多層面的圖像處理系統(tǒng)；

25、（4）具有豐富高效的MATLAB工具箱；（5）具有便利的程序接口（API）;（6）應(yīng)用領(lǐng)域廣泛；（7）嵌入了面向?qū)ο缶幊陶Z(yǔ)言。河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）8第 2 章 IIR 濾波器的設(shè)計(jì)2.1 IIR 濾波器設(shè)計(jì)的常用方法IIR 數(shù)字濾波器的系統(tǒng)函數(shù)為：NkkkMrrrnnzazbznhzXzYzH1001)()()()( （2-1）IIR 濾波器設(shè)計(jì)就是根據(jù)性能指標(biāo)要求，設(shè)計(jì)濾波器的分子和分母多項(xiàng)式。IIR 數(shù)字濾波器的設(shè)計(jì)通常借助于模擬濾波器原型，再將模擬濾波器轉(zhuǎn)換為數(shù)字濾波器。模擬濾波器的理論和設(shè)計(jì)方法已發(fā)展得相當(dāng)成熟，且有一些典型的模擬濾波器供我們選擇，如巴特沃斯(Butterwor

26、th)濾波器、切比雪夫(Chebyshev)濾波器、橢圓(Cauer)濾波器、貝塞爾(Bessel)濾波器等，這些典型的濾波器各有特點(diǎn)。本文采用經(jīng)典設(shè)計(jì)法中的模擬濾波器變換法和完全設(shè)計(jì)函數(shù)法分別來(lái)設(shè)計(jì)數(shù)字濾波器。用 MATLAB 進(jìn)行典型的數(shù)字濾波器的設(shè)計(jì)，一般步驟如下8:(1)按一定規(guī)則將給出的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成模擬低通濾波器的技術(shù)指標(biāo);(2)根據(jù)轉(zhuǎn)換后的技術(shù)指標(biāo)使用濾波器階數(shù)選擇函數(shù)，確定最小階數(shù) N 和固有頻率 Wn，根據(jù)選用的模擬低通濾波器的類型可分別用函數(shù):buttord、cheb1ord、cheb2ord、ellipord 等;(3)運(yùn)用最小階數(shù) N 產(chǎn)生模擬濾波器原型，模

27、擬低通濾波器的創(chuàng)建函數(shù)有:buttap、cheb1ap、cheb2ap、ellipap 等;(4)運(yùn)用固有頻率 Wn 把模擬低通濾波器原型轉(zhuǎn)換成模擬低通、高通、帶通、帶阻濾波器，可分別用函數(shù) lp2lp、lp2hp、lp2bp、lp2bs;(5)運(yùn)用沖激響應(yīng)不變法或雙線性變換法把模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，分別用函數(shù) impinvar 和 bilinear 來(lái)實(shí)現(xiàn)。河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）92.2 雙線性變換法設(shè)計(jì) IIR 數(shù)字濾波器在 MATLAB 中，雙線性變換法的調(diào)用函數(shù)是 bilinear。其調(diào)用格式為：a. zd,pd,kd= bilinear(z,p,k,fs)b. zd,pd

28、,kd= bilinear(z,p,k,fs,fp)c. numd,dend=bilinear(num,den,fs)d. numd,dend=bilinear(num,den,fs,fp)e. Aa,Bb,Cc,Dd=bilinear(A,B,C,D,fs)f. Aa,Bb,Cc,Dd=bilinear(A,B,C,D,fs,fp)zd,pd,kd=bilinear(z,p,k,fs)是把模擬濾波器的零極點(diǎn)模型轉(zhuǎn)換為數(shù)字濾波器的零極點(diǎn)模型，fs 為采樣頻率，z,p,k 分別為濾波器的零點(diǎn)、極點(diǎn)和增益；numd,dend=bilinear(num,den,fs)是把模擬濾波器的傳遞函數(shù)模型轉(zhuǎn)換

29、為數(shù)字濾波器的傳遞模型；Aa,Bb,Cc,Dd=bilinear(A,B,C,D,fs)是把模擬濾波器的狀態(tài)方程模型轉(zhuǎn)換為數(shù)字濾波器狀態(tài)方程模型。例如，用雙線性變換法設(shè)計(jì)一個(gè)巴特沃思數(shù)字低通濾波器，技術(shù)指標(biāo)如下：通帶截止頻率p =24k rad/s，阻帶截止頻率s =28k rad/s，通帶波紋系數(shù)pR =0.3dB, 阻帶波紋系數(shù)sR=50dB,采樣頻率fs=20000Hz。 程序如下：N,Wn=buttord(wp,ws,Rp,Rs,s) %估計(jì)濾波器最小階數(shù)z,p,k=buttap(N);Bap,Aap=zp2tf(z,p,k);b,a=lp2lp(Bap,Aap,Wn);bz,az=b

30、ilinear(b,a,fs)freqz(bz,az,Nn,fs) 程序在 MATLAB 環(huán)境下的運(yùn)行及結(jié)果如圖 2-1 所示：結(jié)果如下：N =11Wn =1.4892e+004bz =河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）10 Columns 1 through 6 0.0110 0.1211 0.6055 1.8166 3.6333 5.0866 Columns 7 through 12 5.0866 3.6333 1.8166 0.6055 0.1211 0.0110az = Columns 1 through 61.0000 2.7098 4.6379 5.2252 4.3685 2.7207

31、Columns 7 through 121.2885 0.4561 0.1181 0.0211 0.0023 0.000105001000150020002500-1000-800-600-400-2000Frequency (Hz)Phase (degrees)05001000150020002500-300-200-1000Frequency (Hz)Magnitude (dB)圖 2-1 巴特沃思數(shù)字低通濾波器幅頻-相頻特性2.3 脈沖響應(yīng)不變法 脈沖響應(yīng)不變法就是要求數(shù)字濾波器的脈沖響應(yīng)序列h(n)與模擬濾波器的脈沖響應(yīng))(tha的采樣值相等，即h(n)=)(thanTt=)(nTha

32、 （2-2） 式中，T為采樣周期。根據(jù)模擬信號(hào)的拉普拉斯變換與離散序列的Z變換之間的關(guān)系，我們知道)(1)(skaezjkSHTzHsT （2-3） 此式表明，)(tha的拉普拉斯變換在s平面上沿虛軸，按照周期s =2/T延拓河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）11后，按式z=STe，進(jìn)行Z變換，就可以將Ha(s)映射為H(z)。事實(shí)上，用脈沖響應(yīng)不變法設(shè)計(jì)IIR濾波器，只適合于Ha(s)有單階極點(diǎn),且分母多項(xiàng)式的階次高于分子多項(xiàng)式階次的情況5。將Ha(s)用部分分式表示： Ha(s)=LTh)(ta=NiiissA1 （2-4） 式中，LT代表拉普拉斯變換，is為的單階極點(diǎn)。將Ha(s)進(jìn)行拉普拉斯

33、反變換，即可得到)(tha=NitSitueAi1)( （2-5） 式中，u(t)是單位階躍函數(shù)。則)(tha的離散序列（對(duì)它進(jìn)行抽樣）：h(n)=)(nTha=NinTSinTueAi1)( （2-6）對(duì)h(n)進(jìn)行z變換之后，可以得到數(shù)字濾波器的系統(tǒng)函數(shù)H(z)H(z)=0)(nnznh= NiTSizeAi111 （2-7） 對(duì)比Ha(s)與H(z)，我們會(huì)發(fā)現(xiàn)：s域中Ha(s)的極點(diǎn)是is，映射到z平面之后，其極點(diǎn)變成了TSie，而系數(shù)沒(méi)有發(fā)生變化，仍為iA。因此，在設(shè)計(jì)IIR濾波器時(shí)，我們只要找出模擬濾波器系統(tǒng)函數(shù)Ha(s)的極點(diǎn)和系數(shù)iA，通過(guò)脈沖響應(yīng)不變法，代入H(z)的表達(dá)式中

34、，即可求出H(z)，實(shí)現(xiàn)連續(xù)系統(tǒng)的離散化8。但是脈沖響應(yīng)不變法只適合于設(shè)計(jì)低通和帶通濾波器，而不適合于設(shè)計(jì)高通和帶阻濾波器。因?yàn)?，如果模擬信號(hào))(tha的頻帶不是介于T/之間，則會(huì)在T/的奇數(shù)倍附近產(chǎn)生頻率混疊現(xiàn)象，映射到z平面后，則會(huì)在附近產(chǎn)生頻率混疊現(xiàn)象11。從而使所設(shè)計(jì)的數(shù)字濾波器不同程度的偏離模擬濾波器在附近的頻率特性，嚴(yán)重時(shí)使數(shù)字濾波器不滿足給定的技術(shù)指標(biāo)。為此，希望設(shè)計(jì)的濾波器是帶限濾波器，如果不是帶限的，例如，高通濾波器、帶阻濾波器，需要在高通濾波器、帶阻濾波器之前加保護(hù)濾波器，濾出高于折疊頻率/T以上的頻帶，以免產(chǎn)生頻率混疊現(xiàn)象9。但這樣會(huì)增加系統(tǒng)的成本和復(fù)雜性。因此，高通與帶

35、阻濾波器不適合用這種方法。在MATLAB中，脈沖響應(yīng)不變法的調(diào)用函數(shù)是impinvar，其調(diào)用格式為：河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）12a. bz,az=impinvar(b,a,fs)b. bz,az=impinvar(b,a)c. bz,az=impinvar(b,a,fs,tol)該函數(shù)的功能是將分子向量為 b、分母向量為 a 的模擬濾波器，轉(zhuǎn)換為分子向量為 bz、分母向量為 az 的數(shù)字濾波器。fs 為采樣頻率，單位為 Hz，默認(rèn)值為1Hz。tol 指誤差容限，表示轉(zhuǎn)換后的離散系統(tǒng)函數(shù)是否有重復(fù)的極點(diǎn)。例如，用脈沖響應(yīng)不變法設(shè)計(jì)一個(gè)契比雪夫 1 型數(shù)字低通濾波器，指標(biāo)要求：通帶截止頻率p

36、 =1000Hz，阻帶截止頻率s =1200Hz，采樣頻率 fs=5000Hz，通帶衰減系數(shù)pR =0.3dB，阻帶衰減系數(shù)sR=40dB。程序如下：wp=1000*2*pi;ws=1200*2*pi;fs=2500;Rp=0.3;Rs=40;N,Wn=cheb1ord(wp,ws,Rp,Rs,s); %估計(jì)濾波器最小階數(shù)z,p,k=cheb1ap(N,Rp); %模擬濾波器函數(shù)引用A,B,C,D=zp2ss(z,p,k); %返回狀態(tài)轉(zhuǎn)移矩陣形式AT,BT,CT,DT=lp2lp(A,B,C,D,Wn); %頻率轉(zhuǎn)換b,a=ss2tf(AT,BT,CT,DT); %返回傳遞函數(shù)形式 bz,a

37、z=impinvar(b,a,fs); %調(diào)用脈沖相應(yīng)不變法H,W=freqz(bz,az); %返回頻率響應(yīng)plot(W*fs/(2*pi),abs(H); %畫圖grid;xlabel(frequency/Hz); ylabel(magnitude);N,WnN =11 Wn =6.2832e+003程序在MATLAB環(huán)境下的運(yùn)行及結(jié)果如圖2-2所示：河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）13020040060080010001200140000.811.21.4frequency/Hzmagnitude圖 2-2 契比雪夫型數(shù)字低通濾波器幅頻特性曲線2.4 完全濾波器的設(shè)計(jì)在 M

38、ATLAB 信號(hào)處理工具箱提供了幾個(gè)直接設(shè)計(jì) IIR 數(shù)字濾波器的函數(shù)，所用到的函數(shù)有 butter、cheby1、cheb2ord、ellipd 以及 besself 等。Butter 用來(lái)直接設(shè)計(jì)巴特沃思數(shù)字濾波器，cheby1 用來(lái)直接設(shè)計(jì)切比雪夫型濾波器，cheb2ord 用來(lái)設(shè)計(jì)切比雪夫型濾波器，ellipd 用來(lái)設(shè)計(jì)橢圓濾波器，besself 用來(lái)設(shè)計(jì)貝塞爾濾波器。完全濾波器的設(shè)計(jì)就是根據(jù)設(shè)計(jì)要求，直接調(diào)用函數(shù)來(lái)設(shè)計(jì)所需要的數(shù)字濾波器。主要函數(shù)的特點(diǎn)如下：1、設(shè)計(jì) Butterworth 濾波器用函數(shù) butter()，可以設(shè)計(jì)低通、高通、帶通和帶阻的數(shù)字和模擬濾波器，其特性是通帶

39、內(nèi)的幅度響應(yīng)最大限度的平滑，但損失了截止頻率處的下降斜度。2、設(shè)計(jì) Chebyshevl 型濾波器用函數(shù) chebyl()?？梢栽O(shè)計(jì)低通、高通、帶通和帶阻的數(shù)字和模擬 Chebyshevl 型濾波器，其通帶內(nèi)為等波紋，阻帶內(nèi)為單調(diào)。Chebyshevl 型濾波器的下降斜度比 II 型大，但其代價(jià)是通帶內(nèi)波紋較大。3、設(shè)計(jì) Chebyshev11 型濾波器用函數(shù) cheby2()。可以設(shè)計(jì)低通、高通、帶通和帶阻的數(shù)字和模擬 ChebyshevII 型濾波器，其通帶內(nèi)為單調(diào)，阻帶內(nèi)等波紋。ChebyshevII 型濾波器的下降斜度比 I 型小，但其阻帶內(nèi)波紋較大。河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）144

40、、設(shè)計(jì)橢圓濾波器用函數(shù) ellip()，與 chebyl、cheby2 類似，可以設(shè)計(jì)低通、高通、帶通和帶阻的數(shù)字和模擬濾波器。與 Butterworth 和 Chebyshev 濾波器相比，ellip 函數(shù)可以得到下降斜度更大的濾波器，得通帶和阻帶均為等波紋。一般情況下，橢圓濾波器能以最低的階實(shí)現(xiàn)指定的性能指標(biāo)。（1）以上幾種函數(shù)在低通濾波器設(shè)計(jì)中的比較：設(shè):Wp=30Hz，Ws=35Hz，F(xiàn)s=100，Rp=0.sdB，Rs=4odB 分別用巴特沃斯(Butterworth)濾波器、切比雪夫(chebyshev)濾波器、橢圓(Cauer)濾波器，程序設(shè)計(jì)如下:巴特沃斯低通濾波器:nl，Wn

41、l=buttord(wp/(Fs/2)，ws/(Fs/2)，rp，rs， z);numl，denl=butter(nl，Wnl);切比雪夫 I 型低通濾波器:n2，Wn2=cheb1ord(wp/(Fs/2)，ws/(Fs/2)，rp，rs， z):num2，den2=ehebyl(n2，rp，Wn2);切比雪夫 II 型低通濾波器:n3，Wn3=ehebZord(wp/(Fs/2)，ws/(Fs/2)，rp，rs， z);num3，den3=cheby2(n3，rp，Wn3):橢圓低通濾波器:n4，Wn4=ellipord(wp/(Fs/2)，ws/(Fs/2)，rp，rs， z):num4

42、，den4=ellip(n3，rp，rs，Wn3);設(shè)計(jì)出的低通濾波器的頻率響應(yīng)如圖 2-3 所示。河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）15圖 2-3 低通數(shù)字濾波器頻率響應(yīng) (a)Butterworth 低通濾波器 (b)Chebyshevl 型低通濾波器 (c)Chebyshevll 型低通濾波器(d)橢圓低通濾波器河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）16（2）幾種類型在高通濾波器設(shè)計(jì)中的比較13:設(shè):Wp=35Hz，Ws=30Hz，F(xiàn)s=100，Rp=0.5dB，Rs=40dB 分別用巴特沃斯(Butterworth)濾波器、切比雪夫(Chebyshev)濾波器、橢圓(Cauer)濾波器。程序設(shè)計(jì)如下:

43、butterworth 高通濾波器:nl，Wnl=buttord(wp/(Fs/2)，ws/(Fs/2)，rp，rs， z);numl，denl=butter(nl，Wnl， high);chebyshevl 高通濾波器:n2，Wn2=eheblord(wp/(Fs/2)，ws/(Fs/2)，rp，rs， z):num2，den2=chebyl(n2，rp，Wn2， high):chebyshevII 高通濾波器:n3，Wn3=Cheb2ord(wp/(Fs/2)，ws/(Fs/2)，rp，rs， z):num3，den3=cheby2(n3，rp，Wn3， high):橢圓高通濾波器:n4，

44、Wn4=ellipord(wp/(Fs/2)，ws/(Fs/2)，rp，rs， z);num4，den4=ellip(n4，rp，rs，Wn4， high):設(shè)計(jì)出的高通濾波器的頻率響應(yīng)如圖 2-4 所示。河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）17圖 2-4 高通數(shù)字濾波器頻率響應(yīng)(a)Butterworth 高通濾波器(b)ChebyshevI 型高通濾波器(c)ChebyshevII 型高通濾波(d)橢圓高通濾波器河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）18從頻率響應(yīng)圖中可以看出：巴特沃斯濾波器具有單調(diào)下降的幅頻特性，通帶內(nèi)平滑;切比雪夫 I 型濾波器的幅頻特性在通帶內(nèi)有波動(dòng)，阻帶內(nèi)單調(diào);chebyshevII

45、 型濾波器的幅頻特性在阻帶內(nèi)有波動(dòng)，通帶內(nèi)單調(diào);橢圓濾波器的選擇性相對(duì)前三種是最好的，下降斜度比較大，通帶和阻帶內(nèi)均為等波紋，同樣的性能指標(biāo)，橢圓濾波器可以最低的階數(shù)來(lái)實(shí)現(xiàn)。這樣根據(jù)不同的要求可以選用不同類型的濾波器。2.5 語(yǔ)音濾波應(yīng)用1語(yǔ)音信號(hào)的采集：用一段語(yǔ)音，在 MATLAB 軟件平臺(tái)下，利用函數(shù) wavread 對(duì)語(yǔ)音信號(hào)進(jìn)行采樣，記住采樣頻率和采樣點(diǎn)數(shù)。z1,fs,bits=wavread(1.wav)運(yùn)行結(jié)果：fs =44100bits =16wavread 函數(shù)調(diào)用格式：y=wavread(file)，讀取 file 所規(guī)定的 wav 文件，返回采樣值放在向量 y 中。y,fs

46、,nbits=wavread(file)，采樣值放在向量 y 中，fs 表示采樣頻率（Hz） ，nbits 表示采樣位數(shù)。y=wavread(file,N)，讀取前 N 點(diǎn)的采樣值放在向量 y 中。y=wavread(file,N1,N2)，讀取從 N1 點(diǎn)到 N2 點(diǎn)的采樣值放在向量 y 中。2語(yǔ)音信號(hào)的頻譜分析首先畫出語(yǔ)音信號(hào)的時(shí)域波形z1=wavread(1.wav);plot(z1);圖像輸出如圖 2-5對(duì)語(yǔ)音信號(hào)進(jìn)行頻譜分析，在 MATLAB 中，可以利用函數(shù) fft 對(duì)信號(hào)進(jìn)行快速付立葉變換，得到信號(hào)的頻譜特性z1=wavread(1.wav);y1=z1(1:8192);Y1=f

47、ft(y1);河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）19n=0:8191;plot(n,Y1);圖像輸出如圖 2-600.511.522.533.544.5x 104-0.8-0.6-0.4-0.60.8 0100020003000400050006000700080009000-600-400-2000200400600圖 2-5 信號(hào)時(shí)域波形 圖 2-6 信號(hào)頻譜分析圖3. 根據(jù)語(yǔ)音信號(hào)的特點(diǎn)，設(shè)計(jì)出不同性能的數(shù)字濾波器，并用 MATLAB 進(jìn)行仿真。根據(jù)語(yǔ)音信號(hào)的特點(diǎn)給出有關(guān)濾波器的性能指標(biāo)：1）低通濾波器性能指標(biāo)，fp=1000Hz，fc=1200 Hz， As=100dB，

48、Ap=1dB；2）高通濾波器性能指標(biāo)，fc=2800 Hz，fp=3000 Hz As=100dB，Ap=1dB；3）帶通濾波器性能指標(biāo)，fp1=1200 Hz，fp2=3000 Hz，fc1=1000 Hz，fc2=3200 Hz，As=100dB，Ap=1dB。用 IIR 濾波器對(duì)其分析，最后，利用 MATLAB 中的函數(shù) freqz 畫出各濾波器的頻率響應(yīng)。設(shè)計(jì)程序及仿真結(jié)果如下：1設(shè)計(jì) IIR 低通濾波器濾波程序如下：clear;close allz1,fs,bits=wavread(1.wav)y1=z1(1:8192);Y1=fft(y1);fp=1000;fc=1200;As=1

49、00;Ap=1; ;Fs=8000;河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）20wc=2*fc/Fs;wb=2*fp/Fs;n,wp=cheb1ord(wc,wb,Ap,As);b,a=cheby1(n,Ap,wp);figure(1);freqz(b,a);x=filter(b,a,z1);X=fft(x,8192);figure(2);subplot(2,2,1);plot(abs(Y1);axis(0,1000,0,1.0);title(濾波前信號(hào)頻譜);subplot(2,2,2);plot(abs(X);axis(0,4000,0,0.03);title(濾波后信號(hào)頻譜);subplot(2,2

50、,3);plot(z1);title(濾波前信號(hào)波形);subplot(2,2,4);plot(x);title(濾波后信號(hào)波形);sound(x,fs,bits);仿真結(jié)果如圖下：01-2000-1500-1000-5000Normalized Frequency ( rad/sample)Phase (degrees)01-600-400-2000Normalized Frequency ( rad/sample)Magnitude (dB)圖 2-7 IIR 數(shù)字低通濾波器幅頻-相

51、頻特性河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）210500100000.81信 信 信 信 信 信 信0100020003000400000.010.020.03信 信 信 信 信 信 信0246x 104-1-0.500.51信 信 信 信 信 信 信0246x 104-1-0.500.51信 信 信 信 信 信 信圖 2-8 濾波前后信號(hào)頻譜和波形對(duì)比2設(shè)計(jì) IIR 高通濾波器濾波程序設(shè)計(jì)如下：clear;close allz1,fs,bits=wavread(1.wav)y1=z1(1:8192);Y1=fft(y1);fc=2800 ;fp=3000 ;As=100;Ap=1;

52、Fs=8000;wc=2*fc/Fs;wb=2*fp/Fs;n,wp=cheb1ord(wc,wb,Ap,As);b,a=cheby1(n,Ap,wp,high);figure(1);freqz(b,a);x=filter(b,a,z1);X=fft(x,8192);figure(2);subplot(2,2,1);plot(abs(Y1);axis(0,1000,0,1.0);title(濾波前信號(hào)頻譜);subplot(2,2,2);plot(abs(X);axis(0,4000,0,0.03);title(濾波后信號(hào)頻譜);subplot(2,2,3);plot(z1);title(濾波

53、前信號(hào)波形);subplot(2,2,4);plot(x);河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）22title(濾波后信號(hào)波形);sound(x,fs,bits);仿真結(jié)果如下：01-2000-1500-1000-5000Normalized Frequency ( rad/sample)Phase (degrees)01-600-400-2000Normalized Frequency ( rad/sample)Magnitude (dB)圖 2-9 IIR 數(shù)字高通濾波器幅頻-相頻特性050

54、0100000.81信 信 信 信 信 信 信0100020003000400000.010.020.03信 信 信 信 信 信 信0246x 104-1-0.500.51信 信 信 信 信 信 信0246x 104-1-0.500.51信 信 信 信 信 信 信圖 2-10 濾波前后信號(hào)頻譜和波形對(duì)比3設(shè)計(jì) IIR 帶通濾波器濾波程序設(shè)計(jì)如下：clear;close allz1,fs,bits=wavread(1.wav)y1=z1(1:8192);Y1=fft(y1);fp1=1200 ;fp2=3000; fc1=1000 ;fc2=3200 ;As=100;Ap=1;

55、 Fs=8000;wc=2*fc1/Fs,2* fc2/Fs;wb=2*fp1/Fs,2*fp2/Fs; n,wp=cheb1ord(wc,wb,Ap,As);河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）23b,a=cheby1(n,Ap,wp,stop);figure(1);freqz(b,a);x=filter(b,a,z1);X=fft(x,8192);figure(2);subplot(2,2,1);plot(abs(Y1);axis(0,1000,0,1.0);title(濾波前信號(hào)頻譜);subplot(2,2,2);plot(abs(X);axis(0,4000,0,0.03);title(濾波

56、后信號(hào)頻譜);subplot(2,2,3);plot(z1);title(濾波前信號(hào)波形);subplot(2,2,4);plot(x);title(濾波后信號(hào)波形);仿真結(jié)果如圖 2-11、圖 2-12 所示。01-15000-10000-500005000Normalized Frequency ( rad/sample)Phase (degrees)01-200-1000100Normalized Frequency ( rad/sample)Magnitude (dB)圖 2-1

57、1 IIR 數(shù)字帶通濾波器幅頻-相頻特性0500100000.81信 信 信 信 信 信 信0100020003000400000.010.020.03信 信 信 信 信 信 信0246x 104-1-0.500.51信 信 信 信 信 信 信0246x 104-1-0.500.51信 信 信 信 信 信 信圖 2-12 濾波前后信號(hào)頻譜和波形對(duì)比河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）24第 3 章 FIR 數(shù)字濾波器的設(shè)計(jì)相對(duì)于無(wú)限沖激響應(yīng)數(shù)字濾波器(IIR)，有限沖激響應(yīng)數(shù)字濾波器的特點(diǎn)是15:(l)具有精確的線性相位;(2)總是穩(wěn)定的;(3)設(shè)計(jì)方式是線性的;(4)硬件容易實(shí)現(xiàn);

58、(5)濾波器過(guò)渡過(guò)程具有有限區(qū)間;(6)相對(duì) IIR 濾波器而言，階次較高,其延遲也要比同樣性能的 IIR 濾波器大得多。線性相位 FIR 濾波器通常采用窗函數(shù)法設(shè)計(jì)。窗函數(shù)法設(shè)計(jì) FIR 濾波器的基本思想是：根據(jù)給定的濾波器技術(shù)指標(biāo)，選擇濾波器長(zhǎng)度 N 和窗函數(shù) （n） ，使其具有最窄寬度的主瓣和最小的旁瓣。其核心是從給定的頻率特性，通過(guò)加窗確定有限長(zhǎng)單位脈沖響應(yīng)序列 h(n)。工程中常用的窗函數(shù)共有 6 種，即矩形窗、巴特利特（Bartlett）窗、漢寧（Hanning）窗、海明（Hamming）窗、布萊克曼（Blackman）窗和凱塞（Kaiser）窗。3.1 典型的窗函數(shù)1、矩形窗(R

59、ectangle Window) )()(nRnwN （3-1）其頻率響應(yīng)和幅度響應(yīng)分別為：21)2/sin()2/sin()(NjjeNeW （3-2） )2/sin()2/sin()(NWR （3-3）2、三角形窗(Bartlett Window) 121,122210,12)(NnNNnNnNnnw （3-4） 其頻率響應(yīng)為： 河南科技大學(xué)畢業(yè)設(shè)計(jì)（論文）25 212)2/sin()4/sin(2)(NjjeNNeW （3-5）3、漢寧(Hanning)窗，又稱升余弦窗 )()12cos(1 21)(nRNnnwN （3-6）其頻率響應(yīng)和幅度響應(yīng)分別為：)12()12(25. 0)(5

60、. 0)()()12()12(25. 0)(5 . 0)()21(NWNWWWeWeNWNWWeWRRRajNjRRRj （3-7）4、漢明(Hamming)窗，又稱改進(jìn)的升余弦窗 )()12cos(46. 054. 0)(nRNnnwN (3-8)其幅度響應(yīng)為：)12()12(23. 0)(54. 0)(NWNWWWRRR (3-9) 5、布萊克曼(Blankman)窗，又稱二階升余弦窗 )()14cos(08. 0)12cos(5 . 042. 0)(nRNnNnnwN (3-10)其幅度響應(yīng)為：)14()14(04. 0)12()12(25. 0)(42. 0)(NWNWNWNWWWRR