Matlab的IIR數(shù)字濾波器設(shè)計

1、封面作者：PanHongliang僅供個人學(xué)習(xí)摘要在現(xiàn)代通信系統(tǒng)中，由于信號中經(jīng)?；煊懈鞣N復(fù)雜成分，所以很多信號分 析都是基于濾波器而進(jìn)行的，而數(shù)字濾波器是通過數(shù)值運(yùn)算實(shí)現(xiàn)濾波，具有處 理精度高、穩(wěn)定、靈活、不存在阻抗匹配問題，可以實(shí)現(xiàn)模擬濾波器無法實(shí)現(xiàn) 的特殊濾波功能。數(shù)字濾波器根據(jù)其沖激響應(yīng)函數(shù)的時域特性，可分為兩種， 即無限長沖激響應(yīng)(IIR)數(shù)字濾波器和有限長沖激響應(yīng)(FIR)數(shù)字濾波器。實(shí)現(xiàn)IIR濾波器的階次較低，所用的存儲單元較少，效率高，精度高，而且能夠保留 一些模擬濾波器的優(yōu)良特性，因此應(yīng)用很廣。Matlab軟件以矩陣運(yùn)算為基礎(chǔ)，把計算、可視化及程序設(shè)計有機(jī)融合到交互式工作環(huán)境

2、中，并且為數(shù)字濾波的 研究和應(yīng)用提供了一個直觀、高效、便捷的利器。尤其是Matlab中的信號處理工具箱使各個領(lǐng)域的研究人員可以直觀方便地進(jìn)行科學(xué)研究與工程應(yīng)用。本文 首先介紹了數(shù)字濾波器的概念，分類以及設(shè)計要求。接著利用MATLA函數(shù)語言編程，用信號處理圖形界面FDAToo來設(shè)計濾波器以及Sptool界面設(shè)計的方法， 并用FDATool模擬IIR數(shù)字濾波器處理信號。重點(diǎn)設(shè)計Chebyshev I型和Chebyshev II型數(shù)字低通濾波器,并介紹最優(yōu)化設(shè)計。【關(guān)鍵字】IIR濾波器FDAToolSptoolSimulinkABSTRACTIn modern communication syste

3、ms,Becauseoften mixed with various signalcomplex components,So many signal analysis is based on filters, and the digital filter isrealized through numerical computation, digital filters filter with high precision, stabilityand flexibility, dont exist, can realize the impedance matching simulating th

4、e special filtercannot achieve filter function. Digital filter according to its impulse response function andcharacteristics of the time can be divided into two kinds, namely the infinite impulseresponse(IIR) digital filter and finite impulse response (FIR digital filters). The order ofrealizing IIR

5、 filter is used, low and high efficiency less storage unit, high precision, andcan keep some simulation characteristics of filter, so it is widely used. Matlab softwarebased on matrix computation, the calculation, visualization and program design oforganic integration to interactiveenvironment for d

6、igital filter, and the research andapplication of provides an intuitive, efficient and convenient tool. Especially in the Matlabsignal processing to all areas of research toolbox personnel can easily for scientificresearch and engineering application. This paper introduces the concept of digital fil

7、ter,classification and design requirements. Then using MATLAB language programming,with functions of signal processing FDATool graphical interface design of interfacedesign and Sptool filter, and FDATool analog signal processing IIR digital filter. Keydesign Chebyshev type I and II digital Chebyshev

8、 lowpass filter, and introducesoptimization design.【Keywords】IIR Filter FDATool Sptool Simulink目錄、八、 _前言5第一章數(shù)字濾波器5第一節(jié)數(shù)字濾波器的概念5第二節(jié)數(shù)字濾波器的分類6第三節(jié)數(shù)字濾波器的設(shè)計要求7第二章IIR數(shù)字濾波器設(shè)計方法7第一節(jié)IIR數(shù)字濾波器的設(shè)計步驟8第二節(jié)用脈沖相應(yīng)不變法設(shè)計IIR數(shù)字濾波器9一、設(shè)計原理9二、脈沖響應(yīng)不變法優(yōu)缺點(diǎn)10第三節(jié)雙線性變換法設(shè)計IIR數(shù)字濾波器11一、設(shè)計原理11二、雙線性變換法優(yōu)缺點(diǎn)12第三章IIR濾波器的MATLAB計13第一節(jié)IIR數(shù)字濾波器

9、的典型設(shè)計法14第二節(jié)IIR數(shù)字濾波器的直接設(shè)計法17第三節(jié)FDATool介紹和界面設(shè)計20第四節(jié)FDATOO設(shè)計IIR數(shù)字濾波器22第五節(jié)SIMULINK仿真IIR濾波器22總結(jié)23致謝24參考文獻(xiàn)24結(jié)束語25、八前言隨著信息時代和數(shù)字世界的到來，數(shù)字信號處理已成為當(dāng)今一門極其重要 的學(xué)科和技術(shù)領(lǐng)域。目前數(shù)字信號處理在通信、語音、圖像、自動控制、雷 達(dá)、軍事、航空航天、醫(yī)療和家用電器等眾多領(lǐng)域得到了廣泛的應(yīng)用。在數(shù)字信號處理中起著重要的作用并已獲得廣泛應(yīng)用的是數(shù)字濾波器(DF，DigitalFilter)，根據(jù)其單位沖激響應(yīng)函數(shù)的時域特性可分為兩類：無限沖激響應(yīng)IIR(Infinite I

10、mpulseResponse) 濾 波 器 和 有 限 沖 激 響 應(yīng)FIR(Finite Impulse Respons&濾波器。與FIR濾波器相比，IIR的實(shí)現(xiàn)采用的是遞歸結(jié)構(gòu)， 極點(diǎn)須 在單位圓內(nèi)， 在相同設(shè)計指標(biāo)下， 實(shí)現(xiàn)IIR濾波器的階次較低，即所用的存儲 單元少，從而經(jīng)濟(jì)效率高。MATLAB是英文MATrixLABoratory(矩陣實(shí)驗(yàn)室)的 縮寫。它是美國的MathWorks公司推出的一套用于科學(xué)計算和圖形處理可視 化、 高性能語言與軟件環(huán)境。MATLAB的信號處理工具箱是專門應(yīng)用于信號處理 領(lǐng)域的專用工具箱，它的兩個基本組成就是濾波器的設(shè)計與實(shí)現(xiàn)部分以及譜分 析部分。

11、工具箱提供了豐富而簡便的設(shè)計，使原來繁瑣的程序設(shè)計簡化成函數(shù) 的調(diào)用。只要以正確的指標(biāo)參數(shù)調(diào)用相應(yīng)的濾波器設(shè)計程序或工具箱函數(shù)，便 可以得到正確的設(shè)計結(jié)果，使用非常方便。第一章 數(shù)字濾波器第一節(jié)數(shù)字濾波器的概念濾波器是指用來對輸入信號進(jìn)行濾波的硬件和軟件。數(shù)字濾波器是對數(shù)字 信號實(shí)現(xiàn)濾波的線性時不變系統(tǒng)。數(shù)字濾波器可以理解為是一個計算程序或算 法，將代表輸入信號的數(shù)字時間序列轉(zhuǎn)化為代表輸出信號的數(shù)字時間序列，并 在轉(zhuǎn)化過程中，使信號按預(yù)定的形式變化。數(shù)字濾波實(shí)質(zhì)上是一種運(yùn)算過程， 實(shí)現(xiàn)對信號的運(yùn)算處理。數(shù)字濾波器和模擬濾波器相比，因?yàn)樾盘柕男问胶蛯?shí) 現(xiàn)濾波的方法不同，數(shù)字濾波器具有比模擬濾波器

