CN109998522A - 一种基于穿戴式动态心电的运动负荷检测装置及方法 - Google Patents
一种基于穿戴式动态心电的运动负荷检测装置及方法 Download PDFInfo
- Publication number
- CN109998522A CN109998522A CN201910146847.3A CN201910146847A CN109998522A CN 109998522 A CN109998522 A CN 109998522A CN 201910146847 A CN201910146847 A CN 201910146847A CN 109998522 A CN109998522 A CN 109998522A
- Authority
- CN
- China
- Prior art keywords
- signal
- value
- vector
- wave
- feature vector
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1118—Determining activity level
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6801—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient specially adapted to be attached to or worn on the body surface
- A61B5/6802—Sensor mounted on worn items
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
- A61B5/7207—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise induced by motion artifacts
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Fuzzy Systems (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Cardiology (AREA)
- Power Engineering (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供了一种可穿戴式的基于动态心电的运动负荷检测装置及方法,装置包括心电信号采集模块、信号调理模块、无线通讯模块、数据存储模块、微处理器,通过采集用户动态心电信号,对信号进行滤波、放大以及AD转换,并传输心电数字信号信息至用户客户端,随后完成包括信号滤波、R波特征检测、特征向量计算、向量归一化、运动负荷检测在内的信号分析过程,最终生成检测报告并发生至监护中心。该装置穿戴方便、舒适,可实现对用户动态心电信号长时间采集与实时分析处理,准确检测用户运动负荷并给出相应建议的功能。
Description
技术领域
本发明属于人体运动过程中心脏健康实时监测领域,具体涉及一种基于穿戴式动态心电 的运动负荷检测装置及相应的检测方法。
背景技术
运动在合理强度下会减少疾病的发生有助于身体健康,但如果运动过度反而会使身体出 现运动损伤的情况,特别是对患有心脏疾病的人群。有数据表明,82%非外伤性猝死与患有 心脏疾病人群运动不当相关,因此科学健身至关重要。
1973年瑞典生理学家Gunnar Bog根据心理学原理在自我感觉和实际运动负荷量评定等 级上设计了主观体力感觉等级表,这是运动对象对自己身体运动负荷的一种自感疲劳分级。 然而这种评估方法结果不准确且对于普通人群来说不适用。随着运动生理学的发展,研究人 员通过测量运动对象的耗氧量、血乳酸浓度、睾酮浓度等生理指标来综合评定运动负荷,这 种方法能较精确的监测到个体对于不同运动负荷的适应情况。然而这些指标的测量需要复杂 的检测设备并且整个过程也十分繁琐,需要强有力的资源支持,目前只适用于专业运动员或 者科研院校研究中,现阶段难以普及。
心电信号(Electrocardiosignal,ECG)能够综合反映人体体内生理的各种变化,并能够 较为准确地描述身体机能对于运动刺激的反应,是对个体综合代谢情况观察的窗口。1842年 意大利生理学家Matteucci观察到心脏搏动时存在电流活动以来,ECG的研究得到了迅猛发 展。早在1990年美国运动医学会就开始提出采用心率和自感疲劳分级主客观信息相结合的方 法确定运动负荷,但那时对于心率的实时测量较为困难,一般采用价格昂贵的心率表。而现 如今随着物联网和可穿戴式心电设备技术的迅速发展,使得实时、不间断远程监控个体心电 信号成为了可能,并广泛应用于及普通民众日常生理活动监测当中。而近几年来,可穿戴式 心电设备技术得到迅猛发展,市面上出现了很多心电监测穿戴式产品,虽然目前利用心率给 出运动建议的方法有很多,但基本上都是基于最佳心率算法演化而来,并没有对动态心电做 进一步的深入分析从多维角度对个体运动情况进行评判,形成完整的运动负荷监测方法,因 此不具有太大的应用参考价值。
发明内容
为解决上述问题,本发明提供了一种可穿戴式的基于动态心电的运动负荷检测装置及方 法。本装置采用I导联心电采集方式采集用户运动过程的动态心电信号,将原始信号滤波、 放大、AD转换后经无线通讯实时上传至手机客户端,最后对信号进行检测、处理以及分析 后得出运动负荷检测结果反馈给用户。
为了达到上述目的,本发明提供如下技术方案:
一种基于穿戴式动态心电的运动负荷检测装置,包括:
心电信号采集模块:与信号调理模块连接,用于采集人体心电信号;
信号调理模块:与微处理器连接,用于进行信号的滤波、放大、AD转换;
无线通讯模块:与微处理器连接,用于通过无线方式根据预先设定的通讯频率将心电数 字信号信息传输至云平台及用户手机客户端进行信号分析;
数据存储模块:与微处理器连接,用于存储用户一段时间内的心电信号数据;
微处理器:与信号调理模块、无线通讯模块和数据存储模块连接,是装置的控制中心, 用于发出命令使各模块正常工作,完成功能的实现;
所述信号分析包括如下步骤:
①信号滤波:除去动态心电信号中存在的工频噪声、基线漂移以及运动干扰;
②R波特征检测:对滤波后的心电信号进行R波特征检测;
③特征向量计算:根据步骤②中检测出的R波峰值点位置得出RR间期时间序列计算出 心率(Heart rate,HR),并进行心率变异性(Heart rate variability,HRV)分析,包括时、频域 及非线性分析得出相关特征向量;
④向量归一化:将步骤③中的特征向量构建成特征矩阵,并对每个特征向量归一化;
⑤运动负荷检测及报告生成:将标准化后的特征矩阵作为输入参数实时输入基于遗传算 法优化的神经网络模型,输出运动负荷检测结果并打印结果报告。
进一步的,所述步骤①信号滤波具体包括如下子步骤:
将采集到的心电信号利用Mallat算法进行小波变换得到不同尺度下信号αj(n),其中 j=1,2,…N,利用小波的多分辨率分析特性对带有工频噪声、基线漂移和运动干扰的心电信号 按照采样频率进行分解,将不同频率组成的含噪心电信号分解成不同频带的子信号;再针对 干扰及噪声所在的不同频带的子信号设定阈值处理进行去噪,对含有心电信号特征波的频段 给予保留,其他的频段阈值化处理进行去噪;之后再利用离散小波逆变换重构出有用心电信 号。
进一步的,所述步骤②R波特征检测具体包括如下子步骤:
根据多尺度分解后的信号确定R波所对应的正-负模极大值对,然后分别确定模极大值阈 值Tmax以及模极小值阈值Tmin;并使用阈值筛选出正极大值大于Tmax,负极大值小于Tmin且 两个极值点之间的距离小于阈值的正-负模极大值对;检测它们的过零点,并使用(2j-1)/2 修正时移;最后使用不应期和回溯的策略对检测出的R波进行修正,得出各R波位置坐标为 (γ1,γ2,…,γk),k为R波点个数。
进一步的,所述步骤③特征向量计算具体包括如下子步骤:
设RR间期时间序列为x()=[RR1,RR2,…,RRn],其中RRj=γj+1-γj(1≤j≤k-1),RR间期均值则心率由步骤①中得出的R波各点位置求出RR间期序列,采用统计学离散趋势法对RR间期时间序列进行HRV时域研究,其中SDNN 表示为RMSSD为相邻RR间期差的均方根值,计算式为 SDSD为相邻RR间期标准差,由计算得到,其中RRSi=RRi+1-RRi;PNN50为相邻RR间期之差大于50ms的个数除以总的RR间期个数的 百分比;
利用基于AR模型的功率谱估计法完成HRV的频域参数求解:AR模型中假定观测的数 据RR时间序列x(n)用方差为σ2的零均值白噪声序列w(n)激励一个全极点的线性系统产生, 模型的阶数为M,系数为aM,k,用差分方程表示x(n)的AR模型为 过程{x(n)}为M阶自回归过程,功率谱密度为使用AR自 回归模型采取Marple算法估算参数aM,k后估算功率谱Pxx(ω)并进行谱分解,计算总功率TP, 低频功率LF以及高频功率HF;利用公式对LF、HF进行归 一化得到参数LFnorm和HFnorm;同时计算LF与HF的比值LF/HF;HRV非线性动力学研 究采用Poincare散点图法,以RR间期时间序列第i个值为作为横坐标,以第(i+1)个值为纵坐 标绘制出Poincare散点图;对Poincare散点图分析得到两个定量分析参数,分别是向量角度 指数SD1和向量长度指数SD2,计算方法分别如下:
进一步的,所述步骤④向量归一化具体包括如下子步骤:
首先构建特征向量J:J=[HR,SDNN,RMSSD,PNN50,LF,HF,LFnorm,HForm,SD1,SD2],对 步骤③得到的每个特征向量采取最大-最小标准化法进行归一化处理,设J为上述特征向量构 成的特征矩阵,minJk为第k个特征向量最小值,maxJk为第k个特征向量最大值;则初始值x 映射到区间[0,1]的值
进一步的,所述步骤⑤中的遗传算法优化的神经网络模型以按误差逆向传播(error Back Propagation,简称BP)算法训练的多层前馈神经网络为基础,训练过程如下:给定训练集 即输入示例由d个属性描述,输出l维实值向量; 对训练例(xk,yk),设神经网络输出则网络在(xk,yk)上的均方误差为 BP算法对于误差Ek及给定学习率η,基于梯度下降策略对层间的权值ω' 和阈值θ'不断更新使训练集D上的累计误差达到最小化的目标。
进一步的,所述心电信号采集模块包括织物干性电极以及信号预处理电路,织物干性电 极包括两个分别置于左右锁骨下方的电极,以及置于右腹部的参考电极,这些电极构成差分 回路用于采集人体心电信号;信号预处理电路包括RC滤波器和右腿驱动电路,用于对原始 心电信号进行简易滤波及去除共模干扰。
进一步的,所述信号调理模块包括低噪声可编程增益放大器、高分辨率的模数转换器、 振荡器、参考电压接口和SPI数据接口,具有高共模抑制比,用于滤除共模信号噪声,并具 有灵活的断电、待机模式。
一种基于穿戴式动态心电的运动负荷检测方法,包括如下步骤:
采集用户的动态心电信号,进行信号的滤波、放大、AD转换,最后进行信号分析;
所述信号分析的过程包括如下步骤:
①信号滤波:除去动态心电信号中存在的工频噪声、基线漂移以及运动干扰;
②R波特征检测:对滤波后的心电信号进行R波特征检测;
③特征向量计算:根据步骤②中检测出的R波峰值点位置得出RR间期时间序列计算出 心率(Heart rate,HR),并进行心率变异性(Heart rate variability,HRV)分析,包括时、频域 及非线性分析得出相关特征向量;
④向量归一化:将步骤③中的特征向量构建成特征矩阵,并对每个特征向量归一化;
⑤运动负荷检测及报告生成:将标准化后的特征矩阵作为输入参数实时输入基于遗传算 法优化的神经网络模型,输出运动负荷检测结果并打印结果报告。
进一步的,所述步骤①信号滤波具体包括如下子步骤:
将采集到的心电信号利用Mallat算法进行小波变换得到不同尺度下信号αj(n),其中 j=1,2,…N,利用小波的多分辨率分析特性对带有工频噪声、基线漂移和运动干扰的心电信号 按照采样频率进行分解,将不同频率组成的含噪心电信号分解成不同频带的子信号;再针对 干扰及噪声所在的不同频带的子信号设定阈值处理进行去噪,对含有心电信号特征波的频段 给予保留,其他的频段阈值化处理进行去噪;之后再利用离散小波逆变换重构出有用心电信 号;
所述步骤②R波特征检测具体包括如下子步骤:
根据多尺度分解后的信号确定R波所对应的正-负模极大值对,然后分别确定模极大值阈 值Tmax以及模极小值阈值Tmin;并使用阈值筛选出正极大值大于Tmax,负极大值小于Tmin且 两个极值点之间的距离小于阈值的正-负模极大值对;检测它们的过零点,并使用(2j-1)/2 修正时移;最后使用不应期和回溯的策略对检测出的R波进行修正,得出各R波位置坐标为 (γ1,γ2,…,γk),k为R波点个数;
所述步骤③特征向量计算具体包括如下子步骤:
设RR间期时间序列为x(n)=[RR1,RR2,…,RRn],其中RRj=γj+1-γj(1≤j≤k-1),RR间期均值则心率由步骤①中得出的R波各点位置求出RR间期序列,采用统计学离散趋势法对RR间期时间序列进行HRV时域研究,其中SDNN 表示为RMSSD为相邻RR间期差的均方根值,计算式为 SDSD为相邻RR间期标准差,由计算得到,其中RRSi=RRi+1-RRi;PNN50为相邻RR间期之差大于50ms的个数除以总的RR间期个数的 百分比;
利用基于AR模型的功率谱估计法完成HRV的频域参数求解:AR模型中假定观测的数 据RR时间序列x(n)用方差为σ2的零均值白噪声序列w(n)激励一个全极点的线性系统产生, 模型的阶数为M,系数为aM,k,用差分方程表示x(n)的AR模型为 过程{x(n)}为M阶自回归过程,功率谱密度为使用AR自 回归模型采取Marple算法估算参数aM,k后估算功率谱Pxx(ω)并进行谱分解,计算总功率TP, 低频功率LF以及高频功率HF;利用公式对LF、HF进行归 一化得到参数LFnorm和HFnorm;同时计算LF与HF的比值LF/HF;HRV非线性动力学研 究采用Poincare散点图法,以RR间期时间序列第i个值为作为横坐标,以第(i+1)个值为纵坐 标绘制出Poincare散点图;对Poincare散点图分析得到两个定量分析参数,分别是向量角度 指数SD1和向量长度指数SD2,计算方法分别如下:
所述步骤④向量归一化具体包括如下子步骤:
首先构建特征向量J:J=[HR,SDNN,RMSSD,PNN50,LF,HF,LFnorm,HForm,SD1,SD2],对 步骤③得到的每个特征向量采取最大-最小标准化法进行归一化处理,设J为上述特征向量构 成的特征矩阵,minJk为第k个特征向量最小值,maxJk为第k个特征向量最大值;则初始值x 映射到区间[0,1]的值
所述步骤⑤中的遗传算法优化的神经网络模型以按误差逆向传播(error BackPropagation, 简称BP)算法训练的多层前馈神经网络为基础,训练过程如下:给定训练集 即输入示例由d个属性描述,输出l维实值向量; 对训练例(xk,yk),设神经网络输出则网络在(xk,yk)上的均方误差为 BP算法对于误差Ek及给定学习率η,基于梯度下降策略对层间的权值ω' 和阈值θ'不断更新使训练集D上的累计误差达到最小化的目标。
与现有技术相比,本发明具有如下优点和有益效果:
1.本发明提供的装置能够完成运动对象心电数据的长程实时采集工作,并针对动态心电信号 设计算法进行滤波去噪工作,为特征检测、提取有效信息进行运动负荷评估提供有利的条件; 分析方法能够无创的完成运动对象的运动负荷评估工作,为用户提供个性化运动指导提供基 础,降低用户运动过疲劳出现突发性心脏疾病或其他运动损伤的风险,让用户进行科学、健 康性的运动。
2.装置使用简单方便,穿戴较为舒适,适用于日常生活场景,可实现对用户动态心电信号长 时间采集与实时分析处理,准确检测用户运动负荷并给出相应建议的功能。
附图说明
图1为本发明一种基于穿戴式动态心电的运动负荷检测装置的结构示意图。
图2为本发明运动负荷检测流程示意图。
图3为本发明对低强度动态心电信号的检测结果图。
图4为本发明对中强度动态心电信号的检测结果图。
图5为本发明对高强度动态心电信号的检测结果图。
附图标记说明:
1、心电信号采集模块,2、信号调理模块,3、微处理器,4、数据存储模块,5、蓝牙模块, 6、手机客户端,7、电源,8、工作指示灯,9、工作按钮。
具体实施方式
以下将结合具体实施例对本发明提供的技术方案进行详细说明,应理解下述具体实施方 式仅用于说明本发明而不用于限制本发明的范围。
如图1所示,本发明提供的一种基于穿戴式动态心电的运动负荷检测装置主要包括心电 信号采集模块1、信号调理模块2、微处理器3、数据存储模块4和无线传输模块5。心电信 号采集模块1与信号调理模块2连接,信号调理模块2、数据存储模块4和无线传输模块5 均与微处理器3连接,同时与微处理器3连接的还有电源7、工作指示灯8和工作按钮9。心电信号采集模块1用于采集心电信号,信号调理模块2将采集到的信号进行滤波、放大、AD转换等工作;随后经微处理器3控制通过无线传输模块5通过Wifi和蓝牙传输方案根据预先设定的通讯频率将心电数字信号传输至云平台和用户手机客户端6完成信号分析。数据存储 模块4用于存储用户一段时间内的心电信号数据;电源7用于给装置各个模块供电;工作指 示灯8用于指示装置是否处于采集信号状态。工作按钮9用于控制装置的工作状态。
具体的说,心电信号采集模块1包括织物干性电极以及信号预处理电路,以双极单导联 方式获取差分心电信号。其中织物干性电极包括置于人体左、右锁骨下方的RA、LA两个织 物干性电极与置于右腹部作为参考电极的RL电极,这些电极构成差分回路用于采集人体心 电信号。信号预处理电路由RC滤波器和右腿驱动电路组成,完成对原始心电信号的简易滤 波及去除共模干扰。
信号调理模块内含低噪声可编程增益放大器以及高分辨率的模数转换器,具有高共模抑 制比,可以很好的滤除共模信号噪声。同时内置振荡器和参考电压,SPI数据接口,具有灵 活的断电、待机模式。
图2给出了本发明运动负荷检测流程示意图。心电信号采集模块中的织物干性电极采集 用户的动态心电信号,然后经过信号调理模块进行信号的滤波、放大、AD转换,最后通过 无线通讯模块上传至客户端完成信号分析,信号分析包括信号滤波、R波特征检测、特征向 量计算、向量归一化、运动负荷检测及报告生成。具体包括如下过程:
①信号滤波:对动态心电信号信号进行滤波,除去动态心电信号中主要存在的工频噪声、 基线漂移以及运动干扰。
将采集到的心电信号按5分钟进行分段,对每段信号利用Mallat算法进行8尺度离散小 波分解,具体的,是将采集到的信号利用Mallat算法进行小波变换得到不同尺度下信号αj(n), 其中j=1,2,…N。由于噪声和信号在多尺度分解下的小波系数具有不同空间分布的特点,利用 小波的多分辨率分析特性对带有工频噪声、基线漂移和运动干扰的心电信号按照采样频率进 行分解,将不同频率组成的含噪心电信号分解成不同频带的子信号。再针对干扰及噪声所在 的不同频带的子信号设定阈值处理进行去噪,对含有心电信号特征波的频段给予保留,其他 的频段阈值化处理进行去噪。之后再利用离散小波逆变换重构出有用ECG信号。
②R波特征检测:对滤波后的心电信号进行R波特征检测。
根据多尺度分解后的信号确定R波所对应的正-负模极大值对,然后分别确定模极大值阈 值Tmax以及模极小值阈值Tmin。并使用阈值筛选出正极大值大于Tmax,负极大值小于Tmin且 两个极值点之间的距离小于200ms的正-负模极大值对。检测它们的过零点,并使用(2j-1)/2 修正时移,其中j为小波分解的阶数。最后使用不应期和回溯的策略对检测出的R波进行修 正,得出各R波位置坐标为(γ1,γ2,…,γk),k为R波点个数。
③特征向量计算:根据步骤②中检测出的R波峰值点位置得出RR间期时间序列计算出 心率(Heart rate,HR),并进行心率变异性(Heart rate variability,HRV)分析,包括时、频域 及非线性分析得出相关特征向量;
具体的,由上述步骤①中得出的R波各点位置求出RR间期序列x(n)= [RR1,RR2,…,RRk-1],其中RRj=γj+1-γj(1≤j≤k-1)。然后计算用于检测运动负荷相关 的特征向量,包括HR、HRV时、频域及非线性分析特征向量。RR间期均值则心采用统计学离散趋势法对RR间期时间序列进行HRV时域研究,计 算RR间期标准差SDNN,该值被认为与HRV的总变化情况有关,表示为 其中为RR间期序列的平均值,N为RR间期序列的长度。
RMSSD为相邻RR间期差的均方根值,计算式为SDSD为相邻 RR间期标准差,由计算得到,其中RRSi=RRi+1-RRi。PNN50为相 邻RR间期之差大于50ms的个数除以总的RR间期个数的百分比,RMSSD与PNN50均被认为与副交感神经(Parasympathetic,PPS)的活性有关。
利用基于AR模型的功率谱分析法完成HRV的频域参数求解。AR模型中假定观测的数 据RR间序列x(n)用方差为σ2的零均值白噪声序列w(n)激励一个全极点的线性系统产生,模 型的阶数为M,系数为aM,k,用差分方程表示x(n)的AR模型为 过程{x(n)}为M阶自回归过程,功率谱密度为其中PM为功 率谱能量,w为对应的角频率。使用AR自回归模型采取Marple算法估算参数aM,k后估算功 率谱Pxx(ω)并进行谱分解,计算总功率TP,低频功率LF(0.04-0.15Hz)以及高频功率HF (0.15-0.40Hz)。利用公式对LF、HF进行归一化得到参数 LFnorm和Hfnorm。同时计算LF与HF的比值LF/HF。HRV非线性动力学研究采用Poincare 散点图法,以上述步骤中得到的RR间期时间序列第i个值为作为横坐标,以第(i+1)个值为纵 坐标绘制出Poincare散点图。对Poincare散点图分析得到两个定量分析参数,分别是向量角 度指数SD1和向量长度指数SD2。画出Poincare散点图后,将离散点的分布近似为椭圆分布, 椭圆的中心位置由RR间期均值确定,椭圆半长轴和半短轴的长度即分别为SD1和SD2值的 大小。其中SD1是由迷走神经传出影响窦房结引起的即刻或短期RR变异性的标准偏差的指 标,SD2是表示全过程中心率的长期或缓慢变化标准偏差的指标。计算方法分别如下:
④向量归一化:将步骤③中的特征向量构建成特征矩阵,并对每个特征向量归一化;
首先构建特征向量J:J=[HR,SDNN,RMSSD,PNN50,LF,HF,LFnorm,HForm,SD1,SD2],对 步骤③得到的每个特征向量采取最大-最小标准化法进行归一化处理,设J为上述特征向量构 成的特征矩阵,minJk为第k个特征向量最小值,maxJk为第k个特征向量最大值。则初始值x 映射到区间[0,1]的值
⑤运动负荷检测及报告生成:将标准化后的特征矩阵作为输入参数实时输入经样本集训 练后的基于遗传算法优化的神经网络模型,输出运动负荷检测结果并打印结果报告。图3-图 5中显示的1分钟的信号均源于为本发明装置分别采集的三段5分钟信号,相应的检测结果 分别为低强度、中强度和高强度运动负荷。最后在用户的手机客户端中实时显示检测结果, 并将结果发送至监护中心结合大数据反馈给用户相应的运动建议。
上述基于遗传算法优化的神经网络模型以按误差逆向传播(error BackPropagation,简称 BP)算法训练的多层前馈神经网络为基础。训练过程如下:给定训练集 即输入示例由d个属性描述,输出l维实值向量。 对训练例(xk,yk),设神经网络输出则网络在(xk,yk)上的均方误差为 BP算法对于误差Ek及给定学习率η,基于梯度下降策略对层间的权值ω' 和阈值θ'不断更新使训练集D上的累计误差达到最小化的目标。然而BP算法 容易陷入局部最优,而遗传算法是一种通过模拟自然选择和遗传学机理进化过程搜索最 优解的方法,具有全局搜索寻优性、鲁棒性强和稳健性高的特点,因此本发明中使用遗传算 法优化初始权值和阈值使BP神经网络更好地逼近全局最小。
本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段,还包括由以上 技术特征任意组合所组成的技术方案。应当指出,对于本技术领域的普通技术人员来说,在 不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的 保护范围。
Claims (10)
1.一种基于穿戴式动态心电的运动负荷检测装置,其特征在于,包括:
心电信号采集模块:与信号调理模块连接,用于采集人体心电信号;
信号调理模块:与微处理器连接,用于进行信号的滤波、放大、AD转换;
无线通讯模块:与微处理器连接,用于通过无线方式根据预先设定的通讯频率将心电数字信号信息传输至云平台及用户手机客户端进行信号分析;
数据存储模块:与微处理器连接,用于存储用户一段时间内的心电信号数据;
微处理器:与信号调理模块、无线通讯模块和数据存储模块连接,是装置的控制中心,用于发出命令使各模块正常工作,完成功能的实现;
所述信号分析包括如下步骤:
①信号滤波:除去动态心电信号中存在的工频噪声、基线漂移以及运动干扰;
②R波特征检测:对滤波后的心电信号进行R波特征检测;
③特征向量计算:根据步骤②中检测出的R波峰值点位置得出RR间期时间序列计算出心率(Heart rate,HR),并进行心率变异性(Heart rate variability,HRV)分析,包括时、频域及非线性分析得出相关特征向量;
④向量归一化:将步骤③中的特征向量构建成特征矩阵,并对每个特征向量归一化;
⑤运动负荷检测及报告生成:将标准化后的特征矩阵作为输入参数实时输入基于遗传算法优化的神经网络模型,输出运动负荷检测结果并打印结果报告。
2.根据权利要求1所述的基于穿戴式动态心电的运动负荷检测装置,其特征在于,所述步骤①信号滤波具体包括如下子步骤:
将采集到的心电信号利用Mallat算法进行小波变换得到不同尺度下信号αj(n),其中j=1,2,...N,利用小波的多分辨率分析特性对带有工频噪声、基线漂移和运动干扰的心电信号按照采样频率进行分解,将不同频率组成的含噪心电信号分解成不同频带的子信号;再针对干扰及噪声所在的不同频带的子信号设定阈值处理进行去噪,对含有心电信号特征波的频段给予保留,其他的频段阈值化处理进行去噪;之后再利用离散小波逆变换重构出有用心电信号。
3.根据权利要求1所述的基于穿戴式动态心电的运动负荷检测装置,其特征在于,所述步骤②R波特征检测具体包括如下子步骤:
根据多尺度分解后的信号确定R波所对应的正-负模极大值对,然后分别确定模极大值阈值Tmax以及模极小值阈值Tmin;并使用阈值筛选出正极大值大于Tmax,负极大值小于Tmin且两个极值点之间的距离小于阈值的正-负模极大值对;检测它们的过零点,并使用(2j-1)/2修正时移,其中j为小波分解的阶数;最后使用不应期和回溯的策略对检测出的R波进行修正,得出各R波位置坐标为(γ1,γ2,...,γk),k为R波点个数。
4.根据权利要求1所述的基于穿戴式动态心电的运动负荷检测装置,其特征在于,所述步骤③特征向量计算具体包括如下子步骤:
设RR间期时间序列为x(n)=[RR1,RR2,...,RRn],其中RRj=γj+1-γj(1≤j≤k-1),RR间期均值则心率由步骤①中得出的R波各点位置求出RR间期序列,采用统计学离散趋势法对RR间期时间序列进行HRV时域研究,其中SDNN表示为其中为RR间期序列的平均值,N为RR间期序列的长度;RMSSD为相邻RR间期差的均方根值,计算式为SDSD为相邻RR间期标准差,由计算得到,其中RRSi=RRi+1-RRi;PNN50为相邻RR间期之差大于50ms的个数除以总的RR间期个数的百分比;
利用基于AR模型的功率谱估计法完成HRV的频域参数求解:AR模型中假定观测的数据RR时间序列x(n)用方差为σ2的零均值白噪声序列w(n)激励一个全极点的线性系统产生,模型的阶数为M,系数为aM,k,用差分方程表示x(n)的AR模型为 过程{x(n)}为M阶自回归过程,功率谱密度为其中PM为功率谱能量,w为对应的角频率;使用AR自回归模型采取Marple算法估算参数aM,k后估算功率谱Pxx(ω)并进行谱分解,计算总功率TP,低频功率LF以及高频功率HF;利用公式对LF、HF进行归一化得到参数LFnorm和HFnorm;同时计算LF与HF的比值LF/HF;HRV非线性动力学研究采用Poincare散点图法,以RR间期时间序列第i个值为作为横坐标,以第(i+1)个值为纵坐标绘制出Poincare散点图;对Poincare散点图分析得到两个定量分析参数,分别是向量角度指数SD1和向量长度指数SD2,计算方法分别如下:
5.根据权利要求1所述的基于穿戴式动态心电的运动负荷检测装置,其特征在于,所述步骤④向量归一化具体包括如下子步骤:
首先构建特征向量J:J=[HR,SDNN,RMSSD,PNN50,LF,HF,LFnorm,HForm,SD1,SD2],对步骤③得到的每个特征向量采取最大-最小标准化法进行归一化处理,设J为上述特征向量构成的特征矩阵,minJk为第k个特征向量最小值,maxJk为第k个特征向量最大值;则初始值x映射到区间[0,1]的值
6.根据权利要求1所述的基于穿戴式动态心电的运动负荷检测装置,其特征在于,所述步骤⑤中的遗传算法优化的神经网络模型以按误差逆向传播(error Back Propagation,简称BP)算法训练的多层前馈神经网络为基础,训练过程如下:给定训练集 即输入示例由d个属性描述,输出l维实值向量;对训练例(xk,yk),设神经网络输出则网络在(xk,yk)上的均方误差为BP算法对于误差Ek及给定学习率η,基于梯度下降策略对层间的权值ω′和阈值θ′不断更新使训练集D上的累计误差达到最小化的目标。
7.根据权利要求1所述的基于穿戴式动态心电的运动负荷检测装置,其特征在于,所述心电信号采集模块包括织物干性电极以及信号预处理电路,织物干性电极包括两个分别置于左右锁骨下方的电极,以及置于右腹部的参考电极,这些电极构成差分回路用于采集人体心电信号;信号预处理电路包括RC滤波器和右腿驱动电路,用于对原始心电信号进行简易滤波及去除共模干扰。
8.根据权利要求1所述的基于穿戴式动态心电的运动负荷检测装置,其特征在于,所述信号调理模块包括低噪声可编程增益放大器、高分辨率的模数转换器、振荡器、参考电压接口和SPI数据接口,具有高共模抑制比,用于滤除共模信号噪声,并具有灵活的断电、待机模式。
9.一种基于穿戴式动态心电的运动负荷检测方法,其特征在于,包括如下步骤:
采集用户的动态心电信号,进行信号的滤波、放大、AD转换,最后进行信号分析;
所述信号分析的过程包括如下步骤:
①信号滤波:除去动态心电信号中存在的工频噪声、基线漂移以及运动干扰;
②R波特征检测:对滤波后的心电信号进行R波特征检测;
③特征向量计算:根据步骤②中检测出的R波峰值点位置得出RR间期时间序列计算出心率(Heart rate,HR),并进行心率变异性(Heart rate variability,HRV)分析,包括时、频域及非线性分析得出相关特征向量;
④向量归一化:将步骤③中的特征向量构建成特征矩阵,并对每个特征向量归一化;
⑤运动负荷检测及报告生成:将标准化后的特征矩阵作为输入参数实时输入基于遗传算法优化的神经网络模型,输出运动负荷检测结果并打印结果报告。
10.根据权利要求9所述的基于穿戴式动态心电的运动负荷检测方法,其特征在于,所述步骤①信号滤波具体包括如下子步骤:
将采集到的心电信号利用Mallat算法进行小波变换得到不同尺度下信号αj(n),其中j=1,2,...N,利用小波的多分辨率分析特性对带有工频噪声、基线漂移和运动干扰的心电信号按照采样频率进行分解,将不同频率组成的含噪心电信号分解成不同频带的子信号;再针对干扰及噪声所在的不同频带的子信号设定阈值处理进行去噪,对含有心电信号特征波的频段给予保留,其他的频段阈值化处理进行去噪;之后再利用离散小波逆变换重构出有用心电信号;
所述步骤②R波特征检测具体包括如下子步骤:
根据多尺度分解后的信号确定R波所对应的正-负模极大值对,然后分别确定模极大值阈值Tmax以及模极小值阈值Tmin;并使用阈值筛选出正极大值大于Tmax,负极大值小于Tmin且两个极值点之间的距离小于阈值的正-负模极大值对;检测它们的过零点,并使用(2j-1)/2修正时移,其中j为小波分解的阶数;最后使用不应期和回溯的策略对检测出的R波进行修正,得出各R波位置坐标为(γ1,γ2,...,γk),k为R波点个数;
所述步骤③特征向量计算具体包括如下子步骤:
设RR间期时间序列为x(n)=[RR1,RR2,...,RRn],其中RRj=γj+1-γj(1≤j≤k-1),RR间期均值则心率由步骤①中得出的R波各点位置求出RR间期序列,采用统计学离散趋势法对RR间期时间序列进行HRV时域研究,其中SDNN表示为其中为RR间期序列的平均值,N为RR间期序列的长度;RMSSD为相邻RR间期差的均方根值,计算式为SDSD为相邻RR间期标准差,由计算得到,其中RRSi=RRi+1-RRi;PNN50为相邻RR间期之差大于50ms的个数除以总的RR间期个数的百分比;
利用基于AR模型的功率谱估计法完成HRV的频域参数求解:AR模型中假定观测的数据RR时间序列x(n)用方差为σ2的零均值白噪声序列w(n)激励一个全极点的线性系统产生,模型的阶数为M,系数为aM,k,用差分方程表示x(n)的AR模型为 过程{x(n)}为M阶自回归过程,功率谱密度为其中PM为功率谱能量,w为对应的角频率;使用AR自回归模型采取Marple算法估算参数aM,k后估算功率谱Pxx(ω)并进行谱分解,计算总功率TP,低频功率LF以及高频功率HF;利用公式对LF、HF进行归一化得到参数LFnorm和HFnorm;同时计算LF与HF的比值LF/HF;HRV非线性动力学研究采用Poincare散点图法,以RR间期时间序列第i个值为作为横坐标,以第(i+1)个值为纵坐标绘制出Poincare散点图;对Poincare散点图分析得到两个定量分析参数,分别是向量角度指数SD1和向量长度指数SD2,计算方法分别如下:
所述步骤④向量归一化具体包括如下子步骤:
首先构建特征向量J:J=[HR,SDNN,RMSSD,PNN50,LF,HF,LFnorm,HForm,SD1,SD2],对步骤③得到的每个特征向量采取最大-最小标准化法进行归一化处理,设J为上述特征向量构成的特征矩阵,minJk为第k个特征向量最小值,maxJk为第k个特征向量最大值;则初始值x映射到区间[0,1]的值
所述步骤⑤中的遗传算法优化的神经网络模型以按误差逆向传播(error BackPropagation,简称BP)算法训练的多层前馈神经网络为基础,训练过程如下:给定训练集 即输入示例由d个属性描述,输出l维实值向量;对训练例(xk,yk),设神经网络输出则网络在(xk,yk)上的均方误差为 BP算法对于误差Ek及给定学习率η,基于梯度下降策略对层间的权值ω′和阈值θ′不断更新使训练集D上的累计误差达到最小化的目标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910146847.3A CN109998522A (zh) | 2019-02-27 | 2019-02-27 | 一种基于穿戴式动态心电的运动负荷检测装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910146847.3A CN109998522A (zh) | 2019-02-27 | 2019-02-27 | 一种基于穿戴式动态心电的运动负荷检测装置及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109998522A true CN109998522A (zh) | 2019-07-12 |
Family
ID=67166138
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910146847.3A Pending CN109998522A (zh) | 2019-02-27 | 2019-02-27 | 一种基于穿戴式动态心电的运动负荷检测装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109998522A (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110807443A (zh) * | 2019-11-27 | 2020-02-18 | 南京信息工程大学 | 基于参数自适应hmm模型的体感网心电信号分段方法 |
CN110916631A (zh) * | 2019-12-13 | 2020-03-27 | 东南大学 | 基于可穿戴生理信号监测的学生课堂学习状态评测系统 |
CN110974212A (zh) * | 2019-11-07 | 2020-04-10 | 曲阜师范大学 | 一种融合心电和肌电特征的康复训练运动状态监测方法及系统 |
CN111557661A (zh) * | 2020-05-15 | 2020-08-21 | 深圳清华大学研究院 | 心电信号处理方法及装置 |
CN111938580A (zh) * | 2020-07-02 | 2020-11-17 | 北京体育大学 | 人体日常运动中摄氧量的测试方法、装置、计算机设备 |
CN112568881A (zh) * | 2020-11-13 | 2021-03-30 | 河北省药品医疗器械检验研究院 | 一种心肺功能评估系统 |
CN112914536A (zh) * | 2021-03-24 | 2021-06-08 | 平安科技(深圳)有限公司 | 运动状态的检测方法、装置、计算机设备和存储介质 |
CN113288168A (zh) * | 2021-05-21 | 2021-08-24 | 天津工业大学 | 一种智能可穿戴疲劳监测及预警系统 |
CN113288158A (zh) * | 2021-05-27 | 2021-08-24 | 河北省科学院应用数学研究所 | 基线漂移和高频噪声的去除方法、装置和设备 |
CN114159072A (zh) * | 2021-12-03 | 2022-03-11 | 信阳市中心医院 | 一种便携式心电检测装置及心电信号分析方法 |
CN114343664A (zh) * | 2022-01-12 | 2022-04-15 | 东南大学 | 一种基于心率刺激-恢复趋势的心肺耐力评估方法与装置 |
WO2022233027A1 (zh) * | 2021-05-07 | 2022-11-10 | 中国科学院深圳先进技术研究院 | 一种基于可穿戴设备的运动性疲劳检测方法和系统 |
CN117435907A (zh) * | 2023-08-14 | 2024-01-23 | 北京体育大学 | 一种基于数据分析的体育记录检测方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104490390A (zh) * | 2014-12-30 | 2015-04-08 | 天津大学 | 基于神经电生理信号联合分析的人体运动能力判别方法 |
CN106073755A (zh) * | 2016-05-27 | 2016-11-09 | 成都信汇聚源科技有限公司 | 一种微型动态心电监测设备上房颤自动识别的实现方法 |
CN108932452A (zh) * | 2017-05-22 | 2018-12-04 | 中国科学院半导体研究所 | 基于多尺度卷积神经网络的心律失常心拍分类方法 |
-
2019
- 2019-02-27 CN CN201910146847.3A patent/CN109998522A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104490390A (zh) * | 2014-12-30 | 2015-04-08 | 天津大学 | 基于神经电生理信号联合分析的人体运动能力判别方法 |
CN106073755A (zh) * | 2016-05-27 | 2016-11-09 | 成都信汇聚源科技有限公司 | 一种微型动态心电监测设备上房颤自动识别的实现方法 |
CN108932452A (zh) * | 2017-05-22 | 2018-12-04 | 中国科学院半导体研究所 | 基于多尺度卷积神经网络的心律失常心拍分类方法 |
Non-Patent Citations (11)
Title |
---|
FIRSTBEAT TECHNOLOGIES LTD: "Stress and Recovery Analysis Method Based on 24-hour Heart Rate Variability", 《HTTPS://ASSETS.FIRSTBEAT.COM/FIRSTBEAT/UPLOADS/2015/11/STRESS-AND-RECOVERY_WHITE-PAPER_20145.PDF》 * |
FIRSTBEAT TECHNOLOGIES LTD: "Stress and Recovery Analysis Method Based on 24-hour Heart Rate Variability", 《HTTPS://ASSETS.FIRSTBEAT.COM/FIRSTBEAT/UPLOADS/2015/11/STRESS-AND-RECOVERY_WHITE-PAPER_20145.PDF》, 4 November 2014 (2014-11-04), pages 5 * |
张坤等: "一种基于小波变换的心电去噪算法", 《现代生物医学进展》 * |
张坤等: "一种基于小波变换的心电去噪算法", 《现代生物医学进展》, vol. 9, no. 19, 15 October 2009 (2009-10-15), pages 1054 - 2 * |
林家瑞: "《微机式医学仪器设计》", 30 September 2004, pages: 303 * |
王小川等: "《MATLAB神经网络43个案例分析》", 31 August 2013, pages: 29 * |
石波等: "静息状态和实时运动状态心率变异性比较研究", 《中国医疗器械杂志》 * |
石波等: "静息状态和实时运动状态心率变异性比较研究", 《中国医疗器械杂志》, vol. 41, no. 3, 30 May 2017 (2017-05-30), pages 168 - 160 * |
董精通: "一种实时心电监护的穿戴式智慧衣", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
董精通: "一种实时心电监护的穿戴式智慧衣", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》, no. 12, 15 December 2017 (2017-12-15), pages 12 - 14 * |
邱雅竹等: "基于Mexican-hat小波的QRS检测新方法", 《生物医学工程学杂志》, vol. 23, no. 6, 31 December 2006 (2006-12-31), pages 1347 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110974212A (zh) * | 2019-11-07 | 2020-04-10 | 曲阜师范大学 | 一种融合心电和肌电特征的康复训练运动状态监测方法及系统 |
CN110807443B (zh) * | 2019-11-27 | 2023-04-25 | 南京信息工程大学 | 基于参数自适应hmm模型的体感网心电信号分段方法 |
CN110807443A (zh) * | 2019-11-27 | 2020-02-18 | 南京信息工程大学 | 基于参数自适应hmm模型的体感网心电信号分段方法 |
CN110916631B (zh) * | 2019-12-13 | 2022-04-22 | 东南大学 | 基于可穿戴生理信号监测的学生课堂学习状态评测系统 |
CN110916631A (zh) * | 2019-12-13 | 2020-03-27 | 东南大学 | 基于可穿戴生理信号监测的学生课堂学习状态评测系统 |
CN111557661A (zh) * | 2020-05-15 | 2020-08-21 | 深圳清华大学研究院 | 心电信号处理方法及装置 |
CN111557661B (zh) * | 2020-05-15 | 2024-03-08 | 深圳清华大学研究院 | 心电信号处理方法及装置 |
CN111938580A (zh) * | 2020-07-02 | 2020-11-17 | 北京体育大学 | 人体日常运动中摄氧量的测试方法、装置、计算机设备 |
CN112568881A (zh) * | 2020-11-13 | 2021-03-30 | 河北省药品医疗器械检验研究院 | 一种心肺功能评估系统 |
CN112914536A (zh) * | 2021-03-24 | 2021-06-08 | 平安科技(深圳)有限公司 | 运动状态的检测方法、装置、计算机设备和存储介质 |
CN112914536B (zh) * | 2021-03-24 | 2023-08-15 | 平安科技(深圳)有限公司 | 运动状态的检测方法、装置、计算机设备和存储介质 |
WO2022233027A1 (zh) * | 2021-05-07 | 2022-11-10 | 中国科学院深圳先进技术研究院 | 一种基于可穿戴设备的运动性疲劳检测方法和系统 |
CN113288168A (zh) * | 2021-05-21 | 2021-08-24 | 天津工业大学 | 一种智能可穿戴疲劳监测及预警系统 |
CN113288158A (zh) * | 2021-05-27 | 2021-08-24 | 河北省科学院应用数学研究所 | 基线漂移和高频噪声的去除方法、装置和设备 |
CN114159072A (zh) * | 2021-12-03 | 2022-03-11 | 信阳市中心医院 | 一种便携式心电检测装置及心电信号分析方法 |
CN114343664A (zh) * | 2022-01-12 | 2022-04-15 | 东南大学 | 一种基于心率刺激-恢复趋势的心肺耐力评估方法与装置 |
CN117435907A (zh) * | 2023-08-14 | 2024-01-23 | 北京体育大学 | 一种基于数据分析的体育记录检测方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109998522A (zh) | 一种基于穿戴式动态心电的运动负荷检测装置及方法 | |
CN108478209B (zh) | 心电信息动态监护方法和动态监护系统 | |
EP3692904B1 (en) | Method and device for self-learning dynamic electrocardiography analysis employing artificial intelligence | |
CN101658425B (zh) | 基于心率变异分析的注意力集中程度的检测装置及方法 | |
Belgacem et al. | ECG based human authentication using wavelets and random forests | |
CN104605841A (zh) | 可穿戴心电信号监测装置及方法 | |
Wang et al. | A wearable ECG monitor for deep learning based real-time cardiovascular disease detection | |
CN110353704B (zh) | 基于穿戴式心电监测的情绪评估方法与装置 | |
CN109288515B (zh) | 基于穿戴式心电信号中早搏信号的周期性监测方法及装置 | |
CN110236508A (zh) | 一种无创血压连续监测方法 | |
Ranganathan et al. | Estimation of heart rate signals for mental stress assessment using neuro fuzzy technique | |
Yu et al. | DDCNN: A deep learning model for AF detection from a single-lead short ECG signal | |
CN114052744A (zh) | 基于脉冲神经网络的心电信号分类方法 | |
CN114648040A (zh) | 生命体征多生理信号提取、融合分析方法 | |
Li et al. | Detection of muscle fatigue by fusion of agonist and synergistic muscle semg signals | |
CN113040738B (zh) | 血压检测装置 | |
CN113812933A (zh) | 基于可穿戴式设备的急性心肌梗死实时预警系统 | |
Zhao et al. | PVC recognition for wearable ECGs using modified frequency slice wavelet transform and convolutional neural network | |
DongseoUniversity | The deployment of novel techniques for mobile ECG monitoring | |
CN106056045B (zh) | 一种三导联重构十二导联心电信号的方法 | |
Jiang et al. | Enhanced machine learning feature selection algorithm for cardiac arrhythmia in a personal healthcare application | |
Kew et al. | Wearable patch-type ECG using ubiquitous wireless sensor network for healthcare monitoring application | |
Rashkovska et al. | Clustering of heartbeats from ECG recordings obtained with wireless body sensors | |
Revanth et al. | Detection of QRS Complexes in ECG Signal | |
CN111543977B (zh) | 基于12导联静息心电图的多级联人工智能排阴方法 |
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: 20190712 |
|
RJ01 | Rejection of invention patent application after publication |