一种基于指端脉搏波的桡动脉脉搏波检测系统及方法【技术领域】，具体涉及一种基于指端脉搏波的桡动脉脉搏波检测系统及方法。[0002]桡动脉是人体若干浅表动脉中的一根，相比较于其他部位的脉搏波，桡动脉脉搏波的检测较方便，因此不论是中医的切脉，或是西医中心血管参数的无创伤检测，往往总是在桡动脉处获取信息。桡动脉处检测到的信息确实可以在一定程度上反映出人体的某些生理、病理特征，这一点早已被中医、藏医等传统医学的脉诊所证实。[0003]关于桡动脉脉搏波的检测，现在常用的方法是在桡动脉处佩戴一款脉搏传感器。可以使用的传感器种类有很多，如压电式脉搏传感器、压阻式脉搏传感器和加速度脉搏传感器等。但是和指端脉搏波的测量相比，桡动脉处的脉搏波测量较复杂，需要找准桡动脉的位置；而且与在手指上佩戴传感器相比，手腕处佩戴传感器对人体影响更大。虽然指端脉搏波采集比较方便，但是和桡动脉脉搏波相比，由于其受到手掌内微循环系统的影响，含有的特征信息量较少。
[0004]针对现有技术 存在的不足，本发明提供一种基于指端脉搏波的桡动脉脉搏波检测系统及方法。[0005]本发明的技术方案是:[0006]一种基于指端脉搏波的桡动脉脉搏波检测系统，包括指端脉搏传感器、桡动脉脉搏传感器、脉搏信号采集模块、中央处理模块和输入输出模块；
[0007]所述脉搏信号采集模块包括指端脉搏信号采集模块和桡动脉脉搏信号采集模块；
[0008]所述指端脉搏传感器、桡动脉脉搏传感器分别安置在指端处和手腕处，指端脉搏传感器的输出端、桡动脉脉搏传感器的输出端分别连接指端脉搏信号采集模块的输入端、桡动脉脉搏信号采集模块的输入端，指端脉搏信号采集模块的输出端、桡动脉脉搏信号采集模块的输出端分别连接中央处理模块的输入端；
[0009]输入输出模块包括按键/触摸屏、IXD屏和SD卡，按键/触摸屏连接中央处理模块的输入端，SD卡和IXD屏分别连接中央处理模块的输出端。
[0010]所述指端脉搏信号采集模块和桡动脉脉搏信号采集模块均包括前置放大电路、脉搏信号获取电路和压力信号获取电路；
[0011]所述脉搏信号获取电路包括高通滤波电路、第一低通滤波电路、第一二级放大电路和电平抬升电路；
[0012]所述压力信号获取电路包括第二低通滤波电路、第二二级放大电路和反相器电路；[0013]指端脉搏传感器的输出端、桡动脉脉搏传感器的输出端分别连接前置放大电路的输入端，前置放大电路的输出端分别连接高通滤波器的输入端和第二低通滤波器的输入端，高通滤波器的输出端连接第一低通滤波器的输入端，第一低通滤波器的输出端连接第一二级放大电路的输入端，第一二级放大电路的输出端连接电平抬升电路的输入端，电平抬升电路的输出端连接中央处理模块，第二低通滤波器的输出端连接第二二级放大电路的输入端，第二二级放大电路的输出端连接反相器电路的输入端，反相器电路的输出端连接中央处理模块。
[0014]采用所述的基于指端脉搏波的桡动脉脉搏波检测系统进行桡动脉脉搏波检测的方法，包括以下步骤:
[0015]步骤1:分别对指端脉搏传感器、桡动脉脉搏传感器进行压力校准，设定指端脉搏传感器检测脉搏波时对指端施加的最佳接触压力和桡动脉脉搏传感器检测脉搏波时对手腕桡动脉处施加的最佳接触压力；
[0016]步骤1.1:指端脉搏传感器、桡动脉脉搏传感器分别采集指端脉搏波和桡动脉脉搏波；
[0017]步骤1.2:中央处理模块控制指端脉搏传感器对指端施加不同的接触压力以及桡动脉脉搏波传感器对手腕桡动脉处施加不同的接触压力，采集不同接触压力状态下的指端脉搏波和烧动脉脉搏波信号；
[0018]步骤1.3:指端脉搏信号采集模块和桡动脉脉搏信号采集模块分别对采集到的指端脉搏波和桡动脉脉搏波进行放大、滤波和电平抬升处理，处理后的信号发送至中央处理模块；
[0019]步骤1.4:根据采集到的不同接触压力下的指端脉搏波信号的幅值和桡动脉脉搏波信号的幅值，确定指端脉搏波信号中幅值最大的脉搏波对应的接触压力和桡动脉脉搏波信号中幅值最大的脉搏波对应的接触压力，将该两个接触压力分别设定为指端脉搏传感器检测脉搏波时对指端施加的最佳接触压力和桡动脉脉搏传感器检测脉搏波时对手腕桡动脉处施加的最佳接触压力；
[0020]步骤2:在设定的指端脉搏传感器检测脉搏波时对指端施加的最佳接触压力下使用指端脉搏传感器采集受试者的指端脉搏波，在设定的桡动脉脉搏传感器检测脉搏波时对手腕桡动脉处施加的最佳接触压力下使用桡动脉脉搏传感器采集受试者的桡动脉脉搏波；
[0021]步骤3:指端脉搏信号采集模块和桡动脉脉搏信号采集模块分别对采集到的指端脉搏波和桡动脉脉搏波进行放大、滤波和电平抬升处理；
[0022]步骤4:建立基于传递函数的桡动脉脉搏波重建模型，该模型的输入为指端脉搏波，输出为烧动脉脉搏波；
[0023]步骤4.1:求解桡动脉脉搏波到指端脉搏波的频域传递函数；
[0024]步骤4.1.1:使用Welch’ s平均周期法估计的指端脉搏波与桡动脉脉搏波的互功率谱密度；
[0025]步骤4.1.2:使用Welch’ s平均周期法估计指端脉搏波的自功率谱密度；
[0026]步骤4.1.3:利用互功率谱密度除以自功率谱密度，即得到指端脉搏波与桡动脉脉搏波的频域传递函数；[0027]步骤4.2:确定基于传递函数的桡动脉脉搏波重建模型的窗宽和阶次；
[0028]步骤4.2.1:分别建立不同窗宽下的不同阶次的基于传递函数的桡动脉脉搏波重建模型；
[0029]步骤4.2.2:利用均方根百分误差和波形匹配度评估桡动脉脉搏波重建模型精度，得到桡动脉脉搏波重建模型的最佳窗宽和最佳阶次，即确定基于传递函数的桡动脉脉搏波重建模型的窗宽和阶次；
[0030]步骤4.3:根据桡动脉脉搏波重建模型的窗宽和阶次建立基于传递函数的桡动脉脉搏波重建模型；
[0031]步骤5:利用基于传递函数的桡动脉脉搏波重建模型和实时采集的指端脉搏波对桡动脉脉搏波进行重建，得到桡动脉脉搏波检测结果，该检测结果存储到SD卡并通过LCD
屏显示。
[0032]有益效果:
[0033]本发明的基于指端脉搏波的桡动脉脉搏波检测系统及方法所得到的桡动脉脉搏波波形与实际测量的桡动脉脉搏波波形的最佳波形匹配度均在70%至92%范围内，该装置适用于桡动脉脉搏波数据的采集，与现有技术中的桡动脉脉搏检测装置相比，使桡动脉脉搏波的检测更加方便。



