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

%本次课程设计采用的谱数据为
iaea-1995 文件夹下 iaearfnwTSTSPEC
%里面的数据。首先来看看 。
中说明了这个谱数据包含的部分
峰的峰位与对应能量如下：
%
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