dsp課程設計論文

1、<p><b>　　摘 要</b></p><p>　　本論文首先介紹了濾波器的濾波原理以及數字濾波器的設計方法及過程。重點介紹了FIR數字濾波器的設計方法。即各種窗函數法和等效最佳一致逼近法。在此基礎上，用DSP虛擬實現任意階FIR數字濾波器。此設計擴展性好，便于調節(jié)濾波器的性能，可以根據不同的要求在DSP上加以實現。</p><p>　　濾波是信號處理中

2、最基本又極為重要的技術，利用濾波器技術可以從復雜的信號中提取出所需要的信號，抑制不需要的信號。絕大多數傳感器輸出的信號，在使用過程中，都必須進行濾波。所以濾波器是具有一定傳輸選擇特性的、對信號進行加工處理的裝置、它允許輸入信號中的一些成分通過，抑制或衰減另一些成分。其功能是將輸入信號變換為人們所需要的輸入信號。本論文主要給出了FIR數字濾波器的設計方法和原理。</p><p>　　關鍵字：DSP；濾波；FIR數字

3、濾波器；</p><p><b>　　Abstract</b></p><p>　　This thesis introduced the design method and process of filtering of filter a principle and number filter first. Particularly introduced the d

4、esign method of FIR number filter. Then the various window function method and etc. effect is the best to unanimously approach a method. On this foundation, carry out the rank FIR number with the DSP conjecture the filte

5、r is arbitrarily.This design expands sex well, easy to regulate the function of filter, can take into to carry out on DSP according to th</p><p>　　Filtering the wave is the most basic and extremely and impor

6、tant technique in the signal processing, making use of the filter technique can withdraw the signal that needs from the complicated signal and repress don't need of signal.The great majority spreads the feeling machi

7、ne output's signal, in the process of using in, have to carry on filtering wave.So filter is have to definitely deliver choice characteristic, carry on to the signal process processed device, it to allow importation

8、signal i</p><p>　　Keywords: DSP；Filter wave；FIR number filter；</p><p><b>　　目 錄</b></p><p><b>　　摘 要I</b></p><p>　　AbstractII</p><p&g

9、t;<b>　　目 錄III</b></p><p><b>　　第一章 緒論1</b></p><p>　　1.1 DSP展發(fā)史1</p><p>　　1.2 DSP的基本特征1</p><p>　　1.3 TMS320C54xx的硬件結構2</p><p>　

10、　1.3.1 C54X的CPU體系結構2</p><p>　　1.3.2 指令緩沖單元（I）3</p><p>　　1.3.3 程序流程單元（P）3</p><p>　　1.3.4 地址程序單元（A）3</p><p>　　1.3.5 數據計算單元（D）4</p><p>　　1.4 DSP的設計

11、4</p><p>　　1.4.1 DSP的設計特點4</p><p>　　1.4.2 DSP系統(tǒng)的設計流程5</p><p>　　第二章 FIR數字濾波器的簡介及基本原理6</p><p>　　2.1 數字濾波器的簡介6</p><p>　　2.2 FIR數字濾波器的結構7</p>&l

12、t;p>　　2.3 FIR數字濾波器的特性9</p><p>　　2.3.1 FIR數字濾波器的相位特性9</p><p>　　2.3.2 線性相位FIR數字濾波器的幅度特性10</p><p>　　第三章 基于MATLAB的FIR數字濾波器設計14</p><p>　　3.1 數字濾波器的設計方法描述14</p>

13、<p>　　3.2 常用窗函數及設計方法14</p><p>　　3.3 基于切比雪夫等效一致逼近法19</p><p>　　3.4 FIR數字濾波器的MATLAB設計20</p><p>　　3.4.1 用各種窗函數設計FIR數字濾波器20</p><p>　　3.4.2 用remez函數設計FIR數字低通濾波器2

14、2</p><p>　　3.4.3 應用FDATOOl設計FIR濾波器26</p><p>　　第四章 FIR數字濾波器的DSP實現（CCS仿真）30</p><p>　　4.1 CSS的特點30</p><p>　　4.2 CCS代碼生成工具及程序流程圖30</p><p>　　4.3 CCS集成開發(fā)環(huán)

15、境的配置32</p><p>　　4.4 CCS 集成開發(fā)環(huán)境應用（與MATLAB聯(lián)合仿真）33</p><p>　　第五章 參 考 文 獻42</p><p><b>　　第一章 緒論</b></p><p>　　1.1 DSP發(fā)展史</p><p>　　自20世紀70年代末80年代初

16、DSP芯片誕生以來DSP芯片得到了飛速的發(fā)展。DSP產品的應用己擴大到人們的學習、工作和生活的各個方面，并逐漸成為電子產品更新?lián)Q代的決定因素。DSP發(fā)展歷程大致分為三個階段:20世紀70年代理論先行，80年代產品普及和90年代的突飛猛進。在DSP出現之前數字信號處理只能依靠微處理器(MPU)來完成。但MPU較低的處理速度無法滿足高速實時的要求。因此，直到20世紀70年代，有人才提出了DSP的理論和算法基礎。隨著大規(guī)模集成電路技術的發(fā)展，

17、1978年AMI公司生產出世界上第一片DSP芯片52811。1979年美國Intel公司發(fā)布的商用可編程器2920是DSP芯片的一個重要里程碑。1980年日本NEC公司推出的月PD7720是第一個具有乘法器的商用DSP芯片。幾年后，第二代基于CMOS工藝的DSP芯片應運而生。80年代后期，第三代DSP芯片問世。90年代DSP發(fā)展最快，相繼出現了第四代和第五代DSP器件。在這之后，最成功的DSP芯片當數美國德州儀器公司 (Texas In

18、struments，簡稱TI)的一系列產品，其DSP市場份額占全世界份額近的50%。目前D</p><p>　　在當今的數字化時代，DSP己成為通信設備、計算機和其它電子產品的基礎器件。數字信號處理器與數字信號處理有著密不可分的關系，我們通常說的“DSP”可以指數字信號處理 (Digital Signal Processing)，也可以代表數字信號處理器(Digital signal Processor)在本文里

