c语言寻峰算法,寻峰算法总结报告

闭合重峰、波峰定位、寻峰算法报告

自动寻峰

由于谱结构的复杂和统计涨落的影响，从谱中正确地找到全部存在的峰是比较困难的。尤其是找到位于很高本底上的弱峰，分辨出相互靠得很近的重峰更为困难。

谱分析对寻峰方法的基本要求如下：

(1) 比较高的重峰分辨能力。能确定相互距离很近的峰的峰位。

(2) 能识别弱峰，特别是位于高本底上的弱峰。

(3) 假峰出现的几率要小。

(4) 不仅能计算出峰位的整数道址，还能计算出峰位的精确值，某些情况下要求峰位的误差小于0.2道。

很多作者对寻峰方法进行了研究，提出了很多有效的寻峰方法。

寻峰的目的 ：

判断有没有峰存在

确定峰位(高斯分布的数学期望)，以便把峰位对应的道址，转换成能量

确定峰边界——为计算峰面积服务(峰边界道的确定，直接影响峰面积的计算)

分为两个步骤：谱变换和峰判定

要求：支持手动/自动寻峰，参数输入，同时计算并显示峰半高宽、精确峰位、峰宽等信息，能够区分康普顿边沿和假峰

感兴区内寻峰

人工设置感兴趣大小，然后在感兴区内采用简单方法寻峰

重点研究：对感兴区内的弱峰寻峰、重峰的分解

对于一个单峰区，当峰形在峰位两侧比较对称时，可以由峰的FWHM计算峰区的左、右边界道址。峰区的宽度取为3FWHM，FWHM的值可以根据峰位mp由测量系统的FWHM刻度公式计算。由于峰形对称，左、右边界道和峰位的距离都是1.5FWHNM。

式中mp是峰位，INT的含义是取整数。

对于存在有低能尾部的峰，其峰形函数描述(参见图)。

，m≥mp－J

，m≤mp－J

式中H为峰高，mp为峰位，是高斯函数的标准偏差，J为接点的道址和峰位之间的距离。在峰位的左侧，有一个接点，其道址为mp－J。在接点的右侧，峰函数是高斯函数。在接点的左侧，峰函数用指数曲线来描述。这时峰区的左、右边界道址为

带有低能尾部的峰函数的图形

全谱自动寻峰

基于核素库法：能量刻度完成后，根据核素库中的能量计算对应的道址，在各个道址附近(左右10道附近)采用简单的寻峰方法(导数法)

方法：

根据仪器选择开发

IF函数法/简单比较法(适于寻找强单峰，速度快)

满足条件： 可认为有峰存在

然后在data i-m至data i+m中找最大值，对应的道值即为峰位

k：找峰阈值，根据高斯分布，