CN107569227A - 一种运动状态下心率的处理方法和监测装置 - Google Patents
一种运动状态下心率的处理方法和监测装置 Download PDFInfo
- Publication number
- CN107569227A CN107569227A CN201711013463.1A CN201711013463A CN107569227A CN 107569227 A CN107569227 A CN 107569227A CN 201711013463 A CN201711013463 A CN 201711013463A CN 107569227 A CN107569227 A CN 107569227A
- Authority
- CN
- China
- Prior art keywords
- ripples
- buffer area
- heart rate
- point
- time
- 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.)
- Pending
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明涉及一种运动状态下心率的处理方法,包括以下步骤:步骤一、采集运动状态下四秒缓存区域内的心电数字信号;步骤二、采用巴特沃斯带通滤波器,滤除了心电数字信号中的基线漂移,工频干扰,肌电噪声,只保留QRS波;步骤三、采用类方差的滤波继续加强的R波;步骤四、确定区间阀值,步骤五、定位R波:步骤六、找出缓冲区域内存在的R波的个数count;步骤七、计算心率;本发明抗干扰能力强,在运动伪差干扰较大的情况下,能够准确的识别出R波,进而准确的计算出心率,有效的解决了运动伪差所造成心率计算的不准确的问题。此外,本发明还提供了运动状态下心率监测的装置。
Description
技术领域
本发明属于运动负荷的评价方法和装置,具体涉及一种运动状态下心率的处理方法和检测装置。
背景技术
心率是反应运动强度与能量消耗的重要指标之一,可以通过心率表收集运动心率和时间,并推算能量消耗水平,因此,监测运动过程中的心率变化状况是反映运动负荷有效的手段之一,其通常的操作是先将体表不同部位的电信号检测出来,再用放大器加以放大,并用记录器描记下来,就可得到心电图形,然后根据心电图形所反映的R波计算出心率值;因此,心率计算的问题实际上也就是正确识别心电信号中R波的问题。但由于运动过程中,由于运动伪差的干扰,R波存在着大量的误检漏检,因此大多数方法对于心率的提取存在不准确的问题。
发明内容
本发明的目的在于克服现有技术的缺陷,提供一种运动状态下心率的处理方法和监测装置,采用该方法和装置,能够准确的识别运动状态下心电信号中R波,并根据R波计算出心率值。
为了实现上述目的,本发明采取的技术方案如下:
技术方案一:
一种运动状态下心率的处理方法,包括以下步骤:
步骤一、设定采样频率和缓冲区域的采样时间,所述缓冲区域的采样时间为0.5秒的整数倍,且大于1.5秒,然后采集运动状态下缓冲区域内的心电数字信号;
步骤二、采用巴特沃斯带通滤波器,滤除了心电数字信号中的基线漂移,工频干扰,肌电噪声以及P波和T波,只保留QRS波;
步骤三、采用类方差滤波计算方法继续加强R波的幅值,得到加强后的心电信号;
步骤四、确定区间阀值:首先将缓存区域的采样时间分成N'个区间,N'≥3,查找每个区间内的采样点幅值最大值V'max,将N'个区间的V'max进行从小到大的排序,去掉排序末位的值,然后将其余V'max做一个平均,得到的平均值即为区间阈值Vth;
步骤五、定位R波:在缓冲区域内从前到后依次对每一个采样点的幅值V'i进行判断,找出高度大于等于区间阈值Vth的第一个采样点,则以第一个采样点为起点,查找向后0.12秒内的最大值点,即峰值点,该峰值点即为缓冲区域内的第一个R波,记录缓冲区域内的第一个R波出现的时间点Tr1;
步骤六、记录下Tr1后,并向后跳过一个不应期,然后继续按照步骤五所述步骤查找高度大于等于阀值的下一个采样点,并找出下一个R波出现的时间Tr,直到缓冲区域的尾端为止,最后计算缓冲区域内存在的R波的个数Rcount,以及缓冲区域内最后一个R波出现的时间点TrR;
步骤七、确定了所有R波的位置后,计算缓冲区域的实时心率,
其中:HR为实时心率值,TrR为缓冲区间内最后一个R波出现的时间点,Tr1为缓冲区间内第一个R波出现的时间点,Rcount为缓冲区域内R波的总个数。
进一步的,所述类方差滤波的计算方法为:采集到的第i个采样点的心电信号幅值为Vi,缓冲区域内总采样点个数N,N≥3,对每个采样点前后取n个采样点进行平均幅值的计算,n≥1;
2n+1个采样点的平均幅值为:Vmean=(Vi-n+Vi-n+1+Vi-n+2+...+Vi+...+Vi+n)/(2n+1);
加强后的心电信号幅值为:V'i=(Vi-n-Vmean)2+(Vi-n+1-Vmean)2+...+(Vi-Vmean)2+...+(Vi+n-Vmean)2;
进一步的,所述不应期的时间为60秒/200,即0.3秒。
技术方案二:
一种运动状态下心率的监测装置,包括用于对体表心电信号采集的心电模拟信号采集与放大模块,通过AD转换模块与所述心电模拟采集与放大模块相连的MCU,以及分别与MCU连接的蓝牙模块和振动马达;所述模拟信号采集与放大模块用于采集体表心电信号后并将所采集的体表心电信号进行放大;所述AD转换模块用于将模拟信号采集与放大模块放大后的心电信号转换为数字信号并传入到MCU,所述MCU用于接收AD转换模块传入的数字信号,并对接收的数字信号进行滤波和分析处理后得到实时心率值;所述蓝牙模块用于与手机端进行通信,传输计算所得实时心率值;振动马达,用于当MCU判断心率发生异常时,进行振动报警。
进一步的,所述实时心率值的获得方法为:
步骤一、设定采样频率和缓冲区域的采样时间,所述缓冲区域的采样时间为0.5秒的整数倍且大于1.5秒,然后采集运动状态下缓冲区域内的心电数字信号;
步骤二、采用巴特沃斯带通滤波器,滤除了心电数字信号中的基线漂移,工频干扰,肌电噪声以及P波和T波,只保留QRS波;
步骤三、采用类方差滤波计算方法继续加强R波的幅值;得到加强后的心电信号;所述类方差滤波的计算方法为:所述类方差滤波的计算方法为:采集到的第i个采样点的心电信号幅值为Vi,缓冲区域内总采样点个数N,N≥3,对每个采样点前后取n个采样点进行平均幅值的计算,n≥1;
2n+1个采样点的平均幅值为:Vmean=(Vi-n+Vi-n+1+Vi-n+2+...+Vi+...+Vi+n)/(2n+1);
加强后的心电信号幅值为:V'i=(Vi-n-Vmean)2+(Vi-n+1-Vmean)2+...+(Vi-Vmean)2+...+(Vi+n-Vmean)2;
步骤四、确定区间阀值:首先将缓存区域的采样时间分成N'个区间,N'≥3,查找每个区间内的采样点幅值最大值V'max,将N'个区间的V'max进行从小到大的排序,去掉排序末位的值,然后将其余V'max做一个平均,得到的平均值即为区间阈值Vth;
步骤五、定位R波:在缓冲区域内从前到后依次对每一个采样点的幅值V'i进行判断,找出高度大于等于区间阈值Vth的第一个采样点,则以第一个采样点为起点,查找向后0.12秒内的最大值点,即峰值点,该峰值点即为缓冲区域内的第一个R波,记录缓冲区域内的第一个R波出现的时间点Tr1;
步骤六、记录下Tr后,并向后跳过0.3秒的不应期,然后继续按照步骤五所述步骤查找高度大于等于阀值的下一个采样点,并找出下一个R波出现的时间Tr,直到缓冲区域的尾端为止,最后计算缓冲区域内存在的R波的个数Rcount,以及缓冲区域内最后一个R波出现的时间点TrR;
步骤七、确定了所有R波的位置后,计算缓冲区域的实时心率,
其中:HR为实时心率值,TrR为缓冲区间内最后一个R波出现的时间点,Tr1为缓冲区间内第一个R波出现的时间点,Rcount为缓冲区域内R波的总个数。
与现有技术相比,本发明所取得的有益效果如下:
1.本发明所需计算量小,可以直接在监测装置中进行计算。
2.本发明抗干扰能力强,在运动伪差干扰较大的情况下,能够准确的识别出R波,进而准确的计算出心率。经过大量运动数据测试,包括上下楼梯,快跑,蹦跳等情况在内,R波抓取准确率均在90%以上,有效的解决了运动伪差所造成心率计算的不准确的问题。
附图说明
图1为原始心电信号图;
图2为类方差滤波加强后的心电信号图;
图3为本发明监测装置的系统结构图;
具体实施方式
以下结合实施例对本发明进行进一步详细的叙述。
如图1~2所示的本发明一种运动状态下心率的处理方法的一个实施例,包括如下步骤:
步骤一、设定采样频率和缓冲区域的采样时间,所述缓冲区域的采样时间为0.5秒的整数倍,且大于1.5秒,然后采集运动状态下缓冲区域内的心电数字信号;(本实施例设定的缓冲区域的采用时间为4s,)
步骤二、采用巴特沃斯带通滤波器,滤除了心电数字信号中的基线漂移,工频干扰,肌电噪声以及P波和T波,只保留QRS波;
步骤三、采用类方差滤波计算方法继续加强R波的幅值,得到加强后的心电信号;所述类方差滤波的计算方法为:采集到的第i个采样点的心电信号幅值为Vi,缓冲区域内总采样点个数N,N≥3,对每个采样点前后取n个采样点进行平均幅值的计算,n≥1;
2n+1个采样点的平均幅值为:Vmean=(Vi-n+Vi-n+1+Vi-n+2+...+Vi+...+Vi+n)/(2n+1);
加强后的心电信号幅值为:V'i=(Vi-n-Vmean)2+(Vi-n+1-Vmean)2+...+(Vi-Vmean)2+...+(Vi+n-Vmean)2;
步骤四、确定区间阀值:首先将缓存区域的采样时间分成N'个区间(本实施例中N'=8,每一个区间0.5秒,),查找每个区间内的采样点幅值最大值V'max,将N'个区间的V'max进行从小到大的排序,去掉排序末位的值,然后将其余V'max做一个平均,得到的平均值即为区间阈值Vth;
步骤五、定位R波:在缓冲区域内从前到后依次对每一个采样点的幅值V'i进行判断,找出高度大于等于区间阈值Vth的第一个采样点,则以第一个采样点为起点,查找向后0.12秒内的最大值点,即峰值点,该峰值点即为缓冲区域内的第一个R波,记录缓冲区域内的第一个R波出现的时间点Tr1;其中0.12秒为人在正常状态下的R波宽度;
步骤六、记录下Tr后,并向后跳过一个不应期,然后继续按照步骤五所述步骤查找高度大于等于阀值的下一个采样点,并找出下一个R波出现的时间Tr,直到缓冲区域的尾端为止,最后计算缓冲区域内存在的R波的个数Rcount(本实施例中找到了5个的R波),以及缓冲区域内最后一个R波出现的时间点TrR;所述不应期指的是在生物对某一刺激发生反应后,在一定时间内,即使再给予刺激,也不发生反应,一般称此期间为不应期;因此,在寻找R波的过程中,找到一个R波后,需要跳过一段对刺激不发生反应的时间,即不应期;所述不应期的之间为60秒/200,即0.3秒;其中,200为60秒内心率的最大支持次数为200次;
步骤七、确定了所有R波的位置后,计算缓冲区域的实时心率,
其中:HR为实时心率值,TrR为缓冲区间内最后一个R波出现的时间点,Tr1为缓冲区间内第一个R波出现的时间点,Rcount为缓冲区域内R波的总个数。
如图3所示的一种运动状态下心率的监测装置的一个实施例,用于对体表心电信号采集的心电模拟信号采集与放大模块,通过AD转换模块与所述心电模拟采集与放大模块相连的MCU,以及分别与MCU连接的蓝牙模块和振动马达;所述模拟信号采集与放大模块用于采集体表心电信号后并将所采集的体表心电信号进行放大;所述AD转换模块用于将模拟信号采集与放大模块放大后的心电信号转换为数字信号并传入到MCU,所述MCU用于接收AD转换模块传入的数字信号,并对接收的数字信号进行滤波和分析处理后得到实时心率值;所述蓝牙模块用于与手机端进行通信,传输计算所得实时心率值;振动马达,用于当MCU判断心率发生异常时,进行振动报警。
进一步的,所述滤波采用巴特沃斯带通滤波器,所述巴特沃斯带通滤波器集成于所述MCU上。
进一步的,所述实时心率值的获得方法为:
步骤一、设定采样频率和缓冲区域的采样时间,所述缓冲区域的采样时间为0.5秒的整数倍且大于1.5秒,然后采集运动状态下缓冲区域内的心电数字信号;
步骤二、采用巴特沃斯带通滤波器,滤除了心电数字信号中的基线漂移,工频干扰,肌电噪声以及P波和T波,只保留QRS波;
步骤三、采用类方差滤波计算方法继续加强R波的幅值;得到加强后的心电信号;所述类方差滤波的计算方法为:所述类方差滤波的计算方法为:采集到的第i个采样点的心电信号幅值为Vi,缓冲区域内总采样点个数N,N≥3,对每个采样点前后取n个采样点进行平均幅值的计算,n≥1;
2n+1个采样点的平均幅值为:Vmean=(Vi-n+Vi-n+1+Vi-n+2+...+Vi+...+Vi+n)/(2n+1);
加强后的心电信号幅值为:V'i=(Vi-n-Vmean)2+(Vi-n+1-Vmean)2+...+(Vi-Vmean)2+...+(Vi+n-Vmean)2;
步骤四、确定区间阀值:首先将缓存区域的采样时间分成N'个区间,N'≥3,查找每个区间内的采样点幅值最大值V'max,将N'个区间的V'max进行从小到大的排序,去掉排序末位的值,然后将其余V'max做一个平均,得到的平均值即为区间阈值Vth;
步骤五、定位R波:在缓冲区域内从前到后依次对每一个采样点的幅值V'i进行判断,找出高度大于等于区间阈值Vth的第一个采样点,则以第一个采样点为起点,查找向后0.12秒内的最大值点,即峰值点,该峰值点即为缓冲区域内的第一个R波,记录缓冲区域内的第一个R波出现的时间点Tr1;
步骤六、记录下Tr后,并向后跳过0.3秒的不应期,然后继续按照步骤五所述步骤查找高度大于等于阀值的下一个采样点,并找出下一个R波出现的时间Tr,直到缓冲区域的尾端为止,最后计算缓冲区域内存在的R波的个数Rcount,以及缓冲区域内最后一个R波出现的时间点TrR;
步骤七、确定了所有R波的位置后,计算缓冲区域的实时心率,
其中:HR为实时心率值,TrR为缓冲区间内最后一个R波出现的时间点,Tr1为缓冲区间内第一个R波出现的时间点,Rcount为缓冲区域内R波的总个数。
以上所述实施方式仅为本发明的优选实施例,而并非本发明可行实施的穷举。对于本领域一般技术人员而言,在不背离本发明原理和精神的前提下对其所作出的任何显而易见的改动,都应当被认为包含在本发明的权利要求保护范围之内。
Claims (5)
1.一种运动状态下心率的处理方法,其特征在于,包括以下步骤:
步骤一、设定采样频率和缓冲区域的采样时间,所述缓冲区域的采样时间为0.5秒的整数倍,且大于1.5秒,然后采集运动状态下缓冲区域内的心电数字信号;
步骤二、采用巴特沃斯带通滤波器,滤除了心电数字信号中的基线漂移,工频干扰,肌电噪声以及P波和T波,只保留QRS波;
步骤三、采用类方差滤波计算方法继续加强R波的幅值,得到加强后的心电信号;
步骤四、确定区间阀值:首先将缓存区域的采样时间分成N'个区间,N'≥3,查找每个区间内的采样点幅值最大值V′max,将N'个区间的V′max进行从小到大的排序,去掉排序末位的值,然后将其余V′max做一个平均,得到的平均值即为区间阈值Vth;
步骤五、定位R波:在缓冲区域内从前到后依次对每一个采样点的幅值Vi'进行判断,找出高度大于等于区间阈值Vth的第一个采样点,则以第一个采样点为起点,查找向后0.12秒内的最大值点,即峰值点,该峰值点即为缓冲区域内的第一个R波,记录缓冲区域内的第一个R波出现的时间点Tr1;
步骤六、记录下Tr1后,并向后跳过一个不应期,然后继续按照步骤五所述步骤查找高度大于等于阀值的下一个采样点,并找出下一个R波出现的时间Tr,直到缓冲区域的尾端为止,最后计算缓冲区域内存在的R波的个数Rcount,以及缓冲区域内最后一个R波出现的时间点TrR;
步骤七、确定了所有R波的位置后,计算缓冲区域的实时心率,
其中:HR为实时心率值,TrR为缓冲区间内最后一个R波出现的时间点,Tr1为缓冲区间内第一个R波出现的时间点,Rcount为缓冲区域内R波的总个数。
2.根据权利要求1所述的一种运动状态下心率的处理方法,其特征在于,
所述类方差滤波的计算方法为:采集到的第i个采样点的心电信号幅值为Vi,缓冲区域内总采样点个数N,N≥3,对每个采样点前后取n个采样点进行平均幅值的计算,n≥1;2n+1个采样点的平均幅值为:Vmean=(Vi-n+Vi-n+1+Vi-n+2+...+Vi+...+Vi+n)/(2n+1);
加强后的心电信号幅值为:
Vi'=(Vi-n-Vmean)2+(Vi-n+1-Vmean)2+...+(Vi-Vmean)2+...+(Vi+n-Vmean)2。
3.根据权利要求1所述的一种运动状态下心率的处理方法,其特征在于,所述不应期的时间为60秒/200,即0.3秒。
4.一种运动状态下心率的监测装置,其特征在于,包括用于对体表心电信号采集的心电模拟信号采集与放大模块,通过AD转换模块与所述心电模拟采集与放大模块相连的MCU,以及分别与MCU连接的蓝牙模块和振动马达;所述模拟信号采集与放大模块用于采集体表心电信号后并将所采集的体表心电信号进行放大;所述AD转换模块用于将模拟信号采集与放大模块放大后的心电信号转换为数字信号并传入到MCU,所述MCU用于接收AD转换模块传入的数字信号,并对接收的数字信号进行滤波和分析处理后得到实时心率值;所述蓝牙模块用于与手机端进行通信,传输计算所得实时心率值;振动马达,用于当MCU判断心率发生异常时,进行振动报警。
5.根据权利要求4所述的一种运动状态下心率的监测装置,其特征在于,所述实时心率值的获得方法为:
步骤一、设定采样频率和缓冲区域的采样时间,所述缓冲区域的采样时间为0.5秒的整数倍且大于1.5秒,然后采集运动状态下缓冲区域内的心电数字信号;
步骤二、采用巴特沃斯带通滤波器,滤除了心电数字信号中的基线漂移,工频干扰,肌电噪声以及P波和T波,只保留QRS波;
步骤三、采用类方差滤波计算方法继续加强R波的幅值;得到加强后的心电信号;所述类方差滤波的计算方法为:所述类方差滤波的计算方法为:采集到的第i个采样点的心电信号幅值为Vi,缓冲区域内总采样点个数N,N≥3,对每个采样点前后取n个采样点进行平均幅值的计算,n≥1;
2n+1个采样点的平均幅值为:Vmean=(Vi-n+Vi-n+1+Vi-n+2+...+Vi+...+Vi+n)/(2n+1);
加强后的心电信号幅值为:
Vi'=(Vi-n-Vmean)2+(Vi-n+1-Vmean)2+...+(Vi-Vmean)2+...+(Vi+n-Vmean)2;
步骤四、确定区间阀值:首先将缓存区域的采样时间分成N'个区间,N'≥3,查找每个区间内的采样点幅值最大值V′max,将N'个区间的V′max进行从小到大的排序,去掉排序末位的值,然后将其余V′max做一个平均,得到的平均值即为区间阈值Vth;
步骤五、定位R波:在缓冲区域内从前到后依次对每一个采样点的幅值Vi'进行判断,找出高度大于等于区间阈值Vth的第一个采样点,则以第一个采样点为起点,查找向后0.12秒内的最大值点,即峰值点,该峰值点即为缓冲区域内的第一个R波,记录缓冲区域内的第一个R波出现的时间点Tr1;
步骤六、记录下Tr后,并向后跳过0.3秒的不应期,然后继续按照步骤五所述步骤查找高度大于等于阀值的下一个采样点,并找出下一个R波出现的时间Tr,直到缓冲区域的尾端为止,最后计算缓冲区域内存在的R波的个数Rcount,以及缓冲区域内最后一个R波出现的时间点TrR;
步骤七、确定了所有R波的位置后,计算缓冲区域的实时心率,
其中:HR为实时心率值,TrR为缓冲区间内最后一个R波出现的时间点,Tr1为缓冲区间内第一个R波出现的时间点,Rcount为缓冲区域内R波的总个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711013463.1A CN107569227A (zh) | 2017-10-26 | 2017-10-26 | 一种运动状态下心率的处理方法和监测装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711013463.1A CN107569227A (zh) | 2017-10-26 | 2017-10-26 | 一种运动状态下心率的处理方法和监测装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107569227A true CN107569227A (zh) | 2018-01-12 |
Family
ID=61037057
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711013463.1A Pending CN107569227A (zh) | 2017-10-26 | 2017-10-26 | 一种运动状态下心率的处理方法和监测装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107569227A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108294737A (zh) * | 2018-01-26 | 2018-07-20 | 深圳还是威健康科技有限公司 | 心率测量方法、装置及智能穿戴设备 |
CN108652611A (zh) * | 2018-05-17 | 2018-10-16 | 广东工业大学 | 一种心率的计算方法、装置、设备及计算机可读存储介质 |
CN108836316A (zh) * | 2018-05-02 | 2018-11-20 | 湖北工业大学 | 一种基于bp神经网络的心电信号r波提取方法 |
CN112006679A (zh) * | 2020-08-28 | 2020-12-01 | 山东省人工智能研究院 | 一种基于窗口方差变换的穿戴式心电信号r波检测方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040215092A1 (en) * | 2002-09-20 | 2004-10-28 | Fischell David R. | Means and method for the detection of cardiac events |
CN203153717U (zh) * | 2013-04-09 | 2013-08-28 | 南方医科大学 | 一种具有跌倒报警功能的移动健康监护装置 |
CN103549950A (zh) * | 2013-11-19 | 2014-02-05 | 上海理工大学 | 移动心电监测差分阈值检测改进算法 |
CN106037720A (zh) * | 2015-12-04 | 2016-10-26 | 贵州大学 | 混合连续信息分析技术在医学上的应用方法 |
-
2017
- 2017-10-26 CN CN201711013463.1A patent/CN107569227A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040215092A1 (en) * | 2002-09-20 | 2004-10-28 | Fischell David R. | Means and method for the detection of cardiac events |
CN203153717U (zh) * | 2013-04-09 | 2013-08-28 | 南方医科大学 | 一种具有跌倒报警功能的移动健康监护装置 |
CN103549950A (zh) * | 2013-11-19 | 2014-02-05 | 上海理工大学 | 移动心电监测差分阈值检测改进算法 |
CN106037720A (zh) * | 2015-12-04 | 2016-10-26 | 贵州大学 | 混合连续信息分析技术在医学上的应用方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108294737A (zh) * | 2018-01-26 | 2018-07-20 | 深圳还是威健康科技有限公司 | 心率测量方法、装置及智能穿戴设备 |
CN108294737B (zh) * | 2018-01-26 | 2020-11-13 | 深圳市元征科技股份有限公司 | 心率测量方法、装置及智能穿戴设备 |
CN108836316A (zh) * | 2018-05-02 | 2018-11-20 | 湖北工业大学 | 一种基于bp神经网络的心电信号r波提取方法 |
CN108652611A (zh) * | 2018-05-17 | 2018-10-16 | 广东工业大学 | 一种心率的计算方法、装置、设备及计算机可读存储介质 |
CN108652611B (zh) * | 2018-05-17 | 2021-03-16 | 广东工业大学 | 一种心率的计算方法、装置、设备及计算机可读存储介质 |
CN112006679A (zh) * | 2020-08-28 | 2020-12-01 | 山东省人工智能研究院 | 一种基于窗口方差变换的穿戴式心电信号r波检测方法 |
CN112006679B (zh) * | 2020-08-28 | 2021-08-03 | 山东省人工智能研究院 | 一种基于窗口方差变换的穿戴式心电信号r波检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107569227A (zh) | 一种运动状态下心率的处理方法和监测装置 | |
CN106344005B (zh) | 一种可移动心电图监测系统 | |
CN104586384B (zh) | 一种心电信号实时心率检测方法及系统 | |
CN102835954B (zh) | 一种心拍波形模板生成方法及模块 | |
CN103654770B (zh) | 移动心电信号qrs波实时波检测方法及装置 | |
CN103099615B (zh) | 一种消除运动心电信号干扰的方法和装置 | |
CN101732050B (zh) | 一种基于光电容积波的呼吸率监测方法 | |
CN103860162B (zh) | 一种起搏信号检测方法、系统和心电检测设备 | |
US8152733B2 (en) | Method for improving recognition rate of respiratory wave | |
CN110664390A (zh) | 基于腕带式ppg和深度学习的心率监测系统及方法 | |
CN108601543B (zh) | 一种ecg信号处理方法及装置 | |
CN103431856A (zh) | 在多导联同步心电信号中选择心电导联的方法及装置 | |
CN103690156A (zh) | 一种心率获取方法及心电信号的处理方法 | |
CN102319063A (zh) | 一种提高胎心率数据加速识别准确性的装置和方法 | |
CN105997043B (zh) | 一种基于腕式可穿戴设备的脉率提取方法 | |
CN112754444A (zh) | 一种基于雷达的非接触式猪只呼吸检测方法 | |
CN109009087B (zh) | 一种心电信号r波的快速检测方法 | |
CN103083011B (zh) | 利用胸阻抗二阶差分图辅助实时定位心电r波峰的方法 | |
CN102613964A (zh) | 获取生理信号周期的方法及系统 | |
CN101176662A (zh) | 心率测量方法 | |
CN106419884A (zh) | 一种基于小波分析的心率计算方法与系统 | |
CN106667474A (zh) | 一种适用于智能穿戴设备的心率测量方法及装置 | |
CN111166322A (zh) | 心电信号特征波提取方法 | |
CN115770028A (zh) | 血压检测方法、系统、装置及存储介质 | |
CN106175755B (zh) | 用于睡眠状态分析的睡眠状态检测器 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180112 |
|
RJ01 | Rejection of invention patent application after publication |