CN111493855B - 个体化心输出量的无创测量***与方法 - Google Patents

个体化心输出量的无创测量***与方法 Download PDF

Info

Publication number
CN111493855B
CN111493855B CN202010317299.9A CN202010317299A CN111493855B CN 111493855 B CN111493855 B CN 111493855B CN 202010317299 A CN202010317299 A CN 202010317299A CN 111493855 B CN111493855 B CN 111493855B
Authority
CN
China
Prior art keywords
pulse wave
artery
module
vector
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
Application number
CN202010317299.9A
Other languages
English (en)
Other versions
CN111493855A (zh
Inventor
肖汉光
黄金锋
任慧娇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chongqing University of Technology
Original Assignee
Chongqing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Technology filed Critical Chongqing University of Technology
Priority to CN202010317299.9A priority Critical patent/CN111493855B/zh
Publication of CN111493855A publication Critical patent/CN111493855A/zh
Application granted granted Critical
Publication of CN111493855B publication Critical patent/CN111493855B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/026—Measuring blood flow
    • A61B5/029—Measuring or recording blood output from the heart, e.g. minute volume
    • 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Physiology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Biophysics (AREA)
  • Hematology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Vascular Medicine (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明公开了个体化心输出量的无创测量***与方法，包括肱动脉血压测量模块、外周浅层动脉脉搏波测量模块与脉搏波传导速度测量模块；脉搏波校正模块根据肱动脉血压参数对外周浅层动脉脉搏波进行校正；脉搏波特征提取模块用于提取脉搏波特征并输出给个体化脉搏特征向量生成模块，生成的个体化脉搏波波特征向量包括脉搏波特征与受测者的个体生理参数；人工神经网络根据个体化脉搏波特征向量识别出相应的每博量，心输出量计算模块，根据每博量计算出心输出量。本发明解决了现有技术中的无创测量方法针对个体的心输出量测量的准确性不高的技术问题，大大缩小了无创测量与金标准法测量精度之间的差距。

Description

个体化心输出量的无创测量***与方法
技术领域
本发明属于心输出量测量技术领域，尤其涉及无创心输出量测量技术领域。
背景技术
心输出量(CO)是指左心室每分钟的射血量，是表征心血管***健康状态最重要的参数，是心脏功能及心血管疾病的重要诊断依据。另外，基于心输出量可以辅助计算出许多其他心血管***参数，因此，心输出量的准确测量在心血管疾病检测和治疗等方面十分关键，具有重要临床意义。
目前，测量心输出量的装置及方法分为有创、微创和无创三大类。
有创的装置及方法是利用介入设备将导管从外周动脉***，通过导管在右心房上部注入一定量生理盐水，与血液混合，同时记录温度变化曲线，通过温度曲线的计算得到CO。该方法称为热稀释法，是测量心输出量的“金标准”。有创法主要用于急救、心血管介入手术和重症监护等领域，具有创伤性大、操作复杂等特点，并且长时间测量容易引起感染或并发症。所以，该方法临床上具有一定的局限性。
微创的装置及方法结合单次热稀释法和脉搏波进行CO测量，该技术称为PICCO技术。FloTrac/Vigileo***是典型的微创装置，它利用FloTrac传感器采集患者股动脉或腋动脉压力波形，通过对动脉压力波形曲线下的面积进行分析，结合患者年龄、性别、身高、体重、体表面积等参数，计算相对每搏量，并且利用单次热稀释法定标，实现心输出量连续测量。虽然微创法的有创性大幅减小，并且大幅降低了感染和并发症的风险，但是，仍然具有创性，无法满足非危重病人和健康人对心输出量测量的需求，因此，一些无创心输出量测量方法被相继提出。
无创的装置及方法主要包括超声法、胸阻抗法、脉搏波波形分析法等。
超声法通过连续多普勒超声波技术，测量主动脉或者肺动脉的血流速度，得到血流速度的时间积分，再乘其管腔截面面积计算出每搏量等指标。代表性仪器有澳大利亚USCOM仪器。基于超声法的CO测量仪器对操作者要求较高，传感器的力度和方向都会影响血流速度波形，从而影响计算结果，超声法与热稀释法的测量结果误差可达40％，所以，其准确性和稳定性有待改进。
胸阻抗法通过测量人体体表电极的电位变化，建立电位与心排量之间的关系，从而实现心排量的测量。该方法的典型代表包括美国的NICOM***和德国的ICON***等。阻抗法因其无创、连续测量、设备简单、安全已普遍应用与临床，但是该技术较依赖于振荡电流通过胸腔的扩散，抗干扰能力较差。
基于脉搏波波形分析法可以通过光学技术获得动脉内血流容量变化引起的透射或反射光的变化。当动脉中血流量随心脏搏动变化时，光源发出的光经过动脉后被探测器接收，其接收光强与血流量同步发生脉动变化，进而利用动脉血流量与CO之间的关系模型计算出CO，从而实现CO测量。光电容积脉搏波测量心输出量时成本较低、信号稳定，但波形形状易受传感器位置和背景光的影响，准确性有待提高。
基于脉搏波波形分析法也可以通过压力传感器获得的外周脉搏波形，然后建立血管数学模型和统计模型，如Windkessel模型、传输线模型、回归模型等，通过模型计算CO。由于外周脉搏波形获取十分容易，且设备廉价，能很好地监测心输出量的相对变化，所以，该技术具有宽广的应用前景，但是，该技术还存在一些不足，如数学模型不够准确，尤其是缺乏对个体差异性的考虑，往往导致心输出量绝对值偏差较大。
综上所述，心输出量的测量正从有创和微创向无创方向发展，基于脉搏波的波形分析法的优势较为突出，但是存在心输出量测量不够准确的问题。
发明内容
针对上述技术的不足，本发明提供了一种个体化心输出量的无创测量方法，解决现有技术中的无创测量方法针对个体的心输出量测量的准确性不高的技术问题。
为解决上述技术问题，本发明的技术方案如下：一种个体化心输出量的无创测量***，包括肱动脉血压测量模块与外周浅层动脉脉搏波测量模块，肱动脉测量血压测量模块用于采集肱动脉的血压参数，包括肱动脉收缩压与肱动脉舒张压；外周浅层动脉脉搏波测量模块用于采集外周浅层动脉脉搏波；
脉搏波校正模块用于输入肱动脉的血压参数与外周浅层动脉脉搏波，利用肱动脉的血压参数对外周浅层动脉脉搏波进行校正，并将校正后的外周浅层动脉脉搏波输出给脉搏波特征提取模块；
脉搏波特征提取模块用于从校正后的外周浅层动脉脉搏波中提取脉搏波特征，并包括以下时域特征：血压特征、时间特征、面积特征与比例特征，并将脉搏波特征输出给个体化脉搏波特征向量生成模块；
个体化脉搏波特征向量生成模块用于根据脉搏波特征与受测者的个体生理参数生成个体化脉搏波特征向量，并输出给已训练完成的人工神经网络；
人工神经网络用于根据输入的个体化脉搏波特征向量，计算出对应的每博量，并输出给心输出量计算模块；
心输出量计算模块用于根据每博量计算出心输出量，计算公式如下：
CO＝SV*N/Fs*60
其中，CO表示心输出量，SV表示每博量，N表示一个心动周期内脉搏波的采样点数，Fs表示***的信号采样频率。
进一步的，还包括用于测量从肱动脉到外周浅层动脉的脉搏波传导速度的脉搏波传导速度测量模块；外周浅层动脉脉搏波为桡动脉血压波、手指动脉血压波或下肢动脉血压波，脉搏波校正模块根据如下公式对外周浅层动脉脉搏波进行校正：
Figure BDA0002459919350000031
其中，
Figure BDA0002459919350000032
是校正前的外周浅层动脉脉搏波；Wr是校正后的外周浅层动脉脉搏波；SBPb是肱动脉收缩压；DBPb是肱动脉舒张压；SBPr是外周浅层动脉脉搏波收缩压；DBPr是外周浅层动脉脉搏波舒张压；α是肱动脉到外周浅层动脉的放大因子，
Figure BDA0002459919350000033
c1、c2均为回归系数，PWV为受测者的脉搏波传导速度，PWV0为脉搏波传导速度参考值。
进一步的，放大因子α的取值范围限制在[1,1.5]。
进一步的，人工神经网络的输入层采用如下公式对个体化脉搏波特征向量进行归一化：
Figure BDA0002459919350000034
其中，
Figure BDA0002459919350000035
表示个体化脉搏波特征向量，
Figure BDA0002459919350000036
表示归一化后的个体化脉搏波特征向量，
Figure BDA0002459919350000037
表示归一化增益向量，
Figure BDA0002459919350000038
表示归一化偏移向量；
归一化后的个体化脉搏波特征向量
Figure BDA0002459919350000039
经过隐藏层与输出层后得到准输出向量
Figure BDA00024599193500000310
准输出向量
Figure BDA00024599193500000311
的计算公式如下：
Figure BDA0002459919350000041
其中，f表示人工神经网络中的结合函数，
Figure BDA0002459919350000042
表示
Figure BDA0002459919350000043
的转置向量，
Figure BDA0002459919350000044
表示权重矩阵，
Figure BDA0002459919350000045
表示权重列向量，
Figure BDA0002459919350000046
与
Figure BDA0002459919350000047
均表示偏移列向量；
输出层将准输出向量
Figure BDA0002459919350000048
进行反归一化，从而得到输出量，即每博量，按如下公式计算：
Figure BDA0002459919350000049
其中，SV表示每博量，
Figure BDA00024599193500000410
表示反归一化增益向量，
Figure BDA00024599193500000411
表示反归一化偏移向量。
进一步的，用于训练人工神经网络的样本数据来源于热稀释法的心输出量、每博量和桡动脉血压波形，每个样本均以饶动脉血压波形的时频特征和个体生理参数组成的特征向量作为输入量，并以与桡动脉血压波形对应的每博量作为标准输出量；样本数据被分成训练样本集与测试样本集，采用训练样本集对人工神经网络进行训练；采用测试样本集对训练后的人工神经网络进行测试，若输出量与标准输出量的误差满足阈值，则表示训练完成，训练完成后的人工神经网络获得能够使输出量与标准输出量的误差满足阈值的参数组合，参数组合中的参数包括归一化增益向量
Figure BDA00024599193500000412
归一化偏移向量
Figure BDA00024599193500000413
权重矩阵
Figure BDA00024599193500000414
权重列向量
Figure BDA00024599193500000415
偏移列向量
Figure BDA00024599193500000416
偏移列向量
Figure BDA00024599193500000417
反归一化增益向量
Figure BDA00024599193500000418
与反归一化偏移向量
Figure BDA00024599193500000419
本发明还提供一种个体化心输出量的无创测量方法，包括以下步骤：
步骤1：获取受测者的肱动脉的血压参数，包括肱动脉收缩压与肱动脉舒张压；获取受测者的外周浅层动脉脉搏波，外周浅层动脉脉搏波为桡动脉血压波、手指动脉血压波或下肢动脉血压波；获取受测者从肱动脉到外周浅层动脉的脉搏波传导速度；
步骤2：利用肱动脉的血压参数对外周浅层动脉脉搏波进行校正，按如下公式：
Figure BDA00024599193500000420
其中，
Figure BDA00024599193500000421
是校正前的外周浅层动脉脉搏波；Wr是校正后的外周浅层动脉脉搏波；SBPb是肱动脉收缩压；DBPb是肱动脉舒张压；SBPr是外周浅层动脉脉搏波收缩压；DBPr是外周浅层动脉脉搏波舒张压；α是肱动脉到外周浅层动脉的放大因子，
Figure BDA00024599193500000422
c1、c2均为回归系数，PWV为受测者的脉搏波传导速度，PWV0为脉搏波传导速度参考值；
步骤3：从校正后的外周浅层动脉脉搏波中提取脉搏波特征，包括时域特征与频域特征；所述时域特征包括血压特征、时间特征、面积特征与比例特征，所述频域特征包括多倍频的谐波幅值；时域特征与频域特征合并组成脉搏波特征；
步骤4：根据脉搏波特征与受测者的个体生理参数生成个体化脉搏波特征向量；
步骤5：将个体化脉搏波特征向量输出给已训练完成的人工神经网络，人工神经网络根据个体化脉搏波特征向量计算出对应的每搏量；
步骤6：根据每博量计算出心输出量，按如下公式：
CO＝SV*N/Fs*60
其中，CO表示心输出量，SV表示每博量，N表示一个心动周期内脉搏波的采样点数，Fs表示***的信号采样频率。
与现有技术相比，本发明具有的优点包括：
1、脉搏波较微弱且容易受到干扰，本发明通过肱动脉血压参数对外周浅层动脉脉搏波进行校正，以降低失真度，是准确提取脉搏波特征的基础。与现有技术采用固定比例因子进行校正所不同的是，本发明在对脉搏波进行校正时，基于动脉硬化与脉搏波放大存在正比例关系，先利用
Figure BDA0002459919350000051
将外周浅层动脉的收缩压与舒张压映射到肱动脉的收缩压与舒张压，再乘以肱动脉到外周浅层动脉的放大因子α，体现放大效应。
2、本发明的放大因子α基于个体的脉搏波传导速度进行计算，使得校正后的脉搏波体现出个体差异性。放大因子α的取值范围限制在[1,1.5]，防止因测量脉搏波传导速度不准确导致α过大的情况发生。同时考虑到个体生理差异对脉搏波传导过程影响的程度不同，在建立脉搏波特征向量时引入个体生理参数(如身高、年龄、体重、性别、病史等)，避免片面的关注波形特征，从而形成多参数融合的个体化脉搏波特征向量。
3、本发明结合人工神经网络对个体化脉搏波特征向量进行每博量识别，避免复杂的数学建模，同时也克服了数学模型表达能力不足的缺陷。训练样本源于热稀释法的测量数据，保证了人工神经网络的准确性，使训练后的人工神经网络能够输出更加接近于热稀释法的测量结果。
4、通过对输入向量进行归一化将向量中的各种参数的数量级调整为一致，便于人工神经网络的训练。反归一化将准输出向量调整至实际范围。
附图说明
图1为个体化心输出量的无创测量***的结构框图；
图2为采用现有技术中的测量模型测量的心输出量与“金标准”方法测量结果的对比图；
图3为采用本具体实施方式中的个体化心输出量的无创测量***测量的心输出量与“金标准”方法测量结果的对比图；
图4为校正前的桡动脉脉搏波的波形图。
具体实施方式
一)***结构
参考图1所示，一种个体化心输出量的无创测量***，包括肱动脉血压测量模块、外周浅层动脉脉搏波测量模块与脉搏波传导速度测量模块，肱动脉测量血压测量模块用于采集肱动脉的血压参数，包括肱动脉收缩压与肱动脉舒张压；外周浅层动脉脉搏波测量模块用于采集外周浅层动脉脉搏波；脉搏波传导速度测量模块用于从肱动脉到外周浅层动脉的脉搏波传导速度；
脉搏波校正模块用于输入肱动脉的血压参数与外周浅层动脉脉搏波，利用肱动脉的血压参数对外周浅层动脉脉搏波进行校正，并将校正后的外周浅层动脉脉搏波输出给脉搏波特征提取模块；
脉搏波特征提取模块用于从校正后的外周浅层动脉脉搏波中提取脉搏波特征，并包括以下时域特征：血压特征、时间特征、面积特征与比例特征，并将脉搏波特征输出给个体化脉搏波特征向量生成模块；
个体化脉搏波特征向量生成模块用于根据脉搏波特征与受测者的个体生理参数生成个体化脉搏波特征向量，并输出给已训练完成的人工神经网络；
人工神经网络用于根据输入的个体化脉搏波特征向量，计算出对应的每博量，并输出给心输出量计算模块；
心输出量计算模块用于根据每博量计算出心输出量，计算公式如下：
CO＝SV*N/Fs*60
其中，CO表示心输出量，SV表示每博量，N表示一个心动周期内脉搏波的采样点数，Fs表示***的信号采样频率。
外周浅层动脉脉搏波为桡动脉血压波、手指动脉血压波或下肢动脉血压波，脉搏波校正模块根据如下公式对外周浅层动脉脉搏波进行校正：
Figure BDA0002459919350000071
其中，
Figure BDA0002459919350000072
是校正前的外周浅层动脉脉搏波；Wr是校正后的外周浅层动脉脉搏波；SBPb是肱动脉收缩压；DBPb是肱动脉舒张压；SBPr是外周浅层动脉脉搏波收缩压；DBPr是外周浅层动脉脉搏波舒张压；α是肱动脉到外周浅层动脉的放大因子，
Figure BDA0002459919350000073
c1、c2均为回归系数(通过测量多个人的PWV和桡动脉相对肱动脉血压的放大倍数，采用回归分析得到回归系数)，PWV为受测者的脉搏波传导速度，PWV0为脉搏波传导速度参考值。
放大因子α的取值范围限制在[1,1.5]。
为了建立能更加准确反映脉搏波的特征向量，脉搏波特征提取模块还从校正后的外周浅层动脉脉搏波中提取包含多倍频的谐波幅值在内的频域特征，时域特征与频域特征合并组成脉搏波特征后输出给个体化脉搏波特征向量生成模块。
人工神经网络的输入层采用如下公式对个体化脉搏波特征向量进行归一化：
Figure BDA0002459919350000074
其中，
Figure BDA0002459919350000075
表示个体化脉搏波特征向量，
Figure BDA0002459919350000076
表示归一化后的个体化脉搏波特征向量，
Figure BDA0002459919350000077
表示归一化增益向量，
Figure BDA0002459919350000078
表示归一化偏移向量；
归一化后的个体化脉搏波特征向量
Figure BDA0002459919350000079
经过隐藏层与输出层后得到准输出向量
Figure BDA00024599193500000710
准输出向量
Figure BDA00024599193500000711
的计算公式如下：
Figure BDA00024599193500000712
其中，f表示人工神经网络中的结合函数，
Figure BDA00024599193500000713
表示
Figure BDA00024599193500000714
的转置向量，
Figure BDA00024599193500000715
表示权重矩阵，
Figure BDA0002459919350000081
表示权重列向量，
Figure BDA0002459919350000082
与
Figure BDA0002459919350000083
均表示偏移列向量；
输出层将准输出向量
Figure BDA0002459919350000084
进行反归一化，从而得到输出量，即每博量，按如下公式计算：
Figure BDA0002459919350000085
其中，SV表示每博量，
Figure BDA0002459919350000086
表示反归一化增益向量，
Figure BDA0002459919350000087
表示反归一化偏移向量。
还包括数据显示模块与数据采集控制模块，所述数据采集控制模块用于控制肱动脉血压测量模块、外周浅层动脉脉搏波测量模块与脉搏波传导速度测量模块的数据采集过程；数据显示模块用于显示心输出量测量过程中的相关数据，包括肱动脉血压参数、外周浅层动脉脉搏波、校正后的外周浅层动脉脉搏波、每博量、心输出量与受测者的个体生理参数。
肱动脉血压测量模块包括袖带、气路导管、气压传感器、充放气电机与血压数据处理模块；电机通过气路导管对袖带进行充放气，气压传感器监测袖带或气路中的气压变化，并将测量的气压信号传递给血压数据处理模块；血压数据处理模块按如下程序计算肱动脉血压参数：对气压信号进行分解，得到振荡信号和线性上升或下降信号；检测出振荡信号的峰值，在峰值左右按峰值幅度的一定比例(如60％或80％)获得对应两个时间点，左右两个时间点对应的上升或下降信号的压力，即为收缩压和舒张压；
外周浅层动脉脉搏波测量模块为桡动脉测量装置、手指动脉测量装置或下肢动脉测量装置。本具体实施方式采用桡动脉血压测量模块：由压电传感器、一路导线和血压信号处理电路组成，其中，压电传感器获得桡动脉的脉动信号，并将信号通过导线传输给血压信号处理电路，该电路完成血压信号的滤波和放大等处理，并将处理后的信号通过数据采集控制模块传递给血压数据处理模块。
脉搏波传导速度测量模块采用肱动脉血压测量模块和外周浅层动脉脉搏波测量模块同步测量肱动脉脉搏波和外周浅层动脉脉搏波，并计算两路波形的波足点时间差PTT(脉搏波传导时间)，然后采用卷尺测量两测量点之间的人体体表距离L，最后利用L除以PTT即为脉搏波传导速度。一般将袖带压升至肱动脉平均压，不能升得过高，这会导致远端外周动脉脉搏波消失，从而使测量失败。
二)、训练人工神经网络
采用误差反向传播法结合梯度下降法对人工神经网络进行训练。用于训练人工神经网络的样本数据来源于热稀释法的心输出量、每博量和桡动脉血压波形，这些数据通过公用数据就可以查询到。每个样本均以饶动脉血压波形的时频特征和个体生理参数组成的特征向量作为输入量，并以与桡动脉血压波形对应的每博量作为标准输出量；样本数据被分成训练样本集与测试样本集，采用训练样本集对人工神经网络进行训练；采用测试样本集对训练后的人工神经网络进行测试，若输出量与标准输出量的误差满足阈值，则表示训练完成，训练完成后的人工神经网络获得能够使输出量与标准输出量的误差满足阈值的参数组合，参数组合中的参数包括归一化增益向量
Figure BDA0002459919350000091
归一化偏移向量
Figure BDA0002459919350000092
权重矩阵
Figure BDA0002459919350000093
权重列向量
Figure BDA0002459919350000094
偏移列向量
Figure BDA0002459919350000095
偏移列向量
Figure BDA0002459919350000096
反归一化增益向量
Figure BDA0002459919350000097
与反归一化偏移向量
Figure BDA0002459919350000098
受测者的部分波形特征参数统计如下表1：
表1病人的部分波形特征参数的统计值
Figure BDA0002459919350000099
图2现有模型ModelRcdecay与“金标准”法的测量结果对比图，横坐标是金标准测量的SV，纵坐标是现有模型ModelRCdecay计算的SV，图中直线表示线性回归直线。
图3为本发明的人工神经网络(ANN)与“金标准”法的测量结果对比图，横坐标是金标准测量的SV，本发明的人工神经网络计算的SV，图中直线表示线性回归直线。
本发明提出人工神经网络(ANN)与其他现有模型的结果比较，如下表2所示，斜率和截距为模型与实测SV的线性回归方程参数，相关系数采用皮尔森相关系数。
表2
Figure BDA00024599193500000910
从上表可以看出，本发明与金标准法的测量结果的标准差远远小于现有模型与金标准法的测量结果之间的标准差，本发明的人工神经网络模型大大提高了无创测量的准确性。
三)、测量过程
采用本具体实施方式中的个体化心输出量的无创测量***对受测者进行检查心输出量测量，具体过程如下。
受测者平躺休息3-5分钟后，通过数据采集控制模块启动和控制肱动脉血压测量模块中气泵的充放气速度和袖带中的气压，完成肱动脉血压的测量。
通过数据采集控制模块启动和控制桡动脉血压测量模块，完成桡动脉血压信号的测量，并通过模拟电路滤波和放大等处理，并将处理后的信号通过数据采集控制模块传递给数据处理模块。
通过数据采集控制模块启动数据处理模块，对其输入的肱动脉气压信号进行处理得到收缩压、平均压和舒张压；对桡动脉血压信号数据进行去噪、去漂移、归一化和校准等处理；然后利用这三个血压对桡动脉血压波形进行校正和一定增益的放大，将校正并放大后的桡动脉脉搏波信号分别发送至脉搏波特征提取模块和显示模块。
在脉搏波特征提取模块中，首先，通过脉搏波的二阶或以上高阶数值差分，检测差分数列的极小值和极大值，从而获得脉搏波下降期的二重波的降中峡和波峰，其特征点参考图4所示。基于这两个特征点和收缩期峰值特征点，计算其他相关的时域特征。
频域特征通过脉搏波的傅里叶变换后获得多倍频的谐波幅值，然后，合并时域和频域特征构建脉搏波特征向量。启动脉搏波特征计算模块，完成脉搏波时域和频域特征的处理和计算，构建脉搏波特征向量，将脉搏波特征向量与受测者的个体生理参数输入个体化脉搏波特征向量生成模块，生成个体化的脉搏波特征向量，并输入至人工神经网络，部分参数发送给数据显示模块显示。
人工神经网络对输入的个体化脉搏波特性向量进行输入归一化、变换和输出反归一化等一系列操作，计算出每搏量并输出给心输出量计算模块，心输出量计算模块根据每博量计算出心输出量。
数据显示模块脉搏波波形、每搏量、心输出量和相关心血管***参数显示在显示屏上，同时生成测量报告。

