基于脈動陣列的FIR濾波器設(shè)計

　　1引 言

　　有限長沖激響應(yīng)(FIR)濾波器在數(shù)字信號處理中是一種基本的處理單元。無限長單位沖激響應(yīng)(IIR)數(shù)字濾波器的優(yōu)點是可以利用模擬濾波器設(shè)計的結(jié)果，但其缺點是不具有線性相位性。圖像處理以及數(shù)據(jù)傳輸都要求信道具有線性相位特性，FIR濾波器可以做成嚴(yán)格的線性相位，避免被處理信號產(chǎn)生相位失真，還可以具有任意的幅度特性。此外，FIR濾波器的單位沖激響應(yīng)是有限長的，因而濾波器一定是穩(wěn)定的。

　　在數(shù)字濾波器的研究中，已經(jīng)提出多種FIR濾波器的設(shè)計和實現(xiàn)方法，如并行結(jié)構(gòu)、流水線結(jié)構(gòu)、分布式結(jié)構(gòu)等[1-3]。FIR濾波器計算量大，且要求實時實現(xiàn)。如何提高速度以滿足信號處理的高效性和實時性一直是人們研究的重點和熱點。脈動陣列是一種新型的流水線結(jié)構(gòu)，所有處理單元是相同且全流水的，并且有模塊化和規(guī)則化的特征，這對于硬件設(shè)計是很重要的。脈動陣列結(jié)構(gòu)可以很好地滿足高速實時信號處理的需要，提高數(shù)據(jù)的通過率及電路的執(zhí)行速度。該結(jié)構(gòu)結(jié)合了并行結(jié)構(gòu)和流水線的優(yōu)勢，使FIR濾波器達(dá)到了更高的處理速度。

　　數(shù)字濾波器現(xiàn)在大多在DSP芯片上實現(xiàn)[4]。隨著集成電路技術(shù)的發(fā)展，現(xiàn)在已有多種專用DSP芯片用于實現(xiàn)濾波、FFT等運算。另外，F(xiàn)PGA內(nèi)部的規(guī)整的邏輯塊陣列，很適合實現(xiàn)脈動陣列這種高度并行的運算。Altera公司新推出的Stratix系列FPGA內(nèi)部有豐富的邏輯陣列資源及高性能嵌入式DSP塊，能夠高效地實現(xiàn)數(shù)字信號處理功能[5]。

　　本文首先介紹了FIR濾波器和脈動陣列的原理，然后設(shè)計了脈動陣列結(jié)構(gòu)的FIR濾波器，畫出電路的結(jié)構(gòu)框圖，并進(jìn)行了時序分析，最后在FPGA上進(jìn)行驗證。結(jié)果表明，脈動陣列的模塊化和高度流水線的結(jié)構(gòu)使FIR濾波器在FPGA上獲得了很好的性能，比串行結(jié)構(gòu)的運算速度更快，呵以更好地滿足數(shù)字信號處理中高效、實時的要求。

　　2 FIR濾波器及脈動陣列原理

　　2.1 FIR濾波器原理

　　數(shù)字濾波器用于改變輸入信號X(n)的頻譜特性以滿足某種特定的設(shè)計要求。一個因果的數(shù)字濾波器可以用它的單位沖激響應(yīng)h(n)、傳輸函數(shù)H(z)或者差分方程來表達(dá)。其中單位沖激響應(yīng)和傳輸函數(shù)描述了系統(tǒng)的時域和頻域性質(zhì)，差分方程則反映了實現(xiàn)該濾波器所需的運算。

　　一個線性時不變因果濾波器可以用式(1)的差分方程描述：

　　如果1≤k≤N時，ak=0，則：

　　這就是一個M階的FIR濾波器，是非遞歸運算。

　　2.2 脈動陣列

