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

1、成却咒z幺程 序设计 课程设计（论文）设计（论文）题目谱的显示、谱光滑、定性分析程序学院名称 核技术与自动化工程学院专业名称辐射防护与环境工程学生姓名蒋卓辰学生学号201 1 060801 1 6任课教师马永红设计（论文）成绩教务处制2015年 1月 10日一、实习目的学习使用各种编程软件，利用计算机工具开展专业信息处理工作。使用面向对象编程思想实现谱数据处理功能以及谱显示功能。二、人员组成及分工本人题目为：谱的显示、谱光滑、定性分析程序具体要求：用算术滑动平均法实现谱光滑一阶导数法寻峰位根据能量刻度系数计算峰位的能量，检索元素三、实习计划1. 首先进行 mfc 绘图区域的学习和了解，实现基本

2、绘图功能2. 建立单文档 mfc 程序，设计一个合理美观的工作区3. 实现文件打开读取功能，并能打开不同文件4. 进行谱数据处理的算法编写5. 编写成谱图形功能并给出处理结论四、基本原理1、重心法选取加权因子和归一化因子，使光滑后的数据成为原来数据的重心。常见的有 5 点和7点光滑。15 点法：datai = (datai +4datai+6datai +4datai+1 + datai42) _ -1 ,7 点法：data = (dataij3 +6data +15dataa +20data +15data卅 +6data书 +data书) 642、简单比较法（极值定峰法、if函数找峰法）a

3、、峰的定义： 满足 datai _m < datai - kdatai > data4m，然后在 data i-m 至 data i+m中找最大信道。b、常用的：5点、7点极大值法。c、一般，用r=n/nbr确定峰是否有意义。n）为净峰幅度与基底之和，n为 基底计数，r为设定值。d、k:找峰阈值，根据高斯统计概率分布，一般 k取值：1 1.5。峰的左右边界道i-l和i+r的确定:datagm -datai k-. dataidata。.加 一 datai . k . datai .3、线性本底法r计算总面积：s=" datai i =l计算本底面积： b =（datal

4、data * （r ；1）计算净峰面积：a=s-b五、主要工作介绍1.建立单文档类mfc可执行程序mfc应用程序向导稣 1ce宝您要创建的应用程序类型是：岸单文档冏多重文档网基本对话框13丁文档度看体系结构支持m您的资源使用的语言是:|中文晡体 中国appwzc忘口了 二上一步 下少 完成 | 取消 |2.建立绘图区函数，初始化绘图设备jiang - microsoft vislmii c+ +忘件但克等嗝亘吞任i疝工程凹咀遭凶工却工if口直i尾和为周宁q。fcjiniiflivicw£- 康 jla)的痴££吧4-口 jinngviw li5 cabouldig

5、 t cjha口加 p "t' cjiiaiigdae hi- cjiangvicw i cmaiinframc 二( glubals/ jidniguirw.lib 二 inttrfact nf thr 匕innqui册w 匚 1上£b /blf mined(>flifx_jlanisuiew_hi_i|ginbt3d_319e_ftcfifi_b：7m_26s7ar2hcf50_ihclude口 tidrfiar hfx jhancuifvf n 庙皿*73。覆 19e 虻的 a7im 26自7正2城f9d thclugedfllf _msc_uefi10

6、04刷mqz oiicffltfldkf / _hsc_uefi > 1000class cjlangview : public cttiemprotected: / create fron serializaitian only cjinguieuf);l>lcl bre_cvhchtol e (cjidnguiew)/ lattriihutfffpubliic:cji-angduc» getooctinentf);/ dp?r.ations publics/ derrides/gereirated virtual function uverrides”(flifxui

7、 rtual( cjianguieui)public: virtual void uirtual idddl pre tbcttd ;(jndraw(cbc* pdc); / overridden to drw this view fiacre at ewinidgwfic rea testiucta cs)：virtual virtudl virtual /bfifk.dddl void voidvirtual0nprtparefrinting(gprintinldfl plnfo);flihbrgiabprintiing(icqc*icpriritififij* |ilnf d)；onen

8、ilpriratingtcdc* pdca cprintlnfo* plnfo);/ inplen&nt.atidnpuli lit)virtyai 飞jiwnguiewlj;tlfdbf _debuguirltiiidl wi d tesrr tuialidovirtual void du>np(counfjcontextk de) const; 知gihhfprotect把配 / generated nessage nap functionsprcte'eted:/1/ftfksgtcjiaingmiewj/ hote - the classyizard uill a

9、dd and renoue nnber fufictions her&a/ do not edit wtit ynu in th廿国 hluich% nf yifnvrtcd eodr t111mkt n3.添加文件打开虚函数4,尝试读取mca文件六、流程及程序1、模块介绍(1)主程序的流程开始bgi初始化调用显示未处理前谱线函数及处理谱函数switch (按键)：各按键作用函数调用end2、程序流程(1)函数及变量定义#define max_element_num 20#define screenheight 480#define screenwidth 640#define lin

10、eheight 370#define linewidth500#define lineleftx(screenwidth-linewidth)/2)#define linerightx(lineleftx+linewidth)#define linetopy40#define linebottomy (linetopy+lineheight)file* datafile;int gdriver,gmode;unsigned long datatable1024,maxval;int dispos500;/* 实际显示的数据*/int oripos500;/* 未处理前的数据*/int curs

