一种基于脉搏波的血压检测装置及方法与流程

本发明属于生物信号检测技术领域，更具体地，涉及一种基于脉搏波的血压检测装置及方法。

背景技术：

血压是指心脏在进行血液循环时对血管壁产生的压力，是反应人体心血管功能的核心参数。如今高血压患者的人数在逐年增多，常规的血压检测必不可少，现有的家用电子血压计测量地点相对固定，不能实时监控血压变化，且需要借助于袖带，在便携性、实时性和舒适度均有所受限，因此研究可穿戴血压检测设备很有意义。基于脉搏波法的血压检测突破了袖带的限制，非常适合可穿戴设备，但是易受干扰，准确度有限，至今还处在探索完善的阶段，尚未进入商业领域。

因此，针对可穿戴血压检测设备的应用需求，研究一种运算量小、准确度高的血压检测装置具有很大的市场价值。

技术实现要素：

针对现有技术的缺陷，本发明的目的在于提供一种基于脉搏波的血压检测装置，旨在解决现有技术提供的血压检测装置占用面积大和准确度不高的问题。

本发明提供了一种基于脉搏波的血压检测装置，血压检测装置包括：ppg信号抗干扰处理模块、ptt提取模块和血压值计算模块；ppg信号抗干扰处理模块用于去除原始ppg信号中的高频噪声、基线漂移和运动伪影，并输出相对干净的ppg信号；ptt提取模块用于利用干净的ppg信号提取特征参数ptt；血压值计算模块用于对特征参数ptt进行曲线拟合并获得sbp值与dbp值。

更进一步地，ppg信号抗干扰处理模块包括：第一低通滤波模块、加权均值滤波模块、自相关运算模块、周期值计算模块和周期平滑滤波模块；第一低通滤波模块的输入端用于接收原始ppg信号ppg0，根据输入端接收的信号输出第一ppg信号ppg1，其中第一ppg信号ppg1是去除了高频噪声后的原始ppg信号ppg0；加权均值滤波模块的输入端用于接收第一ppg信号ppg1，根据输入端接收的信号输出第二ppg信号ppg2；其中第二ppg信号ppg2是去除了基线漂移后的第一ppg信号ppg1；自相关运算模块的输入端用于接收第二ppg信号ppg2，根据输入端接收的信号输出第三ppg信号ppg3；其中第三ppg信号ppg3是自相关运算后的第二ppg信号ppg2；周期值计算模块的输入端用于接收第三ppg信号ppg3，根据输入端接收的信号输出ppg信号的周期值p；周期平滑滤波模块的第一输入端用于接收第二ppg信号ppg2，第二输入端用于接收ppg信号的周期值p，根据输入端接收的信号输出第四ppg信号ppg4，其中第四ppg信号ppg4是去除了运动伪影后的第二ppg信号ppg2。

更进一步地，ptt提取模块包括：二阶微分运算模块、第二低通滤波模块和特征参数提取模块；所述二阶微分运算模块的输入端用于接收去除运动伪影后的信号ppg4，根据输入端接收的信号输出第一apg信号apg1，其中第一apg信号apg1是二阶微分运算后的第四ppg信号ppg4；所述第二低通滤波模块的输入端用于接收第一apg信号apg1，根据输入端接收的信号输出第二apg信号apg2，其中第二apg信号apg2是去除了残余高频噪声后的第一apg信号apg1；所述特征参数提取模块的输入端用于接收第二apg信号apg2，根据输入端接收的信号输出特征参数ptt。

更进一步地，血压值计算模块包括：sbp计算模块和dbp计算模块；所述sbp模块的输入端用于接收特征参数ptt，根据输入端接收的信号输出收缩压sbp；dbp模块的输入端用于接收特征参数ptt，根据输入端接收的信号输出舒张压dbp。

本发明还提供了一种基于脉搏波的血压检测方法，包括下述步骤：

s1：去除原始ppg信号中的高频噪声、基线漂移和运动伪影后，获得相对干净的ppg信号；

s2：利用干净的ppg信号提取特征参数ptt；

s3：将特征参数ptt进行曲线拟合，并获得sbp值与dbp值。

