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

%本次课程设计采用的谱数据为iaea-1995文件夹下iaearfnw TSTSPEC %里面的数据。首先来看看CALIB.ASC。READ_ME.TXT中说明了这个谱数据包含的部分峰的峰位与对应能量如下： % Channel Energy (keV) % 301 122.06 % 1281 511.00 % 1661 661.66 % 2097 834.84 % 2951 1173.24 % 3207 1274.54 % 3353 1332.50 %运行程序，其中参数选择为：选择傅里叶变换法平滑输入3，选择高斯滤波器输入2，然后A=1,FWHM=4，对称零面积法的参数是K=2，H=3,b=1寻出来 %的峰与READ_ME.TXT中说明的部分峰的峰位与对应能量数据相吻合。 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)=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 i=1:count %将平滑好的数据放回原数组，为下一次做好准备。 array_z(i)=array_smooth(i); end end for i=1:count a1(i)=array_z(i); end %***********************重心法平滑结束*************************** %***********************多项式最小二乘移动平滑法***************** elseif(pinghuaxuanze==2) biaoji=2; w=input(input the width of the filter window:); %w为窗口宽度 fwk = zeros(w,1); %存储滤波器系数，产生一个1行，w列的零矩阵；当求平滑之后谱的第i点数据时，先在原始数据第i点的左、右各取m个数据点；也就是形成一个共有w=2m+1个数据点的窗口； for i=1:w