洪 旭 倪師軍 周建斌 馬英杰 劉 易,2 周 偉,2

1（成都理工大學(xué) 核技術(shù)與自動化工程學(xué)院 成都 610059）2（東華理工大學(xué) 放射性地質(zhì)與勘探技術(shù)國防重點實驗室 南昌 330013）

數(shù)字高斯脈沖成形算法仿真研究

洪 旭1倪師軍1周建斌1馬英杰1劉 易1,2周 偉1,2

1（成都理工大學(xué) 核技術(shù)與自動化工程學(xué)院 成都 610059）2（東華理工大學(xué) 放射性地質(zhì)與勘探技術(shù)國防重點實驗室 南昌 330013）

高斯脈沖具有良好的時間、頻率響應(yīng)和較高的信噪比，探測器輸出信號通常被濾波成形為高斯波形或類高斯波形。在高斯脈沖成形算法基礎(chǔ)上，引入截止頻率和品質(zhì)因子，討論成形參數(shù)對成形脈沖幅度、寬度的影響。通過與梯形脈沖成形算法對比，研究高斯脈沖成形算法的濾波效果以及對核信號幅度譜的作用。結(jié)果表明，在相同達(dá)峰時間條件下高斯脈沖成形算法具有更好的噪聲抑制能力，同時可減小核信號幅度譜的半高寬；結(jié)合使用梯形脈沖成形算法可改善高斯脈沖成形算法在處理堆積脈沖中的不足，提高核信號通過率。

高斯脈沖成形，梯形脈沖成形，核信號，半高寬

20世紀(jì)90年代以來，隨著電子技術(shù)的發(fā)展，傳統(tǒng)的模擬核儀器系統(tǒng)逐漸向數(shù)字化方向發(fā)展。數(shù)字核儀器系統(tǒng)結(jié)構(gòu)如圖1所示，探測器輸出信號經(jīng)前端電路濾波成形后，進(jìn)入模數(shù)轉(zhuǎn)換器(Analog-todigital Converter, ADC)采樣，隨后被直接被送入現(xiàn)場可編程門陣列(Field-Programmable Gate Array, FPGA)處理單元進(jìn)行數(shù)字處理及幅度甄別，最后通過微控制單元(Microcontroller Unit, MCU)將譜數(shù)據(jù)傳送至計算機進(jìn)行分析。除在結(jié)構(gòu)上簡化了模擬核儀器系統(tǒng)，數(shù)字核儀器系統(tǒng)在核信號濾波成形方面采用了數(shù)字算法，在實時性、靈活性、穩(wěn)定性等方面都有了很大改進(jìn)。由于高速ADC的使用，疊加在信號中的噪聲也被采集，為減小噪聲對測量結(jié)果的影響，需要對ADC采樣后的數(shù)字信號進(jìn)行數(shù)字濾波成形處理。常用的數(shù)字濾波成形方法包括梯形脈沖成形和高斯脈沖成形。梯形脈沖成形兼顧能量分辨率和計數(shù)率[1-5]，被廣泛應(yīng)用于數(shù)字核儀器系統(tǒng)中。高斯脈沖具有良好的時間響應(yīng)、頻率響應(yīng)和較高的信噪比，在數(shù)字核儀器系統(tǒng)中，前端電路通常將探測器輸出信號濾波成形為高斯波形或類高斯波形[6]。例如，Hallgren等[7-10]采用Sallen-Key (S-K)濾波器將探測器輸出信號成形為高斯波形。另外，利用CR-nRC濾波器可得到類高斯波形[11-13]。除了用電子元器件實現(xiàn)高斯脈沖成形外，部分學(xué)者對高斯脈沖成形數(shù)字實現(xiàn)方法也進(jìn)行了研究。陳世國等[14-15]研究小波分析方法實現(xiàn)了指數(shù)衰減信號高斯脈沖成形；Nakhostin[16]通過CR-nRC濾波器，得到類高斯脈沖成形的遞推算法。本文從模擬高斯脈沖成形電路出發(fā)，在已有的高斯脈沖成形算法基礎(chǔ)上，新引入截止頻率和品質(zhì)因子，研究其成形參數(shù)的選擇及在數(shù)字核儀器系統(tǒng)中的應(yīng)用。

圖1 數(shù)字核儀器系統(tǒng)Fig.1 Digital nuclear spectrometer system.

1 高斯脈沖成形算法

