祝 翔,何劍鋒*,曾家駿,袁兆林

(1. 東華理工大學(xué)信息工程學(xué)院,江西 南昌 330013;2. 江西省放射性地學(xué)大數(shù)據(jù)技術(shù)工程實驗室,江西 南昌 330013;3. 江西省核地學(xué)數(shù)據(jù)科學(xué)與系統(tǒng)工程技術(shù)研究中心,江西 南昌 330013)

1 引言

本文對單指數(shù)核信號梯形成形算法在FPGA的設(shè)計進行了優(yōu)化,在Simulink環(huán)境下,設(shè)計搭建了基于Z變換的一般核脈沖梯形成形平臺和模塊化級聯(lián)型結(jié)構(gòu)核脈沖梯形成形平臺,對兩種平臺結(jié)構(gòu)進行比較分析與抗噪結(jié)構(gòu)性改進,并設(shè)計了一個直接I型FIR低通濾波器,對高速ADC輸入信號進行噪聲濾除,噪聲濾除效果良好,減少了噪聲對測量精度和能量分辨率的影響,最后通過工具自動生成Verilog代碼移植到FPGA上,降低了研究的時間成本。

2 單指數(shù)核信號梯形成形算法公式推導(dǎo)

數(shù)字多道脈沖幅度分析器(DMCA)的設(shè)計中,梯形濾波成形算法是對核信號進行濾波成形處理使用頻率最高的一種方法,因為梯形濾波成形的信號是一個等腰梯形,脈沖前后沿時間相等,脈寬窄,下降速度快,脈沖的上升沿和平頂寬度可以獨立調(diào)節(jié),以適應(yīng)不同成形需求的實驗,能有效抑制彈道虧損[11,12]。理想核脈沖的梯形濾波成形如圖1所示。

圖1 梯形成形示意圖

核輻射信息經(jīng)過閃爍體探測器或半導(dǎo)體探測器等核輻射探測器轉(zhuǎn)換成電信號,這些核脈沖信號幅值不同、隨機分布,但波形相似,其上升速度快、時間很短,下降緩慢、時間較長,頂部尖銳,波形圖形的特征與指數(shù)衰減信號十分類似[13]。

設(shè)前放輸出理想狀態(tài)的單指數(shù)衰減信號,時域表達式為

Vi(t)=Vmax·exp(-t/tao)·u(t)

(1)

式中,Vmax為單指數(shù)衰減信號的幅值;tao為單指數(shù)衰減信號的時間常數(shù);u(t)為標(biāo)準(zhǔn)單位階躍函數(shù)。以Ts為采樣周期進行采樣,則有脈沖序列表達式為

Vi(t)=Vmax·exp(-n·Ts/tao)·u(t)

(2)

對(1)式作Z變換,且令d1=exp(-n·Ts/tao)則有

Vi(z)=z·Vmax/(z·d1)

(3)

圖1中,Vmax為輸出梯形信號的高度,ta為輸出梯形信號的上升沿,tb為輸出脈沖的寬度,d為梯形的平頂寬度,d=tb-ta;用分段函數(shù)表示梯形信號

從高等教育領(lǐng)域來分析，大數(shù)據(jù)在提供教育決策支持、推進教育教學(xué)改革以及深化教育管理變革等方面都起到了顯著的作用。目前各教育主管部門和高校由于重視程度、基礎(chǔ)條件、研究能力以及資金投入不同，大數(shù)據(jù)的應(yīng)用發(fā)展水平也大不相同。如果要加強大數(shù)據(jù)在高等教育領(lǐng)域中的應(yīng)用，充分發(fā)揮教育領(lǐng)域中海量數(shù)據(jù)的真實價值，應(yīng)該需要從以下幾個方面繼續(xù)努力。

(4)

其中,y1(t)=Vmax·t;y2(t)=-y1(t-ta);y3(t)=-y1(t-tb);y4(t)=-y1(t-tc)。把上述子式代入式(4)做單邊Z變換。

則有表達式為

(5)

其中,na=ta/Ts,nb=tb/Ts,nc=tc/Ts。

那么,由表達式(3)和表達式(5)則可以得到梯形濾波成形算法的Z變換函數(shù)

(6)

式中,na、nb為調(diào)節(jié)參數(shù)。

對式(6)進行Z逆變換,則可求得其時域表達式

V0(n)=2·V0(n-1)-V0(n-2)+

+Vi(n-nc-1)-d1·[Vi(n-2)-Vi(n-na-2)-Vi(n-nb-2)+Vi(n-nc-2)]}

(7)

根據(jù)時域表達式(7)就可以在MATLAB中編寫腳本程序進行仿真,并在Simulink環(huán)境下搭建核脈沖梯形成形平臺進行系統(tǒng)級仿真驗證?？梢酝ㄟ^改變na、nb的值在仿真中實現(xiàn)對梯形成形信號的調(diào)節(jié)。