更进一步地，在步骤s1中，将原始ppg信号依次经过低通滤波处理、加权均值滤波处理、自相关运算、周期值计算以及周期平滑滤波处理后获得所述相对干净的ppg信号。

更进一步地，在步骤s1中，所述加权均值滤波处理方法具体为：

(1)获得第一ppg信号ppg1中的第n个数据ppg1(n)，其中n的范围为1～1024；

(2)判断所述信号ppg1(n)是否为ppg1信号中的第1个数据，若是，则进入步骤(3)；若否，则进入步骤(4)；

(3)计算m(n)＝ppg1(1)，其中ppg1(1)为ppg1信号中的第1个数据，m(n)为ppg1(n)的均值信号，并进入步骤(5)；

(4)计算m(n)＝{m(n-1)+[ppg1(n+63)-m(n-1)]}/64，其中m(n-1)为前一个均值信号，ppg1(n+63)为当前数据之后的第63个数据，并进入步骤(5)；

(5)获得ppg2(n)＝ppg1(n)-m(n)，其中ppg2(n)表示第二ppg信号ppg2中的第n个数据；

(6)判断n是否小于1024，若是，则进入步骤(7)；若否，则结束运算；

(7)n自加1，重复步骤(1)～(6)。

加权均值滤波模块运算量小，通过加权均值滤波能有效去除ppg信号中的基线漂移。

更进一步地，在步骤s1中，所述周期值计算方法具体包括：

(1)获得第三ppg信号ppg3中的第m个数据ppg3(m)，其中m的范围为1～1024；

(2)将当前的自相关运算后的信号ppg3(m)与前后两个自相关运算后的信号ppg3(m-1)和ppg3(m+1)分别做对比，判断所述信号ppg3(m)是否为最大值，若是，则进入步骤(3)；若否，则进入步骤(1)；

(3)判定该点为波峰，并记录该点坐标；

(4)判断相邻波峰的间距是否大于等于阈值d，若是，则进入步骤(5)；若否，则进入步骤(1)；

(5)判定该波峰为有效波峰，并记录有效波峰的数量p与坐标xp；

(6)计算np＝p-1，其中np为第三ppg信号ppg3的周期数量；

(7)计算sp＝xp(max)-xp(min)，其中sp为第三ppg信号ppg3的周期总和，xp(max)为最后一个有效波峰的坐标，xp(min)为第一个有效波峰的坐标；

(8)获得p＝sp/np，其中p为第三ppg信号ppg3的平均周期值，也是第二ppg信号ppg2的平均周期值。

通过提取第三ppg信号ppg3的有效波峰可以确定ppg3的平均周期值，也就是第二ppg信号ppg2的平均周期值。

更进一步地，在步骤s2中，所述提取特征参数ptt的方法具体为：

(1)获得第二apg信号apg2中的第n个数据apg2(n)，其中n的范围为1～1024；

(2)将当前第二apg信号apg2(n)与前后两个第二apg信号apg2(n-1)和apg2(n+1)分别做对比，判断所述信号apg2(n)是否为最大值、是否为最小值，若是最大值，则进入步骤(3)；若是最小值，则进入步骤(4)；若既不是最大值也不是最小值，则进入步骤(1)；

(3)判定该点为波峰，记录该点坐标；

(4)判定该点为波谷，记录该点坐标；

(5)判断相邻波谷的间距是否大于等于阈值l且波谷幅值小于等于阈值m，若满足了间距阈值、幅值阈值的条件，则进入步骤(6)；若不满足间距阈值、幅值阈值的条件，则进入步骤(1)；

(6)判定为有效波谷，记录其坐标与数量t；

(7)根据步骤(5)、(6)判断波峰是否为有效波谷两侧最近的波峰，若是，则进入步骤(8)；若否，则进入步骤(1)

(8)判定波峰是有效波谷两侧最近的波峰，记录两波峰的坐标差δxp；

(9)计算sptt＝sptt+δxp，其中sptt为ptt的总和；

(10)计算nptt＝t，其中nptt为ptt的数量；

