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

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

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

3、b=2m+1，為平滑寬度。譜光滑57911Kb3521231429A0177179143A±1126135120A±2-333060A±3-12-55-10A±415-45A±518實(shí)際代碼：/多項(xiàng)式最小二乘法譜光滑（參數(shù)為5） for( i=2;i<1023;i+)datai=1/35*(-3)*datai-2+12*datai-1+17*datai+12*datai+1+(-3)*datai+2); 零面積對(duì)稱法：基本思想：面積為零的“窗”函數(shù)與實(shí)驗(yàn)譜數(shù)據(jù)進(jìn)行褶積變換，且要求“窗”函數(shù)為對(duì)稱函數(shù)。數(shù)學(xué)表達(dá)式：其中，為變換后的譜數(shù)據(jù)，d

4、atai+j為原始實(shí)驗(yàn)譜數(shù)據(jù)，Cj為對(duì)稱零面積變換函數(shù)，W=2m+1為窗寬（變換寬度）。對(duì)稱零面積變換函數(shù)：Gj：實(shí)際代碼：/零面積對(duì)稱法尋峰(m=16,W=11,H=4)for( i=16;i<1007;i+)for(j=0;j<16;j+)Gi+=Gj;d=1/11*Gi;Ci=Gi-d;for( i=16;i<1007;i+)if(yi/dy)>f)fengk+=datai;Wasson（瓦森、沃森）峰面積法： 計(jì)算總面積：計(jì)算本底面積： 計(jì)算凈峰面積：實(shí)際代碼：/ Wasson（瓦森、沃森）峰面積法：for (i=L;i<=R;i+)RL=right-le

5、ft;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ù)能量刻度系數(shù)求峰位能量/*for( i=0;i<1023;i+)Ei=0.*i-0.;for( j=0;j<1024;j+)if(Edaoj<(6.403-6.403*0.1)|Edaoj>(6.403+6.403*0.1)chdaoj="Fe"else

6、if(Edaoj<(7.477-7.477*0.1)|Edaoj>(7.477+7.477*0.1)chdaoj="Ni"else if(Edaoj<(8.047-8.047*0.1)|Edaoj>(8.047+8.047*0.1)chdaoj="Cu"else if(Edaoj<(8.638-8.638*0.1)|Edaoj>(8.638+8.638*0.1)chdaoj="Zn"else if(Edaoj<(10.5-10.5*0.1)|Edaoj>(10.5+10.5*0.1)ch

7、daoj="Ap"else if(Edaoj<(11.725-11.725*0.1)|Edaoj>(11.725+11.725*0.1)chdaoj="As"else if(Edaoj<(12.611-12.611*0.1)|Edaoj>(12.611+12.611*0.1)chdaoj="Ni"*/8. 調(diào)試運(yùn)行五、 體會(huì)在本次程序設(shè)計(jì)實(shí)習(xí)中，我再次學(xué)習(xí)了MFC編程方法。過(guò)程中，我遇到了許多曾經(jīng)都沒(méi)有遇到的問(wèn)題，比如文件的打開(kāi)，文件按鈕虛函數(shù)的添加，繪圖無(wú)效等。在一步一步解決這些問(wèn)題的同時(shí)，漸漸進(jìn)一步的熟悉掌

8、握了譜數(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è)同學(xué)技能的大好機(jī)會(huì)，建議不要輕易放棄，認(rèn)真編寫(xiě)，即使最后沒(méi)有達(dá)到自己預(yù)期的計(jì)劃也是收獲頗豐的。3.在編寫(xiě)的過(guò)程中要有良好的格式及注釋習(xí)慣，以免一段時(shí)間甚至幾分鐘后就看不

9、懂自己的思路。4.在編寫(xiě)算法過(guò)程中，一定要回歸理論，弄明白實(shí)際數(shù)學(xué)過(guò)程，對(duì)算法的編寫(xiě)很有幫助。5.希望老師能夠?qū)Ψ撬惴ú糠纸o予同學(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ì)我在思路上的重要啟發(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)告是本人在指