4 谱数据的定量分析 核数据处理课件.ppt

东华理工大学 张怀强 峰面积的计算 Sterlinski（斯托林斯基）峰面积法 从峰面积计算结果的标准偏差来考虑，科沃尔法峰面积的方差为： 边界道计数data[i0-n], data[i0+n]的权重为(n-1/2)2,对方差的贡献大很多倍。这样，在统计涨落较大的弱峰中，边界道计数的统计涨落会严重影响峰面积计算的精度。 怎么消除这种计数统计上的边叶效应？？？？ 东华理工大学 张怀强 峰面积的计算 Sterlinski（斯托林斯基）峰面积法 基本思想 取几个连续递次增加的n (n=k,k+1,k+2,…,k+l)值。 分别求出各个n值的covell面积。 然后，将此l+1个科沃尔面积再求和以表示“峰面积”。 … 东华理工大学 张怀强 … 峰面积的计算 Sterlinski（斯托林斯基）峰面积法——峰面积计算 分别计算n取不同值的Covell峰面积 …………………………………………………………… 把l+1个面积求和： 东华理工大学 张怀强 峰面积的计算 Sterlinski（斯托林斯基）峰面积法 ——峰面积计算 令n= l+1, 则： 面积统计均方差（标准偏差）： 东华理工大学 张怀强 峰面积的计算 平均总峰面积法 在总峰面积法中，关键是确定峰的边界道l和h，但仍很难消除统计涨落的效应。 基本思想 从峰位置出发分别向两边找，在峰底部第一个最小值作为左、右边界道L、R开始，L、R逐次向外改变一道，即L→L-1, L-2,…；R→R-1, R-2,…，直到其中一边遇到新峰的边界或某一预定道为止。每改变一次边界，都按总峰面积法计算面积，然后由这些面积求出一个平均值，即为平均总峰面积。 东华理工大学 张怀强 峰面积的计算 平均总峰面积法 峰面积的计算 东华理工大学 张怀强 峰面积的计算 平均总峰面积法 标准偏差及总不确定度的计算 由每道上计数的统计涨落引起的。 由于峰边界道的选取引起的。 反映了道计数统计涨落和边界选择引起的总误差。 东华理工大学 张怀强 峰面积的计算 单峰曲线拟合法 基本峰形函数—用高斯函数来描述：（其中，a1为峰的高度，a2为峰位，a3为标准误差） 基底函数（基线函数）——用多项式描述： 则单峰拟合函数：data(i)=G(i)+B(i) 东华理工大学 张怀强 峰面积的计算 单峰曲线拟合法 单峰拟合函数：data(i)=G(i)+B(i) 用峰边界内的点的计数值作为一组已知的数据，代入拟合函数，用非线性最小二乘法，求解各系数（a1, a2, a3, b1, b2, b3） 则峰面积为： a1为峰的高度，a3为标准误差 东华理工大学 张怀强 重峰分析 重峰原因 射线能量相近。如：V的EKβ为5.43keV，Cr的EKα为5.41keV。 统计涨落特性：使每个能量的射线，得到一个高斯分布，即峰。 重峰类型 重合：高斯分布近似重合（能量相近） 组合峰：半边重叠，峰位独立 东华理工大学 张怀强 重峰分析 重峰分析 寻峰，判断是否重峰，并确定重峰个数，及各峰的能量； 选择重峰分析方法，按重峰个数进行解谱，求解各个峰的面积。 非线性最小二乘拟合法 分支比法（剥谱法） 东华理工大学 张怀强 重峰分析 非线性最小二乘拟合法 1.先确定单峰个数n，并确定各单峰的峰位 2.确定非线性最小二乘拟合函数 峰形函数—用高斯函数来描述： 其中，aj1为第j个峰的高度，aj2为第j个峰的峰位，aj3为第j个峰的标准误差，则：aj3= a3，所以，n个峰共有2n+1个系数。 基底函数（基线函数）—用多项式描述： 东华理工大学 张怀强 重峰分析 非线性最小二乘拟合法 2.确定非线性最小二乘拟合函数 则重峰的谱函数： 3.用非线性最小二乘拟合法，求解各系数: b1，b2，b3，aj1，aj2，a3 (j=1,2,…n) 东华理工大学 张怀强 重峰分析 非线性最小二乘拟合法 3.用非线性最小二乘拟合法，求解各系数。 4.则各单峰的面积为： 方法一： 方法二：求重峰总的净峰面积A，然后各个单峰按峰高成比例来分配： 当峰形函数不完全符合高斯分布时误差不大，所以常用 东华理工大学 张怀强 重峰分析 分支比扣除法（剥谱法） 1．如：152Sm(n, γ)135Sm—103.2keV 238U：103.7keV和277.6keV 103.2keV和103.7keV相近，重合 怎么办？？？ 1)用纯U测出谱线，求铀的两个峰103.7keV和277.6keV的分支比(不变)：FZBU=A103.7/A277.6=1.71 2)实测样品的Sm面积为： ASm=A103总-FZBU*A277.6=A103总-1.71*A277.6 东华理工大学 张怀强 重峰分析 分支比扣