频谱数据的寻峰.doc

频谱数据的寻峰
寻峰(Peaksearch)现在还没有一种公认的适用于各种情况的寻峰算法，而人们的眼睛才是最好的寻峰工具。常用的寻峰算法有比较法、导数法、插值法、卷积法、函数拟合法等等。
1.比较法
比较法的原理是根据峰的直观定义，找出比周围数值大得多的（至少大于设定的平均值）的数据值，从而确定峰值。原理上本项目中频谱数据的峰值算法可以采用比较法，即找出所有数据中的最大值，然后确定峰值，但是项目中FFT之后的频谱数据复杂，简单的比较法是无法找到最优值。比较法的优点是程序简单，运算速度快，适用于寻找孤立强峰；缺点是在复杂谱中无法寻找弱峰。
2.导数法
导数法的基本思想是把平滑后的谱线看成是一条连续曲线。通过对谱线上个点求导数，然后依据导数的性质确定峰值。导数法又分为一阶、二阶、三阶三种方法。。以S代表各地址指针数据的连续函数,各阶导数的函数图形如图1所示。

图1各阶导数函数图形
a.一阶导数法，根据图1在峰值处，导数值为零，且随着地址的增加，导数的值由正变负。一阶导数法对孤立的小峰灵敏度高，但是不能分辨重峰。
b.二阶导数法，利用二阶导数寻峰是比较常用的一种寻峰方式。它的工作原理是在峰位置上二阶导数出现极小值。二阶导数寻峰有较好的灵敏度，且对重峰有一定的分辨能力。
c.三阶导数法，在峰值处，三阶导数为零，且随着地址指针的增加，导数的值由负变正。
3.插值法
插值法是人们工程实践与科学研究中经常用到的一种方法。本项目中的FFT之后的心跳频谱和呼吸频谱的寻峰就是采用频谱插值法。常用的插值方法又有拉格朗日插值法、牛顿插值法、样条插值和Hermite插值法。项目中的插值法应该是牛顿一次插值法，该算法编程简单，便于实现。算法中用到的公式：
i=i±x(i+1)xi+1+x(i)
其中，x()FFT之后频谱数据，i为它的地址指针。