核辐射测量数据处理成都理工大学程序设计报告

程序设计课程设计(论文)

设计(论文）题目

谱的显示、谱光滑、定量分析(峰面积）程序

学院名称核技术与自动化工程学院ﻩ

专业名称辐射防护与环境工程

学生姓名袁子程

学生学号20１1０60801１2

任课教师马永红

设计（论文）成绩

教务处制

2015年1月10 日

一、ﻬ实习目的

学习使用各种编程软件，利用计算机工具开展专业信息处理工作。

使用面向对象编程思想实现谱数据处理功能以及谱显示功能。二、人员组成及分工

本人题目为：谱的显示、谱光滑、定量分析(峰面积)程序

具体要求:

多项式最小二乘拟合实现谱光滑

对称零面积法实现寻峰

Wasson(瓦森、沃森)峰面积法实现峰面积计算

寻峰应在5、7、9、１１之间可选,讨论最佳参数

三、实习计划

1.首先进行MFC绘图区域的学习和了解，实现基本绘图功能

2.建立单文档MFＣ程序,设计一个合理美观的工作区

3.实现文件打开读取功能，并能打开不同文件

4.进行谱数据处理的算法编写

5.编写成谱图形功能并给出处理结论

四、主要工作介绍（详细介绍实习中各功能的实现方法、理论公式和计

算原理,程序流程图和程序运行界面截图等成果信息）

1.建立单文档类MFC可执行程序

2.建立绘图区函数,初始化绘图设备

3.添加文件打开虚函数

4.编写读取文件及初步寻峰代码

5.添加绘制谱图代码

6. 尝试读取mca 文件

7. 添加谱数据处理代码

①多项式最小二乘拟合法:

∑-=+=

m

m

j j

i j

b

i data

A K data 1

,ｂ=2m ＋１，为平滑宽度。 谱光滑

5

７

9

11

实际代码：

／/多项式最小二乘法谱光滑(参数为5) ﻩﻩ ｆor( i=2;ｉ<1023；i ＋+) ﻩ ﻩ｛

d ａta[ｉ］=1/35*(（-3）*d ａta[i-2]＋1２*ｄata ［i-１]＋17*d ａta ［i ］+

１2＊ｄa ｔａ[ｉ+1］+(-３)*dat ａ[i+2]); ﻩﻩﻩ ｝

②零面积对称法：

基本思想：面积为零的“窗”函数与实验谱数据进行褶积变换,且要求“窗”函数为对称函数。

数学表达式：j j m

m

j j

m

m

j j

i j

i C C C

data

C y --=-=+===∑∑0

~

其中，i y ~

为变换后的谱数据，ｄata i+j 为原始实验谱数据，Ｃj 为对称零面积变换函数,W=2ｍ＋1为窗宽(变换宽度)。

对称零面积变换函数:∑-==

-=m

m

j j

j j G

W

d d

G C 1

Ｇj ：[]

⎪⎪

⎪⎩⎪⎪⎪

⎨⎧==+=-=)634.2(sec ))2(cos )4))(2ln 4exp )22

2

2

2H j h G d H j G c j H H G b H j G a j j j j 双曲正割函数：

余弦平方函数：柯西函数：高斯函数：π

实际代码:

//零面积对称法寻峰(ｍ=16,W ＝１1,Ｈ=4) ﻩf ｏr （ i=16;ｉ<10０7;ｉ++) ﻩ{ ﻩ

ｆor(ｊ=0；j<１６;j ＋+）

ﻩ ﻩG[i ］+＝G ［ｊ];

ﻩﻩ d=１/11*G[i];

ﻩﻩ ﻩC ［i]＝G ［i ］-d ； ﻩ ﻩ｝

ﻩ f ｏｒ( i=１6;i<1007;ｉ++) ﻩ {

ﻩ ｉf((yi/dy)>f) ﻩﻩ fen ｇ[k ＋+]=da ｔａ［i];

ﻩﻩ }

③Ｗas ｓo ｎ(瓦森、沃森)峰面积法: 计算总面积：∑-=+=

n

n

i i

i data

S 0

计算本底面积:

L L

R n data n L i L R data data b +----=

-)(0

L L

R n data n L i L

R data data b ++---=)(0

)21)((2

1

2)(++=++=--n b b n b b B n n n n

计算净峰面积:B S A -=

实际代码:

/／Wａsson(瓦森、沃森)峰面积法：

ﻩfor(i=L;i<=R；i＋+)

ﻩ{

ﻩﻩﻩＲL=riｇht-leｆt；

ﻩﻩn=(ｒight-leｆｔ)／2;

ﻩﻩim=(riｇht+left)/２;

ﻩﻩS+=ｄata[ｉ];

ﻩﻩBＬ=(daｔa［L］-dａｔa[Ｒ])/RＬ＊(ｉm-Ｌ-ｎ)＋ｄata[L］;

ﻩﻩﻩﻩBＲ=(daｔa［Ｌ]－ｄaｔa[R］)/RL*（im-L+n）+daｔa[L];

ﻩB=(BＬ+BR)＊(n+1/２）;

A=S-Ｂ;

ﻩﻩ｝

计算峰位能量

//根据能量刻度系数求峰位能量

/* ﻩｆｏr（i=0;i<１02３；i+＋)

｛

Ｅ[i]=０.０1758４*i-0.13８１33;

}

ｆor（j＝0;ｊ<1024;j＋+)

{

if(E［daｏ[j]］<（6．４0３-６.４03＊0.1)|｜E[dａo［ｊ]]>(6.403+６．403*