(11)获得ptt＝sptt/nptt，其中ptt为平均ptt值。

通过多重定位的方式对apg信号的波峰和波谷进行筛选，所提取的ptt值更加准确。

通过本发明所构思的以上技术方案，对采集的ppg信号具有较好的抗干扰处理效果，能够准确提取ptt，最终通过曲线拟合测得的sbp误差为1.14±2.15mmhg，dbp误差为0.86±1.38mmhg，均在aami颁布的标准5±8mmhg(平均绝对误差±标准差)范围内。本发明对ppg信号中的高频噪声、基线漂移和运动伪影都进行了滤除，并且在ptt提取时增加了一次低通滤波，提高了检测的准确率。本发明在去除基线漂移时，利用加权均值来替代传统均值滤波，省去了求和步骤，运算量更小；在去除运动伪影时，利用自相关运算来替代传统周期平滑滤波里的内插重构，硬件实现更加方便；这两种方法都更适用于可穿戴血压检测设备。

附图说明

图1是本发明实施例提供的总体结构示意图；

图2是本发明实施例提供的低通滤波模块的算法流程图；

图3是本发明实施例提供的加权均值滤波模块的算法流程图；

图4是本发明实施例提供的自相关运算模块的算法流程图；

图5是本发明实施例提供的周期值计算模块的算法流程图；

图6是本发明实施例提供的周期平滑滤波模块的算法流程图；

图7是本发明实施例提供的二阶微分运算模块的算法流程图；

图8是本发明实施例提供的ptt提取的示意图；

图9是本发明实施例提供的特征参数提取模块的算法流程图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

本发明的目的在于提供一种基于脉搏波的血压检测装置，具有无需袖带，运算量小，准确度较高等特点。它的主要作用是对采集的光电容积脉搏波(photoplethysmograph，ppg)信号进行抗干扰处理，提取特征参数脉搏波传导时间(pulsetransittime，ptt)，最终通过曲线拟合测得收缩压(systolicbloodpressure，sbp)与舒张压(diastolicbloodpressure，dbp)。

为实现上述目的，本发明提供了一种基于脉搏波的血压检测装置。总体结构如图1所示，包括：光电容积脉搏波(photoplethysmograph，ppg)信号抗干扰处理模块、脉搏波传导时间(pulsetransittime，ptt)提取模块和血压值计算模块；ppg信号抗干扰处理模块的输入端作为总输入端，与所采集的原始ppg信号连接，输出端与ptt提模块的输入端连接；用于去除原始ppg信号中的高频噪声、基线漂移和运动伪影，输出相对干净的ppg信号；ptt提取模块的输出端与血压值计算模块的输入端连接；用于利用干净的ppg信号提取特征参数ptt；血压值计算模块的第一输出端和第二输出端均作为总输出端；用于根据ptt进行曲线拟合，输出收缩压(systolicbloodpressure，sbp)值与舒张压(diastolicbloodpressure，dbp)值。

在本发明实施例中，ppg信号抗干扰处理模块包括：第一低通滤波模块、加权均值滤波模块、自相关运算模块、周期值计算模块和周期平滑滤波模块；第一低通滤波模块的输入端作为ppg信号抗干扰处理模块的输入端，输出端与加权均值模块的输入端连接；通过2阶5hz的iir低通滤波器去除ppg信号中的高频噪声；加权均值滤波模块的输出端分别与自相关运算模块输入端、周期平滑滤波模块的第一输入端连接；通过加权均值滤波去除ppg信号中的基线漂移；所述自相关运算模块的输出端与周期值计算模块的输入端连接；用于对去除了高频噪声和基线漂移的ppg信号进行自相关运算；自相关运算前后的信号的频率是相同的，而自相关运算后信号的频率特征非常明显，利用这一特性可确定ppg信号的周期值；所述周期值计算模块的输出端与周期平滑滤波模块的第二输入端连接；通过提取自相关运算后信号的有效波峰来确定自相关运算后信号的平均周期值，也就是自相关运算前ppg信号的周期值；周期平滑滤波的输出端作为ppg信号抗干扰处理模块的输出端；根据周期值去除ppg信号中的运动伪影。

