发明名称 基于支持向量机的液体闪烁体探测器n/γ射线甄别方法
摘要 本发明提供了一种基于支持向量机的液体闪烁体探测器n/γ射线甄别方法,该方法的步骤如下:(1)信号采集,采用液体闪烁体探测器采集由中子源产生的信号,并将所采集的信号转化为脉冲信号;(2)采用电荷比较法甄别中子和γ射线,采用传统电荷比较法对步骤(1)获取的脉冲信号进行甄别,确定电荷比较法对中子和γ射线的甄别能量下限;(3)获取M矩阵;(4)采用LIBSVM软件进行中子和γ射线的分类预测。该方法可以对能量处于传统方法甄别阈值以下的信号进行区分,从而有助于提高中子通量测量的准确性。
申请公布号 CN104536032A 申请公布日期 2015.04.22
申请号 CN201410812137.7 申请日期 2014.12.23
申请人 四川大学 发明人 幸浩洋;余训臻;朱敬军;林兴德;王力;唐昌建;夏玉玺;张才勋;杜强;魏伟伟
分类号 G01T1/204(2006.01)I;G06F19/00(2011.01)I 主分类号 G01T1/204(2006.01)I
代理机构 成都科海专利事务有限责任公司 51202 代理人 刘双兰;郭萍
主权项 基于支持向量机的液体闪烁体探测器n/γ射线甄别方法,其特征在于步骤如下:(1)信号采集采用液体闪烁体探测器采集由中子源产生的信号,并将所采集的信号转化为脉冲信号;(2)采用电荷比较法甄别中子和γ射线采用传统电荷比较法对步骤(1)获取的脉冲信号进行甄别,确定电荷比较法对中子和γ射线的甄别能量下限;(3)获取M矩阵①确定训练集和预测集将能量高于步骤(2)所得甄别能量下限的脉冲信号作为训练点,所有训练点组成训练集,将能量等于或低于步骤(2)所得甄别能量下限的脉冲信号作为预测点,所有预测点组成预测集;②确定M矩阵的特征向量分别扣除步骤(1)采集的各脉冲信号的基线,将扣除基线后的每个脉冲信号分别用M矩阵表示,也就是将各训练点和预测点用M矩阵表示,M=[m<sub>0</sub>,m<sub>1</sub>,…,m<sub>k‑1</sub>,m<sub>k</sub>],其中,m<sub>0</sub>,m<sub>1</sub>,…,m<sub>k‑1</sub>,m<sub>k</sub>分别为M矩阵的特征向量,所述M矩阵的特征向量按如下方法确定:扣除基线后,每个脉冲信号第i个采样点的幅度为g[i],脉冲信号起始点j<sub>start</sub>至脉冲信号第i个采样点的积分值表示为f[i]:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>f</mi><mo>[</mo><mi>i</mi><mo>]</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><msub><mi>j</mi><mi>start</mi></msub></mrow><mi>i</mi></munderover><mi>g</mi><mo>[</mo><mi>i</mi><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000641639010000011.GIF" wi="1101" he="162" /></maths>式(1)中,j<sub>start</sub>为脉冲信号的起始点;M矩阵的特征向量表示如下:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>m</mi><mn>0</mn></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>f</mi><mo>[</mo><mi>i</mi><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000641639010000012.GIF" wi="1102" he="154" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>m</mi><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfrac><mrow><mi>f</mi><mo>[</mo><mi>i</mi><mo>]</mo><mo>&CenterDot;</mo><mi>i</mi></mrow><msub><mi>m</mi><mn>0</mn></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000641639010000013.GIF" wi="1113" he="153" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>m</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><msup><mrow><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfrac><mrow><mi>f</mi><mo>[</mo><mi>i</mi><mo>]</mo><mo>&CenterDot;</mo><msup><mi>i</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msup></mrow><msub><mi>m</mi><mn>0</mn></msub></mfrac><mo>|</mo></mrow><mrow><mn>1</mn><mo>/</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000641639010000014.GIF" wi="1130" he="184" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>m</mi><mi>k</mi></msub><mo>=</mo><msup><mrow><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mfrac><mrow><mi>f</mi><mo>[</mo><mi>i</mi><mo>]</mo><mo>&CenterDot;</mo><msup><mi>i</mi><mi>k</mi></msup></mrow><msub><mi>m</mi><mn>0</mn></msub></mfrac><mo>|</mo></mrow><mrow><mn>1</mn><mo>/</mo><mi>k</mi></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000641639010000015.GIF" wi="1120" he="198" /></maths>式(2)~(5)中,N为整个脉冲信号所包含的采样点的数量,i为脉冲信号第i个采样点,k为≥5的正整数;(4)采用LIBSVM软件进行中子和γ射线的分类预测①分别将步骤(3)得到的各M矩阵进行归一化,使||M||=1;②将归一化后的M矩阵表示的各训练点作为LIBSVM软件的输入信号,使用RBF函数作为核函数对所有训练点组成的训练集进行超平面分类模型训练,得到超平面分类模型;③将归一化后的M矩阵表示的各预测点和超平面分类模型输入LIBSVM软件中,对预测集中的预测点进行分类预测,即完成液体闪烁体探测器n/γ射线甄别。
地址 610065 四川省成都市武侯区一环路南一段24号