基于DSP Builder的14階FIR濾波器的設(shè)計(jì)

時(shí)間：2008-06-16 10:07:00

關(guān)鍵字： DSP BUILDER FIR濾波器 VHDL

手機(jī)看文章

掃描二維碼
隨時(shí)隨地手機(jī)看文章

[導(dǎo)讀]數(shù)字濾波器在數(shù)字信號處理的各種應(yīng)用中發(fā)揮著十分重要的作用，他是通過對采樣數(shù)據(jù)信號進(jìn)行數(shù)學(xué)運(yùn)算處理來達(dá)到頻域?yàn)V波的目的。數(shù)字濾波器既可以是有限長單脈沖響應(yīng)(FIR)濾波器也可以是無限長單脈沖響應(yīng)(IIR)濾波器。在維納濾波器理論發(fā)明的早期，人們使用IIR濾波器，但現(xiàn)在更多是使用FIR濾波器。本文按照Matlab／Simulink／DSP Builder／QuartusⅡ流程，設(shè)計(jì)一個(gè)FIR濾波器。Altera DSP Builder是連接Simulink和QuartusⅡ開發(fā)軟件的DSP開發(fā)工具。在DSP Builder的無縫設(shè)計(jì)流程中，首先在Matlab軟件中進(jìn)行算法設(shè)計(jì)，然后在Simulink軟件中進(jìn)行系統(tǒng)集成，最后將設(shè)計(jì)輸出為硬件描述語言(HDL)文件，以便在QuartusⅡ軟件中使用。

數(shù)字濾波器在數(shù)字信號處理的各種應(yīng)用中發(fā)揮著十分重要的作用，他是通過對采樣數(shù)據(jù)信號進(jìn)行數(shù)學(xué)運(yùn)算處理來達(dá)到頻域?yàn)V波的目的。數(shù)字濾波器既可以是有限長單脈沖響應(yīng)(FIR)濾波器也可以是無限長單脈沖響應(yīng)(IIR)濾波器。在維納濾波器理論發(fā)明的早期，人們使用IIR濾波器，但現(xiàn)在更多是使用FIR濾波器。本文按照Matlab／Simulink／DSP Builder／QuartusⅡ流程，設(shè)計(jì)一個(gè)FIR濾波器。Altera DSP Builder是連接Simulink和QuartusⅡ開發(fā)軟件的DSP開發(fā)工具。在DSP Builder的無縫設(shè)計(jì)流程中，首先在Matlab軟件中進(jìn)行算法設(shè)計(jì)，然后在Simulink軟件中進(jìn)行系統(tǒng)集成，最后將設(shè)計(jì)輸出為硬件描述語言(HDL)文件，以便在QuartusⅡ軟件中使用。

1 FIR濾波器介紹

1.1 FIR濾波器原理

對于一個(gè)FIR濾波器系統(tǒng)，他的沖激響應(yīng)總是有限長的，其系統(tǒng)函數(shù)可以記為：

其中：x(n)是輸入采樣序列；h(n)是濾波器系數(shù)；L是濾波器的階數(shù)；y(n)表示濾波器的輸出序列。也可以用卷積來表示輸出序列y(n)與x(n)，h(n)的關(guān)系。

y(n)=x(n)*h(n)

典型的直接I型FIR濾波器如圖1所示，其輸出序列y(n)滿足下式：

1.2 設(shè)計(jì)要求

數(shù)字濾波器實(shí)際上是一個(gè)采用有限精度算法實(shí)現(xiàn)的線性非時(shí)變離散系統(tǒng)，他的設(shè)計(jì)步驟是先根據(jù)需要確定其性能指標(biāo)，設(shè)計(jì)一個(gè)系統(tǒng)函數(shù)h(n)逼近所需要的性能指標(biāo)，濾波器的系數(shù)的計(jì)算可以借助Matlab強(qiáng)大的計(jì)算能力和現(xiàn)成的濾波器設(shè)計(jì)工具來完成，最后采用有限的精度算法實(shí)現(xiàn)。本系統(tǒng)的設(shè)計(jì)指標(biāo)為，設(shè)計(jì)一個(gè)14階的FIR濾波器。1 MHz和16 MHz的兩個(gè)正弦波的合成波形，經(jīng)過這個(gè)濾波器，濾除16 MHz的成分。

2 DSP Builder設(shè)計(jì)流程

