谱光滑.pptx

谱光滑谱光滑意义由于核衰变及测量的统计性，当计数较小时，计数的统计涨落比较大，计数最多的道不一定是高斯分布的期望值，真正的峰被淹没在统计涨落之中。为了在统计涨落的影响下，能可靠的识别峰的存在，并且准确地确定峰的位置和能量，从而完成定性，就需要进行谱光滑由于散射的影响，峰的两边即峰的边界受统计涨落干扰大，影响峰面积的计算，从而影响准确的定量分析，所以必须对谱数据进行数据光滑。成都理工大学 马英杰谱光滑成都理工大学 马英杰基本思想1）谱数据的光滑处理是逐点进行的2）以待处理点为中心，用其左右m个点的测量数据，对该点作修正，消除统计涨落的影响；3）保留原有意义的特征信息。 峰位、分辨率 能量、峰面积谱光滑成都理工大学 马英杰方法：1）算术滑动平均法（平均移动法）2）重心法3）多项式最小二乘拟合法4）离散函数褶积滑动变换法5）傅立叶变换法其它。。。。。谱光滑算术滑动平均（平均移动法）设x0=i为待光滑点坐标，左右各取m个点(j=i, i±1, i±2, …, i±m)，则共有2m+1个点，用所有2m+1个点的算术平均值作为这点的改正值。此方法：两端各有m个点得不到平滑，称为边沿损失。成都理工大学 马英杰谱光滑算术滑动平均（平均移动法）m=1, 三点法m=2, 五点法 成都理工大学 马英杰谱光滑算术滑动平均（平均移动法）编程：void CMmcaView::SmoothAverage(int data[], int data1[], int Beginch, int Endch, int m){ int temp; for (int i=Beginch; i=Endch; i++) { temp=0; for(int j=-m; j=m; j++) temp=temp+data[i+j]; data1[i]=temp/(2*m+1); }}成都理工大学 马英杰谱光滑成都理工大学 马英杰重心法基本方法：选取加权因子和归一化因子，使光滑后的数据成为原来数据的重心。常见的有3点、5点和7点光滑。谱光滑成都理工大学 马英杰重心法例如：3点法i-1 i-0.5 i i+0.5 i+1第i道和第i+1道计数的重心应是：第i道和第i-1道计数的重心应是：由于道数是整数，没有半道的情况存在，则由上两式再求重心得：谱光滑重心法5点法：成都理工大学 马英杰谱光滑void CMmcaView::SmoothBaryCenter(int data[], int data1[], int Beginch, int Endch, int m){ int temp[2*m+1][m+1]; for (int j=Beginch; j=Endch; j++) { for(int k=0; k=m; k++) { for(int i=k; i=2*m-k; i++) { if(k==0) temp[i][k]=data[j+i-m]; else temp[i][k]=(temp[i-1][k-1]+2*temp[i][k-1]+temp[i+1][k-1])/4; } } data1[j]=temp[m][m]; }}成都理工大学 马英杰重心法7点法：编程一：按给定点数循环计算谱光滑void CMmcaView::SmoothBaryCenter1(int data[], int data1[], int Beginch, int Endch, int m){ for (int i=Beginch; i=Endch; i++) { if(m==1) data1[i]=(data[i-1]+2*data[i]+data[i+1])/4; else if(m==2) data1[i]=(data[i-2]+4*data[i-1]+6*data[i] +4*data[i+1]+data[i+2])/16; else data1[i]=(data[i-3]+6*data[i-2]+15*data[i-1]+20*data[i] +15*data[i+1]+6*data[i+2]+data[i+3])/64; }}成都理工大学 马英杰重心法编程二：直接按点数选择公式计算谱光滑多项式最小二乘拟合法基本思想用一个n次多项式与W=2m+1个谱数据点逐次分段进行拟合，以达到光滑的目的。 设在谱数据中取2m+1个等距点： x -m，x-m+1，…，x-1，x0，x1，…，xm-1，x m对应数据 y- m，y- m+1，…，y-1，y0，y1，…，ym-1，y m令h为xi之间的距离，做变换i=(x-x0)/h，则上述各点将变为： -m，-m+1，…，-1， 0， 1，……，m-1，m成都理工大学 马英杰谱光滑多项式最小二乘拟合法基本思想假定h=1，用一