3 能谱数据平滑 核数据处理课件.ppt

核 数 据 处 理 * * 核 素 指 纹 每种核素辐射γ射线的本征能谱不同 ■ 特征γ射线的数量和能量不同 核素的γ射线指纹 ■ 不同能量γ射线的强度比例不同 只要测出核材料的γ 射线能谱，找出各个特征线的能量和强度，根据每个核素的指纹，就可以决定此和材料中包含那些核素？核素的含量。 由能谱采集系统（γ能谱仪）测量得到的γ 能谱。 不是被测样品的γ 辐射能量分布（本征能谱）。 而是谱仪主要是探测器对于被测样品辐射的γ特征射线的能量响应。 第三章能谱数据平滑 主讲：张怀强 东华理工大学核工程与地球物理学院 在给定时刻，处于特定能态的一定量放射性核素在时间内发生自发核跃迁次数的期望值。 ■ 在一定时间间隔内，核辐射测量中的放射性事件的发生时刻和发生次数都是随机的，即放射性事件的具有统计涨落性。 ■ 放射性活度 特定能态的放射性核素的数目衰减到原数目的一半时，所需要时间的期望值。 ■ 半衰期T1/2 处于特定能态的一定量的放射性核素的平均生存时间。 ■ 放射性核素的平均寿命 在以下情况下 ■ 放射源活度较低 ■ 一定时间内探测器接收到的粒子总数存在统计涨落 ■ 测量时间较短 ■ 待测特征射线分支比较小 ■ 粒子与探测器发生相互作用形式很多 ■ 出现假峰 ■ 丢失弱峰 ■ 峰面积计算误差加大 ■ 能谱数据平滑的作用 减少统计涨落对后续处理的影响。 ■ 平滑前后能谱的特征 ★ 峰位与峰面积尽可能保证稳定变化很小 ★ 尽可能减小重峰出现的几率 3.1 重心法 则在t时刻物体的质心 3.1 重心法 假设能谱相邻两道所受统计涨落影响之和为0，根据质心不变，可知 由于能谱不存在半道，但 三点重心法 3.1 重心法 该方法认为两道所受的统计涨落为0，也就是说频率一定低频，说明过滤了高频。 同理 五点重心法 该滤波器权重因子均为正数，平滑之后不会出现负值，提高平滑数据的可靠性。 2n+1点重心法呢？ 七点重心法？ 2次三点重心法 3.1 重心法 for i=1:count if(i==1) array_smooth(i)=0.5*(array_z(i)+array_z(i+1)); elseif(i1i(count-1)) array_smooth(i)=0.25*array_z(i-1)+0.5*array_z(i)+0.25*array_z(i+1); else array_smooth(i)=0.5*(array_z(count)+array_z(count-1)); end end 能谱左边界做对称镜像处理 能谱右边界做对称镜像处理 for j=1:m for i=1:count array_z(i)=array_smooth(i); end end 判断输入的数是否是奇数，不是则重新输入。 w=input(input the width of the filter window:) while mod(w,2)==0 w=input(input odd number:); end m=floor(w/2); 数据放回原数组，为下一次做好准备。 W点重心法程序 3.2 多项式最小二乘移动平滑法 设 当 3.2 多项式最小二乘移动平滑法 fwk = zeros(1,w); for i=1:w k=i-ceil(w/2); fwk(i)=[1+(15/(w^2-4))*((w^2-1)/12-k^2)]/w; %fwkz=2.5*(w^2-7)*((w^2-1)/12-k^2)-9*((w^2-1)*(3*w^2-7)/240-k^4); %fwk(i)=(1+105/(4*(w^2-4)*(w^2-16))*fwkz)/w; %fwk(i)=(1+(15*(3*w^2-7)/((w^2-4)*((w^2-5)^2+4)))*((w^2-1)*(3*w^2-7)/240-k^4))/w; end 存储滤波器系数 调整，计算采用k=-m:m，窗口的中心点