引言
在信息信號處理過(guò)程中,如對信號的過(guò)濾、檢測、預測等,都要使用濾波器,數字濾波器是數字信號處理(DSP,DigitalSignalProcessing)中使用最廣泛的一種器件。常用的濾波器有無(wú)限長(cháng)單位脈沖響應(ⅡR)濾波器和有限長(cháng)單位脈沖響應(FIR)濾波器兩種[1],其中,FIR濾波器能提供理想的線(xiàn)性相位響應,在整個(gè)頻帶上獲得常數群時(shí)延從而得到零失真輸出信號,同時(shí)它可以采用十分簡(jiǎn)單的算法實(shí)現,這兩個(gè)優(yōu)點(diǎn)使FIR濾波器成為明智的設計工程師的首選,在采用VHDL或VerilogHDL等硬件描述語(yǔ)言設計數字濾波器時(shí),由于程序的編寫(xiě)往往不能達到良好優(yōu)化而使濾波器性能表現一般。而采用調試好的IPCore需要向Altera公司購買(mǎi)。筆者采用了一種基于DSPBuilder的FPGA設計方法,使FIR濾波器設計較為簡(jiǎn)單易行,并能滿(mǎn)足設計要求。
2 FIR濾波器介紹
2.1 FIR濾波器設計的原理
FIR濾波器的數學(xué)表達式可用差分方程(1)來(lái)表示:
其中:r是FIR的濾波器的抽頭數;b(r)是第r級抽頭數(單位脈沖響應);x(n-r)是延時(shí)r個(gè)抽頭的輸入信號。
設計濾波器的任務(wù)就是尋求一個(gè)因果,物理上可實(shí)現的系統函數H(z),使其頻率響應H(ejw)滿(mǎn)足所希望得到的頻域指標。
2.2 設計要求
數字濾波器實(shí)際上是一個(gè)采用有限精度算法實(shí)現的線(xiàn)性非時(shí)變離散系統,它的設計步驟為先根據需要確定其性能指標,設計一個(gè)系統函數H(z)逼近所需要的技術(shù)指標,最后采用有限的精度算法實(shí)現。本系統的設計指標為;設計一個(gè)16階的低通濾波器,對模擬信號的采樣頻率fs為48KHz要求信號的截止頻率fc=10.8kHz輸入序列位寬為9位(最寬位為符號位)。
3 DSPBuilder介紹
DSPbuilder是Altera推出的一個(gè)DSP開(kāi)發(fā)工具,它在QuartusⅡFPGA設計環(huán)境中集成了Mathworks的Matlab和simulinkDSP開(kāi)發(fā)軟件[2]。
以往Matlab工具的使用往往作為DSP算法的建模和基于純數學(xué)的仿真,其數學(xué)模型無(wú)法為硬件DSP應用系統直接產(chǎn)生實(shí)用程序代碼,仿真測試的結果也僅僅是基于數學(xué)算法結構。而以往FPGA所需的傳統的基于硬件描述語(yǔ)言(HDL)的設計由于要考慮FPGA的硬件的δ延時(shí)與VHDL的遞歸算法的銜接,以及補碼運算和乘積結果截取等問(wèn)題,相當繁雜。
對DSP是Builder而言,頂層的開(kāi)發(fā)工具是MatLab/Simulink整個(gè)開(kāi)發(fā)流層幾乎可以在同一環(huán)境中完成,真正實(shí)現了自定向下的設計流程,包括DSP系統的建模、系統級仿真、設計模型向VHDL硬件描述語(yǔ)言代碼的轉換、RTL(邏輯綜合RegisterTransferLevel)級功能仿真測試、編譯適配和布局布線(xiàn)、時(shí)序實(shí)時(shí)仿真直至對DSP目標器件的編程配置,整個(gè)設計流程一氣呵成地將系統描述和硬件實(shí)現有機地融為一體,充分顯示了現代電子設計自動(dòng)化開(kāi)發(fā)的特點(diǎn)與優(yōu)勢。
4 FIR數字濾波器的DSPBuilder設計
4.1 FIR濾波器參數選取
用Matlab提供的濾波器設計的專(zhuān)門(mén)工具箱———FDATool仿真設計濾波器,滿(mǎn)足要求的FlR濾波器幅頻特性如圖1,由于浮點(diǎn)小數FPGA中實(shí)現比較困難,且代價(jià)太大,因而需要將濾波器的系數和輸入數據轉化為整數,其中量化后的系數在Matlab主窗口可直接轉化,對于輸入數據,可乘上一定的增益用Altbus控制位寬轉化為整數輸入。
4.2 FIR濾波器模型建立
根據FIR濾波器原理,可以利用FPGA來(lái)實(shí)現FIR濾波電路,DSPBuilder設計流程的第一步是在Matlab/Simulink中進(jìn)行設計輸入,即在Matlab的Simulink環(huán)境建立一個(gè)MDL模型文件,用圖形方式調用AlteraDSPBuilder和其他的Simulink庫中的圖形模塊,構成系統級或算法級設計框圖(或稱(chēng)Simulink建模),如圖2所示。
4.3 基于DSPBuilder的濾波器仿真
輸入信號分別采用頻率f1=8KHz和f2=16KHz的兩個(gè)正弦信號進(jìn)行疊加。其中的仿真波形如圖3所示,從FIR濾波電路的仿真結果看出,輸入信號通過(guò)濾波器后輸出基本上變成單頻率的正弦信號,進(jìn)一步通過(guò)頻譜儀可看出f2得到了較大的抑制,與條件規定的fc=10.8kHz低通濾波器相符合,至此完成了模型仿真。
4.4 運用Modelsim進(jìn)行功能仿真
在Simulink中進(jìn)行的仿真是屬于系統驗證性質(zhì)的,是對MDL文件進(jìn)行的仿真,并沒(méi)有對生成的VHDL代碼進(jìn)行過(guò)仿真。事實(shí)上,生成VHDL描述是RTL級的,是針對具體的硬件結構的,而在Matlab的Simulink中的模型仿真是算法級(系統級)的,是針對算法實(shí)現的,這二者之間有可能存在軟件理解上的差異,轉換后的VHDL代碼實(shí)現可能與MDL模型描述的情況不完全相符,這就是需要針對生成的RTL級VHDL代碼進(jìn)行功能仿真。
在此,筆者利用Modelsim對生成的VHDL代碼進(jìn)行功能仿真。設置輸入輸出信號均為模擬形式,出現如圖4所示的仿真波形,可以看到這與Simulink里的仿真結果基本一致,即可在QuartusⅡ環(huán)境下進(jìn)行硬件設計。
4.5 在FPGA器件中實(shí)現FIR濾波器
在QuartusⅡ環(huán)境中打開(kāi)DSPBuilder建立的QuartusⅡ項目文件firl.qpf。在QuartusⅡ中進(jìn)行再一次仿真,由此可以看到符合要求時(shí)序波形,然后指定器件引腳并進(jìn)行編譯,最后下載到FPGA器件中,就可以對硬件進(jìn)行測試,加上CLCOK信號和使能信號,用信號發(fā)生器產(chǎn)生所要求的兩個(gè)不同頻率的正弦信號,就可以在示波器上看到濾波以后的結果,需要設計不同的濾波器電路時(shí),僅修改FIR濾波模型文件就可以實(shí)現,這樣不僅避免了繁瑣的VHDL語(yǔ)言編程,而且便于進(jìn)行調整。
5 結束語(yǔ)
在利用FPGA進(jìn)行數字濾波器的開(kāi)發(fā)時(shí),采用DSPBuilder作為設計工具能加快進(jìn)度。當然,在實(shí)際應用中,受精度、速度和器件選擇方面的影響,可以對其轉化的VHDL進(jìn)行進(jìn)一步的優(yōu)化。