Claims (9)

1.一种个体化心输出量的无创测量***，其特征在于，包括肱动脉血压测量模块、外周浅层动脉脉搏波测量模块，肱动脉测量血压测量模块用于采集肱动脉的血压参数，包括肱动脉收缩压与肱动脉舒张压；外周浅层动脉脉搏波测量模块用于采集外周浅层动脉脉搏波；
脉搏波校正模块用于输入肱动脉的血压参数与外周浅层动脉脉搏波，利用肱动脉的血压参数对外周浅层动脉脉搏波进行校正，并将校正后的外周浅层动脉脉搏波输出给脉搏波特征提取模块；
脉搏波特征提取模块用于从校正后的外周浅层动脉脉搏波中提取脉搏波特征，并包括以下时域特征：血压特征、时间特征、面积特征与比例特征，并将脉搏波特征输出给个体化脉搏波特征向量生成模块；
个体化脉搏波特征向量生成模块用于根据脉搏波特征与受测者的个体生理参数生成个体化脉搏波特征向量，并输出给已训练完成的人工神经网络；
人工神经网络用于根据输入的个体化脉搏波特征向量，计算出对应的每博量，并输出给心输出量计算模块；
心输出量计算模块用于根据每博量计算出心输出量，计算公式如下：
CO＝SV*N/Fs*60
其中，CO表示心输出量，SV表示每博量，N表示一个心动周期内脉搏波的采样点数，Fs表示***的信号采样频率。
2.根据权利要求1所述的个体化心输出量的无创测量***，其特征在于，还包括用于测量肱动脉到外周浅层动脉的脉搏波传导速度的脉搏波传导速度测量模块；外周浅层动脉脉搏波为桡动脉血压波、手指动脉血压波或下肢动脉血压波，脉搏波校正模块根据如下公式对外周浅层动脉脉搏波进行校正：
Figure FDA0003925480650000011
其中，
Figure FDA0003925480650000012
是校正前的外周浅层动脉脉搏波；Wr是校正后的外周浅层动脉脉搏波；SBPb是肱动脉收缩压；DBPb是肱动脉舒张压；SBPr是外周浅层动脉脉搏波收缩压；DBPr是外周浅层动脉脉搏波舒张压；α是肱动脉到外周浅层动脉的放大因子，
Figure FDA0003925480650000013
c1、c2均为回归系数，PWV为受测者的脉搏波传导速度，PWV0为脉搏波传导速度参考值。
3.根据权利要求2所述的个体化心输出量的无创测量***，其特征在于，放大因子α的取值范围限制在[1,1.5]。
4.根据权利要求1所述的个体化心输出量的无创测量***，其特征在于，脉搏波特征提取模块还从校正后的外周浅层动脉脉搏波中提取包含多倍频的谐波幅值在内的频域特征，时域特征与频域特征合并组成脉搏波特征后输出给个体化脉搏波特征向量生成模块。
5.根据权利要求1所述的个体化心输出量的无创测量***，其特征在于，人工神经网络的输入层采用如下公式对个体化脉搏波特征向量进行归一化：
Figure FDA0003925480650000021
其中，
Figure FDA0003925480650000022
表示个体化脉搏波特征向量，
Figure FDA0003925480650000023
表示归一化后的个体化脉搏波特征向量，
Figure FDA0003925480650000024
表示归一化增益向量，
Figure FDA0003925480650000025
表示归一化偏移向量；
归一化后的个体化脉搏波特征向量
Figure FDA0003925480650000026
经过隐藏层与输出层后得到准输出向量
Figure FDA0003925480650000027
准输出向量
Figure FDA0003925480650000028
的计算公式如下：
Figure FDA0003925480650000029
其中，f表示人工神经网络中的结合函数，
Figure FDA00039254806500000210
表示
Figure FDA00039254806500000211
的转置向量，
Figure FDA00039254806500000212
表示权重矩阵，
Figure FDA00039254806500000213
表示权重列向量，
Figure FDA00039254806500000214
与
Figure FDA00039254806500000215
均表示偏移列向量；
输出层将准输出向量
Figure FDA00039254806500000216
进行反归一化，从而得到输出量，即每博量，按如下公式计算：
Figure FDA00039254806500000217
其中，SV表示每博量，
Figure FDA00039254806500000218
表示反归一化增益向量，
Figure FDA00039254806500000219
表示反归一化偏移向量。
6.根据权利要求1所述的个体化心输出量的无创测量***，其特征在于，用于训练人工神经网络的样本数据来源于热稀释法的心输出量、每博量和桡动脉血压波形，每个样本均以饶动脉血压波形的时频特征和个体生理参数组成的特征向量作为输入量，并以与桡动脉血压波形对应的每博量作为标准输出量；样本数据被分成训练样本集与测试样本集，采用训练样本集对人工神经网络进行训练；采用测试样本集对训练后的人工神经网络进行测试，若输出量与标准输出量的误差满足阈值，则表示训练完成，训练完成后的人工神经网络获得能够使输出量与标准输出量的误差满足阈值的参数组合，参数组合中的参数包括归一化增益向量
Figure FDA00039254806500000220
归一化偏移向量
Figure FDA00039254806500000221
权重矩阵
Figure FDA00039254806500000222
权重列向量
Figure FDA00039254806500000223
偏移列向量
Figure FDA00039254806500000224
偏移列向量
Figure FDA0003925480650000031
反归一化增益向量
Figure FDA0003925480650000032
与反归一化偏移向量
Figure FDA0003925480650000033
7.根据权利要求6所述的个体化心输出量的无创测量***，其特征在于，采用误差反向传播法结合梯度下降法对人工神经网络进行训练。
8.根据权利要求1所述的个体化心输出量的无创测量***，其特征在于，肱动脉血压测量模块包括袖带、气路导管、气压传感器、充放气电机与血压数据处理模块；电机通过气路导管对袖带进行充放气，气压传感器监测袖带或气路中的气压变化，并将测量的气压信号传递给血压数据处理模块；血压数据处理模块按如下程序计算肱动脉血压参数：对气压信号进行分解，得到振荡信号和线性上升或下降信号；检测出振荡信号的峰值，在峰值左右按峰值幅度的一定比例获得对应两个时间点，左右两个时间点对应的上升或下降信号的压力，即为收缩压和舒张压；
外周浅层动脉脉搏波测量模块为桡动脉测量装置、手指动脉测量装置或下肢动脉测量装置。
9.根据权利要求1所述的个体化心输出量的无创测量***，其特征在于，还包括数据显示模块与数据采集控制模块，所述数据采集控制模块用于控制肱动脉血压测量模块与外周浅层动脉脉搏波测量模块的数据采集过程；数据显示模块用于显示心输出量测量过程中的相关数据，包括肱动脉血压参数、外周浅层动脉脉搏波、校正后的外周浅层动脉脉搏波、每博量、心输出量与受测者的个体生理参数。
CN202010317299.9A 2020-04-21 2020-04-21 个体化心输出量的无创测量***与方法 Active CN111493855B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010317299.9A CN111493855B (zh) 2020-04-21 2020-04-21 个体化心输出量的无创测量***与方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010317299.9A CN111493855B (zh) 2020-04-21 2020-04-21 个体化心输出量的无创测量***与方法