低通Sallen-Key濾波器在核電子學(xué)中被廣泛應(yīng)用于核信號的濾波成形，它可以在較少的級數(shù)下得到類高斯波形的輸出，并且有較大的品質(zhì)因子。圖2是低通S-K濾波器電路原理圖。根據(jù)基爾霍夫電流定律可以建立低通S-K濾波器的電流等式，如式(1)所示。

圖2 低通S-K濾波器電路原理圖Fig.2 Scheme of low-pass S-K filter.

式中：Vi、Vo、Vn、Vp和Vf分別表示對應(yīng)節(jié)點的電壓值；R1、R2、R3、R4為對應(yīng)電阻的阻值；C1、C2為對應(yīng)電容的容值。整理可得：

式中：K為成形脈沖的幅度放大倍數(shù)，其值的選擇與輸出信號的品質(zhì)因子有關(guān)，K=(1+R4/R3)。本文新引入截止頻率fc和品質(zhì)因子Q，研究成形參數(shù)的選擇，即：

為實現(xiàn)截止頻率和品質(zhì)因子單獨控制，對圖2中R1、R2、C1、C2進(jìn)行簡化。令R1=R2=R，C1=C2=C，τ=RC，則式(2)可簡化為：

在原有低通S-K濾波器數(shù)學(xué)模型基礎(chǔ)上[17]，由式(3)、(4)、(5)可得到含有截止頻率和品質(zhì)因子的高斯脈沖成形算法，如式(6) ：

式中：τ為成形脈沖的成形時間，其值的選擇與濾波器的截止頻率有關(guān)。

2 高斯脈沖成形算法實驗測試

2.1成形參數(shù)對輸出信號波形的影響

為研究τ和K值對輸出信號的影響，采用實測核信號進(jìn)行實驗。實測信號由硅漂移半導(dǎo)體探測器(Silicon Drift Detector, SDD)探測，經(jīng)高通濾波和線性放大后被轉(zhuǎn)換速率為20×106s-1的ADC采樣得到。選擇不同的τ和K值進(jìn)行實驗，實驗結(jié)果如圖3所示。

圖3 不同成形時間(a)和不同幅度放大倍數(shù)(b)下高斯脈沖成形Fig.3 Gaussian pulse shaping with different shaping times (a) and different amplification factors (b).

由圖3(a)可以看出，當(dāng)K值不變時(K=1.75，Q=0.8)，隨著τ值增大，成形脈沖下降沿變緩，幅度降低，寬度變寬，極點后移；圖3(b)則說明，當(dāng)τ值不變時(τ=65，fc=50 kHz)，隨著K值增大，成形脈沖波形下降沿變快，幅度增大，寬度變窄，極點前移。

結(jié)合式(6)可知，K值的選擇受輸出信號品質(zhì)因子的約束，實驗證明當(dāng)K接近于3時，輸出信號出現(xiàn)振蕩。在實際應(yīng)用中，K值的選取應(yīng)在2附近，以得到最大品質(zhì)因子的輸出信號。τ值與濾波器的截止頻率有關(guān)，當(dāng)τ值較大時雖然對輸入信號的濾波能力越強，但增加了脈沖寬度。應(yīng)用時可通過比較濾波前后輸入信號的頻譜圖確定最優(yōu)τ值，同時τ值較小時（圖3(a)中τ=25），成形脈沖無下沖。

2.2噪聲抑制

為研究高斯脈沖成形算法在濾波方面的性能，采用常用的梯形脈沖成形算法作對比研究。核信號用負(fù)指數(shù)函數(shù)模擬，并疊加高斯白噪聲。梯形脈沖成形采用Imperiale提出的算法[3]，如式(7)所示。

式中：Vi、Vo分別表示輸入、輸出信號；na、nc分別表示梯形脈沖的上升時間和脈沖寬度；nb表示梯形脈沖上升時間與平頂寬度之和，即nc=na+nb；Ts表示ADC的采樣率；τi為輸入核信號的衰減時間常數(shù)。兩種濾波成形算法得到的成形脈沖具有相同達(dá)峰時間。

模擬核信號的幅度為2V，衰減時間常數(shù)為5μs，分別進(jìn)行高斯脈沖成形和梯形脈沖成形處理，結(jié)果如圖4所示。第一組為核信號和對應(yīng)的高斯成形、梯形成形脈沖結(jié)果；第二組為對應(yīng)信號的頻譜圖。

由圖4中成形脈沖的頻譜圖可以看出，隨著τ、na值的增大，兩種方法的濾波效果均變好；對比兩種方法得到成形脈沖的頻譜圖可得，在相同達(dá)峰時間條件下，高斯脈沖成形算法較梯形脈沖成形算法有更好的噪聲抑制能力。

