心率提取方法及装置、电子设备、压电式传感器及介质与流程

医药医疗技术的改进;医疗器械制造及应用技术1.本技术涉及心率测量技术领域，尤其是涉及一种心率提取方法及装置、电子设备、压电式传感器及计算机可读存储介质。背景技术：2.当前，人们对于自身健康状况越来越关注，而心率作为一项重要的生命体征指标相应受到重视。通过压电式传感器可以检测到人体的生命体征信号，进而从该生命体征信号中可以分离得到人体的呼吸率和心率。在实际应用中，由于心跳信号比较微弱，且容易受到呼吸谐波影响，从而导致对心率提取的准确性低，而这种情况常见于临床心血患者或呼吸不均匀人群中。3.然而，已知的心率提取方法主要存在以下问题，即要么是没有考虑呼吸谐波对心跳信号的影响，导致心率提取的准确性很低；要么是对任何信号都进行呼吸谐波抵消处理，占用较多的运算资源。技术实现要素：4.为解决现有存在的技术问题，本技术提供一种能够灵活且准确地提取心率，运算资源利用率高的心率提取方法及装置、电子设备、压电式传感器及计算机可读存储介质。5.为达到上述目的，本技术实施例的技术方案是这样实现的：6.第一方面，本技术实施例提供一种心率提取方法，所述方法包括：7.获取目标对象的生命体征信号；8.检测所述生命体征信号中是否包含呼吸谐波，得到检测结果；9.根据所述检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率。10.第二方面，本技术实施例提供一种心率提取装置，所述装置包括：11.获取模块，用于获取目标对象的生命体征信号；12.呼吸谐波检测模块，用于检测所述生命体征信号中是否包含呼吸谐波，得到检测结果；13.心率提取模块，用于根据所述检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率。14.第三方面，本技术实施例提供一种电子设备，包括处理器、存储器及存储在所述存储器上并可被所述处理器执行的计算机程序，所述计算机程序被所述处理器执行时实现本技术任一实施例所述的心率提取方法。15.第四方面，本技术实施例提供一种压电式传感器，包括处理器、存储器及存储在所述存储器上并可被所述处理器执行的计算机程序，所述计算机程序被所述处理器执行时实现本技术任一实施例所述的心率提取方法。16.第五方面，本技术实施例提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，所述计算机程序被所述处理器执行时实现本技术任一实施例所述的心率提取方法。17.本技术上述实施例提供的心率提取方法及装置、电子设备及计算机可读存储介质，获取目标对象的生命体征信号，检测所述生命体征信号中是否包含呼吸谐波，得到检测结果，根据所述检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率，如此，通过检测目标对象的生命体征信号中是否包含呼吸谐波，进而根据检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率，从而实现灵活且准确地提取心率，并提高了运算资源利用率。此外，可为生命健康检测提供更准确的参考。附图说明18.图1为适用本技术实施例的应用环境示意图；19.图2为本技术实施例中心率提取方法的流程图；20.图3为本技术一实施例中心率提取装置的示意图；21.图4为本技术另一实施例中心率提取装置的示意图；22.图5为本技术一实施例中自适应滤波示意图；23.图6为本技术一实施例中电子设备的结构示意图；24.图7为本技术一实施例中压电式传感器的结构示意图；25.图8为本技术另一实施例中电子设备的结构示意图；26.图9为本技术一实施例中计算机可读存储介质的示意图。具体实施方式27.以下结合说明书附图及具体实施例对本技术技术方案做进一步的详细阐述。28.除非另有定义，本文所使用的所有的技术和科学术语与属于本技术的技术领域的技术人员通常理解的含义相同。本文中在本技术的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本技术的实现方式。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。29.在本技术的描述中，需要理解的是，术语“中心”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本技术和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本技术的限制。在本技术的描述中，除非另有说明，“多个”的含义是两个或两个以上。30.在本技术的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本技术中的具体含义。31.请参阅图1，图1为适用于本技术实施例的一种应用环境示意图。本技术实施例提供的心率提取方法可以应用于如图1所示的智能家居系统中，该智能家居系统包括压电式传感器101和计算机设备102。32.其中，压电式传感器101可用于在贴附在目标对象的皮肤表面时，采集目标对象的生命体征信号。计算机设备102是指具备对数据进行处理和存储功能的设备，具体可以为个人计算机(personal computer，pc)、手机、平板电脑、穿戴式电子设备、服务器或者其他具有计算能力的电子设备，本技术实施例对此不作限定。计算机设备102与压电式传感器101可以为分离的硬件设备，计算机设备102可以接收压电式传感器101采集的目标对象的生命体征信号，然后作后续的处理，如对生命体征信号进行心率提取。在本技术实施例中，目标对象是指具有心率的人或动物等。33.压电式传感器101与计算机设备102之间可以通过有线或者无线的方式进行通信。例如，压电式传感器101与计算机设备102之间的数据传送可以采用设备到设备(ad-hoc)的方式，也可以在基站或无线访问点(access point，ap)的协调下进行，本技术实施例对此不作限定。34.在一些实施例中，压电式传感器101是独立的，压电式传感器101可以设置于其他电子设备上，例如，压电式传感器101可设置于计算机设备102，也可设置于除计算机设备102外的如呼吸机等其他设备，本技术实施例对此不作限定。在另一些实施例中，压电式传感器101与计算机设备102为集成为一体硬件设备，通过一体的压电式传感器101实施本技术实施例提供的心率提取方法。可以理解的是，以上应用环境仅为方便理解下述实施例作出的示例性描述，并不对本技术实施例实际可实施的应用环境构成任何限定。35.请参阅图2，为本技术一实施例提供的一种心率提取方法，可应用于上述计算机设备，所述心率提取方法包括但不限于s101、s103和s105，具体介绍如下：36.s101，获取目标对象的生命体征信号。37.其中，所述目标对象可以为人体等具有心率的对象。计算机设备获取目标对象的生命体征信号可以是接收压电式传感器发送的其自身所采集的目标对象的生命体征信号，也可以是接收其它终端设备如智能手机等发送的由压电式传感器所采集的目标对象的生命体征信号。需要说明的是，本实施例中以所述生命体征信号为离散时间信号为例，即为通过对原始的生命体征信号进行间隔采样获得的。此外，所述生命体征信号可以为一帧或多帧，本实施例中以获取一帧生命体征信号为例。38.s103，检测所述生命体征信号中是否包含呼吸谐波，得到检测结果。39.可以理解地，因呼吸谐波的存在，包含呼吸谐波的生命体征信号与未包含呼吸谐波的生命体征信号之间将存在一些差异特征，而基于这些差异特征可检测获取的生命体征信号中是否包含呼吸谐波。本实施例中，计算机设备可基于生命体征信号中波峰与波谷的特征检测所述生命体征信号中是否包含呼吸谐波，得到检测结果。40.s105，根据所述检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率。41.其中，根据检测结果的不同，对所述生命体征信号执行心率提取操作也相应不同。具体地，当根据所述检测结果确定所述生命体征信号中不包含呼吸谐波时，可直接对所述生命体征信号进行心率提取，以获得所述目标对象的心率；当根据所述检测结果确定所述生命体征信号中包含呼吸谐波时，需要先消除所述生命体征信号中包含的呼吸谐波，然后再对消除呼吸谐波后的所述生命体征信号进行心率提取，以获得所述目标对象的心率。42.本技术上述实施例中，获取目标对象的生命体征信号，检测所述生命体征信号中是否包含呼吸谐波，得到检测结果，根据所述检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率，如此，通过检测目标对象的生命体征信号中是否包含呼吸谐波，进而根据检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率，从而实现灵活且准确地提取心率，并提高了运算资源利用率。43.在一些实施例中，所述检测所述生命体征信号中是否包含呼吸谐波，得到检测结果，包括：44.对所述生命体征信号进行第一滤波处理，得到呼吸信号；45.检测所述呼吸信号中是否包含呼吸谐波，得到检测结果。46.计算机设备先对所述生命体征信号进行第一滤波处理，从所述生命体征信号中得到呼吸信号，进而检测所述呼吸信号中是否包含呼吸谐波，得到检测结果。其中，所述第一滤波处理可采用带通滤波处理等方式，且当所述第一滤波处理为带通滤波处理时，对应的截止频率可根据呼吸信号的频率大小进行确定。可以理解地，由于呼吸谐波通常包含于呼吸信号中，可通过先对所述生命体征信号进行滤波以得到呼吸信号，进而检测所述呼吸信号中是否包含呼吸谐波，从而实现检测所述生命体征信号中是否包含呼吸谐波，与直接检测所述生命体征信号中是否包含呼吸谐波的方式相比，不仅减少了运算量，提高了处理速度，而且进一步提高了心率提取的准确性。47.在一些实施例中，所述对所述生命体征信号进行第一滤波处理，得到呼吸信号，包括：48.基于预设第一频率区间对获取的目标对象的生命体征信号进行第一带通滤波处理，得到去除基线与高频干扰的生命体征信号；49.基于预设第二频率区间对所述去除基线与高频干扰的生命体征信号进行第二带通滤波处理，得到呼吸信号。50.可以理解地，受到测量环境等因素的干扰，通过压电式传感器采集的目标对象的生命体征信号中除了包括呼吸信号之外，可能还包括其它非必要信号如高频干扰等，而高频干扰等其它非必要信号的存在将影响对呼吸信号的处理，因此，计算机设备可先基于预设第一频率区间对获取的目标对象的生命体征信号进行第一带通滤波处理，以得到去除基线与高频干扰的生命体征信号，然后基于预设第二频率区间对所述去除基线与高频干扰的生命体征信号进行第二带通滤波处理，以得到呼吸信号。其中，所述第一频率区间可根据基线与高频干扰等非必要信号所处的频率范围进行设置，比如可设置第一频率区间为0.1～4hz 等。所述第二频率区间可根据呼吸信号所处的频率范围进行设置，比如可设置第二频率区间为0.1～0.5hz等。在实际应用中，也可能通过对生命体征信号进行一次带通滤波处理操作而得到呼吸信号。如此，通过对目标对象的生命体征信号依次执行两次不同的带通滤波处理，进而获得目标对象的呼吸信号，进一步提高了心率提取的准确性。51.在一些实施例中，所述检测所述呼吸信号中是否包含呼吸谐波，得到检测结果，包括：52.对所述呼吸信号进行峰值检测，得到第一波峰序列和第一波谷序列；53.按照预设规则分别对所述第一波峰序列和所述第一波谷序列进行处理，得到第二波峰序列和第二波谷序列；54.判断基于所述第二波峰序列和所述第二波谷序列所确定的波形是否满足预设对称条件，得到所述呼吸信号中是否包含呼吸谐波的检测结果。55.其中，所述对所述呼吸信号进行峰值检测的目的是查找出所述呼吸信号中所有的波峰和波谷，以得到第一波峰序列和第一波谷序列。具体地，在所述呼吸信号中顺序查找幅值既大于左侧，又大于右侧的点为波峰，而幅值既小于左侧，又小于右侧的点为波谷。需要说明的是，所述第一波峰序列既包括各波峰的幅值，还包括对应位置即时间序号，所述第一波谷序列既包括各波谷的幅值，还包括对应位置即时间序号。这里，基于包含呼吸谐波的呼吸信号与未包含呼吸谐波的呼吸信号之间存在的差异特征，可按照预设规则分别对所述第一波峰序列和所述第一波谷序列进行处理，比如删除明显非呼吸谐波的波峰和波谷等，以得到第二波峰序列和第二波谷序列。优选地，所述按照预设规则分别对所述第一波峰序列和所述第一波谷序列进行处理，得到第二波峰序列和第二波谷序列，包括：删除所述第一波峰序列中幅值与最大波峰幅值的比值小于预设第一比值阈值的波峰、以及位置间隔小于预设第一距离阈值的相邻波峰中幅值小的波峰，并在位置间隔大于预设第二距离阈值的相邻波峰中插入一目标波峰，得到第二波峰序列；其中，所述第二距离阈值大于所述第一距离阈值，所述目标波峰的幅值为相邻波峰幅值的平均值，且所述目标波峰的位置为相邻波峰位置的平均值；删除所述第一波谷序列中幅值与最大波谷幅值的比值小于预设第二比值阈值的波谷、以及位置间隔小于预设第三距离阈值的相邻波谷中幅值小的波谷，并在位置间隔大于预设第四距离阈值的相邻波谷中插入一目标波谷，得到第二波谷序列；其中，所述第四距离阈值大于所述第三距离阈值，所述目标波谷的幅值为相邻波谷幅值的平均值，且所述目标波谷的位置为相邻波谷位置的平均值。这里，所述第一比值阈值和所述第二比值阈值可以根据实际情况需要进行设置，比如可将所述第一比值阈值和所述第二比值阈值都设置为0.5，以将所述第一波峰序列中幅值与最大波峰幅值的比值小于0.5的波峰删除，以及将所述第一波谷序列中幅值与最大波谷幅值的比值小于0.5的波谷删除。所述第一距离阈值、所述第二距离阈值以及所述第三距离阈值和所述第四距离阈值可以根据实际情况需要进行设置，比如可根据所述第一波峰序列中相邻波峰之间的位置间隔平均值设置所述第一距离阈值和所述第二距离阈值，以及根据所述第一波谷序列中相邻波谷之间的位置间隔平均值设置所述第三距离阈值和所述第四距离阈值，当然，也可根据正常生命体征信号中相邻波峰之间的位置间隔设置所述第一距离阈值和所述第二距离阈值，以及根据正常生命体征信号中相邻波谷之间的位置间隔设置所述第三距离阈值和所述第四距离阈值。需要说明的是，所述相邻波峰是指位置间隔最小的两个波峰，所述相邻波谷是指位置间隔最小的两个波谷，其中，位置间隔是指波峰点或波谷点对应的时间序号之差。如此，通过对目标对象的呼吸信号中波峰序列和波谷序列进行删峰和补峰处理，以获得可用于判断波形是否对称的新波峰序列和新波谷序列，操作简便，且提高了心率提取的速度。56.在一些实施例中，所述判断基于所述第二波峰序列和所述第二波谷序列所确定的波形是否满足预设对称条件，得到所述呼吸信号中是否包含呼吸谐波的检测结果，包括：57.获取基于所述第二波峰序列和所述第二波谷序列所确定的波形中每个周期内波峰幅值与波谷幅值的第一比值；其中，相邻的波峰和波谷组成一个周期；58.根据所述每个周期内波峰幅值与波谷幅值的第一比值，获取所述第一比值的平均值；59.在所述第一比值的平均值不满足预设第一对称条件时，确定所述呼吸信号中包含呼吸谐波，否则确定所述呼吸信号中不包含呼吸谐波。60.其中，由于所述第二波峰序列不仅包括各波峰的幅值，还包括对应位置即时间序号，所述第二波谷序列不仅包括各波谷的幅值，还包括对应位置即时间序号，因此，根据各波峰的对应位置和各波谷的对应位置可获知基于所述第二波峰序列和所述第二波谷序列所确定的波形中相邻的波峰和波谷，即可获知每个周期内波峰幅值与波谷幅值。在获取所述每个周期内波峰幅值与波谷幅值的第一比值后，可结合所述第二波峰序列和所述第二波谷序列所确定的波形包含的周期数量，获取所述第一比值的平均值。其中，所述第一对称条件可根据实际情况需要进行设置，比如，所述第一对称条件可设置为所述第一比值的平均值大于1-a1且小于a1，其中，a1可取值为1.3或1.35等。而在所述第一比值的平均值满足预设第一对称条件时，说明基于所述第二波峰序列和所述第二波谷序列所确定的波形上下对称，否则说明基于所述第二波峰序列和所述第二波谷序列所确定的波形上下不对称。在所述第一比值的平均值不满足预设第一对称条件，即基于所述第二波峰序列和所述第二波谷序列所确定的波形上下不对称时，表示检测到呼吸谐波，即所述呼吸信号中包含呼吸谐波，否则确定所述呼吸信号中不包含呼吸谐波。如此，通过检测基于所述第二波峰序列和所述第二波谷序列所确定的波形是否对称而判断所述呼吸信号中是否包含呼吸谐波的检测结果，操作简单且快速，进一步提高了心率提取的准确性和速度。61.在一些实施例中，所述判断基于所述第二波峰序列和所述第二波谷序列所确定的波形是否满足预设对称条件，得到所述呼吸信号中是否包含呼吸谐波的检测结果，包括：62.获取基于所述第二波峰序列和所述第二波谷序列所确定的波形中每个周期内第一幅值与第二幅值的第二比值；其中，相邻的波峰和波谷组成一个周期，所述第一幅值为波峰位置与波谷位置的位置差相对应的幅值，所述第二幅值为波谷位置与波峰位置的位置差相对应的幅值；63.根据所述每个周期内第一幅值与第二幅值的第二比值，获取所述第二比值的平均值；64.在所述第二比值的平均值不满足预设第二对称条件时，确定所述呼吸信号中包含呼吸谐波，否则确定所述呼吸信号中不包含呼吸谐波。65.其中，由于所述第二波峰序列不仅包括各波峰的幅值，还包括对应位置即时间序号，所述第二波谷序列不仅包括各波谷的幅值，还包括对应位置即时间序号，因此，根据各波峰的对应位置和各波谷的对应位置可获知基于所述第二波峰序列和所述第二波谷序列所确定的波形中相邻的波峰和波谷，即可获知每个周期内波峰幅值与波谷幅值以及波峰位置与波谷位置。在获取所述每个周期内第一幅值与第二幅值的第二比值后，可结合所述第二波峰序列和所述第二波谷序列所确定的波形包含的周期数量，获取所述第二比值的平均值。其中，所述第二对称条件可根据实际情况需要进行设置，比如，所述第二对称条件可设置为所述第二比值的平均值大于1-a2且小于a2，其中，a2可取值为1.3或1.4 等。而在所述第二比值的平均值满足预设第二对称条件时，说明基于所述第二波峰序列和所述第二波谷序列所确定的波形左右对称，否则说明基于所述第二波峰序列和所述第二波谷序列所确定的波形左右不对称。在所述第二比值的平均值不满足预设第二对称条件，即基于所述第二波峰序列和所述第二波谷序列所确定的波形左右不对称时，表示检测到呼吸谐波，即所述呼吸信号中包含呼吸谐波，否则确定所述呼吸信号中不包含呼吸谐波。如此，通过检测基于所述第二波峰序列和所述第二波谷序列所确定的波形是否对称而判断所述呼吸信号中是否包含呼吸谐波的检测结果，操作简单且快速，进一步提高了心率提取的准确性和速度。66.在一些实施例中，所述根据所述检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率，包括：67.根据所述检测结果确定所述呼吸信号中包含呼吸谐波时，对所述呼吸信号进行呼吸谐波抵消处理，并对呼吸谐波抵消处理后获得的第一心跳信号进行心率提取，以获得所述目标对象的心率。68.根据所述检测结果确定所述呼吸信号中不包含呼吸谐波时，对所述生命体征信号进行第二滤波处理，并对第二滤波处理后获得的第二心跳信号进行心率提取，以获得所述目标对象的心率。69.具体地，当根据所述检测结果确定所述呼吸信号中包含呼吸谐波时，说明需要消除呼吸谐波以提高心率提取的准确性，此时对所述呼吸信号进行呼吸谐波抵消处理，获得呼吸谐波抵消处理后的第一心跳信号，并对呼吸谐波抵消处理后获得的第一心跳信号进行心率提取，以获得所述目标对象的心率。当根据所述检测结果确定所述呼吸信号中不包含呼吸谐波时，说明无需进行消除呼吸谐波处理，此时可直接对所述生命体征信号进行第二滤波处理，获得第二心跳信号，并对第二滤波处理后获得的第二心跳信号进行心率提取，以获得所述目标对象的心率。其中，所述第二滤波处理可采用带通滤波处理等方式，且当所述第二滤波处理为带通滤波处理时，对应的截止频率可根据心跳信号的频率大小进行确定。需要说明的是，上述对所述第一心跳信号或所述第二心跳信号进行心率提取的方法可采用现有的心率提取方法，如自相关函数方法或者谱估计方法，在此不再详细描述。如此，当目标对象的生命体征信号中不包含呼吸谐波时，无需进行呼吸谐波抵消处理，节省了运算资源；而当目标对象的生命体征信号中包含呼吸谐波时，则进行呼吸谐波抵消处理，以提高心率提取的准确性。70.在一些实施例中，所述对所述呼吸信号进行呼吸谐波抵消处理，包括：71.将所述呼吸信号作为自适应滤波器的输入信号、所述去除基线与高频干扰的生命体征信号作为自适应滤波器的期望信号，通过自适应滤波算法不断调整和更新自适应滤波器的参数，获得呼吸谐波抵消后的第一心跳信号。72.计算机设备通过对所述呼吸信号进行自相关函数功率谱估计，可获得所述呼吸信号的功率谱函数，而该功率谱函数的最大峰值即为所述呼吸信号的基频。基于所述呼吸信号的基频，可建立呼吸谐波信号数学模型，以获取呼吸谐波信号。进而，将所述呼吸信号作为自适应滤波器的主输入、所述去除基线与高频干扰的生命体征信号作为自适应滤波器的期望参考输入、所述呼吸谐波信号作为自适应滤波器的输出，其中，所述自适应滤波器可采用最小均方等算法不断迭代和更新自适应滤波器参数，当所述去除基线与高频干扰的生命体征信号与所述呼吸谐波信号之差的平方最小时，此时所述去除基线与高频干扰的生命体征信号与所述呼吸谐波信号之差即为第一心跳信号。如此，通过自适应滤波对呼吸信号进行呼吸谐波抵消处理，以获得呼吸谐波抵消后的心跳信号，进而基于所述心跳信号提取目标对象的心率，进一步提高了心率提取的准确性。73.本技术实施例另一方面，请参阅图3，还提供一种心率提取装置，所述装置包括：获取模块11，用于获取目标对象的生命体征信号；呼吸谐波检测模块 12，用于检测所述生命体征信号中是否包含呼吸谐波，得到检测结果；心率提取模块13，用于根据所述检测结果对所述生命体征信号执行相应的心率提取操作，以获得所述目标对象的心率。74.在一些实施例中，所述呼吸谐波检测模块12包括：滤波子模块120，用于对所述生命体征信号进行第一滤波处理，得到呼吸信号；检测子模块121，用于检测所述呼吸信号中是否包含呼吸谐波，得到检测结果。75.在一些实施例中，所述滤波子模块120，具体用于：基于预设第一频率区间对获取的目标对象的生命体征信号进行第一带通滤波处理，得到去除基线与高频干扰的生命体征信号；基于预设第二频率区间对所述去除基线与高频干扰的生命体征信号进行第二带通滤波处理，得到呼吸信号。76.在一些实施例中，所述检测子模块121，具体用于：77.对所述呼吸信号进行峰值检测，得到第一波峰序列和第一波谷序列；78.按照预设规则分别对所述第一波峰序列和所述第一波谷序列进行处理，得到第二波峰序列和第二波谷序列；79.判断基于所述第二波峰序列和所述第二波谷序列所确定的波形是否满足预设对称条件，得到所述呼吸信号中是否包含呼吸谐波的检测结果。80.在一些实施例中，所述检测子模块121，具体用于：81.删除所述第一波峰序列中幅值与最大波峰幅值的比值小于预设第一比值阈值的波峰、以及位置间隔小于预设第一距离阈值的相邻波峰中幅值小的波峰，并在位置间隔大于预设第二距离阈值的相邻波峰中插入一目标波峰，得到第二波峰序列；其中，所述第二距离阈值大于所述第一距离阈值，所述目标波峰的幅值为相邻波峰幅值的平均值，且所述目标波峰的位置为相邻波峰位置的平均值；82.删除所述第一波谷序列中幅值与最大波谷幅值的比值小于预设第二比值阈值的波谷、以及位置间隔小于预设第三距离阈值的相邻波谷中幅值小的波谷，并在位置间隔大于预设第四距离阈值的相邻波谷中插入一目标波谷，得到第二波谷序列；其中，所述第四距离阈值大于所述第三距离阈值，所述目标波谷的幅值为相邻波谷幅值的平均值，且所述目标波谷的位置为相邻波谷位置的平均值。83.在一些实施例中，所述检测子模块121，具体用于：84.获取基于所述第二波峰序列和所述第二波谷序列所确定的波形中每个周期内波峰幅值与波谷幅值的第一比值；其中，相邻的波峰和波谷组成一个周期；85.根据所述每个周期内波峰幅值与波谷幅值的第一比值，获取所述第一比值的平均值；86.在所述第一比值的平均值不满足预设第一对称条件时，确定所述呼吸信号中包含呼吸谐波，否则确定所述呼吸信号中不包含呼吸谐波。87.在一些实施例中，所述检测子模块121，具体用于：88.获取基于所述第二波峰序列和所述第二波谷序列所确定的波形中每个周期内第一幅值与第二幅值的第二比值；其中，相邻的波峰和波谷组成一个周期，所述第一幅值为波峰位置与波谷位置的位置差相对应的幅值，所述第二幅值为波谷位置与波峰位置的位置差相对应的幅值；89.根据所述每个周期内第一幅值与第二幅值的第二比值，获取所述第二比值的平均值；90.在所述第二比值的平均值不满足预设第二对称条件时，确定所述呼吸信号中包含呼吸谐波，否则确定所述呼吸信号中不包含呼吸谐波。91.在一些实施例中，所述心率提取模块13，具体用于：92.根据所述检测结果确定所述呼吸信号中包含呼吸谐波时，对所述呼吸信号进行呼吸谐波抵消处理，并对呼吸谐波抵消处理后获得的第一心跳信号进行心率提取，以获得所述目标对象的心率。93.根据所述检测结果确定所述呼吸信号中不包含呼吸谐波时，对所述生命体征信号进行第二滤波处理，并对第二滤波处理后获得的第二心跳信号进行心率提取，以获得所述目标对象的心率。94.在一些实施例中，所述心率提取模块13，具体用于：95.将所述呼吸信号作为自适应滤波器的输入信号、所述去除基线与高频干扰的生命体征信号作为自适应滤波器的期望信号，通过自适应滤波算法不断调整和更新自适应滤波器的参数，获得呼吸谐波抵消后的第一心跳信号。96.在一些实施例中，所述心率提取模块13，具体用于：97.对所述呼吸信号进行自相关函数功率谱估计，将获得的所述呼吸信号的功率谱函数的最大峰值确定为所述呼吸信号的基频；98.基于所述呼吸信号的基频，建立呼吸谐波信号数学模型，以获取呼吸谐波信号；99.将所述呼吸信号作为自适应滤波器的主输入、所述去除基线与高频干扰的生命体征信号作为自适应滤波器的期望参考输入、所述呼吸谐波信号作为自适应滤波器的输出，采用最小均方算法不断迭代和更新自适应滤波器参数，直至当所述去除基线与高频干扰的生命体征信号与所述呼吸谐波信号之差的平方最小时，将所述去除基线与高频干扰的生命体征信号与所述呼吸谐波信号之差确定为第一心跳信号。100.需要说明的是：上述实施例提供的心率提取装置在实现心率提取过程中，仅以上述各程序模块的划分进行举例说明，在实际应用中，可以根据需要而将上述处理分配由不同的程序模块完成，即可将装置的内部结构划分成不同的程序模块，以完成以上描述的全部或者部分方法步骤。另外，上述实施例提供的心率提取装置与心率提取方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。101.为了能够对本技术实施例提供的心率提取方法具有更系统的理解，下面结合一具体示例进行说明，本示例中以所述心率提取装置为压电式传感器或云端服务器为例，请参阅图4，本示例提供的心率提取装置包括如下模块：滤波模块21、呼吸谐波检测模块22、呼吸谐波消除模块23和心率提取模块24；其中，102.首先，滤波模块21对于压电式传感器采集的一帧信号x[n]，n＝1…n(n为离散时间序号，n为一帧数据总长度)，将该信号分别通过截止频率为0.1～4hz 的巴特沃斯带通滤波器(不限于采用其它fir或iir带通滤波器)得到去除基线与高频干扰的生命体征信号xm[n]，以及通过截止频率为0.1～0.5hz的巴特沃斯带通滤波器(不限于采用其它fir或iir带通滤波器)得到呼吸信号xb[n]。[0103]需要说明的是，所述滤波模块21也可通过小波分解、小波去噪等方法对原始信号x[n]进行处理，以得到去除基线与高频干扰的生命体征信号xm[n]和呼吸信号xb[n]等。[0104]对于呼吸谐波检测模块22，所述呼吸谐波检测模块22的作用主要包括峰值检测、删峰与补峰以及对称判断。[0105]关于峰值检测，在呼吸信号xb[n]中顺序查找值既大于左侧，又大于右侧的点为波峰，保存波峰值xp[n]，n＝1λnp和波峰位置lp[n]，n＝1λnp，其中np为波峰数目；并在呼吸信号xb[n]顺序查找值既小于左侧，又小于右侧的点为波谷，保存波谷值xv[n]，n＝1λnv和波谷位置lv[n]，n＝1λnv，其中nv为波谷数目。[0106]关于删峰与补峰，首先求解波峰模值最大值，在波峰点中删除模值小于最大值一半的波峰；并求解波谷模值最大值，在波谷点中删除模值小于最大值一半的波谷；接着，根据相邻两个波峰或波谷之间的位置差的合理范围为最小值lmin和最大值lmax，删除间隔小于lmin的相邻两个波峰中模值较小的波峰和删除间隔小于lmin的相邻两个波谷中模值较小的波谷，在间隔大于lmax的相邻两个波峰之间插入新的波峰，新波峰值为两个波峰值的平均值，新波峰位置为两个波峰位置的平均值，并在间隔大于lmax的相邻两个波谷之间插入新的波谷，新波谷值为两个波谷值的平均值，新波谷位置为两个波谷位置的平均值。经过删峰和补峰处理后，新的波峰值序列为xp’[n]，n＝1λnp’，对应位置序列为lp’[n]，n＝1λnp’，且新的波谷值序列为xv’[n]，n＝1λnv’，对应位置序列为lv’[n]，n＝1λnv’。[0107]关于对称判断，相邻波峰和波谷组成一个完整的周期，计算每一个周期内波峰模值与波谷模值比值τ1，计算一帧数据中所有完整周期的比值平均值若大于α1或小于1-α1(α1值可取1.3等)，则表示波形上下不对称，若大于 1-α1并小于α1，则表示波形上下对称；计算每一个周期内波峰与波谷位置差模值和波谷与波峰位置差模值的比值τ2，计算一帧数据中所有完整周期的比值平均值若大于α2或小于1-α2(α2值可取1.3等)，则表示波形左右不对称，若大于1-α2并小于α2，则表示波形左右对称。[0108]若波形上下不对称或左右不对称，则表示检测到呼吸谐波；否则，表示无呼吸谐波。[0109]对于呼吸谐波消除模块23，所述呼吸谐波消除模块23的作用主要包括呼吸谐波信号数学建模、自适应滤波谐波抵消；其中，[0110]关于呼吸谐波信号数学建模，包括如下过程：[0111]首先对滤波模块21产生的呼吸信号xb[n]做自相关函数的功率谱估计，得到呼吸信号功率谱函数p(f)，找到该功率谱函数的最大峰值即为呼吸信号的基频f0；[0112]接着，建立呼吸谐波信号数学模型[0113][0114]其中，m为呼吸谐波阶数(成人呼吸频段在0.1～0.5hz，心跳频段在0.8～2hz，考虑呼吸谐波影响心跳频段信号，此处阶数可取2～20之间)，ai、和φi分别为呼吸第i阶谐波幅值、角频率和相位。通过三角函数公式变换，该数学模型的幅值和相位形式可转换为同相和正交分量si(n)和sq(n)与系数α(n)和β(n)，其中，si(n)＝[cos(ω1n)，cos(ω2n)，…，cos(ωmn)]，sq(n)＝[sin(ω1n)，sin(ω2n)，…，sin(ωmn)]，α(n)＝[α1，α2，…，αm]，β(n)＝[β1，β2，…，βm]。[0115]关于自适应滤波谐波抵消，包括如下过程：[0116]请参阅图5，为本示例提供的自适应滤波示意图，将滤波模块21得到的生命体征信号xm[n]当做期望信号，滤波模块21得到的呼吸信号xb[n]作为输入信号，呼吸谐波信号作为自适应滤波输出信号，通过期望信号与滤波输出信号的差采用lms(最小均方)算法(也可以采用其它lms类或rls递推最小二乘算法替代)通过不断迭代和更新自适应滤波器参数，当两者差值的平方达到稳态时，两信号的差即为心跳信号。[0117]其中，定义初始化n＝0时刻自适应滤波器参数αi和βi为零向量，采用lms 算法进行迭代和更新的步骤如下：[0118](1)在n时刻，根据期望信号xm[n]和呼吸谐波信号s(n)，计算误差信号 e(n)：e(n)＝xm(n)-s(n)。[0119](2)通过递归更新系数：[0120][0121][0122]其中，μ为迭代步长。[0123](3)重复上述步骤(1)和(2)，直至算法稳态为止，可得到呼吸谐波抵消后的心跳信号xh[n]。[0124]对于心率提取模块24，若没有检测到呼吸谐波，则可将原始信号x[n]通过截止频率为0.8～2hz的巴特沃斯带通滤波器(不限于采用其它fir或iir带通滤波器)得到心跳信号xh[n]；若检测到呼吸谐波，可得到得到呼吸谐波抵消后的心跳信号xh[n]。然后，对心跳信号xh[n]采用自相关函数方法或者谱估计方法即可得到较为准确的心率。[0125]如此，上述实施例提供的心率提取装置，能够消除呼吸高次谐波的影响，得到更为准确的心率；而对没有检测到呼吸谐波的信号，不需要进行谐波抵消处理；此外，可通过呼吸波形是否对称，为生命健康检测提供更多维度的参考。[0126]本技术实施例另一方面，请参阅图6，还提供一种电子设备，包括处理器 211、存储器212及存储在所述存储器212上并可被所述处理器211执行的计算机程序，所述计算机程序被所述处理器211执行时实现本技术任一实施例所提供的心率提取方法，且能达到相同的技术效果，为避免重复，这里不再赘述。所述电子设备可以是智能手机、个人计算机等。[0127]本技术实施例另一方面，请参阅图7，还提供一种压电式传感器，包括处理器111、存储器112及存储在所述存储器112上并可被所述处理器111执行的计算机程序，所述计算机程序被所述处理器111执行时实现本技术任一实施例所提供的心率提取方法，且能达到相同的技术效果，为避免重复，这里不再赘述。所述压电式传感器可以包含与图6所示电子设备相同的硬件组成，且进一步包括采集目标对象的生命体征信号的采集传感组件。[0128]请参阅图8，是本技术实施例提供的电子设备的一种可选的硬件结构框图。如图8所示，电子设备可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器(processingunits，cpu)1110(处理器1110可以包括但不限于微处理器mcu或可编程逻辑器件fpga等的处理装置)、用于存储数据的存储器1130，一个或一个以上存储应用程序1123或数据1122的存储介质1120(例如一个或一个以上海量存储设备)。其中，存储器1130和存储介质1120可以是短暂存储或持久存储。存储在存储介质1120的程序可以包括一个或一个以上模块，每个模块可以包括对心率提取中的一系列指令操作。更进一步地，处理器 1110可以设置为与存储介质1120通信，在电子设备上执行存储介质1120中的一系列指令操作。电子设备还可以包括一个或一个以上电源1160，一个或一个以上有线或无线网络接口1150，一个或一个以上输入输出接口1140，和/或，一个或一个以上操作系统1121，例如windowsservertm，macosxtm， unixtm,linuxtm，freebsdtm等等。[0129]输入输出接口1140可以用于经由一个网络接收或者发送数据。上述的网络具体实例可包括电子设备的通信供应商提供的无线网络。在一个实例中，输入输出接口1140包括一个网络适配器(networkinterfacecontroller，nic)，其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中，输入输出接口1140可以为射频(radiofrequency，rf)模块，其用于通过无线方式与互联网进行通讯。[0130]本领域普通技术人员可以理解，图8所示的结构仅为示意，其并不对上述终端设备的结构造成限定。例如，电子设备还可包括比图8中所示更多或者更少的组件，或者具有与图8所示不同的配置。[0131]本技术实施例还提供一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述心率提取方法实施例的各个过程，且能达到相同的技术效果，为避免重复，这里不再赘述。其中，所述的计算机可读存储介质，如只读存储器(read-onlymemory，简称rom)、随机存取存储器(randomaccessmemory，简称ram)、磁碟或者光盘等。[0132]图9是本技术实施例提供的一种计算机可读存储介质800的结构框图。该计算机可读存储介质800中存储有程序代码，程序代码可被处理器调用执行本技术实施例所提供的智能家居控制方法。计算机可读存储介质800可以是诸如闪存、eeprom(电可擦除可编程只读存储器)、eprom、硬盘或rom之类的电子存储器。可选地，计算机可读存储介质800包括非瞬时性计算机可读介质(non-transitory computer-readable storage medium)。计算机可读存储介质800 具有执行上述方法中的任何方法步骤的程序代码810的存储空间。这些程序代码可以从一个或者多个计算机程序产品中读出或者写入到这一个或者多个计算机程序产品中。程序代码810可以例如以适当形式进行压缩。[0133]以上所述，仅为本技术的具体实施方式，但本技术的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本技术揭露的技术范围之内，可轻易想到变化或替换，都应涵盖在本技术的保护范围之内。因此，本技术的保护范围应以所述权利要求的保护范围为准。