11、orx;float zoom;unsigned int cursorchannel;unsigned char displine;unsignedcharenergymarked,elementloaded,spectrumloaded,peakseeked,cursorpeakseeked;int smoothlevel = 1,peakseeklevel = 9;float a,b;/* 能量刻度方程*/char elementmax_element_num2;unsigned char numofelement,cursorelement;float elementenergymax_e

12、lement_num;unsigned int peakboundmax_element_num2;longunsignedelementareamax_element_num,elementnetareamax_element_num;unsigned int cursorpeakboundleft,cursorpeakboundright;unsigned long cursorpeakarea,cursorpeaknetarea;( 2)程序 #define max_element_num 20#define screenheight 480#define screenwidth640#

13、define lineheight370#define linewidth500#define lineleftx(screenwidth-linewidth)/2)#define linerightx(lineleftx+linewidth)#define linetopy40#define linebottomy(linetopy+lineheight)file* datafile;int gdriver,gmode;unsigned long datatable1024,maxval;int dispos500;/*实际显示的数据*/int oripos500;/* 未处理前的数据*/i

14、nt cursorx; float zoom;unsigned int cursorchannel;unsigned char displine;unsigned char energymarked,elementloaded,spectrumloaded,peakseeked,cursorpeakseeked;int smoothlevel = 1,peakseeklevel = 9;float a,b;/*能量刻度方程 */ char elementmax_element_num2;unsignedcharnumofelement,cursorelement;floatelementene

15、rgymax_element_num;unsigned int peakboundmax_element_num2;unsignedlongelementareamax_element_num,elementnetareamax_element_num;unsigned int cursorpeakboundleft,cursorpeakboundright;unsigned long cursorpeakarea,cursorpeaknetarea;w=input('input the width of the filter window:');for i=1:w k=i-c

16、eil(w/2);end array_z = zeros(count+2*floor(w/2),1);for i=1:count+2*floor(w/2) if(i<=floor(w/2) array_z(i)=array(-i+ceil(w/2);elseif(i>count+floor(w/2) array_z(i)=array(-(i-floor(w/2)+2*count+1);else array_z(i)=array(i-floor(w/2); end end a1=zeros(1,count);for j=1:w smz(j)=array_z(i+j-1);enda1(

17、i)=smz*fwk; enddisp('下面开始输入对称零面积法寻峰)；disp（'下面开始输入对称零面积法的各参数）;disp('如果是方波的话有k=1');k=input(' 请输入参数k=?:n');h=input('请输入参数半宽度 h=?(正奇数):n');m=(2*k+1)*h-1)/2; w=2*m+1; b=input(' 请输入参数b=?:n');a=2*k*b; %k=4; %h=2*k+1; %w=3*h;%b=1; %a=2*k*b; m1=floor(w/2);temporary=zer

18、os(count+2*m1),1);for i=1:count+2*m1 if(i<=m1);temporary(i)=a1(ceil(w/2)-i);elseif(i>(count+m1) temporary(i)=a1(-(i-m1)+2*count+1);else temporary(i)=a1(i-m1); end end a=zeros(count,1);for i=ceil(w/2):count+m1; for j=-(w-1)/2:(w-1)/2;if abs(j)<=(h-1)/2; t=a;else t=-b;end a(i-m1,1)= a(i-m1,1)

19、+t*temporary(i+j);end end for i=1:count; %数据转制ssifenzi(i,1)=a(i,1);end b=zeros(count,1);for i=ceil(w/2):count+m1;for j=-(w-1)/2:(w-1)/2;if abs(j)<=(h-1)/2;t=aa2; else t=ba2;end b(i-m1,1)=b(i-m1,1)+t*temporary(i+j);endend for i=1:count;ssifenmu(i,1)=b(i,1);end for i=1:count;ss(i,1)=ssifenzi(i,1)/s

20、qrt(ssifenmu(i,1);end p=1; q=1; f=30;for i=1:count;if ssifenzi(i)<0;fpdatablow(p,1)=i;fpdatablow(p,2)=ssifenzi(i);p=p+1;elseif ss(i)>f;fpdataup(q,1)=i;fpdataup(q,2)=ssifenzi(i);q=q+1; end end p=1; for i=2:length(fpdataup(:,1)-1;iffpdataup(i,2)>fpdataup(i+1,2)&&fpdataup(i,2)>fpdat

21、aup(i-1,2);mpeak(p,1)=fpdataup(i,1); p=p+1; end end for i=1:length(mpeak(:,1);j=mpeak(i);t=mpeak(i);peak(i)=t+(a1(t+1)-a1(t-1)/(2*a1(t)-a1(t+1)-a1(t-1)/2;disp（'计算峰面积；下面采用总峰面积法）;y1=zeros(1,length(mpeak);for i=1:length(mpeak)for j=(zuobianjic(i)+1):(youbianjic(i)-1)y1(i)=y1(i)+a1(j);ends(i)=y1(i)-(youbianjie(i)-1)*(a1(zuobianjie(i)/2;endendenddisp（&