CN113171067B - 一种基于pwv的血压测量手表 - Google Patents
一种基于pwv的血压测量手表 Download PDFInfo
- Publication number
- CN113171067B CN113171067B CN202110246028.3A CN202110246028A CN113171067B CN 113171067 B CN113171067 B CN 113171067B CN 202110246028 A CN202110246028 A CN 202110246028A CN 113171067 B CN113171067 B CN 113171067B
- Authority
- CN
- China
- Prior art keywords
- pwtt
- pwv
- values
- value
- blood pressure
- 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.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/021—Measuring pressure in heart or blood vessels
- A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
- A61B5/02125—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics of pulse wave propagation time
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Vascular Medicine (AREA)
- Biomedical Technology (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明提供了一种基于PWV的血压测量手表,属于无创连续血压测量技术领域。本发明的PWV测量手表进行PWV的计算,包括采集第一次筛选所需要的指定数据集的初筛,丢弃初始采集到的误差大的数据,同时丢弃后续不在当前总体平均值置信区间内的数据;当上一步初筛后的数据未达到规定的有效值时进行多次循环数据采集,并重新计算当前总体平均值置信区间,直到采集到规定的数据集;最后对采集到的规定数据集进行PWV值的计算,提高了血压测量的精度。
Description
技术领域
本发明涉及无创连续血压测量技术领域,尤其涉及一种基于PWV的血压测量手表。
背景技术
脉搏波传导时间(Pulse Wave Transit Time,PWTT)是动脉血压波沿血管壁传输时间,这个时间是由从心脏开始收缩到大动脉瓣打开再到出现在大动脉血液脉动,经过血管传导到末梢部位时间组成,PWTT参数可以从手腕部的心电和脉搏波来进行获取,目前有大量临床试验可知脉搏波传导时间与血压之间变化关系是生理学上已经明确现象,利用该现象,结合个人性别、身高、体重、臂长和年龄等基础特征,形成特殊算法可以用来监测动脉血压的变化及在一定误差范围内的血压值。脉搏波波速(Pulse Wave Velocity,PWV),是脉搏波在动脉系统的两个既定点间的传播速度。
目前,PWV的测量方法不能根据不同场景下采集到的实际PWTT数据,在指定误差范围和有效时间内计算出PWV值进行灵活动态调整,PWTT的数据筛选通常仅去除前期不稳定的数据,对最终的PWV计算精度影响较大。
因此,有必要提出一种灵活动态的PWV的测量方法,对采集的PWTT进行灵活筛选,从而为不同场景下采集到的实际PWTT数据在有效时间内计算出指定误差范围内的PWV值。
在中国专利申请文献CNCN06618537A中,公开了一种基于脉搏波传导的连续动态血压精确监测法,包括(1)获取被测者个人参数;(2)通过脉搏波传感器和心电信号采集模块连续同步采集心电信号ECG、脉搏波信号s1并发送至控制模块;(3)得到心-桡脉搏波传导时间hrPWTT以及舒张期时长TD;(4)根据心-桡脉搏波传导时间hrPWTT计算心-桡脉搏波传导速度hrPWV;(5)根据心-桡脉搏波传导速度hrPWV计算收缩压PWVSBP;(6)根据收缩压PWVSBP以及舒张期时长TD计算舒张压PWVDBP。PWVDBP采用如下公式计算得到:
现有技术至少存在以下不足:
1.现有技术中在采集PWTT时,通常仅丢掉最初波动期的数据,或者不丢弃,造成PWTT值不准确。
2.现有技术中技术PWV较复杂,通过心-桡脉搏波传导时间,舒张期时长,计算过程未考虑臂长的因素,降低了PWV计算结果精度。
发明内容
本发明的目的在于克服现有算法不能在计算时间和计算误差上灵活动态调整的缺陷,提供基于PWV的血压测量手表,其中对PWTT数据进行灵活筛选,且通过PWTT、臂长与正常人肩膀到心脏的平均距离得到PWV,使得最终计算的PWV精度更高。本发明的方法包括采集第一次筛选所需要的指定数据集的初筛,丢弃初始采集到的误差大的数据,同时丢弃后续不在当前总体平均值置信区间内的数据;当上一步初筛后的数据未达到规定的有效值时进行多次循环数据采集,并重新计算当前总体平均值置信区间,直到采集到规定的数据集;最后对采集到的规定数据集进行PWV值的计算。
本发明提供了一种基于PWV的血压测量手表,所述血压测量手表进行所述PWV的计算,包括如下步骤:
PWTT初步采集及初筛步骤;
S1000:采集N个PWTT值,丢弃前面P个PWTT值;
S2000:计算剩余(N-P)个PWTT值的均值,根据置信度M%计算第一置信区间;其中,5%≤M%≤10%;
S3000:丢弃剩余(N-P)个PWTT值中不在第一置信区间内的PWTT值,如果剩余的PWTT个数满足预设的最低PWTT个数要求,则执行PWV计算步骤,如果剩余的PWTT个数不满足预设的最低PWTT个数要求,则执行PWTT继续采集及筛选步骤;
PWTT继续采集及筛选步骤;
S4000:继续采集m个PWTT值,使得PWTT总个数达到预设的最低PWTT个数;
S5000:计算PWTT值的均值,根据置信度M%计算第二置信区间;每次迭代计算的所述第二置信区间与当前参与置信区间计算的PWTT值有关;
S6000:丢弃PWTT值不在第二置信区间内的PWTT值;
S7000:如果剩余的PWTT值个数满足预设的最低PWTT个数要求,则执行PWV计算步骤,否则,继续执行步骤S4000;
PWV计算步骤;
S8000:计算保留的PWTT值的均值S;
S9000:计算PWV值;其中PWV值采用如下公式进行计算;
其中,
PWV为最终计算得到的脉搏波波速;
L为待测者的臂长;
a为正常人肩膀到心脏的平均距离;
S为最终保留的PWTT值的均值。
优选地,所述初步采集的PWTT值的个数N满足,5<N≤15。
优选地,所述初步采集的PWTT值的个数N=15。
优选地,所述初步筛选时丢弃的PWTT值的个数P满足,5≤P≤9。
优选地,所述初步筛选时丢弃的PWTT值的个数P=5。
优选地,所述预设的最低PWTT个数为N-P。
优选地,置信区间(A1,A2)的计算采用如下公式:
A1=A-A*M%;
A2=A+A*M%;
其中,
A为当前保留的PWTT值的均值;
M%为置信度;
A1为置信区间下限;
A2为置信区间上限。
与现有技术相对比,本发明的有益效果如下:
(1)本发明通过丢弃初始误差较大的数据后,直接提高采集数据的精确度。
(2)本发明通过设置置信区间,丢弃不在当前总体平均值置信区间内的数据,增加采集数据的可信度,从而提高采集数据的精确度。
(3)本发明多次采用置信区间进行PWTT数据筛选,防止了在数据采集的过程中出现异常数据的抖动,而出现误差较大的数据,采用置信区间的方法会筛掉不在置信区间的数据,使整体数据精确度更高,从而使后续最终计算的PWV结果更精确。
(4)本发明在血压计算方法中,添加了心率因素,使得血压计算更准确。
附图说明
图1是本发明的一个实施例的PWV测量方法流程;
图2是本发明的一个实施例中采集PWTT时,PWTT值的波动示例图;
图3是本发明的一个实施例的动态连续血压测量方法流程图。
具体实施方式
下面结合附图1和3,对本发明的具体实施方式作详细的说明。
本发明设置了置信度用于筛选合格数据,由采集的PWTT样本所构造的总体参数的估计区间。置信度展现的是PWTT的真实值有一定概率落在PWTT采集值结果的周围的程度,其给出的是PWTT测量值的可信程度。
根据本发明的一个具体实施方案,基于PWV进行血压测量时的血压计算包括如下步骤:
步骤一,采集多组校准数据,每组校准数据包括如下参量:收缩压SBP、舒张压DBP、脉搏波传导时间PWTT和心率PULSE;
该步的采集数据会用于后续步骤的计算参数,其中,采集的N组数据,理论上是越多对于最终测量血压值结果计算越精确,但是实际中基于存储资源的考量和采集效率及成本的考量,会取一个适当的折中值,从实际测量效果来看,一般第一次校准取7次左右即可计算得到符合国际标准的血压测量值。分别选取4组静止和3组运动状态的数据主要是为了涵盖测试者在静止和运动状态时的生理参数信息,从而为后续获取血压与心率的回归系数做准备。
步骤二,计算各组校准数据中脉搏波波速PWV值,并计算各组校准数据的脉压值PP;在计算之前,对PWTT数值进行筛选和补充,根据保留的PWTT值和补充的PWTT值,在各组校准数据中保留或补充对应于保留或补充的PWTT值的同组的SBP、DBP、PULSE值,其中,脉压值PP=SBP–DBP;
初采集时,数据不稳定,且在测量过程中也会出现数据的拨动,因此在计算PWV时,会去掉一些不稳定造成的不可信PWTT数据,同时,也删除与其对应的同组其他数据;同时为了保证数据量足够,不具有随机性,还要在计算过程中补充采集一些数据,同理,也对应地,将补充的PWTT对应的同组其他数据补充作为校准数据。
血压具有波动性,收缩压(SBP)和舒张压(DBP)分别代表血压波动的两个极端值,收缩压与舒张压之差称为脉压,它反映了一个心动周期中该波动幅度的大小。该步计算的校准脉压值会用于后续测量血压的算法中部分参数;
步骤三,分别计算各组校准数据中SBP与PWV、SBP与PULSE、PWV与PULSE、PP与PWV、PP与PULSE之间的自相关系数;
步骤四,分别计算各组校准数据中SBP与PWV、SBP与PULSE、PP与PWV、PP与PULSE之间的回归系数P1、Q1、P2和Q2;
步骤五,根据校准数据,计算收缩压的回归常数M1和舒张压的回归常数M2;
M1=B-P1*A-Q1*U;
M2=C-P2*A-Q2*U;
其中,
A为各组校准数据的PWV的均值;
B为各组校准数据的收缩压的均值;
C为各组校准数据的脉压值的均值;
U为各组校准数据的心率的均值;
具体来说,回归常数项表示的是未被自变量解释的且长期存在(非随机)的部分,即信息残留与随机误差之和,而随机误差是在自变量解释空间内,预测值和去掉常数项的实际值的误差。在实际血压测量中,当人在运动时,发现测量血压值还与心率呈现一定的线性回归关系,故在计算回归常数时将心率因子也作为自变量添加到常规的血压计算公式中,而回归常数则是通过校准数据由步骤一到步骤三所计算出的结果参数加上以上描述血压计算公式而得出。该回归常数值作为后续测量血压计算公式的参数,将直接影响到后续测量血压值的精确度。
步骤六,实时采集心率值PULSE_rt和PWTT值,根据实时采集的PWTT值与校准数据计算PWV值PWV_rt,根据如下回归方程计算实时的收缩压SBP和舒张压DBP;
SBP=M1+P1*PWV_rt+Q1*PULSE_rt;
DBP=M2+P2*PWV_rt+Q2*PULSE_rt;
其中,
P1、Q1、P2和Q2分别为各组校准数据的SBP与PWV、SBP与PULSE、PP与PWV、PP与PULSE之间的回归系数;
M1为收缩压的回归常数;
M2为舒张压的回归常数;
PWV_rt为根据实时采集的PWTT值与校准数据计算的PWV值;
PULSE_rt为实时采集的心率值。
根据本发明的一个具体实施方案,提供了基于PWTT值或PWV值进行血压测量的方法,其中,基于PWTT值进行血压测量,需根据PWTT值进行PWV的计算,基于PWTT值或PWV值进行血压测量的方法涉及的PWV的计算包括如下步骤:
PWTT初步采集及初筛步骤;
S1000:采集N个PWTT值,丢弃前面P个PWTT值;
当采集到N个PWTT值时,开始PWV的计算,但由于在开始采集数据时,硬件设备会有一个调整阶段而产生数据抖动,从而会使开始采集的数据误差较大,故需要丢弃前P个PWTT值;
S2000:计算剩余(N-P)个PWTT值的均值,根据置信度M%计算第一置信区间;
计算均值后,根据置信度M%,可以得到初步筛选时的置信区间;落在该区间内的PWTT值,为可信数据,将予以保留;不在该区间的数据,为不可信数据,将予以丢弃。这样,使得保留下来的PWTT数据,都是可信数据,保证了后续PWV计算的精确度。
S3000:丢弃剩余(N-P)个PWTT值中不在第一置信区间内的PWTT值,如果剩余的PWTT个数满足预设的最低PWTT个数要求,则执行PWV计算步骤,如果剩余的PWTT个数不满足预设的最低PWTT个数要求,则执行PWTT继续采集及筛选步骤;
用于PWV计算需要有足够的PWTT样本量,因此预设一个最低PWTT个数;该步骤在丢弃初始误差较大的数据后,直接提高采集数据的精确度;其次,通过设置置信区间,丢弃不在当前总体平均值置信区间内的数据,增加采集数据的可信度,从而提高采集数据的精确度。之所以采用置信区间,是为了防止在数据采集的过程中仍然会出现异常数据的抖动而出现误差较大的数据,采用置信区间的方法会使采集到的整体数据精确度更高,从而使后续最终计算的结果更精确。
PWTT继续采集及筛选步骤;
S4000:继续采集m个PWTT值,使得PWTT总个数达到预设的最低PWTT个数;
该步骤主要是当PWTT初步采集及初筛步骤的初步筛选后,PWTT总个数未达到预设的最低PWTT个数时,需要进行多次循环数据采集,直到采集到规定的数据集。
S5000:计算PWTT值的均值,根据置信度M%计算第二置信区间;每次迭代计算的所述第二置信区间与当前参与置信区间计算的PWTT值有关;
补充采集后,数据的整体区间发生变化,因此需要重新计算当前所有PWTT值的均值,并得到新的置信区间,该过程会持续迭代,直至补充数据后再次进行根据置信区间筛选数据后剩余的PWTT值满足预设的最低PWTT个数要求。
S6000:丢弃PWTT值不在第二置信区间内的PWTT值;
S7000:如果剩余的PWTT值个数满足预设的最低PWTT个数要求,则执行PWV计算步骤,否则,继续执行步骤S4000;
该步骤主要是当PWTT初步采集及初筛步骤第一次未采集到规定的有效值时进行多次循环数据采集和当前总体平均值置信区间的计算,直到采集到规定的数据集。之所以采用置信区间,是为了防止在数据采集的过程中仍然会出现异常数据的抖动而出现误差较大的数据,采用置信区间的方法会使采集到的整体数据精确度更高。此时因为参与置信区间计算的数据已经变化,因此,该置信区间与PWTT初步采集及初筛步骤中的不同,且在每次计算时,参与置信区间计算的数据变化时,置信区间都会变化。
PWV计算步骤;
S8000:计算保留的PWTT值的均值T;
S9000:计算PWV值;其中PWV值采用如下公式进行计算;
其中,
PWV为最终计算得到的脉搏波波速;
L为待测者的臂长;
a为正常人肩膀到心脏的平均距离;
T为最终保留的PWTT值的均值。
该步骤主要是通过速度等于距离除以时间的原理来计算最终的PWV值,其中时间是上述步骤采集到的精确PWTT平均值,距离为人体臂长与肩膀到心脏的距离。在计算时通过PWTT、臂长与正常人肩膀到心脏的平均距离得到PWV,使得最终计算的PWV精度更高。
本发明的PWV测量方法用于无创连续血压测量,其充分考虑了与待测量人的臂长以及正常人肩膀到心脏的平均距离,根据保留的可信PWTT数据的均值作为脉搏传播时间,用待测者的臂长与正常人肩膀到心脏的平均距离的和作为脉搏传播距离,根据式(1)得到PWV值。
根据本发明的一个实施方案,所述初步采集的PWTT值的个数N满足,5<N≤15。
根据本发明的一个实施方案,所述初步采集的PWTT值的个数N=15。
根据本发明的一个实施方案,所述初步筛选时丢弃的PWTT值的个数P满足,5≤P≤9。
根据本发明的一个实施方案,所述初步筛选时丢弃的PWTT值的个数P=5。
根据本发明的一个实施方案,所述预设的最低PWTT个数为N-P。
根据本发明的一个实施方案,置信区间(A1,A2)的计算采用如下公式:
A1=A-A*M%;
A2=A+A*M%;
其中,
A为当前保留的PWTT值的均值;
M%为置信度;
A1为置信区间下限;
A2为置信区间上限。
根据本发明的一个具体实施方案,各组校准数据中SBP与PWV、SBP与PULSE、PWV与PULSE、PP与PWV、PP与PULSE之间自相关系数的计算包括如下步骤:
S3100:计算各组校准区数据中各参量的均值、标准差和协方差;
平均值,统计学术语,是表示一组数据集中趋势的量数,是指在一组数据中所有数据之和再除以这组数据的个数。它是反映数据集中趋势的一项指标。该步通过计算各组数据的平均值,通过统计学方法排除极大极小值对结果的影响,使数据更趋于正常值而非极端值,这样会使后续计算结果更加精确。
标准差,是离均差平方的算术平均数(即:方差)的算术平方根,用σ表示。标准差也被称为标准偏差,或者实验标准差,在概率统计中最常使用作为统计分布程度上的测量依据。该步计算的标准差能反映一个数据集的离散程度,该计算结果会作为后续计算血压算法中的一个中间参数,该参数会直接影响到测量血压值的精确度。
协方差在概率论和统计学中用于衡量两个变量的总体误差。而方差是协方差的一种特殊情况,即当两个变量是相同的情况。该步计算的协方差在概率论和统计学中用于衡量两个变量的总体误差,该计算结果会作为参数用于后续血压算法中的相关系数计算和回归系数计算,直接影响到回归相关系数的相关程度,从而最终间接影响到到测量血压值的精确度。
S3200:根据各组校准数据中各参量的均值、标准差和协方差,分别计算SBP与PWV、SBP与PULSE、PWV与PULSE、PP与PWV、PP与PULSE之间的自相关系数。
根据本发明的一个具体实施方案,步骤S3200中自相关系数计算包括如下步骤:
SBP与PWV之间的自相关系数:SBP_Ray=G1/(D*E);
SBP与PULSE之间的自相关系数:SBP_Rby=S1/(E*J);
PWV与PULSE之间的自相关系数:Rab=N1/(D*J);
PP与PWV之间的自相关系数:DIS_Ray=G2/(D*F);
PP与PULSE之间的自相关系数:DIS_Rby=S2/(F*J);
其中,D为各组校准数据中PWV标准差;E为各组校准数据中SBP标准差,F为各组校准数据中PP标准差;J为各组校准数据中心率标准差;各组校准数据的D和E的协方差为G1;各组校准数据的D和F的协方差为G2;各组校准数据的D和J的协方差为N1;各组校准数据的F和J的协方差为S2。
相关系数是最早由统计学家卡尔·皮尔逊设计的统计指标,是研究变量之间线性相关程度的量,计算原理如下:
其中,Cov(X,Y)为X与Y的协方差,Var[X]为X的方差,Var[Y]为Y的方差。
根据本发明的一个具体实施方案,计算各组校准数据的各参量的自相关系数中,需要先计算各组校准数据中各参量的均值、标准差和协方差。
将各组校准数据的各参量的均值、标准差和协方差定义如下:
表1 各参量均值、标准差、协方差符号定义
该步计算的各个相关系数值主要是为了计算两组不同变量之间线性相关程度,会作为参数用于后续回归系数的计算,该参数会直接影响回归系数的相关性,从而最终间接影响到测量血压值的精确度。
根据本发明的一个具体实施方案,各组校准数据中SBP与PWV、SBP与PULSE、PP与PWV、PP与PULSE之间的回归系数P1、Q1、P2和Q2计算包括如下步骤:
P1=(SBP_Ray-SBP_Rby*Rab)/(1-Rab*Rab)*(E/D);
Q1=(SBP_Rby-SBP_Ray*Rab)/(1-Rab*Rab)*(E/J);
P2=(DIS_Ray-DIS_Rby*Rab)/(1-Rab*Rab)*(F/D);
Q2=(DIS_Rby-DIS_Ray*Rab)/(1-Rab*Rab)*(F/J);
其中:D为各组校准数据的PWV标准差;E为各组校准数据的SBP标准差,F为各组校准数据的PP标准差;J为各组校准数据的心率标准差;SBP_Ray为各组校准数据的SBP与PWV之间的自相关系数;SBP_Rby为各组校准数据的SBP与PULSE之间的自相关系数;Rab为各组校准数据的PWV与PULSE之间的自相关系数;DIS_Ray为各组校准数据的PP与PWV之间的自相关系数;DIS_Rby为各组校准数据的PP与PULSE之间的自相关系数。
回归系数在回归方程中表示自变量x对因变量y影响大小的参数。回归系数越大表示x对y影响越大,正回归系数表示y随x增大而增大,负回归系数表示y随x增大而减小。
该步计算的各个回归系数值主要是为了计算两组不同变量之间线性相关程度,会作为参数用于后续回归常数值的计算,该参数会直接影响回归常数的准确性,从而最终间接影响到测量血压值的精确度。
实施例1
根据本发明的一个具体实施方案,下表是同一测试者,固定臂长、N、M和a参数,分别计算丢弃前不同个数的PWTT值,对最终PWV的值的影响。
其中,真实PWV=3.566243,臂长L=630mm,N=15,M=10%,a=200mm。
表2 初筛时丢弃不同个数PWTT值对PWV的影响
真实PWV值采用某品牌AECG100 ECG/PPG及PWTT多功能生理讯号测试仪测得,根据上表1第2列,可得到图2。根据图2可以看到,前5个采集点PWTT数据波动比较大,5个采集点后数据趋于平稳,因此,在P=5-9时,即丢弃前5个到前9个PWTT值时,PWV计算结果更接近与真实值,而在不丢弃或者丢弃少于5个PWTT值时,PWV计算结果偏离真实值更多,且丢弃越少偏离越多。在丢弃5个之后,PWV计算值基本趋于稳定,均接近于真实值。因此将P的取值范围定为5≤P≤9,更优的,为了节省计算时间,P可以为5。
同时,可以看到,在N=6以后,PWTT值基本已经稳定,因此,N取值为5<N≤15,更优地,为保证样本量足够,N可以为15。
实施例2
根据本发明的一个具体实施方案,下表是同一测试者,固定臂长、N、M和a参数,相比于实施例1,改变置信度M值。
其中,真实PWV=3.566243,臂长L=630mm,N=15,M=5%,a=200mm。
表3 置信度为5%时的PWV计算结果
在实施例1的基础上,仅改变置信度M的值,将置信度提高到5%时,根据本发明的PWV计算步骤,舍弃不在置信区间内的数据,同时丢弃后补充采集PWTT数据,最后可以看到,P=0,P=3,P=5,P=7时,相比实施例1置信度M为10%时,PWV最终计算结果更接近于真实值,精度都有所提高。
由于用户在不同硬件环境和开发环境下,采集到的PWTT值的精确度及抖动情况是不一样的,所以用户可以根据实际的硬件环境和开发环境来设置置信区间M值来限定误差范围,从而获取到对应精度的PWV值范围。在追求高精度的时候,可以适当调低置信度M,可以使得PWV最终计算结果更接近于真实值,精度有所提高,但缺点是会消耗更多的计算时间和计算资源;而为了提高PWV的计算速度,可以提高置信度M,但是会降低PWV的精度,用户可以根据实际的需求场景来调节M值对最终结果进行取舍。
实施例3
下面给出传统方法中计算血压未加入心率回归系数时测量的血压与本发明中加入心率回归系数的方案的血压对比。
一、未添加心率回归系数测量的血压与标准血压计测量值对比
表4 测试者不同阶段(校准/测量)以及不同状态(静止/运动)时通过校准血压计所测得的收缩压SBP和舒张压DBP的血压值
表5PWTT血压测量手表中采用的未添加心率计算公式中的m值和p值
未添加心率的血压计算公式为:血压值=m+p*PWV;
表格中:
p为SBP或PP与PWV的回归系数;
m为SBP或PP的回归常数,为SBP或PP的平均值-p×PWV的平均值。
表6 测试者不同状态(静止/运动)时通过PWTT血压测量手表所测得的PWTT值、PWV值、未添加心率的血压计算算法计算的收缩压SBP和舒张压DBP
表7 测试者在不同状态(静止/运动)时通过校准血压计测的血压值与PWTT血压测量手表采用的未添加心率算法计算的血压值之间的误差
采用未添加心率的血压计算方法测得的SBP平均值误差为-4.1mmHg;SBP标准偏差为9.5mmHg;
采用未添加心率的血压计算方法测得的DBP平均值误差为-2.5mmHg;DBP标准偏差为4.1mmHg;
上述数据可以看到在校准后,通过校准数据所计算的系数,再通过校准系数通过未添加心率的计算算法在测量阶段计算收缩压SBP和计算舒张压DBP,而上表显示测量结果标准偏差误差为9.5,超过国际标准8mmHg范围的要求。
一、本发明中添加了心率回归系数测量的血压与标准血压计测量值对比
标准血压计测量数据见表2。
表8 本发明添加了心率计算公式中的m值、p值和q值
SBP | PP | |
p | -9.8588 | -11.0822 |
q | 1.4522 | 0.9990 |
m | 64.5125 | 36.1852 |
将本发明添加心率后的血压计算公式
SBP=M1+P1*PWV_rt+Q1*PULSE_rt;
DBP=M2+P2*PWV_rt+Q2*PULSE_rt;
写成血压计算通用公式为:血压值=m+p*PWV+q*PULSE;
其中,
m为回归常数,分别对应根据各组校准数据得到的收缩压回归常数M1和舒张压回归常数M2;
p为回归系数,分别对应根据各组校准数据得到的SBP与PWV、PP与PWV之间的回归系数P1和P2;
q为回归系数,分别对应根据各组校准数据得到的SBP与PULSE、PP与PULSE之间的回归系数Q1和Q2;
表9 测试者不同阶段(校准/测量)以及不同状态(静止/运动)时通过采用本发明方法的PWTT血压测量手表所测得的PWTT值、PWV值、心率PULSE值、本发明算法计算的收缩压SBP、舒张压DBP和脉压PP
表10测试者在不同状态(静止/运动)时通过校准血压计测的血压值与本发明添加心率后计算的血压值之间的误差,以及最终结果的误差平均值和最终结果的误差标准偏差
采用未添加心率的血压计算方法测得的SBP平均值误差为-1.6mmHg;SBP标准偏差为8mmHg;
采用未添加心率的血压计算方法测得的DBP平均值误差为-2.2mmHg;DBP标准偏差为3.4mmHg;
采用本发明的方法增加心率数据后,从上述数据可以看到在校准后,通过校准数据所计算的系数,再通过校准系数通过添加心率的计算算法在测量阶段计算收缩压SBP和计算舒张压DBP,而上表显示测量结果最终平均值和标准偏差误差均在国际标准8mmHg范围的要求。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均包含在本发明的保护范围之内。
Claims (4)
1.一种基于PWV的血压测量手表,其特征在于,所述血压测量手表进行所述PWV的计算,包括如下步骤:
PWTT初步采集及初筛步骤;
S1000:采集N个PWTT值,丢弃前面P个PWTT值;所述初步采集的PWTT值的个数N满足,5<N≤15;所述初步筛选时丢弃的PWTT值的个数P满足,5≤P≤9;
S2000:计算剩余(N-P)个PWTT值的均值,根据置信度M%计算第一置信区间;置信度根据实际的需求场景进行调节;
S3000:丢弃剩余(N-P)个PWTT值中不在第一置信区间内的PWTT值,如果剩余的PWTT个数满足预设的最低PWTT个数要求,则执行PWV计算步骤,如果剩余的PWTT个数不满足预设的最低PWTT个数要求,则执行PWTT继续采集及筛选步骤;所述预设的最低PWTT个数为N-P;
PWTT继续采集及筛选步骤;
S4000:继续采集m个PWTT值,使得PWTT总个数达到预设的最低PWTT个数;
S5000:计算PWTT值的均值,根据置信度M%计算第二置信区间;每次迭代计算的所述第二置信区间与当前参与置信区间计算的PWTT值有关;
S6000:丢弃PWTT值不在第二置信区间内的PWTT值;
S7000:如果剩余的PWTT值个数满足预设的最低PWTT个数要求,则执行PWV计算步骤,否则,继续执行步骤S4000;
PWV计算步骤;
S8000:计算保留的PWTT值的均值S;
S9000:计算PWV值;其中PWV值采用如下公式进行计算;
其中,
PWV为最终计算得到的脉搏波波速;
L为待测者的臂长;
a为正常人肩膀到心脏的平均距离;
S为最终保留的PWTT值的均值。
2.根据权利要求1所述的基于PWV的血压测量手表,其特征在于,所述初步采集的PWTT值的个数N=15。
3.根据权利要求1所述的基于PWV的血压测量手表,其特征在于,所述初步筛选时丢弃的PWTT值的个数P=5。
4.根据权利要求1所述的基于PWV的血压测量手表,其特征在于,置信区间(A1,A2)的计算采用如下公式:
A1=A-A*M%;
A2=A+A*M%;
其中,
A为当前保留的PWTT值的均值;
M%为置信度;
A1为置信区间下限;
A2为置信区间上限。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110246028.3A CN113171067B (zh) | 2021-03-05 | 2021-03-05 | 一种基于pwv的血压测量手表 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110246028.3A CN113171067B (zh) | 2021-03-05 | 2021-03-05 | 一种基于pwv的血压测量手表 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113171067A CN113171067A (zh) | 2021-07-27 |
CN113171067B true CN113171067B (zh) | 2022-03-15 |
Family
ID=76921947
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110246028.3A Active CN113171067B (zh) | 2021-03-05 | 2021-03-05 | 一种基于pwv的血压测量手表 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113171067B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1442110A (zh) * | 2002-03-01 | 2003-09-17 | 客林公司 | 动脉硬化检查设备 |
CN102397064A (zh) * | 2011-12-14 | 2012-04-04 | 中国航天员科研训练中心 | 连续血压测量装置 |
CN106725396A (zh) * | 2016-12-08 | 2017-05-31 | 上海交通大学 | 一种基于双脉搏波的无创逐拍血压测量装置以及测量方法 |
EP3251588A1 (en) * | 2016-05-31 | 2017-12-06 | Tata Consultancy Services Limited | Method and system for physiological parameter derivation from pulsating signals with reduced error |
CN108366749A (zh) * | 2015-10-12 | 2018-08-03 | 西北大学 | 动态血压与生命体征监测装置、系统和方法 |
CN110403580A (zh) * | 2018-04-28 | 2019-11-05 | 深圳市大耳马科技有限公司 | 一种脉搏波传导参数测量方法和脉搏波传导参数处理设备 |
CN112040849A (zh) * | 2018-04-13 | 2020-12-04 | 深圳市长桑技术有限公司 | 用于确定对象血压的系统和方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20080090194A (ko) * | 2007-04-04 | 2008-10-08 | 엘지전자 주식회사 | 혈압 검출 방법 및 그 장치 |
US20170119264A1 (en) * | 2015-11-04 | 2017-05-04 | Sharp Laboratories of America (SLA), Inc. | Pulse Wave Velocity-to-Blood Pressure Calibration Prompting |
CN106618537B (zh) * | 2016-12-21 | 2020-09-01 | 天津普仁万合信息技术有限公司 | 一种基于脉搏波传导的连续动态血压监测装置和方法 |
CN110037673B (zh) * | 2019-05-13 | 2020-12-01 | 深圳六合六医疗器械有限公司 | 血压个性化区间的统计方法及装置 |
-
2021
- 2021-03-05 CN CN202110246028.3A patent/CN113171067B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1442110A (zh) * | 2002-03-01 | 2003-09-17 | 客林公司 | 动脉硬化检查设备 |
CN102397064A (zh) * | 2011-12-14 | 2012-04-04 | 中国航天员科研训练中心 | 连续血压测量装置 |
CN108366749A (zh) * | 2015-10-12 | 2018-08-03 | 西北大学 | 动态血压与生命体征监测装置、系统和方法 |
EP3251588A1 (en) * | 2016-05-31 | 2017-12-06 | Tata Consultancy Services Limited | Method and system for physiological parameter derivation from pulsating signals with reduced error |
CN106725396A (zh) * | 2016-12-08 | 2017-05-31 | 上海交通大学 | 一种基于双脉搏波的无创逐拍血压测量装置以及测量方法 |
CN112040849A (zh) * | 2018-04-13 | 2020-12-04 | 深圳市长桑技术有限公司 | 用于确定对象血压的系统和方法 |
CN110403580A (zh) * | 2018-04-28 | 2019-11-05 | 深圳市大耳马科技有限公司 | 一种脉搏波传导参数测量方法和脉搏波传导参数处理设备 |
Non-Patent Citations (1)
Title |
---|
"基于智能手机的无袖带血压测量方法";赵冠斌;《中国优秀博硕士学位论文全文数据库(硕士) 医药卫生科技辑》;20190131;第14-15、32-33页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113171067A (zh) | 2021-07-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102429649B (zh) | 连续血压测量装置 | |
CN102397064B (zh) | 连续血压测量装置 | |
CN102488503B (zh) | 连续血压测量装置 | |
CN100528073C (zh) | 基于动脉压力的心血管参数的自动确定 | |
Tanner | The construction of normal standards for cardiac output in man | |
CN106821356A (zh) | 基于Elman神经网络的云端连续血压测量方法及系统 | |
CN112088408A (zh) | 用于睡眠阶段检测的方法、计算设备和可穿戴设备 | |
CN110558960A (zh) | 一种基于ptt和miv-ga-svr的连续血压无创监测方法 | |
CN114788703B (zh) | 高频qrs波形数据分析方法、装置、计算机设备与存储介质 | |
CN104644151B (zh) | 一种基于光电容积脉搏信号的压力脉搏波波形传播预测方法 | |
CN113171067B (zh) | 一种基于pwv的血压测量手表 | |
CN113171070B (zh) | 一种基于pwtt的动态连续血压测量手表 | |
CN103505191A (zh) | 一种利用压脉带的压力脉波震荡讯号以预估中心主动脉脉搏压的方法及其装置 | |
US20200359958A1 (en) | Sleep sufficiency estimation device and sleep sufficiency estimation method | |
CN113470805A (zh) | 建立血压模型的方法 | |
CN113171068B (zh) | 一种基于血压虚拟校准值的动态连续血压测量手表 | |
CN113171071B (zh) | 一种基于pwtt的血压测量手表 | |
CN113180621B (zh) | 基于freeRTOS的连续无创血压测量系统 | |
CN113499048A (zh) | 一种基于CNN-BiLSTM的中心动脉压波形重构的系统及方法 | |
CN111671409A (zh) | 一种脉搏波形的分析方法及装置 | |
CN113171069B (zh) | 一种嵌入式图形血压测量系统 | |
EP2976999B1 (en) | Circulatory organ function arithmetic calculation device | |
CN111387939B (zh) | 心排量的测量方法、装置、计算机设备和存储介质 | |
WO2023190085A1 (ja) | 血圧推定装置、血圧推定方法およびプログラム | |
US20240057878A1 (en) | Continuous self-calibrating blood pressure monitoring device and method for its use |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |