峰面积计算.pptVIP

成都理工大学 马英杰 峰面积的计算 单峰曲线拟合法 单峰拟合函数：data(i)=G(i)+B(i) 用峰边界内的点的计数值作为一组已知的数据，代入拟合函数，用非线性最小二乘法，求解各系数（a1, a2, a3, b1, b2, b3） 则峰面积为： a1为峰的高度，a3为标准误差 成都理工大学 马英杰 峰面积的计算 意义 峰面积的计算是定量分析的基础。 知道了特征峰的净峰面积，就可以计算目标元素的含量。 实测谱中，各特征峰是叠加在环境本底和康普顿散射背景之上的。 总面积S：在一个指定的峰区内，各道计数之和 本底面积B：由环境本底和散射造成的计数总和 净峰面积A：由峰的总面积扣除本底面积即可得出净峰面积 即：峰的总面积—本底面积=净峰面积 S - B = A，所以，计算净峰面积，如何确定B最关键！ 成都理工大学 马英杰 峰面积的计算 A = S – B，关键是如何确定B？ 方法： 线性本底法（总峰面积法，TPA法） Covell（科沃尔）峰面积法 Wasson（瓦森、沃森）峰面积法 Sterlinski(斯托林斯基)峰面积法 平均总峰面积法 单峰曲线拟合法 成都理工大学 马英杰 L R yL yR 峰面积的计算 线性本底法（总峰面积法，TPA法） 确定本底面积计算方法： 左右边界点直线连接即为本底线 线性本底 梯形法计算本底面积 B 成都理工大学 马英杰 峰面积的计算 线性本底法（总峰面积法，TPA法） 1）确定峰的左、右边界L、R 2）计算总面积： 3）计算本底面积： 4）计算净峰面积： L B R yL yR 成都理工大学 马英杰 峰面积的计算 线性本底法——例子 按给定的左右边界道址，用全峰面积法计算该峰面积值 道址 227 228 229 230 231 232 233 234 235 计数 506 373 410 400 481 554 620 763 922 道址 236 237 238 239 240 241 242 243 244 计数 1120 1539 1955 2412 2979 3267 3082 2847 2256 道址 245 246 247 248 249 250 251 252 253 计数 1648 1031 622 343 212 145 92 102 104 成都理工大学 马英杰 峰面积的计算 线性本底法（总峰面积法，TPA法） 编程： int 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; } 成都理工大学 马英杰 峰面积的计算 线性本底法（总峰面积法，TPA法） 编程： int 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; } 计算本底面积 成都理工大学 马英杰 峰面积的计算 线性本底法（总峰面积法，TPA法） 编程： int 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; } 计算总峰面积 成都理工大学 马英杰 峰面积的计算 线性本底法（总峰面积法，TPA法） 编程： int CMmcaView::CalculatePeakArea(int L, int R, ) { int bdarea, zarea=0, area; bdarea = (Data[L] + Data[R]) *