一种四电极人体阻抗获取心率信号电路及方法与流程

本发明属于阻抗测量技术领域，特别涉及一种通过人体阻抗测量获取心率信号的电路及方法。

背景技术：

近年随着人们生活水平的提高，对自身的健康状况日益关注。例如，人们除了关注自己的体重之外，还进一步开始关注人体脂肪的含量，从而可以更加准确的把握自身的胖瘦程度。技术上测量脂肪含量有多种方法，例如经典的水下称重法，ct扫描法，以及人体阻抗测量法。

人体阻抗测量的方法的基本原理是人体脂肪含水极小，因此相对于人体其他组织其阻抗很高；通过测量人体阻抗可以大致判断人体脂肪和非脂肪组织的比例，从而推算脂肪含量。人体阻抗测量大致分为直流阻抗测量和交流阻抗测量。直流阻抗测量通过在人体两个不同点输入直流电流，测量两点之间的电压来获得该两点之间的阻抗，优点是测量电路简单，缺点是测量结果包含了人体皮肤的阻抗；而在直流情况下，人体皮肤阻抗的值很大，且容易受皮肤自身状况和测量电极和皮肤接触状况影响，因此总的人体阻抗测量结果的准确度较差。

随着bia(生物电阻抗)技术的普及与应用，通过bia获取人体成分参数来检测人体健康状况的产品也越来越多和越来越完善，采集阻抗信息的方式也由2电级进化到4电极和8电极，激励电流也从单一频率增加到多频段分析，检测的参数也越来越丰富，不仅可以检测出阻抗值，进而计算获取各种人体成分给健康生活提供更准确更有用的参考，而且还可以通过上肢阻抗信息获取心跳信号，给人体健康监护提供更多一层的支持和保障。但是，原有通过人体阻抗测量心跳的地采集方式需要通过贴电极片到胸腔或腹部等身体部位，采集方式繁琐、不方便。而直接通过手握式的bia获取心跳信号仍然没有广泛普及开来，主要原因是，通过阻抗获取的心跳信号相对干扰较多，无法准确和有效的获取心率。

如专利申请201610168316.0公开了一种基于生物电阻抗技术的心率测量方法及装置，所述方法利用生物电阻抗技术，通过采集受试者手臂桡动脉的电阻抗信号，并进行模数转换后，得到电阻抗数字信号。从所述电阻抗数字信号中提取心率和心率变异性特征。本发明具有简单易行、廉价和操作简便等特点，能够快速准确地采集并分析心率数据，实现对受试者的测量。该方法虽然能通过手臂桡动脉采集电阻抗信号，然而，这种方式获取的心跳信号相对干扰较多，无法准确和有效的获取心率。

技术实现要素：

基于此，因此本发明的首要目地是提供一种四电极人体阻抗获取心率信号电路及方法，该电路及方法能够有效除去基线漂移强烈影响，通过幅度监测(前后峰谷幅度的跳变关系限制)运动信息以及标准心跳模版的建立与比对，能较为准确的排除基线漂移和运动干扰，从而得到较为准确的心率，不仅能准确测量出阻抗及人体成分，还能方便的测量出心跳信号。

本发明的另一个目地在于提供一种四电极人体阻抗获取心率信号电路及方法，该电路及方法通过改进采集阻抗信息的方式改善采集阻抗信息的便捷性和准确性；通过设计阻抗心跳信号提取算法来移植基线漂移和排除运动干扰，从而准确测量出心跳信号。

本发明的再一个目的是提供一种四电极人体阻抗获取心率信号电路及方法，该电路及方法能够通过多种方式(包括但不局限手握式或脚踩式这两种测量方式)获取阻抗的心跳信号，简便实用；另一方面不仅可以获取阻抗值，进而分析人体成份，还可以有效通过心跳算法获取心率，使得健康监测得更为全面和准确；通过设计针对性的算法来抑制极限漂移和排除运动干扰，以准确的检测心率。

为实现上述目的，本发明的技术方案为：

一种四电极人体阻抗获取心率信号电路，其包括有至少四电极，四电极中，n+和n-为激励电流的产生电极和接收电极，i+和i-为测量电极，即从n+和n-间输出不同频率的激励电流，通过测量i+和i-间电压信号，采样后通过采样信号放大电路对采样信号进行放大，其特征在于所述i+和i-两个测量电极与采样信号放大电路间分别串入电容，可以有效滤除信号中无效的直流成分；在n+和i+间、n-和i-间接入电阻，以去除极化电压的影响，解决基线漂移的现象。

