数据处理与能谱分析.doc

数据处理与能谱分析随着计算机技术的发展，利用计算机处理实验数据也越来越常见，随之如来就产生了许多的软件如：Matlab 、 Excel、CAD等等， 但一般这些软件在处理物质放射性衰变时都比较繁琐，因此在处理放射性物质衰变时的数据时，就必须自己依照其规律制作数据处理软件来探究该物质的各种性质，从而来确定该物质的类型以及其运用。我们要使用C+平台和C语言来制作软件及编写对应的响应函数，在进行实验时我们把我们利用专业仪器测得的数据按照要求记录并保存在一个TXT文件中，使用C语言来编写程序对文件里的数据进行读写和操作；例如：调用文件打开函数(fopen)和关闭函数(fclose)语句来对文件进行打开和关闭处理，调用读写语句对文件内的数据进行察看和编写fscanf(fp,”%d %d”,&I,&t)、fprintf(fp,”%d %d”,I,t)。制作EXE软件来对数据进行图谱显示：在相应的操作界面放置相应的显示框图和对应的按钮，如下图然后给框图和按钮赋地址，并添加相应的变量和相应函数。最后组建编译运行。软件运行后先点击数据读写按钮在点击原始图谱按钮最后点击五点平滑按钮就得到了该物质数据图谱。如下图所示：注：横坐标为道址、纵坐标为计数率 添加相应的数据输出框和按钮如：峰值、道值、面积。最后再添加相应的响应函数，再编译运行。得到如下图所示的图谱：至此我们利用C语言对数据处理和能谱分析已经结束，最后我们在根据图谱的形状、峰值、道值以及各个部分的面积来确定物质衰变的性质和物质本身的性质，最终来确定放射性物质的用途。附件：部分响应函数数据读写程序：void CShiyanDlg:OnReadfile() / TODO: Add your control notification handler code hereFILE *fp;int datanum=0; int i;int data1,data2;if(fp=fopen(090623.txt,r)=NULL) printf(Cannot open the file.n); exit(0); while(!feof(fp) fscanf(fp,%d %d,&data1,&data2);datadatanum+=data2;fscanf(fp,n); for(i=0;i2048;i+) if( i2045) data_phi=datai; else data_phi=(datai-2+4*datai-1+6*datai+4*datai+1+datai+2)*1.0/16.0; printf(%d %fn,i+1,data_phi); fclose(fp); if(fp=fopen(out.txt,w)=NULL) printf(file open error.n); exit(0);for(i=0;iGetDC(); pDC-SetMapMode(MM_ANISOTROPIC); CRect rect; double scale; GetDlgItem(IDC_P1)-GetWindowRect(&rect); scale=25.4*10/pDC-GetDeviceCaps(LOGPIXELSX);CDC MemDC; CBitmap MemBitmap; MemDC.CreateCompatibleDC(NULL); MemBitmap.CreateCompatibleBitmap(pDC,rect.Width(),rect.Height(); MemDC.FillSolidRect(0,0,rect.Width(),rect.Height(),RGB(255,255,255); max=data0;for(int i=1;imax ) max=datai; for ( i = 0;iSetPixel(xViewport,rect.Height()-yViewport-1,RGB(255,0,0); pDC-BitBlt(0,0,rect.Width(),rect.Height(),&MemDC,0,0,SRCCOPY); MemBitmap.DeleteObject(); MemDC.DeleteDC(); pDC-DeleteDC(); 五点平滑图谱响应函数：void CShiyanDlg:OnWudianpinghua() / TODO: Add your control notification handler code heredouble xViewport,yViewport;float max;CWnd *pWnd=GetDlgItem(IDC_P2); CDC *pDC =pWnd-GetDC(); pDC-SetMapMode(MM_ANISOTROPIC); CRect rect; double scale; GetDlgItem(IDC_P2)-GetWindowRect(&rect); scale=25.4*10/pDC-GetDeviceCaps(LOGPIXELSX);CDC MemDC; CBitmap MemBitmap; MemDC.CreateCompatibleDC(NULL); MemBitmap.CreateCompatibleBitmap(pDC,rect.Width(),rect.Height(); MemDC.FillSolidRect(0,0,rect.Width(),rect.Height(),RGB(255,255,255); max=data_ph0;for(int i=1;imax ) max=data_phi; for ( i = 0;iSetPixel(xViewport,rect.Height()-yViewport-1,RGB(255,0,0); pDC-BitBlt(0,0,rect.Width(),rect.Height(),&MemDC,0,0,SRCCOPY); MemBitmap.DeleteObject(); MemDC.DeleteDC(); pDC-DeleteDC(); 道值寻峰函数： void CShiyanDlg:OnLButtonUp(UINT nFlags, CPoint point) / TODO: Add your message handler code here and/or call defaultCRect rect1;CWnd *pWnd=GetDlgItem(IDC_P1);pWnd-GetWindowRect(rect1);ScreenToClient(&rect1);CDC *pDC =pWnd-GetDC(); if (rect1.PtInRect(point) pDC-MoveTo(point.x-rect1.left,0); pDC-LineTo(point.x-rect1.left,rect1.Height();if(LButton_Clicked=0)CPen lpen(PS_SOLID, 1, RGB(255,255,255);CPen* pOldPen = pDC-SelectObject(&lpen); pDC-MoveTo(x_old-rect1.left,0); pDC-LineTo(x_old-rect1.left,rect1.Height();lpen.DeleteObject(); x_old=point.x;y_old=point.y;elseLButton_Clicked=0;x_old=point.x;y_old=point.y;CDialog:OnLButtonUp(nFlags, point);void CShiyanDlg:OnButton1() / TODO: Add your control notification handler code hereLButton_Clicked=1;left_flag=1;Left_x=x_old;Left_y=y_old;void CShiyanDlg:OnButton2() / TODO: Add your control notification handler code hereint i,j,xl,xr,max=0,maxi=0;float newx;CRect rect1;CWnd *pWnd=GetDlgItem(IDC_P1);pWnd-GetWindowRect(rect1);ScreenToClient(&rect1);CDC *pDC =pWnd-GetDC(); RButton_Clicked=1;if(left_flag=1)Right_x=x_old; Right_y=y_old; /GetDlgItem(IDC_P1)-GetWindowRect(&rect1); xl=2048*(Left_x-rect1.left)/rect1.Width(); xr=2048*(Right_x-rect1.left)/rect1.Width(); for(i=xl+2;idatai-2) & (datai-1.0*sqrt(datai)datai+2)for(j=i-2;jmax) max=dataj; maxi=j; m_e3=maxi; UpdateData(false);newx=maxi*rect1.Width()/2048;CDC *pDC =pWnd-GetDC(); CPen lpen(PS_SOLID,1,RGB(0,0,255);CPen* pOldPen = pDC-SelectObject(&lpen); /scale=25.4*10/pDC-GetDeviceCaps(LOGPIXELSX);pDC-MoveTo(newx,0); pDC-LineTo(newx,rect1.Height();lpen.DeleteObject(); RButton_Clicked=0;LButton_Clicked=0;void CShiyanDlg:OnZuodaozhi() / TODO: Add your control notification handler code hereCRect rect1;int chx; CWnd *pWnd=GetDlgItem(IDC_P1);pWnd-GetWindowRect(rect1);ScreenToClient(&rect1);chx=2048*(x_old-rect1.left)/rect1.Width();m_e1=chx;UpdateData(false);void CShiyanDlg:OnYoud