3 梯形成形算法MATLAB仿真與模型改進

根據(jù)單指數(shù)衰減脈沖信號式(1)在MATLAB R2017b中編寫梯形成形濾波腳本程序,以單指數(shù)衰減信號來模擬探測器的輸出信號,仿真后得到效果如圖2所示。考慮到實際情況中,探測器輸出的信號會存在一定的噪聲干擾,電子學(xué)噪聲和高速ADC實時采樣高頻噪聲與有用信號疊加后會降低測量精度,影響能量分辨率[14,15]。核電子學(xué)噪聲產(chǎn)生影響的主要是熱噪聲和散粒噪聲,這兩者都是呈高斯正態(tài)分布的,因此,為了使仿真實驗更貼近現(xiàn)實,在單指數(shù)衰減信號上疊加了一個高斯加性白噪聲來模仿現(xiàn)實噪聲所產(chǎn)生的影響,仿真波形如圖3所示。

圖3 加噪聲信號梯形成形MATLAB仿真

經(jīng)過腳本程序進行模擬之后,在Simulink中建立梯形濾波成形平臺,建立了直接型和模塊化級聯(lián)型兩種脈沖梯形成形平臺,讀取疊加噪聲的脈沖信號查看波形,Vmax=5,Vmax=3,Vmax=6,取量化時間ta=50,tb=100,得到濾波成形的波形如圖4所示。腳本程序模擬出的良好波形在Simulink中噪聲更加明顯,對脈沖信號的梯形濾波成形帶來影響,產(chǎn)生了基線誤差。在直接型梯形成形平臺中,梯形成形曲線不夠平滑,且出現(xiàn)了基線漂移,而在初始設(shè)計的模塊化級聯(lián)型梯形成形平臺中,雖然噪聲的影響比直接型平臺的小,但是也產(chǎn)生了基線漂移。這是因為在模型結(jié)構(gòu)上,沒有對噪聲進行優(yōu)化處理,差分運算部分靠前而非差分運算部分在后,這樣的結(jié)構(gòu)導(dǎo)致了濾波成形平臺的整體結(jié)構(gòu)抗噪性弱[16]。

圖4 兩種梯形成形平臺濾波效果對比

綜合以上原因,對上述方法做出了一些改進調(diào)整。首先,在模擬的梯形成形算法上,引入FIR濾波,在采樣頻率為50MHz的情況下,經(jīng)過多次實驗,確定了較好的數(shù)字截止頻率ωn=0.45,基于MATLAB的fir1函數(shù),求出濾波系數(shù),再用MATLAB函數(shù)庫中的filter函數(shù)濾波,設(shè)計了一個長度為32位的低通濾波器,濾除噪聲,降低噪聲產(chǎn)生的負(fù)面影響。其次,在梯形成形平臺設(shè)計中,利用filterDesigner工具箱,選用等波紋的設(shè)計方法,設(shè)置采樣頻率為50MHz,設(shè)計了一個7階等波紋FIR低通濾波器。由于其是有限單位脈沖響應(yīng),且具有很好的選頻特性,在濾除噪聲的同時擁有較好的穩(wěn)定性[17]。經(jīng)過FIR數(shù)字濾波器后,輸入信號的噪聲得到很好的濾除,噪聲濾除前后效果對比如圖5所示。

此外,在梯形成形平臺結(jié)構(gòu)設(shè)計上采用模塊化設(shè)計,設(shè)計了級聯(lián)型濾波平臺的同時,改變了差分結(jié)構(gòu)的位置,將整個負(fù)反饋結(jié)構(gòu)的位置前移,削弱了噪聲累積,使得基線漂移的情況得到改善。整體改進后的級聯(lián)型平臺梯形成形效果如圖6所示?；€漂移得到改善,對輸入信號濾除噪聲后進行的梯形濾波成形,噪聲得到明顯的抑制,濾波成形的梯形效果良好。

3 FPGA硬件實現(xiàn)與分析

采取級聯(lián)型結(jié)構(gòu)實現(xiàn)梯形成形濾波器的設(shè)計與驗證,其本質(zhì)上是由級數(shù)低的濾波器級聯(lián)形成級數(shù)高的濾波器,每個濾波器都是獨立的模塊,并且將前一級的輸出作為后一級的輸入。如果由直接型結(jié)構(gòu)實現(xiàn)梯形成形濾波器的設(shè)計,資源耗費大、結(jié)構(gòu)復(fù)雜,并且數(shù)字處理有限精度的情況下,系統(tǒng)穩(wěn)定性較差。采用級聯(lián)型結(jié)構(gòu)設(shè)計后,降低了資源消耗,優(yōu)化了模型結(jié)構(gòu),提升了運算速度,并且系統(tǒng)更加穩(wěn)定。設(shè)計好的級聯(lián)型梯形成形濾波器模塊化結(jié)構(gòu)如圖7所示。

