一、實(shí)驗(yàn)?zāi)康?/h1>

了解IIR濾波器的特點(diǎn)，掌握Matlab輔助設(shè)計濾波器系數(shù)的方法，并實(shí)現(xiàn)IIR濾波器濾除高頻信號，并在LCD上顯示結(jié)果。

二、實(shí)驗(yàn)原理

IIR濾波器

IIR無限沖激響應(yīng)數(shù)字濾波器（infinite impulse response digital filter）是對單位沖激的輸入信號的響應(yīng)為無限長序列的數(shù)字濾波器?？煞譃橐痪S、二維或多維無限沖激響應(yīng)數(shù)字濾波器。

它的輸出y（n）由當(dāng)前的和過去的輸入信號x(n)及過去的輸出信號共同決定。IIR的幅頻特性精度很高，不是線性相位的，可以應(yīng)用于對相位信息不敏感的音頻信號上。

對于IIR濾波器，沖激響應(yīng)理論上應(yīng)會無限持續(xù)，其輸出不僅取決于當(dāng)前和過去的輸入信號值，也取決于過去的信號輸出值。其系統(tǒng)函數(shù)如下：

IIR濾波器的一種結(jié)構(gòu)和差分方程可見圖示，由IIR的系統(tǒng)函數(shù)出發(fā)，可視作兩個系統(tǒng)的級聯(lián)，并且合并了公共的延時支路。

IIR濾波器特性

無限沖激響應(yīng)數(shù)字濾波器具有以下特性：

（1）系統(tǒng)函數(shù)可以寫成封閉函數(shù)的形式。

（2）采用遞歸型結(jié)構(gòu)，即結(jié)構(gòu)上帶有反饋環(huán)路。IIR濾波器運(yùn)算結(jié)構(gòu)通常由延時、乘以系數(shù)和相加等基本運(yùn)算組成，可以組合成直接型、正準(zhǔn)型、級聯(lián)型、并聯(lián)型四種結(jié)構(gòu)形式，都具有反饋回路。由于運(yùn)算中的舍入處理，使誤差不斷累積，有時會產(chǎn)生微弱的寄生振蕩。

（3）在設(shè)計上可以借助成熟的模擬濾波器的成果，如巴特沃斯、契比雪夫和橢圓濾波器等，有現(xiàn)成的設(shè)計數(shù)據(jù)或圖表可查，其設(shè)計工作量比較小，對計算工具的要求不高。

（4）相位特性不好控制，對相位要求較高時，需加相位校準(zhǔn)網(wǎng)絡(luò)。

IIR濾波器的設(shè)計

設(shè)計IIR數(shù)字濾波器實(shí)質(zhì)上就是找到一個使其頻率響應(yīng)H(Z)滿足給定的通帶截止頻率、通帶衰減系數(shù)、阻帶截止頻率及阻帶衰減系數(shù)的可以在物理上實(shí)現(xiàn)的系統(tǒng)函數(shù)H(Z)。

根據(jù)所要設(shè)計濾波器的參數(shù)去確定一個模擬濾波器的傳輸函數(shù)，然后再根據(jù)這個傳輸函數(shù)，通過雙線性變換、或脈沖響應(yīng)不變法來進(jìn)行數(shù)字濾波器的設(shè)計。它的設(shè)計比較復(fù)雜，復(fù)雜在于它的模擬濾波器傳輸函數(shù)H（s）的確定。這一點(diǎn)可以讓軟件來實(shí)現(xiàn)。具體實(shí)現(xiàn)步驟：

（1）先確定需要一個什么樣的濾波器，巴特沃斯型，切比雪夫型，還是其它什么型的濾波器。

（2）當(dāng)選定一個型號后，就可以根據(jù)設(shè)計參數(shù)和這個濾波器的計算公式來確定其階數(shù)、傳輸函數(shù)的表達(dá)式。通常這個過程中還存在預(yù)扭曲的問題（這只是雙線性變換法所需要注意的問題，脈沖響應(yīng)不變法不存在這種問題）。

（3）確定H（S）后，就可以通過雙線性變換得到其數(shù)字域的差分方程。

MATLAB FIR 濾波器輔助設(shè)計，生成濾波器系數(shù)。

打開Matlab 軟件，并打開"Filter Designed&Analaysis Tool"工具，在彈出的界面中按照所需濾波器修改以下參數(shù)：

（1）濾波器類型：低通Lowpass；

（2）濾波方式：巴特沃斯IIR（Butterworth）；

（3）濾波器階數(shù)：Specify order（4+1階）；

（4）采樣頻率：1000Hz；

（5）通帶截止頻率：150Hz；

程序流程設(shè)計

設(shè)計中首先要進(jìn)行外設(shè)使能配置，然后進(jìn)行LCD管腳復(fù)用配置和LCD中斷配置，接著進(jìn)行LCD顯示和觸摸的初始化。最后產(chǎn)生帶有噪聲的原始信號并進(jìn)行IIR濾波和FFT運(yùn)算，將結(jié)果顯示在LCD屏幕上，在循環(huán)中進(jìn)行觸摸檢測，根據(jù)標(biāo)志位判斷顯示時域或者頻域波形。