12、精度高、穩(wěn)定、體積小、重 量輕、靈活、不要求阻抗匹配等優(yōu)點(diǎn)。輸入數(shù)字信號(數(shù)字序列)通過特定的 運(yùn)算轉(zhuǎn)變?yōu)檩敵龅臄?shù)字序列，因此，數(shù)字濾波器本質(zhì)上是一個完成特定運(yùn)算的 數(shù)字計算過程，也可以理解為是一臺計算機(jī)。描述離散系統(tǒng)輸出與輸入關(guān)系的 卷積和差分方程只是給數(shù)字信號濾波器提供運(yùn)算規(guī)則，使其按照這個規(guī)則完成對輸入數(shù)據(jù)的處理時域離散系統(tǒng)的頻域特性:（式1-1） 其中、分別是數(shù)字濾波器的輸出序列和輸入序列的頻域特性（或稱為頻譜特 性）,是數(shù)字濾波器的單位取樣響應(yīng)的頻譜，又稱為數(shù)字濾波器的頻域響應(yīng)。輸 入序列的頻譜經(jīng)過濾波后,因此，只要按照輸入信號頻譜的特點(diǎn)和處理信號的目 的，適當(dāng)選擇 ， 使得濾波后的

13、滿足設(shè)計的要求，這就是數(shù)字濾波器的濾波原 理。第二節(jié)數(shù)字濾波器的分類按照不同的分類方法，數(shù)字濾波器有許多種類，但總起來可以分成兩大 類：經(jīng)典濾波器和現(xiàn)代濾波器。經(jīng)典濾波器的特點(diǎn)是其輸入信號中有用的頻率 成分和希望濾除的頻率成分占有不同的頻帶，通過一個合適的選頻濾波器濾除 干擾，得到純凈信號，達(dá)到濾波的目的。但是，如果信號和干擾的頻譜相互重 疊，則經(jīng)典濾波器不能有效地濾除干擾，最大限度地恢復(fù)信號，這時就需要現(xiàn) 代濾波器，例如維納濾波器、卡爾曼濾波器、自適應(yīng)濾波器等最佳濾波器?，F(xiàn) 代濾波器是根據(jù)隨機(jī)信號的一些統(tǒng)計特性，在某種最佳準(zhǔn)則下，最大限度地抑 制干擾，同時最大限度地回復(fù)信號，從而達(dá)到最佳濾波

14、的目的。經(jīng)典數(shù)字濾波器從濾波特性上分類，可以分為：低通濾波器、高通濾波 器、帶通濾波器、帶阻濾波器。圖1-1各種理想濾波器的幅頻特性 數(shù)字濾波器根據(jù)其沖激響應(yīng)函數(shù)的時域特性，可分為兩種，即無限長沖激 響應(yīng)（IIR）數(shù)字濾波器和有限長沖激響應(yīng)（FIR）數(shù)字濾波器。IIR數(shù)字濾波器的 特征是，具有無限持續(xù)時間沖激響應(yīng)，需要用遞歸模型來實(shí)現(xiàn)，其差分方程為：（式1-2）系統(tǒng)函數(shù)為：（式1-3）設(shè)計IIR濾波器的任務(wù)就是尋求一個物理上可實(shí)現(xiàn)的系統(tǒng)函數(shù)H（z），使其頻 率響應(yīng)H（z）滿足所希望得到的頻域指標(biāo)，即符合給定的通帶截止頻率、阻帶截 止頻率、通帶衰減系數(shù)和阻帶衰減系數(shù)。第三節(jié)數(shù)字濾波器的設(shè)計要求濾

15、波器的指標(biāo)常常在頻域給出。數(shù)字濾波器的頻響特性函數(shù)一般為復(fù)函 數(shù)，所以通常表示為：（式1-4）其中，|稱為幅頻特性函數(shù)，稱為相頻特性函數(shù)。幅頻特性表示信號 通過該濾波器后各頻率成分的衰減情況，而相頻特性反映各頻率通過濾波器后 在時間上的延時情況。一般IIR數(shù)字濾波器，通常只用幅頻響應(yīng)函數(shù)|來描述 設(shè)計指標(biāo)，相頻特性一般不作要求。IIR濾波器指標(biāo)參數(shù)如下圖所示。圖中，3 p和3s分別為通帶邊界頻率和 阻帶邊界頻率；S1和S2分別為通帶波紋和阻帶波紋；允許的衰減一般用dB數(shù)表示，通帶內(nèi)所允許的最大衰減（dB）和阻帶內(nèi)允許的最小衰減（dB）分別 為ap和口aS表示：（式1-5）（式1-6）一般要求：

16、當(dāng)時，；當(dāng)時，。圖1-2低通濾波器的技術(shù)要求第二章IIR數(shù)字濾波器設(shè)計方法IIR數(shù)字濾波器是一種離散時間系統(tǒng)，其系統(tǒng)函數(shù)為（式2-1）假設(shè)MN,當(dāng)e N寸,系統(tǒng)函數(shù)可以看作一個IIR的子系統(tǒng)和一個（M-N）的FIR子系統(tǒng)的級聯(lián)。IIR數(shù)字濾波器的設(shè)計實(shí)際上是求解濾波器的系數(shù)和，它是 數(shù)學(xué)上的一種逼近問題，即在規(guī)定意義上（通常采用最小均方誤差準(zhǔn)則）去逼 近系統(tǒng)的特性。如果在S平面上去逼近，就得到模擬濾波器；如果在z平面上去逼近，就得到數(shù)字濾波器。第一節(jié)IIR數(shù)字濾波器的設(shè)計步驟IIR數(shù)字濾波器的設(shè)計一般有兩種方法：一個是借助模擬濾波器的設(shè)計方法 進(jìn)行。其設(shè)計步驟是，先設(shè)計模擬濾波器，再按照某種方