进一步，所述2个激励电极间加入释放电阻，以便在没有人体接触时为2个激励电极做一个释放回路。

更进一步，由于测量电极处于悬空时状态不稳定，可能会对人体接触判断造成误判，所述测量电极接入1/2vs，可以在悬空时给一个基准电压，就可以分辨出有无人体接触。

一种四电极人体阻抗获取心率信号方法，其特征在于该方法通过四电极采集人体的阻抗信号，对采集到的信号进行整系数低通滤波处理，然后获取基线值(根据窗口大小)，再将原始的采集到的信号减去基线值，然后进行整系数高通滤波，去除基线漂移；再进行整系数低通滤波除去高频干扰，然后进行滤波处理，滤波后进行凸凹算法获取峰谷信号，根据当前峰谷信息提取心跳模版，然后将后续心跳信号与心跳模版进行分析，提取峰谷信息，再计算心率。

进一步，所述方法中，具体步骤如下：

101、对ad采样信号进行整系数的低通滤波处理，去除高频干扰；

102、获取基线值；

通常情况下，取窗口大小为w0的平滑滤波，保存最新w0个lpdata，除去其中的最大值lmax和最小值lmin，然后剩下的做平均值作为当前的基线值gbasevalue；

103、原始信号去除基线值，将ad采样信号减去基线值，获得缩小信号；

104、对缩小信号进行整系数的高通滤波处理，去除基线漂移；

105、进一步做整系数低通滤波处理，滤除高于心跳频率范围的干扰信号，得到滤波结果信号；

进一步，对滤波后的结果信号进行移动窗口wb的信号幅度监测，并判断是否超过运动信号跳变的阈值，超过阀值则认为是监测到运动，返回步骤101，否则继续106。

106、对滤波后的结果信号进行滤波处理，

进一步，所述滤波处理包括平滑滤波处理和积分滤波处理，其中平滑滤波处理是将w1个数的平滑滤波处理，得平滑滤波处理信号；

然后再对平滑滤波处理信号进行积分处理，积分的作用一方面凸显了心跳信号，减小了小幅度信号的影响，另一方面也起到了平滑的效果，得到心跳波形。

107、用凹凸算法查找出所有的峰和谷的信息，并保存起来；

具体地说，凹凸判断规则如下：

先保存最新ns点gprodata到gfdata[ns]中然后：

峰判断：gfdata[i]<gfdata[i-1]且gfdata[i-1]<＝gfdata[i-2]且gfdata[i-2]>gfdata[i-3]且gfdata[i-3]>gfdata[i-4]；

谷判断：gfdata[i]>gfdata[i-1]且gfdata[i-1]>＝gfdata[i-2]且gprodata[i-2]<gfdata[i-3]且gfdata[i-3]<gfdata[i-4]；

找到波峰或者波谷时，保存到波谷数组gpeakpoint[n]、gpeakvalue[n]或gtroughpoint[n]、gtroughvalue[n]中，并均保准在最后一个位置中，保存前将第一个位置的去掉，后面依次前移。

根据保存的峰谷信息(gpeakpoint[n]、gpeakvalue[n]和gtroughpoint[n]、gtroughvalue[n]，n＝5)，综合判断最新的峰或谷是否是标准的心跳峰谷波形，并根据找到的标准峰或谷模版更新峰谷模版。

108、将最新找到峰或谷信息y与标准峰或谷模版进行匹配分析，即做两者的相关性分析(ρxy)，然后提取峰谷信息，再计算心率。

具体地说，相关性分析公式为：其中n为数据个数，为x的均值，为y的均值；

根据求得的ρxy判断当前峰谷与标准峰谷模版的相关性，在范围[-1，+1]，值越大，相关性越大，设置ρxy>＝gthxy时，为相关，即找到正确峰或谷信息，否则未找到，进入步骤s1。

进一步，根据ρxy结果决定是否保存当前峰谷信息，具有相关性则进行保存；

更进一步，对保存的峰谷信息进行判断，判断方式分幅度及间距的范围限制和前后关系：

幅度：峰谷的幅度在允许的范围内；

间距：在合理范围内(在标准峰谷的最小间距和最大间距之间)；

前后关系：幅度跳变不超过前一幅度的a1倍，间距不超过前已间距的a2倍。