其中，第一低通滤波模块的输入端作为ppg信号抗干扰处理模块的输入端，输出端与加权均值模块的输入端连接；第一低通滤波模块的算法流程如图2所示，具体包括：

(1)将第一ppg信号ppg1的前两个数据ppg1(1)、ppg1(2)置0：

(2)获得原始ppg信号ppg0中的第n个数据ppg0(n)，其中n的范围为1～1024；

(3)判断获得的原始ppg信号ppg0(n)是否大于等于3个，若是，则进入步骤(4)；若否，则进入步骤(2)；

(4)计算b(n)＝ppg0(n)×b(0)+ppg0(n-1)×b(1)+ppg0(n-2)×b(2)，其中b(n)为输入信号的部分和，ppg0(n-1)、ppg0(n-2)分别为前两个获得的原始ppg信号，b(0)为ppg0(n)项的滤波器系数，b(1)为ppg0(n-1)项的滤波器系数，b(2)为ppg0(n-2)项的滤波器系数；

(5)计算a(n)＝ppg1(n-1)×a(1)+ppg1(n-2)×a(2)，其中a(n)为输出信号的部分和，ppg1(n-1)、ppg1(n-2)分别为前两个获得的第一ppg信号，a(1)为ppg1(n-1)项的滤波器系数，a(2)为ppg1(n-2)项的滤波器系数；

(6)计算c(n)＝b(n)+a(n)，其中c(n)为当前输出信号的和值；

(7)获得ppg1(n)＝c(n)/a(0)，其中ppg1(n)为当前的第一ppg信号，a(0)为ppg1(n)项的滤波器系数；

(8)判断n是否小于1024，若是，则进入步骤(7)；若否，则结束运算

(9)n自加1，重复步骤(2)～(8)。

低通滤波模块算法简单，通过低通滤波能有效去除原始ppg信号中的高频噪声。

所述加权均值滤波模块的输出端分别与自相关运算模块输入端、周期平滑滤波模块的第一输入端连接；加权均值滤波模块的算法流程如图3所示，具体包括：

(1)获得第一ppg信号ppg1中的第n个数据ppg1(n)，其中n的取值范围为1～1024；

(2)判断所述信号ppg1(n)是否为ppg1信号中的第1个数据，若是，则进入步骤(3)；若否，则进入步骤(4)；

(3)计算m(n)＝ppg1(1)，其中ppg1(1)为ppg1信号中的第1个数据，m(n)为ppg1(n)的均值信号，并进入步骤(5)；

(4)计算m(n)＝{m(n-1)+[ppg1(n+63)-m(n-1)]}/64，其中m(n-1)为前一个均值信号，ppg1(n+63)为当前数据之后的第63个数据，并进入步骤(5)；

(5)获得ppg2(n)＝ppg1(n)-m(n)，其中ppg2(n)表示第二ppg信号ppg2中的第n个数据；

(6)判断n是否小于1024，若是，则进入步骤(7)；若否，则结束运算；

(7)n自加1，重复步骤(1)～(6)。

加权均值滤波模块运算量小，通过加权均值滤波能有效去除ppg信号中的基线漂移。

所述自相关运算模块的输出端与周期值计算模块的输入端连接；自相关运算模块的算法流程如图4所示，具体包括：

(1)将部分和s置0；

(2)获得第二ppg信号ppg2中的第n个数据ppg2(n)，其中n的范围为1～1024；

(3)计算s(n)＝s(n-1)+ppg2(n)·ppg2(n+m)，其中s(n)为当前部分和，s(n-1)为前一个部分和，ppg2(n+m)为当前数据之后的第m个数据；

(4)判断变量n是否小于1024，若是，则进入步骤(5)；若否，则进入步骤(6)；

(5)n自加1，重复步骤(2)～(4)；

(6)获得ppg3(m)＝s(n)/1024，其中ppg3(m)表示ppg3信号中的第m个数据，m的范围为1～1024；

(7)将部分和s置0、变量n置0；

(8)判断变量m是否小于1024，若是，则进入步骤(9)；若否，则结束运算；

