[发明专利]基于一阶导数特征的峰识别算法

说明书

技术领域

本发明涉及一种信号处理中峰形的检测和识别算法，属于信号处理技术领域。

背景技术

在某些信号处理中，需要检测信号中出现的波峰，并识别出峰的起始点、峰点、结束点。如色谱峰的分割识别。

现有的峰识别算法多通过一阶和二阶导数进行识别。单独使用一阶导数的方法通过设定阈值，将一阶导数超过阈值的点作为峰的起始点。这种方法无法适应基线漂移的情况，对于重叠峰也较难处理。同时使用一阶和二阶导数的方法由于要计算二阶导数使得计算量增加，处理重叠峰的方法也较复杂。

另外还有些峰识别的方法通过高斯波进行匹配，得到峰的位置信息。这种方法只能识别高斯峰，对于峰的起始点和终点位置判断也不够精确。

发明内容

本发明所要解决的技术问题是克服现有技术中的缺陷，提供一种基于一阶导数特征的峰识别算法，基线漂移自适应，无需设置导数阈值，可以准确分割各种重叠峰，也无需计算复杂的二阶导数。

为解决上述技术问题，本发明提供一种基于一阶导数特征的峰识别算法，其特征是，包含以下步骤：

(1)定义一阶导数特征点：分别为局部最大值点、局部最小值点、零点；

(2)基于一阶导数特征点定义完整峰、左肩峰、右肩峰：

完整峰：包含连续的三个一阶导数特征点：局部最大值点、零点、局部最小值点，

左肩峰：包含连续的两个一阶导数特征点：局部最大值点、局部最小值点，并且局部最大值点的导数绝对值大于局部最小值点的导数，

右肩峰：包含连续的两个一阶导数特征点：局部最大值点、局部最小值点，并且局部最大值点的导数绝对值小于局部最小值点的导数；

(3)检测识别出所述完整峰、左肩峰、右肩峰的算法步骤。

检测识别出所述完整峰、左肩峰、右肩峰的算法采用搜索特征点匹配的方式定位峰，完整峰搜索定位起始点、结束点，左肩峰搜索定位起始点，右肩峰搜索定位结束点。

设原始数据为y＝f(x)，其一阶导数为y′＝f(x)，定义三种一阶导数特征点：

1.一阶导数局部最大值点

设x＝x₀，ε＞0，在x₀的邻域d＝(x₀-ε，x₀+ε)中有：f′(x₀)≥f′(x)，x∈d，则称x₀是一阶导数的局部最大值点，记为x₀∈DMAX。DMAX是所有一阶导数局部最大值点的集合。

2.一阶导数局部最小值点

设x＝x₀，ε＞0，在x₀的邻域d＝(x₀-ε，x₀+ε)中有：f′(x₀)≤f′(x)，x∈d，则称x₀是一阶导数的局部最小值点，记为x₀∈DMIN。DMIN是所有一阶导数局部最小值点的集合。

3.一阶导数零点

设x＝x₀，f′(x₀)＝0，则称x0是一阶导数零点，记为x₀∈DZERO。DZERO是所有一阶导数零点的集合。

定义三种峰类型：

1.完整峰

包含上升沿和下降沿的峰。其一阶导数必包含连续的三个特征点：局部最大值点、零点、局部最小值点。

2.左肩峰

只含有上升沿的峰。其一阶导数必包含连续的两个特征点：局部最大值点、局部最小值点，并且局部最大值点的导数绝对值大于局部最小值点的导数。

3.右肩峰

只含有下降沿的峰。其一阶导数必包含连续的两个特征点：局部最大值点、局部最小值点，并且局部最大值点的导数绝对值小于局部最小值点的导数。

检测识别所述完整峰、左肩峰、右肩峰的算法步骤为：

1.输入原始数据，进行滤波去噪，计算一阶导数。

2.根据一阶导数特征点的定义得到所有的一阶导数特征点c_i，i＝1..n，n是特征点个数。

3.在一阶导数特征点中搜索峰：