19、均指數字信號處理器。數字信號處理是一門包括了許多學科并應用于很多領域的學科，是指利用計算機或是專用處理設備，以數字形式對信號進行分析、采集、合成、變換、濾波、估值、壓縮、識別等處理，得到符合要求的信號形式。數字信號處理器是用于處理數字信號的器件，因此它是伴隨著數字信號處理才產生的。</p><p>　　1.2 DSP的基本特征</p><p>　　本設計采用的是TMS320C54x系列DS

20、P，具有改進的哈佛結構、硬件乘法器、流水線結構、高效特殊才旨令集等優(yōu)點，使它的處理速度和容量大大提高，為數字濾波中的復雜算法的實現提供了硬件的保證。</p><p><b>　　1、哈佛總線結構</b></p><p>　　計算機的總線結構可分為兩種。一種是馮·諾依曼結構，其特點是程序和數據共用一個存儲空間。統(tǒng)一編址依靠指令計數器提供的地址來區(qū)分是指令數據還

21、是地址。由于對數據和程序進行分時讀寫，執(zhí)行速度漫。但是半導體工藝的飛速發(fā)展克服了這一缺陷，同時這一結構使計算機結構得到簡化，并成為計算機發(fā)展的一個標準。但由于原理上的特點，這一結構并不適合進行具有高度實時要求的數字信號處理。哈佛結構與馮·諾依曼結構相比，其主要特點是程序和數據具有獨立的存儲空間，有著各自獨立的程序總線和數據總線，大大地提高了數據處理能力，非常適合于實時數字信號處理。TI公司的DSP芯片的結構是改進的哈佛結構，內

22、部采用三個并行的總線，程序讀寫總線、PAB)、讀數據總線 (DABI)、寫數據總線(DABZ)，每組總線連接不同的存儲器空間，使程序存儲器與數據存儲器成為獨立的兩個存儲空間，并且允許數據在程序存儲器和數據存儲器之間進行傳遞，這種結構允許取指令和取操作數并行進行，增加了器件的運算速度。</p><p><b>　　2、硬件乘法器</b></p><p>　　在通用的計算

23、機上，算術邏輯單元(ALU)只能完成兩個操作數的加、減及邏輯運算，而乘法(或除法)則是由加法和移位來實現。因此它們實現乘加運算就比較慢，而在數字濾波算法中需要大量的乘加運算，在TMs320c55xDSP中有一個硬件乘法器，可以在一個指令周期內完成一次乘法和一次加法運算，可以大大提高數字濾波系統(tǒng)的運算速度。</p><p><b>　　3、流水線結構</b></p><p&

24、gt;　　TMS320C54xDSP采用四級流水線指令操作方式，減少了指令的執(zhí)行時間，提高了處理器的吞吐量。任何處理器執(zhí)行一條指令，總是要經過取指令、譯碼、讀操作數和執(zhí)行指令四個階段，需要若干個機器周期才能完成。TMS320C54xDSP系列由于指令和操作數來自不同的空間，同一時刻，可以有四條指令的不同階段在并行處理，盡管每條指令執(zhí)行的時間仍然是幾個機器周期，但由于指令的流水作業(yè)，實現了多條指令的并行執(zhí)行。</p><

25、;p><b>　　4、高效指令集</b></p><p>　　數字信號處理運算的特點之一就是單一運算的重復執(zhí)行，在通用的處理器中，一般采用軟件的方法來解決，基本方法是采用循環(huán)控制或指令重復，循環(huán)控制除了在判斷行環(huán)次數及操作轉向上需要指令開銷外，還會因DSP流水線的頻繁中斷而造成相當大的開銷:指令重復的代價是增加了程序代碼長度，占用更多的存儲空間。而在TM5320C54x DSP中設置了

26、硬件循環(huán)控制電路，提供了重復指令，實現零開銷的循環(huán)控制。所有這些優(yōu)點使得DSP非常適合來做數字濾波器的硬件電路。</p><p>　　1.3 TMS320C54xx的硬件結構 </p><p>　　1.3.1 C54X的CPU體系結構 </p><p>　　C54X有1條32位的程序數據總線（PB），5條16位數據總線（BB、CB、DB、EB、FB）和1條24位的

27、程序地址總線及5條23位地址總線，這些總線分別與CPU相連?？偩€通過存儲單元接口（M）與外部程序總線和數據總線相連，實現CPU對外部存儲器的訪問。這種并行的多總線結構，使CPU能在一個CPU周期內完成1次32位程序代碼讀、3次16位數據讀和兩次16位數據寫。C55X根據功能的不同將CPU分為4個單元，指令緩沖單元（I）、程序流程單元（P）、地址流程單元（A）、和數據計算單元（D）。 </p><p>　　讀程序地

28、址總線（PDA）上傳送24位的程序代碼地址，由讀程序總線（PB）將32位的程序代碼送入指令緩沖單元進行譯碼。 </p><p>　　1.3.2 指令緩沖單元（I）</p><p>　　C54X的指令緩沖單元有指令緩沖隊列IBQ和指令譯碼器組成。在每個CPU周期內，I單元將從程序數據接收的4B程序代碼放入指令緩沖隊列，指令譯碼器從隊列中取6B程序代碼，根據指令的長度可對8位、16位、24位

29、、32位和48位的變長指令進行譯碼，然后把譯碼數據送入P單元、A單元和D單元去執(zhí)行。 </p><p>　　1.3.3 程序流程單元（P）</p><p>　　程序流程單元有程序地址產生電路和寄存器組湊成。程序流程單元產生所有程序空間的地址，并控制指令的讀取順序。</p><p>　　程序地址產生邏輯電路的任務是產生讀取空間的24位地址。一般情況下，它產生的是連續(xù)

