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

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

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

3、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); 零面積對稱法：基本思想：面積為零的“窗”函數(shù)與實(shí)驗(yàn)譜數(shù)據(jù)進(jìn)行褶積變換，且要求“窗”函數(shù)為對稱函數(shù)。數(shù)學(xué)表達(dá)式：其中，為變換后的譜數(shù)據(jù)，dat

4、ai+j為原始實(shí)驗(yàn)譜數(shù)據(jù)，Cj為對稱零面積變換函數(shù)，W=2m+1為窗寬（變換寬度）。對稱零面積變換函數(shù)：Gj：實(shí)際代碼：/零面積對稱法尋峰(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-left

5、;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.017584*i-0.138133;for( j=0;j<1024;j+)if(Edaoj<(6.403-6.403*0.1)|Edaoj>(6.403+6.403*0.1)chdaoj="Fe&

6、quot;else 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+1

7、0.5*0.1)chdaoj="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)行五、 體會在本次程序設(shè)計(jì)實(shí)習(xí)中，我再次學(xué)習(xí)了MFC編程方法。過程中，我遇到了許多曾經(jīng)都沒有遇到的問題，比如文件的打開，文件按鈕虛函數(shù)的添加，繪圖無效等。在一步一步解決這些問題的同時

8、，漸漸進(jìn)一步的熟悉掌握了譜數(shù)據(jù)處理的方法和編程思想。但是依舊有一個問題尚未解決，那就是數(shù)據(jù)類型的變化。這個使得我無法在數(shù)據(jù)處理中正確地獲取應(yīng)有的結(jié)果。令我費(fèi)解的是成譜圖像卻是與原始數(shù)據(jù)非常的吻合。對于這個問題我依然百思不得其解。六、 建議我對本門課程的掌握程度也非常有限，提出以下建議供老師和學(xué)弟學(xué)妹們參考。1.MFC編程時一個系統(tǒng)過程，需要掌握的細(xì)節(jié)可以隨時網(wǎng)上學(xué)習(xí)，但是總體思路要牢記，不容易臨時掌握。2.對于本專業(yè)同學(xué)們而言，這次實(shí)習(xí)是擴(kuò)展我專業(yè)同學(xué)技能的大好機(jī)會，建議不要輕易放棄，認(rèn)真編寫，即使最后沒有達(dá)到自己預(yù)期的計(jì)劃也是收獲頗豐的。3.在編寫的過程中要有良好的格式及注釋習(xí)慣，以免一段時

9、間甚至幾分鐘后就看不懂自己的思路。4.在編寫算法過程中，一定要回歸理論，弄明白實(shí)際數(shù)學(xué)過程，對算法的編寫很有幫助。5.希望老師能夠?qū)Ψ撬惴ú糠纸o予同學(xué)更多的幫助，教授同學(xué)建立軟件的過程，指導(dǎo)同學(xué)完成非專業(yè)的部分，讓同學(xué)著眼于數(shù)學(xué)導(dǎo)算法的翻譯。最后，衷心感謝在本次程序設(shè)計(jì)實(shí)習(xí)過程中，馬永紅老師細(xì)心的指導(dǎo)，您不厭其煩地解答我的各種疑問，并且?guī)椭覍ふ以紨?shù)據(jù)的數(shù)據(jù)格式和數(shù)據(jù)類型。我還要感謝蘆瑞琪、蔣卓辰同學(xué)對我在思路上的重要啟發(fā)，謝謝你們給予我的大力支持！附：程序包附在附件中，項(xiàng)目名稱為37gai學(xué)生學(xué)習(xí)心得在本次程序設(shè)計(jì)實(shí)習(xí)中，我再次學(xué)習(xí)了MFC編程方法。過程中，我遇到了許多曾經(jīng)都沒有遇到的問題，比如文件的打開，文件按鈕虛函數(shù)的添加，繪圖無效等。在一步一步解決這些問題的同時，漸漸進(jìn)一步的熟悉掌握了譜數(shù)據(jù)處理的方法和編程思想。學(xué)生（簽名）：袁子程 2015年 1月10 日誠信承諾本人鄭重聲明所呈交的課