DSP的數(shù)字濾波設(shè)計(jì).docVIP

DSP的數(shù)字濾波器設(shè)計(jì)作者姓名：王珊珊 班級(jí)：200404070526 指導(dǎo)教師：王曉亮摘 要隨著集成電路技術(shù)的發(fā)展，各種新型的大規(guī)模和超大規(guī)模集成電路不斷涌現(xiàn)，集成電路技術(shù)與計(jì)算機(jī)技術(shù)結(jié)合在一起，使得對(duì)數(shù)字信號(hào)處理系統(tǒng)功能的要求越來越強(qiáng)。在數(shù)字信號(hào)處理中數(shù)字濾波占有極其重要的地位。數(shù)字濾波是語音信號(hào)、圖象處理模式識(shí)別和譜分析等領(lǐng)域中的一個(gè)基本的處理技術(shù)。本設(shè)計(jì)基于DSP，完成FIR與IIR數(shù)字濾波器的設(shè)計(jì)與軟硬件仿真。數(shù)字濾波器設(shè)計(jì)的常用方法包括窗函數(shù)法、頻率采樣法和最優(yōu)化法等。本次設(shè)計(jì)在了解了數(shù)字濾波器設(shè)計(jì)原理的基礎(chǔ)上，分別采用了窗函數(shù)法對(duì)不同指標(biāo)的低通、高通、帶通、帶阻以及切比雪夫一型濾波器進(jìn)行設(shè)計(jì)，在matlab軟件下繪畫出各種設(shè)計(jì)方法的幅頻響應(yīng)和相頻響應(yīng)曲線。同時(shí)在CCS環(huán)境以及DSP硬件平臺(tái)上仿真，對(duì)比仿真后波形圖，看得出實(shí)際圖形是否與預(yù)期圖形相符合。關(guān)鍵詞：數(shù)字濾波器；Matlab軟件；窗函數(shù)法；DSPDigital Filter Design Based On DSPAbstract： With the development of the integrated circuit technology, each kind of new large-scale and the ultra large scale integrated circuit emerges increasingly, the integrated circuit technology and the computer technology unified together, it makes digital signal processing system function request to be getting stronger and stronger. The digital filtering holds the extremely important status in the digital signal processing. The digital filtering is a basic processing technology in the area of voice signal, imagery processing domain, and pattern recognition and spectrum analysis.Based on DSP, this article describes the software and hardware simulation of FIR and IIR digital filter. The common method once used on designing digital filter includes window function method, frequency sampling law, most superior reduction and so on. Upon understanding the fundamental principle of designing digital filter, we will use the window function method separately with the different target on low pass, high pass, the band pass, the band stop filter as well as Chebyshev filters to make it work. Drawing the frequency response and the frequency response curve of each method under the matlab software. Meanwhile, with the DSP hardware platform, contrasting the oscillogram after simulation, it will draw the conclusion whether the actual graph matches with the anticipated graph under the CCS environment.Keywords: digital filter; matlab software; window functions; DSP目 錄第1章 前 言11.1 研究背景11.2 研究意義21.3 主要內(nèi)容與要求3第2章 系統(tǒng)開發(fā)平臺(tái)42.1 CCS開發(fā)環(huán)境42.2 硬件平臺(tái)42.3 計(jì)算機(jī)輔助設(shè)計(jì)方法6第3章 數(shù)字濾波器的基本理論73.1 數(shù)字濾波器的基本結(jié)構(gòu)73.1.1 FIR濾波器的基本結(jié)構(gòu)73.1.2 IIR濾波器的基本結(jié)構(gòu)83.2 數(shù)字濾波器的設(shè)計(jì)原理93.3 FIR數(shù)字濾波器的設(shè)計(jì)方法103.3.1 類型，系數(shù)對(duì)稱，N奇數(shù)103.3.2 類型，系數(shù)對(duì)稱，N偶數(shù)113.3.3 類型，系數(shù)反對(duì)稱，N奇數(shù)123.3.4 類型，系數(shù)反對(duì)稱，N偶數(shù)123.3.5 FIR數(shù)字濾波器的窗函數(shù)設(shè)計(jì)方法153.3.6 利用頻率取樣法設(shè)計(jì)FIR濾波器163.3.7 其他幾種窗函數(shù)的介紹193.4 IIR數(shù)字濾波器的設(shè)計(jì)方法283.5 FIR濾波器與IIR濾波器的比較30第4章 利用MATLAB設(shè)計(jì)數(shù)字濾波器324.1 FIR數(shù)字濾波器324.1.1 FIR低通數(shù)字濾波器324.1.2 FIR高通數(shù)字濾波器324.1.3 FIR帶通濾波器334.1.4 FIR帶阻數(shù)字濾波器354.2 IIR數(shù)字濾波器354.2.1 IIR低通數(shù)字濾波器354.2.2 IIR高通數(shù)字濾波器364.2.3 切比雪夫函數(shù)設(shè)計(jì)IIR高通濾波器37第5章 數(shù)字濾波器的DSP實(shí)現(xiàn)445.1 ICETEKF2812-A評(píng)估板技術(shù)指標(biāo)445.2 ICETEKF2812-A評(píng)估板實(shí)物圖445.3 ICETEK-F2812-A器件分布圖455.4 ICETEKF2812A評(píng)估板原理框圖455.5 現(xiàn)場(chǎng)硬件連接圖465.6 DSP實(shí)現(xiàn)465.6.1 有限沖擊響應(yīng)濾波器（FIR）算法465.6.2 無限沖激響應(yīng)濾波器（IIR）算法495.6.3 快速傅立葉變換（FFT）算法525.6.4 語音信號(hào)的FIR濾波55結(jié) 論58致 謝59參考文獻(xiàn)60IV本科畢業(yè)設(shè)計(jì)（論文）第1章 前 言1.1 研究背景隨著集成電路技術(shù)的發(fā)展，各種新型的大規(guī)模和超大規(guī)模集成電路不斷涌現(xiàn)，集成電路技術(shù)與計(jì)算機(jī)技術(shù)結(jié)合在一起， 使得對(duì)數(shù)字信號(hào)處理系統(tǒng)功能的要求越來越強(qiáng)。DSP技術(shù)已在通信、控制、信號(hào)處理、儀器儀表、醫(yī)療、家電等很多領(lǐng)域得到了越來越廣泛的應(yīng)用。在數(shù)字信號(hào)處理中，數(shù)字濾波占有極其重要的地位。數(shù)字濾波在語音信號(hào)、圖象處理、模式識(shí)別和譜分析等領(lǐng)域中是一個(gè)基本的處理技術(shù)。數(shù)字濾波與模擬濾波相比，數(shù)字濾波具有很多突出的優(yōu)點(diǎn)，主要是因?yàn)閿?shù)字濾波器是過濾時(shí)間離散信號(hào)的數(shù)字系統(tǒng)，它可以用軟件（計(jì)算機(jī)程序）或用硬件來實(shí)現(xiàn)，而且在兩種情況下都可以用來過濾實(shí)時(shí)信號(hào)或非實(shí)時(shí)信號(hào)。盡管數(shù)字濾波器這個(gè)名稱一直到六十年代中期才出現(xiàn)，但是隨著科學(xué)技術(shù)的發(fā)展及計(jì)算機(jī)的更新普及，數(shù)字濾波器有著很好的發(fā)展前景。目前，國外有許多院校和科研機(jī)構(gòu)在研究基于FPGA的DSP應(yīng)用，比較突出的有Denmark大學(xué)的研究小組正在從事FPGA實(shí)現(xiàn)數(shù)字濾波器的研究。加州大學(xué)洛杉磯分校的研究小組采用運(yùn)行時(shí)重構(gòu)技術(shù)開發(fā)了一種視頻通訊系統(tǒng)，該系統(tǒng)用一片F(xiàn)PGA可每幀重構(gòu)四次完成視頻圖像壓縮和傳送的操作。此外，他們還在進(jìn)行Mojave項(xiàng)目的開發(fā)工作，力圖采用運(yùn)行時(shí)重構(gòu)技術(shù)來實(shí)現(xiàn)自動(dòng)目標(biāo)識(shí)別應(yīng)用。我國在DSP技術(shù)起步較早，產(chǎn)品的研究開發(fā)成績斐然，基本上與國外同步發(fā)展，而在FPGA方面起步較晚。全國有100來所高等院校從事DSP&FPGA的教學(xué)和科研，除了一部分DSP芯片需要從國外進(jìn)口外，在信號(hào)處理理論和算法方面，與國外處于同等水平.而在FPGA信號(hào)處理和系統(tǒng)方面，有了喜人的進(jìn)展，正在進(jìn)行與世界先進(jìn)國家同樣的研究.西北工業(yè)大學(xué)和國防科學(xué)技術(shù)大學(xué)的ATR實(shí)驗(yàn)室采用了FPGA可重構(gòu)計(jì)算系統(tǒng)進(jìn)行機(jī)載圖像處理和自動(dòng)目標(biāo)識(shí)別，主要是利用該系統(tǒng)進(jìn)行復(fù)雜的卷積運(yùn)算，同時(shí)利用它的可變?nèi)嵝詠磉_(dá)到自適應(yīng)的目的。20世紀(jì)60年代起，由于計(jì)算機(jī)技術(shù)、集成工藝和材料工業(yè)的發(fā)展，濾波器的發(fā)展上了一個(gè)新臺(tái)階，朝著低功耗、高精度、小體積、多功能、穩(wěn)定可靠和價(jià)廉等方向努力，其中高精度、小體積、多功能、穩(wěn)定可靠成為70年代以后的主攻方向，導(dǎo)致數(shù)字濾波器、RC有源濾波器、開關(guān)電容濾波器和電荷轉(zhuǎn)移器等各種濾波器的飛速發(fā)展。到70年代后期，上述幾種濾波器的單片集成己被研制出來并得到應(yīng)用，90年代至現(xiàn)在主要致力于把各類濾波器應(yīng)用于各類產(chǎn)品的開發(fā)和研制。當(dāng)然，對(duì)濾波器本身的研究仍在不斷進(jìn)行。1.2 研究意義幾乎在所有的工程技術(shù)領(lǐng)域中都會(huì)涉及到信號(hào)的處理問題，其信號(hào)表現(xiàn)形式有電、磁、機(jī)械以及熱、光、聲等。信號(hào)處理的目的一般是對(duì)信號(hào)進(jìn)行分析、變換、綜合、估值與識(shí)別等。如何在較強(qiáng)的噪聲背景下提取出真正的信號(hào)或信號(hào)的特征，并將其應(yīng)用于工程實(shí)際是信號(hào)處理的首要任務(wù)。根據(jù)處理對(duì)象的不同，信號(hào)處理技術(shù)分為模擬信號(hào)處理系統(tǒng)和數(shù)字信號(hào)處理系統(tǒng)。數(shù)字信號(hào)處理(Digital signal Processing，DSP)與模擬信號(hào)處理相比有許多優(yōu)點(diǎn)，如相對(duì)于溫度和工藝的變化數(shù)字信號(hào)要比模擬信號(hào)更穩(wěn)健，在數(shù)字表示中，精度可以通過改變信號(hào)的字長來更好地控制，所以DSP技術(shù)可以在放大信號(hào)的同時(shí)去除噪聲和干擾，而在模擬信號(hào)中信號(hào)和噪聲同時(shí)被放大，數(shù)字信號(hào)還可以不帶誤差地被存儲(chǔ)和恢復(fù)、發(fā)送和接收、處理和操縱。許多復(fù)雜的系統(tǒng)可以用高精度、大信噪比和可重構(gòu)的數(shù)字技術(shù)來實(shí)現(xiàn)。數(shù)字信號(hào)處理中一個(gè)非常重要且應(yīng)用普遍的技術(shù)就是數(shù)字濾波。所謂數(shù)字濾波，是指其輸入、輸出均為數(shù)字信號(hào)，通過一定的運(yùn)算關(guān)系改變輸入信號(hào)所含的頻率成分的相對(duì)比例或?yàn)V除某些頻率成分，達(dá)到提取和加強(qiáng)信號(hào)中的有用成份，消弱無用的干擾成份的目的。數(shù)字濾波與模擬濾波相比，有精度高、可靠性高、靈活性好等突出優(yōu)點(diǎn)，可以滿足對(duì)幅度和相位的嚴(yán)格要求，還能降低開發(fā)費(fèi)用，縮短研制到應(yīng)用的時(shí)間，在很多領(lǐng)域逐步代替了傳統(tǒng)的模擬信號(hào)系統(tǒng)。數(shù)字濾波器，在數(shù)字信號(hào)處理中有著廣泛的應(yīng)用，因此，無論是在理論研究上還是在如通訊、HDTV(高清晰度電視)、雷達(dá)、圖象處理、數(shù)字音頻等實(shí)際應(yīng)用上都有著美好的技術(shù)前景和巨大的實(shí)用價(jià)值1 張亞妮.基于MATLAB的數(shù)字濾波器設(shè)計(jì)J. 遼寧工程技術(shù)大學(xué)學(xué)報(bào)，2005.?；趍atlab實(shí)現(xiàn)數(shù)字濾波設(shè)計(jì)，可通過修改濾波器的參數(shù)十分方便地改變?yōu)V波器的特性，因此我們有必要對(duì)濾波器的設(shè)計(jì)方法進(jìn)行研究，理解其工作原理優(yōu)化設(shè)計(jì)方法，設(shè)計(jì)開發(fā)穩(wěn)定性好的濾波器系統(tǒng)。掌握濾波器的設(shè)計(jì)技術(shù)和原理能為在通信領(lǐng)域、信號(hào)處理領(lǐng)域等諸多領(lǐng)域中對(duì)數(shù)字濾波器的設(shè)計(jì)提供技術(shù)和準(zhǔn)備，這不僅具有重要的理論意義同時(shí)還具有重要的現(xiàn)實(shí)意義。1.3 主要內(nèi)容與要求本次畢業(yè)設(shè)計(jì)中主要完成的內(nèi)容包括（1）通過書籍閱讀和文獻(xiàn)搜索學(xué)習(xí)MATLAB和數(shù)字濾波器，并掌握MATLAB的使用和數(shù)字濾波器的設(shè)計(jì)方法。（2）數(shù)字濾波器原理分析與設(shè)計(jì)方法概述。分析數(shù)字濾波器的原理，對(duì)FIR以及IIR數(shù)字濾波器常見設(shè)計(jì)方法的介紹。（3）FIR以及IIR數(shù)字濾波器的matlab軟件實(shí)現(xiàn)（4）設(shè)計(jì)若干類型的FIR以及IIR（低通、高通、帶通、帶阻以及切比雪夫）數(shù)字濾波器，并在matlab上實(shí)現(xiàn)。（5）數(shù)字濾波器的硬件電路實(shí)現(xiàn)選擇LF2812芯片設(shè)計(jì)外圍電路實(shí)現(xiàn)數(shù)字濾波器的硬件設(shè)計(jì)（6）運(yùn)用LF2812實(shí)現(xiàn)IIR以及FIR算法實(shí)驗(yàn)，并以此為基礎(chǔ)，設(shè)計(jì)并實(shí)現(xiàn)FIR語音信號(hào)濾波。（7）各項(xiàng)試驗(yàn)完成后，對(duì)各個(gè)實(shí)驗(yàn)結(jié)果進(jìn)行分析，比對(duì)?？磳?shí)際波形與理想波形的不同之處，比較濾波前后的時(shí)域頻域波形。第2章 系統(tǒng)開發(fā)平臺(tái)2.1 CCS開發(fā)環(huán)境CCS提供了配置、建立、調(diào)試、跟蹤和分析程序的工具，它便于實(shí)時(shí)、嵌入式信號(hào)處理程序的編制和測(cè)試，它能夠加速開發(fā)進(jìn)程，提高工作效率。CCS提供了基本的代碼生成工具，它們具有一系列的調(diào)試、分析能力。CCS支持所示的開發(fā)周期的所有階段。圖 2-1 開發(fā)周期的所有階段2.2 硬件平臺(tái)下面來看一下即將用到的ICETEK DSP 實(shí)驗(yàn)箱，如圖2-2 所示。圖 2-2 實(shí)驗(yàn)箱的組成ICETEK DSP 教學(xué)實(shí)驗(yàn)箱主要由以下幾個(gè)部分組成(1)箱蓋保護(hù)實(shí)驗(yàn)箱設(shè)備；保存教材、使用手冊(cè)、實(shí)驗(yàn)指導(dǎo)書、各種實(shí)驗(yàn)用的連線；可拆卸，在實(shí)驗(yàn)中可從箱體上拆下。(2)箱體裝載實(shí)驗(yàn)箱設(shè)備；左側(cè)外壁上有一個(gè)標(biāo)準(zhǔn)外接電源線插孔；通過固定螺絲與實(shí)驗(yàn)箱底板連為一體。(3)底板固定各模塊；提供電源開關(guān)、實(shí)驗(yàn)用直流電源插座、A/D D/A 輸入輸出插座、各模塊直流供電插座、信號(hào)插座、信號(hào)源輸出插座、測(cè)試點(diǎn)；實(shí)現(xiàn)顯示控制模塊和 DSP 評(píng)估板模塊的信號(hào)互連。(4)信號(hào)源兩組、三路輸出，采用+5V 電源供電；提供切換選擇輸出方波、三角波、正弦波和白噪聲，另可選擇輸出頻率范圍(10Hz-100Hz，100Hz-1KHz，1KHz-10KHz，10KHz-30KHz)，還可進(jìn)行頻率和幅度(0-3.3V)的微調(diào)。(5)仿真器模塊固定 ICETEK 仿真器，支持 PP 型和 USB 型；提供 PP 型仿真器供電+5V 電源插座；仿真器可從底板上拆下單獨(dú)使用或更換。(6)顯示控制模塊通過信號(hào)線連接到底板；從底板提供的+5V 和+12V 直流電源插座輸入電源；提供液晶圖形顯示(128x64 象素)，發(fā)光二極管陣列顯示（8x8 點(diǎn)）（本功能為 V4.3 版本特有)，指示燈(12 只，分為紅、黃、綠三種顏色)，四相步進(jìn)電機(jī)，直流電機(jī)，鍵盤，蜂鳴器。顯示控制模塊可從底板上拆下更換。(7)DSP 評(píng)估板模塊固定各種 DSP 評(píng)估板；提供+5V 直流電源插座(兩個(gè)位置)；34Pin 信號(hào)線插座(4 個(gè))，用于連接 DSP 評(píng)估板和實(shí)驗(yàn)箱底板。DSP 評(píng)估板模塊可從底板上拆下更換。 圖 2-3DSP實(shí)驗(yàn)箱結(jié)構(gòu)圖本文第五章DSP的硬件實(shí)現(xiàn)將要用到上述的CCS開發(fā)環(huán)境，需要借助上述試驗(yàn)箱完成有關(guān)濾波器的相關(guān)軟硬件仿真。后文中將詳細(xì)介紹，在此不再贅述。2.3 計(jì)算機(jī)輔助設(shè)計(jì)方法 Ingle V K,Proakis J G. Digital Processing Using MATLABM. PWS Publishing Company, 2003計(jì)算機(jī)輔助設(shè)計(jì)方法是集電路理論、網(wǎng)絡(luò)圖論、數(shù)值分析、矩陣運(yùn)算、原件建模、優(yōu)化設(shè)計(jì)、高級(jí)計(jì)算機(jī)語言等多交叉學(xué)科于一身的新領(lǐng)域，它把計(jì)算機(jī)的快速、高精度、大存儲(chǔ)量、嚴(yán)格的邏輯判斷和優(yōu)良的數(shù)據(jù)處理能力與人的思維創(chuàng)造能力充分結(jié)合起來，極大的簡化了數(shù)字濾波器的設(shè)計(jì)過程。在優(yōu)秀科技應(yīng)用軟件MATLAB的信號(hào)處理工具箱中，提供了一整套模擬、數(shù)字濾波器的設(shè)計(jì)命令和運(yùn)算函數(shù)，方便準(zhǔn)確，簡單易行，使得設(shè)計(jì)人員除了可以上述傳統(tǒng)設(shè)計(jì)步驟快速地進(jìn)行較復(fù)雜高階選頻濾波器的計(jì)算、分析外，還可以通過原型變換法直接進(jìn)行各種典型數(shù)字濾波器設(shè)計(jì)，即應(yīng)用MATLAB設(shè)計(jì)工具從模擬原型直接變換成滿足原定頻域指標(biāo)要求的數(shù)字濾波器。傳統(tǒng)設(shè)計(jì)方法思路清晰，步驟詳盡，可參閱公式、手冊(cè)循章而行。但是由于計(jì)算繁瑣，手工計(jì)算大多只能用來進(jìn)行簡單低階選頻濾波器的設(shè)計(jì)。而計(jì)算機(jī)輔助設(shè)計(jì)在擁有了傳統(tǒng)的設(shè)計(jì)的優(yōu)勢(shì)上更結(jié)合了一些特定屬于自己的優(yōu)勢(shì)，能進(jìn)行各種相關(guān)的比較復(fù)雜的濾波器設(shè)計(jì)，彌補(bǔ)了傳統(tǒng)上的不足。所以，我最終選擇了計(jì)算機(jī)輔助設(shè)計(jì)的方法，在MATLAB環(huán)境下進(jìn)行數(shù)字濾波器的設(shè)計(jì)和仿真。第3章 數(shù)字濾波器的基本理論3.1 數(shù)字濾波器的基本結(jié)構(gòu)作為線性時(shí)不變系統(tǒng)的數(shù)字濾波器可以用系統(tǒng)函數(shù)來表示，而實(shí)現(xiàn)一個(gè)系統(tǒng)函數(shù)表達(dá)式所表示實(shí)物系統(tǒng)可以用兩種方法一種方法是采用計(jì)算機(jī)軟件來實(shí)現(xiàn)；另一種方法是用加法器、乘法器和延時(shí)器等元件設(shè)計(jì)出專用的數(shù)字硬件系統(tǒng)，即硬件實(shí)現(xiàn)。不論軟件實(shí)現(xiàn)還是硬件實(shí)現(xiàn)，在濾波器設(shè)計(jì)過程中，由同一系統(tǒng)函數(shù)可以構(gòu)成很多不同的運(yùn)算結(jié)構(gòu)。對(duì)于無限精度的系數(shù)和變量，不同結(jié)構(gòu)可能是等效的，與其輸入和輸出特性無關(guān)；但是在系數(shù)和變量精度有限的情況下，不同運(yùn)算結(jié)構(gòu)的性能就有很大的差異。因此，有必要對(duì)離散時(shí)間系統(tǒng)的結(jié)構(gòu)有一些基本知識(shí)。一個(gè)數(shù)字濾波器可以用系統(tǒng)函數(shù)表示為 （3-1）由這樣的系統(tǒng)函數(shù)可以得到表示系統(tǒng)輸入與輸出關(guān)系的常系數(shù)線性差分方程為 （3-2）可見數(shù)字濾波器的功能就是把輸入序列x（n）通過一定的運(yùn)算變換成輸出序列 y（n）。不同的運(yùn)算處理方法決定了濾波器實(shí)現(xiàn)結(jié)構(gòu)的不同。無限沖激響應(yīng)濾波器的單位抽樣響應(yīng)h（n）是無限長的，其差分方程如（3-2）所示，是遞歸式的，即結(jié)構(gòu)上存在著輸出信號(hào)到輸入信號(hào)的反饋，其系統(tǒng)函數(shù)具有（3-1）式的形式，因此在z平面的有限區(qū)間有極點(diǎn)存在。3.1.1 FIR濾波器的基本結(jié)構(gòu)FIR濾波器的單位抽樣響應(yīng)為有線長度，一般采用非遞歸形式實(shí)現(xiàn)。通常的FIR數(shù)字濾波器有橫截性和級(jí)聯(lián)兩種。FIR濾波器實(shí)現(xiàn)的基本結(jié)構(gòu)有(1)FIR濾波器的橫截結(jié)構(gòu)表示系統(tǒng)輸入輸出關(guān)系的差分方程可以寫作 （3-3）若具有對(duì)稱特性，即此FIR濾波器具有線性相位。(2)FIR濾波器的級(jí)聯(lián)型結(jié)構(gòu)將分解成系數(shù)二階因子的乘積形式 （3-4）這時(shí)FIR濾波器可用二階節(jié) 的級(jí)聯(lián)結(jié)構(gòu)來實(shí)現(xiàn)，每個(gè)二階節(jié)用橫截性結(jié)構(gòu)實(shí)現(xiàn)。這種結(jié)構(gòu)的每一節(jié)控制一對(duì)零點(diǎn)，因而在需要控制傳輸零點(diǎn)時(shí)可以采用這種結(jié)構(gòu)。3.1.2 IIR濾波器的基本結(jié)構(gòu)（1）IIR濾波器的直接型結(jié)構(gòu)在實(shí)際中很少采用上述結(jié)構(gòu)實(shí)現(xiàn)高階系統(tǒng)，而是把高階變成一系列不同組合的低階系統(tǒng)來實(shí)現(xiàn)。（2）IIR濾波器的級(jí)聯(lián)型結(jié)構(gòu)特點(diǎn)系統(tǒng)實(shí)現(xiàn)簡單，只需要一個(gè)二階系統(tǒng)通過改變輸入系數(shù)既可以完成;極點(diǎn)位置可單獨(dú)調(diào)整運(yùn)算速度快各二階網(wǎng)絡(luò)的誤差互不影響，總的誤差小，對(duì)字長要求低。缺點(diǎn)不能直接調(diào)整零極點(diǎn)，因多個(gè)二階節(jié)的零點(diǎn)不是整個(gè)系統(tǒng)函數(shù)的零點(diǎn)，當(dāng)需要準(zhǔn)確的傳輸零點(diǎn)時(shí)，級(jí)聯(lián)型是最適合的。（3）IIR濾波器的并聯(lián)型結(jié)構(gòu)優(yōu)點(diǎn)簡化實(shí)現(xiàn)，用一個(gè)二階節(jié)，通過變換系數(shù)就可以實(shí)現(xiàn)整個(gè)系統(tǒng)；極、零點(diǎn)可單獨(dú)使用控制、調(diào)整，調(diào)整、只單獨(dú)調(diào)整了第i對(duì)零點(diǎn)，調(diào)整、則單獨(dú)調(diào)整了第i對(duì)極點(diǎn)；可流水線操作3陳懷琛.數(shù)字信號(hào)處理教程Matlab釋義與實(shí)現(xiàn)M.北京電子工業(yè)出版社，20083.2 數(shù)字濾波器的設(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è)組成部分。理想濾波器是非因果的，物理上不可能實(shí)現(xiàn)。數(shù)字濾波器的設(shè)計(jì)問題就是尋找無限脈沖響應(yīng)（IIR）數(shù)字濾波器的差分方程和對(duì)應(yīng)的系統(tǒng)函數(shù)（兩個(gè)方程）中的一組系和，是的其性能在某種意義上逼近所要求的特性。IIR數(shù)字濾波器的設(shè)計(jì)方法一般分為下述兩類先設(shè)計(jì)一個(gè)合適的模擬濾波器，然后變換成給定指標(biāo)的數(shù)字濾波器。這種方法很方便，這是因?yàn)槟M濾波器的設(shè)計(jì)方法已經(jīng)很成熟，他不僅有完整的設(shè)計(jì)公式，還有完整的圖表供查閱。這種方法的設(shè)計(jì)步驟如下將給定的數(shù)字濾波器的技術(shù)指標(biāo)轉(zhuǎn)換為模擬濾波器的技術(shù)指標(biāo)根據(jù)轉(zhuǎn)換后的技術(shù)指標(biāo)設(shè)計(jì)模擬原型濾波器再按照一定規(guī)則將模擬濾波器轉(zhuǎn)換為數(shù)字濾波器直接在頻域或者時(shí)域中進(jìn)行設(shè)計(jì)，這是一種最優(yōu)化設(shè)計(jì)方法。由于要接聯(lián)立方程，因此需要計(jì)算機(jī)輔助進(jìn)行設(shè)計(jì)。FIR濾波器的設(shè)計(jì)方法與IIR濾波器有很大不同，F(xiàn)IR濾波器的設(shè)計(jì)是通過對(duì)理想濾波器的頻率特性作某種逼近得到的。常用的方法有窗函數(shù)法，頻率取樣法和計(jì)算機(jī)輔助最優(yōu)化設(shè)計(jì)法。FIR濾波器的最大優(yōu)點(diǎn)就是具有線性相位，而設(shè)計(jì)的IIR濾波器的相位都是非線性的。數(shù)字濾波器的設(shè)計(jì)要經(jīng)過三個(gè)步驟確定指標(biāo)在設(shè)計(jì)一個(gè)濾波器前，必須有一些指標(biāo)。這些指標(biāo)要根據(jù)應(yīng)用確定。模型逼近一旦確定了技術(shù)指標(biāo)就可以利用前面的基本原理和關(guān)系式，提出一個(gè)濾波器模型來逼近給定的指標(biāo)體系。這是濾波器設(shè)計(jì)所要研究所的主要問題。實(shí)現(xiàn)上面的兩布的結(jié)果得到的濾波器，通常是以差分方程、系統(tǒng)函數(shù)或者脈沖響應(yīng)來描述的。根據(jù)這個(gè)描述，可以用硬件或計(jì)算機(jī)軟件來實(shí)現(xiàn)它了。3.3 FIR數(shù)字濾波器的設(shè)計(jì)方法線性相位條件要求濾波器分子系數(shù)滿足對(duì)稱性條件。設(shè)濾波器的系數(shù)長度為N，則這些系數(shù)應(yīng)關(guān)于中心點(diǎn)t=n 0=(N-1)/2對(duì)稱。偶對(duì)稱時(shí)，h(n)=h(N-n-1);而奇對(duì)稱時(shí)，h(n)=-h(N-n-1).再考慮到N可以為奇數(shù)或偶數(shù)，總共有四種類型的線性相位FIR濾波器。討論線性相位濾波器頻率響應(yīng)是需要引進(jìn)幅特性正負(fù)號(hào)的概念，以往常設(shè)幅特性為正數(shù)，因?yàn)榉匦缘姆刺?hào)可以用響特性加減來補(bǔ)償。當(dāng)相位特性要求線性，不得隨便增減時(shí)，幅特性就必須分出正負(fù)，稱為符幅特性。為了就算符幅特性，需要直接作解析計(jì)算，求出頻率特性的解析式，并將它分解為幅度和相位特性，即,如果（）滿足線性相位條件，這時(shí)的A()就是符幅特性。對(duì)于長度為N的序列h(n),先考慮兩個(gè)對(duì)稱的系數(shù)h(n)和h（N-1-n）,它們的傅里葉變換為4 高西全，丁玉美.數(shù)字信號(hào)處理. M.陜西西安電子科技大學(xué)出版社，2008. (3-5)3.3.1 類型，系數(shù)對(duì)稱，N奇數(shù)序列h的下標(biāo)n從0算起，到N-1為止，中點(diǎn)位置是L=(N-1)/2.h(L)是一個(gè)孤項(xiàng)，其他可配成L對(duì)系數(shù)。將他們的符幅特性加起來，得到 (3-6) 為了使四種類型濾波器的公式簡明統(tǒng)一，引進(jìn)兩個(gè)參數(shù)。一個(gè)是=（N-1）/2，它就是T值，用以反映延遲的樣本或群延遲，它可以是分?jǐn)?shù)。另一個(gè)是L=(N-1)/2把T值向下取整，因?yàn)橐袻用作下標(biāo)，它必須是整數(shù)。當(dāng)N為奇數(shù)時(shí)，T=L。在推導(dǎo)中，還注意前一個(gè)求和號(hào)是從0到L-1，后一個(gè)求和號(hào)是從0到L，那是把孤項(xiàng)也當(dāng)成成對(duì)相求和，然后再減去一個(gè)孤項(xiàng)h(L)。這樣做的目的也是使四種類型濾波器的公式統(tǒng)一。顯然，為項(xiàng)角項(xiàng)，相角-，它與成嚴(yán)格的線性關(guān)系；A（）就是符幅相。 （3-7)可見，其符幅特性由L+1個(gè)余弦項(xiàng)疊加組成。n=0時(shí)的符幅相為，在頻段間波動(dòng)個(gè)周期余弦函數(shù)在=0，pi和2pi處都不等于零,因此類型線性相位濾波器既可以用作低通濾波器;也可以用作高通濾波器;且可用作帶通和帶阻濾波器。3.3.2 類型，系數(shù)對(duì)稱，N偶數(shù)如果N為偶數(shù)，那么全部系數(shù)都可以配對(duì)，不會(huì)出現(xiàn)中心點(diǎn)的單項(xiàng)，一共有N/2組對(duì)稱出現(xiàn)的系數(shù)，頻率特性就成為 （3-8）此時(shí)T=（N-1）/2將不是整數(shù)，也就是說，對(duì)稱中心將在兩個(gè)樣本點(diǎn)的中間。其相角特性仍為（）=T，它仍與成嚴(yán)格的線性關(guān)系。用整數(shù)L=(N-1)/2=N/2-1來表示求和號(hào)的上限。得出其符幅特性的表達(dá)式 （3-9）可見，它的符幅特性也由L+1個(gè)余弦分量構(gòu)成。這樣，不管N為奇數(shù)或偶數(shù)時(shí)，L都是整數(shù)，求和的項(xiàng)數(shù)都是L+1。由于T=(N-1)/2是分?jǐn)?shù)，所以有余弦分量中都含有一個(gè)0.5。在pi處，就會(huì)出現(xiàn)符幅特性必定為零的情況。所以類型不能用作高通和帶阻濾波器.3.3.3 類型，系數(shù)反對(duì)稱，N奇數(shù)當(dāng)系數(shù)序列反對(duì)稱時(shí)，只要把類型結(jié)果中的cos 換成sin就得到 （3-10）當(dāng)N為奇數(shù)時(shí)，T=L=整數(shù)。相特性為（）=/2-（N-1）=/2-T雖然仍與成線性關(guān)系，但多了一個(gè)常數(shù)項(xiàng)/2。注意反對(duì)稱序列的對(duì)稱中心n=L處，序列的值h（L）=0.所以可以放入求和號(hào)中，得到它的符幅特性為 （3-11）正弦函數(shù)在=0，和2處都等于零，因此類型線性相位濾波器的符幅特性在=0，和處都等于零。它不能用作低通，也不能用作高通濾波器。所以不適宜用于選頻濾波器。3.3.4 類型，系數(shù)反對(duì)稱，N偶數(shù)與類型濾波器相仿，若N為偶數(shù)，則共有N/2組數(shù)值相同反對(duì)稱出現(xiàn)的系數(shù)，不出現(xiàn)單項(xiàng)，頻率特性就成為 （3-12）此處仍有T=（N-1）/2，因此T將不是整數(shù)，帶有小數(shù)部分0.5.故對(duì)稱中心將在兩個(gè)樣本點(diǎn)中間。它的相角特性為()= /2-T,雖然與呈線性關(guān)系，但是多了常數(shù)項(xiàng)/2，故不通過原點(diǎn)。類型的濾波器符幅特性為 （3-13）因?yàn)長是整數(shù)，它的符幅特性也由L+1=N/2個(gè)正弦分量構(gòu)成。因?yàn)閟in(0)=0,在=0處，符幅特性必定為A（0）=0.這樣的特性顯然不能用作低通和帶阻濾波器，實(shí)際上也不適用于選頻濾波器。把這四種類型的線性相位濾波器的相角特性和符幅特性用對(duì)比的方法畫在一張圖上5表3-1 四種類型的線性相位濾波器的相位與符幅特性 Ingle V K，Proakis J G. Digital Signal Processing Using MATL AB M，New YorkPWS Publishing Company,1997 符幅特性舉例2給出序列h0=3,-1,2,-3構(gòu)成的四種類型的系數(shù)向量h1=3,-1,2,-3,5,-3,2,-1,3;h2=3,-1,2,-3,-3,2,-1,3;h3=3,-1,2,-3,0,3,-2,1,-3;h4=3,-1,2,-,3,03,-2,1,-1;分別求出它們的符幅特性曲線，進(jìn)行比較。圖3-1 符幅特性比對(duì)零極點(diǎn)應(yīng)該從FIR濾波器正冪傳遞函數(shù)來判斷。 （3-14）所以它具有與零點(diǎn)數(shù)目相同的M個(gè)極點(diǎn)，它們集中在z平面的原點(diǎn)處，成為M重極點(diǎn)?？梢宰C明，對(duì)于具有對(duì)稱性的系數(shù)向量，如果是上式的分子多項(xiàng)式的根，那么它的倒數(shù)也就是這個(gè)多項(xiàng)式的根。如果是復(fù)數(shù)零點(diǎn)，根據(jù)實(shí)系數(shù)多項(xiàng)式的性質(zhì)，它的共軛復(fù)數(shù)*也是復(fù)數(shù)零點(diǎn)；現(xiàn)在又證明了它的倒數(shù)1/也是零點(diǎn)。因而倒數(shù)的共軛復(fù)數(shù)*也是復(fù)數(shù)零點(diǎn)。四個(gè)復(fù)數(shù)零點(diǎn)必成組出線，圖 3-2 零極點(diǎn)分布圖3.3.5 FIR數(shù)字濾波器的窗函數(shù)設(shè)計(jì)方法宋守鵬.數(shù)字濾波器設(shè)計(jì)及工程應(yīng)用M.江蘇江蘇大學(xué)出版設(shè)，2009.窗函數(shù)設(shè)計(jì)法又稱之為傅里葉技術(shù)設(shè)計(jì)方法，F(xiàn)IR濾波器的單位脈沖響應(yīng)h(n)是有限長的（0nN-1），只有N個(gè)點(diǎn)有值，這N個(gè)值被稱為濾波器的系數(shù)。濾波器的輸出y(n)只與當(dāng)前的輸入x(n)和有限個(gè)過去的輸入x(n-1),x(n-2)x(n-(N-1)有關(guān)，其表達(dá)式為y(n)=h(n)x(n)+h(1)x(n-1) +h(2)x(n-2)+h(N-1)x(n-(N-1)和次表達(dá)式相對(duì)應(yīng)的實(shí)現(xiàn)結(jié)構(gòu)可用圖3-3來表示圖 3-3 實(shí)現(xiàn)窗函數(shù)法的結(jié)構(gòu)圖FIR數(shù)字濾波器的單位沖激響應(yīng)h(n)的Z變換為3 （3-15）利用窗函數(shù)設(shè)計(jì)時(shí)，一般是先給出所要求的理想濾波器的頻率響應(yīng)，要求設(shè)計(jì)一個(gè)FIR濾波器頻率響應(yīng)，來逼近。時(shí)域內(nèi)進(jìn)行的因而先由的傅里葉反變換導(dǎo)出傅里葉系數(shù)即 （3-16）這里的傅里葉系數(shù)就是理想濾波器的單位脈沖響應(yīng)，但是直接應(yīng)用公式（3-15）設(shè)計(jì)數(shù)字濾波器是不可能的，因?yàn)槭菬o限長序列，而且是非因果的。所以要用有限長h(n)來逼近來逼近無限長的，我們可以利用一個(gè)窗函數(shù)序列叫（n）來截取，即 （3-17）然后由h（n）求實(shí)際的，即 （3-18） 然后根據(jù)分析所設(shè)計(jì)出的濾波器的幅頻特性和相頻特性，檢驗(yàn)設(shè)計(jì)結(jié)果是否滿足要求。設(shè)計(jì)過程中，窗函數(shù)將直接影響濾波器的過渡帶帶寬和阻帶衰減。3.3.6 利用頻率取樣法設(shè)計(jì)FIR濾波器以上介紹的窗函數(shù)法是一種時(shí)域設(shè)計(jì)方法，但是在一般情況下，濾波器的技術(shù)指標(biāo)都是由頻率給出，在頻域進(jìn)行設(shè)計(jì)更為直接。頻率取樣法是在頻率域?qū)硐霝V波器取樣，在取樣點(diǎn)上設(shè)計(jì)的濾波幅度值相等，然后根據(jù)濾波器的頻率特性設(shè)計(jì)所需的濾波器。假設(shè)理想濾波器的頻率特性為，對(duì)它在02區(qū)間內(nèi)等間隔的取樣N個(gè)點(diǎn)，即 （3-19）根據(jù)頻域內(nèi)插公式，可得 （3-20）對(duì)應(yīng)的頻率響應(yīng)為 （3-21）其中 （3-22）從（3-17）式中可以看出，根據(jù)理想濾波頻率特性的N個(gè)取樣值，可以通過差值得到實(shí)際設(shè)計(jì)的濾波器的頻率特性，并且所設(shè)計(jì)的濾波器在取樣點(diǎn)上的取值與理想濾波器的值嚴(yán)格相等。這就是頻率取樣法設(shè)計(jì)FIR數(shù)字濾波器的基本原理。在很多情況下要求設(shè)計(jì)的濾波器具有線性相位，前面的討論的濾波器的線性相位條件表述為h(n)= h(N-n-1)。當(dāng)h(n)=h(N-n-1),且N為奇數(shù)時(shí) （3-23）并且幅度函數(shù)H（）應(yīng)該是偶對(duì)稱的，即， （3-25）在0,2pi區(qū)間內(nèi)對(duì)等間隔的取樣N個(gè)點(diǎn)，即 （3-26）其中 （3-27） （3-28）由式（3-26）可得 （3-29）相位條件與N為奇數(shù)時(shí)相同。此時(shí)，奇對(duì)稱要求，即并且當(dāng)h(n)為奇對(duì)稱，按照同樣的方法可以得到幅度和相位應(yīng)滿足的條件。采用頻率取樣法設(shè)計(jì)的FIR數(shù)字濾波器在阻帶內(nèi)的衰減減小，在實(shí)際應(yīng)用中往往達(dá)不到要求。產(chǎn)生這種現(xiàn)象的原因是由于在通帶邊緣取樣點(diǎn)的陡然變換而引起的起伏震蕩。增加阻帶衰減方法是在通帶和阻帶的邊界增加一些過渡的取樣點(diǎn)，從而減小頻帶邊緣的突變，也就減小了起伏震蕩，增大了阻帶最小衰減。為了構(gòu)造特性接近理想低通濾波器的實(shí)際濾波器，應(yīng)當(dāng)使它具備如下特性因果序列有限長度系數(shù)對(duì)稱用長N對(duì)稱中心在對(duì)稱中心在的矩形窗與相乘，截取出實(shí)際系數(shù)序列h(n), （3-30）圖 3-4 實(shí)際系數(shù)序列圖3-5為理想濾波器的脈沖響應(yīng)及其頻譜圖3-5 脈沖響應(yīng)及其頻譜圖3-6 矩形窗函數(shù)及其頻譜矩形窗截取的實(shí)際特性一般濾波器都要求過渡帶愈窄愈好。過渡帶與N成反比，減小過渡帶的措施是增加N；通帶內(nèi)放入波動(dòng)影響濾波器通帶中的信號(hào)通過的均勻性；阻帶內(nèi)的波動(dòng)使某些頻率成分在阻帶內(nèi)不能充分抑制，可能使最小衰減不滿足技術(shù)要求；上述分析說明，增加N對(duì)于通帶波動(dòng)和阻帶衰減幾乎沒有影響。矩形窗帶來的的阻帶衰減為20dB左右，要減小波動(dòng)只有改變窗函數(shù)的形狀。增大N將增加濾波器的復(fù)雜性和計(jì)算時(shí)間。3.3.7 其他幾種窗函數(shù)的介紹 Ingle V K，Proakis J G. Digital Signal Processing Using MATL AB M，New YorkPWS Publishing Company,1997 (1)三角形窗 （3-31）(2)漢寧窗 （3-32）(3)漢明窗 （3-33）Matlab提供的其他窗函數(shù)（程序）N=45;n=1:N;beta=5.658; wdrc=boxcar(N);wdbt=triang(N); wdhn=hanning(N); wdhm=hamming(N);wdbl=blackman(N); wdks=kaiser(N,beta); plot(n,wdrc,wdbt,wdhn,wdhm,wdbl,wdks)Legend(rc,bt,hn,hm,bl,ks)%線性標(biāo)注就可以得到綜合表示幾種窗函數(shù)形狀的圖圖 3-7 窗函數(shù)形狀綜合圖形表3-2 各種窗函數(shù)的過渡帶和波動(dòng) FIR的實(shí)現(xiàn)是非遞歸的，總是穩(wěn)定的；更重要的是，F(xiàn)IR濾波器在滿足幅頻響應(yīng)要求的同時(shí)，可以獲得嚴(yán)格的線性相位特性。線性相位FIR濾波器通常采用窗函數(shù)法設(shè)計(jì)。窗函數(shù)法設(shè)計(jì)FIR濾波器的基本思想是根據(jù)給定的技術(shù)指標(biāo)，選擇濾波器長度N和窗函數(shù)W（n）, 使其具有最窄寬度的主瓣和最小的旁瓣。工程中常用的窗函數(shù)共有五種，即矩形窗。漢寧窗（hanning）。漢明窗（hamming）。布萊克曼窗（blackman）和凱塞窗（kaiser）下面我們運(yùn)用FIR函數(shù)分別運(yùn)用5種不同的窗函數(shù)來設(shè)計(jì)濾波器對(duì)輸入信號(hào)為S(t)=sin(2*pi*50*t)+sin(2*pi*100*t)+sin(2*pi*150*t)進(jìn)行濾波，要求濾波中建的100Hz頻率，采樣頻率為1000Hz。程序 L=30; N=L-1; f1=0,70,80,120,130,500; f2=f1/500; m=0,0,1,1,0,0; b=FIR2(N,f2,m,boxcar(L); figure(1); h,w=freqz(b,1,512,1000); subplot(2,1,1); plot(f1,m,m-,w,abs(h),b-); subplot(2,1,2); plot(w,angle(h); figure(2); freqz(b,1,512,1000); t=0:0.001:1; sgn=sin(2*pi*50*t)+sin(2*pi*100*t)+sin(2*pi*150*t); figure(3); subplot(2,2,1),plot(t(1:100),sgn(1:100); y=fft(sgn,512); Pyy=y.*conj(y)/512; f=1000*(0:155)/512; subplot(2,2,2),plot(f,Pyy(1:156); Yo=filter(b,1,sgn); subplot(2,2,3); plot(t(1:100),Yo(1:100); yo=fft(Yo,512); Pyyo=yo.*conj(yo)/512; fo=1000*(0:155)/512; subplot(2,2,4);plot(f,Pyyo(1:156);程序運(yùn)行后，五種窗函數(shù)的結(jié)果如下(1)矩形窗其相關(guān)參數(shù)，其余。幅度頻率特性過渡帶寬為。旁瓣峰值衰耗為-13dB，最小最帶衰耗為-21dB。運(yùn)行程序結(jié)果如下圖3-8 實(shí)際信號(hào)與理想信號(hào)對(duì)比圖理想結(jié)果與實(shí)際稍有偏差，存在一定失真，這也是無法避免的。圖 3-9濾波器幅相頻特性圖圖 3-10 濾波器前后頻域?qū)Ρ葓D(2)漢寧窗其相關(guān)參數(shù)，其余，幅度頻率特性過渡帶寬為，旁瓣峰值衰耗為31dB,最小阻帶衰耗為-44dB圖 3-11濾波器幅相頻特性圖 圖 3-12 濾波器前后頻域?qū)Ρ葓D由圖可知，沒有完全濾除其他頻率，濾波效果不是很好，可以增加階數(shù)至40.(3)漢明窗其相關(guān)參數(shù)，其余，。幅度頻率特性過渡帶寬為，旁瓣峰值衰耗為-41dB,最小阻帶衰耗為-53dB.圖 3-13濾波器幅相頻特性圖 圖 3-14 濾波器前后頻域?qū)Ρ葓D (4)布萊克曼窗其相關(guān)參數(shù)為，其余，這是一組升余弦窗函數(shù)，升階M是可選參數(shù)。幅度頻率特性過渡帶寬為，旁瓣峰值衰耗為-57dB，最小阻帶衰耗為-75dB.圖 3-15濾波器幅相頻特性圖 圖 3-16 濾波器前后頻域?qū)Ρ葓D (5)凱塞爾窗 相關(guān)參數(shù)，其余，其中為0階貝塞爾函數(shù)，主瓣寬度參數(shù)。圖 3-17 濾波器幅相頻特性圖 圖 3-18 濾波器前后頻域?qū)Ρ葓D由上可以看出，運(yùn)用5種不同窗函數(shù)來對(duì)濾波器進(jìn)行仿真設(shè)計(jì)有著一定的差別。其中，矩形窗函數(shù)和漢寧窗函數(shù)在設(shè)計(jì)上不是很完美，而后的三種傳函數(shù)則效果好很多。3.4 IIR數(shù)字濾波器的設(shè)計(jì)方法IIR數(shù)字濾波器是一種離散時(shí)間系統(tǒng)，其系統(tǒng)函數(shù)為 （3-34）假設(shè)MN，當(dāng)MN時(shí),系統(tǒng)函數(shù)可以看作一個(gè)IIR的子系統(tǒng)和一個(gè)(M-N)的FIR子系統(tǒng)的級(jí)聯(lián)。IIR數(shù)字濾波器的設(shè)計(jì)實(shí)際上是求解濾波器的系數(shù)和 ，它是數(shù)學(xué)上的一種逼近問題，即在規(guī)定意義上（通常采用最小均方誤差準(zhǔn)則）去逼近系統(tǒng)的特性。如果在S平面上去逼近，就得到模擬濾波器；如果在z平面上去逼近，就得到數(shù)字濾波器。利用模擬濾波器來設(shè)計(jì)數(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)式(3-34)中,T是采樣周期。如果令Ha(s)是ha(t)的拉普拉斯變換，H(z)為h(n)的Z變換，利用采樣序列的Z變換與模擬信號(hào)的拉普拉斯變換的關(guān)系得 （3-35）則可看出，脈沖響應(yīng)不變法將模擬濾波器的S平面變換成數(shù)字濾波器的Z平面，這個(gè)從s到z的變換是從S平面變換到Z平面的標(biāo)準(zhǔn)變換關(guān)系式。圖 3-19 S平面轉(zhuǎn)換到Z平面示意圖由（3-35）式，數(shù)字濾波器的頻率響應(yīng)和模擬濾波器的頻率響應(yīng)間的關(guān)系為 （3-36）這就是說，數(shù)字濾波器的頻率響應(yīng)是模擬濾波器頻率響應(yīng)的周期延拓。正如采樣定理所討論的，只有當(dāng)模擬濾波器的頻率響應(yīng)是限帶的，且?guī)抻谡郫B頻率以內(nèi)時(shí)，即 ， （3-37）才能使數(shù)字濾波器的頻率響應(yīng)在折疊頻率以內(nèi)重現(xiàn)模擬濾波器的頻率響應(yīng)，而不產(chǎn)生混疊失真，即 ， （3-38）但是，任何一個(gè)實(shí)際的模擬濾波器頻率響應(yīng)都不是嚴(yán)格限帶的，變換后就會(huì)產(chǎn)生周期延拓分量的頻譜交疊，即產(chǎn)生頻率響應(yīng)的混疊失真，圖3-20。這時(shí)數(shù)字濾波器的頻響就不同于原模擬濾波器的頻響，而帶有一定的失真。當(dāng)模擬濾波器的頻率響應(yīng)在折疊頻率以上處衰減越大、越快時(shí)，變換后頻率響應(yīng)混疊失真就越小。這時(shí)，采用脈沖響應(yīng)不變法設(shè)計(jì)的數(shù)字濾波器才能得到良好的效果。圖 3-20 脈沖不變響應(yīng)示意圖對(duì)某一模擬濾波器的單位沖激響應(yīng)ha(t)進(jìn)行采樣，采樣頻率為fs，若使fs增加，即令采樣時(shí)間間隔（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)，也就是時(shí)域逼近良好，而且模擬頻率和數(shù)字頻率之間呈線性關(guān)系=T。因而，一個(gè)線性相位的模擬濾波器通過脈沖響應(yīng)不變法得到的仍然是一個(gè)線性相位的數(shù)字濾波器。脈沖響應(yīng)不變法的最大缺點(diǎn)是有頻率響應(yīng)的混疊效應(yīng)。所以，脈沖響應(yīng)不變法只適用于限帶的模擬濾波器(例如，衰減特性很好的低通或帶通濾波器)，而且高頻衰減越快，混疊效應(yīng)越小。至于高通和帶阻濾波器，由于它們?cè)诟哳l部分不衰減，因此將