30、地址，如果指令要求讀取非連續(xù)地址的程序代碼時，程序地址產生邏輯電路能夠接收來自I單元的立即數和來自D單元的寄存器值，并將產生的地址傳送到PAB。 </p><p>　　在P單元中使用的寄存器分為5種類型： </p><p>　　a.程序流寄存器：包括程序計數器、返回地址寄存器和控制流程關系寄存器。 </p><p>　　b.塊重復寄存器：包括塊重復寄存器0和1（BR

31、C0、BRC1）BRC1的保存寄存器(BRS1)、塊重復起始地址寄存器0和1以及塊重復結束地址寄存器0和1。 </p><p>　　c.單重復寄存器：包括單重復寄存器和計算單重復寄存器。 </p><p>　　d.中斷寄存器：包括中斷標志寄存器0和1、中斷使能寄存器0和1以及調試中斷使能寄存器0和1。 </p><p>　　e.狀態(tài)奇存期：包括狀態(tài)寄存器0，1，2和

32、3。 </p><p>　　1.3.4 地址程序單元（A） </p><p>　　地址程序單元包括數據地址產生電路、算術邏輯電路和寄存器組構成。 </p><p>　　數據地址產生電路能夠接收來自I單元的立即數和來自A單元的寄存器產生讀取數據空間的地址。對于使用間接尋址模式的指令，有P單元向DAGEN說明采用的尋址模式。 </p><p>

33、　　A單元包括一個16位的算術邏輯單元，它既可以接收來自I單元的立即數也可以與存儲器、I/O空間、A單元寄存器、D單元寄存器和P單元寄存器進行雙向通信。 </p><p>　　A單元包括的寄存器有以下幾種類型： </p><p>　　a.數據頁寄存器：包括數據頁寄存器和接口數據頁寄存器； </p><p>　　b.指針：包括系數數據指針寄存器、堆棧針寄存器和8個輔助

34、寄存器； </p><p>　　c.循環(huán)緩沖寄存器：包括循環(huán)緩沖大小寄存器、循環(huán)緩沖起始地址寄存器; </p><p>　　d.臨時寄存器：包括臨時寄存器。 </p><p>　　1.3.5 數據計算單元（D）</p><p>　　數據計算單元由移位器、算數邏輯電路、乘法累加器和寄存器組構成。D單元包含了CPU的主要運算部件。 </p

35、><p>　　D單元移位器能夠接收來自I單元的立即數，能夠與存儲器、I/O單元、A單元寄存器、D單元寄存器和P單元寄存器進行雙向通信，此外，還可以向D單元的ALU和A單元的ALU提供移位后的數據。移位可以完成以下操作： </p><p>　　a.對40位的累加器可以完成向左最多32位的移位操作，移位數樂意從零食寄存器讀取或由指令中的常數提供； </p><p>　　b.

36、對于16位寄存器、存儲器或I/O空間數據可完成左移31位或32位的移位操作； </p><p>　　c.對于16位立即數可完成向左移最多15位的移位操作。 </p><p>　　1.4 DSP的設計</p><p>　　1.4.1 DSP的設計特點</p><p>　　DSP系統(tǒng)是以數字信號處理為基礎的，因此不但具有數字處理的全部優(yōu)點而且

37、還具有以下特點：</p><p>　　a.接口方便:DSP應用系統(tǒng)與其他以現代；之字技術為基礎的系統(tǒng)或設備都是相互兼容的，這樣的系統(tǒng)接口以實現某種功能要比模擬系統(tǒng)與這些系統(tǒng)接口要容易得多。</p><p>　　b.編程方便：DSP應用系統(tǒng)中的可編程DSP芯片，能靈活方便地進行修改和升級。</p><p>　　c.穩(wěn)定性好：DSP應用系統(tǒng)以數字處理為基礎，受環(huán)境溫度及

38、噪聲的影響較小、可靠性高，無器件老化現象。</p><p>　　d.精度高：16位數字系統(tǒng)可以達到10一級的精度。</p><p>　　e.可重復性好：模擬系統(tǒng)的性能受元器件參數性能變化的影響比較大，而數字系統(tǒng)基本不受影響，因此數字系統(tǒng)便于測試、調試和大規(guī)模生產。</p><p>　　f.集成方便:DSP應用系統(tǒng)中的數字部件有高度的規(guī)范性，便于大規(guī)模集成。</

39、p><p>　　當然，數字信號處理也存在一些缺點。例如，對于簡單信號處理任務，若采用DSP則使成本增加。DSP系統(tǒng)中的高速時鐘可能帶來高頻干擾和電磁泄漏等問題，而且DSP系統(tǒng)消耗的功率也較大。此外，DSP技術更新速度快，對于數學知識要求高，開發(fā)和測試工具還有待進一步完善。</p><p>　　1.4.2 DSP系統(tǒng)的設計流程</p><p>　　一個DSP系統(tǒng)的設計過

40、程大概要有以下幾個步驟：</p><p>　　a.根據系統(tǒng)的任務要求，確定系統(tǒng)處理精度要求、速度要求、實時性要求等性能指標。</p><p>　　b.根據系統(tǒng)的要求進行高級語言的算法模擬，比如使用MATLAB等仿真工具，驗證算法的可行性，得出最佳的處理方法。</p><p>　　c.DSP的系統(tǒng)設計，主要分為硬件設計和軟件設計。硬件設計是指根據系統(tǒng)要求選擇合適的DS