Publications (2)

Publication Number Publication Date
CN111493855A CN111493855A (zh) 2020-08-07
CN111493855B true CN111493855B (zh) 2023-01-06

Family

ID=71865898

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010317299.9A Active CN111493855B (zh) 2020-04-21 2020-04-21 个体化心输出量的无创测量***与方法

Country Status (1)

Country Link
CN (1) CN111493855B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112790748A (zh) * 2020-12-30 2021-05-14 重庆理工大学 一种中心动脉压波形重构***及方法
CN113080907B (zh) * 2021-04-14 2022-10-25 贵州省人民医院 脉搏波信号处理方法及装置
CN113057617B (zh) * 2021-04-30 2022-08-26 重庆理工大学 一种心输出量的无创监测***
CN113143230B (zh) * 2021-05-11 2022-05-20 重庆理工大学 一种外周动脉血压波形重构***
CN113499048B (zh) * 2021-07-22 2022-07-08 重庆理工大学 一种基于CNN-BiLSTM的中心动脉压波形重构的***及方法
CN114027804A (zh) * 2021-12-03 2022-02-11 科思技术(温州)研究院 一种脉象诊断的方法、装置及可读存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5400793A (en) * 1991-01-29 1995-03-28 Nederlandse Organisatie Voor Toegepast-Natuurwetenschappelijk Onderzoek Tno Method of determining the stroke volume and the cardiac output of the human heart
CN103892818A (zh) * 2012-12-28 2014-07-02 吴健康 一种无创中心动脉血压测量方法和设备
CN104244814A (zh) * 2012-05-15 2014-12-24 皇家飞利浦有限公司 心输出量的监测
CN104323768A (zh) * 2012-11-20 2015-02-04 深圳市理邦精密仪器股份有限公司 一种心输出量连续监测的参数校准方法

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5913826A (en) * 1996-06-12 1999-06-22 K-One Technologies Wideband external pulse cardiac monitor
US7822470B2 (en) * 2001-10-11 2010-10-26 Osypka Medical Gmbh Method for determining the left-ventricular ejection time TLVE of a heart of a subject
US7815578B2 (en) * 2003-02-10 2010-10-19 Massachusetts Institute Of Technology Methods and apparatus for determining cardiac output
US20070197924A1 (en) * 2004-03-05 2007-08-23 O'rourke Michael F Method and apparatus for determination of cardiac output from the arterial pressure pulse waveform
CN101006919A (zh) * 2007-01-26 2007-08-01 北京工业大学 一种用于大压差下的心输出量检测方法及其装置
CN101176663B (zh) * 2007-12-06 2010-10-27 山东大学 基于脉搏波的无创心输出量检测方法和装置
US10405762B2 (en) * 2009-04-22 2019-09-10 Vital Metrix, Inc. System and method for noninvasively measuring ventricular stroke volume and cardiac output
CN102499669B (zh) * 2011-10-26 2014-12-24 中国科学院深圳先进技术研究院 心脏参数的测量装置
JP5985355B2 (ja) * 2012-10-30 2016-09-06 日本光電工業株式会社 血液量測定方法および測定装置
US9949696B2 (en) * 2013-03-14 2018-04-24 Tensys Medical, Inc. Apparatus and methods for computing cardiac output of a living subject via applanation tonometry
US10542961B2 (en) * 2015-06-15 2020-01-28 The Research Foundation For The State University Of New York System and method for infrasonic cardiac monitoring
CN105726000B (zh) * 2016-01-29 2019-03-22 北京工业大学 一种基于四肢血压脉搏的心脏血管功能参数的装置
WO2018017542A1 (en) * 2016-07-19 2018-01-25 Mayo Foundation For Medical Education And Research Non-invasive cardiac output assessment
CN109512412B (zh) * 2018-11-21 2021-10-08 重庆理工大学 一种中心动脉血压测量装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5400793A (en) * 1991-01-29 1995-03-28 Nederlandse Organisatie Voor Toegepast-Natuurwetenschappelijk Onderzoek Tno Method of determining the stroke volume and the cardiac output of the human heart
CN104244814A (zh) * 2012-05-15 2014-12-24 皇家飞利浦有限公司 心输出量的监测
CN104323768A (zh) * 2012-11-20 2015-02-04 深圳市理邦精密仪器股份有限公司 一种心输出量连续监测的参数校准方法
CN103892818A (zh) * 2012-12-28 2014-07-02 吴健康 一种无创中心动脉血压测量方法和设备

