一种基于四肢血压脉搏的心脏血管功能参数的计算方法
【技术领域】
[0001] 本发明属于生物力学技术，涉及计算机程序，电路设计等，可用于心脏血管功能参 数的无创计算。
【背景技术】
[0002] 心血管疾病已经成为危害全球和我国成人健康并导致死亡的首要病症，但是心血 管疾病早期临床症状不明显，直至并发症恶化才会对心血管疾病进行诊断。目前，心血管疾 病的诊断技术主要分为直接法和间接法两种。直接测量法可提供精度较高的、可记录的连 续数据，但是这些参数数据的获取是通过有创检测手段得到，该方法不仅存在着一定的风 险性，而且对医生的临床水平或医疗设备的技术都有当高的要求。间接测量法在临床上应 用较多的是通过心血管仪器诊断，如超声心动图法、心阻抗法、螺旋CT、核医学、磁共振 (MRI)等方法;虽然是无创检测，但是检测环境要求严格、不能随身携带等缺点，被检查者只 能在专门的医疗场所、由专业的医疗人员进行检测。虽然目前的心血管疾病的检测技术在 临床上已经结合应用并提供了好的预测方法，但是仍然有很多高危患者未被发现。
[0003] 国内外学者一直致力于基于脉搏波研究心脏血管系统的无创检测技术，已经建立 设计了无创检测模型、方法和仪器装置，但仍存在一些问题，如:评价参数单一、模型中涉及 较多的假设和猜测及产品技术壁皇等问题。所以心脏血管系统的无创检测方法有待于进一 步改进与优化。
[0004] 此外，由于现代人们的健康意识大大提高，针对上述心功能的检测方法成本高、操 作复杂且只适用于专门的医疗机构，离心血管疾病的健康普检及家庭保健还存在很大差 距。本发明提出了基于四肢血压脉搏的心脏血管的无创计算方法，该方法将为心脏血管功 能检测提供一种经常性、无创伤、易操作、快输出、价格低廉的检测方法和手段;同时也为心 脏血管系统的无创检测方法提供理论基础。
【发明内容】
[0005] 本发明的目的是提供一种基于四肢血压脉搏波的心脏血管系统功能参数的计算 方法，该方法可以无创同步检测人体的四肢肱动脉和踝动脉的血压和脉搏波的8路信息，依 据8路信息和非线性脉搏波理论同步检测人体的心脏功能、动脉网硬化程度、血液流变学指 标等心脏血管功能的状态参数。其目的在于，提高无创检测计算心功能参数的准确性，避免 现有的有创检测手段的有创性及高成本等问题。
[0006] 所述的四肢肱动脉和踝动脉搏波需进行波形平均，其平均的是四肢压力传感器提 取的压力脉搏波分量，波形平均的步骤，包括:首先读取波形数据，标定完整的波形求出周 期，然后对波形进行周期和幅值归一化处理，使脉搏波单波数据具有统一的幅值和长度，并 求得平均波形。
[0007] 所述的非线性脉搏波理论，建立四肢血压、血流量两个基本血流动力学参量和四 肢脉搏波之间的关系式，并根据人体心脏血管系统的固有特征参数或者几何参数，推导出 四肢血压脉搏数据、血流量与心脏血管功能参数的关系式并进行计算。
[0008] 所述的计算方法，主要利用非线性脉搏波理论和四肢血压脉搏波数据，建立压力 梯度与血流量，血流量与心脏血管功能参数的关系模型；利用数值分析计算出一个心动周 期内不同时刻的四肢血流量的数值解，利用统计学的回归方法分别推算四肢血流量的回归 系数，修正得出人体血流量及平均血流量;在血流量和四肢脉搏血压数据信息基础上，计算 求解心脏血管功能参数参量。算法中将通过误差分析对计算结果进行判断是否达到收敛， 若收敛，则输出结果;否则将进行迭代计算。
[0009] 所述的8路信息的检测装置，包括脉搏波传感器、血压传感器，A/D转换电路、接口 电路及附属电路等，对同步检测并采集的四肢血压及脉搏波8路数据信息，通过USB集线器 传输到PC端，通过上位机的程序控制下位机，并对检测到的血压波形数据进行处理、存储和 输出。该硬件部分，操作简单，易检测，成本低，用时短。
[0010] 所述的8路信息的检测状态，包括静息状态、定量运动和运动后再次静息状态下的 血压和脉搏波信号。其目的在于，用静态参数表示个体某一时刻的生理状态，用动态参数表 示个体经历运动时的自身调节功能及运动结束后个体做出的自我调节反应。
[0011] -种基于四肢血压脉搏的心脏血管功能参数的计算方法，可以为人体的心脏血管 系统提供一种经常性、无创伤、易操作、快输出、价格低廉的检测手段;可以适用于社区医 疗、康复保健、运动健身、休闲疗养等医疗保健机构的健康普检。
[0012] 本发明的特征在于，在上位机和下位机中依次按以下步骤实现的：
[0013] 步骤⑴，上位机和下位机初始化：
[0014] 上位机输入:被测者的个人信息，包括性别、年龄、身高、体重、肩宽既往病史等信 息。
[0015] 下位机检测:被测者肱动脉和踝动脉四肢的脉搏波和血压。
[0016] 上位机记录:被测者肱动脉和踝动脉四肢的脉搏波波形和血压压力值。
[0017] 步骤(2)，上位机自动检测、记录、保存四肢的脉搏波波形图，对得到的波形图利用 周期和幅值归一化处理进行波形平均，然后按下式计算分别得到一个心动周期内被测者的 四肢平均动脉压P mf，平均动脉压为一个心动周期内四肢分别对应的动脉压力的平均值。
[0019] T为一个心动周期，f为四肢的标号，€=1，2,3,4分别表示右上肢、左上肢、右下肢、 左下肢，P(t)f为一个心动周期内不同时刻各采样点对应的脉压值，（f = 1，2，3，4)。
[0020] 步骤(3)，按下式计算波形因子kf，（f =1，2,3,4)。
[0022] Pdf为舒张压，Psf为收缩压，均为四肢肱动脉和踝动脉的实际测量值，（f=l，2,3， 4)〇
[0023]步骤(4)，按下式计算肱动脉和踝动脉压力脉搏波的传播速度c(t)f:
[0025] P(t)f为一个心动周期内不同时刻各采样点对应的脉压值，（f = 1，2，3，4);
[0026] p为血液密度，p = l .05*10-3kg/cm3;
[0027] bf为脉搏波波形参数，按下式计算，（f=l，2,3,4):
[0029] a为与速度剖面有关的参数，a = 0.57;
[0030] 步骤(5)，按下式计算肱动脉和踝动脉压力脉搏波的空间压力梯度F(t)f，从而得 到一个心动周期内不同时刻的压力梯度，（f = l，2,3,4)。
[0032]
为一个心动周期内不同时时刻压力对时间的导数，clf(f=l，2,3,4)为四肢 动脉内的血液脉搏波的传播速度，ω 四肢动脉内脉搏波的角频率：
四肢动脉内脉搏波的系数m = 1，2，3……，按下式递推得出：
[0039] Pm*平均动脉压，P为血压在动脉内坐标Z和时间t的函数，其傅立叶形式为：
幅值和频率等由实测波形获得。
[0040] 步骤(6)，按下式计算一个心动周期内不同时刻肱动脉和踝动脉四肢的血流量Q (t)f，（f=l，2,3,4)〇
[0042] X(t)f为血流量随时间变化的一次幂系数；e(t)f为血流量随时间变化的二次幂系 数;A(t) f为一个心动周期内血管半径随时间变化的函数，主要描述一个心动周期内血管弹 性收缩扩张的动态过程，Q(t) f为一个心动周期内不同时刻肱动脉和踝动脉各采样点对应 的血流量值(f = l，2,3,4)。
[0043] A(t)f用下式表示，（f=l，2,3,4):
[0045] α为生理条件下人体修正系数，是计算A(t)f时涉及一个心动周期内血流量随时间 变化时的非线性脉搏波传播系数，α = 0.57;
[0046] β为生理条件下人体修正系数，是计算A(t)f时一个心动周期内血管半径随时间 变化时的非线性脉搏波传播系数，与血管半径弹性扩张有关，β=〇. 1.173;
[0047] γ为血液