利用MATLAB生成2×105個帶高斯白噪聲的負(fù)指數(shù)信號模擬原始核信號，分別做兩組實驗：

1) 對原始核信號先進(jìn)行高斯脈沖成形，再做幅度甄別，成形參數(shù)τ=15，K=2（Q=1，使得輸出信號品質(zhì)最佳）；

2) 對原始核信號先進(jìn)行相同達(dá)峰時間的梯形脈沖成形，再做幅度甄別，成形參數(shù)為na=30，nb=50。

實驗所得幅度譜如圖5所示。從圖5中可以看出，在不考慮堆積脈沖的情況下，采用高斯脈沖成形算法所得幅度譜的半高寬(Full Width at Half Maximum, FWHM)較采用梯形脈沖成形算法所得幅度譜的半高寬值提高了7.9%，高斯脈沖成形算法適用于核信號濾波處理。

圖4 高斯脈沖成形與梯形脈沖成形(a) τ=25，K=1.5，na=50，nb=100，(b) τ=50，K=1.75，na=100，nb=150，(c) τ=100，K=2，na=200，nb=250Fig.4 Gaussian pulse shaping and trapezoidal pulse shaping. (a) τ=25, K=1.5, na=50, nb=100, (b) τ=50, K=1.75, na=100, nb=150, (c) τ=100, K=2, na=200, nb=250

圖5 高斯脈沖成形與梯形脈沖成形半高寬對比Fig.5 Comparison of Gaussian pulse shaping and trapezoidal pulse shaping on FWHM.

2.3堆積脈沖分離

計數(shù)率是衡量核測量系統(tǒng)另一個重要指標(biāo)，其主要受堆積脈沖的影響。通過在第一個核信號結(jié)束前疊加相同的核信號，模擬堆積脈沖。分別利用高斯脈沖成形和梯形脈沖成形處理堆積脈沖。模擬核信號的幅度為2V，衰減時間常數(shù)為5μs，兩種方法得到的成形脈沖達(dá)峰時間均為2.5μs，結(jié)果如圖6所示。

由圖6可以看出，相同達(dá)峰時間條件下，高斯濾波成形不能分離堆積脈沖；梯形脈沖成形算法在分離堆積脈沖的同時，恢復(fù)第二脈沖的真實幅度。

圖6 高斯脈沖成形(a)和梯形脈沖成形(b)分離堆積脈沖Fig.6 Pile-up pulse separation with Gaussian pulse shaping (a) and trapezoidal pulse shaping (b).

在原始核信號中增加10萬個堆積脈沖，采用與§2.2相同方法研究高斯脈沖成形算法在分離堆積脈沖中的應(yīng)用。高斯脈沖成形參數(shù)τ=15，K=2；梯形脈沖成形參數(shù)na=30，nb=30，所得幅度譜如圖7所示。由圖7可以看出，采用高斯脈沖成形算法所得幅度譜的總計數(shù)明顯低于采用梯形脈沖成形算法所得到的幅度譜，梯形脈沖成形算法適用于分離堆積脈沖。

圖7 高斯脈沖成形與梯形脈沖成形總計數(shù)對比Fig.7Comparison of Gaussian pulse shaping and trapezoidal pulse shaping on total counts.

2.4先濾波再成形

由§2.2和2.3可以看出，高斯脈沖成形與梯形脈沖成形分別在濾波和堆積脈沖分離中各具優(yōu)勢。結(jié)合兩者的優(yōu)點，先利用高斯脈沖成形提高核信號信噪比，然后利用梯形脈沖成形，分離堆積脈沖。從§2中可知，雖然τ值越大對噪聲的抑制能力越好，但成形脈沖的寬度增大，易產(chǎn)生堆積脈沖。

模擬30萬個帶高斯白噪聲的核信號作為原始核信號（其中含有10萬個堆積脈沖），先做高斯脈沖成形處理得到高斯脈沖，然后再對高斯脈沖做梯形脈沖成形處理，最后對梯形脈沖做幅度甄別得到幅度譜。高斯脈沖成形參數(shù)τ=8，K=2；梯形脈沖成形參數(shù)na=30，nb=30，所得幅度譜與單獨采用高斯脈沖成形和梯形脈沖成形結(jié)果對比如圖8所示。

圖8 先高斯脈沖成形再梯形脈沖成形Fig.8 Performance of Gaussian pulse shaping before trapezoidal pulse shaping.