41、P芯片然后設計相應的外圍電路。軟件設計主要是指根據系統(tǒng)的要求和選用的DSP芯片編寫相應的程序。程序的編寫可以使用匯編語言，匯編語言編寫的程序效率高，但比較煩雜；也可采用C語言，DSP的C語言基本上是標準C語言，編寫比較簡單，但效率低。在實際系統(tǒng)開發(fā)時往往是兩種語言結合編寫，在算法運算量大的地方使用匯編語言，在運算量小的地方使用C語言，這樣既能縮短軟件的開發(fā)周期，提高程序的可讀性和可移植性，又滿足了系統(tǒng)的實時性要求。</p>

42、<p>　　第二章 FIR數字濾波器的簡介及基本原理</p><p>　　2.1 數字濾波器的簡介</p><p>　　隨著信息時代和數字世界的到來，數字信號處理已成為一門極其重要的學科和技術領域。數字信號處理在通信、語音、圖像、自動控制、雷達、軍事、航空航天、醫(yī)療和家用電器等眾多領域得到了廣泛的應用。數字信號處理（DSP）包括兩重含義：數字信號處理技術（Digital Sig

43、nal Processing ）和數字信號處理器（ Digital Signal Processor ）。數字信號處理（DSP）是利用計算機或專用處理設備，以數值計算的方法、對信號進行采集、濾波、增強、壓縮、估值和識別等加工處理，借以達到提取信息和便于應用的目的，其應用范圍涉及幾乎所有的工程技術領域。</p><p>　　在數字信號處理的應用中，數字濾波器很重要而且得到了廣泛的應用。按照數字濾波器的特性，它可以被

44、分為線性與非線性、因果與非因果、無限長沖擊響應（IIR）與有限長沖擊響應（FIR）等等。其中，線性時不變的數字濾波器是最基本的類型；而由于數字系統(tǒng)可以對延時器加以利用，因此可以引入一定程度的非因果性，獲得比傳統(tǒng)的因果濾波器更靈活強大的特性；IIR 濾波器的特征是具有無限持續(xù)時間沖激響應，這種濾波器一般需要用遞歸模型來實現，因而有時也稱之為遞歸濾波器，而FIR濾波器的沖激響應只能延續(xù)一定時間，在工程實際中可以采用遞歸的方式實現，也可以采用

45、非遞歸的方式實現，但其結構主要還是是非遞歸結構，沒有輸出到輸入的反饋，并且FIR濾波器很容易獲得嚴格的線性相位特性，避免被處理信號產生相位失真，而線性相位體現在時域中僅僅是h( n)在時間上的延遲，這個特點在圖像信號處理、數據傳輸等波形傳遞系統(tǒng)中是非常重要的，且不會發(fā)生阻塞現象，能避免強信號淹沒弱信號，因此特別適合信號強弱相差懸殊的情況。相對于IIR濾波器，FIR濾波器有著易于實現和系統(tǒng)絕對穩(wěn)定的優(yōu)勢，因此得到廣泛的應用；對于時變系統(tǒng)濾

46、波器的研究則導致了以卡爾曼濾</p><p>　　一個數字濾波器可以用系統(tǒng)函數表示為：</p><p>　　由此式可得到表示輸入輸出關系的常系數線性差分方程為：</p><p>　　可見數字濾波器的功能就是把輸入序列x(n)通過一定的運算變換成輸出序列y(n)。不同的運算處理方法決定了濾波器實現結構的不同。</p><p>　　數字濾波器的運

47、算結構有兩種表示方法：方框圖和信號流圖法，如圖所示：最常見的3個基本運算單元：加法器、單位延時器和常數乘法器。</p><p>　　信號流圖 方框圖</p><p><b>　　單位延時 </b></p><p><b>　　乘常數 </b></p>

48、;<p><b>　　相加 </b></p><p>　　研究濾波器實現結構的意義：</p><p>　?。?）濾波器的基本特性，如有線長沖激響應與無限長沖激響應，決定了結構上有不同的特點；</p><p>　?。?）不同結構所需的存儲單元及乘法次數不同，前者影響復雜性，后者影響運算速度；</p><

49、;p>　?。?）在有限精度（有限字長）實現情況下，不同運算結構的誤差及穩(wěn)定性不同；</p><p>　?。?）好的濾波器結構應該易于控制濾波器性能，適合于模塊化實現，便于時分復用；</p><p>　　2.2 FIR數字濾波器的結構</p><p>　　設 h ( n) ( n = 0，1， 2 ?N - 1)為濾波器的沖激響應，輸入信號為 x ( n) ，

50、則 FIR 濾波器就是要實現下列差分方程: 數字濾波器具有一下差分方程：</p><p><b>　　(1)</b></p><p>　　式中，為輸入序列，為輸出序列，和為濾波器系數，N是濾波器階數。若所有的均為0，則有：</p><p><b>　　(2)</b></p><p>　　(2)式即為

51、FIR的差分方程。為了不失一般性，用下式來表示FIR的濾波器差分</p><p><b>　　方程：</b></p><p><b>　　(3)</b></p><p>　　將(2)式進行Z變換，整理后可得FIR濾波器的傳遞函數：</p><p><b>　　(4)</b>&l

52、t;/p><p>　　FIR濾波器實質上就是一個分節(jié)的延遲線，把每一節(jié)的輸出用濾波器系數進行加權累加，便得到濾波器的輸出結果，它總是穩(wěn)定并且可實現的。在一些工程實際應用（如：圖像處理、數據調制解調）中，往往對相位要求較高。FIR濾波器可以實現嚴格的線性相位，從而得到了廣泛應用。它的差分方程數學表達式為：</p><p><b>　　(5)</b></p>&

53、lt;p>　　(5)式中，N是FIR濾波器的抽頭數，x(n)表示在n 時刻輸入的信號樣值，h(n)表示濾波器的第n級抽頭系數。橫截型FIR濾波器的結構如下圖所示：</p><p>　　(FIR濾波器的結構圖)</p><p>　　2.3 FIR數字濾波器的特性</p><p>　　2.3.1 FIR數字濾波器的相位特性</p><p>

