核辐射测量数据处理成都理工大学程序设计介绍.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)) ch[dao[j]]=Ap; els