不同成形方法所得幅度譜的半高寬與核信號總計數(shù)如表1所示。由表1可以看出，采用先高斯脈沖成形再做梯形脈沖成形的方法得到的核信號總計數(shù)較單獨采用高斯成形方法所得總計數(shù)提高了2.3倍。在對能量分辨率要求相對較低，而對計數(shù)率要求相對較高的應(yīng)用場合可采用兩種濾波成形方法相結(jié)合的處理方式。

表1 不同方法得到的半高寬與總計數(shù)Table 1 Total counts and energy resolution with different shaping algorithms.

3 結(jié)語

在高斯脈沖成形算法基礎(chǔ)上，引入了截止頻率和品質(zhì)因子，對成形參數(shù)的選擇進(jìn)行了研究。τ與濾波器截止頻率有關(guān)，τ值越大，對輸入信號的高頻噪聲抑制能力越強，同時增加了成形脈沖寬度；K與輸出脈沖品質(zhì)因子有關(guān)，K值越大，成形脈沖品質(zhì)越高，但尾部出現(xiàn)下沖；當(dāng)K=2時，輸出信號具有最大品質(zhì)因子。通過與梯形脈沖成形算法對比研究，高斯脈沖成形算法在濾波方面具有更好的噪聲抑制能力；梯形脈沖成形算法在堆積脈沖分離方面更具優(yōu)勢。在數(shù)字核儀器系統(tǒng)設(shè)計時，對采樣后的核信號先進(jìn)行高斯脈沖成形再進(jìn)行梯形脈沖成形處理可改善系統(tǒng)的計數(shù)率，適用于計數(shù)率要求相對較高的應(yīng)用場合。

1 Jordanov V T, Glenn F. Knoll digital synthesis of pulse shapes in real time for high resolution radiation spectroscopy[J]. Nuclear Instruments and Methods in Physics Research A, 1994,345: 337-345. DOI: 10.1016/ 0168-9002(94)91011-1

2 Warburton W K, Momayezi M, Hubbard-Nelson B, et al. Digital pulse processing: new possibilities in nuclear spectroscopy[J]. Applied Radiation and Isotopes, 2000,53: 913-920. DOI: 10.1016/S0969-8043(00)00247-5

3 Imperiale C, Imperiale A. On nuclear spectrometry pulse digital shaping and processing[J]. Measurement, 2001,30: 49-73. DOI: 10.1016/S0263-2241(00)00057-9

4 Esmaeili-sani V, Moussavi-zarandi A, Akbar-ashrafi N, et al. Neutron-gamma discrimination based on bipolar trapezoidal pulse shaping using FPGAs in NE213[J]. Nuclear Instruments and Methods in Physics Research A, 2012,694: 113-118. DOI: 10.1016/j.nima.2012.08.025

5 Regadio A, Sanchez-Prieto S, Prieto M, et al. Implementation of real-time adaptive digital shaping for nuclear spectroscopy[J]. Nuclear Instruments and Methods in Physics Research A, 2014,735: 297-303. DOI: 10.1016/j.nima.2013.09.063

6 Fernandes A M, Pereira R C, Sousa J, et al. Real time algorithm for digital pulse processing applied to gamma-ray and hard X-ray spectroscopy[J]. Fusion Engineering and Design, 2012,87: 2156-2160. DOI: 10.1016/j.fusengdes.2012.03.033

7 Hallgren B, Bal F, Barr G, et al. The Na48 LKr calorimeter digitizer electronic chain[J]. Nuclear Instruments and Methods in Physics Research A, 1998,419: 680-685. DOI: 10.1016/S0168-9002(98)00846-8

8 Pauly S W. Development of a shaping amplifier, gated integrator and 16K channel ADC for a portable spectroscopy system[J]. Nuclear Instruments and Methods in Physics Research A, 1999,422: 379-384. DOI: 10.1016/S0168-9002(98)00987-5

9 Gevin O, Lugiez F, Limousin O, et al. IDeF-X V1.0: a new 16-channel low-noise analog front-end for Cd(Zn)Te detectors[J]. Nuclear Instruments and Methods in Physics Research A, 2006,567: 140-144. DOI: 10.1016/j.nima. 2006.05.227

10 Bisello D, Candelori A, Giubilato P, et al. Position sensitive detectors for ion electron emission microscopy[J]. Nuclear Instruments and Methods in Physics Research A, 2007,573: 23-26. DOI: 10.1016/ j.nima.2006.10.268

11 Kalinin A I, Bednyakov V A. Pulse shaping for Ge-spectrometers optimized for ballistic deficit and electronic noise[J]. Nuclear Instruments and Methods in Physics Research A, 2005,538: 718-722. DOI: 10.1016/ j.nima.2004.08.126

