一种基于相关性的脉搏波传导时间获取方法与流程

本发明涉及一种脉搏波传导时间获取方法。特别是涉及一种基于相关性的脉搏波传导时间获取方法。

背景技术：

脉搏波传导时间是人体动脉脉压波在动脉血管中的传导时间。脉搏波传导时间的获取方法一般有两大类，一种是常规方法，计算不同位置的两路脉搏波信号，计算脉搏波传导时间，另一种测量方法是利用心电信号和脉搏波信号测量脉搏波传导时间。

常规测量方法就是在脉搏波传播过程中，同步地测量两个不同位置的动脉脉搏波，通过测量某个相对应点的传播时间得到脉搏波传导时间，但由于脉搏波是一个具有多重次谐波组合而成的复合波，其波形复杂并且在传播过程中又在不停的发生变化，所以常规算法并不能保证测量的准确性。

目前常用的脉搏波传导时间测量方法的就是通过心电信号和光电脉搏波信号的时间差来获得脉搏波传导时间。其主要是选择一些特征点，如选择人体心电信号的r波峰值点与光电脉搏波信号的某些特征点作为脉搏波传导时间的起点与终点，从而获得脉搏波传导时间。该方法可能由于个别特征点的无法准确识别，导致最后结果的不准确性。基于相关性的脉搏波传导时间获取方法，则是选取心电信号和光电脉搏波信号的一段数据，避免了由于个别点引起的误差。基于相关性的脉搏波传导时间获取方法的关键在于计算选取片段内心电信号和光电脉搏波信号的相关系数，查找相关系数最大值时光电脉搏波信号的平移量。

技术实现要素：

本发明所要解决的技术问题是，提供一种能够准确获取脉搏波传导时间的基于相关性的脉搏波传导时间获取方法。

本发明所采用的技术方案是：一种基于相关性的脉搏波传导时间获取方法，包括如下步骤：

1)利用电极片和光电脉搏波传感器同步采集人体心电信号和光电脉搏波信号；

2)确定在设定长度中人体心电信号数据与光电脉搏波信号数据的相关系数公式：

式中，r0为相关系数，设定长度的起点为loc_1，设定长度的终点为loc_2，xi为人体心电信号在起点loc_1和终点loc_2之间的任意一个值，ex为人体心电信号在起点loc_1和终点loc_2之间的均值，yi为光电脉搏波信号在起点loc_1和终点loc_2之间的任意一个值，ey为光电脉搏波信号在起点loc_1和终点loc_2之间的均值，n＝loc_2-loc_1；

3)在时域内，将光电脉搏波信号向左移动1个单位的偏移量，即ppg(loc_1+1，loc_2+1)，利用相关系数公式计算偏移后ecg(loc_1，loc_2)与ppg(loc_1+1，loc_2+1)的相关系数，记为r1，其中，ecg为人体心电信号，ppg为光电脉搏波信号；

4)重复步骤3)，计算ecg(loc_1，loc_2)与ppg(loc_1+i，loc_2+i)的相关系数，直至得到m个相关系数：r1、r2……rm，其中m≥i；

5)选择m个相关系数r1、r2……rm中最大的相关系数ri，将相关系数ri所对应的光电脉搏波移动的偏移量用i表示，从而得到脉搏波传导时间ptt，ptt＝i×0.5×10³。

步骤2)包括：找到人体心电信号的一个r波，以r波为起点loc_1，选取一段设定长度的数据，终点位置记为loc_2，设定xi为人体心电信号在起点loc_1和终点loc_2之间的任意一个值，ex为人体心电信号在起点loc_1和终点loc_2之间的均值，yi为光电脉搏波信号在起点loc_1和终点loc_2之间的任意一个值，ey为光电脉搏波信号在起点loc_1和终点loc_2之间的均值，n＝loc_2-loc_1，所述相关系数记为r0，则：

在matlab软件中使用r0＝corrcoef(ecg(loc_1，loc_2)，ppg(loc_1，loc_2))。