(9)m自加1，重复步骤(2)～(8)；

对去除了高频噪声和基线漂移的ppg信号进行自相关运算，自相关运算前后的信号的频率是相同的，而自相关运算后信号的频率特征非常明显，利用这一特性可确定去除了高频噪声和基线漂移的ppg信号的周期值。

所述周期值计算模块的输出端与周期平滑滤波模块的第二输入端连接；周期值计算模块的算法流程如图5所示，具体包括：

(1)获得第三ppg信号ppg3中的第m个数据ppg3(m)，其中m的范围为1～1024；

(2)将当前的自相关运算后的信号ppg3(m)与前后两个自相关运算后的信号ppg3(m-1)和ppg3(m+1)分别做对比，判断所述信号ppg3(m)是否为最大值，若是，则进入步骤(3)；若否，则进入步骤(1)；

(3)判定该点为波峰，记录该点坐标；

(4)判断相邻波峰的间距是否大于等于阈值d，若是，则进入步骤(5)；若否，则进入步骤(1)；

(5)判定该波峰为有效波峰，记录有效波峰的数量p与坐标xp；

(6)计算np＝p-1，其中np为第三ppg信号ppg3的周期数量；

(7)计算sp＝xp(max)-xp(min)，其中sp为第三ppg信号ppg3的周期总和，xp(max)为最后一个有效波峰的坐标，xp(min)为第一个有效波峰的坐标；

(8)获得p＝sp/np，其中p为第三ppg信号ppg3的平均周期值，也是第二ppg信号ppg2的平均周期值。

通过提取自相关运算后信号的有效波峰可以确定自相关运算后信号的平均周期值，也就是自相关运算前的信号ppg2(n)的平均周期值。

所述周期平滑滤波的输出端作为ppg信号抗干扰处理模块的输出端；周期平滑滤波模块的算法流程如图6所示，具体包括：

(1)获得第二ppg信号ppg2的平均周期值p；

(2)获得第二ppg信号ppg2中的第n个数据ppg2(n)，其中n的范围为1～1024；

(3)计算d(n)＝ppg2(n)+ppg2(n+p)+ppg2(n+2p)+ppg2(n+3p)，其中d(n)为输入信号的部分和，ppg2(n+p)、ppg2(n+2p)、ppg2(n+3p)分别为当前数据之后的第p个、第2p个、第3p个数据；

(4)获得ppg4(n)＝d(n)/4，其中ppg4(n)为第四ppg信号；

(5)判断变量n是否小于1024，若是，则进入步骤(6)；若否，则结束运算；

(6)n自加1，重复步骤(2)～(5)。

周期平滑滤波运算简单，通过周期平滑滤波能有效去除ppg信号中的运动伪影。

在本发明实施例中，ptt提取模块包括：二阶微分运算模块、第二低通滤波模块和特征参数提取模块；二阶微分运算模块的输入端作为ptt提取模块的输入端，输出端与第二低通滤波模块连接；用于对干净的ppg信号进行二阶微分运算，获取apg信号；第二低通滤波模块的输出端与特征参数提取模块的输入端连接；二阶微分运算会放大相对干净的ppg信号中仍存在的少量高频噪声，需要通过2阶5hz的iir低通滤波器去除apg信号中的残余的高频噪声；特征参数提取模块的输出端作为ptt提取模块的输出端；通过多重定位的方式对apg信号的波峰和波谷进行筛选，提取ptt值。

其中，二阶微分运算模块的输入端作为ptt提取模块的输入端，输出端与第二低通滤波模块连接；二阶微分运算模块的算法流程如图7所示，具体包括：

(1)获得第四ppg信号ppg4中的第n个数据ppg4(n)，其中n的范围为1～1024；

(2)获得apg1(n)＝ppg4(n)+ppg4(n-2)-ppg4(n-1)×2，其中apg1(n)为第一apg信号，ppg4(n-1)、ppg4(n-2)分别为前两个第四ppg信号ppg4。

(3)判断变量n是否小于1024，若是，则进入步骤(4)；若否，则结束运算；

