计算物理 ›› 2017, Vol. 34 ›› Issue (5): 603-610.

• 研究论文 • 上一篇    下一篇

中子深度分析反演算法研究

杨鑫, 李润东, 王冠博, 窦海峰, 周刚   

  1. 核物理与化学研究所, 四川 绵阳 621900
  • 收稿日期:2016-06-21 修回日期:2016-09-05 出版日期:2017-09-25 发布日期:2017-09-25
  • 作者简介:杨鑫(1988-),男,博士,助理研究员,从事反应堆物理及应用技术研究,E-mail:yangx05@126.com
  • 基金资助:
    国家自然科学基金(11475152)资助项目

Study on Unfolding Algorithms for Neutron Depth Profiling

YANG Xin, LI Rundong, WANG Guanbo, DOU Haifeng, ZHOU Gang   

  1. Institution of Nuclear Physics and Chemistry, Mianyang Sichuan 621900, China
  • Received:2016-06-21 Revised:2016-09-05 Online:2017-09-25 Published:2017-09-25

摘要: 针对中子深度分析技术,研究几类反演算法:概率迭代法、奇异值分解求解最小二乘法、线性正则化方法及约束线性正则化法.开展相关数值实验,分别在欠定方程、超定方程情况下,对反演算法的结果进行比较.一般情况下几种算法都能得到较为理想的结果,概率迭代法和约束线性正则化法由于有迭代过程,因此在离子源强度发生阶跃处得到的结果并不理想.在超定方程情况下,当选取的反演范围较为随意时,线性正则化方法可能使得解的稳定性较差,而约束线性正则化法能够很好地抑制求解过程中的不稳定性.对实际测量能谱进行反演时,由于能谱有随机误差,线性正则化方法不能很好地抑制误差的影响,反演结果振荡较强,其余方法的结果与参考值符合很好.

关键词: 中子深度分析, 反演算法, 欠定方程, 超定方程

Abstract: Unfolding algorithms for neutron depth profiling:probability iteration, SVD solving least square, linear regularization(LR) and constrained linear reqularization(CLR) are studied. All algorithms are applied for both under-determined and over-determined equations, and results are compared and discussed. Due to iterative processes, probability iteration and CLR could not work well in the case that sources intensities change sharply. LR is unstable as unfolding range is chosen arbitrarily, which can be constrained by CLR. A practical spectrum of NDP experiment is unfolded by algorithms. LR could not restrain stochastic errors caused by statistical fluctuation, while other algorithms show good unfolding results and agree well with references.

Key words: neutron depth profiling, unfolding algorithms, under-determined equations, over-determined equations

中图分类号: