能谱谱数据分解方法研究程序核数据处理课程设计.doc

% 本次课程设计采用的谱数据为 iaea-1995 文件夹下 iaearfnw TSTSPEC % 里面的数据。首先来看看 。 中说明了这个谱数据包含的部分峰的峰位与对应能量如下: % Channel Energy (keV) % 301 % 1281 % 1661 % 2097 % 2951 % 3207 % 3353 % 运行程序,其中参数选择为:选择傅里叶变换法平滑输入 3 ,选择高斯滤波器输入 2 ,然后 A=1,FWHM=4 ,对称零面积法的参数是 K=2 , H=3,b=1 寻出来% 的峰与 中说明的部分峰的峰位与对应能量数据相吻合。 clc; clear; [Filename,Pathname]=uigetfile('*.*',' 选择谱数据'); fid=fopen([Pathname Filename],'r') %fid 为文件指针, r 表示读操作[array,count]=fscanf(fid,'%d',[1 inf]); % 指定格式转换后返回给矩阵 array ,同时返回成功的读出的数据数量 count ,1 表示读出一个元素到一个列向量, inf 表示读到文件结束返回一个与文件数据元素相同的列向量 fclose(fid); %%%% 下面开始能谱平滑%%%%%%% pinghuaxuanze=input(' 请选择平滑方法:\n 输入 1 选择重心法平滑\n 输入 2 选择多项式最小二乘移动平滑法\n 输入 3 选择傅里叶变换法\n 输入 4 选择小波变换:\n'); %************************ 重心法平滑**************************** if (pinghuaxuanze==1) biaoji=1; for i=1:count array_z(i)=array(i); end w=input('input the width of the filter window:'); %w 表示 w 点平滑公式 while mod(w,2)==0 % 判断输入的数是否是奇数,不是则重新输入。 w=input('input odd number:'); end m=floor(w/2); for j=1:m for i=1:count if(i==1) array_smooth(i)=*(array_z(i)+array_z(i+1)); % 能谱左边界做对称镜像处理 elseif(i>1&&i<(count-1)) array_smooth(i)=*array_z(i-1)+*array_z(i)+*array_z(i+1); else array_smooth(i)=*(array_z(count)+array_z(count-1)); % 能谱右边界做对称镜像处理 end end for i=1:count % 将平滑好的数据放回原数组,为下一次做好准备。 array_z(i)=array_smooth(i); end end for i=1:count a1(i)=arra