CN105662345B - 心跳信号处理方法、装置和系统 - Google Patents

心跳信号处理方法、装置和系统 Download PDF

Info

Publication number
CN105662345B
CN105662345B CN201610005670.1A CN201610005670A CN105662345B CN 105662345 B CN105662345 B CN 105662345B CN 201610005670 A CN201610005670 A CN 201610005670A CN 105662345 B CN105662345 B CN 105662345B
Authority
CN
China
Prior art keywords
waveform
heartbeat
curve
heart rate
rising edge
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610005670.1A
Other languages
English (en)
Other versions
CN105662345A (zh
Inventor
刘�文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Shuliantianxia Intelligent Technology Co Ltd
Original Assignee
Shenzhen H&T Intelligent Control Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shenzhen H&T Intelligent Control Co Ltd filed Critical Shenzhen H&T Intelligent Control Co Ltd
Priority to CN201610005670.1A priority Critical patent/CN105662345B/zh
Publication of CN105662345A publication Critical patent/CN105662345A/zh
Priority to PCT/CN2016/101249 priority patent/WO2017118127A1/zh
Application granted granted Critical
Publication of CN105662345B publication Critical patent/CN105662345B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/02Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024Detecting, measuring or recording pulse rate or heart rate
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16ZINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00Subject matter not provided for in other main groups of this subclass

Abstract

本发明公开了一种心跳信号处理方法,该方法包括:获取心跳信号原始波形,得到所述原始波形的波峰和波谷位置信息;根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形;根据所述心跳波形,计算得到心率值。本发明还公开了一种心跳信号处理装置和一种心跳信号处理系统。本发明通过拟合呼吸曲线,实现了去除原始波形上呼吸引起的干扰,解决了心跳信号受呼吸信号影响而导致误差大的问题,并根据心跳信号计算得到的了精准心率值。

Description

心跳信号处理方法、装置和系统
技术领域
本发明涉及信号处理技术领域,尤其涉及一种心跳信号处理方法、装置和系统。
背景技术
当前,睡眠监控主要是通过雷达红外、超声波检测、压电传感器,获取人们的心跳信号、呼吸信号等,进行处理得到人们睡眠时的心率值等各项睡眠参数。其中,雷达红外和超声检测手段成本较高,并且会产生红外或者超声辐射,检测效果也不理想。
通过压电传感器监测睡眠参数的原理是,通过灵敏度较高的压电薄膜传感器监测受测者由呼吸引起的胸腔收扩运动和心跳震颤引起的机械信号,进而分析信号得到呼吸、心跳等参数。
因而,使用压电传感器监测心率时,容易受到各种信号干扰,例如:较强烈的呼吸干扰、人体生理信号(包括体动等)干扰、外界环境干扰等。其中,呼吸干扰影响最大,因为呼吸信号较为强烈,基本淹没了心跳信号。另一方面,加上体动和外界环境干扰,以及同心率频率段的信号干扰,使得原本就微弱的心跳信号更加弱小,加大处理难度。如何在干扰下进行心跳信号处理,得到有效的心跳信号,成为当前心率监测的重要研究课题。
发明内容
本发明的主要目的在于提供一种心跳信号处理方法、装置和系统,旨在解决心跳信号处理准确性差、误差大的技术问题。
为实现上述目的,本发明提供一种心跳信号处理方法,所述心跳信号处理方法包括以下步骤:
获取心跳信号原始波形,得到所述原始波形的波峰和波谷位置信息;
根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形;
根据所述心跳波形,计算得到心率值。
优选的,所述根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形的步骤包括:
根据所述原始波形的波峰和波谷位置信息,得到波形趋势;
根据所述原始波形的波峰和波谷位置信息、所述波形趋势,拟合呼吸波形得到呼吸曲线;
根据所述呼吸曲线,去除所述呼吸曲线对所述原始波形的干扰,并进行滤波,得到心跳波形。
优选的,所述呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线,所述心跳波形包括上升沿心跳曲线和下降沿心跳曲线,所述根据所述呼吸曲线,去除所述呼吸曲线对所述原始波形的干扰,并进行滤波,得到心跳波形的步骤包括:
根据所述原始波形的波峰和波谷位置信息、所述波形趋势,得到原始波形的上升沿波形曲线和下降沿波形曲线;
在所述上升沿波形曲线上去除所述上升沿呼吸曲线,并进行滤波,得到上升沿心跳曲线;
在所述下降沿波形曲线上去除所述下降沿呼吸曲线,并进行滤波,得到下降沿心跳曲线。
优选的,所述根据所述心跳波形,计算得到心率值的步骤包括:
分别根据所述上升沿心跳曲线和所述下降沿心跳曲线,计算得到对应的上升沿心率和下降沿心率;
根据预设的处理规则,去除所述上升沿心率和下降沿心率中的异常数据;
根据所述处理后的上升沿心率和下降沿心率,得到心率值。
此外,为实现上述目的,本发明还提供一种心跳信号处理装置,所述心跳信号处理装置包括:
获取模块,用于获取心跳信号原始波形,得到所述原始波形的波峰和波谷位置信息;
波形处理模块,用于根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形;
计算模块,用于根据所述心跳波形,计算得到心率值。
优选的,所述波形处理模块包括:
波形趋势单元,用于根据所述原始波形的波峰和波谷位置信息,得到波形趋势;
拟合单元,用于根据所述原始波形的波峰和波谷位置信息、所述波形趋势,拟合呼吸波形得到呼吸曲线;
去干扰单元,用于根据所述呼吸曲线,去除所述呼吸曲线对所述原始波形的干扰,并进行滤波,得到心跳波形。
优选的,所述呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线,所述心跳波形包括上升沿心跳曲线和下降沿心跳曲线,所述去干扰单元包括:
波形定位子单元,用于根据所述原始波形的波峰和波谷位置信息、所述波形趋势,得到原始波形的上升沿波形曲线和下降沿波形曲线;
去除子单元,用于在所述上升沿波形曲线上去除所述上升沿呼吸曲线,并进行滤波,得到上升沿心跳曲线;在所述下降沿波形曲线上去除所述下降沿呼吸曲线,并进行滤波,得到下降沿心跳曲线。
优选的,所述计算模块包括:
计算单元,用于分别根据所述上升沿心跳曲线和所述下降沿心跳曲线,计算得到对应的上升沿心率和下降沿心率;
异常处理单元,用于根据预设的处理规则,去除所述上升沿心率和下降沿心率中的异常数据;
综合处理单元,用于根据所述处理后的上升沿心率和下降沿心率,得到心率值。
此外,为实现上述目的,本发明还提供一种心跳信号处理系统,其特征在于,所述心跳信号处理系统包括压电传感器、模数转换器和心率监测CPU,其中:
所述压电传感器用于,采集压电模拟信号;
所述模数转换器用于,将所述压电模拟信号转换为数字信号,得到心跳信号原始波形;
所述心率监测CPU,包括获取模块、波形处理模块和计算模块。
本发明实施例提出的一种心跳信号处理方法、装置和系统,通过获取心跳信号原始波形,得到原始波形的波峰和波谷位置信息;然后,根据原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并根据呼吸曲线对原始波形去除干扰,得到心跳波形;然后,根据得到的心跳波形,计算得到心率值。本发明实施例通过拟合呼吸曲线,实现了去除原始波形上呼吸引起的干扰,解决了心跳信号受呼吸信号影响而导致误差大的问题,根据心跳信号计算得到的了精准心率值。
附图说明
图1为本发明心跳信号处理方法第一实施例的流程示意图;
图2为本发明心跳信号处理方法第二实施例的流程示意图;
图3为本发明心跳信号处理方法第三实施例的流程示意图;
图4为本发明实施例的一种心跳信号原始波形处理流程示意图;
图5为本发明心跳信号处理方法第四实施例的流程示意图;
图6为本发明实施例中的一种心跳波形示意图;
图7为本发明心跳信号处理装置第一实施例的功能模块示意图;
图8为本发明心跳信号处理装置第二实施例的功能模块示意图;
图9为本发明心跳信号处理装置第三实施例的功能模块示意图;
图10为本发明心跳信号处理装置第四实施例的功能模块示意图;
图11为本发明心跳信号处理系统第一实施例的模块示意图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明实施例的主要解决方案是:获取心跳信号原始波形,得到原始波形的波峰和波谷位置信息;然后,根据原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并根据呼吸曲线对原始波形去除干扰,得到心跳波形;然后,根据得到的心跳波形,计算得到心率值。
由于现有技术不能有效的去除呼吸信号的干扰,根据心跳信号得到的心率值误差大。
本发明提供一种解决方案,通过拟合呼吸曲线,实现了去除原始波形上呼吸引起的干扰,解决了心跳信号受呼吸信号影响而导致误差大的问题,根据心跳信号计算得到的了精准心率值。
参照图1,本发明心跳信号处理方法第一实施例提供一种心跳信号处理方法,所述心跳信号处理方法包括:
步骤S10、获取心跳信号原始波形,得到所述原始波形的波峰和波谷位置信息。
本实施例主要用于心跳信号的处理,得到心率。
具体的,作为一种实施方式,本实施例通过压电传感器获取用户呼吸和心跳的身体振动信号,例如,将压电传感器置于床垫中或床单下,获取人体睡觉时,呼吸引起的胸腔收括振动信号和心跳引起的有规律的身体震动信号。
然后,压电传感器将得到的模拟信号发送给A/D转换器(Analog to DigitalConverter,模数转换器)。
A/D转换器收到模拟信号后,将模拟线号转换为数字信号,得到心跳信号原始波形,并将得到的心跳信号原始波形发送给心率监测CPU,由心率监测CPU进行心跳信号的处理。
心率监测CPU接收A/D转换器发送的采样波(A/D Sample),得到心跳信号的原始波形后,首先获取原始波形的波峰和波谷位置信息。
作为一种实施方式,心率监测CPU可直接进行差分阈值,得到原始波形的波峰(Sample-Peak)和波谷(Sample-Trough)位置信息并存储。
原始波形的波峰位置信息包括原始波形中各波峰的位置和振幅,原始波形的波谷位置信息包括原始波形中各波谷的位置和振幅。
作为一种实施方式,首先,心率监测CPU通过高低通滤波,得到呼吸波形。
需要说明的是,由于人体正常的呼吸周期约为1.5秒到5s一次,因此,可根据呼吸周期灵活设置高低通滤波的频率截止范围。例如,本实施例设置呼吸波形的频率范围为0.1—0.5Hz,则通过高低通滤波得到的波形即为0.1—0.5Hz范围内的波形。
在得到呼吸波形后,心率监测CPU进行差分阈值,得到呼吸波形的波峰(Breath-Peak)和波谷(Breath-Trough)位置信息。
然后,分别根据呼吸波形的波峰和波谷的位置,查找原始波形上对应的位置的波峰和波谷。
例如,以当前呼吸波形波峰的位置为基准,设置区间范围,查找原始波形对应位置区间范围内的振幅最大点,以此点为原始波形对应位置的波峰。同样的,以当前呼吸波形波谷的位置为基准,设置区间范围,查找原始波形对应位置区间范围内的振幅最小点,以此点为原始波形对应位置的峰谷。
根据呼吸波形的波峰和波谷的位置,查找原始波形上对应的位置的波峰和波谷,可以避免直接获取原始波形的波峰和波谷时,由于干扰因素导致获取的波峰、波谷位置不是呼吸波形的准确位置,从而导致拟合的呼吸曲线不准确
需要说明的是,预设的区间范围可以根据实际需要灵活设置,例如取基准位置前后2秒内采集到的数据点为预设的区间范围。
由此,心率监测CPU根据呼吸波形的波峰和波谷位置信息得到了原始波形的波峰和波谷位置信息。
步骤S20、根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形。
在获取原始波形的波峰和波谷位置信息后,心率监测CPU根据得到的原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除呼吸曲线对原始波形的干扰,得到心跳波形。
具体的,作为一种实施方式,首先,心率监测CPU根据原始波形的波峰和波谷位置信息,判断波形趋势。
例如,根据得到的原始波形的波峰和波谷位置信息,若原始波形上,一个波峰的前一个点为波谷,则次波峰的波形趋势为上升沿(Wave-S);一个波谷的前一个点为波峰,则次波峰的波形趋势为下降沿(Wave-X)。
然后,根据原始波形相邻波峰和波谷的位置信息,及波形趋势,在原始波形上拟合呼吸波形,得到呼吸曲线。
然后,根据呼吸曲线去除原始波形上的呼吸干扰,得到去除呼吸干扰后的波形。
对得到的去除呼吸干扰后波形进行滤波,根据人体正常心跳周期,设置心跳波形的频率范围为0.8—1.6Hz,则通过高低通滤波得到的波形即为0.8—1.6Hz范围内的波形。
由此,得到心跳波形。
步骤S30、根据所述心跳波形,计算得到心率值。
在得到心跳波形后,心率监测CPU根据心跳波形计算心率值。
具体的,作为一种实施方式,首先,去除心跳波形的开始部分和结束部分,用于去除开始部分和结束部分采集的不规律心跳信号。
然后,根据心跳波形的相邻波峰采样点距离得到两个波峰的采样间隔时间,根据两个波峰的间隔时间,计算得到采样时间节点的心率值。同理,也可以根据波谷的位置信息计算得到心率值。
然后,对计算得到心率值进行取舍、范围限制等优化。例如,根据人体正常心率值,设置心率范围为每分钟40-120次,根据此范围去掉由起床等原因引起的不合理的过快心跳或慢心跳;若得到的某个采样时间节点的心率值,与此时间节点的前一个时间节点得到的心率值相比,小于前一个心率值的三分之一或大于前一个心率值的三分之四,则判定得到的心率值浮动过大,去除此心率值。
然后,平滑处理得到的心率值,取连续的数个心率值的平均值,得到的心率值即为采样时间内的心率值。
在本实施例中,心率监测CPU获取心跳信号原始波形,得到原始波形的波峰和波谷位置信息;然后,根据原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并根据呼吸曲线对原始波形去除干扰,得到心跳波形;然后,根据得到的心跳波形,计算得到心率值。本实施例通过拟合呼吸曲线,实现了去除原始波形上呼吸引起的干扰,解决了心跳信号受呼吸信号影响而导致误差大的问题,根据心跳信号计算得到了精准心率值。
进一步的,参照图2,本发明心跳信号处理方法第二实施例提供一种心跳信号处理方法,基于上述图1所示的实施例,所述步骤S20包括:
步骤S21、根据所述原始波形的波峰和波谷位置信息,得到波形趋势。
在获取原始波形的波峰和波谷位置信息后,心率监测CPU根据得到的原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除呼吸曲线对原始波形的干扰,得到心跳波形。
作为另一种实施方式,首先,心率监测CPU根据原始波形的波峰和波谷位置信息,判断波形趋势,定位原始波形的上升沿和下降沿。
例如,根据得到的原始波形的波峰和波谷位置信息,若原始波形上,一个波峰的前一个点为波谷,则此波峰的波形趋势为上升沿;一个波谷的前一个点为波峰,则此波谷的波形趋势为下降沿。
步骤S22、根据所述原始波形的波峰和波谷位置信息、所述波形趋势,拟合得到呼吸曲线。
在获取原始波形的波峰和波谷位置信息后,心率监测CPU根据原始波形的波峰和波谷位置信息、波形趋势分别拟合呼吸波形。
将原始波形相邻波峰和波谷位置,作为起始点和结束点,对上升沿波形拟合呼吸波形,得到上升沿呼吸曲线;根据下降沿波形拟合呼吸波形,得到下降沿呼吸曲线。
由此,连接相邻的、有一端点相同的上升沿呼吸曲线和下降沿呼吸曲线,得到原始波形上拟合的呼吸曲线。
步骤S23、根据所述呼吸曲线,去除所述呼吸曲线对所述原始波形的干扰,并进行滤波,得到心跳波形。
在得到拟合的呼吸曲线后,心率监测CPU去除呼吸曲线对原始波形的干扰,得到心跳波形。
具体的,作为一种实施方式,由于得到的原始波形与呼吸曲线位置对应,在原始波形上,减去对应位置的呼吸曲线,通过波形上相同位置,采样点振幅的相减,进行呼吸曲线的去除或削弱,得到的波形即为去除呼吸干扰后的波形。
然后,对得到的去除呼吸干扰后波形进行高低通滤波,根据人体正常心跳周期,设置心跳波形的频率范围为0.8—1.6Hz,则通过高低通滤波得到的波形即为0.8—1.6Hz范围内的波形。
由此,得到心跳波形。
在本实施例中,心率监测CPU根据原始波形的波峰和波谷位置信息,得到波形趋势;然后,根据原始波形的波峰和波谷位置信息、波形趋势,拟合呼吸波形得到呼吸曲线;然后,根据得到的呼吸曲线,去除呼吸曲线对原始波形的干扰,得到心跳波形。本实施例根据原始波形的波峰和波谷位置信息得到波形趋势,拟合得到更准确的呼吸曲线,然后将原始波形减去对应位置的呼吸曲线去除呼吸干扰,并进行滤波去除杂波干扰,得到准确的、无干扰的心跳波形。
进一步的,参照图3,本发明心跳信号处理方法第三实施例提供一种心跳信号处理方法,基于上述图2所示的实施例,所述呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线,所述心跳波形包括上升沿心跳曲线和下降沿心跳曲线,所述步骤S23包括:
步骤S231、根据所述原始波形的波峰和波谷位置信息、所述波形趋势,得到原始波形的上升沿波形曲线和下降沿波形曲线。
在得到拟合的呼吸曲线后,心率监测CPU去除呼吸曲线对原始波形的干扰,得到心跳波形。
具体的,作为一种实施方式,心率监测CPU根据原始波形的波峰和波谷位置信息,针对上升沿波形和下降沿波形分别拟合,拟合得到的呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线。
心率监测CPU根据原始波形的波峰和波谷位置信息、波形趋势,获取原始波形的上升沿波形曲线和下降沿波形曲线。例如,若相邻两点的波形趋势为上升沿波形,则根据这两点的位置信息,得到这段原始波形为上升沿呼吸曲线。
根据波峰和波谷位置,原始波形的上升沿波形曲线与拟合得到的上升沿呼吸曲线对应,原始波形的下降沿波形曲线与拟合得到的下降沿呼吸曲线。
步骤S232、在所述上升沿波形曲线上去除所述上升沿呼吸曲线,并进行滤波,得到上升沿心跳曲线;在所述下降沿波形曲线上去除所述下降沿呼吸曲线,并进行滤波,得到下降沿心跳曲线。
心率监测CPU去除原始波形的呼吸信号干扰时,具体的,作为一种实施方式,在原始波形上,心率监测CPU分别对上升沿波形曲线和下降沿波形曲线进行处理。
参照图4,取上升沿波形曲线为Ss,下降沿波形曲线为Sx,上升沿呼吸曲线为ys,下降沿呼吸曲线为wx。由于原始波形与呼吸曲线位置对应,上升沿波形曲线Ss对应上升沿呼吸曲线ys,下降沿波形曲线Sx对应下降沿呼吸曲线wx。
则上升沿波形曲线减去上升沿呼吸曲线,得到曲线Hs,曲线Hs=Ss-ys。曲线Hs上点的振幅,为相同位置的上升沿波形曲线点与上升沿呼吸曲线点的振幅之差。
然后,对曲线Hs进行高低通滤波,根据人体正常心跳周期,设置心跳波形的频率范围为0.8—1.6Hz,则通过高低通滤波得到的波形即为0.8—1.6Hz范围内的上升沿心跳曲线。
同理,下降沿波形曲线减去下降沿呼吸曲线,得到曲线Hx,曲线Hx=Sx-wx。下降沿心跳曲线上点的振幅,为相同位置的下降沿波形曲线点与下降沿呼吸曲线点的振幅之差。
然后,对曲线Hx进行高低通滤波,根据人体正常心跳周期,设置心跳波形的频率范围为0.8—1.6Hz,则通过高低通滤波得到的波形即为0.8—1.6Hz范围内的上升沿心跳曲线。
由此,心率监测CPU得到上升沿心跳曲线和下降沿心跳曲线,也即得到心跳波形。
本实施例中,呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线,心率监测CPU首先得到原始波形的上升沿波形曲线和下降沿波形曲线;然后,在上升沿波形曲线上去除上升沿呼吸曲线,并进行滤波,得到上升沿心跳曲线;在下降沿波形曲线上去除下降沿呼吸曲线,并进行滤波,得到下降沿心跳曲线,由此,得到心跳波形。本实施例通过波形趋势,定位上升沿波形和下降沿波形,根据上升沿波形和下降沿波形分别拟合得到上升沿呼吸曲线和下降沿呼吸曲线,然后对原始波形的上升沿波形和下降沿波形的分别进行去干扰处理,避免呼吸信号的波峰与心跳信号的波峰一致而导致的去干扰产生的误差,实现了心跳信号的准确获取。
进一步的,参照图5,本发明心跳信号处理方法第四实施例提供一种心跳信号处理方法,基于上述图3所示的实施例,所述步骤S30包括:
步骤S31、分别根据所述上升沿心跳曲线和所述下降沿心跳曲线,计算得到对应的上升沿心率和下降沿心率。
在获取上升沿心跳曲线和下降沿心跳曲线后,心率监测CPU分别计算心率值。
具体的,作为一种实施方式,首先,心率监测CPU获取上升沿心跳曲线和下降沿心跳曲线的波峰、波谷位置信息。
然后,心率监测CPU根据上升沿心跳曲线相邻两个波峰或波谷采样点位置的距离,得到两个波峰或波谷的采样间隔时间,根据两个波峰或波谷的间隔时间,计算得到采样时间节点的心率值。由此,得到上升沿心跳曲线对应的上升沿心率值HR1。
同理,心率监测CPU根据下降沿心跳曲线相邻两个波峰或波谷采样点位置的距离,得到两个波峰或波谷的采样间隔时间,根据两个波峰或波谷的间隔时间,计算得到采样时间节点的心率值。由此,得到下降沿心跳曲线对应的下降沿心率值HR2。
步骤S32、根据预设的处理规则,去除所述上升沿心率和下降沿心率中的异常数据。
在得到上升沿心率值HR1和下降沿心率值HR2后,心率监测CPU根据预设的规则对得到心率值进行处理,去除异常数据。
具体的,作为一种实施方式,预设的处理规则包括取舍、范围限制、平滑处理。
首先,心率监测CPU舍去开始进行心跳信号采样时、一定时间范围内的初始波形,舍去结束心跳采样时、一定时间范围内的结束波形。初始波形和结束波形由于受到采样仪器的初始化等外部原因,会引起波形异常。
然后,根据预先设置的心率范围,去除过大或过小的心率值。根据人体正常心率值,设置心率范围为每分钟40-120次,根据此范围去掉由起床等原因引起的不合理的过快心跳或慢心跳。
另外,若得到的某个采样时间节点的心率值,与此时间节点的前一个时间节点得到的心率值相比,小于前一个心率值的三分之一或大于前一个心率值的三分之四,则判定得到的心率值浮动过大,去除此心率值。
另外,若一个上升沿心跳曲线或下降沿心跳曲线对应的心跳数量小于2个,则舍弃此上升沿心跳曲线或下降沿心跳曲线对应的心率值。
需要说明的是,预设的处理规则还可以包括其他异常数据处理方式,可根据实际需要灵活设置。
由此,心率监测CPU得到去除异常数据后的上升沿心率数据和下降沿心率数据。
步骤S33、根据所述处理后的上升沿心率和下降沿心率,得到心率值。
在得到去除异常数据后的上升沿心率数据和下降沿心率数据后,心率监测CPU综合处理上升沿心率和下降沿心率,得到原始波形对应的心率值。
具体的,作为一种实施方式,心率监测CPU对得到的心率值平均值,得到的心率值即为原始波形采样时间内的心率值。
或者,作为另一种实施方式,心率监测CPU将得到的上升沿心率和下降沿心率按对应的采样时间节点排序,并按时间顺序划分为数个采样区。然后,分别求各采样区的心率平均值,再综合处理得到原始波形采样时间内的心率值。
本发明实施例采用段处理,采集心跳信号时,将得到的预设长度的原始波形作为一段原始波形进行处理,计算得到心率此段原始波形对应的心率值。其中,预设的长度可根据实际需要灵活设置。
例如,每采集70个点后,将这70个点的波形作为一段原始波形进行处理,得到对应的心跳波形。参照图6所示的心跳波形图,心率监测CPU得到多段心跳波形。
由于开始心跳信号采样和结束心跳信号采样时,心跳信号误差较大,因此心率监测CPU得到多段心跳波形后,去除开始采样和结束采样时预设时间范围内的心跳波形。
由此,心率监测CPU计算得到多个心率值。
然后,心率监测CPU对得到的多个心率值进行综合处理,去除异常数据,计算多段心跳波形的平均心率值,得到最终的心率值。
在本实施例中,心率监测CPU分别根据上升沿心跳曲线和下降沿心跳曲线,计算得到对应的上升沿心率和下降沿心率;然后,根据预设的处理规则,去除上升沿心率和下降沿心率中的异常数据;然后,根据处理后的上升沿心率和下降沿心率,得到原始波形对应的心率值。本实施例根据上升沿心跳曲线和下降沿心跳曲线分别计算得到心率之后,在对得到的心率值进行综合处理得到最终的心率值,实现了对心跳信号的处理,得到心率值。
参照图7,本发明心跳信号处理装置第一实施例提供一种心跳信号处理装置,所述心跳信号处理装置包括:
获取模块100,用于获取心跳信号原始波形,得到所述原始波形的波峰和波谷位置信息。
本实施例主要用于心跳信号的处理,得到心率。
具体的,作为一种实施方式,本实施例通过压电传感器获取用户呼吸和心跳的身体振动信号,例如,将压电传感器置于床垫中或床单下,获取人体睡觉时,呼吸引起的胸腔收括振动信号和心跳引起的有规律的身体震动信号。
然后,压电传感器将得到的模拟信号发送给A/D转换器(Analog to DigitalConverter,模数转换器)。
A/D转换器收到模拟信号后,将模拟线号转换为数字信号,得到心跳信号原始波形,并将得到的心跳信号原始波形发送给获取模块100,由获取模块100进行心跳信号的处理。
获取模块100接收A/D转换器发送的采样波(A/D Sample),得到心跳信号的原始波形后,首先获取原始波形的波峰和波谷位置信息。
作为一种实施方式,获取模块100可直接进行差分阈值,得到原始波形的波峰(Sample-Peak)和波谷(Sample-Trough)位置信息并存储。
原始波形的波峰位置信息包括原始波形中各波峰的位置和振幅,原始波形的波谷位置信息包括原始波形中各波谷的位置和振幅。
作为一种实施方式,首先,获取模块100通过高低通滤波,得到呼吸波形。
需要说明的是,由于人体正常的呼吸周期约为1.5秒到5s一次,因此,可根据呼吸周期灵活设置高低通滤波的频率截止范围。例如,本实施例设置呼吸波形的频率范围为0.1—0.5Hz,则通过高低通滤波得到的波形即为0.1—0.5Hz范围内的波形。
在得到呼吸波形后,心率监测CPU进行差分阈值,得到呼吸波形的波峰(Breath-Peak)和波谷(Breath-Trough)位置信息。
然后,分别根据呼吸波形的波峰和波谷的位置,查找原始波形上对应的位置的波峰和波谷。
例如,以当前呼吸波形波峰的位置为基准,设置区间范围,查找原始波形对应位置区间范围内的振幅最大点,以此点为原始波形对应位置的波峰。同样的,以当前呼吸波形波谷的位置为基准,设置区间范围,查找原始波形对应位置区间范围内的振幅最小点,以此点为原始波形对应位置的峰谷。
根据呼吸波形的波峰和波谷的位置,查找原始波形上对应的位置的波峰和波谷,可以避免直接获取原始波形的波峰和波谷时,由于干扰因素导致获取的波峰、波谷位置不是呼吸波形的准确位置,从而导致拟合的呼吸曲线不准确
需要说明的是,预设的区间范围可以根据实际需要灵活设置,例如取基准位置前后2秒内采集到的数据点为预设的区间范围。
由此,获取模块100根据呼吸波形的波峰和波谷位置信息得到了原始波形的波峰和波谷位置信息。
波形处理模块200,用于根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形。
在获取原始波形的波峰和波谷位置信息后,波形处理模块200根据得到的原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除呼吸曲线对原始波形的干扰,得到心跳波形。
具体的,作为一种实施方式,首先,波形处理模块200根据原始波形的波峰和波谷位置信息,判断波形趋势。
例如,根据得到的原始波形的波峰和波谷位置信息,若原始波形上,一个波峰的前一个点为波谷,则次波峰的波形趋势为上升沿(Wave-S);一个波谷的前一个点为波峰,则次波峰的波形趋势为下降沿(Wave-X)。
然后,根据原始波形相邻波峰和波谷的位置信息,及波形趋势,在原始波形上拟合呼吸波形,得到呼吸曲线。
然后,根据呼吸曲线去除原始波形上的呼吸干扰,得到去除呼吸干扰后的波形。
对得到的去除呼吸干扰后波形进行滤波,根据人体正常心跳周期,设置心跳波形的频率范围为0.8—1.6Hz,则通过高低通滤波得到的波形即为0.8—1.6Hz范围内的波形。
由此,得到心跳波形。
计算模块300,用于根据所述心跳波形,计算得到心率值。
在得到心跳波形后,计算模块300根据心跳波形计算心率值。
具体的,作为一种实施方式,首先,去除心跳波形的开始部分和结束部分,用于去除开始部分和结束部分采集的不规律心跳信号。
然后,根据心跳波形的相邻波峰采样点距离得到两个波峰的采样间隔时间,根据两个波峰的间隔时间,计算得到采样时间节点的心率值。同理,也可以根据波谷的位置信息计算得到心率值。
然后,对计算得到心率值进行取舍、范围限制等优化。例如,根据人体正常心率值,设置心率范围为每分钟40-120次,根据此范围去掉由起床等原因引起的不合理的过快心跳或慢心跳;若得到的某个采样时间节点的心率值,与此时间节点的前一个时间节点得到的心率值相比,小于前一个心率值的三分之一或大于前一个心率值的三分之四,则判定得到的心率值浮动过大,去除此心率值。
然后,平滑处理得到的心率值,取连续的数个心率值的平均值,得到的心率值即为采样时间内的心率值。
本实施例中,获取模块100获取心跳信号原始波形,得到原始波形的波峰和波谷位置信息;然后,波形处理模块200根据原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并根据呼吸曲线对原始波形去除干扰,得到心跳波形;然后,计算模块300根据得到的心跳波形,计算得到心率值。本实施例通过拟合呼吸曲线,实现了去除原始波形上呼吸引起的干扰,解决了心跳信号受呼吸信号影响而导致误差大的问题,本实施例根据心跳信号计算得到的心率值非常精准,接近于实际心率值。
进一步的,参照图8,本发明心跳信号处理装置第二实施例提供一种心跳信号处理装置,基于上述图7所示的实施例,所述波形处理模块200包括:
波形趋势单元210,用于根据所述原始波形的波峰和波谷位置信息,得到波形趋势。
在获取原始波形的波峰和波谷位置信息后,波形处理模块200根据得到的原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除呼吸曲线对原始波形的干扰,得到心跳波形。
作为另一种实施方式,首先,波形趋势单元210根据原始波形的波峰和波谷位置信息,判断波形趋势,定位原始波形的上升沿和下降沿。
例如,根据得到的原始波形的波峰和波谷位置信息,若原始波形上,一个波峰的前一个点为波谷,则此波峰的波形趋势为上升沿;一个波谷的前一个点为波峰,则此波谷的波形趋势为下降沿。
拟合单元220,用于根据所述原始波形的波峰和波谷位置信息、所述波形趋势,拟合呼吸波形得到呼吸曲线。
在获取原始波形的波峰和波谷位置信息后,拟合单元220根据原始波形的波峰和波谷位置信息、波形趋势分别拟合呼吸波形。
将原始波形相邻波峰和波谷位置,作为起始点和结束点,对上升沿波形拟合呼吸波形,得到上升沿呼吸曲线;根据下降沿波形拟合呼吸波形,得到下降沿呼吸曲线。
由此,连接相邻的、有一端点相同的上升沿呼吸曲线和下降沿呼吸曲线,得到原始波形上拟合的呼吸曲线。
去干扰单元230,用于根据所述呼吸曲线,去除所述呼吸曲线对所述原始波形的干扰,并进行滤波,得到心跳波形。
在得到拟合的呼吸曲线后,去干扰单元230去除呼吸曲线对原始波形的干扰,得到心跳波形。
具体的,作为一种实施方式,由于得到的原始波形与呼吸曲线位置对应,在原始波形上,减去对应位置的呼吸曲线,通过波形上相同位置,采样点振幅的相减,进行呼吸曲线的去除或削弱,得到的波形即为去除呼吸干扰后的波形。
然后,对得到的去除呼吸干扰后波形进行高低通滤波,根据人体正常心跳周期,设置心跳波形的频率范围为0.8—1.6Hz,则通过高低通滤波得到的波形即为0.8—1.6Hz范围内的波形。
由此,得到心跳波形。
在本实施例中,波形趋势单元210根据原始波形的波峰和波谷位置信息,得到波形趋势;然后,拟合单元220根据原始波形的波峰和波谷位置信息、波形趋势,拟合呼吸波形得到呼吸曲线;然后,去干扰单元230根据得到的呼吸曲线,去除呼吸曲线对原始波形的干扰,得到心跳波形。本实施例根据原始波形的波峰和波谷位置信息得到波形趋势,拟合得到更准确的呼吸曲线,然后将原始波形减去对应位置的呼吸曲线去除呼吸干扰,并进行滤波去除杂波干扰,得到准确的、无干扰的心跳波形。
进一步的,参照图9,本发明心跳信号处理装置第三实施例提供一种心跳信号处理装置,基于上述图8所示的实施例,所述呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线,所述心跳波形包括上升沿心跳曲线和下降沿心跳曲线,所述去干扰单元230包括:
波形定位子单元231,用于根据所述原始波形的波峰和波谷位置信息、所述波形趋势,得到原始波形的上升沿波形曲线和下降沿波形曲线。
在得到拟合的呼吸曲线后,波形定位子单元231去除呼吸曲线对原始波形的干扰,得到心跳波形。
具体的,作为一种实施方式,波形定位子单元231根据原始波形的波峰和波谷位置信息,针对上升沿波形和下降沿波形分别拟合,拟合得到的呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线。
波形定位子单元231根据原始波形的波峰和波谷位置信息、波形趋势,获取原始波形的上升沿波形曲线和下降沿波形曲线。例如,若相邻两点的波形趋势为上升沿波形,则根据这两点的位置信息,得到这段原始波形为上升沿呼吸曲线。
根据波峰和波谷位置,原始波形的上升沿波形曲线与拟合得到的上升沿呼吸曲线对应,原始波形的下降沿波形曲线与拟合得到的下降沿呼吸曲线。
去除子单元232,用于在所述上升沿波形曲线上去除所述上升沿呼吸曲线,并进行滤波,得到上升沿心跳曲线;在所述下降沿波形曲线上去除所述下降沿呼吸曲线,并进行滤波,得到下降沿心跳曲线。
去除子单元232去除原始波形的呼吸信号干扰时,具体的,作为一种实施方式,在原始波形上,去除子单元232分别对上升沿波形曲线和下降沿波形曲线进行处理。
参照图4,取上升沿波形曲线为Ss,下降沿波形曲线为Sx,上升沿呼吸曲线为ys,下降沿呼吸曲线为wx。由于原始波形与呼吸曲线位置对应,上升沿波形曲线Ss对应上升沿呼吸曲线ys,下降沿波形曲线Sx对应下降沿呼吸曲线wx。
则上升沿波形曲线减去上升沿呼吸曲线,得到曲线Hs,曲线Hs=Ss-ys。曲线Hs上点的振幅,为相同位置的上升沿波形曲线点与上升沿呼吸曲线点的振幅之差。
然后,对曲线Hs进行高低通滤波,根据人体正常心跳周期,设置心跳波形的频率范围为0.8—1.6Hz,则通过高低通滤波得到的波形即为0.8—1.6Hz范围内的上升沿心跳曲线。
同理,下降沿波形曲线减去下降沿呼吸曲线,得到曲线Hx,曲线Hx=Sx-wx。下降沿心跳曲线上点的振幅,为相同位置的下降沿波形曲线点与下降沿呼吸曲线点的振幅之差。
然后,对曲线Hx进行高低通滤波,根据人体正常心跳周期,设置心跳波形的频率范围为0.8—1.6Hz,则通过高低通滤波得到的波形即为0.8—1.6Hz范围内的上升沿心跳曲线。
由此,去除子单元232得到上升沿心跳曲线和下降沿心跳曲线,也即得到心跳波形。
本实施例中,呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线,波形定位子单元231得到原始波形的上升沿波形曲线和下降沿波形曲线;然后,去除子单元232在上升沿波形曲线上去除上升沿呼吸曲线,并进行滤波,得到上升沿心跳曲线;在下降沿波形曲线上去除下降沿呼吸曲线,并进行滤波,得到下降沿心跳曲线,由此,得到心跳波形。本实施例通过波形趋势,定位上升沿波形和下降沿波形,根据上升沿波形和下降沿波形分别拟合得到上升沿呼吸曲线和下降沿呼吸曲线,然后对原始波形的上升沿波形和下降沿波形的分别进行去干扰处理,避免呼吸信号的波峰与心跳信号的波峰一致而导致的去干扰产生的误差,实现了心跳信号的准确获取。
进一步的,参照图10,本发明心跳信号处理装置第四实施例提供一种心跳信号处理装置,基于上述图9所示的实施例,所述计算模块300包括:
计算单元310,用于分别根据所述上升沿心跳曲线和所述下降沿心跳曲线,计算得到对应的上升沿心率和下降沿心率。
在获取上升沿心跳曲线和下降沿心跳曲线后,计算单元310分别计算心率值。
具体的,作为一种实施方式,首先,计算单元310获取上升沿心跳曲线和下降沿心跳曲线的波峰、波谷位置信息。
然后,计算单元310根据上升沿心跳曲线相邻两个波峰或波谷采样点位置的距离,得到两个波峰或波谷的采样间隔时间,根据两个波峰或波谷的间隔时间,计算得到采样时间节点的心率值。由此,得到上升沿心跳曲线对应的上升沿心率值HR1。
同理,计算单元310根据下降沿心跳曲线相邻两个波峰或波谷采样点位置的距离,得到两个波峰或波谷的采样间隔时间,根据两个波峰或波谷的间隔时间,计算得到采样时间节点的心率值。由此,得到下降沿心跳曲线对应的下降沿心率值HR2。
异常处理单元320,用于根据预设的处理规则,去除所述上升沿心率和下降沿心率中的异常数据。
在得到上升沿心率值HR1和下降沿心率值HR2后,异常处理单元320根据预设的规则对得到心率值进行处理,去除异常数据。
具体的,作为一种实施方式,预设的处理规则包括取舍、范围限制、平滑处理。
首先,异常处理单元320舍去开始进行心跳信号采样时、一定时间范围内的初始波形,舍去结束心跳采样时、一定时间范围内的结束波形。初始波形和结束波形由于受到采样仪器的初始化等外部原因,会引起波形异常。
然后,根据预先设置的心率范围,去除过大或过小的心率值。根据人体正常心率值,设置心率范围为每分钟40-120次,根据此范围去掉由起床等原因引起的不合理的过快心跳或慢心跳。
另外,若得到的某个采样时间节点的心率值,与此时间节点的前一个时间节点得到的心率值相比,小于前一个心率值的三分之一或大于前一个心率值的三分之四,则判定得到的心率值浮动过大,去除此心率值。
另外,若一个上升沿心跳曲线或下降沿心跳曲线对应的心跳数量小于2个,则舍弃此上升沿心跳曲线或下降沿心跳曲线对应的心率值。
需要说明的是,预设的处理规则还可以包括其他异常数据处理方式,可根据实际需要灵活设置。
由此,异常处理单元320得到去除异常数据后的上升沿心率数据和下降沿心率数据。
综合处理单元330,用于根据所述处理后的上升沿心率和下降沿心率,得到心率值。
在得到去除异常数据后的上升沿心率数据和下降沿心率数据后,综合处理单元330综合处理上升沿心率和下降沿心率,得到原始波形对应的心率值。
具体的,作为一种实施方式,综合处理单元330对得到的心率值平均值,得到的心率值即为原始波形采样时间内的心率值。
或者,作为另一种实施方式,综合处理单元330将得到的上升沿心率和下降沿心率按对应的采样时间节点排序,并按时间顺序划分为数个采样区。然后,分别求各采样区的心率平均值,再综合处理得到原始波形采样时间内的心率值。
本发明实施例采用段处理,采集心跳信号时,将得到的预设长度的原始波形作为一段原始波形进行处理,计算得到心率此段原始波形对应的心率值。其中,预设的长度可根据实际需要灵活设置。
例如,每采集70个点后,将这70个点的波形作为一段原始波形进行处理,得到对应的心跳波形。参照图6所示的心跳波形图,计算模块300得到多段心跳波形。
由于开始心跳信号采样和结束心跳信号采样时,心跳信号误差较大,因此计算模块300得到多段心跳波形后,去除开始采样和结束采样时预设时间范围内的心跳波形。
由此,计算模块300计算得到多个心率值。
然后,计算模块300对得到的多个心率值进行综合处理,去除异常数据,计算多段心跳波形的平均心率值,得到最终的心率值。
在本实施例中,计算单元310分别根据上升沿心跳曲线和下降沿心跳曲线,计算得到对应的上升沿心率和下降沿心率;然后,异常处理单元320根据预设的处理规则,去除上升沿心率和下降沿心率中的异常数据;然后,综合处理单元330根据处理后的上升沿心率和下降沿心率,得到原始波形对应的心率值。本实施例根据上升沿心跳曲线和下降沿心跳曲线分别计算得到心率之后,在对得到的心率值进行综合处理得到最终的心率值,实现了对心跳信号的处理,得到心率值。
参照图11,本发明心跳信号处理系统提供一种心跳信号处理系统,所述心跳信号处理系统包括压电传感器A、模数转换器B和心率监测CPU C,其中:
所述压电传感器A用于,采集压电模拟信号。
本实施例通过压电传感器A获取用户呼吸和心跳的身体振动信号,例如,将压电传感器置A于床垫中或床单下,获取人体睡觉时,呼吸引起的胸腔收括振动信号和心跳引起的有规律的身体震动信号。
例如,可以将压电传感器A设置为1.2米长、1-2毫米厚、5厘米宽,放置于人体平躺时靠近胸腔的位置,持续采集压电信号。
压电传感器A采集到的信号为模拟信号。
然后,压电传感器A将得到的模拟信号发送给模数转换器B。
所述模数转换器B用于,将所述压电模拟信号转换为数字信号,得到心跳信号原始波形。
在本实施例中,模数转换器B可以采用A/D转换器(Analog to DigitalConverter,模数转换器)。
A/D转换器收到压电传感器A发送的模拟信号后,将模拟线号转换为数字信号,持续得到多段心跳信号原始波形。
然后,A/D转换器将得到的多段心跳信号原始波形,按照时间顺序,依次发送给心率监测CPU C。
所述心率监测CPU C,包括获取模块100、波形处理模块200和计算模块300。
心率监测CPU C分别对得到的多段心跳信号原始波形进行处理。
具体的,作为一种事实方式,首先,获取模块100接收模数转换器B发送的心跳信号原始波形,并通过高低通滤波提取出频率范围为0.1-0.5Hz的呼吸波形。
然后,获取模块100进行差分阈值,得到呼吸波形的波峰和波谷位置。然后,根据呼吸波形的波峰和波谷位置查找原始波形对应位置区间范围内的波峰和波谷。其中,预设的区间范围可根据实际需要灵活设置。
由此,获取模块100得到原始波形的波峰和波谷位置信息。
然后,波形处理模块200根据原始波形的波峰和波谷位置信息,判断波形趋势。然后,根据波峰和波谷的位置,将相邻的波峰和波谷作为起始点或结束点,并根据波形趋势,分别拟合得到上升沿呼吸曲线和下降沿呼吸曲线。
然后,波形处理模块200得到原始波形的上升沿波形曲线和下降沿波形曲线。原始波形的上升沿波形曲线与上升沿呼吸曲线位置对应,下降沿波形曲线与下降沿呼吸曲线位置对应。
然后,波形处理模块200将上升沿波形曲线减去对应位置的上升沿呼吸曲线,并对得到的波形进行高低通滤波,得到上升沿心跳曲线;将下降沿波形曲线减去对应位置的下降沿呼吸曲线,并对得到的波形进行高低通滤波,得到下降沿心跳曲线。
然后,计算模块300根据得到的上升沿心跳曲线和下降沿呼吸曲线,分别计算得到对应的上升沿心率和下降沿心率。
然后,计算模块300对得到的上升沿心率和下降沿心率分别去除异常数据,得到人体心率合理范围内的上升沿心率数据和下降沿心率数据。
然后,计算模块300对上升沿心率数据和下降沿心率数据取平均值,得到当前原始波形的心率值。
由于本实施例包括多段心跳信号原始波形,在分别对每段原始波形进行处理,得到对应的心率值后,计算模块300舍去异常数据,例如心跳数小于2个的心跳波形对应的心率值,进行平滑处理,得到最终当前检测时间范围内的心率值。
本实施例中,心跳信号处理系统包括压电传感器A、模数转换器B和心率监测CPUC,其中压电传感器A采集压电模拟信号并发送给模数转换器B;模数转换器B将压电模拟信号转换为数字信号,得到心跳信号原始波形,并发送给心率监测CPU C;心率监测CPU C接收模数转换器B发送的心跳信号原始波形并进行处理,去除呼吸信号的干扰,计算得到心率值。本实施例通过压电传感器采集信号,心率监测CPU去除心跳信号的严重干扰,根据得到的心跳波形计算并优化处理,得到了最终的心率值,准确性高。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (5)

1.一种心跳信号处理方法,其特征在于,所述心跳信号处理方法包括以下步骤:
获取心跳信号原始波形,得到所述原始波形的波峰和波谷位置信息;
根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形;
根据所述心跳波形,计算得到心率值;
所述根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形的步骤包括:
根据所述原始波形的波峰和波谷位置信息,得到波形趋势;
根据所述原始波形的波峰和波谷位置信息、所述波形趋势,拟合呼吸波形得到呼吸曲线;
根据所述呼吸曲线,去除所述呼吸曲线对所述原始波形的干扰,并进行滤波,得到心跳波形;
所述呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线,所述心跳波形包括上升沿心跳曲线和下降沿心跳曲线,所述根据所述呼吸曲线,去除所述呼吸曲线对所述原始波形的干扰,并进行滤波,得到心跳波形的步骤包括:
根据所述原始波形的波峰和波谷位置信息、所述波形趋势,得到原始波形的上升沿波形曲线和下降沿波形曲线;
在所述上升沿波形曲线上去除所述上升沿呼吸曲线,并进行滤波,得到上升沿心跳曲线;
在所述下降沿波形曲线上去除所述下降沿呼吸曲线,并进行滤波,得到下降沿心跳曲线。
2.如权利要求1所述的心跳信号处理方法,其特征在于,所述根据所述心跳波形,计算得到心率值的步骤包括:
分别根据所述上升沿心跳曲线和所述下降沿心跳曲线,计算得到对应的上升沿心率和下降沿心率;
根据预设的处理规则,去除所述上升沿心率和下降沿心率中的异常数据;
根据所述处理后的上升沿心率和下降沿心率,得到心率值。
3.一种心跳信号处理装置,其特征在于,所述心跳信号处理装置包括:
获取模块,用于获取心跳信号原始波形,得到所述原始波形的波峰和波谷位置信息;
波形处理模块,用于根据所述原始波形的波峰和波谷位置信息拟合得到呼吸曲线,并去除所述呼吸曲线对所述原始波形的干扰,得到心跳波形;
计算模块,用于根据所述心跳波形,计算得到心率值;
所述波形处理模块包括:
波形趋势单元,用于根据所述原始波形的波峰和波谷位置信息,得到波形趋势;
拟合单元,用于根据所述原始波形的波峰和波谷位置信息、所述波形趋势,拟合呼吸波形得到呼吸曲线;
去干扰单元,用于根据所述呼吸曲线,去除所述呼吸曲线对所述原始波形的干扰,并进行滤波,得到心跳波形;
所述呼吸曲线包括上升沿呼吸曲线和下降沿呼吸曲线,所述心跳波形包括上升沿心跳曲线和下降沿心跳曲线,所述去干扰单元包括:
波形定位子单元,用于根据所述原始波形的波峰和波谷位置信息、所述波形趋势,得到原始波形的上升沿波形曲线和下降沿波形曲线;
去除子单元,用于在所述上升沿波形曲线上去除所述上升沿呼吸曲线,并进行滤波,得到上升沿心跳曲线;在所述下降沿波形曲线上去除所述下降沿呼吸曲线,并进行滤波,得到下降沿心跳曲线。
4.如权利要求3所述的心跳信号处理装置,其特征在于,所述计算模块包括:
计算单元,用于分别根据所述上升沿心跳曲线和所述下降沿心跳曲线,计算得到对应的上升沿心率和下降沿心率;
异常处理单元,用于根据预设的处理规则,去除所述上升沿心率和下降沿心率中的异常数据;
综合处理单元,用于根据所述处理后的上升沿心率和下降沿心率,得到心率值。
5.一种心跳信号处理系统,其特征在于,所述心跳信号处理系统包括压电传感器、模数转换器和心率监测CPU,其中:
所述压电传感器用于,采集压电模拟信号;
所述模数转换器用于,将所述压电模拟信号转换为数字信号,得到心跳信号原始波形;
所述心率监测CPU,包括如权利要求3-4任一项所述的心跳信号处理装置。
CN201610005670.1A 2016-01-05 2016-01-05 心跳信号处理方法、装置和系统 Active CN105662345B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201610005670.1A CN105662345B (zh) 2016-01-05 2016-01-05 心跳信号处理方法、装置和系统
PCT/CN2016/101249 WO2017118127A1 (zh) 2016-01-05 2016-09-30 一种心跳信号处理方法、装置和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610005670.1A CN105662345B (zh) 2016-01-05 2016-01-05 心跳信号处理方法、装置和系统

Publications (2)

Publication Number Publication Date
CN105662345A CN105662345A (zh) 2016-06-15
CN105662345B true CN105662345B (zh) 2018-11-16

Family

ID=56299020

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610005670.1A Active CN105662345B (zh) 2016-01-05 2016-01-05 心跳信号处理方法、装置和系统

Country Status (2)

Country Link
CN (1) CN105662345B (zh)
WO (1) WO2017118127A1 (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105662345B (zh) * 2016-01-05 2018-11-16 深圳和而泰智能控制股份有限公司 心跳信号处理方法、装置和系统
CN108697352B (zh) * 2017-06-29 2021-04-20 深圳和而泰智能控制股份有限公司 生理信息测量方法及生理信息监测装置、设备
CN107981841A (zh) * 2017-10-27 2018-05-04 深圳和而泰智能控制股份有限公司 一种信号处理方法、装置、设备及介质
CN109745026A (zh) * 2017-11-07 2019-05-14 深圳欧德蒙科技有限公司 一种心率测量方法和系统
CN110115574A (zh) * 2018-02-07 2019-08-13 普天信息技术有限公司 心率监测的方法和装置
CN109259750B (zh) * 2018-11-12 2023-04-11 浙江清华柔性电子技术研究院 心率计算方法、装置、计算机设备和存储介质
CN111358435B (zh) * 2020-03-13 2023-02-28 珠海向量科技有限公司 一种提高深度神经网络精度的数据增强方法
CN112464794A (zh) * 2020-11-25 2021-03-09 易方达基金管理有限公司 基于图像的波动趋势识别方法、装置、计算机设备和介质
CN114027822B (zh) * 2021-04-19 2022-11-25 北京超思电子技术有限责任公司 一种基于ppg信号的呼吸率测量方法及装置
CN116610221B (zh) * 2023-07-19 2024-03-12 深圳市巨朋电子科技有限公司 一种可监测用户生理指标的智能键盘及计算机系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4781201A (en) * 1984-12-27 1988-11-01 American Home Products Corporation (Del.) Cardiovascular artifact filter
JP2001252251A (ja) * 2000-03-10 2001-09-18 Yokohama Rubber Co Ltd:The 心負担評価方法および心負担評価装置
WO2004098409A1 (ja) * 2003-05-07 2004-11-18 Seijiro Tomita 心拍や呼吸等の生体信号の抽出法及び装置
CN104605829A (zh) * 2015-02-03 2015-05-13 南京理工大学 一种非接触式生命体征监护系统的心跳信号优化算法
CN105105739A (zh) * 2015-04-01 2015-12-02 上海兆观信息科技有限公司 短距离无线心率及心率变异性检测方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3521654B2 (ja) * 1996-11-07 2004-04-19 セイコーエプソン株式会社 脈拍計
KR100462182B1 (ko) * 2002-04-15 2004-12-16 삼성전자주식회사 Ppg 기반의 심박 검출 장치 및 방법
CN102440768B (zh) * 2010-10-13 2013-08-28 兰州理工大学 脉搏波形特征点提取方法
CN102393870B (zh) * 2011-06-15 2014-05-14 西安电子科技大学 脉搏波数据基线不平的修正方法
CN102626307A (zh) * 2012-01-16 2012-08-08 兰州理工大学 动态脉搏信号实时检测系统及检测方法
CN102885616A (zh) * 2012-07-17 2013-01-23 桂林电子科技大学 一种去除脉搏波信号基线漂移的方法
CN103020472B (zh) * 2012-12-27 2015-12-09 中国科学院深圳先进技术研究院 基于约束估计的生理信号质量评估方法和系统
CN103027692B (zh) * 2012-12-27 2014-11-26 天津大学 一种基于不确定度的动态光谱数据处理方法
CN103169456A (zh) * 2013-03-29 2013-06-26 深圳职业技术学院 一种脉搏波信号的处理方法及处理系统
CN104027095A (zh) * 2014-06-25 2014-09-10 哈尔滨工业大学 一种脉搏数据的预处理方法
CN104605830A (zh) * 2015-02-03 2015-05-13 南京理工大学 基于非接触式生命体征监护系统的运动趋势项消除算法
CN105662345B (zh) * 2016-01-05 2018-11-16 深圳和而泰智能控制股份有限公司 心跳信号处理方法、装置和系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4781201A (en) * 1984-12-27 1988-11-01 American Home Products Corporation (Del.) Cardiovascular artifact filter
JP2001252251A (ja) * 2000-03-10 2001-09-18 Yokohama Rubber Co Ltd:The 心負担評価方法および心負担評価装置
WO2004098409A1 (ja) * 2003-05-07 2004-11-18 Seijiro Tomita 心拍や呼吸等の生体信号の抽出法及び装置
CN104605829A (zh) * 2015-02-03 2015-05-13 南京理工大学 一种非接触式生命体征监护系统的心跳信号优化算法
CN105105739A (zh) * 2015-04-01 2015-12-02 上海兆观信息科技有限公司 短距离无线心率及心率变异性检测方法

Also Published As

Publication number Publication date
WO2017118127A1 (zh) 2017-07-13
CN105662345A (zh) 2016-06-15

Similar Documents

Publication Publication Date Title
CN105662345B (zh) 心跳信号处理方法、装置和系统
KR101656611B1 (ko) 무구속적으로 측정한 생체신호를 이용하여 산소탈포화지수를 획득하는 방법
US20200260962A1 (en) System and methods for acquisition and analysis of health data
CA2809764C (en) Systems and methods for respiratory event detection
JP5153770B2 (ja) いびき検知及び確認のためのシステム及び方法
US20210219925A1 (en) Apparatus and method for detection of physiological events
US6626842B2 (en) Heart-sound analyzing apparatus
JP5090155B2 (ja) 呼吸速度、心拍速度、及び無呼吸に関する非侵襲性モニタリング
CN105769122B (zh) 一种睡眠呼吸暂停低通气指数的测量设备
WO2019096175A1 (zh) 一种生命体征信号分析处理方法和生命体征监测设备
US20080082018A1 (en) Systems and methods for respiratory event detection
CN109414204A (zh) 用于确定针对对象的呼吸信息的方法和装置
CN107920766A (zh) 睡眠现象监测
KR20070048201A (ko) 비접촉 마이크로폰을 사용한 호흡 패턴 판별 장치 및 방법
CN108888271B (zh) 一种生理参数测量系统及配备该测量系统的智能座椅
CN103153183A (zh) 用于诊断阻塞性睡眠呼吸暂停的装置和方法
CN112155560B (zh) 基于实时心冲击信号的呼吸暂停检测方法及系统
CN106308801A (zh) 一种利用智能手机检测人体呼吸频率的方法
JP2628690B2 (ja) 呼吸数モニタ装置
CN110448280A (zh) 基于压电薄膜传感器实现的瑜伽冥想质量检测方法及系统
CN110115583A (zh) 呼吸监测的方法和装置
WO2017042350A1 (en) Method and system for monitoring ventilatory parameter
CN115397318A (zh) 用于确定呼吸努力的系统和方法
AU2011203044C1 (en) Systems and methods for respiratory event detection
JP2018130513A (ja) 呼吸成分抽出方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20191030

Address after: 518000, building 10, building ten, building D, Shenzhen Institute of Aerospace Science and technology, 6 hi tech Southern District, Nanshan District, Shenzhen, Guangdong 1003, China

Patentee after: Shenzhen heertai Home Furnishing Online Network Technology Co. Ltd.

Address before: 518057, building 10, building ten, building D, Shenzhen Institute of Aerospace Science and technology, 6 hi tech Southern District, Nanshan District, Shenzhen, Guangdong 1010-1011, China

Patentee before: ShenZhen H&T Intelligent Control CO., LTD.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200409

Address after: 1706, Fangda building, No. 011, Keji South 12th Road, high tech Zone, Yuehai street, Nanshan District, Shenzhen City, Guangdong Province

Patentee after: Shenzhen shuliantianxia Intelligent Technology Co., Ltd

Address before: 518000, building 10, building ten, building D, Shenzhen Institute of Aerospace Science and technology, 6 hi tech Southern District, Nanshan District, Shenzhen, Guangdong 1003, China

Patentee before: SHENZHEN H & T HOME ONLINE NETWORK TECHNOLOGY Co.,Ltd.