17、法轉(zhuǎn)換成數(shù)字濾波 器。這種方法比較容易一些，因?yàn)槟M濾波器的設(shè)計方法已經(jīng)非常成熟，不僅 有完整的設(shè)計公式，還有完善的圖表供查閱；另外一種直接在頻率或者時域內(nèi) 進(jìn)行，由于需要解聯(lián)立方程，設(shè)計時需要計算機(jī)做輔助設(shè)計。其設(shè)計步驟是： 先設(shè)計過渡模擬濾波器得到系統(tǒng)函數(shù)，然后將按某種方法轉(zhuǎn)換成數(shù)字濾波器的 系統(tǒng)函數(shù)。這是因?yàn)槟M濾波器的設(shè)計方法已經(jīng)很成熟，不僅有完整設(shè)計公 式，還有完善的圖表和曲線供查閱；另外，還有一些典型的優(yōu)良濾波器類型可 供我們使用。為了保證轉(zhuǎn)換后的穩(wěn)定且滿足技術(shù)指標(biāo)要求，對轉(zhuǎn)換關(guān)系提出兩點(diǎn)要求:(1)因果穩(wěn)定的模擬濾波器轉(zhuǎn)換成數(shù)字濾波器，仍是因果穩(wěn)定的。(2)數(shù)字濾波器的頻率相應(yīng)模

18、仿模擬濾波器的頻響特性，s平面的虛軸映射為z平面的單位圓，相應(yīng)的頻率之間呈線性關(guān)系。利用模擬濾波器成熟的理論設(shè)計IIR數(shù)字濾波器的過程是：(1)確定數(shù)字低通濾波器的技術(shù)指標(biāo)：通帶邊界頻率、通帶最大衰減、阻帶 截止頻率、阻帶最小衰減。(2)將數(shù)字低通濾波器的技術(shù)指標(biāo)轉(zhuǎn)換成相應(yīng)的模擬低通濾波器的技術(shù)指 標(biāo)。(3)按照模擬低通濾波器的技術(shù)指標(biāo)設(shè)計過渡模擬低通濾波器。(4)用所選的轉(zhuǎn)換方法，將模擬濾波器轉(zhuǎn)換成數(shù)字低通濾波器系統(tǒng)函數(shù)。IIR數(shù)字濾波器的設(shè)計流程圖如下：圖2-1IIR數(shù)字濾波器的設(shè)計步驟流程圖成熟的模擬濾波器設(shè)計方法主要有脈沖響應(yīng)不變法和雙線性變換法。第二節(jié)用脈沖相應(yīng)不變法設(shè)計IIR數(shù)字濾

19、波器一、設(shè)計原理利用模擬濾波器來設(shè)計數(shù)字濾波器，也就是使數(shù)字濾波器能模仿模擬濾波 器的特性，這種模仿可以從不同的角度出發(fā)。脈沖響應(yīng)不變法是從濾波器的脈 沖響應(yīng)出發(fā)，使數(shù)字濾波器的單位脈沖響應(yīng)序列h(n)模仿模擬濾波器的沖激響應(yīng)ha(t)，即將ha(t)進(jìn)行等間隔采樣，使h(n)正好等于ha(t)的采樣值，滿足h(n)=ha(nT)式中,T是采樣周期。如果令HL(s)是ha(t)的拉普拉斯變換，H(z)為h(n)的Z變換，利用采樣序 列的Z變換與模擬信號的拉普拉斯變換的關(guān)系得(式2-2)則可看出，脈沖響應(yīng)不變法將模擬濾波器的S平面變換成數(shù)字濾波器的Z平面，這個從s到z的變換z=esT是從S平面變

20、換到Z平面的標(biāo)準(zhǔn)變換關(guān)系式。圖2-2脈沖響應(yīng)不變法的映射關(guān)系由(2-2)式，數(shù)字濾波器的頻率響應(yīng)和模擬濾波器的頻率響應(yīng)間的關(guān)系為(式2-3)這就是說，數(shù)字濾波器的頻率響應(yīng)是模擬濾波器頻率響應(yīng)的周期延拓。正 如采樣定理所討論的，只有當(dāng)模擬濾波器的頻率響應(yīng)是限帶的，且?guī)抻谡郫B頻率以內(nèi)時，即（式2-4）才能使數(shù)字濾波器的頻率響應(yīng)在折疊頻率以內(nèi)重現(xiàn)模擬濾波器的頻率響應(yīng)，而不產(chǎn)生混疊失真，即（式2-5）但是，任何一個實(shí)際的模擬濾波器頻率響應(yīng)都不是嚴(yán)格限帶的，變換后就 會產(chǎn)生周期延拓分量的頻譜交疊，即產(chǎn)生頻率響應(yīng)的混疊失真。這時數(shù)字濾波 器的頻響就不同于原模擬濾波器的頻響，而帶有一定的失真。當(dāng)模擬濾波器

21、的 頻率響應(yīng)在折疊頻率以上處衰減越大、越快時，變換后頻率響應(yīng)混疊失真就越 小。 這時， 采用脈沖響應(yīng)不變法設(shè)計的數(shù)字濾波器才能得到良好的效果。圖2-3脈沖響應(yīng)不變法中的頻響混疊現(xiàn)象對某一模擬濾波器的單位沖激響應(yīng)ha（t）進(jìn)行采樣，采樣頻率為fs,若使fs增加，即令采樣時間間隔（T=1/fs）減小，則系統(tǒng)頻率響應(yīng)各周期延拓分量之 間相距更遠(yuǎn)，因而可減小頻率響應(yīng)的混疊效應(yīng)。二、脈沖響應(yīng)不變法優(yōu)缺點(diǎn)從以上討論可以看出，脈沖響應(yīng)不變法使得數(shù)字濾波器的單位脈沖響應(yīng)完 全模仿模擬濾波器的單位沖激響應(yīng)，也就是時域逼近良好，而且模擬頻率Q和數(shù)字頻率3之間呈線性關(guān)系3=QT。因而，一個線性相位的模擬濾波器（例如

