基于FDATool的FIR濾波器設(shè)計(jì)方法（二）

　　脈動(dòng)型(Systolic)FIR濾波器設(shè)計(jì)#e#

　　脈動(dòng)型FIR濾波器是對(duì)直接型的升級(jí)，在每個(gè)操作后都加入流水線級(jí)，每個(gè)動(dòng)作都打一拍，就跟心臟跳動(dòng)一樣，因此稱為脈動(dòng)型，這種結(jié)構(gòu)非常適用于高速數(shù)據(jù)流的處理。如圖1所示為脈動(dòng)型FIR濾波器結(jié)構(gòu)。

　　圖1

　　與直接型結(jié)構(gòu)不同的是，輸入數(shù)據(jù)到下一個(gè)處理單元都需要打2拍，這是為了使乘法后的累加數(shù)據(jù)同步，下面推導(dǎo)驗(yàn)證：

　　x(n)為輸入數(shù)據(jù)，yt(n)為直接型結(jié)構(gòu)的輸出

　　yt(n)=x(n)h(0)+x(n-1)h(1)+x(n-2)h(2)…x(n-10)h(10)

　　ys(n)為脈動(dòng)型結(jié)構(gòu)的輸出，如圖1中有P1、P2…P10共10個(gè)節(jié)點(diǎn)

　　P1=x(n-4)h(0)

　　P2=(P1 + x(n-5)h(1))*Z-1=x(n-5)h(0) + x(n-6)h(1)

　　…

　　P10=(P9 + x(n-23)h(10))*Z-1

　　ys(n)=x(n-14)h(0) + x(n-15)h(1) + … + x(n-23)h(9)+ x(n-24)h(10)

　　由ys(n)和yt(n)的表達(dá)式，可以推導(dǎo)出ys(n)=yt(n-14)

　　因此脈動(dòng)型FIR濾波器的延遲較大

　　如圖2所示為11抽頭系數(shù)脈動(dòng)型FIR濾波器FPGA實(shí)現(xiàn)結(jié)構(gòu)(實(shí)例與前幾節(jié)相同)，穿了一層“衣服”，采用Xilinx FPGA中的DSP48E1 實(shí)現(xiàn)，基本處理單元中的操作都可在一個(gè)DSP48E1中完成，輸入數(shù)據(jù)經(jīng)過(guò)DSP48E1中寄存2拍后通過(guò)ACOUT輸出，直接連接到下一個(gè) DSP48E1中的ACIN端口，累加輸出PCOUT直接連接到下一個(gè)DSP48E1中的PCIN端口，這些連接都沒(méi)有經(jīng)過(guò)FPGA的Fabric連線邏輯，而是通過(guò)DSP Block的內(nèi)部走線連接，這樣實(shí)現(xiàn)能夠縮短路徑的延時(shí)。

　　圖2

　　編寫了相關(guān)代碼，綜合結(jié)果如下：

　　Number of Slice Registers： 4

　　Number of Slice LUTs： 19

　　Number of DSP48E1s： 11

　　Minimum period： 3.006ns{1} (Maximum frequency： 332.668MHz)

　　線性相位實(shí)現(xiàn)：

　　與前幾節(jié)相同，由于FIR濾波器的線性相位特性，相對(duì)應(yīng)有線性相位的實(shí)現(xiàn)結(jié)構(gòu)，如圖3所示，利用DSP48E1中預(yù)加器實(shí)現(xiàn)乘法前的加法操作。對(duì)于脈動(dòng)型 FIR濾波器的線性相位結(jié)構(gòu)有很多注意點(diǎn)，其中預(yù)加器數(shù)據(jù)的配對(duì)，常規(guī)情況下，此例中應(yīng)是x(n)和x(n-10)、x(n-1)和x(n-9)、 x(n-2)和x(n-8)、x(n-3)和x(n-7)、x(n-4)和x(n-6)，而圖3中結(jié)構(gòu)，加入了延時(shí)11的移位寄存器，預(yù)加器配對(duì)的數(shù)據(jù)為 x(n-2)和n(n-12)、x(n-4)和x(n-12)、x(n-6)和x(n-12)、x(n-8)和x(n-12)、x(n-10)和x(n- 12)，可以發(fā)現(xiàn)預(yù)加器配對(duì)數(shù)據(jù)中有一個(gè)數(shù)據(jù)始終是x(n-12)，但是每一個(gè)配對(duì)數(shù)據(jù)的相對(duì)延時(shí)與常規(guī)情況下相同：10、8、6、4和2。

