一种肱动脉内皮功能的电阻抗检测分析装置及方法
【技术领域】
[0001] 本发明属于医疗仪器技术领域，具体涉及一种血管功能的无创伤评估方法，特别 涉及一种肱动脉内皮功能的无创伤检测分析装置及方法。
【背景技术】
[0002] 近年来，因心脑血管疾病发病和死亡人数急剧增多，治疗费用也急剧增加，给患者 家庭和社会带来了沉重的负担。大量研究表明，动脉硬化程度增加是心脑血管疾病的独立 危险因子之一。动脉硬化作为一种全身性的病理变化过程，会直接导致重要脏器的损坏，威 胁人类健康。
[0003] 动脉内皮功能受损或者功能障碍是动脉硬化早期的主要表现。动脉内皮作为血管 流动和血管壁之间的功能性屏障，不仅维持血管结构的完整，还在调节血管细胞生长、调节 抗凝及纤溶系统、介导炎症与免疫、调节白细胞与血小板在血管内皮粘附、调节脂质氧化、 调节血管通透性等方面扮演十分重要的角色。如果内皮结构破坏、功能紊乱，血液中的脂肪 便会进入血管壁，与中膜相互作用，最终形成动脉粥样硬化性斑块。
[0004] 肱动脉不仅是测量血压的部位，也是检测全身动脉功能的一个重要"窗口"。研究 表明，肱动脉内皮功能与心脏冠状动脉的功能呈正相关，因此，正确评价肱动脉内皮功能， 有助于及时发现动脉功能异常，从而为心脑血管疾病的早期诊断和治疗提供一定的参考。 胧动脉血流介导性扩张(flow-mediated dilation，FMD)反映了动脉内皮在血流突然增加 时的舒张能力。临床上，一般使用高空间分辨率的超声多普勒技术无创伤检测肱动脉血流 介导性扩张引起的动脉管径扩张率来评价动脉的内皮功能。这种动脉管径扩张率定义为：
，Dbase3line3表示静息状态下动脉搏动管径的平均值，D peak 表示动脉扩张后搏动管径平均值的最大峰值。
[0005] 虽然利用高空间分辨率的超声多普勒技术测量肱动脉FMD评估内皮功能被美国心 脏学会认为是一项"有前途的技术"，但用该方法检测动脉管径的变化需要对动脉管径进行 精确定位，不仅操作程序烦琐、耗时，而且相关的仪器设备笨重、价格昂贵，不便于推广到社 区或者家庭简单使用。
【发明内容】
[0006] 用气袖无创伤阻断肱动脉血流一定时间（一般3-5分钟）后释放，将引起下游前臂 动脉的血流动力学环境特别是作用于动脉内皮的壁面剪应力发生改变，改变的剪应力将诱 发前臂动脉内皮释放血管舒张因子N0,进一步引起前臂动脉的管径发生改变，这一过程称 为肱动脉血流介导性扩张(FMD)。在该过程中动脉血流量发生改变，动脉血液的电导特性发 生相应变化，引起前臂动脉的总电阻抗发生改变。通过检测静息和FMD实施过程电阻抗的变 化情况可以反映 FMD引起的动脉管径变化，从而达到间接评估动脉内皮功能的目的。
[0007] 前臂可看成一个电导体，其电阻抗由前臂内所有动脉、静脉、皮肤、肌肉、骨骼和其 他组织的电导特性共同决定。由于FMD过程中只有动脉血流的电导特性才发生明显变化，因 此所测量的前臂电阻抗的改变就是前臂内动脉段电阻抗的改变。将前臂内的动脉等效于一 个纵向长度L不变的圆筒形腔室，其半径R则在动脉正常搏动和FMD过程中发生改变，该动脉 段的电阻抗Z a表达为
[0009] 式中，σ为血液电导率。
[0010] 假定除动脉外的组织和器官总阻抗为常量ΖΡ，则前臂总电阻抗Ζ为心和23的并联， 表示为
[0012]进一步假定外围组织和器官的阻抗2(5在?1?)实施前后保持不变，则对公式(2)两边 微分可得到
[0014] 由于ZP?Za，可近似为
[0015] Δ Z ? Δ Za (4)
[0016] 即前臂电阻抗的改变就是前臂内动脉电阻抗的改变。
[0017] 本发明的前臂电阻抗检测装置的技术方案如下：
[0018] 本发明的装置包括前臂上双电极、前臂下双电极、高频低幅恒流源发生器、交流电 压检测电路、信号放大器、滤波电路、A/D转换电路和计算机。前臂上双电极置于待测手臂的 上部，前臂下双电极置于待测手臂的下部；高频低幅恒流源发生器和一对前臂上下电极构 成电流回路；另一对前臂上下电极和交流电压检测电路构成回路检测电压变化。如图1所 示，高频低幅恒流源发生器给前臂上电极和前臂下电极加一个一定频率ω的高频低幅的恒 流电源，使前臂通过一个强度恒定的交流电流I，这样在前臂长度为L段两端就建立起一定 的交流电压U。交流电压检测电路检测到的交流电压U经放大器、滤波，A/D转换将模拟的电 压信号转换成数字信号输入计算机，并根据如下公式(5)计算电阻抗Ζ随时间的变化曲线：
[0020]由公式(5)计算的电阻抗信号存入计算机进一步进行分析。
[0021 ]本发明的前臂电阻抗信号分析方法及提供的特征参数如下：
[0022]分别检测FMD实施前静息状态和实施后的前臂电阻抗信号Zo(t)和ZF(t)。对前臂电 阻抗信号Zo(t)和ZF(t)进行小波分析获得特征参数：
[0023]对于电阻抗信号2(〇46[以32]，该信号的连续小波变换被定义为：
[0025]其中，平5/是基函数Ws,t的共辄，基函数Ws,t是母小波Φ⑴经过时间平移t和尺度 伸缩s变换得到：
[0027]通过小波变换后，信号在时间尺度上的能量密度被定义为：
[0029] 选择用Morlet母小波进行时频分析，其中，Morlet小波的时域表达式如下：
[0031] 其中，ω Q = 23T，C的值取决于母小波的形状，
是Ws,t(u)的傅里叶变换。
[0032] 对连续的电阻抗波形Z(t)进行小波分析，可得到具有五个特征峰的频谱图，频谱 区间为0.0005-2.0Hz。其中，频率区间0.005-0.021Hz对应动脉内皮的生理活动，计算频谱 图中该频率区间的特征频率f。、绝对能量E。、平均振幅A。：
[0033] (1)特征频率乜
[0035] 在频率区间 0·005-0·021Ηz 内取 A(f)的极大值 A(fpl)，A(fp2)，A(fp3)···，若 fc 满足
[0036] A (fc)=max( A (fpi), A (fp2), A (fp2), . . .) (11)
[0037] 则fc为频率区间[fcl，fc2]内的特征频率。
[0042]其中，E。(f。1，f。2)表示信号在频率间隔[f。1，f。 2]和时间间隔[t1，t2]内的能量。A。 (fcl，fc2 )表不在频率间隔[fcl，fc2