[发明专利]一种γ射线能谱中重峰区域全能峰基底的计算方法及系统

权利要求书

1.一种γ射线能谱中重峰区域全能峰基底的计算方法，其特征在于，包括：

S1、对γ射线能谱中重峰区域的每个全能峰，将每个全能峰的净计数部分表示为高斯函数，将每个全能峰的的基底部分表示为余误差函数，将每个全能峰表示为对应的高斯函数与对应的余误差函数之和；

S2、对所述重峰区域的能谱进行数据拟合，计算得到每个高斯函数和每个余误差函数中的各待定参数值；

S3、根据每个高斯函数和每个余误差函数中的各待定参数值，计算得到所述重峰区域的各道基底计数。

2.根据权利要求1所述的一种γ射线能谱中重峰区域全能峰基底的计算方法，其特征在于，步骤S1包括：

如果所述重峰区域共有N个全能峰，则第i个全能峰表示为：

N个全能峰表示为：

其中，表示第i个全能峰的净计数部分，L_ierfc(w_i(E-E_i))表示第i个全能峰的基底部分；H_i为高斯函数幅值，L_i为余误差函数幅值，w_i为高斯函数和余误差函数的宽度参数，E_i为高斯函数和余误差函数的中心位置，E是γ射线能谱上的数据点的横坐标；H_i、L_i、w_i、E_i(i＝1,…,N)为待定参数；erfc为余误差函数；i、N为正整数。

3.根据权利要求2所述的一种γ射线能谱中重峰区域全能峰基底的计算方法，其特征在于，步骤S2包括：

基于式(2)，对所述重峰区域的能谱进行数据拟合，即求下式的最小值：

其中，y_j为第j道实测计数，f(E_j)为第j个全能峰，M为所述重峰区域的总道数，即总数据点个数,j、M为正整数。

4.根据权利要求3所述的一种γ射线能谱中重峰区域全能峰基底的计算方法，其特征在于，通过迭代算法求解式(3)的最小值，得到各待定参数H_i、L_i、w_i、E_i(i＝1,…,N)的值。

5.根据权利要求4所述的一种γ射线能谱中重峰区域全能峰基底的计算方法，其特征在于，步骤S3包括：

基于计算得到的H_i、L_i、w_i、E_i(i＝1,…,N)，通过下式计算得到所述重峰区域的各道基底计数B(E_j)为：

其中，余误差函数erfc表示为：

其中，x为自变量，在式(1)、(2)和(4)中x＝w_i(E-E_i)；t为积分变量。

6.一种γ射线能谱中重峰区域全能峰基底的计算系统，其特征在于，包括：

表示模块，用于对γ射线能谱中重峰区域的每个全能峰，将每个全能峰的净计数部分表示为高斯函数，将每个全能峰的的基底部分表示为余误差函数，将每个全能峰表示为对应的高斯函数与对应的余误差函数之和；

数据拟合模块，用于对所述重峰区域的能谱进行数据拟合，计算得到每个高斯函数和每个余误差函数中的各待定参数值；

计算模块，用于根据每个高斯函数和每个余误差函数中的各待定参数值，计算得到所述重峰区域的各道基底计数。