峰面积计算.ppt

成都理工大学马英杰峰面积的计算峰面积的计算?意义?峰面积的计算是定量分析的基础。?知道了特征峰的净峰面积,就可以计算目标元素的含量。?实测谱中,各特征峰是叠加在环境本底和康普顿散射背景之上的。?总面积S:在一个指定的峰区内,各道计数之和?本底面积B:由环境本底和散射造成的计数总和?净峰面积A:由峰的总面积扣除本底面积即可得出净峰面积即:峰的总面积—本底面积=净峰面积S - B = A,所以,计算净峰面积,如何确定B最关键!1成都理工大学马英杰峰面积的计算峰面积的计算?A = S – B,关键是如何确定B??方法:?线性本底法(总峰面积法,TPA法)?Covell(科沃尔)峰面积法?Wasson(瓦森、沃森)峰面积法?Sterlinski(斯托林斯基)峰面积法?平均总峰面积法?单峰曲线拟合法2成都理工大学马英杰LRyLyR峰面积的计算峰面积的计算?线性本底法(总峰面积法,TPA法)?确定本底面积计算方法:?左右边界点直线连接即为本底线?线性本底?梯形法计算本底面积2)1(*)(2)1(*)(????????LRyyLRdatadataBRLRLB3成都理工大学马英杰峰面积的计算峰面积的计算?线性本底法(总峰面积法,TPA法)?1)确定峰的左、右边界L、R?2)计算总面积:?3)计算本底面积:?4)计算净峰面积:???RLiidataSBSA??2)1(*)(????LRdatadataBRLLBRyLyR4成都理工大学马英杰峰面积的计算峰面积的计算?线性本底法——例子?按给定的左右边界道址,用全峰面积法计算该峰面积值道址227228229230231232233234235计数506373410400481554620763922道址236237238239240241242243244计数112015391955241229793267308228472256道址245246247248249250251252253计数16481031622343212145921021042970092145...400410?????????RLiidataS57732/23*)92410(2)1(*)(???????LRdatadataBRL23927577329700???A5成都理工大学马英杰峰面积的计算峰面积的计算?线性本底法(总峰面积法,TPA法)?编程:???RLiidataSBSA??2)1(*)(????LRdatadataBRLint CMmcaView::CalculatePeakArea(int L, int R ){int bdarea, zarea=0, area;bdarea = (Data[L] + Data[R]) * (R - L + 1) / 2;//计算总面积for (int i = L; i<= R; i++)zarea = Data[i] + zarea;//计算净峰面积area = zarea - bdarea;return area;}6成都理工大学马英杰峰面积的计算峰面积的计算?线性本底法(总峰面积法,TPA法)?编程:???RLiidataSBSA??2)1(*)(????LRdatadataBRLint CMmcaView::CalculatePeakArea(int L, int R, ){int bdarea, zarea=0, area;bdarea = (Data[L] + Data[R]) * (R - L + 1) / 2;//计算总面积for (int i = L; i<= R; i++)zarea = Data[i] + zarea;//计算净峰面积area = zarea - bdarea;return area;}计算本底面积7成都理工大学马英杰峰面积的计算峰面积的计算?线性本底法(总峰面积法,TPA法)?编程:???RLiidataSBSA??2)1(*)(????LRdatadataBRLint CMmcaView::CalculatePeakArea(int L, int R, ){int bdarea, zarea=0, area;bdarea = (Data[L] + Data[R]) * (R - L + 1) / 2;//计算总面积for (int i = L; i<= R; i++)zarea = Data[i] + zarea;//计算净峰面积area = zarea - bdarea;return area;}计算总峰面积8成都理工大学马英杰峰面积的计算峰面积的计算?线性本底法(总峰面积法,TPA法)?编程:???RLiidataSBSA??2)1(*)(????LRdatadataBRLint CMmcaView::CalculatePea