《電子技術應用》
您所在的位置:首頁 > 可編程邏輯 > 業界動態 > FIR數字濾波器分布式算法的原理及FPGA實現

FIR數字濾波器分布式算法的原理及FPGA實現

2008-08-22
作者:畢占坤 吳伶錫

??? 摘 要: 在利用FPGA實現數字信號處理" title="數字信號處理">數字信號處理方面,分布式算法" title="分布式算法">分布式算法發揮著關鍵作用,與傳統的乘積-和結構相比,具有并行處理的高效性特點。詳細研究了基于FPGA、采用分布式算法實現FIR數字濾波器" title="數字濾波器">數字濾波器的原理和方法,并通過Xilinx ISE在Modelsim下進行了仿真。
??? 關鍵詞: 分布式算法? DALUT? FPGA??? FIR

?

??? 數字濾波器正在迅速地代替傳統的由R、L、C元件和運算放大器組成的模擬濾波器并且日益成為DSP的一種主要處理環節。FPGA也在逐漸取代ASIC和PDSP,用作前端數字信號處理的運算(如:FIR濾波、CORDIC算法或FFT)。乘累加運算是實現大多數DSP算法的重要途徑,而分布式算法則能夠大大提高乘累加運算的效能。
1 傳統的乘累加結構FIR數字濾波器基本理論
??? FIR濾波器被稱為有限長脈沖響應" title="脈沖響應">脈沖響應濾波器,與IIR數字濾波器相對應,它的單位脈沖響應h(n)只有有限個數據點。輸入信號經過線性時不變系統輸出的過程是一個輸入信號與單位脈沖響應進行線性卷積的過程,即:
???

式中,x(n)是輸入信號,y(n)是卷積輸出,h(n)是系統的單位脈沖響應。可以看出,每次采樣y(n)需要進行L次乘法和L-1次加法操作實現乘累加之和,其中L是濾波器單位脈沖響應h(n)的長度。可以發現,當L很大時,每計算一個點,則需要很長的延遲時間。
2 乘累加運算的位寬分配
??? DSP算法最主要的就是進行乘累加運算。假設采樣信號的位寬用N來表示,則N位與N位的乘積結果需要2N位的寄存器來保存;如果兩個操作數都是有符號數,則乘積只有2N-1個有效位,因為產生了兩個符號位。
??? 為了使累加器的結果不產生溢出,需要對累加器進行冗余設計,也就是說要在累加器2N的位寬上多設計出K位,累加器的長度M計算方式如下(L為濾波器的長度):
??? 對于無符號數:M=2N+K=2N+log2L
????對于有符號數:M=2N+K=2N+log2L-1
3 乘累加運算的分布式算法原理分析
??? 得益于Xilinx FPGA查找表" title="查找表">查找表結構的潛能,分布式算法在濾波器設計方面顯示出了很高的效率,自20世紀90年代初以來越來越受到人們的重視。分布式算法是基于查找表的一種計算方法,在利用FPGA實現數字信號處理方面發揮著重要的作用,可以大大提高信號的處理效率。它主要應用于數字濾波、頻率轉換等數字信號處理的乘累加運算。
??? 分布式算法推導如下:
??? 設Ak是已知常數(如濾波器系數、FFT中的正弦/余弦基本函數等),xk(n)是變量,可以看作是n時刻的第k個采樣輸入數據,y(n)代表n時刻的系統響應。那么它們的內積為:
???

其中,xk(n)變量可以寫成下面的格式:
???

式中,B為數據格式的字長;xkb是變量的二進制位,只有“0”和“1”兩種狀態。將(2)式代入(1)式得:
???

4 FPGA實現過程中查找表的構造方法
??? 根據以上論述,括號中的每一乘積項代表著輸入變量的某一位與常量的二進制“與”操作,加號代表著算術和操作,指數因子對括號中的值加權。如果事先構造一個查找表,該表存儲著括號中所有可能的組合值,就可以通過所有輸入變量相對應位的組合向量(xNb,x(N-1)b,…x1b)對該表進行尋址,該查找表稱為DALUT。DALUT的構造規則如表1所示。

?


??? 根據以上分析,可以得出FPGA實現分布式算法的硬件結構,如圖1所示。

?


5 采用分布式算法實現FIR數字濾波器
?? ?為了說明問題,以一個三個系數的FIR數字濾波器為例設計分布式算法,字寬也設置為三位。設FIR數字濾波器系數為:h(0)=5,h(1)=2,h(2)=3。相對應的ADLUT表如表2所示。

?


