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

%本次课程设计采用的谱数据为iaea-1995文件夹下iaearfnw TSTSPEC %里面的数据。首先来看看CALIB.ASCo READ_ME.TXT屮说明了这个谱数据包含的部分 峰的峰位与对应能量如下： % Chan nel Energy (keV) % 301 122.06 % 1281 511.00 % 1661 661.66 % 2097 834.84 % 2951 1173.24 % 3207 1274.54 % 3353 1332.50 %运行程序，其屮参数选择为：选择傅里叶变换法平滑输入3,选择高斯滤波器输入 后 A=1,FWHM=4, 对称零面积法的参数是K=2, H=3,b=l寻出来 %的峰与READ_ME.TXT屮说明的部分峰的峰位与对应能暈数据相吻介。 2,然 clc; clear; [Filename,Pathname]=uigetfile(,*.*,,i^ 择谱数据); fid=fopen([Pathname Filename],r) %fid 为文件指针，i■表示读操作 [array,count]=fscanf(fid,,%d,,[l inf]); %指定格式转换后返冋给矩阵 array,同 时返冋成功的读出的数据数量count, 1表示读出一个元素到一个列向量，inf表示读到文件 结朿返冋一个与文件数据元素相同的列向量 fclose(fid); %%%%ifii 开始能谱平 zH-%%%%%%% pinghuaxuanze=i叩ut(*请选择平滑方法:\n输入1选择重心法平滑\n输入2选择多项式最小二 乘移动平滑法\n输入3选择傅里叶变换法\n输入4选择小波变换:W); 1 ?￡? ?土 ?士 ?士 ?￡? ?土 ?士 ■一|、 ■丿? I ? \ I / % J f Ia ?士 ?土 ?士 ?士 1 v y. 1 ? ? w 1A ? ~ if (pinghuaxuanze== 1) biaoji=l; for i=l:count array_z(i)=array(i); end w=i叩ut「i叩ut the width of the filter window:); %w 表示 w 点平滑公式 while mod(w,2)==0 %判断输入的数是否是奇数，不是则重新输入。 w=input(4nput odd number/); end m=floor(w/2); for j=l:m for i=l:count if(i==l) array_smooth(i)=0.5*(array_z(i)+array_z(i+1)); %能谱左边界做对称镜像处理 elseif(i 1 i(count-l)) array_smooth(i)=0.25*array_z(i-l)+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=l:count %将平滑好的数据放冋原数组，为下一次做好准备。 array_z( i )=array_smooth( i); end end for i=l:count al(i)=array_z(i); end 逐 ****** ****** ***********重心 法平 滑结束 ********** ********** ******* 轻 ******* ****** *** ******* 多项式最小二乘移动平滑法 **** *** ** ***** *** elseif(pi nghuaxua nze==2) biaoji=2; w=input(*i叩ut the width of the filter window:); %w 为窗口 宽度 fwk = zeros(w,l); %存储滤波器系数，产生一个1行，w列的零矩阵；当求平滑 之后谱的笫i点数据时，先在原始数据第i点的左、右各取m个数据点；也就是形成一个共 有w=2m+l个数据点的窗口； for i=l :w k=i-ceil(w/2); %调整，计算采用k=- m:m,窗口的中心点为ceil( w/2)点 fwk(i)=[ 1 +(15/(wA2-4))*((wA2-1)/12-kA2)]/w; %2 次或 3 次滤波器 %fwkz=2.5*(vvA2-7)*((wA2-1)/12-kA2)-9*((wA2-1 )*(3*wA2-7)/240-kA4); %4 次或5次滤波器 %fwk(i)=( 1 +105/(4*(wA2-4)*(wA2-16))*fwkz)/w; %fwk(i)=( 1+(15 *(3 * w