使用DSP Builder完成設(shè)計(jì)時(shí)，首先在Matlab／Simulink軟件中建立模型文件(.mdl)，DSP Builder SignalCom-piler模塊讀取由DSP Builder和MegaCore模塊構(gòu)建的Simulink建模文件(.mdl)，生成VHDL文件和工具命令語言(Tcl)腳本，進(jìn)行綜合、硬件實(shí)施和仿真。

3 FIR數(shù)字濾波器的DSP Builder設(shè)計(jì)

3.1 Matlab／Simulink建模

根據(jù)FIR數(shù)字濾波器的算法，在Matlab的Simulink環(huán)境中建立一個(gè)MDL模型文件。用圖形方式調(diào)用Alter-aDSPBuilder和其他的Simulink庫中的圖形模塊。

3.2 在Simulink中仿真并生成VHDL代碼

完成模型設(shè)計(jì)之后，可以先在Simulink中對模型進(jìn)行仿真，可以通過Simulink中的示波器模塊查看各步驟的中間結(jié)果。輸入信號采用頻率為1 MHz和16 MHz的兩個(gè)正弦波的疊加。仿真結(jié)果如圖2所示。從仿真波形可以看出，經(jīng)過FIR濾波器之后，16 MHz的高頻信號被很好地濾除了。

3.3 使用Modelsim進(jìn)行RTL仿真

在Simulink中進(jìn)行的仿真是屬于系統(tǒng)驗(yàn)證性質(zhì)的，是對mdl文件進(jìn)行仿真，并沒有對生成的VHDL代碼進(jìn)行過仿真。事實(shí)上，生成VHDL描述的是RTL級的，是針對具體的硬件結(jié)構(gòu)的，而在Matlab的Simulink中的模型仿真是算法級的，兩者之間有可能存在軟件理解上的差異。轉(zhuǎn)換后的VHDL代碼實(shí)現(xiàn)可能與mdl模型描述的情況不完全相符。這就需要針對生成的RTL級VHDL代碼進(jìn)行功能仿真，如圖3所示。

仿真結(jié)果說明，本濾波器可以較好地起到濾波的效果，由于輸入正弦波經(jīng)過數(shù)字化過程之后在Modelsim中還原出的模擬波形出現(xiàn)失真現(xiàn)象。

3.4 使用QuartusⅡ進(jìn)行時(shí)序仿真和綜合編譯

Modelsim完成的RTL級仿真只是功能仿真，其仿真結(jié)果并不能精確地反映電路的全部硬件特性，進(jìn)行門級的時(shí)序仿真仍然是十分重要的。時(shí)序仿真即針對具體硬件芯片，在仿真時(shí)加入芯片內(nèi)部的時(shí)間信息包括：建立時(shí)間，保持時(shí)間，關(guān)鍵路徑延時(shí)等。綜合編譯后，tsu=4.622 ns，tco=6.925 ns，th=0.415 ns，時(shí)間分析符合實(shí)際情況。

4 濾波器設(shè)計(jì)在FPGA上的實(shí)現(xiàn)

由于從Simulink建模仿真，到Modelsim RTL仿真和QuartusⅡ時(shí)序仿真，使用的輸入正弦波都是仿真信號，而不是實(shí)際的信號源。在硬件實(shí)際運(yùn)行時(shí)，可以從外部信號源接入芯片內(nèi)部或者在芯片內(nèi)部存儲(chǔ)正弦波的數(shù)據(jù)。這里采用的是后者，即在頂層文件中引入LPM_ROM宏模塊，在其中存入正弦波數(shù)據(jù)的mif文件(存儲(chǔ)初始化文件)，F(xiàn)IR濾波器模塊直接從ROM中讀取數(shù)據(jù)，經(jīng)過濾波處理之后數(shù)據(jù)存入另一ROM中，可以通過In-systemmemory content editor工具從中讀出數(shù)據(jù)，再利用Matlab還原出波形，從而在實(shí)際硬件中驗(yàn)證了濾波器的功能。

實(shí)現(xiàn)的頂層圖如圖4所示。

5 結(jié) 語

本文采用Matlab／Simulink／DSP Builder／QuartusⅡ的設(shè)計(jì)流程，實(shí)現(xiàn)了14階FIR濾波器。通過系統(tǒng)仿真，RTL仿真，時(shí)序仿真和實(shí)際硬件測試(這里采用了一種不需要使用SignaltapⅡ軟件邏輯分析儀的方法)對濾波器的功能、性能等進(jìn)行了分析，從而驗(yàn)證了采用DSP，Builder實(shí)現(xiàn)濾波器設(shè)計(jì)的硬件化的獨(dú)特優(yōu)勢。