54、　　IIR數字濾波器能夠保留一些模擬濾波器的優(yōu)良特性，比如具有良好的幅頻特性，但是其相位是非線性的。FIR數字濾波器可以設計成嚴格線性相位的，避免被處理信號產生相位失真。</p><p>　　FIR數字濾波器設計就是用多項式： （1） </p><p>　　來逼近所要求的頻率特性指標。由于它的單位沖激響應是有限長的，所以FIR數字濾波器是穩(wěn)定的。由式（1）可以得到FIR數字

55、濾波器的頻率響應：</p><p><b>　　（2）</b></p><p>　　其中，是幅頻特性，是相頻特性。如果要求FIR數字濾波器具有嚴格線性相位，即相位不失真時，其相位和頻率呈正比，即相頻特性滿足：</p><p><b>　?。?） </b></p><p>　　其中，為群延時。式（3）

56、說明系統(tǒng)對信號中所有頻率分量都具有相同的時間延遲。對上述條件降低一點的要求是相位和頻率呈線性關系，即</p><p><b>　　（4）</b></p><p>　　雖然的存在使相位呈非線性，但是它的群延時仍保持常數。</p><p>　　FIR數字濾波器的沖激響應是實數，當h(n)是偶對稱的，即</p><p><

57、;b>　?。?）</b></p><p>　　其對稱中心為： （6）</p><p><b>　　根據式（2），得到</b></p><p><b>　　（7）</b></p><p>　　其中，幅頻特性為：

58、 （8）</p><p>　　相頻特性為： （9）</p><p>　　滿足式（3）的條件。</p><p>　　當h(n)為奇對稱時，即 （10）</p><p>　　其對稱中心為。同理，可得其幅度特性為：</p><p

59、><b>　　（11）</b></p><p>　　相頻特性為： （12）</p><p>　　滿足式（4）的條件。</p><p>　　綜上所述，FIR數字濾波器具有線性相位的充要條件是：</p><p>　　或者 （13）</p><p>　　2.3.

60、2 線性相位FIR數字濾波器的幅度特性</p><p>　　如果濾波器的系數h(n)的長度為N，且這些系數是關于對稱的，根據h(n)的奇偶對稱性和N的奇偶性，線性相位FIR數字濾波器可以分為4種類型，下面分別介紹這4種類型濾波器的頻率響應。</p><p>　　1）I型濾波器，系數h(n)為偶對稱，N為奇數</p><p>　　當系數h(n)為偶對稱，N為奇數時，根

61、據式（8），該類型濾波器的幅度特性函數為： （14）</p><p>　　其中， （15）</p><p>　　濾波器的幅度響應對呈偶對稱。</p><p>　　2）II型濾波器，系數h(n)為偶對稱，N為偶數</p><p>　　當系數h(n)為偶對稱，N為偶數時

62、，根據式（8），該類型濾波器的幅度特性函數為： （16）</p><p>　　其中， （17）</p><p>　　濾波器的幅度響應對呈奇對稱。但是，由于在時等于零，不能用這種方式實現在有頻率響應的頻率特性，比如高通濾波器和帶阻濾波器。</p><p>　　3）III

63、型濾波器，系數h(n)為奇對稱，N為奇數</p><p>　　當系數h(n)為奇對稱，N為奇數時，根據式（11），該類型濾波器的幅度特性函數為： （18）</p><p>　　其中， （19）</p><p>　　濾波器的幅度響應對奇對稱。但是，由于在時等于零，不能用

64、這種方式實現低通濾波器、高通濾波器和帶阻濾波器，只能用做帶通濾波器。</p><p>　　4）IV型濾波器，系數h(n)為奇對稱，N為偶數</p><p>　　當系數h(n)為奇對稱，N為偶數時，根據式（11），該類型濾波器的幅度特性</p><p>　　函數為： （20）</p><p>　　其中，

65、 （21）</p><p>　　濾波器的幅度響應對呈奇對稱，對呈偶對稱。但是，由于在時等于零，不能用這種方式實現低通濾波器和帶阻濾波器。</p><p>　　2.3.3 線性相位FIR數字濾波器的零極點特性</p><p>　　FIR數字濾波器的零點是其系數多項式的根，它的極點與原點數目相同，集中在Z平面的

66、原點處。</p><p>　　由于線性相位FIR數字濾波器的單位脈沖響應具有對稱性，即，可得 （22）</p><p>　　由上式可以看出，如果是該濾波器的一個零點，則也是它的零點。又由于h(n)是實數，H(z)的零點必定共軛成對出現，則和也是零點。所以，線性相位FIR數字濾波器的零點必是互為倒數的共軛對。</p><p>　　根據4種類型線性相位F

67、IR數字濾波器的特點，可以得到它們零點特性的主要區(qū)別是在z=1處和z=-1處的零點數量，即</p><p>　　1）I型線性相位FIR數字濾波器在z=1和z=-1處有偶數個零點或者沒有零點。</p><p>　　2）II型線性相位FIR數字濾波器在z=1處有偶數個零點或者沒有零點，在z=-1處有奇數個零點。</p><p>　　3）III型線性相位FIR數字濾波器在

