CN101176660A - 一种心血管系统联合变异性指标的检测方法和装置 - Google Patents
一种心血管系统联合变异性指标的检测方法和装置 Download PDFInfo
- Publication number
- CN101176660A CN101176660A CNA2007101151437A CN200710115143A CN101176660A CN 101176660 A CN101176660 A CN 101176660A CN A2007101151437 A CNA2007101151437 A CN A2007101151437A CN 200710115143 A CN200710115143 A CN 200710115143A CN 101176660 A CN101176660 A CN 101176660A
- Authority
- CN
- China
- Prior art keywords
- variability
- time series
- signal
- cardiovascular system
- layer element
- 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
本发明提供了一种心血管系统联合变异性指标的检测方法和装置,检测方法包括以下步骤:(1)采集受检测者的心电、心音、挠动脉脉搏信号;(2)对三路信号进行模拟-数字转换,形成信号波形图;(3)识别和提取三路信号各自的特征点;(4)构造联合变异性指标的心电间期时间序列RR,电机械延迟时间序列和脉搏波传播时间序列;(5)检验各时间序列的有效性;(6)计算心率变异性、电机械延迟变异性和脉搏波传播时间变异性;(7)计算心血管联合变异性指标AV。检测装置包含心电信号检测模块、心音信号检测模块、挠动脉脉搏信号检测模块、模拟-数字转换装置和计算机。本发明能够较好地反映人体心血管系统的不同生理、病理状态,区分度高,具有广泛的临床应用价值。
Description
技术领域
本发明涉及一种人体心血管系统联合变异性指标AV的检测方法和装置,特别是对人体进行无创心血管系统功能状态检测的方法和装置,属于心血管系统功能状态检测技术领域。
背景技术
人体心血管系统是一个复杂的非线性系统,其中的诸多信号如心电信号、心音信号、脉搏信号蕴含着丰富的非线性动力学信息。非线性动力学信息的多少与人体心血管系统的功能状态密切相关,同时反应了自主神经对心血管系统的调控功能,也反应了呼吸因素、体液因素等与心血管系统的相互作用。
心血管系统信号的变异性是其非线性动力学信息的一个主要表现,研究信号的变异性是分析心血管系统非线性的一个重要方法,以此可以评价心血管系统的不同生理、病理状态。因此,对心血管系统信号的变异性研究引起了相关人员极大的兴趣,并有人尝试将研究成果应用到临床。连续窦性心搏间瞬时心率的微小涨落称为心率变异性(HRV),有关HRV的研究开展的相对较早,HRV信号中已被证实蕴含了有关心血管系统及神经、呼吸、体液调节等的大量信息,对这些信息的提取分析可以定量评估心交感神经和迷走神经的紧张性和均衡性及其对心血管活动的影响,HRV分析对一些疾病的早期诊断、治疗及预后评估等均具有重要意义。现已证明,许多疾病如心肌梗死、充血性心力衰竭、心源性猝死、冠心病、隐性心肌缺血、糖尿病、高血压等疾病的发生,在HRV分析结果上都表现出明显的特征。
由于心脏电活动和机械活动的非同步性,存在电机械延迟,在同一心动周期内,心脏电机械延迟可以由第一心音(S1)的二尖瓣关闭成分(A1)出现时间减去心电Q波起始点出现时间获得,这样连续的心脏电机械延迟构成了电机械延迟时间序列QA1,QA1的变异性反映了心脏细胞的传导性能和心脏自律细胞的协调性能。由于脉搏波由心脏传播到挠动脉需要一定的时间,这个时间称为脉搏波传播时间PWTT(Pulse Wave Translation Time),在同一心动周期内,PWTT可由挠动脉脉搏波切迹点In出现时间减去第二心音(S2)的主动脉瓣关闭成分(A2)出现时间获得,这样连续的脉搏波传播时间构成PWTT序列,PWTT的变异性反映了血管的硬化度和顺应性。
从心率变异性角度评价心脏功能的研究已开展多年,但到目前为止,还未见有利用心率变异性、电机械延迟变异性、脉搏波传播时间变异性三个方面的信息来全面评价人体心血管系统功能状态的研究报道。
变异性分析常用的方法有两类:线性分析方法和非线性分析方法,其中线性分析方法又分为时域分析方法、频域分析方法和时频分析方法等。(1)时域分析指标:主要包括检测时间序列(RR序列、QA1序列、PWTT序列)的标准差SDNN、每5min时间序列平均值的标准差SDANN、时间序列相邻间期差值的均方根值RMSSD、时间序列间期直方图、时间序列间期分布图等;(2)频域分析指标:主要包括检测时间序列的总功率、极低频、低频、高频成分含量、低频与高频成分含量比值等;(3)时频分析指标:主要是检测时间序列谱峰的动态特性。非线性分析方法目前有散点图分析法、分形维数分析法、复杂度分析法、近似熵分析法、模式熵分析法和符号动力学分析法等,但是真正在临床上有较多应用的是散点图分析方法。线性分析方法由于并不符合心血管系统本身是非线性系统的事实,因此分析的结果特异性比较低,并对外界噪声比较敏感。散点图分析方法主要依靠医生肉眼观察,得出一些定性或半定量指标,指标物理意义不是很明确,难以进行准确诊断。其它的非线性分析方法分别从不同侧面显示了其理论价值,但是距离应用到临床短时心血管系统的变异性分析上还有差距。
发明内容
本发明针对现有从心率变异性角度评价心脏功能技术存在的不足,提供一种心血管系统联合变异性指标的检测方法,该检测方法综合考虑心率变异性、电机械延迟变异性、脉搏波传播时间变异性三个方面信息,全面评价人体心血管系统的功能状态。同时提供一种心血管系统联合变异性指标的检测装置。
本发明的心血管系统联合变异性指标的检测方法包括以下步骤:
(1)在同一心动周期内同步采集受检测者的心电、心音、挠动脉脉搏信号;
(2)对三路信号进行模拟-数字转换,形成信号波形图;
(3)识别和提取三路信号各自的特征点:心电信号的P波、QRS波群和T波,心音信号的第一心音的二尖瓣关闭成分、第二心音的主动脉瓣关闭成分,挠动脉脉搏信号的脉搏起始点和重搏波切迹点;
(4)根据提取的特征点构造联合变异性指标的时间序列:心电信号波形中相邻的R波峰之间的心电间期时间序列RR,用同一心动周期内第一心音的二尖瓣关闭成分出现时间减去心电Q波起始点出现时间获取的电机械延迟时间序列QA1,用同一心动周期内挠动脉脉搏波切迹点出现时间减去第二心音的主动脉瓣关闭成分出现时间获取的脉搏波传播时间序列PWTT;
(5)检验各时间序列的有效性:若心电间期时间序列RR、电机械延迟时间序列QA1、脉搏波传播时间序列PWTT中相同位置的值RRi,QA1i,PWTTi(i=1,2,…,N)都是同一心动周期内的时间间隔信息,则获取的时间序列为有效序列,否则时间序列无效,需重新获取时间序列;
(6)根据获取的三个时间序列,采用尺度熵算法计算心率变异性HRV、电机械延迟变异性QA1V和脉搏波传播时间变异性PWTTV;
(7)根据心率变异性HRV、电机械延迟变异性QA1V和脉搏波传播时间变异性PWTTV计算心血管联合变异性指标AV。
下面给出本发明计算心血管系统联合变异性指标AV的具体过程。
首先,在同一心动周期内,同步采集一受检测者的心电、心音、挠动脉脉搏信号,采集时间持续5~10分钟,三路信号采用相同的采样率进行模拟-数字转换。其次,利用各自的特征提取算法由计算机自动识别出信号特征点,构造联合变异性指标AV的时间序列,具体有三个:1.心电信号相邻的R波峰之间的心电间期时间序列RR,RR={RR1,RR2,…,RRN};2.在同一心动周期内,用第一心音(S1)的二尖瓣关闭成分(A1)出现时间减去心电Q波起始点出现时间获取的电机械延迟时间序列QA1,QA1={QA11,QA12,…,QA1N};3.在同一心动周期内,用挠动脉脉搏波切迹点In出现时间减去第二心音(S2)的主动脉瓣关闭成分(A2)出现时间获取的脉搏波传播时间序列PWTT,PWTT={PWTT1,PWTT2,…,PWTTN}。再次,对获取的三个时间序列进行有效性检验:若心电间期时间序列RR、电机械延迟时间序列QA1、脉搏波传播时间序列PWTT中相同位置的值RRi,QA1i,PWTTi(i=1,2,…,N)都是同一心动周期内的时间间隔信息,则获取的时间序列为有效序列,否则时间序列无效,需重新获取时间序列。最后,采用尺度熵算法计算三个时间序列各自的变异性HRV、QA1V、PWTTV。下面以心电信号为例说明如何利用尺度熵算法计算HRV。
(1)构造m维相空间:对于心电间期时间序列RR={RR1,RR2,…,RRN},取m个连续点组成-m维矢量
X(i)=[RRi,RRi+1,…,RRi+m-1]
i=1,2,…,N-m+1,共有N-m+1个m维矢量;
(2)矢量符号化:对于每一个矢量X(i),首先计算其均值u(i)和均方差σ(i),即
根据均值u(i)和均方差σ(i)将矢量X(i)转换为符号矢量S(i)=[si,si+1,…,si+m-1],
其中i=1,2,…,N-m+1,j=0,1,…,m-1,λ为一个参数;
(3)计算尺度熵:每一个m维符号矢量S(i)都代表了心电间隔序列RR的一种波动模式,统计N-m+1个符号矢量S(i)中的不同波动模式为{S(k)},k=1,2,…,K,计算每一种模式S(k)的分布概率P(k)
其中Nk为N-m+1个符号矢量中出现波动模式S(k)的次数;
计算RR序列的尺度熵为:
采用H(m)值来代表心率变异性HRV。
同理,电机械延迟变异性QA1V和脉搏波传播时间变异性PWTTV也采用上述类似的步骤获得。
利用计算得出的三个变异性HRV、QA1V、PWTTV,由下面公式计算联合变异性指标AV:
其中a,b,c为权系数,代表HRV、QA1V、PWTTV三种变异性对联合变异性指标AV贡献度的大小。权系数a,b,c采用层析分析法确定,具体步骤为:
(1)以HRV、QA1V、PWTTV为对象层元素,总数为N;以高血压、冠心病、心肌梗塞、心率不齐、动脉硬化、年龄、性别为准则层元素,总数为K;以AV为目标层元素;
(2)判断第i个和第j个对象层元素对准则层元素k的重要程度相对比值aij,其中i,j=1,2,…,N,k=1,2,…,K,构成成对比较矩阵Ak,并对矩阵Ak进行一致性检验:计算矩阵Ak中最大特征值λkmax,若λkmax≥0.9×N,则通过一致性检验,转向步骤(3);否则重新获取成对比较矩阵Ak;
(3)计算Ak中对应最大特征值λkmax的特征向量wk,wk为行向量,表示对象层元素对准则层元素k的权重,依次获取的K个wk构成权矩阵w,w=[w1;w2;…;wK];
(4)判断第s个和第p个准则层元素对目标层元素的重要程度相对比值bsp,构成成对比较矩阵B,并对矩阵B进行一致性检验:计算矩阵B中最大特征值λmax,若λmax≥0.9×K,则通过一致性检验,转向(5);否则重新获取成对比较矩阵B;
(5)计算B中对应最大特征值λmax的特征向量v,v为行向量,表示准则层元素对目标层元素的权重;
本发明的心血管系统联合变异性指标的检测装置采用以下技术解决方案:
该心血管联合变异性指标检测装置包含心电信号检测模块、心音信号检测模块、挠动脉脉搏信号检测模块、模拟-数字转换装置和计算机;心电信号检测模块、心音信号检测模块和挠动脉脉搏信号检测模块均与模拟-数字转换装置连接,模拟-数字转换装置与计算机连接;心电信号检测模块包括采集电极和滤波、去噪、放大电路,用于采集受检测者的心电信号,并完成信号的滤波、去噪、放大;心音信号检测模块包括心音传感器和滤波、去噪、放大电路,用于采集受检测者的心音信号,并完成对信号的滤波、去噪、放大;挠动脉脉搏信号检测模块包括脉搏波传感器、固定支架和滤波、去噪以及放大电路,用于采集受检测者的挠动脉脉搏信号,脉搏波传感器安装在固定支架上。
各模块中的滤波、去噪、放大电路采用现有通用电路。
固定支架由固定座和滑动块组成,固定座上设有导轨,滑动块安装在固定座的导轨上,滑动块上部安装有进给螺栓,脉搏波传感器安放在滑动块上。滑动块带动传感器沿固定座导轨水平运动,进给螺栓与脉搏波传感器相连,进给螺栓转动使脉搏波传感器上下运动。脉搏波传感器的水平运动使其准确定位在检测点正上方,上下运动使其与受检测者皮肤保持适当压力。
模拟-数字转换装置采用现有四通道同步数据采集卡。
与计算机相连有一打印机,用于打印分析报告。
与计算机相连有一组指示灯,至少包含三个指示灯,分别指示装置处于准备就绪状态、正在采集信号状态、采集结束状态。
与计算机相连有一启动开关,闭合启动开关,实现对检测装置的上电和预热;一工作按钮,闭合工作按钮,实现信号的整个自动采集过程。
本发明的心血管联合变异性指标检测方法和装置是使用一种联合变异性指标AV对心血管系统功能状态进行评价,AV包含心率变异性、电机械延迟变异性和脉搏波传播时间变异性三个有机组成部分,全面反映了自主神经系统对心血管系统的调控作用、心脏细胞的传导性能、心脏自律细胞的协调性能、血管的顺应性和硬化度,能够实现对心血管系统功能状态的整体评价。经过了严格的临床试验,验证了联合变异性指标AV对人体心血管系统的功能状态有较好的综合评价效果。
附图说明
图1是本发明的心血管系统联合变异性指标的检测方法的流程图。
图2是在同一心动周期内同步采集到的心电、心音、挠动脉脉搏信号的波形示意图。
图3是根据提取的特征点获取的用于构造联合变异性指标的心电间期时间序列RR、心脏电机械延迟时间序列QA1和脉搏波传播时间序列PWTT的示意图。
图4是实现本发明方法的心血管系统联合变异性指标的检测装置的结构原理示意图。
图5是本发明的检测装置中挠动脉脉搏传感器固定支架的结构示意图。
图6是本发明的检测装置的操作流程图。
其中:1、心电信号检测模块,2、心音信号检测模块,3、挠动脉脉搏检测模块,4、模拟-数字转换装置,5、计算机,6、打印机,7、一组指示灯,8、启动开关、9、工作按钮,10、右腕电极,11、右腿电极,12、左腿电极,13、心音传感器,14、脉搏波传感器,15、固定座,16、滑动块,17、卡隼,18、进给螺栓,19、挂带柱,20、固定绷带,21a、轴向定位标志,21b、轴向定位标志,22、径向定位标志。
具体实施方式
图1给出了本发明心血管联合变异性指标AV的检测方法的流程。首先在同一心动周期内同步采集受检测者的心电、心音、挠动脉脉搏信号,对三路信号进行滤波、去噪、放大以及模拟-数字转换,形成如图2所示的信号波形图,图2中由上到下依次为心电信号、心音信号、挠动脉脉搏信号,并于图上标注出了心电信号的P波、QRS波群和T波,心音信号的第一心音(S1)的二尖瓣关闭成分(A1)、第二心音(S2)的主动脉瓣关闭成分(A2),挠动脉脉搏信号的脉搏起始点u和重搏波切迹点In,将上述标注点作为信号的特征点。各信号波形是利用图4所示实现本发明的装置采集到的。根据图2提取的特征点获取用于构造联合变异性指标分析的心电间期时间序列RR、电机械延迟时间序列QA1和脉搏波传播时间序列PWTT如图3所示。具体构造方法为:心电信号相邻的R波峰之间的时间间隔构成心电间期时间序列RR,用同一心动周期内第一心音(S1)的二尖瓣关闭成分(A1)出现时间减去心电Q波起始点出现时间获取电机械延迟时间序列QA1;用同一心动周期内挠动脉脉搏波切迹点In出现时间减去第二心音(S2)的主动脉瓣关闭成分(A2)出现时间获取脉搏波传播时间序列PWTT。
检验各时间序列的有效性:若心电间期时间序列RR、电机械延迟时间序列QA1、脉搏波传播时间序列PWTT中相同位置的值RRi,QA1i,PWTTi(i=1,2,…,N)都是同一心动周期内的时间间隔信息,则获取的时间序列为有效序列,否则时间序列无效,需重新获取时间序列。
根据获取的三个时间序列,采用尺度熵算法计算心率变异性HRV、电机械延迟变异性QA1V和脉搏波传播时间变异性PWTTV。根据心率变异性HRV、电机械延迟变异性QA1V和脉搏波传播时间变异性PWTTV计算心血管联合变异性指标AV,由下面的公式计算:
其中a,b,c为权系数,代表HRV、QA1V、PWTTV三种变异性对联合变异性指标AV贡献度的大小。
图4给出了实现本发明方法的装置的结构原理图,该心血管联合变异性指标检测装置主要包括心电信号检测模块1、心音信号检测模块2、挠动脉脉搏检测模块3、模拟-数字转换装置4和计算机5,与计算机5连接有打印机6、一组指示灯7、启动开关8和工作按钮9。心电检测模块1主要包含滤波、去噪、放大电路和三个检测电极:置于受检测者右手腕关节处的右腕电极10,置于受检测者右腿踝关节处的右腿电极11,置于受检测者左腿踝关节处的左腿电极12,三个检测电极构成心电II导联检测方式。心音信号检测模块2主要包含一心音传感器13和滤波、去噪、放大电路。挠动脉脉搏检测模块3主要包含一置于受检测者左手挠动脉处的脉搏波传感器14和滤波、去噪、放大电路。挠动脉脉搏检测模块3中的脉搏波传感器14安装在如图5所示的固定支架上。该支架由固定座15和滑动块16组成,脉搏波传感器14固定在滑动块16内,滑动块16两侧分别设有一个卡隼17,按住卡隼17可以使滑动块16沿固定座15上的导轨水平方向滑动。在固定座15一侧设有三个定位标志:轴向定位标志21a和21b,径向定位标志22。检测时,调节支架的位置,使挠动脉轴向与固定座15上的轴向定位标志21a和21b连线对齐,脉搏波搏动点径向投影对准径向定位标志22。在滑动块16的上方设有一进给螺栓18,转动该螺栓可以改变传感器14的探头与受检测者皮肤的接触压力。固定座15上最外侧设有挂带柱19,其上安装有固定绷带20。模拟-数字转换装置4与计算机5连接,计算机5接受来自模拟-数字转换装置4的三路数字信号,将其存储在计算机5中的RAM存储器中,并通过调用计算机5的ROM存储器中的程序实现下述操作:开始信号采集、取消信号采集、暂停信号采集、存储采集数据、分析存储数据、显示分析结果,通过与计算机5连接的打印机6打印报告。通过闭合启动开关8实现对检测装置的上电和预热;通过闭合工作按钮9实现信号的整个自动采集过程。在检测过程中,与计算机5连接的一组指示灯7指示检测装置所处状态,至少包含三个指示灯,分别指示装置处于准备就绪、正在采集信号和采集结束三个状态。
图6给出了本发明的心血管系统联合变异性指标检测装置的操作流程。按照该操作流程,装置的整个检测过程持续5~10分钟,实现无创、快速、简便地对受检测者进行检测。
本发明提供的心血管联合变异性指标AV的检测方法和装置具有下述优点:
1、综合考虑心率变异性、电机械延迟变异性、脉搏波传播时间变异性三个方面的信息,提出一种心血管联合变异性指标AV,来全面评价人体心血管系统的功能状态。
2、同一心动周期内,用第一心音(S1)的二尖瓣关闭成分(A1)出现时间减去心电P波起始点出现时间获取电机械延迟时间序列,计算电机械延迟变异性,评价心脏细胞的传导性能和心脏自律细胞的协调性能;
3、同一心动周期内,用挠动脉脉搏波切迹点In出现时间减去第二心音(S2)的主动脉瓣关闭成分(A2)出现时间获取脉搏波传播时间序列,计算脉搏波传播时间变异性,评价血管的顺应性和硬化度;
4、采用尺度熵算法计算序列的变异性,该算法针对短时程时间序列(5分钟左右)即可收敛,得出的指标重复性好,对疾病的特异性高,适宜于临床应用;
5、基于联合变异性指标AV的心血管系统功能状态检测是一种无创、实时、简便的检测技术,本发明所提供装置检测得出的信号,不仅局限于检测联合变异性指标AV,还可应用于检测心脏等容收缩期、心脏射血时间等其它心血管系统功能状态参数,有较强的临床应用价值。
Claims (9)
1.一种心血管系统联合变异性指标的检测方法,包括以下步骤:
(1)在同一心动周期内同步采集受检测者的心电、心音、挠动脉脉搏信号;
(2)对三路信号进行模拟-数字转换,形成信号波形图;
(3)识别和提取三路信号各自的特征点:心电信号的P波、QRS波群和T波,心音信号的第一心音的二尖瓣关闭成分、第二心音的主动脉瓣关闭成分,挠动脉脉搏信号的脉搏起始点和重搏波切迹点;
(4)根据提取的特征点构造联合变异性指标的时间序列:心电信号波形中相邻的R波峰之间的心电间期时间序列RR,用同一心动周期内第一心音的二尖瓣关闭成分出现时间减去心电Q波起始点出现时间获取的电机械延迟时间序列QA1,用同一心动周期内挠动脉脉搏波切迹点出现时间减去第二心音的主动脉瓣关闭成分出现时间获取的脉搏波传播时间序列PWTT;
(5)检验各时间序列的有效性:若心电间期时间序列RR、电机械延迟时间序列QA1、脉搏波传播时间序列PWTT中相同位置的值RRi,QA1i,PWTTi,i=1,2,…,N,都是同一心动周期内的时间间隔信息,则获取的时间序列为有效序列,否则时间序列无效,需重新获取时间序列;
(6)根据获取的三个时间序列,采用尺度熵算法计算心率变异性HRV、电机械延迟变异性QA1V和脉搏波传播时间变异性PWTTV;
(7)根据心率变异性HRV、电机械延迟变异性QA1V和脉搏波传播时间变异性PWTTV计算心血管联合变异性指标AV。
2.根据权利要求1所述的心血管系统联合变异性指标的检测方法,其特征在于:所述采用尺度熵算法计算心率变异性HRV的步骤是:
(1)构造m维相空间:对于心电间期时间序列RR={RR1,RR2,…,RRN},取m个连续点组成一m维矢量
X(i)=[RRi,RRi+1,…,RRi+m-1]
i=1,2,…,N-m+1,共有N-m+1个m维矢量;
(2)矢量符号化:对于每一个矢量X(i),首先计算其均值u(i)和均方差σ(i),即
根据均值u(i)和均方差σ(i)将矢量X(i)转换为符号矢量S(i)=[si,si+1,…,si+m-1],
其中i=1,2,…,N-m+1,j=0,1,…,m-1,λ为一个参数;
(3)计算尺度熵:每一个m维符号矢量S(i)都代表了心电间隔序列RR的一种波动模式,统计N-m+1个符号矢量S(i)中的不同波动模式为{S(k)},k=1,2,…,K,计算每一种模式S(k)的分布概率P(k)
其中Nk为N-m+1个符号矢量中出现波动模式S(k)的次数;
计算RR序列的尺度熵为:
采用H(m)值来代表心率变异性HRV。
4.根据权利要求3所述的心血管系统联合变异性指标的检测方法,其特征在于:所述权系数a,b,c采用层析分析法确定,具体步骤为:
(1)以HRV、QA1V、PWTTV为对象层元素,总数为N;以高血压、冠心病、心肌梗塞、心率不齐、动脉硬化、年龄、性别为准则层元素,总数为K;以AV为目标层元素;
(2)判断第i个和第j个对象层元素对准则层元素k的重要程度相对比值aij,其中i,j=1,2,…,N,k=1,2,…,K,构成成对比较矩阵Ak,并对矩阵Ak进行一致性检验:计算矩阵Ak中最大特征值λkmax,若λkmax≥0.9×N,则通过一致性检验,转向步骤(3);否则重新获取成对比较矩阵Ak;
(3)计算Ak中对应最大特征值λkmax的特征向量wk,wk为行向量,表示对象层元素对准则层元素k的权重,依次获取的K个wk构成权矩阵w,w=[w1;w2;…;wK];
(4)判断第s个和第p个准则层元素对目标层元素的重要程度相对比值bsp,构成成对比较矩阵B,并对矩阵B进行一致性检验:计算矩阵B中最大特征值λmax,若λmax≥0.9×K,则通过一致性检验,转向(5);否则重新获取成对比较矩阵B;
(5)计算B中对应最大特征值λmax的特征向量v,v为行向量,表示准则层元素对目标层元素的权重;
5.一种心血管系统联合变异性指标的检测装置,其特征在于:包含心电信号检测模块、心音信号检测模块、挠动脉脉搏信号检测模块、模拟-数字转换装置和计算机;心电信号检测模块、心音信号检测模块和挠动脉脉搏信号检测模块均与模拟-数字转换装置连接,模拟-数字转换装置与计算机连接;心电信号检测模块包括采集电极和滤波、去噪、放大电路,用于采集受检测者的心电信号,并完成信号的滤波、去噪、放大;心音信号检测模块包括心音传感器和滤波、去噪、放大电路,用于采集受检测者的心音信号,并完成对信号的滤波、去噪、放大;挠动脉脉搏信号检测模块包括脉搏波传感器、固定支架和滤波、去噪以及放大电路,用于采集受检测者的挠动脉脉搏信号,脉搏波传感器安装在固定支架上。
6.根据权利要求5所述的心血管系统联合变异性指标的检测装置,其特征在于:所述固定支架由固定座和滑动块组成,固定座上设有导轨,滑动块安装在固定座的导轨上,滑动块上部安装有进给螺栓,脉搏波传感器安放在滑动块上。
7.根据权利要求5所述的心血管系统联合变异性指标的检测装置,其特征在于:与所述计算机相连有一打印机,用于打印分析报告。
8.根据权利要求5所述的心血管系统联合变异性指标的检测装置,其特征在于:与所述计算机相连有一组指示灯,至少包含三个指示灯,分别指示装置处于准备就绪状态、正在采集信号状态、采集结束状态。
9.根据权利要求5所述的心血管系统联合变异性指标的检测装置,其特征在于:与所述计算机相连有一启动开关,闭合启动开关,实现对检测装置的上电和预热;一工作按钮,闭合工作按钮,实现信号的整个自动采集过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101151437A CN100515327C (zh) | 2007-12-06 | 2007-12-06 | 一种心血管系统联合变异性指标的检测方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2007101151437A CN100515327C (zh) | 2007-12-06 | 2007-12-06 | 一种心血管系统联合变异性指标的检测方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101176660A true CN101176660A (zh) | 2008-05-14 |
CN100515327C CN100515327C (zh) | 2009-07-22 |
Family
ID=39403190
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2007101151437A Expired - Fee Related CN100515327C (zh) | 2007-12-06 | 2007-12-06 | 一种心血管系统联合变异性指标的检测方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100515327C (zh) |
Cited By (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101862201A (zh) * | 2010-06-23 | 2010-10-20 | 南京邮电大学 | 手机智能心音看诊器中心音图的表示方法 |
CN102334986A (zh) * | 2011-07-08 | 2012-02-01 | 重庆大学 | 动脉压信号中重搏切迹点识别方法 |
CN102387742A (zh) * | 2009-02-11 | 2012-03-21 | 爱德华兹生命科学公司 | 心输出量相关波形的参数检测 |
CN102551690A (zh) * | 2011-12-29 | 2012-07-11 | 乐普(北京)医疗器械股份有限公司 | 人体信号自适应分析方法 |
CN102727245A (zh) * | 2012-07-16 | 2012-10-17 | 山东大学 | 一种心音信号采集调理传输装置 |
CN103006207A (zh) * | 2013-01-11 | 2013-04-03 | 山东师范大学 | 基于形态学的t波交替散点图法的心电信号分析方法 |
CN104414679A (zh) * | 2013-09-04 | 2015-03-18 | 北京航空航天大学 | 一种中心动脉压重建中外周压力波延迟时间的确定方法 |
CN104665798A (zh) * | 2015-01-26 | 2015-06-03 | 周常安 | 血压管理装置及方法 |
CN105078424A (zh) * | 2015-08-26 | 2015-11-25 | 深圳诺康医疗设备有限公司 | 信号采集方法、脉搏监测和主监测装置及综合监测系统 |
CN105310724A (zh) * | 2015-09-28 | 2016-02-10 | 何宗彦 | 一种测量脉搏波传播速度的方法及系统 |
CN105496402A (zh) * | 2015-11-20 | 2016-04-20 | 北京理工大学 | 基于散点图和符号动力学的心电特征分析方法 |
CN107126203A (zh) * | 2016-02-28 | 2017-09-05 | 南京全时陪伴健康科技有限公司 | 人体健康信息实时测量显示系统 |
CN107174197A (zh) * | 2016-03-10 | 2017-09-19 | 深圳市理邦精密仪器股份有限公司 | 监护设备及其波形显示方法、装置 |
CN107361753A (zh) * | 2017-08-29 | 2017-11-21 | 哈尔滨理工大学 | 基于脉搏波形特征点的人体健康状态监护方法 |
CN108185996A (zh) * | 2017-12-27 | 2018-06-22 | 中国科学院深圳先进技术研究院 | 动脉血管年龄估算模型构建方法和装置 |
CN108577883A (zh) * | 2018-04-03 | 2018-09-28 | 上海交通大学 | 一种冠心病筛查装置、筛查系统以及信号特征提取方法 |
CN109044299A (zh) * | 2018-09-25 | 2018-12-21 | 北方民族大学 | 一种心血管健康检测的生理信号采集设备及方法 |
CN109222964A (zh) * | 2018-07-20 | 2019-01-18 | 广州视源电子科技股份有限公司 | 房颤检测装置及存储介质 |
CN110200617A (zh) * | 2019-06-26 | 2019-09-06 | 慧言科技(天津)有限公司 | 一种利用心率变异性定量评价疼痛的系统 |
CN110490176A (zh) * | 2019-09-03 | 2019-11-22 | 南京邮电大学 | 一种基于多时窗心音信号的多元特征提取方法 |
CN110710960A (zh) * | 2019-09-12 | 2020-01-21 | 东莞市康助医疗科技有限公司 | 一种动脉血管年龄推算方法及装置 |
CN110956406A (zh) * | 2019-12-07 | 2020-04-03 | 中国科学院心理研究所 | 一种基于心率变异性的团队协同能力的评估模型 |
CN111543978A (zh) * | 2020-05-11 | 2020-08-18 | 中北大学 | 基于心电、心音和脉搏同步诊断的方法及系统 |
CN112842350A (zh) * | 2021-02-20 | 2021-05-28 | 无锡市中健科仪有限公司 | 基于正反传播算法获得心拍的主导心率的方法 |
CN113598724A (zh) * | 2021-09-03 | 2021-11-05 | 上海市高血压研究所 | 用于评估心血管自主神经功能的装置 |
CN113712566A (zh) * | 2020-05-12 | 2021-11-30 | 深圳市科瑞康实业有限公司 | 一种生成心搏间期差值数据序列的方法和装置 |
CN113712564A (zh) * | 2020-05-12 | 2021-11-30 | 深圳市科瑞康实业有限公司 | 一种心电信号分类设备和方法 |
CN114916923A (zh) * | 2022-05-19 | 2022-08-19 | 南京中医药大学 | 一种电机械互联的心电脉搏信号分析方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1188635A (zh) * | 1997-01-22 | 1998-07-29 | 张才德 | 心功能电脑检测诊断仪 |
CN2363654Y (zh) * | 1999-01-19 | 2000-02-16 | 华南理工大学 | 心功能参数检测仪 |
CN100413464C (zh) * | 2006-05-26 | 2008-08-27 | 中国人民解放军空军航空医学研究所 | 在脉搏波法动脉血压连续测量中的脉搏波传导时间的获取方法和装置 |
-
2007
- 2007-12-06 CN CNB2007101151437A patent/CN100515327C/zh not_active Expired - Fee Related
Cited By (43)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8771197B2 (en) | 2009-02-11 | 2014-07-08 | Edwards Lifesciences Corporation | Detection of parameters in cardiac output related waveforms |
CN105243273A (zh) * | 2009-02-11 | 2016-01-13 | 爱德华兹生命科学公司 | 心输出量相关波形的参数检测 |
CN102387742A (zh) * | 2009-02-11 | 2012-03-21 | 爱德华兹生命科学公司 | 心输出量相关波形的参数检测 |
CN102387742B (zh) * | 2009-02-11 | 2015-11-25 | 爱德华兹生命科学公司 | 心输出量相关波形的参数检测 |
CN105243273B (zh) * | 2009-02-11 | 2018-03-27 | 爱德华兹生命科学公司 | 心输出量相关波形的参数检测 |
CN101862201A (zh) * | 2010-06-23 | 2010-10-20 | 南京邮电大学 | 手机智能心音看诊器中心音图的表示方法 |
CN101862201B (zh) * | 2010-06-23 | 2011-08-31 | 南京邮电大学 | 手机智能心音看诊器中心音图的表示方法 |
CN102334986A (zh) * | 2011-07-08 | 2012-02-01 | 重庆大学 | 动脉压信号中重搏切迹点识别方法 |
CN102334986B (zh) * | 2011-07-08 | 2013-09-25 | 重庆大学 | 动脉压信号中重搏切迹点识别方法 |
CN102551690B (zh) * | 2011-12-29 | 2013-10-09 | 乐普(北京)医疗器械股份有限公司 | 人体信号自适应分析方法 |
CN102551690A (zh) * | 2011-12-29 | 2012-07-11 | 乐普(北京)医疗器械股份有限公司 | 人体信号自适应分析方法 |
CN102727245A (zh) * | 2012-07-16 | 2012-10-17 | 山东大学 | 一种心音信号采集调理传输装置 |
CN103006207A (zh) * | 2013-01-11 | 2013-04-03 | 山东师范大学 | 基于形态学的t波交替散点图法的心电信号分析方法 |
CN103006207B (zh) * | 2013-01-11 | 2016-03-02 | 山东师范大学 | 基于形态学的t波交替散点图法的心电信号分析方法 |
CN104414679A (zh) * | 2013-09-04 | 2015-03-18 | 北京航空航天大学 | 一种中心动脉压重建中外周压力波延迟时间的确定方法 |
CN104665798A (zh) * | 2015-01-26 | 2015-06-03 | 周常安 | 血压管理装置及方法 |
CN105078424A (zh) * | 2015-08-26 | 2015-11-25 | 深圳诺康医疗设备有限公司 | 信号采集方法、脉搏监测和主监测装置及综合监测系统 |
CN105310724A (zh) * | 2015-09-28 | 2016-02-10 | 何宗彦 | 一种测量脉搏波传播速度的方法及系统 |
CN105496402A (zh) * | 2015-11-20 | 2016-04-20 | 北京理工大学 | 基于散点图和符号动力学的心电特征分析方法 |
CN105496402B (zh) * | 2015-11-20 | 2018-03-02 | 北京理工大学 | 基于散点图和符号动力学的心电特征分析方法 |
CN107126203A (zh) * | 2016-02-28 | 2017-09-05 | 南京全时陪伴健康科技有限公司 | 人体健康信息实时测量显示系统 |
CN107174197A (zh) * | 2016-03-10 | 2017-09-19 | 深圳市理邦精密仪器股份有限公司 | 监护设备及其波形显示方法、装置 |
CN107361753A (zh) * | 2017-08-29 | 2017-11-21 | 哈尔滨理工大学 | 基于脉搏波形特征点的人体健康状态监护方法 |
CN108185996A (zh) * | 2017-12-27 | 2018-06-22 | 中国科学院深圳先进技术研究院 | 动脉血管年龄估算模型构建方法和装置 |
CN108577883A (zh) * | 2018-04-03 | 2018-09-28 | 上海交通大学 | 一种冠心病筛查装置、筛查系统以及信号特征提取方法 |
CN109222964B (zh) * | 2018-07-20 | 2021-02-09 | 广州视源电子科技股份有限公司 | 房颤检测装置及存储介质 |
CN109222964A (zh) * | 2018-07-20 | 2019-01-18 | 广州视源电子科技股份有限公司 | 房颤检测装置及存储介质 |
CN109044299A (zh) * | 2018-09-25 | 2018-12-21 | 北方民族大学 | 一种心血管健康检测的生理信号采集设备及方法 |
CN109044299B (zh) * | 2018-09-25 | 2021-03-30 | 北方民族大学 | 一种心血管健康检测的生理信号采集设备及方法 |
CN110200617A (zh) * | 2019-06-26 | 2019-09-06 | 慧言科技(天津)有限公司 | 一种利用心率变异性定量评价疼痛的系统 |
CN110490176A (zh) * | 2019-09-03 | 2019-11-22 | 南京邮电大学 | 一种基于多时窗心音信号的多元特征提取方法 |
CN110490176B (zh) * | 2019-09-03 | 2022-08-26 | 南京邮电大学 | 一种基于多时窗心音信号的多元特征提取方法 |
CN110710960A (zh) * | 2019-09-12 | 2020-01-21 | 东莞市康助医疗科技有限公司 | 一种动脉血管年龄推算方法及装置 |
CN110956406A (zh) * | 2019-12-07 | 2020-04-03 | 中国科学院心理研究所 | 一种基于心率变异性的团队协同能力的评估模型 |
CN110956406B (zh) * | 2019-12-07 | 2022-11-18 | 中国科学院心理研究所 | 一种基于心率变异性的团队协同能力的评估方法 |
CN111543978A (zh) * | 2020-05-11 | 2020-08-18 | 中北大学 | 基于心电、心音和脉搏同步诊断的方法及系统 |
CN113712566A (zh) * | 2020-05-12 | 2021-11-30 | 深圳市科瑞康实业有限公司 | 一种生成心搏间期差值数据序列的方法和装置 |
CN113712564A (zh) * | 2020-05-12 | 2021-11-30 | 深圳市科瑞康实业有限公司 | 一种心电信号分类设备和方法 |
CN113712564B (zh) * | 2020-05-12 | 2023-09-01 | 深圳市科瑞康实业有限公司 | 一种心电信号分类设备和方法 |
CN113712566B (zh) * | 2020-05-12 | 2024-02-06 | 深圳市科瑞康实业有限公司 | 一种生成心搏间期差值数据序列的方法和装置 |
CN112842350A (zh) * | 2021-02-20 | 2021-05-28 | 无锡市中健科仪有限公司 | 基于正反传播算法获得心拍的主导心率的方法 |
CN113598724A (zh) * | 2021-09-03 | 2021-11-05 | 上海市高血压研究所 | 用于评估心血管自主神经功能的装置 |
CN114916923A (zh) * | 2022-05-19 | 2022-08-19 | 南京中医药大学 | 一种电机械互联的心电脉搏信号分析方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN100515327C (zh) | 2009-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100515327C (zh) | 一种心血管系统联合变异性指标的检测方法和装置 | |
CN101732040B (zh) | 一种无创多路脉搏波检测方法 | |
CN102302361B (zh) | 一种心率及身体状态监测装置及方法 | |
Stein et al. | Association of Holter-based measures including T-wave alternans with risk of sudden cardiac death in the community-dwelling elderly: the Cardiovascular Health Study | |
CN103690156B (zh) | 一种心率获取方法及心电信号的处理方法 | |
CN101427917B (zh) | 基于固有趋势子序列模式分解的出现新的心脏活动趋势的ecg信号获取方法 | |
Altan et al. | A new approach to early diagnosis of congestive heart failure disease by using Hilbert–Huang transform | |
KR102008196B1 (ko) | 심전도 데이터를 이용한 혈중 칼륨농도 예측모델 생성장치 및 그 방법 | |
US20100113893A1 (en) | Method for measuring physiological stress | |
CN101357062A (zh) | 基于容积脉搏波检测的生命体征参数测评装置 | |
CN109691994A (zh) | 一种基于心电图的心率监测分析方法 | |
CN102462494A (zh) | 一种新型智能心电测试保健器 | |
Mai et al. | Non-contact heartbeat detection based on ballistocardiogram using UNet and bidirectional long short-term memory | |
Bhoi et al. | QRS Complex Detection and Analysis of Cardiovascular Abnormalities: A Review. | |
CN100342820C (zh) | 一种心率变异性的预测度指标检测分析方法及仪器 | |
Pang et al. | Real time heart ischemia detection in the smart home care system | |
CN1422591A (zh) | 能由颈部同时测量心电、脉搏和声波信号的传感器 | |
CN104027097A (zh) | 血管功能无创检测方法及装置 | |
Amanipour et al. | The effects of blood glucose changes on frequency-domain measures of HRV signal in type 1 diabetes | |
Palacios et al. | Changes of autonomic information flow due to idiopathic dilated cardiomyopathy | |
CN102940484B (zh) | 可远程通讯的家用心血管检测和诊断装置 | |
US8165665B2 (en) | Chip for sensing a physiological signal and method for sensing the same | |
Chen et al. | Beat-to-beat heart rate detection based on seismocardiogram using BiLSTM network | |
KR102574331B1 (ko) | 멀티스케일 누적 잔여 분산 엔트로피 기반 향상된 심박변이율 정량화 방법 | |
CN114886403A (zh) | 一种基于脉搏主波间期的恶性心律失常识别与预测系统 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090722 Termination date: 20171206 |