亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        n-γ核信號的matlab仿真及其數(shù)字波形甄別方法研究

        2016-12-25 08:53:52曹真?zhèn)?/span>顏擁軍付德順王慶震易凌帆
        核科學(xué)與工程 2016年2期
        關(guān)鍵詞:信號方法

        曹真?zhèn)?顏擁軍,付德順,王慶震,易凌帆

        (南華大學(xué)核科學(xué)技術(shù)學(xué)院,湖南 衡陽421001)

        n-γ核信號的matlab仿真及其數(shù)字波形甄別方法研究

        曹真?zhèn)?顏擁軍,付德順,王慶震,易凌帆

        (南華大學(xué)核科學(xué)技術(shù)學(xué)院,湖南 衡陽421001)

        由于中子射線中常?;煊笑蒙渚€,所以在探測器探測到中子信號的同時也可能會探測到γ信號,而對不同的核信號進行區(qū)分具有重大意義。本文利用matlab軟件開展了中子和γ粒子數(shù)字化核信號波形仿真研究,并利用幾種核脈沖波形甄別方法對仿真波形進行甄別,最后對不同方法的甄別效果展開了分析討論。結(jié)果表明,電荷比較法、上升時間法和脈沖梯度分析法都是穩(wěn)定的甄別方法,可有效實現(xiàn)n-γ甄別,其中電荷比較法的甄別效果相對較好。

        matlab;數(shù)字核信號;仿真;波形甄別

        基于模擬電路的核脈沖波形甄別發(fā)展至今已比較成熟,常用的方法有電荷比較法、上升時間法等等[1-2]。其中,電荷比較法是一種較為經(jīng)典的方法,常常用于n-γ、α-β(γ)甄別等;上升時間法通過對信號進行積分起到了平滑處理信號的作用,可以得到一定的甄別效果。然而隨著技術(shù)要求的提高,模擬電子電路過于復(fù)雜,處理能力有限和不穩(wěn)定性等弊端越發(fā)地明顯,這正好促進了數(shù)字化核脈沖波形甄別的發(fā)展,特別是在以高速模數(shù)轉(zhuǎn)換器(ADC)為代表的數(shù)字器件得到廣泛應(yīng)用之后,數(shù)字核脈沖波形甄別的研究取得了較大進步[3]。

        中子射線中往往混有γ射線,這兩類射線射入閃爍體時會產(chǎn)生大量的次級帶電粒子,中子和γ粒子的質(zhì)量阻止本領(lǐng)不同,從而能量沉積密度也就不同,于是中子和γ粒子的從光電倍增管中出來的電流信號脈沖在形狀上就會有差異,這就是n-γ波形甄別的依據(jù)。

        本研究根據(jù)核脈沖信號的特征,利用matlab軟件編寫了n-γ數(shù)字核脈沖波形數(shù)值仿真程序,并將仿真的核脈沖分別用電荷比較法、上升時間法和脈沖梯度分析法進行波形甄別,最后比較了這三種方法的甄別效果。

        1 核信號的特征和仿真

        1.1 核信號數(shù)學(xué)模型

        閃爍體探測器探測效率高,對射線靈敏,時間分辨好,環(huán)境適應(yīng)性強,有著良好的n-γ分辨能力。在有機閃爍體中,中子和γ粒子激發(fā)的熒光脈沖含有快、慢兩種成分。實驗證明,有機閃爍體內(nèi)中子和γ粒子所產(chǎn)生的熒光脈沖可表示為:

        其中:τf——快成分的衰減時間;

        τs——慢成分的衰減時間;

        If——快成分強度;

        Is——慢成分強度。

        往往τs比τf大一個量級,而If和Is的相對大小與入射粒子類型有關(guān),一般而言,由中子引起的次級帶電粒子能量沉積密度較γ粒子的大,這就使得在退激發(fā)光過程中,中子的Is不小于If,而γ粒子的情況則相反[4]。因此,中子與γ粒子在閃爍體中產(chǎn)生的熒光脈沖的形狀就存在差異,如圖1所示。

        當光電倍增管工作在線性范圍時,從它的

        圖1 中子與γ粒子對應(yīng)的熒光脈沖形狀Fig.1 The Fluorescence Pulse of Neutron and Gamma Particle

        最后一級倍增極引出的電流脈沖形狀應(yīng)當反映閃爍體發(fā)射的光脈沖形狀,近似為[4]:

        式中:if——電流脈沖中與熒光脈沖相對應(yīng)的

        快成分最大值;

        is——電流脈沖中與熒光脈沖相對應(yīng)的

        慢成分最大值。

        本研究以(2)式作為仿真核信號數(shù)學(xué)模型。

        1.2 核信號統(tǒng)計特性

        (1)脈沖時間間隔E(t)滿足指數(shù)分布規(guī)律:

        其中:C——脈沖的平均計數(shù)率。

        在本研究對仿真核脈沖信號的甄別處理中,由于脈沖時間間隔對不會影響到甄別效果,故為簡化程序,仿真核脈沖以固定時間間隔出現(xiàn)。

        (2)脈沖幅度分布H(A)滿足正態(tài)分布:

        式中:A——平均幅度;

        σA——幅度標準差,與探測器固有能量分辨率Rd有如下關(guān)系:

        其中:E0——能譜峰中心值;

        FWHM ——峰半高寬。

        (3)脈沖白噪聲分布h(a)滿足正態(tài)分布規(guī)律:

        式中:a——噪聲平均幅度;

        σa——噪聲幅度標準差。

        1.3 核脈沖信號仿真

        本研究的數(shù)字核脈沖信號仿真,為了對比電荷比較法、上升時間法和脈沖梯度分析法的n-γ甄別性能,設(shè)置n、γ核脈沖信號的幅度均值、幅度標準偏差、噪聲標準偏差相同,而其他波形參數(shù)見表1[5],仿真出的一組幅度不等的濾波之前和濾波之后的n、γ波形圖例如圖2、圖3所示,仿真程序流程見圖4。

        圖2 濾波前n-γ仿真波形Fig.2 n-γSimulation Pulse Shape before Filtering

        圖3 濾波后n-γ仿真波形Fig.3 n-γSimulation Pulse Shapeafter Filtering

        表1 中子射線和γ射線部分波形參數(shù)Table 1 Partial Parameter of Neutron and GammaPulse Shape

        2 數(shù)字波形甄別方法

        2.1 電荷比較法

        通過上述可知,光電倍增管出來的電流脈沖也包含快慢成分,那么在閃爍體探測器等效

        圖4 仿真程序流程圖Fig.4 Simulation Program Flow Chart

        輸出回路的電容上積累的電荷也由快、慢成分組成:

        其中:Q f——電荷的快成分;

        Qs——電荷的慢成分。

        由于中子和γ粒子所產(chǎn)生的核脈沖波形不同,就導(dǎo)致了中子和γ粒子的Q f/Q、Qs/Q、Q s/Q f等比值有較大差異。因此可以通過比較不同帶電粒子的這些比值中的任何一個值來鑒別粒子。

        本研究選擇比值Q f/Q來作為中子和γ粒子的甄別依據(jù),對于快成分電荷,采用從電流波形的起始到波形下降至峰值的25%處這一段波形的積分作為快成分電荷Q f;總電荷積分Q則采用全電流波形積分。

        具體的程序原理:在matlab中建立一個單元型變量,將濾波后的各個波形作為該單元型變量的各個元素存儲在該單元型變量中,再對每一個波形進行尋峰,求出對應(yīng)的峰值時刻和峰值大小,再根據(jù)峰值求出波形下降至峰值25%時的對應(yīng)時刻,將此時刻設(shè)為t1,波形起始時刻作為t0,波形基本衰減完全的對應(yīng)時刻為t2。t0-t1時間對應(yīng)的電流脈沖積分作為波形的快成分積分Q f,而t0-t2時間對應(yīng)的電流脈沖積分作為波形的總積分Q,選擇比值Q f/Q作為中子和γ粒子的電荷比較法波形甄別因子。

        波形甄別效果的好壞往往用品質(zhì)因子FOM來衡量:

        式中 Δt——甄別譜上γ粒子峰位和中子峰位之間的距離;

        FWHMn——甄別譜上中子峰的半高寬;

        FWHMγ——甄別譜上γ粒子峰的半高寬。

        FOM值越大,則甄別效果越好。

        本次仿真分別產(chǎn)生中子波形和γ粒子波形各2 500個,運行程序所得到的仿真波形的甄別因子Q f/Q散點圖如圖5所示,其中橫軸為仿真波形的序號,縱軸為Q f/Q數(shù)值,甄別譜如圖6所示,得到甄別品質(zhì)因子FOM=1.35。

        圖5 電荷比較法甄別因子散點圖Fig.5 Scatter Chart of Discrimination Factor for Charge Comparison Method

        圖6 電荷比較法甄別譜Fig.6 Discrimination Spectrum of Charge Comparison Method

        2.2 上升時間法

        圖1 所示的核信號波形是有個先上升后下降,最后下降衰減至零的過程。這就意味著如果對這個波形求積分,那么其積分值F(t)最后應(yīng)該會接近一個常值,設(shè)為F。設(shè)常數(shù)a、b滿足:a<b,且a∈[0,1],b∈[0,1],在積分值增加至F的過程中存在一小一大的兩個不同時刻T a、T b,在這兩個時刻時分別有F(T a)=aF,F(T b)=bF,定義ΔT=T b-T a為上升時間。

        分別對中子和γ粒子的核信號波形進行積分。對于中子核信號積分,設(shè)a F n增加至b F n所對應(yīng)的時間ΔTn為其上升時間;對于γ粒子核信號積分,再設(shè)a Fγ增加至b Fγ所對應(yīng)的時間ΔTγ為其上升時間,由于γ粒子核信號衰減較中子核信號快,所以ΔTγ要比ΔTn小[6]。

        在仿真中取a=0.15,b=0.85,分別仿真中子和γ粒子核信號各2 500個,上升時間散點圖和甄別譜分別如圖7、8所示,得到甄別品質(zhì)因子FOM=1.20。

        圖7 上升時間法甄別因子散點圖Fig.7 Scatter Chart of Discrimination Factor for Risetime Method

        圖8 上升時間法甄別譜Fig.8 Discrimination Spectrum of Risetime Method

        2.3 脈沖梯度分析法

        脈沖梯度分析法又稱為PGA(Pulse Gradient Analysis)方法,是隨著數(shù)字技術(shù)的發(fā)展而開發(fā)出的一種數(shù)字波形甄別方法[6]。其原理在于:中子核信號比γ粒子核信號衰減得更慢,這也表明中子核信號下降部分的梯度比γ粒子核信號的更小,通過不同核脈沖梯度的差異便可實現(xiàn)波形甄別[7]。

        在仿真中,設(shè)脈沖梯度估計值:

        其中:yp——脈沖峰值;

        ys——另一個采樣點的幅值,在波形下降沿中選取;

        ΔT——兩個采樣點的時間間隔。

        分別仿真中子和γ粒子核信號各2 500個,ys在峰值大小15%的對應(yīng)點上取,G值散點圖和甄別譜分別如圖9、10所示,得到甄別品質(zhì)因子FOM=0.99。

        圖9 脈沖梯度分析法散點圖Fig.9 Scatter Chart of PGA Method

        圖10 脈沖梯度分析法甄別譜Fig.10 Discrimination Spectrum of PGA Method

        2.4 三種方法的甄別效果對比

        先后分別仿真中子和γ波形目各1 500個、2 000個、2 500個,求出這些波形目分別在上述三種方法下的甄別品質(zhì)因子,如圖11所示。

        圖11 三種方法甄別效果比較圖Fig.11 Dicrimination Effect Comparison Chart of the Three Different Methods

        由圖11可以看出,隨著仿真波形數(shù)目的增多,各個方法的FOM值近似線性地增大,但增幅很小,并沒有因為波形的增多而產(chǎn)生較為劇烈的波動,說明在其他條件不變的情況下,仿真波形的數(shù)目不會對三種甄別方法的效果產(chǎn)生較大影響,這些方法都較為穩(wěn)定。在甄別效果上電荷比較法較好,上升時間法次之,脈沖梯度分析法較差,這是因為脈沖梯度分析法在原理上只利用了核信號波形中的兩個有效的采樣點,對完整的波形的利用率不高,而電荷比較法比起另外兩種方法更多地利用了全部脈沖所包含的信息,因此FOM更高,這一結(jié)論是與相關(guān)文獻[9]的結(jié)論相似的。

        3 結(jié)論

        本研究基于matlab軟件先根據(jù)核信號波形的統(tǒng)計特性進行了閃爍體探測器的n-γ核脈沖信號波形仿真,再利用電荷比較法、上升時間法和脈沖梯度分析法對仿真信號開展甄別研究。這三種波形甄別方法看似不同,但實際上都是尋找一個能將不同粒子的核信號在波形衰減上的差異表征出來的量作為甄別依據(jù)進行粒子鑒別。三種甄別方法的結(jié)果表示,三種方法都是穩(wěn)定的甄別方法,可有效實現(xiàn)n-γ甄別,其中電荷比較法由于更多地利用了波形信號而得到相對較好的甄別效果。

        [1] Brooks.Nucl.Instr.and Method,1959,4:151.

        [2] Roush M L,Wilson M A,Hornyak W F.Nucl.Instr.a(chǎn)nd Method,1964,31:112.

        [3] Simoes J.Basilio,Correia Carlos M.B.A,Pulse processing architectures,Nucl.Instr.a(chǎn)nd Method 1972,99:525-539.

        [4] 凌球,郭蘭英.核輻射探測[M].北京:原子能出版社,2002.

        [5] 左廣霞,何彬,等.n-γ脈沖信號的仿真及其波形甄別技術(shù)研究[J].核技術(shù),2012,35(5):365-368.

        [6] Mellow B D,et a1.Digital discrimination of neutrons and γ-rays in liquid scintillators using pulse gradient analysis[J].Nuclear Instruments and Methods A,2007,578(1):191-197.

        [7] S.D.Jastaniah,P.J.Sellin.Digital techniques for n/γ pulse shape discrimination and capture-gated neutron spectroscopy using liquid scintillators[J].NIM A,2004,517:202-210.

        [8] 陳小輝.基于氘化液閃探測器的數(shù)字化γ射線探測研究[D].綿陽:中國工程物理研究院,2013:15-16.

        [9] 黃林.基于數(shù)字化技術(shù)的芪晶體探測器脈沖形狀甄別方法研究[J].核電子學(xué)與探測技術(shù),2012,32(9):993-997.

        The Research on Matlab Simulation and Digital Pulse Shape Discrimination for n-γNuclear Signal

        CAO Zhen-wei,YAN Yong-jun,FU De-shun,WANG Qing-zhen,YI Ling-fan
        (University of South China,Hengyang of Hunan Prov.421001,China)

        Neutron radiation usually mingle withγradiation,as a result of it,the neutron signal of radiation detector usually mixγsignal.It is significant to discriminate different nuclear signal.In this study,we use matlab to simulate n-γnuclear pulse shape digital signal and to discriminate them by pulse shape discrimination methods,then,we discuss the discrimination effect of different method.The result shows that charge comparison method,rising time method and pulse gradient analysis are all stable and reliable for n-γdiscrimination,and charge comparison method is relatively better.

        Matlab;Digital Nuclear Signal;Simulation;Pulse Shape Discrinimation

        TL812

        A

        0258-0918(2016)01-0293-06

        2015-07-27

        湖南省研究生科研創(chuàng)新項目(2015SCX23)

        曹真?zhèn)?1989—),男,湖南郴州人,南華大學(xué)核科學(xué)技術(shù)學(xué)院碩士研究生,研究方向為核測控

        顏擁軍:yan_jason@163.com

        猜你喜歡
        信號方法
        信號
        鴨綠江(2021年35期)2021-04-19 12:24:18
        完形填空二則
        學(xué)習(xí)方法
        孩子停止長個的信號
        可能是方法不對
        用對方法才能瘦
        Coco薇(2016年2期)2016-03-22 02:42:52
        基于LabVIEW的力加載信號采集與PID控制
        一種基于極大似然估計的信號盲抽取算法
        四大方法 教你不再“坐以待病”!
        Coco薇(2015年1期)2015-08-13 02:47:34
        賺錢方法
        激情一区二区三区视频| 国产啪精品视频网站| 九九在线视频| 亚洲精彩视频一区二区| 精品国产一区二区三区av麻| 色欲aⅴ亚洲情无码av| 日韩成人无码一区二区三区| 欧美日韩亚洲国产无线码| 日本淫片一区二区三区| 中国妇女做爰视频| a级毛片毛片免费观看久潮喷| 在线国产视频精品视频| 日本人妻高清免费v片| 国产高清在线观看av片| www国产精品内射熟女| 国产丝袜免费精品一区二区 | 国产激情一区二区三区| 色悠久久久久综合欧美99| 男女高潮免费观看无遮挡| 国产精品一区二区三区在线观看 | 一区二区韩国福利网站| 国产精品又污又爽又色的网站| 久久久99精品免费视频| 图片区小说区激情区偷拍区| 91人妻无码成人精品一区91| 国产一区二区三区不卡在线播放| 18禁在线永久免费观看| 在线va免费看成| 毛片一级精油按摩无码| 国内自拍视频一区二区三区| 香港台湾经典三级a视频| 日本午夜福利| 国产免费一区二区三区三| 亚洲av无码日韩av无码网站冲| 亚洲不卡中文字幕无码| 日韩精品精品一区二区三区| 成人影院视频在线免费观看| 精品水蜜桃久久久久久久| 一区二区三区不卡在线| 国内偷拍精品一区二区| 亚洲成a人无码|