　　脈動陣列結(jié)構(gòu)是1978年由Carneige-Mellon大學(xué)的H.T.Kung等人發(fā)展起來的一種專用處理器設(shè)計模型，是一種新的流水線結(jié)構(gòu)，表示一種有節(jié)奏地汁算并通過系統(tǒng)傳輸數(shù)據(jù)的處理單元網(wǎng)絡(luò)。這些處理單元規(guī)則地泵入泵出數(shù)據(jù)以維持規(guī)則的數(shù)據(jù)流[6]。由于其數(shù)據(jù)運算與傳送方式類似于心臟或脈搏有節(jié)奏地跳動，因此被稱為心動陣列或脈動陣列。利用脈動陣列可設(shè)計出一系列適合數(shù)字信號處理應(yīng)用的模塊化、規(guī)則且有效的運算結(jié)構(gòu)。脈動陣列結(jié)構(gòu)的基本特性如下：

　　(1)脈動性

　　數(shù)據(jù)在一個全局時鐘的精確控制下，有節(jié)奏地經(jīng)過計算和傳遞過程，最后通過網(wǎng)絡(luò)。

　　(2)模塊性和規(guī)則性

　　陣列由模塊化的處理單元組成，各處理單元之間的互連方式均勻一致，并且計算網(wǎng)絡(luò)可任意擴(kuò)大。

　　(3)空間局部性和時間局部性

　　表現(xiàn)在數(shù)據(jù)從一個計算節(jié)點傳遞到下一個計算節(jié)點至少需要一個單元時間。

　　(4)有效的流水線結(jié)構(gòu)

　　脈動陣列免除了形成數(shù)據(jù)流所需的控制開銷。陣列內(nèi)處理單元間的局部連接方式使陣列中負(fù)載均勻、連線極短，最大限度地減小了系統(tǒng)內(nèi)部的通信延時，提高了處理單元的利用率，使整個陣列的系統(tǒng)性能得到充分發(fā)揮。

　　脈動結(jié)構(gòu)是將線性映射技術(shù)用于規(guī)則依賴圖上進(jìn)行設(shè)計的。依賴圖的邊表示優(yōu)先約束。脈動陣列設(shè)計中的基本矢量有：

　　(1)投影矢量：

　　(2)處理器矢量：PT=(p1 p2)

　　(3)調(diào)度矢量：ST=(S1 S2)

　　(4)硬件利用率：HUE=1/│STd │

　　3 FlR濾波器的脈動結(jié)構(gòu)設(shè)計

　　3.1 FIR濾波器的脈動陣列結(jié)構(gòu)

　　用調(diào)度不等式選擇可行的調(diào)度矢量ST，并根據(jù)脈動結(jié)構(gòu)的可行性限制條件選擇投影矢量d和處理器空間矢量PT，然后用線性映射技術(shù)設(shè)計脈動列。

　　選擇投影矢量、處理器矢量和調(diào)度矢量如下：

　　于是，節(jié)點IT=(i，j)被處理為：

　　節(jié)點IT=(i，j)的執(zhí)行時間為：

　　硬件利用率：

　　以五階FIR濾波器為例：

　　FIR濾波器的脈動結(jié)構(gòu)設(shè)計框圖如圖1所示。

　　其中，d0,d1，d2和d3作為延時單元，m0，m1，…，m4為固定乘數(shù)的乘法器，乘數(shù)分別為ω0，ω1，…，ω4。乘法器m0和加法器a0構(gòu)成處理器P0;乘法器m1和加法器a1構(gòu)成處理器P1;類似的，m4和a4構(gòu)成處理器P4。5個處理器同時工作，x(n)輸入后同時到達(dá)5個處理器，分別和權(quán)重ωi相乘后，經(jīng)不同的延時相加得到結(jié)果y(n)。

　　該設(shè)計框圖很好地體現(xiàn)了2.2節(jié)提到的脈動陣列的4個基本特性。數(shù)據(jù)在網(wǎng)絡(luò)中有節(jié)奏的輸入、計算、傳遞并輸出。

