核数据处理课程设计--能谱谱数据分解方法研究程序.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:cou