文献类型：专利

专利类型：发明专利

是否失效：否

是否授权：否

申 请 号：CN201610101592.5

申 请 日：20160225

发 明 人：艾宪芸 肖无云 李京伦 马新华 王广柱 史志兰 梁卫平 陈小猛 张羽中 张磊 张斌

申 请 人：中国人民解放军63973部队

申请人地址：102205 北京市1044信箱200号

公 开 日：20190212

公 开 号：CN105785424B

代 理 人：刘永盛;肇文丽

代理机构：11046 中国人民解放军防化研究院专利服务中心

语　　种：中文

摘　　要：本发明涉及一种碲锌镉探测器伽玛谱全能峰非线性拟合算法。该算法提供了碲锌镉探测器伽玛能谱全能峰拟合函数，在拟合函数的非线性拟合过程中，借鉴工业上物理退火的物理过程和Metropolis准则，结合最小二乘Levenberg‑Marquardt(LM)迭代算法，采用内、外两个循环，进行非线性拟合，既提高了迭代计算的速度，又达到全局收敛的目的，克服了迭代结果对拟合参数初始值过度依赖的缺点。

主 权 项：1.一种碲锌镉CZT探测器伽玛全能峰非线性拟合算法，其特征在于该拟合算法具体步骤如下：第一步：确定CZT探测器伽玛谱全能峰的拟合函数式为：f(i，p)＝G(i)+B(i)+S(i)+D(i)式中：G(i)＝H_g×exp[-(i-i₀)²/(2σ²)]；B(i)＝A₅+A₆×i；S(i)＝H_s×H_g×erfc[(i-i₀)/(σ2^1/2)]；D(i)＝H_t×H_g×exp[(i-i₀)/(T_sσ)]×erfc[(i-i₀)/(σ2^1/2)+1/(T_s2^1/2)]；i_o＝A₁+A₂E_peak；σ²＝A₃+A₄E_peak； 式中：H_s为本底台阶高度；H_g为高斯峰的峰高；i_o为高斯峰峰位道址，i_o与E_peak线性相关，E_peak是伽玛谱全能峰峰位道址的测量值；σ²为高斯峰的方差，σ²与E_peak线性相关；H_t为尾巴高度；T_s为指数尾巴的反斜率；i为道址；第二步：确定需要拟合的参数共10个：H_s、H_g、A₁、A₂、A₃、A₄、A₅、A₆、H_t、T_s；第三步：设定CZT探测器伽玛谱全能峰的测量谱为y(i)，i表示道址；第四步：设定每个拟合参数的取值范围，其极大、极小值分别存入数组AH和AL中；第五步：温度初始化，对每个拟合参数在其取值范围内均匀抽样得到一组拟合参数的初始值，该初始参数向量记作p₁，把p₁带入拟合函数得到f(i，p₁)，计算残差平方和：