核辐射测量数据处理成都理工大学程序设计报告..doc
文本预览下载声明
程序设计 课程设计(论文)
设计（论文）题目
谱的显示、谱光滑、定量分析（峰面积）程序
学院名称 核技术与自动化工程学院
专业名称 辐射防护与环境工程
学生姓名 袁子程
学生学号 201106080112
任课教师 马永红
设计（论文）成绩
教务处 制
2015年 1 月 10 日
实习目的
学习使用各种编程软件，利用计算机工具开展专业信息处理工作。
使用面向对象编程思想功能以及谱显示功能。题目为：
具体要求：
首先进行MFC绘图区域的学习和了解，实现基本绘图功能
MFC程序，设计一个合理美观的工作区
文件打开读取功能，并能打开不同文件
谱数据处理的算法编写
谱图形并给出处理结论单文档类MFC可执行程序
绘图区函数，初始化绘图设备
打开虚函数
文件及初步寻峰代码
谱图代码
读取mca文件
处理代码
，b=2m+1，为平滑宽度。
谱光滑 5 7 9 11 Kb 35 21 231 429 A0 17 7 179 143 A±1 12 6 135 120 A±2 -3 3 30 60 A±3 -12 -55 -10 A±4 15 -45 A±5 18
实际代码：
为） for( i=2;i1023;i++)
{
data[i]=1/35*((-3)*data[i-2]+12*data[i-1]+17*data[i]+12*data[i+1]+(-3)*data[i+2]);
}
②零面积对称法：
数学表达式：
其中，为变换后的谱数据，datai+j为原始实验谱数据，Cj为对称零面积变换函数，W=2m+1为窗宽（变换宽度）。
对称零面积变换函数：
Gj：
实际代码
//零面积对称法m=16,W=11,H=4)
for( i=16;i1007;i++)
{
for(j=0;j16;j++)
G[i]+=G[j];
d=1/11*G[i];
C[i]=G[i]-d;
}
for( i=16;i1007;i++)
{
if((yi/dy)f)
feng[k++]=data[i];
}
③Wasson（瓦森、沃森）峰面积法：
计算总面积：
计算本底面积：
计算净峰面积：
实际代码：
for (i=L;i=R;i++)
{
RL=right-left;
n=(right-left)/2;
im=(right+left)/2;
S+=data[i];
BL=(data[L]-data[R])/RL*(im-L-n)+data[L];
BR=(data[L]-data[R])/RL*(im-L+n)+data[L];
B=(BL+BR)*(n+1/2);
A=S-B;
}
计算峰位能量
/* for( i=0;i1023;i++)
{
E[i]=0.017584*i-0.138133;
}
for( j=0;j1024;j++)
{
if(E[dao[j]](6.403-6.403*0.1)||E[dao[j]](6.403+6.403*0.1))
ch[dao[j]]=Fe;
else if(E[dao[j]](7.477-7.477*0.1)||E[dao[j]](7.477+7.477*0.1))
ch[dao[j]]=Ni;
else if(E[dao[j]](8.047-8.047*0.1)||E[dao[j]](8.047+8.047*0.1))
ch[dao[j]]=Cu;
else if(E[dao[j]](8.638-8.638*0.1)||E[dao[j]](8.638+8.638*0.1))
ch[dao[j]]=Zn;
else if(E[dao[j]](10.5-10.5*0.1)||E[dao[j]](10.5+10.5*0.1))
显示全部