CN1887223A - 心脏状态的实时趋势动力学特征分析方法 - Google Patents
心脏状态的实时趋势动力学特征分析方法 Download PDFInfo
- Publication number
- CN1887223A CN1887223A CN 200610052490 CN200610052490A CN1887223A CN 1887223 A CN1887223 A CN 1887223A CN 200610052490 CN200610052490 CN 200610052490 CN 200610052490 A CN200610052490 A CN 200610052490A CN 1887223 A CN1887223 A CN 1887223A
- Authority
- CN
- China
- Prior art keywords
- characteristic
- time
- parameter
- scatterplot
- heart
- 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.)
- Granted
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种心脏状态的实时趋势动力学特征分析方法。由高输入阻抗、低噪声、高共模抑制比多通道心电放大器,12位A/D转换板获得心电波形;采用时空相关技术和散点图技术,形成时间序列相关心率变异时空散点图;自动和人工交互提取时序相关特征参数,短时-实时特征转变图谱参数,提取特征图谱时空量化参数指标;对上述图谱及其量化参数指标进行人工神经网络分类:参数指标在两个分别由9个准静态时序离散度特征参数集合和5个动态趋势特征参数集合构成的9维与5维空间来刻画时序心率变异的动力学特征和相关不变性特征,存在两个独立曲面及其交界来表征心脏功能的正常,异常和关键的临界状态。用于早期心脏病理动力学趋势的分类预测。
Description
技术领域
本发明涉及一种用于诊断目的的探测测量方法。尤其涉及一种基于时序相关散点图技术的心脏疾病早期分析和心脏动力学特征识别方法。
背景技术
反映心率变异性的Poincare散点图的图形特征参数与许多心脏疾病包括心室传导阻滞、心肌缺血、心肌损伤、交感-副交感调节和体液调节能力有着密切关系。用心率变异性检测评估自律心率控制能力的方法已经成为心脏血管功能的关键信息来源。特别是,它已经作为急性心肌缺血患者的阴性预后的独立判断标识。多年来国内外学者通过时域分析和频域分析,已经证实心率变异性对心脏疾病动态趋势变化具有很强和简洁的敏感性表征,因此在临床心卒和猝死具有早期分析的预测的重要意义。但是没有明确地引入时间信息标尺,不能划分短时心脏活动的状态,不能以时域尺度刻画心脏状态的稳定行和心脏状态的动力学适应性。因而不能全面的描述心脏的当前的生理/病理特征和具有预后行的趋势特征。
小波变换(Wavelet transform,WT)是一种线性运算,它对信号进行不同尺度的分解,能将各种交织在一起的不同频率组成的混合信号分解成不同频率的块信号,可有效地应用于信号分离和ECG复杂波的识别,提高时频两域的分辨率等。将它与时域极值方法结合起来分析可以获得精确的RR间期参数。
神经网络分析系统,以模仿人脑功能,完成类似人脑的工作,是生理学上的真实人脑神经网络的结构和功能,以及若干基本特性的某种理论抽象、简化和模拟而构成的一种信息处理系统。从系统观点来看,这种人工神经网络是有大量神经元通过极其丰富和完善的联接而构成的自适应非线性动态系统。由于神经元之间有着不同的联接方式,所以可组成不同结构形态的神经网络系统。误差反向传播算法(BP算法)网络以及自适应谐振理论算法(ART2)网络均有很好的系统辨识与分类能力。
发明内容
本发明的目的在于:提供一种基于时间序列的心率变异离散图,实现心脏状态的实时趋势动力学特征分析方法。
本发明采用的技术方案如下:
1、用体表电极和高输入阻抗、低噪音、高共模抑制比以及宽带多通道心电放大器得到心电模拟信号,再由12位A/D转换器获得数字心电波形;采用小波变换和极值获得QRS复杂波形位置,并通过RR波峰的时间间隔形成RR(n)/RR(n+1)的表达心率变异率的散点图;再从散点图中生成不同时间间隔的短时散点图形成时序散点图,以分钟作为短时散点图的最小间隔,通过时序散点图提取丰富的心脏动力学状态特征参数和趋势特征参数;时间序列节点内离散度参数;时间序列节点间离散度参数;时间序列节点内离散度描述了心脏动力学状态的稳定度;时间序列节点间离散度距离描述了心脏动力学状态的迁移转换特性;状态稳定度和迁移转换特性表征了心脏的适应性以及体液、交感-副交感神经系统活动影响。
2、对上述时间序列节点内离散度参数,时间序列节点间离散度参数进行人工神经网络的分类分析,分析结果由14个特征参数分别构成9维和5维空间,综合决定心率变异离散度的性质,存在一个14维的曲面,14维曲面将分为3个部分:心脏功能的正常、异常以及过渡区域;BP网络中确定网络隐层为一层;网络的输入层的神经元个数N由上述特征参数的个数选定了总共9个特征参数,这样输入层神经元需要9个;隐层神经元个数M的选择依据经验公式
其中P是训练库的样本个数,隐层神经元个数与训练的样本个数相关,式中M与P均跟网络拟合的非线性关系的复杂程度相关;网络的输出层神经元个数K由需要分类的种数决定;
本发明的特点是:将多通道心电波形的检测、先进的极值定位、小波变换、独特的自动时序心率变异散点图形识别及其参数化提取与人工神经网络的分类分析方法进行有机结合,多个时序心率变异离散图参数在多维空间中拟合多维曲面进行分类的创新方法对临床多种心脏病,特别是心肌缺血心脏病的分析、状态变迁趋势和预测具有重要的意义。用于早期心脏病理动力学趋势的分类预测。
附图说明
图1.POINCARE心率变异散点图及其轮廓表示;
图2.心率变异散点图构造方法;
图3.散点集合准静态图形参数描述图;
图4.准静态参数BP逻辑图;
图5.动态参数BP逻辑图;
图6.ART2神经网络逻辑图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
1.RR心电图间期的获得。心电波形由如下装置获得,装置由体表电极和高输入阻抗、低噪声、高共模抑制比心电放大器;12位A/D转换器;XSCAL ARM嵌入式处理系统;
2.时序心率离散变异散点分析方法以及提取心脏动力学特征参数和趋势特征参数的分析方法:
a)小波变换从3/5导联ECG的QRS波群的提取出RR间期。如图2右子图所示计算RR间期,如图2左子图所示计算得到散点图坐标点,当前RR间隔(RRn+1)对前一个RR间期(RRn+1)构成一个序偶二维坐标(RRn+1,RRn),一系列的序偶二维坐标点集合构成的反映心率变动性(HRV)的可视化散点图又称为Poincare散点图。
b)构造生成不同时间间隔的短时散点图形成时序散点图(以分钟作为短时散点图的最小间隔),并且如图1和图2所示可采用全自动方法和人机交互方法计算确定散点集合轮廓特征参数(其轮廓线可以按照图1A的包络曲线或者按照图1B的平行四边形周线计算):散点集合质心位置、正交轴向角、轮廓正交轴向顶点、基于质心正交轴向离散度、散点集合密度(轮廓闭合曲线面积)、最大/最小RR、最大/最小ΔRR等9个主要参数以及如图1和图3的轮廓内疏密点集子聚合区形状、位置、数量等多个辅助参数。其中每个点集子聚合区也可以计算9个主要准静态特征参数。9个主要准静态特征参数在多维空间是相互独立(正交必要但不充分完备);
c)序心率变异散点图的表示方法:颜色维标记时间序列与二维心率变异散点图构成颜色三维时序心率变异散点图;或者与二维心率变异散点图坐标正交的时间轴构成空间三维时序心率变异散点图。其中散点集合可以用集合质心与集合轮廓表示成简洁三维时序心率变异散点图。
d)提取心脏动力学状态特征参数及其趋势特征参数:二维空间质心运动参数(方向θ、速度ν)、散点集合密度变化率,散点集合轮廓轴向伸缩运动参数(扩张率、收缩率)、散点集合轮廓质心为圆点旋转运动参数(角速度ω)以及质心运动轨迹图谱。
3.人工神经网络分析特征参数获得病理特征分类和趋势特征分类结果。
从临床的原始动态心电数据出发得到病理特征分类和动态特征分类以及趋势预测的心电数据集合。如图4~5,准静态参数和时序动态参数每个分类各包含64个训练数据集作为人工神经网络输入数据集合。这两类相同结构都利用ART2自适应学习网络(如图6)在实际应用中分别对应于病理准静态特征分类和动态趋势特征进行分类分析。ART2神经网络实际上按照训练数据在模式空间中相聚的距离远近进行分类,如图在训练过程中,包括准静态特征参数描述和动力学特征及趋势特征参数描述的两个训练集合,分别进行训练得到两个特征空间集。这两个空间集合就是分类的结果。
在实际分类分析工作中,如果出现准静态特征参数描述误差,通过匹配器触发个体差异性自学习事件;如果出现准动力学特征及趋势特征参数描述误差,触发指导教师存在的人工神经网络自学习过程,自学习过程结果获得病理特征演化的可视化预测信息并进行可视化病理变迁预测。
Claims (2)
1、一种心脏状态的实时趋势动力学特征分析方法,其特征在于:用体表电极和高输入阻抗、低噪音、高共模抑制比以及宽带多通道心电放大器得到心电模拟信号,再由12位A/D转换器获得数字心电波形;采用小波变换和极值获得QRS复杂波形位置,并通过RR波峰的时间间隔形成RR(n)/RR(n+1)的表达心率变异率的散点图;再从散点图中生成不同时间间隔的短时散点图形成时序散点图,以分钟作为短时散点图的最小间隔,通过时序散点图提取丰富的心脏动力学状态特征参数和趋势特征参数;时间序列节点内离散度参数;时间序列节点间离散度参数;时间序列节点内离散度描述了心脏动力学状态的稳定度;时间序列节点间离散度距离描述了心脏动力学状态的迁移转换特性;状态稳定度和迁移转换特性表征了心脏的适应性以及体液、交感—副交感神经系统活动影响。
2、根据权利要求1所述的一种心脏状态的实时趋势动力学特征分析方法,其特征在于:对上述时间序列节点内离散度参数,时间序列节点间离散度参数进行人工神经网络的分类分析,分析结果由14个特征参数分别构成9维和5维空间,综合决定心率变异离散度的性质,存在一个14维的曲面,14维曲面将分为3个部分:心脏功能的正常、异常以及过渡区域;BP网络中确定网络隐层为一层;网络的输入层的神经元个数N由上述特征参数的个数选定了总共9个特征参数,这样输入层神经元需要9个;隐层神经元个数M的选择依据经验公式
其中P是训练库的样本个数,隐层神经元个数与训练的样本个数相关,式中M与P均跟网络拟合的非线性关系的复杂程度相关;网络的输出层神经元个数K由需要分类的种数决定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100524905A CN100415159C (zh) | 2006-07-17 | 2006-07-17 | 心脏状态的实时趋势动力学特征分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100524905A CN100415159C (zh) | 2006-07-17 | 2006-07-17 | 心脏状态的实时趋势动力学特征分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1887223A true CN1887223A (zh) | 2007-01-03 |
CN100415159C CN100415159C (zh) | 2008-09-03 |
Family
ID=37576504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100524905A Expired - Fee Related CN100415159C (zh) | 2006-07-17 | 2006-07-17 | 心脏状态的实时趋势动力学特征分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100415159C (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101427917B (zh) * | 2008-09-08 | 2011-02-09 | 电子科技大学 | 基于固有趋势子序列模式分解的出现新的心脏活动趋势的ecg信号获取方法 |
CN101675456B (zh) * | 2007-05-04 | 2013-02-13 | 皇家飞利浦电子股份有限公司 | 心脏轮廓繁殖 |
CN103549949A (zh) * | 2013-10-21 | 2014-02-05 | 华南理工大学 | 基于确定学习理论的心肌缺血辅助检测方法 |
CN104318061A (zh) * | 2014-09-25 | 2015-01-28 | 北京国双科技有限公司 | 用于散点图的数据展示处理方法及装置 |
CN104970789A (zh) * | 2014-04-04 | 2015-10-14 | 中国科学院苏州纳米技术与纳米仿生研究所 | 心电图分类方法及系统 |
CN106073755A (zh) * | 2016-05-27 | 2016-11-09 | 成都信汇聚源科技有限公司 | 一种微型动态心电监测设备上房颤自动识别的实现方法 |
CN107260161A (zh) * | 2017-07-18 | 2017-10-20 | 上海图灵医疗科技有限公司 | 一种心电动力学数据量化分析方法 |
CN107438399A (zh) * | 2015-04-08 | 2017-12-05 | 皇家飞利浦有限公司 | 心血管恶化预警评分 |
CN107898454A (zh) * | 2017-11-13 | 2018-04-13 | 湖北科技学院 | 一种非线性洛伦兹散点图形态学计算方法 |
CN108577831A (zh) * | 2018-03-19 | 2018-09-28 | 武汉海星通技术股份有限公司 | 单导心贴数据长程监控诊断系统及其处理方法 |
CN108784680A (zh) * | 2018-03-19 | 2018-11-13 | 武汉海星通技术股份有限公司 | 基于散点图的心电图智能分析方法及系统 |
CN109431488A (zh) * | 2018-09-28 | 2019-03-08 | 上海优加利健康管理有限公司 | 一种心电图动态实时分析数据中典型数据的识别方法 |
CN109493933A (zh) * | 2018-08-08 | 2019-03-19 | 浙江大学 | 一种基于注意力机制的不良心血管事件的预测装置 |
CN110133629A (zh) * | 2018-08-31 | 2019-08-16 | 深圳迈睿智能科技有限公司 | 微波探测器和微波探测方法及智能设备 |
CN113273986A (zh) * | 2021-05-19 | 2021-08-20 | 于杨 | 一种心率变异性分析方法、装置、智能终端及存储介质 |
CN114767123A (zh) * | 2022-03-30 | 2022-07-22 | 合肥工业大学 | 一种基于多模态神经网络的心电散点图异常识别分类方法 |
-
2006
- 2006-07-17 CN CNB2006100524905A patent/CN100415159C/zh not_active Expired - Fee Related
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101675456B (zh) * | 2007-05-04 | 2013-02-13 | 皇家飞利浦电子股份有限公司 | 心脏轮廓繁殖 |
CN101427917B (zh) * | 2008-09-08 | 2011-02-09 | 电子科技大学 | 基于固有趋势子序列模式分解的出现新的心脏活动趋势的ecg信号获取方法 |
CN103549949A (zh) * | 2013-10-21 | 2014-02-05 | 华南理工大学 | 基于确定学习理论的心肌缺血辅助检测方法 |
CN103549949B (zh) * | 2013-10-21 | 2015-04-29 | 华南理工大学 | 基于确定学习理论的心肌缺血辅助检测方法 |
CN104970789A (zh) * | 2014-04-04 | 2015-10-14 | 中国科学院苏州纳米技术与纳米仿生研究所 | 心电图分类方法及系统 |
CN104318061B (zh) * | 2014-09-25 | 2018-02-02 | 北京国双科技有限公司 | 用于散点图的数据展示处理方法及装置 |
CN104318061A (zh) * | 2014-09-25 | 2015-01-28 | 北京国双科技有限公司 | 用于散点图的数据展示处理方法及装置 |
CN107438399A (zh) * | 2015-04-08 | 2017-12-05 | 皇家飞利浦有限公司 | 心血管恶化预警评分 |
CN106073755A (zh) * | 2016-05-27 | 2016-11-09 | 成都信汇聚源科技有限公司 | 一种微型动态心电监测设备上房颤自动识别的实现方法 |
CN107260161A (zh) * | 2017-07-18 | 2017-10-20 | 上海图灵医疗科技有限公司 | 一种心电动力学数据量化分析方法 |
WO2019015017A1 (zh) * | 2017-07-18 | 2019-01-24 | 上海图灵医疗科技有限公司 | 一种心电动力学数据量化分析方法 |
CN107260161B (zh) * | 2017-07-18 | 2019-02-19 | 上海图灵医疗科技有限公司 | 一种心电动力学数据量化分析方法 |
CN107898454A (zh) * | 2017-11-13 | 2018-04-13 | 湖北科技学院 | 一种非线性洛伦兹散点图形态学计算方法 |
CN108784680A (zh) * | 2018-03-19 | 2018-11-13 | 武汉海星通技术股份有限公司 | 基于散点图的心电图智能分析方法及系统 |
CN108577831A (zh) * | 2018-03-19 | 2018-09-28 | 武汉海星通技术股份有限公司 | 单导心贴数据长程监控诊断系统及其处理方法 |
CN108577831B (zh) * | 2018-03-19 | 2021-02-09 | 武汉海星通技术股份有限公司 | 单导心贴数据长程监控诊断系统及其处理方法 |
CN109493933A (zh) * | 2018-08-08 | 2019-03-19 | 浙江大学 | 一种基于注意力机制的不良心血管事件的预测装置 |
CN110133629A (zh) * | 2018-08-31 | 2019-08-16 | 深圳迈睿智能科技有限公司 | 微波探测器和微波探测方法及智能设备 |
CN110187397A (zh) * | 2018-08-31 | 2019-08-30 | 深圳迈睿智能科技有限公司 | 用于车辆的微波探测系统、装置及其应用 |
CN110133629B (zh) * | 2018-08-31 | 2022-03-22 | 深圳迈睿智能科技有限公司 | 微波探测器和微波探测方法及智能设备 |
CN109431488A (zh) * | 2018-09-28 | 2019-03-08 | 上海优加利健康管理有限公司 | 一种心电图动态实时分析数据中典型数据的识别方法 |
CN113273986A (zh) * | 2021-05-19 | 2021-08-20 | 于杨 | 一种心率变异性分析方法、装置、智能终端及存储介质 |
CN114767123A (zh) * | 2022-03-30 | 2022-07-22 | 合肥工业大学 | 一种基于多模态神经网络的心电散点图异常识别分类方法 |
CN114767123B (zh) * | 2022-03-30 | 2024-04-26 | 合肥工业大学 | 一种基于多模态神经网络的心电散点图异常识别分类方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100415159C (zh) | 2008-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100415159C (zh) | 心脏状态的实时趋势动力学特征分析方法 | |
CN109583346A (zh) | 基于lstm-fc的脑电信号特征提取与分类识别方法 | |
CN104166548B (zh) | 基于运动想象脑电数据的深度学习方法 | |
CN106293057A (zh) | 基于bp神经网络的手势识别方法 | |
CN109846469A (zh) | 一种基于卷积神经网络的非接触式心率测量方法 | |
CN110974214A (zh) | 一种基于深度学习的自动心电图分类方法、系统及设备 | |
CN111476158B (zh) | 一种基于pso-pca-svm的多通道生理信号体感手势识别方法 | |
CN104035563A (zh) | 基于w-pca和无监督ghsom的脑电信号识别方法 | |
CN110826625A (zh) | 一种基于表面肌电信号的手指手势分类方法 | |
CN109033415A (zh) | 一种基于拉普拉斯特征映射的多维神经信号的降维和可视化方法 | |
CN106529570A (zh) | 基于深度脊波神经网络的图像分类方法 | |
CN116503673B (zh) | 一种基于心电图的心律失常识别检测方法及系统 | |
CN114159079B (zh) | 基于特征提取和gru深度学习模型的多类型肌肉疲劳检测方法 | |
CN114925734A (zh) | 一种基于神经拟态计算的在线神经元分类方法 | |
CN106845348A (zh) | 一种基于手臂表面肌电信号的手势识别方法 | |
CN114767123B (zh) | 一种基于多模态神经网络的心电散点图异常识别分类方法 | |
CN109948465A (zh) | 一种基于多目标的表面肌电信号分类方法 | |
CN111783669B (zh) | 一种面向个人用户的表面肌电信号分类识别方法 | |
Zhao et al. | Artificial neural networks based on fractal growth | |
Lata et al. | Disease classification using ECG signals based on R-peak analysis with ABC and ANN | |
CN114343679A (zh) | 基于迁移学习的表面肌电信号上肢动作识别方法及系统 | |
Wang et al. | Study on intelligent syndrome differentiation in traditional Chinese medicine based on information fusion technology | |
Zhou et al. | Gesture Recognition Based on Attention and Dilated Causal Convolution | |
Krishnapriya et al. | Surface Electromyography based Hand Gesture Signal Classification using 1D CNN | |
Shatnyi et al. | Electrocardiogram Processing System Design with Parallel Computing and Memory Transferring Using Fuzzy ART Neural Network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080903 Termination date: 20130717 |