??? 在進行FPGA設計時,該表以組件Component形式構建,設置為ROM結構,提供輸入尋址端口table_in[2..0],輸出端口table_out[3..0]。FPGA算法的結構圖如圖2所示。

?


??? 算法實現中的幾個關鍵問題為:
??? (1)采用狀態機實現分布式算法的狀態轉移
??? 狀態機的實現如圖3所示,設置三個狀態s0、s1、s2。狀態s0完成數據的裝入,數據寄存器需要成對出現,一個完成數據的延遲,另一個完成數據的移位,并將狀態轉移到s1;狀態s1完成查找表功能、數據移位和分布式算法的乘累加運算,數據移位一個數據寬度后將狀態轉移到s2;狀態s2完成數據的輸出,并將狀態轉移到s0。利用狀態機可以條理清楚地簡化計算過程,在算法實現時發揮著關鍵的作用。
??? (2)系統時鐘與數據輸入時鐘的關系
??? 根據上述的狀態轉移關系,可以得出:每輸入一個數據,在下一次數據輸入之前,需要在狀態s1停留一個數據寬度(三位)的時鐘時間,在s2停留一個時鐘的數據輸出時間。也就是說,系統時鐘頻率應是數據輸入頻率的5倍,即fclkock=5fxin
??? (3)分布式算法中的乘累加公式推導及核心代碼實現
?? ?設B是數據的字寬,Pn是分布式算法第n位的結果,則有:
???

??? 有了該關系式,就可以通過for...loop循環,使用一條語句完成分布式乘累加算法。具體如下:
??? for n in 0 to B-1 loop
???? ??P:=P/2+tableout(n)*2B-1;
??? end loop;
6 算法仿真驗證與結論
?? ?本文實現的FIR濾波器在Xilinx的集成開發環境ISE下利用ModelSim進行了仿真。當輸入數據為7,3,1...時,仿真輸出依次為35,29,32,16...,與乘累加方式FIR濾波算法得出的結果完全一致。假設查找表和PDSP的通用乘法器延遲時間相同,分布式算法的等待時間是B?子,通用乘法器的等待時間是NI。可見,對于位寬較小的數據來說,分布式算法的執行速度遠高于乘累加運算。可見,利用FPGA實現分布式計算大大提高了計算的速度,在高速信號處理中發揮著重要作用。
參考文獻
1劉 偉,孟憲元.利用FPGA實現小型聲納的片上系統集成.電子技術應用,2000;26(10)
2 盧 毅. VHDL與數字電路設計.北京:科學出版社, 2002
3 褚振勇. FPGA設計與應用. 西安:西安電子科技大學出版社,2002
4 程佩青.數字信號處理教程.北京:清華大學出版社,2000

本站內容除特別聲明的原創文章之外,轉載內容只為傳遞更多信息,并不代表本網站贊同其觀點。轉載的所有的文章、圖片、音/視頻文件等資料的版權歸版權所有權人所有。本站采用的非本站原創文章及圖片等內容無法一一聯系確認版權者。如涉及作品內容、版權和其它問題,請及時通過電子郵件或電話通知我們,以便迅速采取適當措施,避免給雙方造成不必要的經濟損失。聯系電話:010-82306118;郵箱:aet@chinaaet.com。
主站蜘蛛池模板: 最近中文国语字幕在线播放视频 | 国产无套粉嫩白浆在线精品 | 亚洲综合男人的天堂色婷婷 | 欧美日韩视频一区三区二区 | 99久久伊人 | 成人自拍网站 | 九九99视频在线观看视频观看 | 另类视频欧美 | 巨尻在线播放 | 在线观看亚洲天堂 | 亚洲综合网址 | 精品久久久久久中文字幕202o | 蜜臀91精品国产免费观看 | 精品国产欧美一区二区三区成人 | 成年网站免费视频黄 | 午夜视频一区二区 | 国产成人综合高清在线观看 | 嘿嘿嘿视频免费网站在线观看 | 亚洲视频99| 亚洲精品欧美日本中文字幕 | 日韩免费视频观看 | 日本欧美一区二区三区免费不卡 | 青草91| 黄网站免费视频 | videoxxoo欧美老师 | 久久精品亚洲综合 | 国产va免费精品观看 | 欧洲视频一区 | 国产精品推荐天天看天天爽 | 911精品国产亚洲日本美国韩国 | 欧美一级片在线观看 | 成人免费激情视频 | 私人午夜影院 | 日韩一中文字幕 | 午夜视频免费观看黄 | 亚洲国产成a人v在线观看 | 国产高清大尺度一区二区不卡 | 一级午夜免费视频 | 免费看色网站 | 福利在线观看视频 | 国外成人在线视频网站 |