Also Published As

Publication number Publication date
CN111493855A (zh) 2020-08-07

Similar Documents

Publication Publication Date Title
CN111493855B (zh) 个体化心输出量的无创测量***与方法
US5265615A (en) Method and apparatus for continuous measurement of cardiac output and SVR
US7220230B2 (en) Pressure-based system and method for determining cardiac stroke volume
AU704699B2 (en) Method and device for determining brachial arterial pressure wave on the basis of noninvasively measured finger blood pressure wave
US4873987A (en) Noninvasive continuous monitor of arterial blood pressure waveform
US7361148B2 (en) Cuff volumetric pulse wave obtaining apparatus, cuff volumetric pulse wave analyzing apparatus, pressure pulse wave obtaining apparatus, and pressure pulse wave analyzing apparatus
US6120459A (en) Method and device for arterial blood pressure measurement
CN109512412B (zh) 一种中心动脉血压测量装置
CN101765398B (zh) 前负荷依赖性和流体反应性的评估
US20020028998A1 (en) Method and instrument to measure vascular impedance
EP0651970A1 (en) Method and apparatus for assessing cardiovascular performance
US20060247543A1 (en) High resoution bio-impedance device
EP0818175A1 (en) Living body condition measuring apparatus
CN114145724A (zh) 基于ecg和ppg多生理特征参数动态监测血压的方法
WO1992006633A1 (en) Method and apparatus for measuring cardiac output
EP1150604B1 (en) Method and device for continuous analysis of cardiovascular activity of a subject
WO2008156377A1 (en) Method and apparatus for obtaining electronic oscillotory pressure signals from an inflatable blood pressure cuff
Zhang et al. A LabVIEW based measure system for pulse wave transit time
CN210095711U (zh) 一种无创连续血压测量设备
Casacanditella et al. Indirect measurement of the carotid arterial pressure from vibrocardiographic signal: Calibration of the waveform and comparison with photoplethysmographic signal
CN114587307A (zh) 一种基于电容耦合电极的非接触血压检测仪及方法
Wang et al. An improved algorithm for noninvasive blood pressure measurement
RU2107456C1 (ru) Способ измерения давления в глубокой вене плеча и измерительная система для его реализации
RU2327414C1 (ru) Способ определения артериального давления по объемной компрессионной осциллограмме
Sidhu et al. Comparison of artificial intelligence based oscillometric blood pressure estimation techniques: a review paper

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