喬俊松　李琪　陳誠(chéng)

摘要：在大型吊裝設(shè)備工作時(shí)，需要采集設(shè)備的振動(dòng)信號(hào)和電機(jī)轉(zhuǎn)速信號(hào)作為故障診斷的數(shù)據(jù)來(lái)源，但由于不同電機(jī)轉(zhuǎn)速下混入的高頻噪聲信號(hào)頻率不相同，導(dǎo)致傳統(tǒng)低通濾波器無(wú)法得到信噪比較高的信號(hào)，對(duì)故障診斷準(zhǔn)確率造成影響，急需設(shè)計(jì)可靈活修改各種參數(shù)的低通濾波器。使用MATLAB設(shè)計(jì)了FIR數(shù)字低通濾波器，而后使用QuartusⅡ軟件在FPGA上實(shí)現(xiàn)該濾波器。使用Modelsim軟件對(duì)FIR濾波器進(jìn)行波形仿真，并把數(shù)據(jù)導(dǎo)出至MATLAB分析。結(jié)果表明，所設(shè)計(jì)的FIR低通濾波器具有良好的濾波效果，滿足了預(yù)期的設(shè)計(jì)要求，也可方便修改各種濾波參數(shù)，達(dá)到了預(yù)期的設(shè)計(jì)目的。

關(guān)鍵詞：低通濾波;FIR數(shù)字濾波器;FPGA;MATLAB

中圖分類(lèi)號(hào)：TU745

文獻(xiàn)標(biāo)識(shí)碼：A

DOI： 10.15913/j.cnki.kjycx.2019.09.001

1 引言

軍事裝備中常用各種大型吊裝設(shè)備來(lái)完成武器裝備的吊裝和安放，設(shè)備的工作運(yùn)行狀態(tài)需要實(shí)時(shí)監(jiān)控，通過(guò)狀態(tài)監(jiān)控和故障診斷預(yù)判提高識(shí)別問(wèn)題的概率，確保武器裝備的安全。傳統(tǒng)故障診斷設(shè)備采集設(shè)備開(kāi)車(chē)后的電機(jī)轉(zhuǎn)速和振動(dòng)信號(hào)，通過(guò)故障診斷算法完成判斷，而采集到的振動(dòng)信號(hào)包含大量的高頻干擾信號(hào)，傳統(tǒng)濾波器設(shè)置10倍工頻作為截止頻率，在轉(zhuǎn)速達(dá)到6 000 r/min以上時(shí)可較好地濾除高頻干擾，而在轉(zhuǎn)速在2 000 -6 000 r/min時(shí)由于截止頻率過(guò)高，無(wú)法達(dá)到濾除效果。而數(shù)字FIR濾波器如能方便修改參數(shù)矩陣，就可以實(shí)現(xiàn)改變各項(xiàng)濾波參數(shù)靈活配置濾波器的要求，為達(dá)到濾波要求，本文摸索出一種方便利用MATLAB和QUARTUS II設(shè)計(jì)FIR濾波器的方法和流程，實(shí)踐證明，能夠改善振動(dòng)信號(hào)去噪能力，可提高信噪比。

2 FIR濾波器的原理及設(shè)計(jì)方法

2.1 FIR濾波器數(shù)學(xué)模型

數(shù)字濾波是將輸入的信號(hào)序列按規(guī)定的算法處理，得到所期望的輸出序列。一個(gè)線性時(shí)不變系統(tǒng)的輸出序列v（n）與輸入序列x（n）之間應(yīng)滿足常系數(shù)線性差分方程：窗、Kaiser窗等，前三種窗函數(shù)都是以增加主瓣寬度為代價(jià)換取一定程度的旁瓣抑制。只有Kaiser窗是通過(guò)調(diào)整參數(shù)值來(lái)折中選擇主瓣寬度和旁瓣衰減，因此具有很高的靈活性。

3 基于MATLAB&QUARTUS¨的FIR濾波器設(shè)計(jì)

3.1

數(shù)字濾波器總體設(shè)計(jì)方案 FIR濾波器設(shè)計(jì)流程如圖2所示。

MATLAB數(shù)學(xué)工具優(yōu)勢(shì)在于不僅可以完成基本的數(shù)學(xué)運(yùn)算，對(duì)信號(hào)進(jìn)行處理，而且根據(jù)各專(zhuān)門(mén)領(lǐng)域中的特殊需要提供了許多可選的工具箱，本文采用其設(shè)計(jì)濾波器專(zhuān)用的Filter Design& Analysis To01（ FDA）工具箱，該工具箱使用方便，操作簡(jiǎn)單，可輸出需要的濾波器參數(shù)數(shù)組。

