基于STM32的心率计（二）动态阈值算法获取心率值

[导读]前言上一篇文章：基于uFUN开发板的心率计（一）DMA方式获取传感器数据，介绍了如何获取PulseSensor心率传感器的电压值，并对硬件电路进行了计算分析。心率计，重要的是要获取到心率值，本篇文章将介绍一种采样数据处理算法——动态阈值算法，来获取心率值，这种算法来自于一位网友：...

前言

上一篇文章：基于uFUN开发板的心率计（一）DMA方式获取传感器数据，介绍了如何获取PulseSensor心率传感器的电压值，并对硬件电路进行了计算分析。心率计，重要的是要获取到心率值，本篇文章将介绍一种采样数据处理算法——动态阈值算法，来获取心率值，这种算法来自于一位网友：玩的就是心跳 —— 使用 PulseSensor 脉搏传感器测量心率(http://www.shaoguoji.cn/2017/06/29/use-pulsesensor-with-stm32/)，本文部分内容摘自这篇文章。

IBI和BPM

心率，指的是一分钟内的心跳次数，得到心率最笨的方法就是计时一分钟后数有多少次脉搏。但这样的话每次测心率都要等上个一分钟才有一次结果，效率极低。另外一种方法是，测量相邻两次脉搏的时间间隔，再用一分钟除以这个间隔得出心率。这样的好处是可以实时计算脉搏，效率高。
IBI：相邻两次脉搏的时间间隔，单位：s。
BPM：心率，一分钟内的心跳次数。即 BPM=6000/IBI

例如，在这张心率传感器输出信号的波形图中，可以计算出，两次波峰之间的时间为：0.685s，心率值为：60/0.685 = 87。

从网上找来的 arduino 开源算法复杂的一匹，看了一遍感觉一头雾水（反正我暂时没看懂）。由上面的分析可以得出，我们的最终目的就是要求出 IBI 的值，并通过 IBI 计算出实时心率。既然知道原理了那就自己来把算法实现吧。有兴趣研究官方算法的朋友，可以下载：PulseSensor资料.rar。

核心操作 —— 识别一个脉搏信号

无论是采用计数法还是计时法，只有能识别出一个脉搏，才能数出一分钟内脉搏数或者计算两个相邻脉搏之间的时间间隔。那怎么从采集的电压波形数据判断是不是一个有效的脉搏呢？

显然，可以通过检测波峰来识别脉搏。最简单粗暴的方法是设定一个阈值，当读取到的信号值大于此阈值时便认为检测一个脉搏。似乎用一个 if 语句就轻轻松松解决。但，事情真的有那么简单么？

其实这里存在两个问题。

问题一：阈值的选取

作为判断的参考标尺，阈值该选多大？10？100？还是1000？我们不得而知，因为波形的电压范围是不确定的，振幅有大有小并且会改变，根本不能用一个写死的值去判断。就像下面这张图一样：

可以看出，两个形状相同波形的检测结果截然不同 —— 同样是波峰，在不同振幅的波形中与阈值比较的结果存在差异。实际情况正是如此：传感器输出波形的振幅是在不断随机变化的，想用一个固定的值去判定波峰是不现实的。

既然固定阈值的方法不可取，那自然想到改变阈值 —— 根据信号振幅调整阈值，以适应不同信号的波峰检测。通过对一个周期内的信号多次采样，得出信号的最高与最低电压值，由此算出阈值，再用这个阈值对采集的电压值进行判定，考虑是否为波峰。也就是说电压信号的处理分两步，首先动态计算出参考阈值，然后用用阈值对信号判定、识别一个波峰。