68、z=1和z=-1處有奇數個零點。</p><p>　　4）IV型線性相位FIR數字濾波器在z=1處有奇數個零點，在z=-1處有偶數個零點或者沒有零點。</p><p>　　例：根據給出4種類型濾波器的系數，分別畫出其零極點圖。</p><p>　　h1=[-4,1,-1,-2,5,6,5,-2,-1,1,-4]; h2=[-4,1,-1,-2,5,6,6,5,-2,

69、-1,1,-4];</p><p>　　h3=[-4,1,-1,-2,5,0,-5,2,1,-1,4]; h4=[-4,1,-1,-2,5,6,-6,-5,2,1,-1,4];</p><p>　　clear all; close all; clc;</p><p>　　h1=[-4,1,-1,-2,5,6,5,-2,-1,1,-4];h2=[-4,1,-1,

70、-2,5,6,6,5,-2,-1,1,-4];</p><p>　　h3=[-4,1,-1,-2,5,0,-5,2,1,-1,4];h4=[-4,1,-1,-2,5,6,-6,-5,2,1,-1,4];</p><p>　　subplot(2,2,1);zplane(h1,1);title('I型零極點')</p><p>　　subplot(2,2

71、,2);zplane(h2,1);title('II型零極點')</p><p>　　subplot(2,2,3);zplane(h3,1);title('III型零極點')</p><p>　　subplot(2,2,4);zplane(h4,1);title('IV型零極點')</p><p>　　第三章 基于MA

72、TLAB的FIR數字濾波器設計</p><p>　　3.1 數字濾波器的設計方法描述</p><p>　　數字濾波器有多種設計方法，如雙線性變換法、窗函數設計法、插值逼近法和Chebyshev逼近法等等，但是通常采用窗函數法設計。窗函數法設計FIR濾波器的基本思想是：根據給定的濾波器技術指標選擇濾波器長度N和窗函數，使其具有最窄寬度的主瓣和最小的旁瓣，其核心是從給定的頻率特性，通過加窗確定

