语音信号的数字滤波处理(五).docVIP

目 錄 1 緒論 1 1.1 課程設計背景 1 1.2 課程設計目的 1 2 課程設計預習 2 2.1 卷積運算的演示 2 2.2 採樣定理的演示 4 3 程式設計的基本原理 6 3.1 IIR濾波器的設計原理 7 3.2 FIR濾波器的設計原理 7 4 設計程式的調試和運行結果 8 4.1 濾波器的設計 9 4.1.1設計方案 10 4.1.2設計步驟 11 4.2 IIR濾波器的設計 12 4.2.1巴特沃斯低通程式的設計 13 4.2.2巴特沃斯高通程式的設計 14 4.2.3巴特沃斯帶通程式的設計 15 4.3 FIR濾波器的設計 16 4.3.1 Kaiser窗低通程式的設計 17 4.3.2 Kaiser窗高通程式的設計 19 4.3.3 Kaiser窗帶通程式的設計 20 4.3 回放語音信號 22 5 心得與體會 23 參考文獻 24 附錄 25  1 緒論 1.1課程設計背景 數字濾波器是數字信號處理的基礎，用來對信號進行過濾、檢測與參數估計等處理，在通信、圖像、語音、雷達等許多領域都有著十分廣泛的應用。尤其在圖像處理、數據壓縮等方面取得了令人矚目的進展和成就。由此可見，數字濾波器的設計尤為重要。 目前數字濾波器的設計有許多現成的高級語言設計程式，但他們都存在設計效率較低，不具有可視圖形，不便於修改參數等缺點，而Matlab為數字濾波的研究和應用提供了一個直觀、高效、便捷的利器。他以矩陣運算為基礎，把計算、可視化、程式設計融合到了一個互動式的工作環境中。尤其是Matlab工具箱使各個領域的研究人員可以直觀方便地進行科學研究與工程應用。其中的信號處理工具箱、圖像處理工具箱、小波工具箱等更是為數字濾波研究的蓬勃發展提供了可能。 本次課程設計主要介紹基於Matlab環境的IIR數字濾波器和FIR數字濾波器的設計，然後將設計結果應用於語音信號的處理，針對處理結果進一步分析了濾波器性能的。 1.2課程設計目的 (1)掌握數字信號處理的基本概念，基本理論和基本方法。 (2)熟悉離散信號和系統的時域特性。 (3)掌握序列快速傅裏葉變換方法。 (4)學會MATLAB的使用，掌握MATLAB的程式設計方法。 (5)掌握利用MATLAB對語音信號進行頻譜分析。 (6)掌握濾波器的網路結構。 (7)掌握MATLAB設計IIR、FIR數字濾波器的方法和對信號進行  2 課程設計預習 2.1卷積運算的演示 1. 線性卷積 序列x1(n)=[2 0 1 2 5 7 0 5 0 2 1 3，序列x2(n)=[ 2 0 1 1 9 9]。動態演示兩個序列進行線性卷積x1(n)﹡x2(n)的翻轉、移位、乘積、求和的過程。其中翻轉採用fliplr[1]，程式見附錄A。 圖2.1 線性卷積圖示  圖2.2 動態演示N=2和N=7 2.1.2 迴圈卷積 序列x1(n)=[ 2 0 1 2 5 7 0 5 0 2 1 3]，序列x2(n)=[ 2 0 1 1 9 9 ]，N=12。動態演示兩個序列進行迴圈卷積x1(n)⊙x2?(n)的翻轉、移位、乘積、求和的過程,程式見附錄B。 圖2.3 迴圈卷積圖示 當N12時，迴圈卷積發生混疊；當N≥20時，迴圈卷積等於線性卷積。 2.2 採樣定理的演示 信號進行採樣，A為幅度因數，a為衰減因數，為模擬角頻率，調試結果如圖 1所示，程式見附錄C 圖2.4 採樣定理圖示  3 課程設計的基本原理 3.1 IIR濾波器的設計原理 IIR數字濾波器的設計一般是利用目前已經很成熟的模擬濾波器的設計方法來進行設計，通常採用模擬濾波器原型有巴特沃斯濾波器、切比雪夫濾波器、橢圓濾波器等。 IIR數字濾波器單位取樣回應是時寬無限的h(n)，即n0＜n＜∞。 數字濾波器按照實現的方法和結構形式分為遞歸型或非遞歸型兩類。遞歸型數字濾波器的當前輸出y(n)是輸入x(n)的當前值和以前各輸入值x(n)，x(n–1)，….，及以前各輸出值y(n)，y(n–1)，….的函數。 一個N階遞歸型數字濾波器(IIR濾波器)的差分方程如公式2-1所示： 2-1 其中，式中的係數至少有一項不為零。說明必須將延時的輸出序列進行回饋。 遞歸系統的傳統函數定義如公式2-2所示： 2-2 遞歸系統的傳遞函數H(z)在Z平面上不僅有零點，而且有極點。非遞歸型數字濾波器當前的輸出值y(n)僅為當前及以前的輸入序列的函數，而與以前的各個輸出值無關，因此從結構上看非遞歸系統沒有回饋環路。 3.2 FIR濾波器的設計原理 FIR數字濾波器的單位取樣回應是時寬有限的h(