本发明的一种基于相关性的脉搏波传导时间获取方法，具有如下有益效果：

1、本发明选取的心电信号和光电脉搏波一段区域，而不是单个的特征点计算脉搏波传导时间，提高了脉搏波传导时间测量的准确性；

2、本发明在加入高斯白噪声和基线漂移的情况下同样具有较好的效果，具有更大的潜在应用前景。

附图说明

图1是本发明的一种基于相关性的脉搏波传导时间获取方法流程图；

图2是本发明中脉搏波信号向左平移示意图。

具体实施方式

下面结合实施例和附图对本发明的一种基于相关性的脉搏波传导时间获取方法做出详细说明。

本发明的一种基于相关性的脉搏波传导时间获取方法，是一种通过相关性分析采集的人体心电信号和光电脉搏波信号获取脉搏波传导时间ptt的方法。

如图1所示，本发明的一种基于相关性的脉搏波传导时间获取方法包括如下步骤：

1)利用电极片和光电脉搏波传感器同步采集人体心电信号(ecg)和光电脉搏波信号(ppg)；

2)确定在设定长度中人体心电信号数据与光电脉搏波信号数据的相关系数公式：

式中，r0为相关系数，设定长度的起点为loc_1，设定长度的终点为loc_2，xi为人体心电信号在起点loc_1和终点loc_2之间的任意一个值，ex为人体心电信号在起点loc_1和终点loc_2之间的均值，yi为光电脉搏波信号在起点loc_1和终点loc_2之间的任意一个值，ey为光电脉搏波信号在起点loc_1和终点loc_2之间的均值，n＝loc_2-loc_1；

相关系数公式的确定包括：找到人体心电信号的一个r波，以r波为起点loc_1，选取一段设定长度的数据，终点位置记为loc_2，设定xi为人体心电信号在起点loc_1和终点loc_2之间的任意一个值，ex为人体心电信号在起点loc_1和终点loc_2之间的均值，yi为光电脉搏波信号在起点loc_1和终点loc_2之间的任意一个值，ey为光电脉搏波信号在起点loc_1和终点loc_2之间的均值，n＝loc_2-loc_1，所述相关系数记为r0，则：

在matlab软件中使用r0＝corrcoef(ecg(loc_1，loc_2)，ppg(loc_1，loc_2))。

3)在时域内，将光电脉搏波信号向左移动1个单位的偏移量，即ppg(loc_1+1，loc_2+1)，利用相关系数公式计算偏移后ecg(loc_1，loc_2)与ppg(loc_1+1，loc_2+1)的相关系数，记为r1，其中，ecg为人体心电信号，ppg为光电脉搏波信号；

4)重复步骤3)，计算ecg(loc_1，loc_2)与ppg(loc_1+i，loc_2+i)的相关系数，直至得到m个相关系数：r1、r2……rm，其中m≥i；

5)选择m个相关系数r1、r2……rm中最大的相关系数ri，将相关系数ri所对应的光电脉搏波移动的偏移量用i表示，从而得到脉搏波传导时间ptt，ptt＝i×0.5×10³。

下面给出实例：

1、利用analogdiscovery采集卡同步采集心电信号(ecg)和光电脉搏波信号(ppg)，采样频率为2khz，将数据(长度为8000)上传到pc机上；

2、在ecg信号中找到一个r波，然后选取一段数据，同样在ppg信号中选取同样长度的数据，计算ecg信号与ppg信号的相关系数；

3、将光电脉搏波信号向左移动1个单位(0.5×10^-3s)，然后计算与上述心电信号的相关系数，判断是否为最大值；

4、如图2所示，依次向左移动1,2,3，……，m个单位，分别计算心电信号和光电脉搏波的相关系数；

5、如图2所示，选择相关系数最大的一个值，记录此时的平移量i(i≤m)，由平移量计算此时的ptt，ptt＝i×0.5×10³。