核輻射測(cè)量數(shù)據(jù)處理成都理工大學(xué)程序設(shè)計(jì)報(bào)告資料(共12頁(yè))

1、 程序設(shè)計(jì)(chn x sh j) 課程(kchng)設(shè)計(jì)(shj)(論文)設(shè)計(jì)（論文）題目 譜的顯示、譜光滑、定量分析（峰面積）程序 學(xué)院名稱 核技術(shù)與自動(dòng)化工程學(xué)院 專業(yè)名稱 輻射防護(hù)與環(huán)境工程 學(xué)生姓名 袁子程 學(xué)生學(xué)號(hào) 201106080112 任課教師 馬永紅 設(shè)計(jì)（論文）成績(jī) 教務(wù)處 制2015年 1 月 10 日實(shí)習(xí)目的學(xué)習(xí)使用各種編程軟件，利用計(jì)算機(jī)工具(gngj)開(kāi)展專業(yè)信息處理工作。使用(shyng)面向?qū)ο缶幊趟枷?sxing)實(shí)現(xiàn)譜數(shù)據(jù)處理功能以及譜顯示功能。人員組成及分工本人題目為：譜的顯示、譜光滑、定量分析（峰面積）程序具體要求：多項(xiàng)式最小二乘擬合實(shí)現(xiàn)譜光滑對(duì)稱零面

2、積法實(shí)現(xiàn)尋峰Wasson（瓦森、沃森）峰面積法實(shí)現(xiàn)峰面積計(jì)算尋峰應(yīng)在5、7、9、11之間可選，討論最佳參數(shù)實(shí)習(xí)計(jì)劃首先進(jìn)行MFC繪圖區(qū)域的學(xué)習(xí)和了解，實(shí)現(xiàn)基本繪圖功能建立單文檔MFC程序，設(shè)計(jì)一個(gè)合理美觀的工作區(qū)實(shí)現(xiàn)文件打開(kāi)讀取功能，并能打開(kāi)不同文件進(jìn)行譜數(shù)據(jù)處理的算法編寫(xiě)編寫(xiě)成譜圖形功能并給出處理結(jié)論主要工作介紹（詳細(xì)介紹實(shí)習(xí)中各功能的實(shí)現(xiàn)方法、理論公式和計(jì)算原理，程序流程圖和程序運(yùn)行界面截圖等成果信息）建立單文檔類MFC可執(zhí)行程序建立(jinl)繪圖區(qū)函數(shù)(hnsh)，初始化繪圖設(shè)備添加(tin ji)文件打開(kāi)虛函數(shù)編寫(xiě)(binxi)讀取文件(wnjin)及初步尋峰代碼添加(tin ji)

3、繪制譜圖代碼嘗試(chngsh)讀取mca文件(wnjin)添加(tin ji)譜數(shù)據(jù)處理代碼多項(xiàng)式最小二乘擬合法：，b=2m+1，為平滑寬度。譜光滑57911Kb3521231429A0177179143A1126135120A2-333060A3-12-55-10A415-45A518實(shí)際(shj)代碼(di m)：/多項(xiàng)式最小二乘法(chngf)譜光滑（參數(shù)為5） for( i=2;i1023;i+)datai=1/35*(-3)*datai-2+12*datai-1+17*datai+12*datai+1+(-3)*datai+2); 零面積對(duì)稱法：基本思想：面積為零的“窗”函數(shù)與實(shí)驗(yàn)

4、譜數(shù)據(jù)進(jìn)行褶積變換，且要求“窗”函數(shù)為對(duì)稱函數(shù)。數(shù)學(xué)表達(dá)式：其中，為變換后的譜數(shù)據(jù)，datai+j為原始實(shí)驗(yàn)譜數(shù)據(jù)，Cj為對(duì)稱零面積變換函數(shù)，W=2m+1為窗寬（變換寬度）。對(duì)稱零面積變換函數(shù)：Gj：實(shí)際代碼：/零面積(min j)對(duì)稱(duchn)法尋峰(m=16,W=11,H=4)for( i=16;i1007;i+)for(j=0;j16;j+)Gi+=Gj;d=1/11*Gi;Ci=Gi-d;for( i=16;if)fengk+=datai;Wasson（瓦森、沃森）峰面積(min j)法： 計(jì)算總面積：計(jì)算本底面積： 計(jì)算凈峰面積：實(shí)際代碼：/ Wasson（瓦森、沃森）峰面積法：