　　圖3

　　而各節(jié)點(diǎn)P1、P2、P3、P4、P5和y(n)的表達(dá)式如下

　　P1=x(n-5)h(0) + x(n-15)h(0)

　　P2=( P1 + (x(n-6)h(1) + x(n-14)h(1)) )Z-1=x(n-6)h(0) + x(n-16)h(0) + x(n-7)h(1) + x(n-15)h(1)

　　P3=( P2 + (x(n-8)h(2) + x(n-14)h(2)) )Z-1=x(n-7)h(0) + x(n-17)h(0) + x(n-8)h(1) + x(n-16)h(1) + x(n-9)h(2) + x(n-15)h(2)

　　P4=( P3 + (x(n-10)h(3) + x(n-14)h(3)) )Z-1=x(n-8)h(0) + x(n-18)h(0) + x(n-9)h(1) + x(n-17)h(1) + x(n-10)h(2) + x(n-16)h(2) + x(n-11)h(3) + x(n-15)h(3)

　　P5=( P4 + (x(n-12)h(4) + x(n-14)h(4)) )Z-1=x(n-9)h(0) + x(n-19)h(0) + x(n-10)h(1) + x(n-18)h(1) + x(n-11)h(2) + x(n-17)h(2) + x(n-12)h(3) + x(n-16)h(3) + x(n-13)h(4) + x(n-15)h(4)

　　y(n)=(P5 + x(n-14)h(5))Z-1= x(n-10)h(0) + x(n-20)h(0) + x(n-11)h(1) + x(n-19)h(1) + x(n-12)h(2) + x(n-18)h(2) + x(n-13)h(3) + x(n-17)h(3) + x(n-14)h(4) + x(n-16)h(4) + x(n-15)h(5)

　　因抽頭系數(shù)對(duì)稱，由h(0)=h(10)，h(1)=h(9)，h(2)=h(8)，h(3)=h(7)，h(4)=h(6)可得

　　y(n)= x(n-10)h(0) + x(n-11)h(1) + x(n-12)h(2) + x(n-13)h(3) + x(n-14)h(4) + x(n-15)h(5) + x(n-16)h(6) + x(n-17)h(7) + x(n-18)h(8) + x(n-19)h(9) + x(n-20)h(10)

　　驗(yàn)證得到y(tǒng)(n)=yt(n-10)，比普通脈動(dòng)結(jié)構(gòu)延時(shí)小，但是相比于其他結(jié)構(gòu)的FIR濾波器延時(shí)還是較大的。

　　編寫了相關(guān)代碼，綜合結(jié)果如下：

　　Number of Slice Registers： 84

　　Number of Slice LUTs： 99

　　Number of DSP48E1s： 6

　　Minimum period： 3.256ns{1} (Maximum frequency： 307.125MHz)

　　在DSP in FPGA： FIR濾波器設(shè)計(jì)(一)、(二)中分別講解了直接型、轉(zhuǎn)置型和脈動(dòng)型結(jié)構(gòu)FIR濾波器的實(shí)現(xiàn)方法，這三種結(jié)構(gòu)是FPGA實(shí)現(xiàn)中比較常用的方法，以下對(duì)這三種結(jié)構(gòu)做一個(gè)比較：

　　(1) 直接型：方法簡(jiǎn)單易實(shí)現(xiàn)，但是使用加法樹(shù)優(yōu)化后增加了功耗

　　(2) 轉(zhuǎn)置型：關(guān)鍵路徑延時(shí)較小，時(shí)序易滿足，但是輸入數(shù)據(jù)扇出較大，不適用于階數(shù)較高的濾波器實(shí)現(xiàn)

　　(3) 脈動(dòng)型：適用于高速數(shù)據(jù)處理，但是延時(shí)相比于其它結(jié)構(gòu)較大