(4)n自加1，重复步骤(1)～(3)。

通过对去除了高频噪声、基线漂移和运动伪影之后的ppg信号进行二阶微分运算，所获取的apg信号更加准确。

所述第二低通滤波模块的输出端与特征参数提取模块的输入端连接；二阶微分运算会放大相对干净的ppg信号中仍存在的少量高频噪声，需要通过2阶5hz的iir低通滤波器去除apg信号中的残余的高频噪声，第二低通滤波模块的算法流程与第一低通滤波模块一致；

基于apg信号进行ptt提取流程如图8所示；具体包括：b为主波谷，a、c分别为主波谷两侧最近的波峰，a与c的时间差即为所求的ptt值。

所述特征参数提取模块的输出端作为ptt提取模块的输出端；特征参数提取模块的算法流程如图9所示，具体包括：

(1)获得第二apg信号apg2中的第n个数据apg2(n)，其中n的范围为1～1024；

(2)将当前第二apg信号apg2(n)与前后两个第二apg信号apg2(n-1)和apg2(n+1)分别做对比，判断所述信号apg2(n)是否为最大值、是否为最小值，若是最大值，则进入步骤(3)；若是最小值，则进入步骤(4)；若既不是最大值也不是最小值，则进入步骤(1)；

(3)判定该点为波峰，记录该点坐标；

(4)判定该点为波谷，记录该点坐标；

(5)判断相邻波谷的间距是否大于等于阈值l且波谷幅值小于等于阈值m，若满足了间距阈值、幅值阈值的条件，则进入步骤(6)；若不满足间距阈值、幅值阈值的条件，则进入步骤(1)；

(6)判定为有效波谷，记录其坐标与数量t；

(7)根据步骤(5)、(6)判断波峰是否为有效波谷两侧最近的波峰，若是，则进入步骤(8)；若否，则进入步骤(1)

(8)判定波峰是有效波谷两侧最近的波峰，记录两波峰的坐标差δxp；

(9)计算sptt＝sptt+δxp，其中sptt为ptt的总和；

(10)计算nptt＝t，其中nptt为ptt的数量；

(11)获得ptt＝sptt/nptt，其中ptt为平均ptt值。

通过多重定位的方式对apg信号的波峰和波谷进行筛选，所提取的ptt值更加准确。

在本发明实施例中，血压值计算模块包括：sbp计算模块和dbp计算模块；sbp计算模块的输入端作为血压值计算模块的输入端，输出端作为血压值计算模块的第一输出端；根据拟合公式sbp＝a*ptt+b获取收缩压sbp，其中a、b分别为线性公式的一次项系数和常数项系数；

dbp计算模块的输入端作为血压值计算模块的输入端，输出端作为血压值计算模块的第一输出端；根据拟合公式dbp＝c*ptt+d+e*ptt^-2获取舒张压dbp，其中c、d、e分别为线性公式的一次项系数、常数项系数和负二次项系数。

本发明还提供了一种基于脉搏波的血压检测方法，包括下述步骤：

s1：去除原始ppg信号中的高频噪声、基线漂移和运动伪影后，获得相对干净的ppg信号；

s2：利用干净的ppg信号提取特征参数ptt；

s3：将特征参数ptt进行曲线拟合，并获得sbp值与dbp值。

其中，在步骤s1中，将原始ppg信号依次经过低通滤波处理、加权均值滤波处理、自相关运算、周期值计算以及周期平滑滤波处理后获得所述相对干净的ppg信号。

作为本发明的一个实施例，在步骤s1中，所述加权均值滤波处理方法具体为：

(1)获得第一ppg信号ppg1中的第n个数据ppg1(n)，其中n的范围为1～1024；

(2)判断所述信号ppg1(n)是否为ppg1信号中的第1个数据，若是，则进入步骤(3)；若否，则进入步骤(4)；

(3)计算m(n)＝ppg1(1)，其中ppg1(1)为ppg1信号中的第1个数据，m(n)为ppg1(n)的均值信号，并进入步骤(5)；