图片声明：本站部分配图来自人工智能系统AI生成,觅知网授权图片,PxHere摄影无版权图库。本站只作为美观性配图使用,无任何非法侵犯第三方意图,一切解释权归图片著作权方,本站不承担任何责任。如有恶意碰瓷者,必当奉陪到底严惩不贷!




内容声明：本文中引用的各种信息及资料（包括但不限于文字、数据、图表及超链接等）均来源于该信息及资料的相关主体（包括但不限于公司、媒体、协会等机构）的官方网站或公开发表的信息。部分内容参考包括:(百度百科,百度知道,头条百科,中国民法典,刑法,牛津词典,新华词典,汉语词典,国家院校,科普平台)等数据,内容仅供参考使用,不准确地方联系删除处理！本站为非盈利性质站点,发布内容不收取任何费用也不接任何广告!




免责声明：我们致力于保护作者版权，注重分享，被刊用文章因无法核实真实出处，未能及时与作者取得联系，或有版权异议的，请联系管理员，我们会立即处理，本文部分文字与图片资源来自于网络，部分文章是来自自研大数据AI进行生成,内容摘自(百度百科,百度知道,头条百科,中国民法典,刑法,牛津词典,新华词典,汉语词典,国家院校,科普平台)等数据,内容仅供学习参考,不准确地方联系删除处理!的,若有来源标注错误或侵犯了您的合法权益，请立即通知我们，情况属实，我们会第一时间予以删除，并同时向您表示歉意,谢谢!