核数据处理b_003.pdf

1、成都理工大学 马英杰核数据处理第三章第三章 谱数据处理谱数据处理-3 3成都理工大学成都理工大学 核自学院核自学院成都理工大学 马英杰谱数据处理谱数据处理谱数据处理谱数据处理定量分析定量分析定量分析定量分析?定量分析的内容定量分析的内容?1）峰边界道的确定；）峰边界道的确定；?2）峰面积计算；）峰面积计算；?3）重峰分析；）重峰分析；?4）含量计算。）含量计算。成都理工大学 马英杰峰边界道的确定峰边界道的确定峰边界道的确定峰边界道的确定?目的目的?峰位对应的单个计数小，精度不够高；而且整个高斯分布都是同一能量射线作用的，所以常用整个分布的计数和作评价计算用。峰位对应的单个计数小，精度不够高；而

2、且整个高斯分布都是同一能量射线作用的，所以常用整个分布的计数和作评价计算用。?峰边界道：整个高斯分布的起始道、终止道。起始道就是峰的左边界，终止道就是峰的右边界。峰边界道：整个高斯分布的起始道、终止道。起始道就是峰的左边界，终止道就是峰的右边界。?峰边界道的确定，直接影响峰面积的计算。峰边界道的确定，直接影响峰面积的计算。成都理工大学 马英杰峰边界道的确定峰边界道的确定峰边界道的确定峰边界道的确定?方法方法?1）根据观察谱线，直接输入左右边界道址）根据观察谱线，直接输入左右边界道址?2）各寻峰法中确定峰边界的方法来确定）各寻峰法中确定峰边界的方法来确定?简单比较法简单比较法?导数法导数法?对称

3、零面积法等等对称零面积法等等?3）用峰的全宽度确定峰边界道址）用峰的全宽度确定峰边界道址成都理工大学 马英杰峰边界道的确定峰边界道的确定峰边界道的确定峰边界道的确定?简单比较法简单比较法?峰的左右边界道i-L和i+R的确定峰的左右边界道i-L和i+R的确定?当峰形出现 左、右当峰形出现 左、右“拖尾拖尾”时，应限制时，应限制 L、R 的长度的长度半高宽半高宽ririmrililimlidatakdatadatadatakdatadata+)()(FWTMchchFWTMchchpRLp成都理工大学 马英杰峰边界道的确定峰边界道的确定峰边界道的确定峰边界道的确定?导数法导数法?1）一阶导数法）一

4、阶导数法?一阶导数由负变正一阶导数由负变正=0处处?2）二阶导数法）二阶导数法?二阶导数正极大值点二阶导数正极大值点?3）三阶导数法）三阶导数法?三阶导数由正变负三阶导数由正变负=0处处成都理工大学 马英杰峰边界道的确定峰边界道的确定峰边界道的确定峰边界道的确定?对称零面积变换法对称零面积变换法?找峰的准则：找峰的准则：?为了使找峰的灵敏度与实验谱的统计性无关，一般用变换谱的标准偏差为单位找峰，即当变换谱与它的标准偏差之比出现正极值，且此极值超过一给定的数值为了使找峰的灵敏度与实验谱的统计性无关，一般用变换谱的标准偏差为单位找峰，即当变换谱与它的标准偏差之比出现正极值，且此极值超过一给定的数值

5、f时，就认为是峰。时，就认为是峰。?峰边界的确定：峰边界的确定：SSi的正峰两边相邻的负峰址。的正峰两边相邻的负峰址。fdataCdataCyySSmmjjijmmjjijiii=+=+212成都理工大学 马英杰峰边界道的确定峰边界道的确定峰边界道的确定峰边界道的确定?用峰的全宽度确定峰边界道址用峰的全宽度确定峰边界道址?全宽度：全宽度：1/10峰高处对应的道址峰高处对应的道址?方法：方法：?1.选择不受干扰的单峰；选择不受干扰的单峰；?2.分别在峰左、峰右找峰底：分别在峰左、峰右找峰底：?1)确定峰高h，h=datap-datab;?2)计算的值；?3)确定处的ch值chL、chR。RLhh

