摘要：本文主要介紹基于DSP的數字濾波器的設計，使用CCS5000Simulator實現FTSK數據輸入,使用FIR濾波器對FTSK調制信號進行處理，輸出需要的波形與頻譜。文中采用線性緩沖區(qū)和帶移位雙操作尋址的方法實現FIR濾波器。


前言

一個實際的應用系統(tǒng)中，總存在各種干擾。使用DSP進行數字信號處理時，可以從噪聲中提取信號，即對一個具有噪聲和信號的混合源進行采樣，然后經過一個數字濾波器，濾除噪聲，提取有用信號；數字濾波器是DSP最基本的應用領域，也是熟悉DSP應用的重要環(huán)節(jié)。在系統(tǒng)設計中，濾波器的好壞將直接影響系統(tǒng)的性能。

數字濾波器的基本理論和設計

對于數字濾波器的系統(tǒng)函數可以表示為:


可以直接寫成表示輸出與輸入的關系,即常系數線性差分方程:以下是IIR濾波器的表達式


當全部,即系統(tǒng)函數和單位抽樣響應，則系統(tǒng)是FIR。

其濾波結構圖如圖一所示——橫向濾波結構

圖一橫向濾波器結構圖

FIR濾波算法實際上是一種乘法累加運算。它不斷輸入樣本，經延時，作乘法累加，再輸出濾波結果y（n）。在這里使用FIR濾波器，它有以下幾個特點：
（1）系統(tǒng)的單位沖激響應h（n）在有限個n值處不為零；
（2）系統(tǒng)函數H（z）在|z|>0處收斂，在|z|>0處只有零點，有限z平面只有零點，而全部極點都在z=0處；
（3）結構主要是非遞歸結構，沒有輸出到輸入的反饋。

這本次設計中FTSK輸入數據中包含頻率為800HZ，1200HZ，1600HZ，2021HZ，中心頻率為1600HZ，提取該頻率的信號。利用Matlab設計一個帶通濾波器。具體參數為：采樣頻率為22050HZ，通帶寬度為250HZ，則Fpass1=1475HZ，Fpass2=1725HZ，衰減1db，過渡帶為200HZ則Fstop1=1275HZ，Fstop2=1925HZ，阻帶衰減為30db。運行Matlab獲得126階的帶通濾波器，并提取系數。

用線性緩沖區(qū)和帶移位雙操作數尋址方法實現FIR濾波器

在這里介紹用線性緩沖區(qū)法實現，其特點是：
（1）對于N級的FIR濾波器，在數據存儲區(qū)中開辟一個稱之為滑窗的N個單元的緩沖區(qū)，存放最新的N個輸入樣本。
（2）從最老的樣本開始，每讀一個樣本后，將此樣本向下移位。讀完最后一個樣本后，輸入最新樣本至緩沖區(qū)的頂部。
在這里N=126，y（n）=
在數據存儲區(qū)中存放系數，并設置線性緩沖區(qū)存放輸入數據。

具體實現程序如下：
.title“simfir126.asm”
..global_firinit;兩個匯編程序，可作為C語言的調用，一個是初始化濾波器
.global_asmfir；另一個是濾波器執(zhí)行程序
.mmregs
X.usect“data1”，127
.bssdata1，1
SIZE.set127
.data
COEF；從matlab中獲得的濾波器系數
；濾波器初始化
_firinit：
SSBXFRCT
STM#X+（SIZE-1），AR2
STM#（SIZE-2）,AR0
RET
；濾波器執(zhí)行：
_asmfir：
STLA，*AR2+0
RPTZA，#（SIZE-2）
MACD*AR2-，COEF，A
STHA，*AR2
LD*AR2+，A
NOP
NOP
RETD
NOP
NOP
.END

Simulator仿真結果分析與硬件調試

仿真在系統(tǒng)調試中起著重要作用，TI公司也提供了軟件仿真器（Simulator）來調試程序。其中提供的探測點（ProbePoint）功能非常強大，它是一個開發(fā)算法的工具，將計算機文件數據傳送到目標板的buffer提供DSP軟件應用，同時可以將計算結果輸出到計算機文件中供分析，也可以通過CCS提供的圖形窗口觀察輸入輸出數據情況。

在本次設計中利用CCS提供的斷點和探測點，指定FTSK數據文件的輸入點，進行相關設定，同時利用CCS提供的圖形窗口觀察輸入和輸出的波形與頻譜。運行程序，分別得到輸入波形和頻譜圖（圖二），輸出波形和頻譜圖（圖三）下面分別對這兩個圖形進行分析。

輸入波形和頻譜圖


圖二輸入信號的波形和頻譜圖

由圖二輸入信號的波形圖可以看出，輸入是