　　該結(jié)構(gòu)由5個相同的模塊化處理器構(gòu)成，容易擴(kuò)展，可實現(xiàn)任意階的FIR濾波器。數(shù)據(jù)的輸入、和權(quán)重的相乘、中間結(jié)果的寄存、數(shù)據(jù)的輸出各需要一個單元時問，但這幾種操作可同時進(jìn)行，不會互相干擾。該脈動陣列結(jié)構(gòu)以高度的流水線方式運行。

　　需要注意的是，在乘法器輸出的時候需要對輸出的數(shù)據(jù)進(jìn)行一位擴(kuò)展，以避免加法器的溢出問題。

　　3.2 時序分析

　　該結(jié)構(gòu)的空問時間表示如圖2所示。橫軸為時間軸，縱軸為處理器軸。圖中的每一行對應(yīng)一個處理器?？梢钥吹皆跁r間軸的某一點上，輸入數(shù)據(jù)x(n)在同一時間劍達(dá)所有處理器，即輸入"廣播";在處理器軸的某一點上，權(quán)重ωi在各處理器中處于相同的空間坐標(biāo)，因此權(quán)重"保持";而輸出數(shù)據(jù)y(n)在不同的空間和時間得到，所以輸出"移動"。該FIR濾波器為輸入廣播、權(quán)重保持、輸出移動型的脈動陣列結(jié)構(gòu)。

　　借助于脈動陣列技術(shù)，輸入和輸出之間的處理可以同時進(jìn)行，因此，總的執(zhí)行時間變得最小。對該結(jié)構(gòu)進(jìn)行時序分析，如圖3所示。在第一個時鐘周期x(0)送到各處理器，和權(quán)重分別相乘后在第二個時鐘周期得到y(tǒng)(0)=ω0x(0)，在此同時數(shù)據(jù)x(1)已經(jīng)輸入;第三個時鐘周期輸入x(2)，在處理器P0得到ω0x(1)，此時處理器P1的結(jié)果ω1x(0)經(jīng)一級延時后和P0的結(jié)果相加得到y(tǒng)(1)=ω0(1)+ω1x(O)。同理，第四個周期得到y(tǒng)(2)=ω0x(2)+ω1x(1)+ω2x(0)，第五個周期得到y(tǒng)(3)=ω0x(3)+ω1x(2)+ω2x(1)+ω3x(O)…。

　　4 FPGA驗證

　　以一個輸入輸出均為8 b的五階FIR濾波器為例，對該脈動陣列結(jié)構(gòu)的FIR濾波器在FPGA上進(jìn)行驗證。目標(biāo)器件選用Stratix系列器件，編譯軟件為QuartusⅡ4.2。

　　Stratix系列器件是Altera公司推出的新型FPGA，其內(nèi)部有豐富的嵌入式存儲器、優(yōu)化的數(shù)字信號處理(DSP)塊和高性能I/O能力。Stratix器件是設(shè)計復(fù)雜的高性能系統(tǒng)的理想選擇。StratiX器件使用DSP塊實現(xiàn)大計算量應(yīng)用所需的大數(shù)據(jù)吞吐量。DSP塊由硬件乘法器、加法器、減法器、累加器和流水線寄存器組成[6]。從圖1中可以看出，每個處理器單元由一個乘法器和一個加法器構(gòu)成，這種結(jié)構(gòu)很適合用FPGA內(nèi)部的DSP塊實現(xiàn)。

　　實驗結(jié)果如表1所示。為了便于比較，使用相同的器件，還給出了串行結(jié)構(gòu)FIR濾波器的實驗結(jié)果。在串行結(jié)構(gòu)中，每次計算1個x(n)，需要3個時鐘周期。脈動陣列結(jié)構(gòu)每次計算5個x(n)，需要5個時鐘周期。每個x(n)為8 b數(shù)據(jù)。運算速度可以由下面的公式得出：

　　運算速度=(每次處理數(shù)據(jù)長度/運算耗用的時鐘周期數(shù))