核数据处理课程设计

《核数据处理课程设计》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、核数据处理一、实习目的1、掌握图形程序的编写技能与技巧；2、掌握文件读写函数的使用；3、理解整个谱数据处理的流程；4、进一步掌握对常用核数据处理的基本方法5、掌握软件设计所需的基本动手能力，解决实际问题。二、实习内容1、读谱数据文件、保存光滑后的谱数据文件，文件名可选2、显示谱数据，各算法要以函数方式实现3、显示分析谱段范围内，寻得的各峰位道址、能量、对应的元素名称4、显示寻峰及峰面积计算结果5、光滑、寻峰宽度2k+1应在5、7之间可选，讨论最佳参数。三、需求分析1、输入、输出的内容和形式（1）输入：a、打开应用软件，弹出要求输入谱数

2、据文件名对话框。b、需要刻度谱线时，分别按e键和c键输入元素文件名和能量刻度文件名。c、人机交互时，采用键盘按键方式将命令传给程序（具体操作：在运行程序后按h键即可弹出操作说明对话框）。（2）输出：a、屏幕显示谱线；光标位置及光标处的道址、能量和计数；屏幕上方显示“Softwareinstructionspleaseinput'h'”。b、当谱线被刻度好后，按n键会弹出相应元素对应的能量、面积和净峰面积。c、按h键弹出应用程序操作说明对话框。2、预期功能、界面（1）显示输入谱文件名对话框；在文件名输入正确时，能够正确读取文件中的数据。

3、（2）正确显示谱线以及在刻度后显示道址、能量及计数。（3）显示相应元素对应的能量、面积和净峰面积对话框。（4）达到谱数据光滑的效果。（5）采用键盘输入命令，使屏幕：a、即时显示光标的位置以及对应的道址、能量及计数。b、即时显示不同幅度放大倍数时的谱线形状。c、即时显示不同光滑次数处理后的谱线微细结构（只有在幅度放大倍数很大时光滑效果才明显）。d、即时显示谱线的点和线的切换。3、测试数据（1）文件读入是否正确。（2）谱数据段的选择功能是否实现。（3）move和curse的边界值是否设置妥当。（4）谱数据光滑时前后几个数据（不能进行光滑）

4、是否与光滑次数的关系对应好。（5）道址和能量的换算是否正确。（6）寻峰是否正确。（7）最小能量差的计算是否正确。（8）根据能量查找核素是否正确。（9）净峰面积的计算是否有缺陷。（10）检查谱线幅度放大时是否会出现异常。四、概要设计1、开发工具：TurboC。2、主要算法：（1）、重心法：选取加权因子和归一化因子，使光滑后的数据成为原来数据的重心。常见的有5点和7点光滑。5点法：7点法：（2）、简单比较法（极值定峰法、IF函数找峰法）：a、峰的定义：满足，然后在datai-m至datai+m中找最大值道。b、常用的：5点、7点极大值法。

5、c、一般，用R=N0/Nb≥Rc确定峰是否有意义。N0为净峰幅度与基底之和，Nb为基底计数，Rc为设定值。d、k：找峰阈值，根据高斯统计概率分布，一般k取值：1—1.5。峰的左右边界道i-L和i+R的确定：（3）、线性本底法：计算总面积：计算本底面积：计算净峰面积：3、主程序的流程及各模块之间的关系。开始BGI初始化调用显示未处理前谱线函数及处理谱函数NY主程序流程：谱文件是否打开？等待Switch（按键）：各按键作用函数调用End主程序开始第一次画谱线函数originalLine()标注峰drawPeak()画坐标drawCoord

6、inate()文件读写函数寻峰elementPeakArea()计算峰面积calculatePeakArea()找峰边界findPeakBoundary()指定道值找峰findPeakByChannel()光标当前位置寻峰cursorPeakArea()能量刻度energyToChannel()计算峰面积calculatePeakArea()计算峰面积calculatePeakArea()Key？缩放谱线w（s）移动光标a（d）快速移动光标z（x）载入元素文件e载入标定文件c计算峰面积nKey？谱光滑m光标位置寻峰p还原谱线到刚刚载入

7、的状态r将谱线数据以ASCII码文件存储f点线切换0(9)退出程序q各模块之间的关系：等待有键按下？设定键被按下未设定键被按下五、详细设计1、函数及变量定义：#defineMAX_ELEMENT_NUM20#defineScreenHeight480#defineScreenWidth640#defineLineHeight370#defineLineWidth500#defineLineLeftX((ScreenWidth-LineWidth)/2)#defineLineRightX(LineLeftX+LineWidth)#def

8、ineLineTopY40#defineLineBottomY(LineTopY+LineHeight)FILE*dataFile;intgdriver,gmode;unsignedlongdataTable[1024]