程序功能：使用MATLAB 輔助設(shè)計濾波器系數(shù)實(shí)現(xiàn)IIR濾波，并在LCD上顯示結(jié)果。

數(shù)字信號處理庫

DSPLIB 包含優(yōu)化的、C語言可調(diào)用的通用信號處理例程，用于計算密集型實(shí)時應(yīng)用程序。 調(diào)用這些例程的運(yùn)行速度比直接用C語言編寫的等效代碼快得多。使用DSPLIB可以縮短應(yīng)用程序開發(fā)時間。

DSPLIB 3_4_0_0包括適用于 C64x+ 或 C66x 或 C674x 處理器的 Windows 或 Linux 安裝可執(zhí)行文件。 每個可執(zhí)行文件安裝一個組件包存儲庫、一個文檔目錄、一個 Eclipse 插件目錄和一個擴(kuò)展的組件目錄結(jié)構(gòu)，其中包含組件庫、頭文件和測試示例。

TMS320C6748處理器使用的是dsplib_c674x_3_4_0_0。

函數(shù)源碼

程序使用DSPLIB 的庫來進(jìn)行IIR濾波，調(diào)用的程序源碼和使用說明可以安裝DSPLIB后查看。調(diào)用的IIR濾波函數(shù)中，

第一個參數(shù)是延遲元素值。

第二個參數(shù)是指向輸入數(shù)組的指針。

第三個參數(shù)是指向輸出數(shù)組的指針。

第四個參數(shù)是自回歸濾波器系數(shù)。

第五個參數(shù)是動態(tài)平均濾波器系數(shù)。

第六個參數(shù)是輸出樣本數(shù)。

程序使用DSPLIB 的庫來進(jìn)行FFT運(yùn)算，調(diào)用的程序源碼和使用說明可以安裝DSPLIB后查看。調(diào)用的FFT函數(shù)中，

第一個參數(shù)是樣本中FFT 的長度，

第二個參數(shù)是指向數(shù)據(jù)輸入的指針。

第三個參數(shù)是指向復(fù)雜旋轉(zhuǎn)因子的指針。

第四個參數(shù)是指向復(fù)雜輸出數(shù)據(jù)的指針。

第五個參數(shù)是指向包含64 個條目的位反轉(zhuǎn)表的指針。如果樣本的FFT長度可以表示為 4 的冪，

第六個參數(shù)是4，否則第六個參數(shù)是2 。

第七個參數(shù)是從主FFT開始的樣本中的子 FFT偏移索引 。

第八個參數(shù)是樣本中主FFT的大小。

程序使用DSPLIB 的庫來進(jìn)行FFT逆變換，調(diào)用的程序源碼和使用說明可以安裝DSPLIB后查看。調(diào)用的IFFT函數(shù)中，

第一個參數(shù)是樣本中FFT 的長度。

第二個參數(shù)是指向數(shù)據(jù)輸入的指針。

第三個參數(shù)是指向復(fù)雜旋轉(zhuǎn)因子的指針。

第四個參數(shù)是指向復(fù)雜輸出數(shù)據(jù)的指針。

第五個參數(shù)是指向包含64 個條目的位反轉(zhuǎn)表的指針 。如果樣本的FFT長度可以表示為 4 的冪，

第六個參數(shù)是4，否則第六個參數(shù)是 2 。

第七個參數(shù)是從主FFT開始的復(fù)雜樣本中的子FFT偏移索引 。

第八個參數(shù)是樣本中主FFT的大小。

三、操作現(xiàn)象

硬件連接

（1）連接仿真器和電腦的USB接口。

（2）將撥碼開關(guān)撥到DEBUG模式01111，連接實(shí)驗(yàn)箱電源，撥動電源開關(guān)上電。

操作現(xiàn)象

導(dǎo)入工程，選擇Demo文件夾下的對應(yīng)工程

編譯工程，生成可執(zhí)行文件

將CCS連接實(shí)驗(yàn)箱并加載程序

程序加載完成后點(diǎn)擊運(yùn)行程序

運(yùn)行程序后，LCD會顯示IIR濾波前后的時域波形，上方為100Hz信號與450Hz噪聲組成的混合信號，下方為濾波后得到的結(jié)果（保留100Hz信號）。

點(diǎn)擊一下LCD屏幕，會切換頻域波形，上方顯示原始頻域信號+噪聲，下方顯示濾波后的頻域波形，輸入波形為一個低頻率的正弦波與一個高頻的正弦波疊加而成，通過觀察頻域圖可知輸入波形中的低頻波形通過了濾波器，而高頻部分被濾除。

實(shí)驗(yàn)結(jié)束后，先點(diǎn)擊黃色按鈕暫停程序運(yùn)行，再點(diǎn)擊紅色按鈕退出CCS與實(shí)驗(yàn)板的連接，最后實(shí)驗(yàn)箱斷電即可。