5、for (i=L;i=R;i+)RL=right-left;n=(right-left)/2;im=(right+left)/2;S+=datai;BL=(dataL-dataR)/RL*(im-L-n)+dataL;BR=(dataL-dataR)/RL*(im-L+n)+dataL;B=(BL+BR)*(n+1/2);A=S-B;計(jì)算(j sun)峰位能量(nngling)/根據(jù)能量刻度(kd)系數(shù)求峰位能量/*for( i=0;i1023;i+)Ei=0.017584*i-0.138133;for( j=0;j1024;j+)if(Edaoj(6.403+6.403*0.1)chdaoj

6、=Fe;else if(Edaoj(7.477+7.477*0.1)chdaoj=Ni;else if(Edaoj(8.047+8.047*0.1)chdaoj=Cu;else if(Edaoj(8.638+8.638*0.1)chdaoj=Zn;else if(Edaoj(10.5+10.5*0.1)chdaoj=Ap;else if(Edaoj(11.725+11.725*0.1)chdaoj=As;else if(Edaoj(12.611+12.611*0.1)chdaoj=Ni;*/調(diào)試(dio sh)運(yùn)行(ynxng)體會(huì)(thu)在本次(bn c)程序設(shè)計(jì)(chn x sh j)實(shí)

7、習(xí)中，我再次學(xué)習(xí)了MFC編程方法。過(guò)程中，我遇到了許多曾經(jīng)都沒(méi)有遇到的問(wèn)題，比如文件的打開(kāi)，文件按鈕虛函數(shù)的添加，繪圖無(wú)效等。在一步一步解決這些問(wèn)題的同時(shí)，漸漸進(jìn)一步的熟悉掌握了譜數(shù)據(jù)處理的方法和編程思想。但是依舊有一個(gè)問(wèn)題尚未解決，那就是數(shù)據(jù)類型的變化。這個(gè)使得我無(wú)法在數(shù)據(jù)處理中正確地獲取應(yīng)有的結(jié)果。令我費(fèi)解的是成譜圖像卻是與原始數(shù)據(jù)非常的吻合。對(duì)于這個(gè)問(wèn)題我依然百思不得其解。建議我對(duì)本門(mén)課程的掌握程度也非常有限，提出以下建議供老師和學(xué)弟學(xué)妹們參考。1.MFC編程時(shí)一個(gè)系統(tǒng)過(guò)程，需要掌握的細(xì)節(jié)可以隨時(shí)網(wǎng)上學(xué)習(xí)，但是總體思路要牢記，不容易臨時(shí)掌握。2.對(duì)于本專業(yè)同學(xué)們而言，這次實(shí)習(xí)是擴(kuò)展我專業(yè)

8、同學(xué)技能的大好機(jī)會(huì)，建議不要輕易放棄，認(rèn)真編寫(xiě)，即使最后沒(méi)有達(dá)到自己預(yù)期的計(jì)劃也是收獲頗豐的。3.在編寫(xiě)的過(guò)程中要有良好的格式及注釋習(xí)慣，以免一段時(shí)間甚至幾分鐘后就看不懂自己的思路。4.在編寫(xiě)算法過(guò)程中，一定要回歸(hugu)理論，弄明白實(shí)際數(shù)學(xué)過(guò)程，對(duì)算法的編寫(xiě)很有幫助。5.希望(xwng)老師能夠?qū)Ψ撬惴ú糠纸o予(jy)同學(xué)更多的幫助，教授同學(xué)建立軟件的過(guò)程，指導(dǎo)同學(xué)完成非專業(yè)的部分，讓同學(xué)著眼于數(shù)學(xué)導(dǎo)算法的翻譯。最后，衷心感謝在本次程序設(shè)計(jì)實(shí)習(xí)過(guò)程中，馬永紅老師細(xì)心的指導(dǎo)，您不厭其煩地解答我的各種疑問(wèn)，并且?guī)椭覍ふ以紨?shù)據(jù)的數(shù)據(jù)格式和數(shù)據(jù)類型。我還要感謝蘆瑞琪、蔣卓辰同學(xué)對(duì)我在思路上的

9、重要啟發(fā)，謝謝你們給予我的大力支持！附：程序包附在附件中，項(xiàng)目名稱為37gai學(xué)生學(xué)習(xí)心得在本次程序設(shè)計(jì)實(shí)習(xí)中，我再次學(xué)習(xí)了MFC編程方法。過(guò)程中，我遇到了許多曾經(jīng)都沒(méi)有遇到的問(wèn)題，比如文件的打開(kāi)，文件按鈕虛函數(shù)的添加，繪圖無(wú)效等。在一步一步解決這些問(wèn)題的同時(shí)，漸漸進(jìn)一步的熟悉掌握了譜數(shù)據(jù)處理的方法和編程思想。學(xué)生（簽名）：袁子程 2015年 1月10 日誠(chéng)信承諾本人鄭重聲明所呈交的課程報(bào)告是本人在指導(dǎo)教師指導(dǎo)下進(jìn)行的研究工作及取得的研究成果。據(jù)我所知，除了文中特別加以標(biāo)注的地方外，論文中不包含其他人已經(jīng)發(fā)表或撰寫(xiě)過(guò)的研究成果。與我一同工作的同學(xué)對(duì)本文研究所做的貢獻(xiàn)均已在報(bào)告中作了明確的說(shuō)明并表