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

1、京#Wz大挈程序設(shè)計(jì) 課程設(shè)計(jì)（論文）設(shè)計(jì)（論文）題目譜的顯示、譜光滑、定量分析（峰面積）程序?qū)W院名稱 核技術(shù)與自動(dòng)化工程學(xué)院專業(yè)名稱 輻射防護(hù)與環(huán)境工程學(xué)生姓名袁子程學(xué)生學(xué)號(hào)2任課教師馬永紅設(shè)計(jì)（論文）成績(jī)教務(wù)處制一、實(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、7、9、11之間可選，討論最佳參數(shù)三、實(shí)習(xí)計(jì)劃1 .首先進(jìn)行MFC繪

2、圖區(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ù)處理的算法編寫5 .編寫成譜圖形功能并給出處理結(jié)論四、主要工作介紹（詳細(xì)介紹實(shí)習(xí)中各功能的實(shí)現(xiàn)方法、理論公式和計(jì)算原理，程序流程圖和程序運(yùn)行界面截圖等成果信息）1 .建立單文檔類MFC可執(zhí)行程序2 .建立繪圖區(qū)函數(shù)，初始化繪圖設(shè)備 -SI. I * |； T V « MuM"d +二產(chǎn)啊才.-l MH.工 * 口1f/ff "一 """""J MM. .打打 MF3NE

3、”¥Z： 7N3N*ff niirgpplrn irjBritH2叫.11匕1 HI i',phid iR " IHM - 4M1 Wafa- *吁 41*1 ari fW*11!M LlW-3G irh.1所-利啪jHliHihC"WWLd.* -J, *- EBrtE f-Kt； |旨"0IFLb«iIi 1 ray - md -tI/M ml .Irl I1： Li * - r MBVafi r Ml 木*r'i-i* + i3 .添加文件打開(kāi)虛函數(shù)New Virtual Override for cla« CM

4、y37giDocNew Virtual FunctionsExisting virtual function averridesCwnCloseFfamEOnNewDocumentDelete Contents G etFirstViewPo sition GetNextViewDnChs ngedViewList OnCloseDocument OnCmdMsg OnFinalRelease OnSaveDocument Rep ortS ave LoadExceplion SaveModilied SetPalhNameSctnticOnOnenDocumentSerialize確定Add

5、 HandlerAdd and EditEdit ExistingOnOpenDocumentQ: Called 君s part of the File Open command4,編寫讀取文件及初步尋峰代碼2 TIRE VtGI Z*E nt= :* *«£正 i 餐 餐工»¥ *-».*«* fJ *- »M 1 >» -*- - *- r!p>n ,»*111, I It 4ff fWC : Airi jwr tri al i frI n«wrid-H ri4f “f llij

6、Hdiiaiii嚕。M.'rypilTiaE13rqiilf 町/修.4*3 tail Mjiq*« « fi 3UEB也即F* / i百* * CMf TgiiW i.ilf -HI 明*FWA，i!F"» * C”-i，舊尸，k r- r k j n -她勒411H Ud 11 電 %fi MM .llE.“T I警也” Ie*科 .I kuT4.l-i|cim J «"l MU| ：）L，EJ7叼 I -HITrilWt4'i.'->«MJL3 j-!-U-«JI l-H.-J

7、'-£ir«l J-J!LIE 巾 F”1liik4 tl I4M d| I ：|"rf* a'|rij B ." I |if串 Im v-fa./C單i!鼻不主 卜尸口鼻/以3!rHkB'«Kil| .匕 I - i *l|tb .1I - «n “-D：“ H MWMl |i =4«, r4 vlr I .戶 H4Hmii n I箝I 1-ME"*I 口 LF開(kāi) 21 一打 t山/藝*工¥>«本柳*產(chǎn)T iLr p> /* si j1 lbI A>a

8、- pi i i .jt Cm i na薩 i i iMtmi >. MfM* » I*I ' / f/ . T-nnaif .肝r/n _5.添加繪制譜圖代碼6,嘗試讀取mca文件7.添加譜數(shù)據(jù)處理代碼多項(xiàng)式最小二乘擬合法：-1 m ,data, Ajdataj j , b=2m+1 ,為平滑范度。Kbj m譜光滑Kb3521231429A0177179143A ±1126135120A +2-333060A ±3-12-55-10A ±415-45A ±518實(shí)際代碼:多項(xiàng)式最小二乘法譜光滑（參數(shù)為 5）for( i=2;i&

9、lt;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á)式:myi Cjdatai jj mmCj 0 Cj C j j m其中，區(qū)為變換后的譜數(shù)據(jù),data+j為原始實(shí)驗(yàn)譜數(shù)據(jù)，Cj為對(duì)稱零面積變換函數(shù)，W=2m+1為窗寬（變換寬度）。對(duì)稱零面積變換函數(shù):CjmGjma）高斯函數(shù)：Gjb）柯西函數(shù)：Gj Gj:C）余弦平方函數(shù):d）雙曲正割函數(shù):exp_HFGj , Gj41n 2(j

10、. H)224j2 cos2( j. 2H)sech(2.634j H)實(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 (瓦森、沃森)峰面積法：n計(jì)算總面積：S datai0 i i n計(jì)算本底面積：dataR dataLb n r L(i0 L n) dataL R Lu dataR dataL /bn(i0 L n) dataLR

11、L2n 11B (bn bn)2- (b n bn)(n12)計(jì)算凈峰面積：A S B實(shí)際代碼：/Wasson (瓦森、沃森)峰面積法：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ù)能量刻度系數(shù)求峰位能量/* for( i=0;i<1023;i+)(Ei=0.017584*i-

12、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"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.63

13、8-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)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=&

14、quot;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)一步的熟悉掌握了譜數(shù)據(jù)處理的方法和編程思想。但是依舊有一個(gè)問(wèn)題尚未解決，那就是數(shù)據(jù)類型的變化。這個(gè) 使得我無(wú)法在數(shù)據(jù)處理中正確地獲取應(yīng)有的結(jié)果。令我費(fèi)解的是成 譜圖像卻是與原始數(shù)據(jù)非常的吻合。對(duì)于這個(gè)問(wèn)題我依然百思不得 其解。六、建議我對(duì)本門課程的掌握程度也非常有限，提出以下建議供老師和 學(xué)弟學(xué)妹們參考。1 .MFC編程時(shí)一個(gè)系統(tǒng)過(guò)程，需要掌握的細(xì)節(jié)可以隨時(shí)網(wǎng)上

15、學(xué) 習(xí)，但是總體思路要牢記，不容易臨時(shí)掌握。2 .對(duì)于本專業(yè)同學(xué)們而言，這次實(shí)習(xí)是擴(kuò)展我專業(yè)同學(xué)技能的大 好機(jī)會(huì)，建議不要輕易放棄，認(rèn)真編寫，即使最后沒(méi)有達(dá)到自己預(yù) 期的計(jì)劃也是收獲頗豐的。3 .在編寫的過(guò)程中要有良好的格式及注釋習(xí)慣， 以免一段時(shí)間甚 至幾分鐘后就看不懂自己的思路。4 .在編寫算法過(guò)程中，一定要回歸理論，弄明白實(shí)際數(shù)學(xué)過(guò)程,對(duì)算法的編寫很有幫助。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)題，比如文件