(4)计算m(n)＝{m(n-1)+[ppg1(n+63)-m(n-1)]}/64，其中m(n-1)为前一个均值信号，ppg1(n+63)为当前数据之后的第63个数据，并进入步骤(5)；

(5)获得ppg2(n)＝ppg1(n)-m(n)，其中ppg2(n)表示第二ppg信号ppg2中的第n个数据；

(6)判断n是否小于1024，若是，则进入步骤(7)；若否，则结束运算；

(7)n自加1，重复步骤(1)～(6)。

本发明中，加权均值滤波模块运算量小，通过加权均值滤波能有效去除ppg信号中的基线漂移。

作为本发明的一个实施例，在步骤s1中，所述周期值计算方法具体包括：

(1)获得第三ppg信号ppg3中的第m个数据ppg3(m)，其中m的范围为1～1024；

(2)将当前的自相关运算后的信号ppg3(m)与前后两个自相关运算后的信号ppg3(m-1)和ppg3(m+1)分别做对比，判断所述信号ppg3(m)是否为最大值，若是，则进入步骤(3)；若否，则进入步骤(1)；

(3)判定该点为波峰，并记录该点坐标；

(4)判断相邻波峰的间距是否大于等于阈值d，若是，则进入步骤(5)；若否，则进入步骤(1)；

(5)判定该波峰为有效波峰，并记录有效波峰的数量p与坐标xp；

(6)计算np＝p-1，其中np为第三ppg信号ppg3的周期数量；

(7)计算sp＝xp(max)-xp(min)，其中sp为第三ppg信号ppg3的周期总和，xp(max)为最后一个有效波峰的坐标，xp(min)为第一个有效波峰的坐标；

(8)获得p＝sp/np，其中p为第三ppg信号ppg3的平均周期值，也是第二ppg信号ppg2的平均周期值。

本发明通过提取第三ppg信号ppg3的有效波峰可以确定ppg3的平均周期值，也就是第二ppg信号ppg2的平均周期值。

作为本发明的一个实施例，在步骤s2中，所述提取特征参数ptt的方法具体为：

(1)获得第二apg信号apg2中的第n个数据apg2(n)，其中n的范围为1～1024；

(2)将当前第二apg信号apg2(n)与前后两个第二apg信号apg2(n-1)和apg2(n+1)分别做对比，判断所述信号apg2(n)是否为最大值、是否为最小值，若是最大值，则进入步骤(3)；若是最小值，则进入步骤(4)；若既不是最大值也不是最小值，则进入步骤(1)；

(3)判定该点为波峰，记录该点坐标；

(4)判定该点为波谷，记录该点坐标；

(5)判断相邻波谷的间距是否大于等于阈值l且波谷幅值小于等于阈值m，若满足了间距阈值、幅值阈值的条件，则进入步骤(6)；若不满足间距阈值、幅值阈值的条件，则进入步骤(1)；

(6)判定为有效波谷，记录其坐标与数量t；

(7)根据步骤(5)、(6)判断波峰是否为有效波谷两侧最近的波峰，若是，则进入步骤(8)；若否，则进入步骤(1)

(8)判定波峰是有效波谷两侧最近的波峰，记录两波峰的坐标差δxp；

(9)计算sptt＝sptt+δxp，其中sptt为ptt的总和；

(10)计算nptt＝t，其中nptt为ptt的数量；

(11)获得ptt＝sptt/nptt，其中ptt为平均ptt值。

本发明通过多重定位的方式对apg信号的波峰和波谷进行筛选，所提取的ptt值更加准确。

本发明通过第一低通滤波模块去除ppg信号中的高频噪声，通过加强均值滤波模块去除ppg信号中的基线漂移，通过自相关运算模块、周期值计算模块与周期平滑滤波模块去除ppg信号中运动伪影，整体运算量小，且具有较好的抗干扰处理效果；本发明还通过二阶微分运算模块产生加速度apg信号，通过第二低通滤波模块去除apg信号中的残余高频噪声，通过特征参数提取模块提取特征参数ptt，进一步的去噪处理与多重定位的提取方式提高了提取ptt的准确度；本发明最终通过曲线拟合来获取sbp值与dbp值，实现简单。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。