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

％本次课程设计采用的谱数据为ｉａｅａ－19９５文件夹下iaearfnwﻩＴＳＴＳPEC
%里面的数据。首先来看看CＡLIB．ASC。ＲＥAD_ＭＥ．TXT中说明了这个谱数据包含的部分峰的峰位与对应能量如下:
%　 　Channｅl 　 Ｅnerｇy (ｋeＶ)
％ 30１ 　 　　　
％ﻩ １２81 　
％ﻩ　 1661 　　　 　
％ﻩ ２0９7　 　 8３４．84
％ﻩ 　 2９５1　　 　 　1173．２４
%ﻩ 　 32０7
%ﻩ 　　３353 　 　
%运行程序,其中参数选择为:选择傅里叶变换法平滑输入3,选择高斯滤波器输入2,然后Ａ=1,FWHM=4,对称零面积法的参数是K=2,Ｈ=３,ｂ=1寻出来
%。
clｃ;
cleaｒ;
[]=uigetfilｅ（＇*.*＇,'选择谱数据＇）;
fid=ｆｏpen（[Pathnaｍｅ ],'r'）　 　　 　 　 　 %ｆｉd为文件指针,r表示读操作 　
［arraｙ,ｃount]=fscaｎf（fｉd,'％d',［1 iｎf］);　 　 　　　%指定格式转换后返回给矩阵array,同时返回成功的读出的数据数量cｏunt,1表示读出一个元素到一个列向量,ｉｎf表示读到文件结束返回一个与文件数据元素相同的列向量
fｃｌosｅ(fｉd）;
%%%%下面开始能谱平滑%％%％％%％
pｉｎghuaxuanze＝inｐut('请选择平滑方法:\n输入１选择重心法平滑\n输入2选择多项式最小二乘移动平滑法＼ｎ输入3选择傅里叶变换法\n输入4选择小波变换:\n');
%**＊*********＊＊**＊**＊**＊*重心法平滑**＊***＊****＊**＊***＊*＊＊*****＊
if (pinghuaｘuanze＝=1)
　biaｏｊi=１;
　 　foｒ　i=１:count 　 　
　　 　 aｒraｙ＿ｚ(i）=ａｒray（i);
　　 end
w=iｎput(＇ｉnｐuｔ　tｈｅ widｔh　ｏf ｔｈe fｉｌter wiｎdow:'); 　 　　 %w表示ｗ点平滑公式
　while　mod（w,2）==0 　　 　 %判断输入的数是否是奇数,不是则重新输入。
　 　 w=inpｕt（＇inpｕt　ｏdd ｎumber:');　
　 　　eｎd
m=flooｒ(w/2);
for j=1:m 　
foｒ i=1:ｃｏuｎt
　iｆ(i==1)
　　 arｒａｙ＿smooth(i)=*(arraｙ_z(ｉ)+arraｙ_z(ｉ＋1）); ％能谱左边界做对称镜像处理
　 　 eｌｓeｉｆ（i>1&&i<(coｕnt－１）)
　 　 　 　 arｒay＿sｍｏotｈ(ｉ)=0．2５＊arｒaｙ_ｚ(i-1)+0．５*arｒay_z(i)＋