12 Fleury J, Taille C L, Martin-Chassard G. Front-end ASICs development for W-Si calorimeter at ILC (CALICE collaboration)[J]. Nuclear Instruments and Methods in Physics Research A, 2007,572: 371-375. DOI: 10.1016/ j.nima.2006.10.223

13 Gómez-Galán J A, López-Ahumada R, Sánchez-Rodríguez T, et al. High speed low power FEE for silicon detectors in nuclear physics applications[J]. Nuclear Instruments and Methods in Physics Research A, 2013,714: 155-162. DOI: 10.1016/j.nima.2013.03.002

14 陳世國, 吉世印, 劉萬松. 基于小波分析的指數(shù)衰減信號高斯脈沖成形[J]. 物理學(xué)報, 2008,57(5): 2882-2887. DOI: 10.7498/aps.57.2882

CHEN Shiguo, JI Shiyin, LIU Wansong. Gaussian pulse shaping of exponential decay signal based on wavelet analysis[J]. Acta Physica Sinica, 2008,57(5): 2882-2887. DOI: 10.7498/aps.57.2882

15 陳世國, 吉世印, 劉萬松, 等. 基于小波分析的高斯脈沖成形的遞歸實現(xiàn)[J]. 物理學(xué)報, 2009,58(5): 3041-3046. DOI: 10.7498/aps.58.3041

CHEN Shiguo, JI Shiyin, LIU Wansong, et al. Recursive implementation of Gaussian pulse shaping based on wavelet analysis[J]. Acta Physica Sinica, 2009,58(5): 3041-3046. DOI: 10.7498/aps.58.3041

16 Nakhostin M. Recursive algorithm for real-time digital CR-(RC)n pulse shaping[J]. IEEE Transactions on Nuclear Science, 2011,58(5): 2378-2381. DOI: 10.1109/ TNS.2011.2164556

17 Zhou J B, Zhou W, Hong X. Improvement of digital S-K filter and its application in nuclear signal processing[J]. Nuclear Science and Techniques, 2013,24(6): 060401. DOI: 10.13538/j.1001-8042/nst.2013.06.020

Simulation study on Gaussian pulse shaping algorithm

HONG Xu1NI Shijun1ZHOU Jianbin1MA Yingjie1LIU Yi1,2ZHOU Wei1,2
1(CollegeofNuclearTechnologyand Automation Engineering,ChengduUniversity of Technology,Chengdu610059,China)
2(Fundamental Science on Radioactive Geology and Exploration Technology Laboratory,East China Institute of Technology,Nanchang 330013,China)

Background:The Gaussian pulse possesses the advantages of time response, frequency response and signal-to-noise ratio (SNR), thereby nuclear signal is usually shaped to Gaussian or Semi-Gaussian pulse shape.Purpose:This work aims to simulation study on Gaussian pulse shaping algorithm and discuss the influence of shaping parameters on the amplitude and width of the shaped pulse.Methods:The cutoff frequency and quantify factor are introduced into Gaussian pulse shaping algorithm. The filtering performance of the algorithm and its effects on amplitude spectrum are evaluated and compared with the trapezoidal pulse shaping algorithm.Results:The simulation results suggest that Gaussian pulse shaping algorithm has better performance than trapezoidal pulse shaping algorithm in filtering when the peaking times of the shaped pulses are equal. Also, Gaussian pulse shaping algorithm can reduce the full width at half maximum (FWHM) of nuclear signal amplitude spectrum.Conclusion:Combined with the trapezoidal pulse shaping algorithm, the defect of Gaussian pulse shaping algorithm in pile-up pulse separation can be improved to increase output stream of nuclear signal.

Gaussian pulse shaping, Trapezoidal pulse shaping, Nuclear signal, Full width at half maximum

HONG Xu, male, born in1989, graduated from Chengdu University of Technology with a master’s degree in 2015, doctoral student, engaged in nuclear signal acquisition and processing research

TL822+.6

10.11889/j.0253-3219.2016.hjs.39.110403

國家自然科學(xué)基金(No.11475036)、放射性地質(zhì)與勘探技術(shù)國防重點實驗室開放基金(No.RGET1201)資助

洪旭，男，1989年出生，2015年于成都理工大學(xué)獲碩士學(xué)位，現(xiàn)為博士研究生，從事核信號獲取與處理研究

Supported by National Natural Science Foundation of China (No.11475036), National Defense Key Laboratory of Radioactive Geology and

Exploration Technology, National Natural Science Foundation of China (No.RGET1201)

2016-05-09，

2016-09-19