更进一步，根据获取的有效峰谷间距值gpeaksample[n]、gtroughsample[n]来计算心率gbreathvalue，并做结果的四舍五入和平滑处理。

更进一步，所述平滑处理为：保存最新rn个gbreathvalue到gresultdata[rn]中，然后去除其中最大max和最小min值，余下的做平均作为最终的gbreathvalue。

所述四电极采用手握接触式或脚踩接触式测量阻抗信息和心跳信息。

本发明所实现的四电极人体阻抗获取心率信号电路及方法，能够有效除去基线漂移强烈影响，通过幅度监测(前后峰谷幅度的跳变关系限制)运动信息以及标准心跳模版的建立与比对，能较为准确的排除基线漂移和运动干扰，从而得到较为准确的心率，不仅能准确测量出阻抗及人体成分，还能方便的测量出心跳信号。

且本发明实现方便快捷，通过电路设计和优化心跳提取算法能很好的抑制强烈的基线漂移和排除运动干扰，从而准确和有效地计算出心率。

该电路及方法能够通过多种方式(包括但不局限手握式或脚踩式这两种测量方式)获取阻抗的心跳信号，简便实用；另一方面不仅可以获取阻抗值，进而分析人体成份，还可以有效通过心跳算法获取心率，使得健康监测得更为全面和准确；通过设计针对性的算法来抑制极限漂移和排除运动干扰，以准确的检测心率。

附图说明

图1是本发明所实施四电级测人体阻抗心跳的电路原理图。

图2是本发明所实施的心跳提取流程图。

图3是本发明所实施测量的阻抗心跳原始波形图。

图4是本发明所实施的阻抗心跳去基线漂移波形图。

图5是本发明所实施测量的心跳运动监测示意图。

图6是图5的去基线及处理后心跳波形图。

图7是本发明所实施的查找波谷及分析筛选波形图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

参见附图1所示，为本发明所实现的四电极采集阻抗信息的电路原理图，图中所示，其中，采用四电极手握接触式或脚踩接触式测量阻抗信息和心跳信息。n+和n-为激励电流的产生电极和接收电极，i+和i-为测量电极，即从n+和n-间输出不同频率的激励电流，通过测量i+和i-间电压信号，然后转换成阻抗信号，进而提取阻抗值和心跳信号。

为了解决测量中的基线漂移，电极极化干扰，接触判断等难点。电路做了如下处理。无人体接触时，两个激励电极和测量电极处于悬空状态，在2个激励电极间加入释放电阻(r1)，以便在没有人体接触时为2个激励电极做一个释放回路。由于测量电极处于悬空时状态不稳定，可能会对人体接触判断造成误判，为测量电极接入1/2vs可以在悬空时给一个基准电压，就可以分辨出有无人体接触，方便后续处理。在测量电极中串入电容(c1、c2)可以有效滤除信号中无效的直流成分。在n+和i+间、n-和i-间接入电阻(r2、r3)以去除极化电压的影响，解决基线漂移的现象。

本发明所实现的四电极人体阻抗获取心率信号方法，通过四电极采集人体的阻抗信号，对采集到的信号进行整系数低通滤波处理，然后获取基线值(根据窗口大小)，再将原始的采集到的信号减去基线值，然后进行整系数高通滤波，去除基线漂移；再进行整系数低通滤波除去高频干扰，然后进行滤波处理，滤波后进行凸凹算法获取峰谷信号，根据当前峰谷信息提取心跳模版，然后将后续心跳信号与心跳模版进行分析，提取峰谷信息，再计算心率。

心跳提取算法具体步骤如图2所示：

步骤s1：对ad采样信号sampledatas(如图3所示，原始获得的采样信号有较强的基线漂移干扰)进行整系数的低通滤波处理，去除高频干扰lpdata。

步骤s2：窗口大小为w0的平滑滤波，保存最新w0个lpdata，除去其中的最大值lmax和最小值lmin，然后剩下的做平均值作为当前的基线值gbasevalue。

步骤s3：原始信号去除基线值，即：谷gbasedata＝sampledatas-gbasevalue，这样将sampledatas从较大的整数统一缩小在0附近的值gbasedata，而不影响原来信号的波动情况。

步骤s4：对gbasedata进行整系数的高通滤波处理，目的是去除基线漂移，如图4所示为去除基线后的波形图。

