3 谱数据的定性分析 核数据处理课件.ppt

东华理工大学 张怀强 寻峰 简单比较法 1)满足条件： 可认为有峰存在 2)然后在data i-m至data i+m中找最大值，对应的道 值即为峰位 k：找峰阈值，根据高斯分布，一般k取值1—1.5 一般，用 R=N0 / Nb ≥ R0确定峰是否有意义 R为峰谷比， R0为设定值 （经验值） N0为净峰幅度与基底之和 Nb为基底计数 东华理工大学 张怀强 寻峰 简单比较法 3)峰的左右边界道i-L和i+R的确定 当峰形出现 左、右 “拖尾” 时，应限制 L、 R 的长度——半高宽 东华理工大学 张怀强 寻峰 简单比较法 方法简单，只能用于简单明显的峰，不能用于寻找重峰 如果寻峰宽度因子小，即m小，寻峰灵敏度高，但可能会把统计涨落引起的假峰也作为真峰，所以，用R=N0 / Nb ≥ R0确定峰是否有意义 如果寻峰宽度因子大，即m大，寻峰灵敏度低，可能丢掉弱峰 常用的5点、7点极大值法。即 m 取 2 、3 东华理工大学 张怀强 寻峰 简单比较法——例子 请用五点简单比较法判断该段谱数据是否有峰存在？若有，请求出该峰的峰位和左、右边界的道址。 道址 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 m=2，k=1 从左向右依次计算判断 L、R:从峰位开始分别向左、右依次计算判断 东华理工大学 张怀强 寻峰 简单比较法 编程——寻峰位 int CMmcaView::SearPeakCampare(int Beginch, int Endch, int m, float k) { int n1=0, Maxtemp; for (int j = Beginch; j = Endch; j++) { double f = data[j] - k * sqrt(double(data[j])); if ((data[j-m] f) || (data[j+m] f)) continue ; n1=j-m; Maxtemp=data[j-m]; for (int i=j-m; i=j+m; i++) { if (data[i]Maxtemp) { Maxtemp=data[i]; n1=i; } } } return n1; } 找最大值 东华理工大学 张怀强 寻峰 简单比较法 编程——寻峰边界——左边界 int CMmcaView::SearLBorderCampare(int Beginch, int peak, int m, float k) { int n1=beginch; for (int j = Peak; j = Beginch; j--) { double f = data[j] + k * sqrt(double(data[j])); if ((data[j-m] = f))) { n1=j; break; } } return n1; } 东华理工大学 张怀强 寻峰 简单比较法 编程——寻峰边界——右边界 int CMmcaView::SearRBorderCampare(int Endch, int peak, int m, float k) { int n1=Endch; for (int j = Peak; j = Endch; j++) { double f = data[j] + k * sqrt(double(data[j])); if ((data[j+m] = f))) { n1=j; break; } } return n1; } 东华理工大学 张怀强 寻峰 高斯乘积函数找峰法 基本思路 描述谱峰形状的函数主要是高斯函数： 则由相邻的数据点定义一个新的函数（第一高斯乘积函数，只与H有关）： 理论上峰的判断： 东华理工大学 张怀强 寻峰 高斯乘积函数找峰法 基本思路 同理，定义一个更一般的高斯乘积函数： 其中，m是步长（用道表示），是高斯乘积函数的阶数，则Pm(i)