6、)(101)(101RLhh)(101)(101成都理工大学 马英杰峰面积的计算峰面积的计算峰面积的计算峰面积的计算?意义意义?峰面积的计算是定量分析的基础。峰面积的计算是定量分析的基础。?知道了特征峰的净峰面积，就可以计算目标元素的含量。知道了特征峰的净峰面积，就可以计算目标元素的含量。?实测谱中，各特征峰是叠加在环境本底和康普顿散射背景之上的。实测谱中，各特征峰是叠加在环境本底和康普顿散射背景之上的。?总面积总面积S：在一个指定的峰区内，各道计数之和：在一个指定的峰区内，各道计数之和?本底面积本底面积B：由环境本底和散射造成的计数总和：由环境本底和散射造成的计数总和?净峰面积净峰面积A：由

7、峰的总面积扣除本底面积即可得出净峰面积：由峰的总面积扣除本底面积即可得出净峰面积即：峰的总面积即：峰的总面积本底面积本底面积=净峰面积净峰面积S-B =A，所以，计算净峰面积，如何确定，所以，计算净峰面积，如何确定B最关键！最关键！成都理工大学 马英杰峰面积的计算峰面积的计算峰面积的计算峰面积的计算?A=S B，关键是如何确定，关键是如何确定B？?方法：方法：?线性本底法（总峰面积法，线性本底法（总峰面积法，TPA法）法）?Covell（科沃尔）峰面积法（科沃尔）峰面积法?Wasson（瓦森、沃森）峰面积法（瓦森、沃森）峰面积法?Sterlinski(斯托林斯基斯托林斯基)峰面积法峰面积法?平

8、均总峰面积法平均总峰面积法?单峰曲线拟合法单峰曲线拟合法成都理工大学 马英杰LRyLyR峰面积的计算峰面积的计算峰面积的计算峰面积的计算?线性本底法（总峰面积法，线性本底法（总峰面积法，TPA法）法）?确定本底面积计算方法：确定本底面积计算方法：?左右边界点直线连接即为本底线左右边界点直线连接即为本底线?线性本底线性本底?梯形法计算本底面积梯形法计算本底面积2)1(*)(2)1(*)(+=+=LRyyLRdatadataBRLRLB成都理工大学 马英杰峰面积的计算峰面积的计算峰面积的计算峰面积的计算?线性本底法（总峰面积法，线性本底法（总峰面积法，TPA法）法）?1）确定峰的左、右边界）确定峰

9、的左、右边界L、R?2）计算总面积：）计算总面积：?3）计算本底面积：）计算本底面积：?4）计算净峰面积：）计算净峰面积：=RLiidataSBSA=2)1(*)(+=LRdatadataBRLLBRyLyR成都理工大学 马英杰峰面积的计算峰面积的计算峰面积的计算峰面积的计算?线性本底法（总峰面积法，线性本底法（总峰面积法，TPA法）法）?编程：编程：=RLiidataSBSA=2)1(*)(+=LRdatadataBRLint CMmcaView:CalculatePeakArea(int L,int R)int bdarea,zarea=0,area;bdarea=(DataL+DataR

10、)*(R-L+1)/2;/计算总面积计算总面积for(int i=L;i=R;i+)zarea=Datai+zarea;/计算净峰面积计算净峰面积area=zarea-bdarea;return area;成都理工大学 马英杰峰面积的计算峰面积的计算峰面积的计算峰面积的计算?线性本底法（总峰面积法，线性本底法（总峰面积法，TPA法）法）?编程：编程：=RLiidataSBSA=2)1(*)(+=LRdatadataBRLint CMmcaView:CalculatePeakArea(int L,int R,)int bdarea,zarea=0,area;bdarea=(DataL+DataR

11、)*(R-L+1)/2;/计算总面积计算总面积for(int i=L;i=R;i+)zarea=Datai+zarea;/计算净峰面积计算净峰面积area=zarea-bdarea;return area;计算本底面积计算本底面积成都理工大学 马英杰峰面积的计算峰面积的计算峰面积的计算峰面积的计算?线性本底法（总峰面积法，线性本底法（总峰面积法，TPA法）法）?编程：编程：=RLiidataSBSA=2)1(*)(+=LRdatadataBRLint CMmcaView:CalculatePeakArea(int L,int R,)int bdarea,zarea=0,area;bdarea=

12、(DataL+DataR)*(R-L+1)/2;/计算总面积计算总面积for(int i=L;i=R;i+)zarea=Datai+zarea;/计算净峰面积计算净峰面积area=zarea-bdarea;return area;计算总峰面积计算总峰面积成都理工大学 马英杰峰面积的计算峰面积的计算峰面积的计算峰面积的计算?线性本底法（总峰面积法，线性本底法（总峰面积法，TPA法）法）?编程：编程：=RLiidataSBSA=2)1(*)(+=LRdatadataBRLint CMmcaView:CalculatePeakArea(int L,int R,)int bdarea,zarea=0,

13、area;bdarea=(DataL+DataR)*(R-L+1)/2;/计算总面积计算总面积for(int i=L;in时，用最小二乘拟合法求解各系数。时，用最小二乘拟合法求解各系数。a0，a1，a2，an，成都理工大学 马英杰含量计算含量计算含量计算含量计算?多标法多标法最小二乘拟合法最小二乘拟合法(工作曲线法工作曲线法)?计算含量：计算含量：?工作曲线为工作曲线为n次多项式：（常取次多项式：（常取n=2）wx=a0+a1Ax+a2Ax2+anAxn。?在相同条件下测量样品，得对应峰的面积计数在相同条件下测量样品，得对应峰的面积计数Ax，?带入上多项式方程即可得含量值。带入上多项式方程即可得含量值。