QUARTUSⅡ是Altera公司的綜合性FPGA開(kāi)發(fā)軟件，內(nèi)嵌自帶的綜合器和仿真器，可以完成從設(shè)計(jì)輸入到硬件配置的完整FPGA設(shè)計(jì)流程。QUARIUSⅡ支持Altera的IP核，用戶可以充分利用成熟的模塊，簡(jiǎn)化設(shè)計(jì)的復(fù)雜性，加快設(shè)計(jì)速度。

設(shè)計(jì)時(shí)將二者結(jié)合使用，發(fā)揮各自的優(yōu)勢(shì)，先利用MATLAB的Filter Design& Analysis To01（ FDA）工具箱進(jìn)行數(shù)字低通濾波器的設(shè)計(jì)，再將濾波器參數(shù)導(dǎo)人QUARTUSⅡ中進(jìn)行編程，整個(gè)設(shè)計(jì)流程如圖2所示。

3.2 基于MATLAB的數(shù)字低通濾波器參數(shù)設(shè)計(jì)

利用MATLAB中的數(shù)字濾波器設(shè)計(jì)模塊FDATOOL進(jìn)行低通濾波器的設(shè)計(jì)，步驟如下。在MATLAB指令窗口中輸入“FDATOOL”，彈出窗口如圖3所示。

根據(jù)整體設(shè)計(jì)要求，將濾波器參數(shù)設(shè)置如下。

類(lèi)型（Filer Type）：低通（Low Pass）。

設(shè)計(jì)方法（Design Method）： FIR。

采用窗函數(shù)法（ Window）：Kaiser窗。

采樣頻率Fs為2 000 Hz。

截止頻率Fc為500 Hz。

濾波器的階數(shù)取16，因此Filter order取15。

阻帶衰減不小于50 dB，beta值取4.5。

數(shù)字濾波器分析方面，F(xiàn)IR濾波器幅頻響應(yīng)如圖4所示。

設(shè)計(jì)好濾波器后，通過(guò)FDA工具中的Analysis按鈕進(jìn)行濾波器分析，啟動(dòng)幅頻響應(yīng)分析如圖4所示，其中x軸為頻率，y軸為幅度值（單位為dB）。

由圖4可以得到，頻率在200 Hz以下的信號(hào)基本沒(méi)有發(fā)生衰減，而頻率在500 Hz及以上的信號(hào)則發(fā)生了較大幅度的衰減，并且不小于50 dB。指標(biāo)達(dá)到了濾波器的設(shè)計(jì)要求。FIR濾波器相頻響應(yīng)分析如圖5所示。由圖5可知設(shè)計(jì)的FIR濾波器在通帶內(nèi)為線性相位響應(yīng)，即該濾波器是一個(gè)線性相位的濾波器。線性相位意味著通帶信號(hào)保真，滿足了設(shè)計(jì)要求。

導(dǎo)出濾波器系數(shù)方面，利用MATLAB的FDA工具箱設(shè)計(jì)完成FIR低通濾波器并且滿足設(shè)計(jì)要求后，再將濾波器系數(shù)，即一個(gè)包含16個(gè)濾波器參數(shù)的數(shù)組導(dǎo)出到TXT文本中供QUARTUSⅡ在設(shè)計(jì)FIR IP核時(shí)調(diào)用。

3.3 基于FPGA的數(shù)字低通濾波器設(shè)計(jì)

在完成MATLAB軟件仿真及數(shù)據(jù)準(zhǔn)備之后，開(kāi)始進(jìn)行FPGA設(shè)計(jì)工作，流程如圖6所示。

建立工程，設(shè)置FIR核參數(shù)，如圖7所示。根據(jù)需要將FIR核的參數(shù)設(shè)置濾波器系數(shù)位寬（Bit Width）為12 bit，目標(biāo)器件（ Device Family）為CycloneⅣE，流水線級(jí)數(shù)（ Pipline Level）為1，輸入數(shù)據(jù)及濾波器系數(shù)的存儲(chǔ)資源保持默認(rèn)值為L(zhǎng)ogic Cells，濾波器結(jié)構(gòu)（Structure）為DistributedArithmetic：Fully parallel。濾波器系數(shù)通過(guò)外部文件（FirCoe.txt）導(dǎo)人，如圖8所示。

編寫(xiě)Verilog HDL程序，實(shí)現(xiàn)濾波功能，主要程序如下：

module FirIPDa （reset_n， clk， Sin， S_out）;

input reset_n;

//復(fù)位信號(hào)，低電平有效