22、 貝塞爾濾波器）通過脈沖響應(yīng)不變法得到的仍然是一個線性相位的數(shù)字濾波 器。脈沖響應(yīng)不變法的最大缺點(diǎn)是有頻率響應(yīng)的混疊效應(yīng)。所以，脈沖響應(yīng)不變法只適用于限帶的模擬濾波器（例如，衰減特性很好的低通或帶通濾波器）， 而且咼頻衰減越快，混疊效應(yīng)越小。至于咼通和帶阻濾波器，由于它們在咼頻 部分不衰減，因此將完全混淆在低頻響應(yīng)中。如果要對高通和帶阻濾波器采用 脈沖響應(yīng)不變法，就必須先對高通和帶阻濾波器加一保護(hù)濾波器，濾掉高于折 疊頻率以上的頻率，然后再使用脈沖響應(yīng)不變法轉(zhuǎn)換為數(shù)字濾波器。當(dāng)然這樣 會進(jìn)一步增加設(shè)計復(fù)雜性和濾波器的階數(shù)。第三節(jié)雙線性變換法設(shè)計IIR數(shù)字濾波器一、設(shè)計原理脈沖響應(yīng)不變法的主要缺

23、點(diǎn)是產(chǎn)生頻率響應(yīng)的混疊失真。這是因?yàn)閺腟平 面到Z平面是多值的映射關(guān)系所造成的。為了克服這一缺點(diǎn)，可以采用非線性 頻率壓縮方法，將整個頻率軸上的頻率范圍壓縮到-n/Tn/T之間，再用z=esT轉(zhuǎn)換到Z平面上。也就是說，第一步先將整個S平面壓縮映射到S平面的-n/Tn/T一條橫帶里；第二步再通過標(biāo)準(zhǔn)變換關(guān)系z=e3T將此橫帶變換到整個Z平面上去。這樣就使S平面與Z平面建立了一一對應(yīng)的單值關(guān)系，消除 了多值變換性，也就消除了頻譜混疊現(xiàn)象，映射關(guān)系如圖2-3所示。圖2-4雙線性變換的映射關(guān)系為了將S平面的整個虛軸jQ壓縮到S1平面jQ1軸上的-n/T到n/T段 上，可以通過以下的正切變換實(shí)現(xiàn)（式2-

24、6）式中,T仍是采樣間隔。當(dāng)Q1由-n/T經(jīng)過0變化到n/T時，Q由-x經(jīng)過o變化到+*，也即映射 了整個jQ軸。將式（2-6）寫成（式2-7）將此關(guān)系解讀延拓到整個S平面和S1平面，令jQ=s，jQ1=s1,則得（式2-8）再將S1平面通過以下標(biāo)準(zhǔn)變換關(guān)系映射到Z平面（式2-9）從而得到S平面和Z平面的單值映射關(guān)系為：（式2-10）（式2-11）式（2-10）與式（2-11）是S平面與Z平面之間的單值映射關(guān)系，這種變 換都是兩個線性函數(shù)之比，因此稱為雙線性變換式（2-6）與式（2-10）的雙線性變換符合映射變換應(yīng)滿足的兩點(diǎn)要求。 首先,把z=ejS可得式2-12)即S平面的虛軸映射到Z平面的

25、單位圓。其次，將S=c+jQ代入式（2-12），得（式2-13）因此（式2-14）由此看出，當(dāng)（70時，|z|0時，|z|1。也就是說，S平面的左 半平面映射到Z平面的單位圓內(nèi)，S平面的右半平面映射到Z平面的單位圓 外，S平面的虛軸映射到Z平面的單位圓上。 因此， 穩(wěn)定的模擬濾波器經(jīng)雙線 性變換后所得的數(shù)字濾波器也一定是穩(wěn)定的。二、雙線性變換法優(yōu)缺點(diǎn)雙線性變換法與脈沖響應(yīng)不變法相比，其主要的優(yōu)點(diǎn)是避免了頻率響應(yīng)的 混疊現(xiàn)象。 這是因?yàn)镾平面與Z平面是單值的一一對應(yīng)關(guān)系。S平面整個jQ軸 單值地對應(yīng)于Z平面單位圓一周，即頻率軸是單值變換關(guān)系。這個關(guān)系如式（2-10）所示，重寫如下：（式2-15）

26、上式表明，S平面上Q與Z平面的3成非線性的正切關(guān)系，如圖2-4所示。由圖2-4看出，在零頻率附近，模擬角頻率Q與數(shù)字頻率3之間的變換關(guān)系接近于線性關(guān)系；但當(dāng)Q進(jìn)一步增加時，3增長得越來越慢，最后當(dāng)時，3終止在折疊頻率3=冗處，因而雙線性變換就不會出現(xiàn)由于高頻部分超過 折疊頻率而混淆到低頻部分去的現(xiàn)象，從而消除了頻率混疊現(xiàn)象。圖2-5雙線性變換法的頻率變換關(guān)系 但是雙線性變換的這個特點(diǎn)是靠頻率的嚴(yán)重非線性關(guān)系而得到的，如式（2-12）及圖2-4所示。由于這種頻率之間的非線性變換關(guān)系，就產(chǎn)生了新的 問題。首先，一個線性相位的模擬濾波器經(jīng)雙線性變換后得到非線性相位的數(shù) 字濾波器，不再保持原有的線性相

27、位了；其次，這種非線性關(guān)系要求模擬濾波 器的幅頻響應(yīng)必須是分段常數(shù)型的，即某一頻率段的幅頻響應(yīng)近似等于某一常 數(shù) （這正是一般典型的低通、 高通、 帶通、 帶阻型濾波器的響應(yīng)特性） ，不然 變換所產(chǎn)生的數(shù)字濾波器幅頻響應(yīng)相對于原模擬濾波器的幅頻響應(yīng)會有畸變， 如圖2-5所示。圖2-6雙線性變換法幅度和相位特性的非線性映射 對于分段常數(shù)的濾波器，雙線性變換后，仍得到幅頻特性為分段常數(shù)的濾 波器，但是各個分段邊緣的臨界頻率點(diǎn)產(chǎn)生了畸變，這種頻率的畸變，可以通 過頻率的預(yù)畸變來加以校正。也就是將臨界模擬頻率事先加以畸變，然后經(jīng)變 換后正好映射到所需要的數(shù)字頻率上。第三章IIR濾波器的MATLAB設(shè)計

28、MATLAB是矩陣實(shí)驗(yàn)室(Matrix Laboratory)之意。除具備卓越的數(shù)值計 算能力外，它還提供了專業(yè)水平的符號計算，文字處理，可視化建模仿真和實(shí) 時控制等功能。MATLAB的基本數(shù)據(jù)單位是矩陣，它的指令表達(dá)式與數(shù)學(xué),工程 中常用的形式十分相似,故用MATLAB來解算問題要比用C,FORTRA等語言完相 同的事情簡捷得多在新的版本中也加入了對C,FORTRAN,c+ JAVA的支持.可 以直接調(diào)用,用戶也可以將自己編寫的實(shí)用程序?qū)氲組ATLAB函數(shù)庫中方便自己以后調(diào)用，此外許多的MATLA愛好者都編寫了一些經(jīng)典的程序，用戶可以直 接進(jìn)行下載就可以用，非常的方便。模擬濾波器的理論和設(shè)

