技术简介:
本专利针对传统血压测量方法存在侵入性、准确性不足的问题,提出一种基于高灵敏传感器的无创连续监测方案。通过采集脉搏波并进行滤波处理,提取主波幅度、降中峡幅度等7项关键特征参数,结合线性回归模型建立血压计算公式,实现无需袖带的实时血压监测,显著提升测量便捷性与舒适度。
关键词:无创血压监测,脉搏波分析,线性回归模型
一种无创连续血压的测量方法及装置制造方法
【专利摘要】本发明适用于血压测量【技术领域】,提供了一种无创连续血压的测量方法及装置,包括:获取传感器采集的脉搏数据,所述脉搏数据包括脉搏波;对获取的脉搏波进行滤波;根据滤波后的脉搏数据提取脉搏波的特征参数,所述特征参数包含主波幅度h1、主动脉扩张降压点h2、降中峡幅度h3、重搏波幅度h4、左心室的快速射血期t1、左心室的收缩期t2、左心室舒张期t3与脉动周期t;将脉搏波的特征参数代入线性回归方程:Y=b0+b1x1+b2x2+b3x3+……+b7x7,计算出相应的偏回归系数;根据脉搏波的数据与计算出的偏回归系数,计算出用户的血压值。增加了血压测量的便捷性,从而提高了测量血压的舒适度与测量时间。
【专利说明】
【技术领域】
[0001] 本发明属于血压测量【技术领域】,尤其涉及一种无创连续血压的测量方法及装置。 一种无创连续血压的测量方法及装置
【背景技术】
[0002] 血压是指血液在血管中流动时对血管壁所施加的侧压力,它是反映心血管功能的 重要生理指标,在疾病诊断、治疗效果观察、进行预后判断等都有着重要的意义。心室收缩, 血液从心室流入动脉,此时血液对动脉的压力最高,称为收缩压。心室舒张,动脉血管弹性 回缩,血液仍慢慢继续向前流动,但血压下降,此时的压力称为舒张压。
[0003] 目前,测量血压的方法大致分为两类,一种称为直接测量法,另一种称为间接测量 法。
[0004] 直接测量法是将连接压力传感器的导管经皮直接插入大动脉或心脏检查血压信 号,能够进行连续测量。其优点是数据准确,缺点是会在使用者的身体上留有一定的创伤, 不但会对患者造成痛苦,还可能引起感染;而且,介入会引起机体的应激反应,身体处于非 自然状态,一定程度上影响血压值。
[0005] 间接测量法通过检测动脉管壁的搏动,血管容积变化等参数间接得到血压,这种 方法优点简便易行,缺点是在临床上应用广泛;缺点:只能测出人体某一时刻的血压值,且 测量的血压往往不在同一心动周期内,测量的精度容易受到环境的影响。
[0006] 然而,现有的无创动态血压采用的是示波法与容积补偿法,产品化的无创动态血 压测量仪大多使用这两种方法设计而成,在血压测量过程中始终有气囊对人体的束缚,而 且设备及测量过程较为复杂,不便于长时间的操作和使用。
【发明内容】
[0007] 本发明实施例的目的在于提供一种无创连续血压测量方法,旨在解决不需要阻断 人体动脉血管内血液的流动,只需使用接触式传感器轻微压迫人体腕部桡动脉,即可测得 人体腕部动脉血压值,并将人体腕部动脉压力数值经过计算,换算出人体动脉血压的问题。
[0008] 本发明的技术方案这样实现的,一种无创连续血压测量方法,包括:
[0009] 步骤1 :获取传感器采集的脉搏数据,所述脉搏数据包括脉搏波;
[0010] 步骤2 :对获取的脉搏波进行滤波;
[0011] 步骤3 :根据滤波后的脉搏数据提取脉搏波的特征参数,所述特征参数包含主波 幅度4、主动脉扩张降压点h2、降中峡幅度h 3、重搏波幅度h4、左心室的快速射血期&、左心 室的收缩期t2、左心室舒张期t 3与脉动周期t ;
[0012] 步骤4 :将脉搏波的特征参数代入线性回归方程,计算出相应的偏回归系数;
[0013] Y = b0+b1x1+b2x2+b 3x3+b4x4+b5x5+b 6x6+b7x7,
[0014] 式中Y为收缩压或舒张压,k为常数项,bp b2, b3,, b4, b5, b6, b7为偏回归系数, x2, x3, x4, x5, x6, x7为被选入的脉搏特征参数;
[0015] 步骤5 :根据脉搏波的数据与计算出的偏回归系数,计算出用户的血压值。
[0016] 进一步,所述步骤1中获取所述脉搏数据的具体步骤:
[0017] 所述脉搏数据还包括粗略血压值,通过血压计得到用户相应粗略血压值,使用脉 搏传感器提取用户桡动脉的脉搏波。
[0018] 进一步,所述步骤2中对获取的脉搏波进行滤波方式为采用巴特斯滤波滤去高次 谐波。
[0019] 进一步,所述步骤4中计算偏回归系数具体包括:
[0020] 脉搏特征参数 Xp x2, x3, x4, x5, x6, x7 分别对应的为 hi, lvV lvV W V,Z, beat, 其中,V为上升沿平均斜率h/tp Z为h# (l+t2/t3), beat为脉率1/(n/fs) *60,其中,n为单 位周期采样点个数,fs为采样频率。
[0021] 进一步,根据同一完整周期T内的所述脉搏波中的最大值c_aBP_max与最小值c_ aBP_min,通过满足公式:(A^-Ai) * (Ai-Ai+1)〈0
[0022] 其中Ai是第i点的幅值的点,即是极值点,根据脉搏波中的最大幅值与极值点幅 值之差设置一个阈值w,所述阈值w = Τ/2 π,若在极值点满足阈值w条件,则找出主波幅度 h ;根据主波幅度h左边的最小极值点找出脉搏周期的起始点B,将主波峰左边距离最近 的极小值点确认为脉搏波周期的起始点B ;分别设置横向阈值Wl、w2,若有极值点与它前面 的主波峰时间间隔满足小于阈值A,则确认该极值点是降中峡E,其幅度值为h 3 ;若有极值 点与它前面的主波峰时间间隔满足小于阈值《2,则确认该极值点是重搏波点F,其幅度值为 h 4 ;在主波峰h和降中峡h3之间找到二阶差分最大值点D,其幅度值为h2。
[0023] 本发明的另一目的在于提供一种无创连续血压测量装置,包括:
[0024] 采集单元,用于获取传感器采集的脉搏数据,所述脉搏数据为脉搏波;
[0025] 滤波单元,用于对获取的脉搏波进行滤波;
[0026] 提取单元,用于对滤波后的脉搏数据提取脉搏波的特征参数,所述特征参数包含 主波幅度4、主动脉扩张降压点h 2、降中峡幅度h3、重搏波幅度h4、左心室的快速射血期心、 左心室的收缩期t 2、左心室舒张期t3与脉动周期t ;
[0027] 模型单元,将脉搏波的特征参数代入线性回归方程:Y = x5+b6x6+b7x7,计算出相应的偏回归系数;
[0028] 计算单元,用于根据脉搏波的数据与计算出的偏回归系数,套入模型单元计算出 用户的血压值。
[0029] 进一步,所述采集单元包含通过血压计采集用户粗略的血压值与使用脉搏传感器 提取用户桡动脉的脉搏波。
[0030] 进一步,所述滤波单元采用巴特斯滤波方式滤去脉搏波的高次谐波。
[0031] 进一步,所述模型单元包括偏回归系数子单元,用于将脉搏特征参数Xl,χ 2, χ3, χ4, x5, x6, x7分别对应的为hp hyhD lv^,h/hp V,Z, beat,套入线性回归方程计算出偏回归系 数bQ至b7。
[0032] 进一步,所述提取单元包括:
[0033] 检测子单元,用于检测在同一周期内对满足公式所有的极值点;
[0034] 第一定位子单元,用于根据脉搏波中的最大值与最小值,查询对应极值点定位起 始点B与主波峰C,幅度值h ;
[0035] 第二定位子单元,用于根据主波峰C设置横向阈值Wl、w2,若有极值点与主波峰C 的时间间隔小于阈值Wi,查找出降中峡E,幅度值h3 ;若有极值点与主波峰C的时间间隔小 于阈值w2,查找出重搏波点F,幅度值h4 ;
[0036] 第三定位子单元,用于根据主波峰h和降中峡1!3之间找到二阶差分最大值点即降 压值点D,幅度值h 2。
[0037] 在本发明中,通过灵敏度很高的压力传感器获取到用户的脉搏波,对所述脉搏波 进行滤波,并提取脉搏波的特征参数,将特征参数套入公式Y中,计算出用户的血压值;这 种血压测量方法以使用户在无创且连续的状态下测量血压,同时避免了无法准确血压的间 接测量法,增加了血压测量的便捷性,从而提高了测量血压的舒适度与测量时间。
【专利附图】
【附图说明】
[0038] 图1是本发明实施例提供的一种无创连续血压测量方法的实现流程图;
[0039] 图2是本发明实施例提供的文一种无创连续血压测量滤波前的图;
[0040] 图3是本发明实施例提供的一种无创连续血压测量滤波后的图;
[0041] 图4是本发明实施例提供的脉搏波形图;
[0042] 图5是本发明实施例提供的脉搏波特征参数列表;
[0043] 图6是本发明实施例提供的一种无创连续血压测量装置的结构框图。
【具体实施方式】
[0044] 为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对 本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并 不用于限定本发明。
[0045] 图1是本发明实施例提供的一种无创连续血压测量方法的实现流程图,详述如 下:
[0046] 步骤1 :获取传感器采集的脉搏数据,所述脉搏数据包括脉搏波;
[0047] 在步骤1中的实施例中,首先利用欧姆龙电子血压计得到相应的粗略的血压值, 然后通过脉搏波传感器HK2000C提取桡动脉的脉搏波,传感器的采样频率为200Hz,其中传 感器每秒钟1秒收到200个数据,每次采集800个数据进行处理,大约4到5个脉搏波周期, 采集到脉搏波。
[0048] 步骤2 :对获取的脉搏波进行滤波;
[0049] 在步骤2中的实施例中,由于人体生物电信号具有高阻抗、信号弱频率低的特点, 并且处于较大的噪声背景之中,基于这些特性可知原始波形中包含有许多的噪声(即高次 谐波),首先进行巴特沃斯滤波(一种常见的低通滤波)滤去高次噪声,滤去高次谐波。
[0050] 步骤3 :根据滤波后的脉搏数据提取脉搏波的特征参数,所述特征参数包含主波 幅度4、主动脉扩张降压点h2、降中峡幅度h 3、重搏波幅度h4、左心室的快速射血期&、左心 室的收缩期t2、左心室舒张期t 3与脉动周期t ;
[0051] 在步骤3中的实例中,通过对脉搏波的中波形图分析,结合公式(2)中的极值点分 析对比,依次定位出满足极值点的B、C、E、F、D对应找出相应的特征参数。
[0052] 步骤4 :将脉搏波的特征参数代入线性回归方程,计算出相应的偏回归系数;
[0053] Y = b0+b1x1+b2x2+b 3x3+b4x4+b5x5+b 6x6+b7x7,
[0054] 式中Y为收缩压或舒张压,K为常数项,bp b2, b3,, b4, b5, b6, b7为偏回归系数,Xp x2, x3, x4, x5, x6, x7为被选入的脉搏特征参数;
[0055] 步骤5 :根据脉搏波的数据与计算出的偏回归系数,计算出用户的血压值。
[0056] 优选的,所述步骤1中获取所述脉搏数据的具体步骤:
[0057] 通过血压计得到用户相应粗略血压值,使用脉搏传感器提取用户桡动脉的脉搏 波。
[0058] 如图2至图3所示,为脉搏波滤波前和滤波后的图。
[0059] 优选的,所述步骤2中对获取的脉搏波进行滤波方式为采用巴特斯滤波滤去高次 谐波。
[0060] 如图5所示,为本发明实施例提供的脉搏波特征参数列表。
[0061] 在本实例中,根据滤波后脉搏波的波形特点选取有效的特征参数:
[0062] h :主波幅度,为主波峰到脉搏波图基线的高度。
[0063] h2:主动脉扩张降压点,为主波峰到降中峡的二阶差分最大值。
[0064] h3 :降中峡幅度,为降中峡谷底到脉搏波图基线的高度。
[0065] h4:重搏波幅度,为重搏波峰到脉搏波图基线的高度。
[0066] t1:左心室的快速射血期,即脉搏波形起始点到主波峰的时间间隔。
[0067] t2:左心室的收缩期,即脉搏波形起始点到降中峡的时间间隔。
[0068] t3:左心室舒张期,即降中峡到脉搏波形终止点(下一周期脉搏波的起始点)之间 的时间间隔。
[0069] T :脉动周期,即脉搏波形起点到终点之间的时间间隔。
[0070] 优选的,所述步骤4中计算偏回归系数具体包括:
[0071] 脉搏特征参数 XpXrXrXpXyXvXT 分别对应的为 hi, lvV lvV W V,Z, beat, 其中,V为上升沿平均斜率h/tp Z为h# (l+t2/t3), beat为脉率1/(n/fs) *60,其中,n为单 位周期采样点个数,fs为采样频率。
[0072] 如图4所示,为本发明实施例提供的脉搏波形图。
[0073] 优选的,根据同一完整周期T内的所述脉搏波中的最大值c_aBP_max与最小值c_ aBP_min,通过满足公式:(A^-Ai) * (Ai-Ai+1)〈0
[0074] 其中Ai是第i点的幅值的点,即是极值点,根据脉搏波中的最大幅值与极值点幅 值之差设置一个阈值w,所述阈值w = Τ/2 π,若在极值点满足阈值w条件,则找出主波幅度 h ;根据主波幅度h左边的最小极值点找出脉搏周期的起始点B,将主波峰左边距离最近 的极小值点确认为脉搏波周期的起始点B ;分别设置横向阈值Wl、w2,若有极值点与它前面 的主波峰时间间隔满足小于阈值A,则确认该极值点是降中峡E,其幅度值为h 3 ;若有极值 点与它前面的主波峰时间间隔满足小于阈值《2,则确认该极值点是重搏波点F,其幅度值为 h 4 ;在主波峰h和降中峡h3之间找到二阶差分最大值点D,其幅度值为h2。
[0075] 通过上述步骤将脉搏波的特征参数以及各种偏回归系数代入线性回归方程:Y = b0+b1x1+b2x2+b 3x3+......+b7x7,算出用户在一个周其月内连续状态的血压值。为了减少个体差异 性,提升测量血压精确性,因为血压测量设备会面临不同的人使用,刚刚建立起来的辨识参 数就需要更新,来消除个体差异引起的血压测量误差,因此这就涉及到校正部分。该部分主 要是为了解决不同的人,或者当设备电量不足时所引起的误差问题。所以算法里必须要体 现不断的更新辨识数据库来重新获取新的辨识参数。
[0076] 图6所示,为本发明实施例提供的一种无创连续血压测量装置的结构框图。
[0077] 采集单元1,用于获取传感器采集的脉搏数据,所述脉搏数据为脉搏波;
[0078] 滤波单元2,用于对获取的脉搏波进行滤波;
[0079] 提取单元3,用于对滤波后的脉搏数据提取脉搏波的特征参数,所述特征参数包含 主波幅度4、主动脉扩张降压点h2、降中峡幅度h3、重搏波幅度h4、左心室的快速射血期心、 左心室的收缩期t 2、左心室舒张期t3与脉动周期t ;
[0080] 模型单元4,将脉搏波的特征参数代入线性回归方程:Y = bd+bA+bj^+t^xflv^+b 5x5+b6x 6+b7x7,计算出相应的偏回归系数;
[0081] 计算单元5,用于根据脉搏波的数据与计算出的偏回归系数,套入模型单元4计算 出用户的血压值。
[0082] 进一步,所述采集单元1包含通过血压计采集用户粗略的血压值与使用脉搏传感 器提取用户桡动脉的脉搏波。
[0083] 进一步,所述滤波单元2采用巴特斯滤波方式滤去脉搏波的高次谐波。
[0084] 进一步,所述模型单元4包括偏回归系数子单元41,用于将脉搏特征参数Xl,x 2, X3,......χ7分别对应的为hp hyhD lv^,h/hp V,Z, beat,套入公式(1)计算出偏回归系数 b0 至 b7。
[0085] 进一步,所述提取单元3包括:
[0086] 检测子单元31,用于检测在同一周期内对满足公式(2)所有的极值点;
[0087] 第一定位子单元32,用于根据脉搏波中的最大值与最小值,查询对应极值点定位 起始点B与主波峰C,幅度值h ;
[0088] 第二定位子单元33,用于根据主波峰C设置横向阈值Wl、w2,若有极值点与主波峰 C的时间间隔小于阈值 Wl,查找出降中峡E,幅度值h3;若有极值点与主波峰C的时间间隔小 于阈值w2,查找出重搏波点F,幅度值h 4 ;
[0089] 第三定位子单元34,用于根据主波峰h和降中峡h3之间找到二阶差分最大值点 即降压值点D,幅度值h 2。
[0090] 在本发明实施例中,通过灵敏度很高的压力传感器获取到用户的脉搏波,对所述 脉搏波进行滤波,并提取脉搏波的特征参数,将特征参数套入公式Y中,计算出用户的血压 值;这种血压测量方法以使用户在无创且连续的状态下测量血压,同时避免了无法准确血 压的间接测量法,增加了血压测量的便捷性,从而提高了测量血压的舒适度与测量时间。
[0091] 以上所述,仅为本发明的【具体实施方式】,但本发明的保护范围并不局限于此,任何 熟悉本【技术领域】的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应 涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
【权利要求】
1. 一种无创连续血压的测量方法,其特征在于,包括:步骤1 :获取传感器采集的脉搏数据,所述脉搏数据包括脉搏波;步骤2 :对获取的脉搏波进行滤波;步骤3 :根据滤波后的脉搏数据提取脉搏波的特征参数,所述特征参数包含主波幅度 4、主动脉扩张降压点h2、降中峡幅度h3、重搏波幅度h4、左心室的快速射血期&、左心室的 收缩期t 2、左心室舒张期t3与脉动周期t ;步骤4 :将脉搏波的特征参数代入线性回归方程,计算出相应的偏回归系数;Y = b0+b1x1+b2x2+b 3x3+b4x4+b5x5+b 6x6+b7x7,式中Y为收缩压或舒张压,k为常数项,bp b2, b3,, b4, b5, b6, b7为偏回归系数,Χρ x2, x3, x4, x5, x6, x7为被选入的脉搏特征参数;步骤5 :根据脉搏波的数据与计算出的偏回归系数,计算出用户的血压值。
2. 根据权利要求1所述的一种无创连续血压的测量方法,其特征在于,所述步骤1中获 取所述脉搏数据的具体步骤: 所述脉搏数据还包括粗略血压值,通过血压计得到用户相应粗略血压值,再使用脉搏 传感器提取用户桡动脉的脉搏波。
3. 根据权利要求1所述的一种无创连续血压的测量方法,其特征在于,所述步骤2中对 获取的脉搏波进行滤波方式为采用巴特斯滤波滤去高次谐波。
4. 根据权利要求1所述的一种无创连续血压的测量方法,其特征在于,所述步骤4中计 算偏回归系数具体包括: 脉搏特征参数Xp X2, x3, x4, x5, xe,分别对应的为lvV lvV W V,Z, beat,其 中,V为上升沿平均斜率h/tp Z为hfQ+iv%) ,beat为脉率l/(n/fs)*60,其中,n为单位 周期采样点个数,fs为采样频率。
5. 根据权利要求1或4所述的一种无创连续血压的测量方法,其特征在于,根据同一 完整周期T内的所述脉搏波中的最大值c_aBP_max与最小值c_aBP_min,通过满足公式: (A i_1-Ai)*(Ai-Ai+1)<0 其中Ai是第i点的幅值的点,即是极值点,根据脉搏波中的最大幅值与极值点幅值之 差设置一个阈值w,所述阈值w = Τ/2 π,若在极值点满足阈值w条件,则找出主波幅度h ; 根据主波幅度h左边的最小极值点找出脉搏周期的起始点B,将主波峰左边距离最近的极 小值点确认为脉搏波周期的起始点B ;分别设置横向阈值Wl、w2,若有极值点与它前面的主 波峰时间间隔满足小于阈值A,则确认该极值点是降中峡E,其幅度值为h 3 ;若有极值点与 它前面的主波峰时间间隔满足小于阈值w2,则确认该极值点是重搏波点F,其幅度值为h 4 ; 在主波峰h和降中峡h3之间找到二阶差分最大值点D,其幅度值为h2。
6. -种无创连续血压的测量装置,其特征在于,包括: 采集单元,用于获取传感器采集的脉搏数据,所述脉搏数据包括脉搏波;滤波单元,用于对获取的脉搏波进行滤波;提取单元,用于对滤波后的脉搏数据提取脉搏波的特征参数,所述特征参数包含主波 幅度4、主动脉扩张降压点h2、降中峡幅度h3、重搏波幅度h4、左心室的快速射血期&、左心 室的收缩期t 2、左心室舒张期t3与脉动周期t ;模型单元,将脉搏波的特征参数代入线性回归方程:Y = bd+bA+bj^+t^xflv^+bgXg+bj;x6+b7x7,计算出相应的偏回归系数;计算单元,用于根据脉搏波的数据与计算出的偏回归系数,套入模型单元计算出用户 的血压值。
7. 根据权利要求6所述的一种无创连续血压的测量装置,其特征在于,所述采集单元 包含通过血压计采集用户粗略的血压值与使用脉搏传感器提取用户桡动脉的脉搏波。
8. 根据权利要求6所述的一种无创连续血压的测量装置,其特征在于,所述滤波单元 采用巴特斯滤波方式滤去脉搏波的高次谐波。
9. 根据权利要求6所述的一种无创连续血压的测量装置,其特征在于,所述模型单元 包括偏回归系数子单元,用于将脉搏特征参数Xp X2, x3, Xy x;5, Xe,分别对应的为hph;/ hp h/hp V,Z, beat,套入线性回归方程计算出偏回归系数至b7。
10. 根据权利要求6或9所述的一种无创连续血压的测量装置,其特征在于,所述提取 单元包括: 检测子单元,用于检测在同一周期内对满足公式(Ag-AiMAi-AwXO的所有极值点;第一定位子单元,用于根据脉搏波中的最大值与最小值,查询对应极值点定位起始点B 与主波峰C,其幅度值为h ;第二定位子单元,用于根据主波峰C设置横向阈值Wl、w2,若有极值点与主波峰C的时 间间隔小于阈值^,查找出降中峡E,其幅度值为h3 ;若有极值点与主波峰C的时间间隔小 于阈值w2,查找出重搏波点F,其幅度值为h4 ;第三定位子单元,用于根据主波峰4和降中峡1!3之间找到二阶差分最大值点即降压值 点D,其幅度值为h2。
【文档编号】A61B5/021GK104116503SQ201410339206
【公开日】2014年10月29日 申请日期:2014年7月16日 优先权日:2014年7月16日
【发明者】陈智华, 潘林强, 何龙, 顾超, 梅晰山, 程昆仑 申请人:华中科技大学