input

clk; //FPGA系統(tǒng)時(shí)鐘2kHz

input signed [11：0]S_in; //數(shù)據(jù)輸入頻率為2kHz

output signed [24：0]S_out.//濾波后的輸出數(shù)據(jù)

wire ast sink_valid， astsourceready，

wire ast_ source_valid;

wire[1：0]ast_sourceerror;

wire [1：0]ast_sink_error;

assign astsink_valid-1fb1;

assign ast_sourceread～l'bl;

assign astsink_error2'd0;

fir firinst（.clk（clk）， .reset_n（reset_n），

.ast sink data（S_in），

.ast_ sink valid（ast_sink_valid），

.ast_ sourceready（ast_sourceready），

.ast sink_ error（ast_sink_error），

.ast source_data（Sout），

.ast sink ready（ast_sink_ready），

.ast_ source_valid（astsource_valid），

.ast_ sourceerror（ast_sourceerror》;

endmodule

利用Modelsim進(jìn)行仿真并將仿真結(jié)果輸出至外部文件中，再使用MATLAB軟件對(duì)測(cè)試結(jié)果進(jìn)行分析對(duì)比，對(duì)比情況如圖9.圖10所示。

圖10輸入信號(hào)與輸出信號(hào)的時(shí)域波形對(duì)比

從圖9中可以看出，800 Hz附近的噪聲信號(hào)衰減都大于50 dB，而200 Hz附近的信號(hào)幾乎沒(méi)有改變，滿足了設(shè)計(jì)要求。從圖10中可以看出，輸入信號(hào)有兩個(gè)峰值，經(jīng)過(guò)濾波后變?yōu)榱藛畏?，并且保留下?lái)那個(gè)信號(hào)的幅度也有所加強(qiáng)，滿足了設(shè)計(jì)要求。

4 結(jié)束語(yǔ)

本文完整地論述了基于MATLAB和FPGA的FIR低通濾波器的設(shè)計(jì)思路和仿真，達(dá)到了對(duì)FIR濾波器設(shè)計(jì)的要求，仿真實(shí)驗(yàn)驗(yàn)證了設(shè)計(jì)的正確性，使得設(shè)備故障診斷信號(hào)源的信噪比大大提高，為后期故障診斷算法運(yùn)行提供高質(zhì)量的信號(hào)。這一方法可以應(yīng)用到各種采用FIR濾波器的數(shù)字系統(tǒng)中，提高信號(hào)濾波效果。

參考文獻(xiàn)：

[l]何蘊(yùn)良，耿淑琴，汪金輝.基于Verilog的FIR數(shù)字濾波器設(shè)計(jì)與仿真[J].現(xiàn)代電子技術(shù)，2016（ 10）： 1-4.

[2]解亞南.基于FPGA的數(shù)字濾波器的設(shè)計(jì)研究[D].青島：青島大學(xué)，2016.

[3]楊峰.基于FPGA的FIR數(shù)字濾波器設(shè)計(jì)與仿真[J].四川文理學(xué)院學(xué)報(bào)，2016，26 （5）： 33-35.

[4]杜勇.數(shù)字濾波器的MATLAB與FPGA實(shí)現(xiàn)（Altera/Verilog版）[M].北京：電子工業(yè)出版社，2015.

[5]曹振吉，何敏.基于FPGA和Matlab的FIR數(shù)字濾波器[J].現(xiàn)代電子技術(shù)，2015，38 （1）： 98-102.

[6]劉東華.Altera系列FPGA芯片口核詳解[M].北京：電子工業(yè)出版社，2014.

[7]汪麗娜.音頻信號(hào)采集系統(tǒng)中數(shù)字濾波器的研究與設(shè)計(jì)[D].蘭州：蘭州交通大學(xué)，2014.

[8]王香，張莉莉.基于FPGA的16階FIR數(shù)字濾波器的設(shè)計(jì)[J].電子世界，2013（ 16）： 148-149。

[9]單文軍，周雪純，李文華.基于FPGA的FIR數(shù)字濾波器設(shè)計(jì)與實(shí)現(xiàn)[J].現(xiàn)代電子技術(shù)，2013，36（14）： 123-126.

[10]張景芝.基于FPGA的數(shù)字濾波器設(shè)計(jì)與實(shí)現(xiàn)[D].石家莊：河北經(jīng)貿(mào)大學(xué)，2012.

[11]劉慶良，盧榮軍，李建清.FIR數(shù)字濾波器的FPGA實(shí)現(xiàn)研究[J].電子設(shè)計(jì)工程，2010（3）：59-61，64.