29、計方法已發(fā)展得相當(dāng)成熟，且有一些典型的模擬濾 波器供我們選擇，如巴特沃斯(Butterworth )濾波器、切比雪夫(Chebyshev)濾 波器、橢圓(Cauer )濾波器、貝塞爾(Bessel )濾波器等，這些典型的濾波器各 有特點(diǎn)。用MATLA進(jìn)行數(shù)字濾波器的設(shè)計的步驟如下：(1) 將設(shè)計指標(biāo)歸一化處理。如果采用雙線性變換法，還需進(jìn)行預(yù)畸變。(2)根據(jù)歸一化頻率，確定最小階數(shù)N和頻率參數(shù)Wn可供選用的階數(shù)擇函 數(shù)有：buttord，cheblord，cheb2ord，ellipord等。(3) 運(yùn)用最小階數(shù)N設(shè)計模擬低通濾波器原型。模擬低通濾波器的創(chuàng)建函數(shù)有:buttap,cheblap

30、,cheb2ap,ellipap和besselap，這些函數(shù)輸出的是零極點(diǎn)式形式，還要用zp2tf函數(shù)轉(zhuǎn)換成分子分母多項(xiàng)式形式。如果想 根 據(jù) 最 小 階數(shù) 直 接 設(shè) 計 模 擬 低 通 濾 波 器 原 型 ， 可 用Wn設(shè)為1。(4)根據(jù)第2步的頻率參數(shù)Wn模擬低通濾波原型轉(zhuǎn)換模擬低通、高通、帶 通、帶阻濾波器，可用函數(shù)分別是：lp21p,lp2hp,lp2bp,lp2bs。(5) 運(yùn)用脈沖響應(yīng)不變法或雙線性變法把模濾波器轉(zhuǎn)數(shù)字濾波器，調(diào)用的函 數(shù)是impinvar和bilinear。脈沖響應(yīng)不變法適用于采樣頻率大于4倍 截止頻率的銳截止低通帶通濾波器，而雙線性變換法適合于相位特性要 求不

31、高的各型濾波器。(6)根據(jù)輸出的分子分母系數(shù)，調(diào)用函數(shù)buttord計算N和wc,有系數(shù)向量 可以寫出數(shù)字濾波器系統(tǒng)函數(shù)乙再用freqz函數(shù)驗(yàn)證設(shè)計結(jié)果。第一節(jié)IIR數(shù)字濾波器的典型設(shè)計法設(shè)計巴特沃斯數(shù)字低通濾波器和橢圓數(shù)字低通濾波器,要求通帶邊界頻fp=2.1kHZ，通帶最大衰減Rp=0.5dB;阻帶邊界頻率fs=8kHZ,阻帶最小衰減Rs=30dB采樣頻率為Fs=20kHZ低通巴特沃斯濾波器設(shè)計步驟如下：(1) 確定階數(shù)N。取N=4(2) 求極點(diǎn)? ? ? ? 歸一化低通原型系統(tǒng)函數(shù)為由N=4直接查表得到：極點(diǎn)：歸一化低通濾波器系統(tǒng)函數(shù)為式中，0.0000,0.0999,0.1914,0.

32、0252(3) 將去歸一化最終得到butter,chebyl,cheby2,ellip,bessel等函數(shù)，只是注意要將函數(shù)中的通過計算可以總結(jié)出過程太麻煩，而且容易出錯，結(jié)果不直觀。下面用M程序設(shè)計來實(shí)現(xiàn)用脈沖響應(yīng)不變法設(shè)計的巴特沃斯數(shù)字低通濾波器的M程序如下:fp=2100。fs=8000。Fs=20000。Rp=0.5。Rs=30。T=1/Fs。%設(shè)計指標(biāo)W1p=fp/Fs*2。W1s=fs/Fs*2。求歸一化頻率N,Wn=buttord(W1p,W1s,Rp,Rs,s)。%確定butterworth的最小介數(shù)N和頻率參數(shù)Wn設(shè)計模擬低通原型的零極點(diǎn)增益參數(shù)將零極點(diǎn)增益轉(zhuǎn)換成分子分母參數(shù)