步骤s5：对步骤s4处理结果进一步做整系数低通滤波处理，目的是滤除高于心跳频率范围的干扰信号，得到滤波结果gfiltedata。

步骤s6：对滤波后的结果gfiltedata进行移动窗口wb的信号幅度监测，并判断是否超过运动信号跳变的阈值gsportth，当监测到运动(如图5所示，方形为监测到了运动)返回步骤s1，否则继续s7。

步骤s7：对gfiltedata进行w1个数的平滑滤波处理处理处理，得：gsmoothdata。

步骤s8：在对gsmoothdata进行积分处理，积分窗口大小w2，积分的作用一方面凸显了心跳信号，减小了小幅度信号的影响，另一方面也起到了平滑的效果，得到心跳波形gprodata，如图6所示，下方为去基线后波形，上方为预处理后的心跳波形。

步骤s9：用凹凸算法查找出所有的峰和谷的信息，并保存起来。

凹凸判断规则如下：

先保存最新ns点gprodata到gfdata[ns]中然后：

峰判断：gfdata[i]<gfdata[i-1]且gfdata[i-1]<＝gfdata[i-2]且gfdata[i-2]>gfdata[i-3]且gfdata[i-3]>gfdata[i-4]；

谷判断：gfdata[i]>gfdata[i-1]且gfdata[i-1]>＝gfdata[i-2]且gprodata[i-2]<gfdata[i-3]且gfdata[i-3]<gfdata[i-4]；

找到波峰或者波谷时，保存到波谷数组gpeakpoint[n]、gpeakvalue[n]或gtroughpoint[n]、gtroughvalue[n]中，并均保准在最后一个位置中，保存前将第一个位置的去掉，后面依次前移。

步骤s10：根据保存的峰谷信息(gpeakpoint[n]、gpeakvalue[n]和gtroughpoint[n]、gtroughvalue[n]，n＝5)，综合判断最新的峰或谷是否是标准的心跳峰谷波形，找到标准峰或谷模版后，更新峰谷模版。

步骤s11：最新找到峰或谷信息y与标准峰或谷模版x进行匹配分析，即做两者的相关性分析(ρxy)，公式：其中n为数据个数，为x的均值，为y的均值；

根据求得的ρxy判断当前峰谷与标准峰谷模版的相关性，在范围[-1，+1]，值越大，相关性越大，设置ρxy>＝gthxy时，为相关，即找到正确峰或谷信息，否则未找到，进入步骤s1。

步骤s12：根据ρxy结果决定是否保存当前峰谷信息，如图7所示最上方为积分处理后波形，中间为找到的波谷示意图，最下面为筛选后的波谷点。

步骤s13：对保存的峰谷信息进行判断，判断方式分幅度及间距的范围限制和前后关系：

幅度：峰谷的幅度在合理范围内(由硬件采集决定)；

间距(gpeaksample[n]/gtroughsample[n])：在合理范围内(在标准峰谷的最小间距和最大间距之间)；

前后关系：幅度跳变不超过前一幅度的a1倍，间距不超过前已间距的a2倍。

步骤s14：根据获取的有效峰谷间距值gpeaksample[n]、gtroughsample[n]来计算心率gbreathvalue，并做结果的四舍五入和平滑处理。

平滑处理：保存最新rn个gbreathvalue到gresultdata[rn]中，然后去除其中最大max和最小min值，余下的做平均作为最终的gbreathvalue。

总之，本发明所实现的四电极人体阻抗获取心率信号电路及方法，能够有效除去基线漂移强烈影响，通过幅度监测(前后峰谷幅度的跳变关系限制)运动信息以及标准心跳模版的建立与比对，能较为准确的排除基线漂移和运动干扰，从而得到较为准确的心率，不仅能准确测量出阻抗及人体成分，还能方便的测量出心跳信号。

且本发明实现方便快捷，通过电路设计和优化心跳提取算法能很好的抑制强烈的基线漂移和排除运动干扰，从而准确和有效地计算出心率。

该电路及方法能够通过多种方式(包括但不局限手握式或脚踩式这两种测量方式)获取阻抗的心跳信号，简便实用；另一方面不仅可以获取阻抗值，进而分析人体成份，还可以有效通过心跳算法获取心率，使得健康监测得更为全面和准确；通过设计针对性的算法来抑制极限漂移和排除运动干扰，以准确的检测心率。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。