基于人体微动信号的呼吸率和心率的提取方法

技术领域

本发明属于人体生物测量技术领域，尤其涉及一种基于人体微动信号的呼吸率和心率的提取方法。

背景技术

心跳和呼吸是人体生命的基本体征，通过对心跳和呼吸的检测，特别是在安静或睡眠情境下的检测，可以探知身体活动状况，分析评估健康质量，提供改善健康的建议。

人体微动信号来源于，人在呼吸时胸腔的扩张、收缩运动以及心脏的跳动，引起的身体其它部位的运动；因此人体一直处于运动状态之中，但是相对于肢体运动，呼吸或心跳引起的身体运动非常微弱。

呼吸信号、心跳信号强度幅值随时间变化形成波形。现有技术是通过比较相邻两个点的数值大小来确定波形处于上升阶段或下降阶段，同时对上升波形和下降波形计数。接着统计波峰与波谷出现的次数，找出呼吸波的波峰值与波谷值，当波峰值或波谷值达到预定义好的阈值时，认为这是一个正常的呼吸事件。

当存在噪音或其他因素干扰时，呼吸波形不规整，现有技术容易造成误判，出现计数错误，导致计算结果不准确。另外，预定义好的阈值是固定值，当被测者的状态或所处环境变化时，固定的阈值就不适合了。

发明内容

本发明的目的在于针对现有技术的不足，提供一种基于人体微动信号的呼吸率和心率的提取方法。本发明在安静或睡眠情境下，通过对呼吸和心跳引起的人体微弱体动信号数据的计算处理，获得被测者稳定、准确的呼吸率、心率值。

本发明的目的是通过以下技术方案来实现的：一种基于人体微动信号的呼吸率和心率的提取方法，在安静或睡眠情境下，利用传感器采集人体微动信号，滤除噪声得到混合信号；从所述混合信号中去除心率对应频率的信号，得到呼吸信号，根据呼吸信号持续获得呼吸率；从所述混合信号中去除呼吸率对应频率的信号，得到心跳信号，根据心跳信号持续获得心率值。

进一步地，包括如下步骤：

1.参数初始化

参数初始化包括设置采样频率、采样持续时间、单个样本的采样时长、样本推移的时间步长、初始呼吸率值、初始心率值；

2.采样

2.1接收来自传感器的数据包，每个数据包包含了一个样本推移的时间步长的采样数据。

2.2每个数据包的采样数据按照时序依次拼接，直到达到成单个样本的采样时长，构成一个完整的样本；之后每接收一个新的数据包，删除最先的数据包并将新数据包按时序拼接在上一个样本的末尾，得到新的样本。

3.滤波

滤除高于心跳频率的高频信号；所述高频信号为大于等于预设阈值的信号。经过高频滤波后的样本数据主要是呼吸与心跳信号的混合信号。

4.提取呼吸率值

4.1呼吸波形的过滤计算。利用步骤6.1确定的心跳频率，其中首次执行时利用初始心率值，滤除步骤3得到的混合信号中对应频率的信号获得呼吸信号。

4.2呼吸计次。根据步骤6.1确定的呼吸频率，其中首次执行时根据初始呼吸率值，确定步骤4.1获得的呼吸信号的分段个数，为单个样本的采样时长内呼吸的次数；在步骤4.1得到的呼吸信号中取一个分段，找到分段内所有的极大值，将其中的最大极大值计为一次呼吸，合计单个样本内的呼吸信号出现的最大极大值个数。

5.提取心率值

5.1心跳波形的过滤计算。利用步骤6.1确定的呼吸频率，其中首次执行时利用初始呼吸率值，过滤步骤3得到的混合信号中对应频率的信号获得心跳信号。

5.2心跳计次。利用步骤6确定的心跳频率，其中首次执行时利用初始心率值，确定步骤5.1获得的心跳信号的分段个数，为单个样本的采样时长内心跳的次数；在步骤5.1得到的心跳信号中取一个分段，找到分段内所有的极大值，将其中的最大极大值计为一次心跳，合计单个样本内的心跳信号出现的最大极大值个数。

6.过滤参数优化

6.1根据步骤4、步骤5中呼吸信号、心跳信号的最大极大值个数，确定呼吸频率和心跳频率。

6.2查看呼吸率、心率是否有变化。如果当前样本获得的呼吸率、心率和上一样本不同，则基于当前样本重复执行步骤4～6.1；如果重复执行三次后仍有变化，则直接进入步骤7；如果没变化，则得到当前样本的准确呼吸率、心率值，进入步骤7。

7.按时间步长推移样本

跳到步骤2获得新的样本，循环执行步骤2～6，持续获得呼吸率、心率值，直至达到步骤1设置的采样持续时间，完成计算处理的全部过程。

进一步地，所述传感器为压电传感器。

进一步地，所述采样频率为25Hz；所述采样持续时间为8h；所述单个样本的采样时长为60s；所述样本推移的时间步长为5s。

进一步地，步骤3中的预设阈值取2～3Hz。

进一步地，初始心率值设为75次/分钟。

进一步地，初始呼吸率值设为20次/分钟。

本发明的有益效果是：

(1)本发明能够处理呼吸、心跳引起的人体微弱运动混合在一起的信号数据，消减干扰信号；

(2)本发明利用采样数据连续更新样本，分段计算，获得稳定、客观的结果；

(3)本发明的滤波参数持续优化，计算结果趋向准确。

附图说明

图1是本发明流程图；

图2是睡觉时人体微动的原始采样数据示意图；其中，横坐标为采样点，纵坐标是电压值(单位：毫伏)；