73、有限長單位脈沖響應序列h(n)。一般采用以下五種窗函數：</p><p>　　矩形窗(Rectangular Window)、三角形窗(Triangular Window)、漢寧窗（Hanning Window）、哈明窗（Hamming Window）、布萊克曼窗(Blackman Window）。</p><p>　　目前FIR濾波器的實現方法大致可分為三種：利用單片通用數字濾波器集成電

74、路、DSP器件和可編程邏輯器件實現。單片通用數字濾波器使用方便，但由于字長和階數的規(guī)格較少，不能完全滿足實際需要，使用以串行運算為主導的通用DSP芯實現要簡單，是一種實時、快速、特別適合于實現各種數字信號處理運算的微處理器，借助于通用數字計算機按濾波器的設計算法編出程序進行數字濾波計算 。由于它具有豐富的硬件資源、改進的哈佛結構、高速數據處理能力和強大的指令系統(tǒng)，因此在通信、航空、航天、雷達、工業(yè)控制網絡及家用電器等各個領域得到廣泛應用

75、。</p><p>　　3.2 常用窗函數及設計方法</p><p><b>　　1、矩形窗</b></p><p>　　矩形窗的時域表達式為： （23）</p><p>　　它的頻域表達式為： （24）</p><p>　　矩形窗的主

76、瓣寬度為,第一旁瓣比主瓣低13dB，阻帶最小衰減為21dB。</p><p>　　在MATLAB中，矩形窗函數為boxcar,調用格式為： w=boxcar(N)</p><p>　　其中，N是窗函數的長度；返回值w是一個長度為N的矩形窗序列。</p><p><b>　　2、三角窗</b></p><p>　　三角窗的

77、時域表達式為以下幾種。</p><p>　　窗長為奇數時： （25）</p><p>　　窗長為偶數時： （26）</p><p>　　它的頻域表達式為： （27）</p><p>　　三角窗的主瓣寬度為,第一旁瓣比主瓣低26dB，阻帶最小衰減為25dB。&

78、lt;/p><p>　　在MATLAB中，三角窗函數為triang,調用格式為： w=triang(N)</p><p>　　其中，N是窗函數的長度；返回值w是一個長度為N的三角窗序列。</p><p><b>　　3、漢寧窗</b></p><p>　　漢寧窗函數又稱升余弦函數，是余弦平方函數。</p>&l

79、t;p>　　它的時域表達式為： （28）</p><p><b>　　它的頻域表達式為：</b></p><p><b>　?。?9）</b></p><p>　　其中，為矩形窗的幅度函數。漢寧窗的主瓣寬度為,第一旁瓣比主瓣低31dB，阻帶最小衰減為44dB。</p><p>

80、;　　在MATLAB中，漢寧窗函數為hanning,調用格式為： w=hanning(N)</p><p>　　其中，N是窗函數的長度；返回值w是一個長度為N的漢寧窗序列。</p><p><b>　　4、哈明窗</b></p><p>　　哈明窗函數是一種改進的升余弦函數，或者稱為改進的漢寧窗。</p><p>　　它

81、的時域表達式為： （30）</p><p>　　哈明窗在保證與漢寧窗相同的主瓣寬度的條件下，使大約99.96%的能量集中在主瓣，第一旁瓣比主瓣低41dB，阻帶最小衰減為53dB。</p><p>　　在MATLAB中，哈明窗函數為hamming,調用格式為： w=hamming(N)</p><p>　　其中，N是窗函數的長度；返回值w是一個長

82、度為N的哈明窗序列。</p><p><b>　　5、布萊克曼窗</b></p><p>　　布萊克曼窗函數又稱三階升余弦函數，它對升余弦函數加上一個二次諧波余弦分量，從而進一步降低旁瓣。</p><p>　　它的時域表達式為：（31）</p><p><b>　　它的頻域表達式為：</b><

83、/p><p><b>　?。?2）</b></p><p>　　布萊克曼窗的主瓣寬度為,第一旁瓣比主瓣低57dB，阻帶最小衰減為74dB。</p><p>　　在MATLAB中，布萊克曼窗函數為blackman,調用格式為： w=blackman(N)</p><p>　　其中，N是窗函數的長度；返回值w是一個長度為N的布萊

84、克曼窗序列。上述各種窗函數的性能比較如表所示: </p><p>　　6、窗函數的MATLAB實現</p><p>　　繪出長度為64的矩形窗、三角窗、漢寧窗、哈明窗和布萊克曼窗函數的時域和頻域幅度特性曲線。</p><p>　　程序代碼如下：clear all;close all;clc;</p><p>　　N=64; w1=b

85、oxcar(N);w2=triang(N);</p><p>　　w3=hanning(N);w4=hamming(N);w5=blackman(N);</p><p>　　wvtool(w1);wvtool(w2);wvtool(w3);</p><p>　　wvtool(w4);wvtool(w5);</p><p><b>　

86、　結果如下圖所示：</b></p><p><b>　　矩形窗</b></p><p><b>　　三角窗</b></p><p><b>　　漢寧窗</b></p><p><b>　　哈明窗</b></p><p>

87、<b>　　布萊克曼窗</b></p><p>　　7、 MATLAB信號處理工具箱提供了基于窗函數法的FIR濾波器的設計函數fir1和fir2,它們能使濾波器的設計更加簡單。</p><p><b>　　1）fir1</b></p><p>　　功能：基于窗函數的FIR濾波器設計——標準頻率響應形狀。</p>

88、<p>　　格式：b=fir1(n,Wn,’ftype’,window)</p><p>　　其中，b為濾波器系數；n為濾波器階數；Wn為截止頻率，0<Wn<1，1對應于采樣頻率的一半；ftype為當指定ftype時可設計高通和帶阻濾波器，ftype=high時設計高通FIR濾波器，ftype=stop時設計帶阻FIR濾波器，低通和帶通FIR濾波器無需輸入ftype參數；window為窗

89、函數，窗函數的長度應等于FIR濾波器系數的個數，即階數n+1。</p><p>　　說明：標準頻率響應應指所設計的濾波器的預期特性為理想頻率響應，包括低通、帶通、高通或帶阻特性。</p><p><b>　　2）fir2</b></p><p>　　功能：基于窗函數的FIR濾波器設計——任意頻率響應形狀。</p><p>

90、;　　格式：b=fir2(n,f,m,window)</p><p>　　其中，b為濾波器系數；n為濾波器階數；f為頻率點矢量，0<f<1，f=1時對應的頻率為采樣頻率的一半，矢量f按升序排列，且第一個必須為0，最后一個必須為1，并允許出現相同的頻率值；m為幅度矢量，按長期共存矢量包含與f相對應的期望濾波器響應幅度，矢量f與矢量m的長度必須相同；window為窗函數，窗函數的長度應等于FIR濾波器系數

91、的個數，即階數n+1。</p><p>　　說明：fir2函數用于設計具有任意頻率響應形狀的加窗線性相位FIR數字濾波器，其幅頻特性由頻率點向量f和幅度值向量m給出，0 f 1,要求f為單增向量，而且從0開始，以1結束，1表示數字頻率w= .m與f等長度，m(k)表示頻點f(k)的幅頻響應曲線。</p><p>　　3.3 基于切比雪夫等效一致逼近法</p><p>

92、;　　工具箱采用remez算法實現線性相位FIR數字濾波器的等波紋最佳一致逼近法設計。與其他設計法相比，其優(yōu)點是，設計指標相同時，使濾波器階數最低，或者階數相同時，使通帶最平坦，阻帶最小，衰減最大；通帶和阻帶均為等波紋形式，最適合設計片段常數特性的濾波器。其調用格式如下：</p><p>　　b=remez(n,f,m,w,’ftype’)</p><p>　　其中，w和ftype可默認。

93、b為濾波器系數向量，調用參數n,f,m的含義與函數fir2中類同，但這里有一點不同，期望逼近的頻幅響應值位于f(k)與f(k+1)（k為奇數）之間的頻段上，而f(k+1)與f(k+2)之間為無關區(qū)。w為加權向量，其長度為f的一半。w(k)為對m中第k個常數片段的逼近精度加權值，w值越大逼近精度越高。ftype用于指定濾波器類型。</p><p>　　Remezord函數用于估算FIR數字濾波器的等波紋最佳一致逼近

94、設計的最低階數n，從而使濾波器在滿足指標的前提下造價最低?；菊{用格式如下：</p><p>　　[n,fo,mo,w]=remezord(f,m,dev,Fs)</p><p>　　其返回參數供remez函數使用。設計的濾波器可以滿足由參數f,m,dev和Fs指定的指標。F和m與remez中所用的類似，這里f可以是模擬頻率(Hz)或歸一化數字頻率，但必須以0開始，以Fs/2(用歸一化頻率

95、時為1)結束，而且其中省略了0和Fs/2兩個頻點。Fs為采樣頻率，省略時默認為2Hz。</p><p>　　dev為各逼近頻段允許的幅頻響應偏差（波紋振幅）。</p><p>　　remez函數可直接調用remezord返回的參數，使用格式如下：</p><p>　　b=remez(n,fo,mo,w)。</p><p>　　3.4 FIR

96、數字濾波器的MATLAB設計</p><p>　　MATLAB設計的程序流程圖： </p><p>　　3.4.1 用各種窗函數設計FIR數字濾波器</p><p>　　1）窗函數設計FIR數字濾波器的步驟</p><p>　　實際的濾波器設計指標通常包括通帶截止頻率、阻帶起始頻率、通帶最大衰減和阻帶最小衰減。其步驟如下：</p>

97、;<p>　　（1）構造希望逼近的頻率響應函數。</p><p>　?。?）根據阻帶最小衰減選定窗函數類型，由“過渡帶寬=窗函數主瓣寬度”確定窗長N。為保證所設計的FIR數字濾波器具有第一類線性相位，一般選取N為奇數。根據N寫出窗函數表達式。</p><p>　?。?）求出理想沖激響應，理想濾波器的截止頻率。</p><p>　　（4）求出實際濾波器的

98、沖激響應：，其中,。</p><p>　　2）例題： 分別用矩形窗、漢寧窗、哈明窗和布萊克曼窗設計FIR數字低通濾波器，要求：通帶截止頻率，窗長N=21。</p><p>　　程序代碼： clear all; close all; clc;</p><p>　　wd=0.2*pi;N=21;M=(N-1)/2;</p><p>　　n

99、n=-M:M;n=nn+eps;hd=sin(wd*n)./(pi*n); </p><p>　　w1=boxcar(N)';h1=hd.*w1; w2=hanning(N)';h2=hd.*w2;</p><p>　　w3=hamming(N)';h3=hd.*w3; w4=blackman(N)';h4=hd.*w4;</p>&l

100、t;p>　　H1=20*log10(abs(fft(h1,1024)));H2=20*log10(abs(fft(h2,1024)));</p><p>　　H3=20*log10(abs(fft(h3,1024)));H4=20*log10(abs(fft(h4,1024)));</p><p>　　HH1=[H1(513:1024) H1(1:512)];</p>

101、<p>　　HH2=[H2(513:1024) H2(1:512)];</p><p>　　HH3=[H3(513:1024) H3(1:512)];</p><p>　　HH4=[H4(513:1024) H4(1:512)];</p><p>　　w=(-512:511)/512;plot(w,HH1,w,HH2,':',w,HH3,&

102、#39;-',w,HH4,'--');</p><p>　　axis([-1.2 1.2 -150 20]);</p><p>　　legend('矩形窗','漢寧窗','哈明窗','布萊克曼窗');</p><p>　　xlabel('\omega/pi');s

103、et(gcf,'color','w'); </p><p>　　4種窗函數設計的低通濾波器分貝幅度特性</p><p>　　從圖中可以看出，用矩形窗設計時得到的過渡帶最窄，但阻帶衰減最小；布萊克曼窗設計時得到的過渡帶最寬，但換來的是阻帶衰減最大。通常情況下，濾波器設計的兩個重要指標很難同時滿足，設計過程中應該折中考慮。</p>

104、<p>　　調用fir2函數設計一個FIR數字低通濾波器：其在0到pi/8的幅度響應為1，在pi/8到2pi/8幅度響應為1/2，在2pi/8到4pi/8幅度響應為1/4，在4pi/8到6pi/8幅度響應為1/6，在6pi/8到pi幅度響應為1/8，并且濾波器的階數為60。畫出理想濾波器和設計得到的濾波器的幅度頻率響應進行比較。程序代碼如下： </p><p>　　f=[0 0.125 0.125 0.

105、25 0.25 0.5 0.5 0.75 0.75 1];</p><p>　　m=[1 1 0.5 0.5 0.25 0.25 1/6 1/6 0.125 0.125]; b=fir2(60,f,m);</p><p>　　[h,w]=freqz(b);plot(f,m,w/pi,abs(h));grid on; </p><p>　　legend('‘理

106、想濾波器',’設計濾波器'); 結果如下圖所示：</p><p>　　3.4.2 用remez函數設計FIR數字低通濾波器</p><p>　　1）設計濾波器，使逼近低通濾波特性 。要求通帶波紋 ，阻帶衰減 ，并用最小階數實現。繪出設計的FIR數字濾波幅頻特性曲線。</p><p><b>　　程序代碼： </b>&

107、lt;/p><p>　　%用remez函數設計FIR低通濾波器</p><p>　　clear;close all; fc=1/4;fs=5/16; %輸入給定指標 Rp=3;As=60;Fs=2; f=[fc,fs];m=[1,0]; dev=[(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-As/20)];</p><p>　　%計算re

108、mezord函數所需要參數f,m,dev</p><p>　　[N,fo,mo,w]=remezord(f,m,dev,Fs);  %確定remez函數所需參數</p><p>　　hn=remez(N,fo,mo,w); %調用remez函數進行設計</p><p>　　hw=fft(hn,512);   

109、0;       %求設計出的濾波器頻率特性</p><p>　　w=[0:511]*2/512;plot(w,20*log10(abs(hw)));grid; %畫對數幅頻特性曲線</p><p>　　axis([0,max(w)/2,-90,5]);xlabel('w/pi'); ylabel('Ma

110、gnitude(dB)');</p><p>　　title('幅頻特性'); line([0,0.4],[-3,-3]);    %畫線檢驗設計結果</p><p>　　line([1/4,1/4],[-90,5]); line([5/16,5/16],[-90,5]);</p><p><b>　

111、　結果如圖所示：</b></p><p>　　結果分析：如圖所示，圖中橫線為-3dB，兩條豎線分別位于頻率 1/4和5 /16。顯然，通帶指標稍有富裕，過渡帶寬度和阻帶最小衰減剛好滿足指標要求。</p><p>　　2）用remez函數設計FIR數字濾波器（在濾波器設計中，技術指標越高，實現濾波器的階數也就越高，在remez函數調用格式b=remez(N,f,m,w)。）<

112、;/p><p>　　程序代碼如下：%用remez函數設計數字低通濾波器</p><p>　　Fs=1000;num=8192;t=(1:10000)/Fs;</p><p>　　s1=sin(2*pi*t*100);s2=sin(2*pi*t*150);s3=sin(2*pi*t*200);</p><p>　　s=s1+s2+s3;</p

113、><p>　　figure(1);plot(t,s);axis([0,0.1,-3,3]);grid;</p><p>　　xlabel('time/s');ylabel('amplitude');</p><p>　　%設計一個低通濾波器 </p><p>　　fedge=[110 140];mval=[1 0]

114、;dev=[0.0559 0.01];</p><p>　　[N,fpts,mag,wt]=remezord(fedge,mval,dev,Fs);</p><p>　　b=remez(N,fpts,mag,wt);disp(b);</p><p>　　[h,w]=freqz(b,1,256);</p><p>　　figure(2);plo

115、t(w/pi,20*log10(abs(h)));grid;</p><p>　　xlabel('\omega/pi');ylabel('Gain(dB)');</p><p><b>　　%對信號進行濾波</b></p><p>　　sf=filter(b,1,s);</p><p>