[0034]图1为本发明
[0035]图2为本发明【具体实施方式】的指端桡动脉脉搏信号采集模块信号走向示意图；
[0036]图3为本发明【具体实施方式】的前置放大电路原理图；
[0037]图4为本发明【具体实施方式】的高通滤波电路原理图；
[0038]图5为本发明【具体实施方式】的第一低通滤波电路原理图；
[0039]图6为本发明【具体实施方式】的第一二级放大电路原理图；
[0040]图7为本发明【具体实施方式】的电平抬升电路原理图；
[0041]图8为本发明【具体实施方式】的第二低通滤波电路原理图；
[0042]图9为本发明【具体实施方式】的第二二级放大电路原理图；
[0043]图10为本发明【具体实施方式】的反相器电路原理图；
[0044]图11为本发明【具体实施方式】的桡动脉脉搏波检测的方法流程图；
[0045]图12为本发明【具体实施方式】的SD卡电路原理图；
[0046]图13为本发明【具体实施方式】的传递函数建模训练和测试过程框图；
[0047]图14为本发明【具体实施方式】的不同窗宽下的最佳模型估计结果的PRD和FIT分布图；
[0048]图15为本发明【具体实施方式】的一位被试者窗宽为1024时各阶次模型估计结果的PRD和FIT分布图；
[0049]图16为本发明【具体实施方式】的系统工作流程图。
【具体实施方式】[0050]下面结合附图对本发明的【具体实施方式】做详细说明。
[0051]如图1所示，基于指端脉搏波的桡动脉脉搏波检测系统，包括指端脉搏传感器、桡动脉脉搏传感器、脉搏信号采集模块、中央处理模块和输入输出模块；
[0052]脉搏信号采集模块包括指端脉搏信号采集模块和桡动脉脉搏信号采集模块；
[0053]指端脉搏传感器、桡动脉脉搏传感器分别安置在指端处和手腕处，指端脉搏传感器的输出端、桡动脉脉搏传感器的输出端分别连接指端脉搏信号采集模块的输入端、桡动脉脉搏信号采集模块的输入端，指端脉搏信号采集模块的输出端、桡动脉脉搏信号采集模块的输出端分别连接中央处理模块的输入端；
[0054]输入输出模块包括按键/触摸屏、IXD屏和SD卡，按键/触摸屏连接中央处理模块的输入端，SD卡和IXD屏分别连接中央处理模块的输出端。
[0055]指端脉搏信号采集模块和桡动脉脉搏信号采集模块均包括前置放大电路、脉搏信号获取电路和压力信号获取电路；
[0056]脉搏信号获取电路包括高通滤波电路、第一低通滤波电路、第一二级放大电路和电平抬升电路；
[0057]压力信号获取电路包括第二低通滤波电路、第二二级放大电路和反相器电路；
[0058]指端脉搏传感器的输出端、桡动脉脉搏传感器的输出端通过接插件(如航空插头、DB插头等)分别连 接前置放大电路的输入端，如图2所示，前置放大电路的输出端分别连接高通滤波器的输入端和第二低通滤波器的输入端，高通滤波器的输出端连接第一低通滤波器的输入端，第一低通滤波器的输出端连接第一二级放大电路的输入端，第一二级放大电路的输出端连接电平抬升电路的输入端，电平抬升电路的输出端连接中央处理模块，第二低通滤波器的输出端连接第二二级放大电路的输入端，第二二级放大电路的输出端连接反相器电路的输入端，反相器电路的输出端连接中央处理模块，其中，电平抬升电路输出脉搏波信号，反相器电路输出压力信号。
[0059]前置放大电路主芯采用仪表运算放大器AD620，其电路如图3所示。本实施方式为了分别得到脉搏信号(交流信号)和压力信号(直流信号)，引入脉搏信号获取电路和压力信号获取电路，该两个电路以TI公司的0PA2277为主芯片。
[0060]脉搏信号获取电路包括高通滤波电路、第一低通滤波电路、第一二级放大电路和电平抬升电路。信号走向为:前置放大输出信号首先接高通滤波电路(如图4所示)滤除信号中的直流成分；然后连接第二低通滤波电路(如图5所示)滤除脉搏信号中的高频噪声；然后接入第一二级放大电路(如图6所示)将脉搏信号放大；最后接入电平抬升电路(如图7所示)后得到满足AD采样条件的脉搏波信号。
[0061]脉搏传感器压力信号获取电路包括:第二低通滤波电路、第二二级放大电路和反相器电路。信号走向为:前置放大输出信号首先接第二低通滤波电路(如图8所示)滤除信号中的交流成分；然后接入第二二级放大电路(如图9所示)将压力信号放大，最后接入反相器电路(如图10所示)得到满足AD采样条件的压力信号。
[0062]中央处理模块采用ARM处理器，指端脉搏信号采集模块和桡动脉脉搏信号采集模块的输出端接ARM处理器的AD采样模块模拟信号输入端，由ARM中内置AD模块完成AD采样；因为脉搏信号频率范围集中在0-40Hz范围内，根据采样定理，本实施方式中将AD采样频率设定为IOOHz，采样精度为12位。[0063]SD卡电路如图12所示，SD卡与ARM处理器间通过SPI接口实现通信，采用弹出式SD卡插槽，电路中“C/D”脚用于检测插槽内是否插入SD卡。在进行SD卡初始化时，SPI输出时钟频率设置为330Khz，SD卡初始化成功后，对SD卡数据进行读写操作时，将SPI输出时钟频率设置为16.5MHz ο
[0064]在进行数据存储时，首先调用SD卡内的FAT16文件系统完成对数据文件的创建，本实施方式中数据采集文件和结果文件的格式均为“dat”格式，命名方式分别为“fpwavOOO”和“rpwavOOO”，文件名后三位数字分别代表被试者检查编号的个位、十位和百位，便于数据的记录。
[0065]本实施方式中使用按键或触摸屏实现对系统各个功能模块间的切换。系统中包含5个(Keyl、Key2、Key3、Key4、Reset)按键,通过外中断的方式,在各个按键中断服务子程序中实现相应的功能；
[0066]本实施方式中央处理模块采用三星公司型号为S3C6410的ARM板，其内部含有触摸屏接口和LCD屏接口，使用传统的四线触屏方式，通过对点击触摸屏动作产生的电压值的采样结果的判断，得到坐标信息，在触摸屏中断服务子程序，根据该坐标信息，执行相应的操作。
[0067]本实施方式的IXD屏选用群创公司型号为AT070TN92(7寸屏)的IXD液晶屏，采用18位(6: 6: 6) RGB真彩色模式。
[0068]本实施方式中可以通过IXD屏显示数据采集和数据分析请假，实时显示当前被试者的编号和实测的脉搏数据波形、用于数据分析的实测数据和检测结果数据和重建的桡动脉脉搏波与实测的桡 动脉脉搏波间的波形匹配度(FIT)值和均方根百分误差(PRD)值。
[0069]LCD屏显示的字符为16*16的点阵，字符取模方式为左高位横向取模，数据排列方式为从左到右，从上到下；数据采集界面的波形显示区域背景颜色设置为黑色，网格线条为红色虚线(点与点之间隔三个像素点)，横坐标值代表时间，打点间隔为10ms，每100个像素点长度代表ls，纵坐标代表采样幅值，单位为“伏特”，两通路脉搏波(桡动脉和指端脉搏)采样结果值分别映射为区域I和区域2的纵坐标值，映射关系如公式(I)和(2)所示:
[0070]Y1 = 51+(4095^)/22(I)
[0071]Y2 = 259+(4095-V3)/22(2)
[0072]其中，Y1和Y2分别代表显示区域I和显示区域2的纵坐标值，V1和V2分别代表两通路数据采样结果值。
[0073]采用基于指端脉搏波的桡动脉脉搏波检测系统进行桡动脉脉搏波检测的方法，其流程如图11所示，包括以下步骤:
[0074]步骤1:分别对指端脉搏传感器、桡动脉脉搏传感器进行压力校准，设定指端脉搏传感器检测脉搏波时对指端施加的最佳接触压力和桡动脉脉搏传感器检测脉搏波时对手腕桡动脉处施加的最佳接触压力；
[0075]步骤1.1:指端脉搏传感器、桡动脉脉搏传感器分别采集指端脉搏波和桡动脉脉搏波；
[0076]步骤1.2:中央处理模块控制指端脉搏传感器对指端施加不同的接触压力以及桡动脉脉搏波传感器对手腕桡动脉处施加不同的接触压力，采集不同接触压力状态下的指端脉搏波和烧动脉脉搏波信号；[0077]步骤1.3:指端脉搏信号采集模块和桡动脉脉搏信号采集模块分别对采集到的指端脉搏波和桡动脉脉搏波进行放大、滤波和电平抬升处理，处理后的信号发送至中央处理模块；
[0078]步骤1.4:根据采集到的不同接触压力(50g、60g、...140g、150g)下的指端脉搏波信号的幅值和桡动脉脉搏波信号的幅值，确定指端脉搏波信号中幅值最大的脉搏波对应的接触压力和桡动脉脉搏波信号中幅值最大的脉搏波对应的接触压力,将该两个接触压力分别设定为指端脉搏传感器检测脉搏波时对指端施加的最佳接触压力和桡动脉脉搏传感器检测脉搏波时对手腕桡动脉处施加的最佳接触压力；
[0079]本实施方式的指端脉搏传感器和桡动脉脉搏传感器可以选用电机加压式/充气加压式传感器。首先，系统对两个传感器进行调零操作，调零操作成功后，分别完成对指端脉搏传感器和桡动脉脉搏传感器压力的校准，使其工作在最佳压力状态下。
[0080]步骤2:在设定的指端脉搏传感器检测脉搏波时对指端施加的最佳接触压力下使用指端脉搏传感器采集受试者的指端脉搏波，在设定的桡动脉脉搏传感器检测脉搏波时对手腕桡动脉处施加的最佳接触压力下使用桡动脉脉搏传感器采集受试者的桡动脉脉搏波；
[0081]步骤3:指端脉搏信号采集模块和桡动脉脉搏信号采集模块分别对采集到的指端脉搏波和桡动脉脉搏波进行放大、滤波和电平抬升处理；
[0082]步骤4:建立基 于传递函数的桡动脉脉搏波重建模型，该模型的输入为指端脉搏波，输出为烧动脉脉搏波；
[0083]本实施方式中求解指端脉搏波到桡动脉脉搏波频域传递函数Txy (f)方法的原理为:使用输入输出信号的互功率谱密度除以输入信号的功率谱密度，其表达式为:
[0084]Txy (f) = (Pyx (f))/(Pxx (f))(3)
[0085]其中，Pyx(f)为X与y间的互功率谱密度，Pxx(f)为x的功率谱密度，二者均由Welch’s平均周期法估计。将桡动脉脉搏波Pr作为输入，指端脉搏波Pf作为输出，可求得桡动脉脉搏波到指端脉搏波的频域传递函数Trf。根据求得的Trf的幅值和相位信息即可求得不同阶次的时域FIR模型。
[0086]步骤4.1:求解桡动脉脉搏波到指端脉搏波的频域传递函数；
[0087]步骤4.1.1:使用Welch’ s平均周期法估计的指端脉搏波与桡动脉脉搏波的互功率谱密度；
[0088]步骤4.1.2:使用Welch’ s平均周期法估计指端脉搏波的自功率谱密度；
[0089]步骤4.1.3:利用互功率谱密度除以自功率谱密度，即得到指端脉搏波与桡动脉脉搏波的频域传递函数；
[0090]从上述步骤中可知求解桡动脉脉搏波到指端脉搏波的频域传递函数的核心在于Welch’s平均周期法的使用。Welch’s平均周期法是为了克服周期法的随机起伏问题提出的。随机起伏是指在使用周期法计算信号的功率谱时，随着所取的信号序列长度不同，得到的周期图不同。平均周期法是对信号进行分段加窗，首先得到每段信号的周期图，然后对所得的周期图取平均，得到信号的功率谱。即对信号加窗处理可以减小随机起伏。因此，窗宽的选择对信号功率谱的求解影响很大，进而对输入输出信号的频域传递函数的求解和FIR模型的精度起到决定性作用。[0091]步骤4.2:确定基于传递函数的桡动脉脉搏波重建模型的窗宽和阶次；
[0092]步骤4.2.1:分别建立不同窗宽下的不同阶次的基于传递函数的桡动脉脉搏波重建模型；
[0093]步骤4.2.2:利用均方根百分误差和波形匹配度评估桡动脉脉搏波重建模型精度，得到桡动脉脉搏波重建模型的最佳窗宽和最佳阶次，即确定基于传递函数的桡动脉脉搏波重建模型的窗宽和阶次；
[0094]本实施方式为了客观地估计模型算法的性能，采用在心电信号压缩领域常用的评价指标:均方根百分误差(PRD)和波形匹配度(FIT)。二者的计算公式为: