CN110420014A - 基于脉搏波智能分解的可穿戴心血管功能评价系统及方法 - Google Patents
基于脉搏波智能分解的可穿戴心血管功能评价系统及方法 Download PDFInfo
- Publication number
- CN110420014A CN110420014A CN201910509555.1A CN201910509555A CN110420014A CN 110420014 A CN110420014 A CN 110420014A CN 201910509555 A CN201910509555 A CN 201910509555A CN 110420014 A CN110420014 A CN 110420014A
- Authority
- CN
- China
- Prior art keywords
- wave
- radial artery
- pulse wave
- pulse
- 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.)
- 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/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/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- 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
-
- 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/022—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers
- A61B5/0225—Measuring pressure in heart or blood vessels by applying pressure to close blood vessels, e.g. against the skin; Ophthalmodynamometers the pressure being controlled by electric signals, e.g. derived from Korotkoff sounds
-
- 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/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/6813—Specially adapted to be attached to a specific body part
- A61B5/6824—Arm or wrist
-
- 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/683—Means for maintaining contact with the body
- A61B5/6831—Straps, bands or harnesses
-
- 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
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Signal Processing (AREA)
- Vascular Medicine (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Power Engineering (AREA)
- Ophthalmology & Optometry (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明提供一种基于脉搏波智能分解的可穿戴心血管功能评价系统及方法,涉及生物医学电子健康领域。包括桡动脉血流阻断单元、脉搏波与前向波采集单元、控制单元、传输单元、反射波计算单元和心血管功能评估单元;脉搏波与前向波采集单元采集脉搏波和前向波波形;控制单元用于控制双袖带分别达到血流阻断或前向波测量的要求;传输单元用于将采集到的数据传输到反射波计算单元;反射波计算单元通过前向波计算得到脉搏波的反射波;心血管功能评估单元根据反射波计算反映心血管功能的相关参数和波形的实时显示,并根据参数指标评估心血管功能。本发明实现对脉搏波进行分解和反射波提取,计算前向波与反射波相关指标从而对心血管功能进行有效评估。
Description
技术领域
本发明涉及生物医学电子健康技术领域,尤其涉及一种基于脉搏波智能分解的可穿戴心血管功能评价系统及方法。
背景技术
心血管疾病负担日渐加重,已成为重大的公共卫生难题,防治心血管疾病刻不容缓。采用心血管功能评价可以有效早期诊断和干预心血管疾病。动脉硬化程度、顺应性等动脉相关参数可以反映心血管系统功能,因此早期筛查动脉硬化等疾病是心血管病预防的根本途径。
基于脉搏波波形分析的血流动力学指标检测,是一种常用的无创心血管功能评估方法,利用这种原理的医疗仪器均存在价格昂贵、操作过程烦琐、且对操作技能要求较高等应用瓶颈,难以在大众健康领域中普及推广。因此,建立适用于群体性普查和个体化动态检测的动脉硬化检测方法和技术十分必要。
专利号为201811213852.3的发明专利提出了一种基于脉搏波特征参数对人体健康状况的检测方法,通过便携的腕式脉搏检测装置对人体脉搏波进行检测,并提取特征,对脉搏波特征进行样本训练,从而判断人体的健康状况。但这项发明专利只在时、频域对脉搏波信号进行了特征提取,并没有对脉搏波进行分解并根据反射波评估动脉硬化情况;专利号为 201010558798.3的发明专利利用小波变换来提取脉搏波的起始点,计算脉搏波波速对血管硬化程度进行检测,但该专利只完成了脉搏波分解中的起始点提取没有提取出脉搏波的反射波进行深入计算;专利号201510955386.6为的发明专利提出了一种基于FMD的血管内皮功能评估检测系统及方法,采集左指尖和右指尖的血容积脉搏波,对手臂的肱动脉进行加压阻断,但该方法只可以计算血管内皮功能指数,没有对袖带加压的方式分解脉搏波这一方面进行分析。
发明内容
本发明要解决的技术问题是针对上述现有技术的不足,提供一种基于脉搏波智能分解的可穿戴心血管功能评价系统及方法,实现更有效评估动脉及心血管的功能。
为解决上述技术问题,本发明所采取的技术方案是:一方面,本发明提供一种基于脉搏波智能分解的可穿戴心血管功能评价系统,包括桡动脉血流阻断单元、脉搏波与前向波采集单元、控制单元、传输单元、反射波计算单元和心血管功能评估单元;所述桡动脉血流阻断单元包括桡动脉阻断袖带和压力反馈传感器;所述桡动脉血流阻断袖带佩戴于桡动脉,用于阻断桡动脉血流,并通过压力反馈传感器测量反馈桡动脉处压力状态;所述脉搏波与前向波采集单元用于采集不同压力下脉搏波和前向波;所述控制单元用于控制桡动脉阻断袖带压力,当压力超过自适应阈值,或压力脉搏波波动幅值小于自适应阈值时,控制桡动脉阻断袖带停止加压,达到血流阻断或前向波测量的要求;所述传输单元用于将脉搏波与前向波采集单元采集到的数据传输到反射波计算单元和心血管功能评估单元;所述反射波计算单元用于将测得的脉搏波数据进行计算获得脉搏波的反射波形;所述心血管功能评估单元用于根据反射波计算反映心血管功能的参数及反射波形的实时显示,并根据前向波与反射波的综合分析评估心血管功能;
所述桡动脉阻断袖带,用于阻断血流和脉搏波向手臂末端传导,从而改变脉搏波终端反射的位置,在桡动脉处实现脉搏波的全反射;所述压力反馈传感器,用于精确测量袖带内部的气体压力,并将压力值实时反馈给控制单元;
所述脉搏波与前向波采集单元包括脉搏波与前向波测量模块和信号调理模块;所述脉搏波与前向波测量模块包括可调压力测量袖带和前向波测量传感器;可调压力测量袖带紧贴放在桡动脉阻断袖带的近心端,用于分时精确采集:桡动脉阻断袖带不加压时的脉搏波波形,及桡动脉阻断袖带加压后的前向波及其在桡动脉阻断袖带的全反射波形的叠加,实现不同压力下的测量;所述前向波测量传感器,用于精准测量不同脉象压力下的前向波,并将压力波形实时传输到信号调理模块;所述信号调理模块,用于将前向波测量传感器得到的前向波波形进行信号滤波与调理;
所述反射波计算单元,通过正常状态下的脉搏波信号与桡动脉阻断状态下测量的前向波及其在桡动脉阻断袖带的全反射波形的叠加,计算得到脉搏波的反射波;
所述心血管功能评估单元包括参数计算与波形显示模块和心血管功能评价模块;所述参数计算与波形显示模块,根据反射波计算单元得到的反射波计算心血管相关参数,包括动脉反射指数RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr和动脉反射波增强指数AIx,并对直接采集到的脉搏波和前向波进行实时显示;所述心血管功能评价模块通过参数计算与波形显示模块得到的心血管相关参数,综合分析前向波与反射波,对心血管功能和健康程度进行评估和显示。
优选地,所述桡动脉血流阻断单元和脉搏波与前向波测量模块,均包括两条充气袖带、微型气泵和电磁阀;所述充气袖带由高弹性橡胶材质制成,袖带包括充气气囊、充气口和细软管,微型气泵和电磁阀通过细软管连接充气口为充气气囊充气,保证充气气囊内压力处处相等,两条充气袖带通过魔术贴绑在桡动脉处,桡动脉血流阻断袖带在桡动脉处,脉搏波与前向波测量模块的可调压力测量袖带在紧挨桡动脉血流阻断袖带的近心端;所述微型气泵和电磁阀连接控制单元,根据控制单元设定的压力要求为充气气囊充气。
优选地,所述信号调理模块,包括滤波电路和AD转换装置;所述滤波电路和AD转换装置均与桡动脉血流阻断单元的输出端相连,为桡动脉血流阻断单元的输出信号进行调理和 AD转换,AD转换后的信号经传输单元传输到反射波计算单元和心血管功能评估单元。
优选地,所述控制单元采用微控制器,根据压力反馈传感器反馈的桡动脉阻断袖带内实时压力值判断是否达到阻断的要求,从而控制微型气泵是否继续充气,实现精确控制并调节加压压力稳定测量袖带内压力,并在微控制器的I/O口上连接四个三极管元件实现控制信号与供电电压隔离,加速微型气泵的充气速率;同时控制单元通过控制桡动脉血流阻断单元的桡动脉阻断袖带加压时间,先后分时测量同一位置的脉搏波和阻断后的前向波。
优选地,所述传输模块包括存储模块和蓝牙收发模块;所述存储模块保存脉搏波与前向波采集单元采集的脉搏波和前向波,并将这些波形数据通过蓝牙收发模块传输到反射波计算单元。
优选地,所述反射波计算单元包括第一处理模块、第二处理模块、第三处理模块和第四处理模块;所述第一处理模块读取采集的测量者脉搏波和前向波,并进行波形去噪;所述第二处理模块对脉搏波和前向波信号分别去除其基线漂移,并为了反射波波形的准确计算而将脉搏波和前向波信号对齐在同一水平;所述第三处理模块将脉搏波和前向波按照第一个采集到的有效波谷位置对齐,避免因相位差导致的反射波不准确;所述第四处理模块,按照如下公式求解反射波:
Pb(t)=P(t)-Pf(t)
其中,Pb(t)为反射波波形,P(t)和Pf(t)分别为采集的脉搏波波形和前向波波形,t为时间。
优选地,所述参数计算与波形显示模块,通过分解得到的反射波,提取反射波、脉搏波和前向波相关波形特征,计算心血管动脉反射指数RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr和动脉反射波增强指数AIx。
另一方面,本发明还提供一种基于脉搏波智能分解的可穿戴心血管功能评价方法,包括以下步骤:
步骤1、调节可调压力测量袖带和桡动脉阻断袖带压力进行脉搏波采集和分解,得到脉搏波和前向波,具体方法为:
步骤1.1、控制单元控制脉搏波与前向波测量模块的可调压力测量袖带加压至70mmHg,脉搏波与前向波测量模块的前向波测量传感器采集70mmHg压力下的脉搏波,即为脉搏波波形,保存并传输至信号调理模块;
步骤1.2、控制单元控制桡动脉血流阻断单元的桡动脉阻断袖带加压至200mmHg时,同时桡动脉血流阻断单元的压力反馈传感器将桡动脉阻断袖带中的压力波形实时反馈给控制单元,当压力波形幅值小于阈值时,停止加压,此时血流在桡动脉血流阻断单元的桡动脉阻断袖带处实现完全阻断,桡动脉的前向波在此处发生完全反射,在脉搏波与前向波测量模块袖带处采集脉搏波前向波的等量叠加;
步骤1.3、脉搏波与前向波测量单元采集桡动脉血流阻断时的脉搏波和前向波等量叠加的压力波形,即两倍前向波,并进行信号调理;
步骤2、反射波计算单元对脉搏波与前向波测量模块采集到的脉搏波和两倍前向波进行处理和计算,得到反射波,具体方法为:
步骤2.1、读取采集的测量者脉搏波和两倍前向波各十个连续周期,并进行波形去噪和去基线漂移处理,并为了反射波波形的准确计算将两路波形信号对齐在同一水平,将两路波形信号各十个周期分别进行叠加并取平均,进一步保证反射波波形精准计算;
步骤2.2、将脉搏波和两倍前向波按照波谷位置水平在时间域对齐,避免因相位差导致的反射波不准确;
步骤2.3、按照如下公式求解反射波:
Pb(t)=P(t)-Pf(t)
其中,Pb(t)为反射波波形,P(t)为采集到的脉搏波波形,Pf(t)为直接采集到的两倍前向波波形除以2,即为脉搏波的前向波;
步骤3、分解得到的反射波,通过参数计算与波形显示模块,提取反射波和脉搏波、前向波相关波形特征,计算心血管相关参数:桡动脉反射指数RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr、动脉反射波增强指数AIx;
所述桡动脉反射指数RI如下公式所示:
其中,Pb为反射波压力,Pf为前向波压力;
所述桡动脉反射强度RM如下公式所示:
其中,Pb为反射波压力,Pf为前向波压力;
所述桡动脉反射波增强指数AIr为是桡动脉反射波和前向波叠加后波形的压力峰值P2与收缩早期即反射波到达前压力峰值P1的比值,如下公式所示:
所述动脉反射波增强指数AIx如下公式所示:
其中,Ps为脉搏波波形收缩期峰值,Pd脉搏波波形舒张期最低值,Pi为反射波波形中的反射波峰值;
步骤4、心血管功能评价模块根据计算得到的心血相关参数对心血管功能和健康程度进行评估和显示。
采用上述技术方案所产生的有益效果在于:本发明提供的基于脉搏波智能分解的可穿戴心血管功能评价系统及方法,采用袖带加压阻断桡动脉血流的物理方式分解脉搏波,直接测量和计算人体桡动脉处的反射波波形,提高反射波分解的准确性;微控制器根据压力自动控制袖带充放气的时间,并可以实现在不同压力下脉搏波波形的采集和分析,实现全自动采集与反射波计算;以反射波波形为依据计算相关参数,反应动脉硬化程度,通过前向波与反射波的综合分析,进行个性化心血管功能智能评估与深度解析。
附图说明
图1为本发明实施例提供的基于脉搏波智能分解的可穿戴心血管功能评价系统的结构框图;
图2为本发明实施例提供的佩戴示意图;
图3为本发明实施例提供的压力袖带结构示意图,其中,(a)为袖带正面结构示意图, (b)为袖带反面结构示意图;
图4为本发明实施例提供的压力袖带充气状态立体结构示意图;
图5为本发明实施例提供的微控制器的电路连接图,其中,(a)为微控制器的电路连接, (b)为第一个三极管与微控制器的连接,(c)为第二个三极管与微控制器的连接,(d)为第三个三极管与微控制器的连接,(e)为第四个三极管与微控制器的连接;
图6为本发明实施例提供的基于脉搏波智能分解的可穿戴心血管功能评价方法的流程图;
图7为本发明实施例提供的本发明方法分解得到的反射波与基于波形分析得到的反射波对比结果图,其中,(a)为本发明方法分解得到的反射波,(b)为基于波形分析得到的反射波。
图中,201、可调压力充气袖带;202、细软管;203、桡动脉血流阻断模块;204、脉搏波与前向波测量模块;301、袖带充气气囊;302、魔术贴;303、袖带充/放气口。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
本实施例中,基于脉搏波智能分解的可穿戴心血管功能评价系统,如图1所示,包括桡动脉血流阻断单元、脉搏波与前向波采集单元、控制单元、传输单元、反射波计算单元和心血管功能评估单元;桡动脉血流阻断单元包括桡动脉阻断袖带和压力反馈传感器;桡动脉血流阻断袖带佩戴于桡动脉,用于阻断桡动脉血流,并通过压力反馈传感器测量反馈桡动脉处压力状态;脉搏波与前向波采集单元用于采集不同压力下脉搏波和前向波;控制单元用于控制桡动脉阻断袖带压力,当压力超过自适应阈值,或压力脉搏波波动幅值小于自适应阈值时,控制桡动脉阻断袖带停止加压,达到血流阻断或前向波测量的要求;传输单元用于将脉搏波与前向波采集单元采集到的数据传输到反射波计算单元和心血管功能评估单元;反射波计算单元用于将测得的脉搏波数据进行计算获得脉搏波的反射波形;心血管功能评估单元用于根据反射波计算反映心血管功能的参数及反射波形的实时显示,并根据前向波与反射波的综合分析评估心血管功能;
桡动脉阻断袖带,用于阻断血流和脉搏波向手臂末端传导,从而改变脉搏波终端反射的位置,在桡动脉处实现脉搏波的全反射;所述压力反馈传感器,用于精确测量袖带内部的气体压力,并将压力值实时反馈给控制单元;
脉搏波与前向波采集单元包括脉搏波与前向波测量模块和信号调理模块;脉搏波与前向波测量模块包括可调压力测量袖带和前向波测量传感器;可调压力测量袖带紧贴放在桡动脉阻断袖带的近心端,用于分时精确采集:桡动脉阻断袖带不加压时的脉搏波波形,及桡动脉阻断袖带加压后的前向波及其在桡动脉阻断袖带的全反射波形(相当于两个前向波)的叠加,实现不同压力下的测量;前向波测量传感器,用于精准测量不同脉象压力下的前向波,并将压力波形实时传输到信号调理模块;信号调理模块,用于将前向波测量传感器得到的前向波波形进行信号滤波与调理;
反射波计算单元,通过正常状态下的脉搏波信号与桡动脉阻断状态下测量的前向波及其在桡动脉阻断袖带的全反射波形的叠加,计算得到脉搏波的反射波;
心血管功能评估单元包括参数计算与波形显示模块和心血管功能评价模块;参数计算与波形显示模块,根据反射波计算单元得到的反射波计算心血管相关参数,包括动脉反射指数 RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr和动脉反射波增强指数AIx,并对直接采集到的脉搏波和前向波进行实时显示;心血管功能评价模块通过参数计算与波形显示模块得到的心血管相关参数,综合分析前向波与反射波,对心血管功能和健康程度进行评估和显示。
桡动脉血流阻断单元和脉搏波与前向波测量模块,如图2所示,均包括两条充气袖带201、微型气泵和电磁阀;充气袖带由高弹性橡胶材质制成,袖带包括充气气囊、充气口和细软管 202,微型气泵和电磁阀通过细软管连接充气口为充气气囊充气,保证充气气囊内压力处处相等,两条充气袖带通过魔术贴绑在桡动脉处,桡动脉血流阻断袖带在桡动脉处,脉搏波与前向波测量模块的可调压力测量袖带在紧挨桡动脉血流阻断袖带的近心端;微型气泵和电磁阀连接控制单元,根据控制单元设定的压力要求为充气气囊充气。
本实施例中,袖带的正反面结构如图3(a)、(b)所示,袖带的立体结构如图4所示,301 为袖带充气气囊,302为魔术贴,303为袖带充/放气口,袖带上的魔术贴302将桡动脉血流阻断模块203和脉搏波与前向波测量模块204的两条袖带绑在桡动脉处,桡动脉血流阻断模块203的袖带固定在远心端;脉搏波与前向波测量模块204的袖带固定在近心端;微型气泵和电磁阀连接控制电路,根据设定压力要求为气囊充气。
信号调理模块,包括滤波电路和AD转换装置;滤波电路和AD转换装置均与桡动脉血流阻断单元的输出端相连,为桡动脉血流阻断单元的输出信号进行调理和AD转换,AD转换后的信号经传输单元传输到反射波计算单元和心血管功能评估单元。
控制单元采用微控制器,根据压力反馈传感器反馈的桡动脉阻断袖带内实时压力值判断是否达到阻断的要求,从而控制微型气泵是否继续充气,实现精确控制并调节加压压力并稳定测量袖带内压力,并在微控制器的I/O口上连接四个三极管元件如图5所示,实现控制信号与供电电压隔离,加速微型气泵的充气速率;同时控制单元通过控制桡动脉血流阻断单元的桡动脉阻断袖带加压时间,先后分时测量同一位置的脉搏波和阻断后的前向波。
传输模块包括存储模块和蓝牙收发模块;存储模块保存脉搏波与前向波采集单元采集的脉搏波和前向波,并将这些波形数据通过蓝牙收发模块传输到反射波计算单元。
反射波计算单元包括第一处理模块、第二处理模块、第三处理模块和第四处理模块;第一处理模块读取采集的测量者脉搏波和前向波,并进行波形去噪;第二处理模块对脉搏波和前向波信号分别去除其基线漂移,并为了反射波波形的准确计算而将脉搏波和前向波信号对齐在同一水平;第三处理模块将脉搏波和前向波按照第一个采集到的有效波谷位置对齐,避免因相位差导致的反射波不准确;第四处理模块,按照如下公式求解反射波:
Pb(t)=P(t)-Pf(t)
其中,Pb(t)为反射波波形,P(t)和Pf(t)分别为采集的脉搏波波形和前向波波形,t为时间。
参数计算与波形显示模块,通过分解得到的反射波,提取反射波、脉搏波和前向波相关波形特征,计算心血管动脉反射指数RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr 和动脉反射波增强指数AIx。
基于脉搏波智能分解的可穿戴心血管功能评价方法,如图6所示,包括以下步骤:
步骤1、调节可调压力测量袖带和桡动脉阻断袖带压力进行脉搏波采集和分解,得到脉搏波和前向波,具体方法为:
步骤1.1、控制单元控制脉搏波与前向波测量模块的可调压力测量袖带加压至70mmHg,脉搏波与前向波测量模块的前向波测量传感器采集70mmHg压力下的脉搏波,即为脉搏波波形,保存并传输至信号调理模块;
步骤1.2、控制单元控制桡动脉血流阻断单元的桡动脉阻断袖带加压至200mmHg时,同时桡动脉血流阻断单元的压力反馈传感器将桡动脉阻断袖带中的压力波形实时反馈给控制单元,当压力波形幅值小于阈值时,停止加压,此时血流在桡动脉血流阻断单元的桡动脉阻断袖带处实现完全阻断,桡动脉的前向波在此处发生完全反射,在脉搏波与前向波测量模块袖带处采集脉搏波前向波的等量叠加;
步骤1.3、脉搏波与前向波测量单元采集桡动脉血流阻断时的脉搏波和前向波等量叠加的压力波形,即两倍前向波,并进行信号调理;
步骤2、反射波计算单元对脉搏波与前向波测量模块采集到的脉搏波和两倍前向波进行处理和计算,得到反射波,具体方法为:
步骤2.1、读取采集的测量者脉搏波和两倍前向波各十个连续周期,并进行波形去噪和去基线漂移处理,并为了反射波波形的准确计算将两路波形信号对齐在同一水平,将两路波形信号各十个周期分别进行叠加并取平均,进一步保证反射波波形精准计算;
步骤2.2、将脉搏波和两倍前向波按照波谷位置水平在时间域对齐,避免因相位差导致的反射波不准确;
步骤2.3、按照如下公式求解反射波:
Pb(t)=P(t)-Pf(t)
其中,Pb(t)为反射波波形,P(t)为采集到的脉搏波波形,Pf(t)为直接采集到的两倍前向波波形除以2,即为脉搏波的前向波;
步骤3、分解得到的反射波,通过参数计算与波形显示模块,提取反射波和脉搏波、前向波相关波形特征,计算心血管相关参数:桡动脉反射指数RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr、动脉反射波增强指数AIx;
所述桡动脉反射指数RI如下公式所示:
其中,Pb为反射波压力,Pf为前向波压力;
所述桡动脉反射强度RM如下公式所示:
其中,Pb为反射波压力,Pf为前向波压力;
本实施例中,根据桡动脉反射强度RM的健康人群判断标准为RM数值小于0.5;
所述桡动脉反射波增强指数AIr为是桡动脉反射波和前向波叠加后波形的压力峰值P2与收缩早期即反射波到达前压力峰值P1的比值,如下公式所示:
本实施例中,根据桡动脉反射波增强指数AIr的健康人群判断标准如下:
对于18岁至30岁的测试者,判断AIr数值大于31并且小于88;
对于30岁至45岁的测试者,判断AIr数值大于50并且小于104;
对于45岁至80岁的测试者,判断AIr数值大于58并且小于105;
所述动脉反射波增强指数AIx如下公式所示:
其中,Ps为脉搏波波形收缩期峰值,Pd脉搏波波形舒张期最低值,Pi为反射波波形中的反射波峰值;
步骤4、心血管功能评价模块根据计算得到的心血相关参数对心血管功能和健康程度进行评估和显示。
本实施例还根据采集得到的脉搏波波形,通过波形分析的方法,将脉搏波分解为前向波与反射波,并将反射波计算单元通过调节袖带压力方法分解的反射波与基于波形分析得到的反射波进行比对分析;
首先读取脉搏波并对脉搏波进行平滑去噪及去除基线漂移预处理,增强脉搏波信号的质量;然后进行脉搏波特征点的识别与提取;脉搏波特征点包括脉搏波每个周期的起始点、主波峰值、切迹点和反射点;脉搏波每个周期的起始点与主波峰值通过识别脉搏波曲线的波峰波谷得到,脉搏波的切迹点与反射点通过计算脉搏波的二阶微分极值点来定位,每周期脉搏波的二阶微分第一个极大值为反射点,第二个极大值为切迹点,根据以上算法实现了脉搏波主要特征点的识别定位。最后通过得到的脉搏波特征点构造三角形以拟合血液流速波,三角形的底边长为起始点至切迹点的长度,在反射点定位成功的情况下,三角形顶点横坐标位于反射点处,在反射点找不准的情况下,三角形顶点位于底边长30%处,三角形幅值无需校准,建立三角形流速波后,根据以下公式:
Pf(t)=[P(t)+Zc×F(t)]/2
Pb(t)=[P(t)-Zc×F(t)]/2
分别计算出前向波与反射波的波形;
其中,F(t)为构造的拟合血液流速波的三角形,Zc为动脉特征阻抗;
动脉特征阻抗的计算如下公式所示:
ZC=(Pi-Pd)/Peakflow
其中,Pd为舒张末期压力值,Pi为压力曲线反射点压力值即流速波峰值对应位置,在程序中取为三角形顶点横坐标,Peakflow为三角形的高(幅值)。
对比公式可以得知,ZC是压力增量与流速增量的比值因而使用三角形拟合血流速波形时幅值无需校准,在本实施例中默认构造高为100的三角形;
本实施例中,通过本发明的反射波计算单元的调节袖带压力方法分解得到的反射波与基于波形分析得到的反射波的结果如图7所示。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (9)
1.一种基于脉搏波智能分解的可穿戴心血管功能评价系统,其特征在于:包括桡动脉血流阻断单元、脉搏波与前向波采集单元、控制单元、传输单元、反射波计算单元和心血管功能评估单元;所述桡动脉血流阻断单元包括桡动脉阻断袖带和压力反馈传感器;所述桡动脉血流阻断袖带佩戴于桡动脉,用于阻断桡动脉血流,并通过压力反馈传感器测量反馈桡动脉处压力状态;所述脉搏波与前向波采集单元用于采集不同压力下脉搏波和前向波;所述控制单元用于控制桡动脉阻断袖带压力,当压力超过自适应阈值,或压力脉搏波波动幅值小于自适应阈值时,控制桡动脉阻断袖带停止加压,达到血流阻断或前向波测量的要求;所述传输单元用于将脉搏波与前向波采集单元采集到的数据传输到反射波计算单元和心血管功能评估单元;所述反射波计算单元用于将测得的脉搏波数据进行计算获得脉搏波的反射波形;所述心血管功能评估单元用于根据反射波计算反映心血管功能的参数及反射波形的实时显示,并根据前向波与反射波的综合分析评估心血管功能;
所述桡动脉阻断袖带,用于阻断血流和脉搏波向手臂末端传导,从而改变脉搏波终端反射的位置,在桡动脉处实现脉搏波的全反射;所述压力反馈传感器,用于精确测量袖带内部的气体压力,并将压力值实时反馈给控制单元;
所述脉搏波与前向波采集单元包括脉搏波与前向波测量模块和信号调理模块;所述脉搏波与前向波测量模块包括可调压力测量袖带和前向波测量传感器;可调压力测量袖带紧贴放在桡动脉阻断袖带的近心端,用于分时精确采集:桡动脉阻断袖带不加压时的脉搏波波形,及桡动脉阻断袖带加压后的前向波及其在桡动脉阻断袖带的全反射波形的叠加,实现不同压力下的测量;所述前向波测量传感器,用于精准测量不同脉象压力下的前向波,并将压力波形实时传输到信号调理模块;所述信号调理模块,用于将前向波测量传感器得到的前向波波形进行信号滤波与调理;
所述反射波计算单元,通过正常状态下的脉搏波信号与桡动脉阻断状态下测量的前向波及其在桡动脉阻断袖带的全反射波形的叠加,计算得到脉搏波的反射波;
所述心血管功能评估单元包括参数计算与波形显示模块和心血管功能评价模块;所述参数计算与波形显示模块,根据反射波计算单元得到的反射波计算心血管相关参数,包括动脉反射指数RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr和动脉反射波增强指数AIx,并对直接采集到的脉搏波和前向波进行实时显示;所述心血管功能评价模块通过参数计算与波形显示模块得到的心血管相关参数,综合分析前向波与反射波,对心血管功能和健康程度进行评估和显示。
2.根据权利要求1所述的基于脉搏波智能分解的可穿戴心血管功能评价系统,其特征在于:所述桡动脉血流阻断单元和脉搏波与前向波测量模块,均包括两条充气袖带、微型气泵和电磁阀;所述充气袖带由高弹性橡胶材质制成,袖带包括充气气囊、充气口和细软管,微型气泵和电磁阀通过细软管连接充气口为充气气囊充气,保证充气气囊内压力处处相等,两条充气袖带通过魔术贴绑在桡动脉处,桡动脉血流阻断袖带在桡动脉处,脉搏波与前向波测量模块的可调压力测量袖带在紧挨桡动脉血流阻断袖带的近心端;所述微型气泵和电磁阀连接控制单元,根据控制单元设定的压力要求为充气气囊充气。
3.根据权利要求1所述的基于脉搏波智能分解的可穿戴心血管功能评价系统,所述信号调理模块,包括滤波电路和AD转换装置;所述滤波电路和AD转换装置均与桡动脉血流阻断单元的输出端相连,为桡动脉血流阻断单元的输出信号进行调理和AD转换,AD转换后的信号经传输单元传输到反射波计算单元和心血管功能评估单元。
4.根据权利要求1所述的基于脉搏波智能分解的可穿戴心血管功能评价系统,所述控制单元采用微控制器,根据压力反馈传感器反馈的桡动脉阻断袖带内实时压力值判断是否达到阻断的要求,从而控制微型气泵是否继续充气,实现精确控制并调节加压压力稳定测量袖带内压力,并在微控制器的I/O口上连接四个三极管元件实现控制信号与供电电压隔离,加速微型气泵的充气速率;同时控制单元通过控制桡动脉血流阻断单元的桡动脉阻断袖带加压时间,先后分时测量同一位置的脉搏波和阻断后的前向波。
5.根据权利要求1所述的基于脉搏波智能分解的可穿戴心血管功能评价系统,所述传输模块包括存储模块和蓝牙收发模块;所述存储模块保存脉搏波与前向波采集单元采集的脉搏波和前向波,并将这些波形数据通过蓝牙收发模块传输到反射波计算单元。
6.根据权利要求1所述的基于脉搏波智能分解的可穿戴心血管功能评价系统,所述反射波计算单元包括第一处理模块、第二处理模块、第三处理模块和第四处理模块;所述第一处理模块读取采集的测量者脉搏波和前向波,并进行波形去噪;所述第二处理模块对脉搏波和前向波信号分别去除其基线漂移,并为了反射波波形的准确计算而将脉搏波和前向波信号对齐在同一水平;所述第三处理模块将脉搏波和前向波按照第一个采集到的有效波谷位置对齐,避免因相位差导致的反射波不准确;所述第四处理模块,按照如下公式求解反射波:
Pb(t)=P(t)-Pf(t)
其中,Pb(t)为反射波波形,P(t)和Pf(t)分别为采集的脉搏波波形和前向波波形,t为时间。
7.根据权利要求1所述的基于脉搏波智能分解的可穿戴心血管功能评价系统,所述参数计算与波形显示模块,通过分解得到的反射波,提取反射波、脉搏波和前向波相关波形特征,计算心血管动脉反射指数RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr和动脉反射波增强指数AIx。
8.一种基于脉搏波智能分解的可穿戴心血管功能评价方法,采用权利要求1所述的基于脉搏波智能分解的可穿戴心血管功能评价系统进行评价,包括以下步骤:
步骤1、调节可调压力测量袖带和桡动脉阻断袖带压力进行脉搏波采集和分解,得到脉搏波和前向波,具体方法为:
步骤1.1、控制单元控制脉搏波与前向波测量模块的可调压力测量袖带加压至70mmHg,脉搏波与前向波测量模块的前向波测量传感器采集70mmHg压力下的脉搏波,即为脉搏波波形,保存并传输至信号调理模块;
步骤1.2、控制单元控制桡动脉血流阻断单元的桡动脉阻断袖带加压至200mmHg时,同时桡动脉血流阻断单元的压力反馈传感器将桡动脉阻断袖带中的压力波形实时反馈给控制单元,当压力波形幅值小于阈值时,停止加压,此时血流在桡动脉血流阻断单元的桡动脉阻断袖带处实现完全阻断,桡动脉的前向波在此处发生完全反射,在脉搏波与前向波测量模块袖带处采集脉搏波前向波的等量叠加;
步骤1.3、脉搏波与前向波测量单元采集桡动脉血流阻断时的脉搏波和前向波等量叠加的压力波形,即两倍前向波,并进行信号调理;
步骤2、反射波计算单元对脉搏波与前向波测量模块采集到的脉搏波和两倍前向波进行处理和计算,得到反射波,具体方法为:
步骤2.1、读取采集的测量者脉搏波和两倍前向波各十个连续周期,并进行波形去噪和去基线漂移处理,并为了反射波波形的准确计算将两路波形信号对齐在同一水平,将两路波形信号各十个周期分别进行叠加并取平均,进一步保证反射波波形精准计算;
步骤2.2、将脉搏波和两倍前向波按照波谷位置水平在时间域对齐,避免因相位差导致的反射波不准确;
步骤2.3、按照如下公式求解反射波:
Pb(t)=P(t)-Pf(t)
其中,Pb(t)为反射波波形,P(t)为采集到的脉搏波波形,Pf(t)为直接采集到的两倍前向波波形除以2,即为脉搏波的前向波;
步骤3、分解得到的反射波,通过参数计算与波形显示模块,提取反射波和脉搏波、前向波相关波形特征,计算心血管相关参数:桡动脉反射指数RI、桡动脉反射强度RM、桡动脉反射波增强指数AIr、动脉反射波增强指数AIx;
步骤4、心血管功能评价模块根据计算得到的心血相关参数对心血管功能和健康程度进行评估和显示。
9.根据权利要求8所述的基于脉搏波智能分解的可穿戴心血管功能评价方法,其特征在于:
所述桡动脉反射指数RI如下公式所示:
其中,Pb为反射波压力,Pf为前向波压力;
所述桡动脉反射强度RM如下公式所示:
其中,Pb为反射波压力,Pf为前向波压力;
所述桡动脉反射波增强指数AIr为是桡动脉反射波和前向波叠加后波形的压力峰值P2与收缩早期即反射波到达前压力峰值P1的比值,如下公式所示:
所述动脉反射波增强指数AIx如下公式所示:
其中,Ps为脉搏波波形收缩期峰值,Pd脉搏波波形舒张期最低值,Pi为反射波波形中的反射波峰值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910509555.1A CN110420014A (zh) | 2019-06-13 | 2019-06-13 | 基于脉搏波智能分解的可穿戴心血管功能评价系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910509555.1A CN110420014A (zh) | 2019-06-13 | 2019-06-13 | 基于脉搏波智能分解的可穿戴心血管功能评价系统及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110420014A true CN110420014A (zh) | 2019-11-08 |
Family
ID=68407606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910509555.1A Pending CN110420014A (zh) | 2019-06-13 | 2019-06-13 | 基于脉搏波智能分解的可穿戴心血管功能评价系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110420014A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110974193A (zh) * | 2019-12-30 | 2020-04-10 | 北京中科芯健医疗科技有限公司 | 一种人体功能状态评估装置 |
CN111358498A (zh) * | 2020-03-23 | 2020-07-03 | 贵州民族大学 | 去除超声脉搏波上冲段反射波和量化噪声的方法及系统 |
CN113017576A (zh) * | 2021-02-03 | 2021-06-25 | 西安宏源视讯设备有限责任公司 | 一种智能脉搏仪 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1991019451A1 (en) * | 1990-06-19 | 1991-12-26 | Mitsuei Tomita | Device for measuring blood flow velocity and quantity in aorta |
US5095912A (en) * | 1989-05-19 | 1992-03-17 | Mitsuei Tomita | Arterial distensibility measuring apparatus |
EP1179318A1 (en) * | 2000-08-11 | 2002-02-13 | Colin Corporation | Central-artery-pressure-waveform estimating apparatus |
EP1402814A1 (en) * | 2002-09-27 | 2004-03-31 | Colin Corporation | Arteriosclerosis-degree evaluating apparatus |
CN1909829A (zh) * | 2004-02-18 | 2007-02-07 | 伊里斯·米科诺斯 | 测量血液动力学参数的装置及方法 |
CN101156771A (zh) * | 2007-09-28 | 2008-04-09 | 天津市先石光学技术有限公司 | 基于脉搏波频谱分析提高血管硬度测量精度的方法及装置 |
CN201533838U (zh) * | 2009-08-11 | 2010-07-28 | 秦皇岛市康泰医学系统有限公司 | 动态血压记录仪 |
CN102462493A (zh) * | 2010-11-19 | 2012-05-23 | 伍霆杰 | 一种血压测量用测量带、装置、其制备方法及用途 |
CN103140165A (zh) * | 2010-09-28 | 2013-06-05 | 欧姆龙健康医疗事业株式会社 | 血压信息测定装置及用该装置的动脉硬化度指标计算方法 |
CN103260503A (zh) * | 2010-12-08 | 2013-08-21 | 欧姆龙健康医疗事业株式会社 | 血压信息测定装置及在该装置中的动脉硬化度的指标的计算方法 |
-
2019
- 2019-06-13 CN CN201910509555.1A patent/CN110420014A/zh active Pending
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5095912A (en) * | 1989-05-19 | 1992-03-17 | Mitsuei Tomita | Arterial distensibility measuring apparatus |
WO1991019451A1 (en) * | 1990-06-19 | 1991-12-26 | Mitsuei Tomita | Device for measuring blood flow velocity and quantity in aorta |
EP1179318A1 (en) * | 2000-08-11 | 2002-02-13 | Colin Corporation | Central-artery-pressure-waveform estimating apparatus |
EP1402814A1 (en) * | 2002-09-27 | 2004-03-31 | Colin Corporation | Arteriosclerosis-degree evaluating apparatus |
CN1909829A (zh) * | 2004-02-18 | 2007-02-07 | 伊里斯·米科诺斯 | 测量血液动力学参数的装置及方法 |
CN101156771A (zh) * | 2007-09-28 | 2008-04-09 | 天津市先石光学技术有限公司 | 基于脉搏波频谱分析提高血管硬度测量精度的方法及装置 |
CN201533838U (zh) * | 2009-08-11 | 2010-07-28 | 秦皇岛市康泰医学系统有限公司 | 动态血压记录仪 |
CN103140165A (zh) * | 2010-09-28 | 2013-06-05 | 欧姆龙健康医疗事业株式会社 | 血压信息测定装置及用该装置的动脉硬化度指标计算方法 |
CN102462493A (zh) * | 2010-11-19 | 2012-05-23 | 伍霆杰 | 一种血压测量用测量带、装置、其制备方法及用途 |
CN103260503A (zh) * | 2010-12-08 | 2013-08-21 | 欧姆龙健康医疗事业株式会社 | 血压信息测定装置及在该装置中的动脉硬化度的指标的计算方法 |
Non-Patent Citations (1)
Title |
---|
王宇: "基于有效反射距离理论及单路脉波的动脉硬化评价", 《中国优秀硕士学位论文全文数据库医药卫生科技辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110974193A (zh) * | 2019-12-30 | 2020-04-10 | 北京中科芯健医疗科技有限公司 | 一种人体功能状态评估装置 |
CN111358498A (zh) * | 2020-03-23 | 2020-07-03 | 贵州民族大学 | 去除超声脉搏波上冲段反射波和量化噪声的方法及系统 |
CN111358498B (zh) * | 2020-03-23 | 2021-08-24 | 贵州民族大学 | 去除超声脉搏波上冲段反射波和量化噪声的方法及系统 |
CN113017576A (zh) * | 2021-02-03 | 2021-06-25 | 西安宏源视讯设备有限责任公司 | 一种智能脉搏仪 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100413463C (zh) | 利用脉搏品质值拒绝赝象的设备 | |
CN110420014A (zh) | 基于脉搏波智能分解的可穿戴心血管功能评价系统及方法 | |
CN102008296B (zh) | 基于脉搏波信号和心电信号测量动脉血压装置及测量方法 | |
US20120157791A1 (en) | Adaptive time domain filtering for improved blood pressure estimation | |
EP0651970A1 (en) | Method and apparatus for assessing cardiovascular performance | |
US20020038090A1 (en) | Central-artery-pressure-waveform estimating apparatus | |
CN110141196B (zh) | 基于双三角血流模型的外周动脉血管弹性评价系统 | |
US20050187483A1 (en) | Oscillometric determination of blood pressure | |
CN111493855B (zh) | 个体化心输出量的无创测量系统与方法 | |
CN108024733B (zh) | 评估内皮功能的系统和方法 | |
CN106388789B (zh) | 一种脉搏波测量装置及方法 | |
EP3644935B1 (en) | Monitoring vital parameters of a compression garment wearer | |
US20130060152A1 (en) | Apparatus and method for continuous oscillometric blood pressure measurement | |
WO2010058169A1 (en) | Apparatus and method for measuring blood pressure | |
CN104042200A (zh) | 一种逐拍动脉血压的无创监测装置及其方法 | |
US7887491B2 (en) | Impedance based device for non-invasive measurement of blood pressure and ankle-brachial index | |
EP3457929B1 (en) | Non-invasive system and method for measuring blood pressure variability | |
CN107960998A (zh) | 血压测量方法以及血压测量装置 | |
CN103767694B (zh) | 一种准确提取袖带压力震荡波的方法 | |
CN101224106A (zh) | 一种人体动脉顺应性的检测方法及装置 | |
CN101810474B (zh) | 一种实现无创测量血压和脉搏波速度的装置 | |
CN210095711U (zh) | 一种无创连续血压测量设备 | |
US8690786B2 (en) | System and method for a non-invasive blood pressure monitor | |
US9974449B2 (en) | Method and apparatus for arterial blood pressure measurement and individualized rectifying technology | |
WO2019222923A1 (zh) | 一种脉象仪及脉象仪系统 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191108 |