圖7 梯形成形濾波器模塊化結(jié)構(gòu)

輸入信號從Input進入后先經(jīng)過FIR低通濾波器,濾除噪聲后進入第一級模塊做延時處理,然后將第一級模塊的輸出作為后一級模塊的輸入,以此類推直至第五級模塊結(jié)束,最后由Output輸出梯形成形結(jié)果,成形波形可添加一個示波器模塊Scope進行觀察。

FIR低通濾波器模塊是對輸入信號進行噪聲濾除,其輸出信號的噪聲大大降低,噪聲的濾除效果如圖5所示。FIR低通濾波器輸出的信號作為第一級模塊的的輸入,第一級模塊對濾除噪聲的核脈沖信號進行nb個延時,輸出信號就是從低通濾波器輸出的信號減去nb個延時后的信號,輸出波形是在濾除噪聲信號波形后nb處將濾除噪聲信號進行翻轉(zhuǎn)。第二級模塊將第一級模塊的輸出作為輸入,并將這個信號延時一個時鐘周期再乘以d＿1,輸出信號則由將第一級模塊輸出的信號減去延時一個時鐘周期再乘d＿1的信號所得到,仿真輸出的波形則是第一級模塊輸出波形的幅度峰值。第三級模塊的結(jié)構(gòu)與第一級模塊相同,輸出信號是第二級模塊輸出的信號與將第二級模塊輸出的信號作na個延時后的信號做差得到的結(jié)果。第三級模塊的輸出作為第四級模塊的輸入,輸出信號則是輸入信號加上延時一個周期輸入信號的兩倍增益與輸入信號延時兩個周期的差,第四級模塊是一個負(fù)反饋結(jié)構(gòu),不能夠輕易調(diào)換位置,否則會產(chǎn)生基線漂移。第五級模塊則是將第四級模塊輸出的信號除以na,輸出波形就是將本季模塊的輸入信號波形縮小na倍,也就是最后成形的梯形。

在Simulink設(shè)計好梯形濾波成形平臺后,在Code選項卡調(diào)用HDL Workflow Advisor生成工程文件,用Quartus和Model-sim仿真軟件打開工程文件對梯形成形效果作仿真驗證。在生成工程時,硬件語言可以選擇為VHDL或者Verilog HDL,可以根據(jù)實驗需求進行選擇。選用Verilog HDL語言生成工程后,由Quartus開發(fā)工具對整個工程進行編譯、綜合、布局布線,直至在FPGA上完成整個工作流程。用Model-sim仿真得到梯形成形波形如圖8所示,得到的波形與MATLAB仿真測試波形的效果基本一致,FIR低通濾波器對噪聲濾除效果好。

圖8 梯形濾波成形Modelsim仿真

模擬驗證通過后,在Quartus Prime中調(diào)用工具欄的Signal Tap Logic Analyzer 進行板級調(diào)試,調(diào)試時設(shè)置好所需參數(shù),添加捕捉的波形,并設(shè)置觸發(fā)條件,設(shè)置完成后進行一次綜合編譯,最后將FPGA硬件開發(fā)板上電,選定設(shè)備為開發(fā)板,將綜合后的程序下載到FPGA硬件,傳輸數(shù)據(jù)并捕捉波形。本次實驗選用的是Altera公司的EP4CE10F17C8芯片完實現(xiàn)多道數(shù)字濾波成形算法的設(shè)計。采用的是CIT-3000SY X 熒光元素錄井儀實測的部分原始數(shù)據(jù)進行實際測試,通過Signal Tap對梯形成形濾波器輸入輸出的數(shù)據(jù)進行捕捉和觀察驗證,捕捉的波形如圖9所示。波形捕捉準(zhǔn)確,最終完成了數(shù)字多道梯形成形算法在FPGA的實現(xiàn)與改進。

圖9 梯形濾波成形Signal Tap 時序驗證圖

4 結(jié)論

本文對單指數(shù)核脈沖的級聯(lián)型結(jié)構(gòu)梯形濾波成形平臺改進后,提升了結(jié)構(gòu)抗噪性,設(shè)計的直接型FIR低通濾波器對噪聲能進行有效濾除,降低了噪聲的干擾;對單指數(shù)核脈沖作梯形成形仿真,確定了合適的成形參數(shù),能對成形脈沖進行調(diào)節(jié)從而適應(yīng)需求;使用Simulink的HDL Code工具后,整體的設(shè)計使得FPGA開發(fā)流程得以簡化,硬件資源消耗降低,通過Signal Tap Logic Analyzer 內(nèi)嵌邏輯分析儀對實際數(shù)據(jù)的成形波形進行觀察分析,結(jié)果與仿真結(jié)果基本一致,取得了不錯的測試結(jié)果,對數(shù)字多道脈沖幅度分析器后續(xù)的設(shè)計工作具有一定幫助。