33、。%將低通原型轉(zhuǎn)換為模擬低通。%用脈沖響應(yīng)不變法進(jìn)行模數(shù)變換sys=tf(bz,az,T)。%H,W=freqz(bz,az,512,Fs)。subplot(2,1,1)。plot(W,20*log10(abs(H)。gridon。% xlabel(頻率/Hz)。ylabel(振幅/dB)。subplot(2,1,2)。plot(W,abs(H)。grid on。xlabel(頻率/Hz)。ylabel(振幅/H)。運(yùn)行后的波形如下:圖3-1典型濾波器在Matlab上運(yùn)行波形運(yùn)行結(jié)果：N =4 bz = 0.0000 0.0999 0.1914 0.0252az=1.0000 -1.4336

34、1.0984 -0.4115 0.0627可以得出：只需編程，結(jié)果非常直觀。 雙線性變換法設(shè)計步驟如下：z,p,k=buttap(N)。%bp,ap=zp2tf(z,p,k)。%bs,as=lp2lp(bp,ap,Wn*pi*Fs)bz,az=impinvar(bs,as,Fs)給出傳輸函數(shù)H(Z)%生成頻率響應(yīng)參數(shù)%繪制幅頻響應(yīng)加坐標(biāo)網(wǎng)格（1）首先寫出該濾波器的系統(tǒng)函數(shù)（2）利用雙線性變換法轉(zhuǎn)換，數(shù)字濾波器的系統(tǒng)函數(shù)為用雙線性變換法設(shè)計橢圓數(shù)字低通濾波器的M程序如下:fs=20000。wp=2*pi*2100/fs。ws=2*pi*8000/fs。Rp=0.5。Rs=30。Ts=1/fs。W

35、p=2/Ts*tan(wp/2)。Ws=2/Ts*tan(ws/2)。%按頻率轉(zhuǎn)換公式進(jìn)行轉(zhuǎn)換N,Wn=ellipord(Wp,Ws,Rp,Rs,s)。%計算模擬濾波器的最小階數(shù)z,p,k=ellipap(N,Rp,Rs)。%設(shè)計模擬原型濾波器Bap,Aap=zp2tf(z,p,k)。%零點(diǎn)極點(diǎn)增益形式轉(zhuǎn)換為傳遞函數(shù)形式b,a=lp2lp(Bap,Aap,Wn)。%低通轉(zhuǎn)換為低通濾波器的頻率轉(zhuǎn)化bz,az=bilinear(b,a,fs)。%運(yùn)用雙線性變換法得到數(shù)字濾波器傳遞函數(shù)H,f=freqz(bz,az,512,fs)subplot(2,1,1)。plot(f,20*log10(abs(

36、H)title(N=2頻率響應(yīng))grid on。xlabel(頻率/Hz)。ylabel(振幅/dB)。subplot(2,1,2)。plot(f,abs(H)。grid on。xlabel(頻率/Hz)。ylabel(振幅/H)。圖3-2雙線性變換法設(shè)計的橢圓數(shù)字低通濾波器運(yùn)行波形 運(yùn)行結(jié)果：N=2bz= 0.1213 0.1662 0.1213az= 1.0000 -0.9889 0.4218第二節(jié)IIR數(shù)字濾波器的直接設(shè)計法除了典型設(shè)計以外，MATLAB言號處理工具箱提供了幾個直接設(shè)計IIR數(shù)字濾波器的函數(shù)，直接調(diào)用就可以設(shè)計濾波器，這為設(shè)計通用濾波器提供了方 便。主要有以下幾種方法：1

37、.零極點(diǎn)累試法 這種方法也稱為零極點(diǎn)累試法。在確定零極點(diǎn)位置時要注意：(1)極點(diǎn)必須位于z平面單位圓內(nèi)，保證數(shù)字濾波器因果穩(wěn)定；(2)復(fù)數(shù)零極點(diǎn)必須共軛成對，保證系統(tǒng)函數(shù)有理式的系數(shù)是實(shí)的。2.在頻域利用幅度平方誤差最小法直接設(shè)計IIR數(shù)字濾波器3.在時域直接設(shè)計IIR數(shù)字濾波器設(shè)計Butterworth濾波器用函數(shù)butter()，可以設(shè)計低通、高通、帶通和 帶阻的數(shù)字和模擬濾波器，其特性是通帶內(nèi)的幅度響應(yīng)最大限度的平滑，但損 失了截止頻率處的下降斜度。設(shè)計Chebyshev I型濾波器用函數(shù)chebyl()?？梢栽O(shè)計低通、高通、帶通 和帶阻的數(shù)字和模擬Chebyshev I型濾波器，其通帶

38、內(nèi)為等波紋，阻帶內(nèi)為單 調(diào)。Chebyshev I型濾波器的下降斜度比II型大，但其代價目是通帶內(nèi)波紋較 大。設(shè)計Chebyshev II型濾波器用函數(shù)cheby2()?？梢栽O(shè)計低通、高通、帶 通和帶阻的數(shù)字和模擬Chebyshev II型濾波器，其通帶內(nèi)為單調(diào)，阻帶內(nèi)等波 紋。Chebyshev II型濾波器的下降斜度比I型小，但其阻帶內(nèi)波紋較大。設(shè)計橢圓濾波器用函數(shù)ellip()，與chebyl, cheby2類似，可以設(shè)計低 通、高通、帶通和帶阻的數(shù)字和模擬濾波器。與Butterworth和chebyshev濾 波器相比，ellip函數(shù)可以得到下降斜度更大的濾波器，得通帶和阻帶均為等 波

39、紋。一般情況下，橢圓濾波器能以最低的階實(shí)現(xiàn)指定的性能指標(biāo)。在使用各類濾波器函數(shù)時應(yīng)當(dāng)注意以下重點(diǎn):(1)、階數(shù)和固有頻率的選擇:N,Wn=buttord(Wp,Ws,Rp,Rs)可得到符合要 求性質(zhì)的濾波器的最小階數(shù)N以及數(shù)字Butterworth濾波器的固有頻率Wn即3dB)。設(shè)計的要求是在通帶內(nèi)的衰減不超過Rp,在阻帶內(nèi)的衰減不小于Rs,通帶和阻帶有截止頻率分別是Wp, Ws,它們是歸一化的頻率，范圍是0, 1，對應(yīng)n弧度。(2)、有關(guān)濾波器設(shè)計當(dāng)中的頻率歸一化問題:信號處理工具箱中經(jīng)常使用的 頻率是Nyquist頻率，它被定義為采樣頻率的一半，在濾波器的階數(shù)選擇和設(shè) 計中的截止頻率均使用

40、Nyquist頻率進(jìn)行歸一化處理。例如對于一個采樣頻率為1000 Hz的系統(tǒng)，400Hz的歸一化即為400/500=0.8。歸一化頻率的范圍在0, 1之間。如果要將歸一化頻率轉(zhuǎn)換為角頻率，則將歸一化頻率乘以n。如果要將歸一化頻率轉(zhuǎn)換為Hz,則將歸一化頻率乘以采樣頻率的一半。(3)、設(shè)計一個N階的低通Butterworth濾波器使用函數(shù)B,A=butter(N,Wn),返回濾波器系數(shù)矩陣B,A。其中固有頻率Wn必須是歸一化頻率。它的最 大值是采樣頻率的一半。Fs缺省時默認(rèn)為2Hz。如果Wn=Wl,W2是一個兩元素 的向量，則函數(shù)將設(shè)計出一個2N階的帶通濾波器，通帶為W1,W2。設(shè)計Chebysh

41、ev I型和Chebyshev II型數(shù)字低通濾波器，要求通帶邊界頻 率fp=2.1kHZ，通帶最大衰減Rp=0.5dB;阻帶邊界頻率fs=8kHZ，阻帶最小衰 減Rs=30dB采樣頻率為Fs=20kHZChebyshev I型的M程序如下：Fs=20000。%抽樣頻率20KHzFlp=2100。運(yùn)行后的波形如下：圖3-3Chebyshev I型運(yùn)行波形 運(yùn)行結(jié)果：N=2 b= 0.1007 0.2014 0.1007a= 1.0000 -0.9872 0.4140Chebyshev II型M程序如下:Fs=20000。Flp=2100。Fls=8000。Wp=2*Flp/Fs。Ws=2*Fl

42、s/Fs。Rp=0.5。Rs=30。N,Wn=cheb2ord(Wp,Ws,Rp,Rs)。%返回最小階數(shù)和截止頻率Fls=8000Wp=2*Flp/Fs%歸一化的通帶截止頻率Ws=2*Fls/Fs%歸一化的阻帶截止頻率Rp=0.5。%通帶最大衰減（單位：dB）Rs=30。%阻帶最小衰減（單位：dB）N,Wn=cheb1ord（Wp,Ws,Rp,Rs）。%返回最小階數(shù)和截止頻率b,a=cheby1(N,Rp,Wn)。%返回H（z）的分子分母系數(shù)hw,w=freqz(b,a)。subplot(2,1,1)。plot(w/pi,20*log10(abs(hw)。grid on。xlabel(s /

43、nylabel(幅度(dB)I型幅頻響應(yīng)）。subplot(2,1,2)。plot(w/pi,abs(hw)。gridon。xlabel(s/nylabel(幅度(H)。%抽樣頻率20KHz%歸一化的通帶截止頻率%歸一化的阻帶截止頻率%通帶最大衰減（單位:dB）b,a=cheby2(N,Rs,Wn)。%返回H(z)的分子分母系數(shù)hw,w=freqz(b,a)。subplot(2,1,1)。plot(w/pi,20*log10(abs(hw)。grid on。xlabel(S /n)。ylabel(幅度(dB)title(切比雪夫II型幅頻響應(yīng))。subplot(2,1,2)。plot(w/pi

44、,abs(hw)。grid on。xlabel(S /n)。ylabel(幅度(H)。運(yùn)行后的波形：圖3-4 Chebyshev II型運(yùn)行波形運(yùn)行結(jié)果：N=2b= 0.2357 0.42410.2357a= 1.0000 -0.2996 0.1950從頻率響應(yīng)圖中可以看出:巴特沃斯濾波器具有單調(diào)下降的幅頻特性，通帶 內(nèi)平滑。切比雪I型濾波器的幅頻特性在通帶內(nèi)有波動，阻帶內(nèi)單調(diào)。chebyshev II型濾波器的幅頻特性在阻帶內(nèi)有波動，通帶內(nèi)單調(diào)。橢圓濾波器 的選擇性相對前三種是最好的，下降斜度比較大，通帶和阻帶內(nèi)均為等波紋， 同樣的性能指標(biāo)，橢圓濾波器可以最低的階數(shù)來實(shí)現(xiàn)。這樣根據(jù)不同的要求

45、可 以選用不同類波器。第三節(jié)FDATool介紹和界面設(shè)計FDATool(Filter Design&Analysis Tool)是MATLAB言號處理工具箱里專用 的濾波器設(shè)計分析工具，MATLAB 6.0以上的版本還專門增加了濾波器設(shè)計工具 箱(Filter DesignToolbox)。FDATool可以設(shè)計幾乎所有的常規(guī)濾波器，包括FIR和IIR的各種設(shè)計方法。它操作簡單，方便靈活。FDAToo。界面總共分兩大部分，一部分是Design Filter。在界面的下半部，用來設(shè)置濾波器的設(shè)計參數(shù)；另一部分則是特性區(qū)，在界面的上半部分， 用來顯示濾波器的各種特性。Design Filt

46、er部分主要分為：Filter Type（濾 波器類型）選項(xiàng)，包括Lowpass（低通）、Highpass（高通）、Bandpass（帶通）、Bandstop（帶阻）和特殊的FIR濾波器。Design Method（設(shè)計方法）選項(xiàng)，包括IIR濾波器的Butterwotth（巴特沃 思）法、Chebyshev Type I（切比雪夫I型）法、Chebyshev Type II（切比雪夫II型）法、Elliptic（橢圓濾波器）法和FIR濾波器的Equiripple法、Least- Squares（最小乘方）法、Window（窗函數(shù)）法。Filter Order（濾波器階數(shù)）選項(xiàng)，定義濾波器的階數(shù)

47、，包括SpecifyOrder（指定階數(shù)）和Minimum Order（最小階數(shù)）。在Specify Order中填入所要 設(shè)計的濾波器的階數(shù)（N階濾波器，Specify Order=N-1）。如果選擇Minimum Order，則MATLA根據(jù)所選擇的濾波器類型自動使用最小階數(shù)。Frequency Specifications選項(xiàng)，可以詳細(xì)定義頻帶的各參數(shù)，包括采樣 頻率和頻帶的截止頻率。它的具體選項(xiàng)由Filter Type選項(xiàng)和Design Method選項(xiàng)決定。例如Bandpass（帶通）濾波器需要定義Fstop1 （下阻帶截止頻率）、Fpass1 （通帶下限截止頻率）、Fpass2（通

48、帶上限截止頻率）、Fstop2（上阻帶截 止頻率），而Lowpass（低通）濾波器只需要定義Fstop1、Fpass1。采用窗函數(shù)設(shè) 計濾波器時，由于過渡帶是由窗函數(shù)的類型和階數(shù)所決定，所以只需定義通帶 截止頻率，而不必定義阻帶參數(shù)。Magnitude Specifications選項(xiàng)，可以定義幅值衰減的情況。例如設(shè)計帶 通濾波器時，可以定義Wstop1（頻率Fstop1處的幅值衰減）、Wpass （通帶范圍 內(nèi)的幅值衰減）、Wstop2（頻率Fstop2處的幅值衰減）。當(dāng)采用窗函數(shù)設(shè)計時， 通帶截止頻率處的幅值衰減固定為6db，所以不必定義。Window Specifications選項(xiàng)，當(dāng)

49、選取采用窗函數(shù)設(shè)計，該選項(xiàng)可定義， 它包含了各種窗函數(shù)?；贛ATLABS數(shù)的濾波器設(shè)計完成后，需要對已設(shè)計的濾波器的頻率響應(yīng)要進(jìn)行校核。要得到幅頻、相頻響應(yīng)特性，運(yùn)算量也是很大的。而利用MATLAB強(qiáng)大的信號處理界面工具進(jìn)行計算機(jī)輔助設(shè)計，可以快速有效地設(shè)計數(shù)字濾波 器，大大地簡化了計算量。第四節(jié)FDATOOL設(shè)計IIR數(shù)字濾波器我們將以一個IIR濾波器的設(shè)計實(shí)例來具體說明使用MATLAB工具箱的方 便。要求設(shè)計Chebyshev1型低通濾波器滿足指標(biāo)，首先在Response Type中選 擇Bandpass高通濾波器，然后在下面的Desigh Method中選擇IIR類型，并且 指定Fil

50、ter Order項(xiàng)中的階數(shù)SpecifyOrder=10，由于是設(shè)計Chebyshev濾波 器，其下面Option就不必選擇。然后在FrequencySpecifications中選擇Unit為Hz,給出采樣頻率Fs=1000,;最后在Magnitude Specifications中選擇Unit為dB,Apass=1。設(shè)置完成后點(diǎn)擊Design Filter即可得到所設(shè)計的IIR濾波器。通過菜單選項(xiàng)Analysis可以在特性區(qū)看到所設(shè)計的幅頻響應(yīng)、相頻響 應(yīng)、沖擊響應(yīng)和零極點(diǎn)配置等特性。設(shè)計完成后將結(jié)果保存,封裝為filterl.fda文件,下面是運(yùn)用FDATool設(shè)計濾波圖3-5FDA設(shè)

51、計主界面圖3-6相頻特性圖3-7沖擊響應(yīng)特性從以上這些我們可以清晰明了的看到設(shè)計的該濾波器的幅頻,相頻和沖擊 響應(yīng)特性：我們能夠看到，圖3-3中Chebshev1型濾波器具有等紋波的通頻帶、 阻帶衰減的特性，總之，使用FDATOOLT具包設(shè)計和分析濾波器，是非常方便 易行的,而且交互性良好,不需要極其復(fù)雜的程序編制就可以實(shí)現(xiàn)。而且我們 在最后的仿真模型設(shè)計中將載入這一模塊進(jìn)行信號濾波。第五節(jié)SIMULINK仿真IIR濾波器SIMULIN!是一個進(jìn)行動態(tài)系統(tǒng)建模、仿真和綜合分析的集成軟件包。它可 以處理的系統(tǒng)包括：線性、非線性系統(tǒng);離散、連續(xù)及混合系統(tǒng);單任務(wù)、多 任務(wù)離散事件系統(tǒng)。本章通過調(diào)用

52、Simulink中的功能模塊構(gòu)成IIR數(shù)字濾波器的仿真框圖，在仿 真過程中，可以雙擊各功能模塊，隨時改變參數(shù)，獲得不同狀態(tài)下的仿真結(jié) 果。首先建立信號傳輸仿真模塊，先打開SIMULINK工具箱，并且建立一個Model，在這個空白Model中進(jìn)行環(huán)節(jié)庫及框圖的建立，在SIMULINK菜單下找到Source，雙擊Source圖標(biāo)，將正弦信號源和常量信號源拉到Model中，然后分別找到DSP離散正弦信號源,乘法 器和加法器以及示波器，然后把 需要用來設(shè)計IIR數(shù)字濾波器的模塊FDATOOL都拉到Model中，把環(huán)節(jié) 都布好后，把各環(huán)節(jié)的端口按框圖連接起來。下面就來仿真3.5節(jié)中信號處理的模型，以框圖

53、的方式直觀地現(xiàn)實(shí)信號處理 的過程，信號源為f= sin(2*pi*t*40)+sin(3*pi*t*50)+sin(2*pi*t*100)，在DSP離 散信號模塊中設(shè)定參數(shù)分別為40HZ,75HZ,100HZ,在simulink環(huán)境中，將3.5中封裝好的filterl.fda濾波器文件載入FDATOOL中，生成仿真框圖和濾波效果 如圖所示。圖3-8仿真模型圖3-9濾波效果結(jié)果顯示，仿真模型中顯示40HZ, 75H Z, 100HZ三個不同頻率的正弦信相加 后的波形，干擾后的混合信號是一個夾雜多頻率的信號，我們所需要的有用信 號(40HZ)在傳輸過程中不可避免的受到其他頻率信號干擾，我們的目的就

54、是 得到所需要的40HZ頻段信號，盡可能濾除干擾信號，濾波效果的圖中就是濾波 后信號與原始抽樣信號的的波形圖，經(jīng)濾波后，信號周期約為0.025s即頻率為40HZ,說明濾波效果很好。通過分析，說明SIMULINK中各種非常有用的工具 箱不僅對于設(shè)計IIR數(shù)字濾波器非常有用，而且對于整個信號仿真處理具有相當(dāng) 可視化的效果，從仿真的效果來看，是達(dá)到了技術(shù)指標(biāo)的要求?？偨Y(jié)通過這幾個實(shí)驗(yàn)，對設(shè)計數(shù)字濾波器的整個過程有了很好的掌握。其中對 雙線性變換法，巴特沃斯設(shè)計模擬濾波器的運(yùn)用，也更加清楚了。通過對數(shù)字帶通濾波器的設(shè)計，熟悉了MATLAB的運(yùn)行環(huán)境，初步掌握了MATLAB語言在數(shù)字信號處理中一些基本庫

55、函數(shù)的調(diào)用和編寫基本程序等應(yīng)用； 熟悉了濾波器設(shè)計的一般原理，對濾波器有了一個感性的認(rèn)識；學(xué)會了數(shù)字高 通濾波器設(shè)計的一般步驟；加深了對濾波器設(shè)計中產(chǎn)生誤差的原因以及雙線性 變換法優(yōu)缺點(diǎn)的理解和認(rèn)識?？傊?，使理論聯(lián)系了實(shí)際，鞏固并深化了對課本 基本知識的認(rèn)識和理解，使理論得以升華。致謝本論文是在雷老師的親切關(guān)懷和悉心指導(dǎo)下完成的。他嚴(yán)肅的科學(xué)態(tài)度， 嚴(yán)謹(jǐn)?shù)闹螌W(xué)精神，精益求精的工作作風(fēng)，深深地感染和激勵著我。從課題的選 擇到工程的最終完成，雷老師都始終給予我細(xì)心的指導(dǎo)和不懈的支持。一直以 來，雷老師不僅在學(xué)業(yè)上給我以精心指導(dǎo)，同時還在思想、生活上給我以無微 不至的關(guān)懷，在此謹(jǐn)向雷老師致以誠摯的謝

56、意和崇高的敬意。在此，我還要感謝在一起愉快的度過畢業(yè)設(shè)計生活的同學(xué)，正是由于你們 的幫助和支持，我才能克服一個一個的困難和疑惑，直至本文的順利完成。特 別感謝我的同學(xué)，她們對本課題做了不少工作，給予我不少的幫助。在論文即將完成之際，我的心情無法平靜，從開始進(jìn)入課題到論文的順利 完成，有多少可敬的師長、同學(xué)、朋友給了我無言的幫助，在這里請接受我誠 摯的謝意!最后我還要感謝培，謝謝你們!參考文獻(xiàn)1倪養(yǎng)華數(shù)字信號處理與實(shí)現(xiàn)上海：上海交通大學(xué)出版社，19982飛思科技產(chǎn)品研發(fā)中心 MATLAB7輔助信號處理技術(shù)與應(yīng)用 北京：電子 工業(yè)出版社，20053樓順天，李博苗 基于MATLAB的系統(tǒng)分析與設(shè)計一信號處理 西安：西 安電子科技大學(xué)出版社，1998.4張葛祥，李娜 MATLAB仿真技術(shù)與應(yīng)用 北京：清華大學(xué)出版社，20035胡廣書數(shù)字信號處理、理論、算法與實(shí)現(xiàn)M 北京：清華大學(xué)出版社，1997陳希林，肖明清 一種LabWindows/CVI與MATLAB混合編程的實(shí)現(xiàn)方法J.微計算機(jī)信息，20057劉波.MATLAB信號處理.北京：電子工業(yè)出版社，200