图3是滤除心跳信号后的呼吸波形图；其中，横坐标为采样点，纵坐标是电压值(单位：毫伏)；

图4是呼吸计次图；其中，横坐标为采样点，纵坐标是电压值(单位：毫伏)；

图5是滤除呼吸信号后的心跳波形图；其中，横坐标为采样点，纵坐标是电压值(单位：毫伏)；

图6是心跳计次图；其中，横坐标为采样点，纵坐标是电压值(单位：毫伏)。

具体实施方式

本发明一种基于人体微动信号的呼吸率和心率的提取方法，在安静特别是睡眠状态下，利用灵敏传感器(比如压电传感器)，检测呼吸和心跳引起的人身体的微弱运动，采集人体微动信号数据，通过对采样数据样本放大、滤波、提取等计算处理，获知受测者每分钟的呼吸次数(呼吸率)、每分钟心跳次数(心率)。对人体微动的检测往往是长时间进行的，比如检测一个晚上的睡眠，持续数小时。检测结果呼吸率、心率需要在开始检测后的较短时间内显示，比如60秒。所以检测分段进行，处理计算传感器采集的一段数据，获得结果。检测过程持续数小时地长时间进行，不断地采集数据，在开始检测后的较短时间内不断地获得计算结果。检测过程不对被测者的生理活动过程产生干扰，容易实施检测，获得稳定准确的结果，有利于给出真实客观评估，有效帮助提高健康质量。

下面以某一种实施形式为例，结合图1至图6进一步说明本发明。

如图1所示，本实施例包括以下步骤：

1.参数初始化。

参数初始化包括设置采样频率(比如25Hz)，设置采样持续时间(比如8h)，设置单个样本的采样时长(比如60s)，设置样本推移的时间步长(比如5s)，同时，呼吸率、心率值归零。

2.采样。

2.1接收来自传感器的数据包，每个数据包包含了一个样本推移的时间步长的采样数据。

2.2每个数据包的采样数据按照时序依次拼接，直到达到成单个样本的采样时长，构成一个完整的样本，如图2所示；之后每接收一个新的数据包，删除最先的数据包并将新数据包按时序拼接在上一个样本的末尾，得到新的样本，实现样本推移。

3.滤波。

样本数据中包含有呼吸信号、心跳信号和噪声。正常成人呼吸频率在12～20次/分钟的范围内，也就是频率为0.2～0.33Hz。正常成人心跳频率在50～100次/分钟的范围内，也就是频率为0.83～1.67Hz。

滤除噪声，主要是滤除明显高于心跳频率的高频信号；所述高频信号为大于等于预设阈值的信号，预设阈值一般取2～3Hz，比如电源线路的50Hz信号。经过高频滤波后的样本数据主要包含的是呼吸与心跳混合在一起的混合信号。

为了提取呼吸率与心率，把滤波后的样本数据分别处理，分别执行步骤4、步骤5两个过程。

4.提取呼吸率值。

4.1呼吸波形的过滤计算。利用步骤6.1确定的心跳频率，滤除步骤3得到的混合信号中对应频率的信号获得呼吸信号，如图3所示。初始时，心率值预估为75次/分钟。

4.2呼吸计次。根据步骤6.1确定的呼吸频率，确定步骤4.1获得的呼吸信号的分段个数。初始时，呼吸率值预估为20次/分钟。分段个数为单个样本的采样时长内呼吸的次数；将步骤4.1得到的呼吸信号分段，在呼吸信号中取一个数据分段，找到分段内所有的极大值，将其中的最大极大值计为一次呼吸，并记录其出现的时间。在整个样本中，逐个数据地往前推移数据分段，得到全部的数据分段最大极大值及其出现的时间。合计单个样本内的呼吸信号出现的最大极大值个数，即得到呼吸率值，如图4所示，虚线连接的各圆点即是该样本数据内出现的最大极大值。

5.提取心率值

5.1心跳波形的过滤计算。利用步骤6.1确定的呼吸频率，过滤步骤3得到的混合信号中对应频率的信号获得心跳信号，如图5所示。初始时，呼吸率值预估为20次/分钟。

5.2心跳计次。利用步骤6确定的心跳频率，确定步骤5.1获得的心跳信号的分段个数。初始时，心率值预估为75次/分钟。分段个数为单个样本的采样时长内心跳的次数；将步骤5.1得到的心跳信号分段，在心跳信号中取一个数据分段，找到分段内所有的极大值，将其中的最大极大值计为一次心跳，并记录其出现的时间。在整个样本中，逐个数据地往前推移数据分段，得到全部的数据分段最大极大值及其出现的时间。合计单个样本内的心跳信号出现的最大极大值个数，即得到心率值，如图6所示，虚线连接的各圆点即是该样本数据内出现的最大极大值。

6.过滤参数优化

6.1根据步骤4、步骤5确定呼吸频率和心跳频率。

6.2查看呼吸率、心率是否有变化。如果当前单个样本获得的呼吸率、心率和上一个样本不同，则不采用当前呼吸率、心率的监测值，基于当前样本重复执行步骤4～6.1；如果重复执行三次后仍与上一个样本不同，则直接进入步骤7；如果没变化，则已稳定得到当前样本的准确呼吸率、心率值，进入步骤7。

7.按时间步长推移样本

跳到步骤2获得新的样本，循环执行步骤2～6，持续获得呼吸率、心率值；直至达到步骤1设置的采样持续时间，完成计算处理的全部过程。

本说明书中所描述的具体实施例仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代，但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。