具体实施方式

下面将结合附图和实施例，对本发明技术方案做优选详细描述，显然所描述的实施例是本发明一部分实施例，而不是全部的实施例，本发明的实施方式并不限于此。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例

如图1所示，一种基于ECG信号和PPG信号的血压监测系统，包括生理信号采集模块、数据传输模块、数据预处理模块、特征提取模块、异常值检测模块、特征选择模块、血压回归预测模块。通过部署在手表或手环装置的生理信号采集模块，同步采集佩戴者的肢体导联的心电图(ElectroCardioGram)ECG信号和光电容积脉搏波(PhotoPlethysmoGraphy)PPG信号，使用简单、方便携带、可随时进行血压监测；通过数据预处理模块对生理信号数据进行预处理；特征提取模块提取ECG信号和PPG信号的特征，数据异常值检测模块根据提取的生理信号特征进行异常值检测，并剔除相关的异常样本，特征选择模块选择出与血压相关的生理信号特征，血压回归预测模块使用迁移学习算法结合建立的数据库训练用户的个人血压预测模型，对提取的生理信号特征进行血压预测，使血压预测模型更加适配于每个用户，确保血压监测设备的精度。

生理信号采集模块，包括ECG传感器和PPG传感器，用于同步采集ECG信号、PPG信号这两种生理信号数据，得到采样频率为250Hz的ECG信号和PPG信号。优选地，ECG传感器和PPG传感器集成于一台智能手表或手环装置上，ECG传感器采用的是肢体单导联三电极，其中两个电极集成在智能手表或手环装置底部，另外一个电极集成在智能手表或手环装置右侧；PPG传感器采用的是光电传感器，集成于智能手表或手环装置底部，位于智能手表或手环装置上ECG传感器底部双电极的中央，将智能手表或手环装置佩戴于左手手腕处，开启采集生理信号数据时，右手食指轻触智能手表或手环装置右侧电极，智能手表或手环装置的主控程序开启250Hz的定时器中断，通过信号放大电路、模数转换电路等读取传感器中的ECG信号和PPG信号，同步采集采样频率为250Hz的ECG信号和PPG信号。

数据传输模块，用于将生理信号数据ECG信号和PPG信号通过蓝牙通信从传感器端传输到手机移动端，手机移动端再通过网络将生理信号数据上传至服务器。本实施例中，部署在智能手表端的蓝牙模块与手机移动端的蓝牙相连接，同步更新智能手表端上的时间信息，再将采集到的生理信号数据打包传输到手机移动端。数据包中以一秒的数据为一帧，每一帧数据均携带手表上的时间信息，然后手机移动端再通过网络将数据上传至服务器，减少手表终端的计算量同时能够保证整个血压预测模型的运算速度。

数据预处理模块，用于对采集到的ECG信号和PPG信号进行时间对齐、数据切分、信号去噪、归一化处理、峰值检测、数据剔除操作等预处理，剔除不满足条件的生理信号数据片段。如图2所示，数据预处理模块整体流程，数据传输模块将生理信号数据包传送至手机移动端，生理信号数据包中以1秒为一帧，帧头为当前时间信息，数据部分则为采集到的ECG信号和PPG信号。通过生理信号数据包中的时间信息，对齐两种生理信号，然后以T秒间隔进行数据切分，优选地T为5。ECG信号去噪首先设计200ms和600ms的中值滤波器，把心电信号的QRS波群、P波和T波滤除，得到ECG信号基线，在将原始信号和基线相减，得到消除基线漂移后的心电信号；然后设计50Hz的数字陷波器去除工频噪声干扰；最后设计截至频率为45Hz的低通滤波器滤除肌电噪声的影响，得到纯净的心电信号。PPG信号去噪首先使用小波变换分解出噪声细节分量，再通过小波逆变换进行信号重构，从而消除基线漂移，然后设计50Hz的数字陷波器去除工频噪声干扰以及设计截止频率为0.5Hz-4Hz的带通滤波器滤除高频噪声。先对滤波后的ECG信号和PPG信号进行归一化处理，然后再进行峰值检测，设置峰值检测的幅值阈值以及时间间隔阈值，避免误识别极大值。根据峰值检测结果进行数据剔除，由于正常静息状态下同时采集的ECG信号和PPG信号波峰是交替出现的，因此在5秒的生理信号片段中，检查ECG和PPG信号波峰的数量以及波峰位置是否交替出现，不满足条件的生理信号数据片段将被剔除。

特征提取模块，用于对ECG信号和PPG信号进行特征提取，提取的特征包括ECG信号的时频域特征、PPG信号的时频域特征以及两种生理信号的在时间上的联系，进行多特征融合以全面提取与血压相关的生理信息。其中，ECG信号的时域特征包括提取ECG信号一阶小波变换后获得的细节信号和近似信号的方差、均值、最大值、最小值、排列熵、样本熵，频域特征包括频带能量和频带能量占比；PPG信号的时域特征包括提取PPG原始信号的心率特征、脉冲与脉冲之间的皮尔逊相关系数、PPG原始信号、一阶导信号、二阶导信号的最大振幅、最小振幅、周期、最小振幅/最大振幅、脉冲幅度为10％、25％、33％、50％、66％、75％、100％时信号的收缩期时间、舒张期时间、收缩期时间/舒张期时间、收缩期时间+舒张期时间、收缩期信号面积、舒张期信号面积、收缩期信号面积或舒张期信号面积、收缩期信号面积和舒张期信号面积；PPG信号经过一阶、二阶和三阶小波变换后获得的细节信号和近似信号的方差、均值、最大值、最小值、排列熵、样本熵、频带能量以及频带能量占比；以及PPG原始信号的基频、基波振幅、二次谐波频率、二次谐波振幅、三次谐波频率、三次谐波振幅。两种生理信号在时间上的联系指同步采集的ECG和PPG波峰之间的时间差，称为脉冲到达时间。

具体地、如图3和图4所示，数据预处理模块将生理信号切分成5秒一段的生理信号数据片段，提取生理信号数据片段的如下特征参数：

(1)ECG时频域特征参数包括：小波变换获得的细节信号(cD)、近似信号(cA)的最小值(E_Min_cD1、E_Min_cD2、E_Min_cD3、E_Min_cA3)、最大值(E_Max_cD1、E_Max_cD2、E_Max_cD3、E_Max_cA3)、方差(E_Var_cD1、E_Var_cD2、E_Var_cD3、E_Var_cA3)、均值(E_Mean_cD1、E_Mean_cD2、E_Mean_cD3、E_Mean_cA3)、频带能量(E_Energy_cD1、E_Energy_cD2、E_Energy_cD3、E_Energy_cA3)、频带能量占比(E_Ratio_Energy_cD1、E_Ratio_Energy_cD2、E_Ratio_Energy_cD3、E_Ratio_Energy_cA3)、m＝1样本熵(E_SampEn_1_cD1、E_SampEn_1_cD2、E_SampEn_1_cD3、E_SampEn_1_cA3)、m＝2样本熵(E_SampEn_2_cD1、E_SampEn_2_cD2、E_SampEn_2_cD3、E_SampEn_2_cA3)、排列熵(E_PE_cD1、E_PE_cD2、E_PE_cD3、E_PE_cA3)；

(2)如图5所示，PPG时频域特征参数包括：脉冲与脉冲之间的皮尔逊相关系数(pearson)、最大振幅(Ih)、脉冲最小振幅(Il)、脉冲信号周期(T_pi)、最小振幅/最大振幅(Il/Ih)、脉冲幅度为10％、25％、33％、50％、66％、75％、100％时收缩期时间(St10、St25、St33、St50、St66、St75、St)、舒张期时间(Dt10、Dt25、Dt33、Dt50、Dt66、Dt75、Dt)、收缩期和舒张期之比(St10/Dt10、St25/Dt25、St33/Dt33、St50/Dt50、St66/Dt66、St75/Dt75、、St/Dt)、收缩期和舒张期之和(St10+Dt10、St25+Dt25、St33+Dt33、St50+Dt50、St66+Dt66、St75+Dt75、St+Dt)、收缩期面积(A_St10、A_St25、A_St33、A_St50、A_St66、A_St75、、A_St)、舒张期面积(A_Dt10、A_Dt25、A_Dt33、A_Dt50、A_Dt66、A_Dt75、A_Dt)、收缩期和舒张期面积之比(A_St10/A_Dt10、A_St25/A_Dt25、A_St33/A_Dt33、A_St50/A_Dt50、A_St66/A_Dt66、A_St75/A_Dt75、A_St/A_Dt)、收缩期和舒张期面积之和(A_St10+A_Dt10、A_St25+A_Dt25、A_St33+A_Dt33、A_St50+A_Dt50、A_St66+A_Dt66、A_St75+A_Dt75、A_St+A_Dt)；一阶、二阶导信号峰值(Ih_1d、Ih_2d)、谷值(Ih_1d、Ih_2d)、谷值/峰值(Il_1d/Ih_1d、Il_2d/Ih_2d)、上升期时间(St_1d、St_2d)、下降期时间(Dt_1d、Dt_2d)、从开始到谷值时间(St_1d+Dt_1d、St_2d+Dt_2d)、下降期时间/上升期时间(Dt_1d/St_1d、Dt_2d/St_2d)、上升期面积(A_St_1d、A_St_2d)、下降期面积(A_Dt_1d、A_Dt_2d)、上升期和下降期面积之和(A_St_1d+A_Dt_1d、A_St_2d+A_Dt_2d)、上升期和下降期面积之比(A_St_1d/A_Dt_1d、A_St_2d/A_Dt_2d)、基频(f_base)、基波振幅(sp_mag_base)、二次谐波频率(f_2)、二次谐波振幅(sp_mag_2)、三次谐波频率(f_3)、三次谐波振幅(sp_mag_3)、心率(HR)、小波变换获得的细节信号(cD)和近似信号(cA)的最小值(P_Min_cD1、P_Min_cD2、P_Min_cD3、P_Min_cA3)、最大值(P_Max_cD1、P_Max_cD2、P_Max_cD3、P_Max_cA3)、方差(P_Var_cD1、P_Var_cD2、P_Var_cD3、P_Var_cA3)、均值(P_Mean_cD1、P_Mean_cD2、P_Mean_cD3、P_Mean_cA3)、频带能量(P_Energy_cD1、P_Energy_cD2、P_Energy_cD3、P_Energy_cA3)、频带能量占比(P_Ratio_Energy_cD1、P_Ratio_Energy_cD2、P_Ratio_Energy_cD3、P_Ratio_Energy_cA3)、m＝1样本熵(P_SampEn_1_cD1、P_SampEn_1_cD2、P_SampEn_1_cD3、P_SampEn_1_cA3)、m＝2样本熵(P_SampEn_2_cD1、P_SampEn_2_cD2、P_SampEn_2_cD3、P_SampEn_2_cA3)、排列熵(P_PE_cD1、P_PE_cD2、P_PE_cD3、P_PE_cA3)；

优选地，通过计算脉冲与脉冲之间的皮尔逊相关系数，获得PPG信号在某一时间段内脉冲与脉冲之间的相关性矩阵，从而反应短时间内PPG信号脉冲的变化程度。从5秒的PPG信号提取n个脉冲周期S₁、S₂、S₃…S_n，脉冲与脉冲之间的计算公式如下：

其中，S_i、S_j表示5秒PPG信号片段中分割出来的两个不同的脉冲信号周期，脉冲与脉冲之间的相关性矩阵如下：

计算相关性矩阵的中所有元素的均值、方差、最小值，均值反应了脉冲与脉冲之间相关性的平均水平，方差反应了脉冲与脉冲之间相关性的变化程度，最小值反应了脉冲与脉冲之间变化的最大程度。

(3)如图6所示，ECG信号和PPG信号之间时间上的相关的特征参数为ECG信号波峰和PPG信号波峰之间的时间差(PAT)。

异常值检测模块，用于检测并剔除在特征空间中分布在比较稀疏区域且距离密度高的群体比较远的异常样本，使得剔除后纯净的样本数据更能表征当前用户的血压信息。通过孤立森林算法对样本在特征空间中进行异常值检测，在特征空间中分布在比较稀疏区域且距离密度高的群体比较远的样本点被认为是异常样本，孤立异常样本并将异常样本剔除。所采集的4分钟生理信号数据切分成5秒一段的生理信号数据片段，然后进行滤波以及根据波峰位置、波峰大小、以及波峰数量剔除不合格样本，剩余的样本数据片段依然存在着难以滤除的强噪声和运动伪迹，从这些样本中提取的特征无法表征当前用户的血压情况，因此，通过孤立森林算法将切分的样本点映射至高维度的特征空间中，根据随机选择的每一维特征切分出不同的子空间，直到子空间中只有一个样本点，计算每个样本点的异常分数s，设置异常样本比例值，优选异常样本比例值为0.2，最后剔除在特征空间中分布在比较稀疏区域且距离密度高的群体比较远的异常样本。如图7所示，实施例中孤立森林算法效果图。

优选地，通过孤立森林算法对样本在特征空间中进行异常值检测的具体步骤包括：

步骤1：随机从剩余的样本数据片段中选取X个样本点并将其映射至高维度的特征空间中；

步骤2：随机选取一维特征，通过设置随机阈值，产生一个超平面，将特征空间中的样本数据切分为两个子空间，即产生一个孤立树的两个分支；

步骤3：在切割出来的子空间中不断重复步骤2，直到切分的子空间中只有一个样本点或者孤立树的分支数量到达上限；

步骤4：产生多个孤立树，综合每棵树的结果，计算每个剩余样本数据片段的异常分数s，计算公式如下：

其中n(x)为每棵树中节点的个数，C(X)为X个样本的平均路径。

步骤5：设置异常值分数阈值，当样本的异常分数评分超过该阈值时，则该样本将被视为异常样本，从剩余的样本数据片段中剔除。

特征选择模块，用于通过遗传算法搜索全局最优特征子集，筛选出与血压相关的生理信号特征以及剔除冗余特征，减少血压回归模型的计算量同时提高血压回归预测精度，主要包括以下步骤：初始化种群数量；设置适应度函数；随机从种群中选出特定数量的个体；交换个体之间的染色体，形成新的子代个体；按照一定概率对子代个体的染色体进行变异。

该遗传算法基于sklearn-genetic0.2工具包进行开发，设置参数如下：man_population＝100，crossover_proba＝0.5，mutation_proba＝0.3，n_generations＝100；选择的特征数量为5至40个；使用5折交叉验证的方法避免选择的特征具有局限性和特殊性。

血压回归预测模型模块，用于结合已建立的生理信号血压数据库，使用迁移学习算法Tradaboost针对每个用户进行个性化模型训练，结合线性回归算法对选取的生理信号特征进行血压回归预测，输出其对应的收缩压(SBP)和舒张压(DBP)的预测值；通过调整模型参数，而使得血压预测模型更加适配于每一位用户。

如图8所示，使用迁移学习算法Tradaboost针对每个用户进行个性化模型训练的具体步骤包括：

步骤1、建立群体的生理信号血压数据库，在隔音测试间中进行数据采集，招募受试者，采集受试者连续五天中8:00-10:00、14:00-16:00、18:00-20:00三个时段的真实血压以及对应的4分钟生理信号，作为源域数据；

步骤2、使用电子血压计和生理信号采集模块采集一天中8:00-10:00、14:00-16:00、18:00-20:00三个时段的真实血压以及对应的4分钟生理信号作为目标域数据；

步骤3、使用Tradaboost迁移学习算法，在训练过程中，通过计算源域数据中血压预测值和血压真实值的误差更新对应的样本权重，预测值与真实值越接近，其对应权重越大；反之，权重越小；通过不断降低源域数据中与目标域数据相矛盾的样本数据权重，从而过滤源域中与目标域不匹配的数据，使得模型更加适配于每位用户。

优选地，调整模型参数的具体训练步骤包括：

(1)、初始化单个样本权重，其计算公式如下：

其中，i代表单个样本，n代表源域样本数量，m代表目标域样本数量，则所有样本权重向量为：

(2)、使用回归算法得到每一个生理信号样本的血压预测值y_p，利用生理信号样本的真实血压值y_i，计算误差公式如下：

其中，i代表单个样本，D_t代表样本的真实血压值和血压预测值y_p的最大绝对误差；

(3)、更新样本权重，公式如下：

其中，i代表单个样本，n代表源域样本数量，m代表目标域样本数量，为真实血压值和血压预测值的误差，Z_t为归一化常数，N为设置的迭代次数；

(4)、迭代次数N为10，重复步骤(2)、(3)，不断更新样本权重，预测值与真实值越接近，其对应权重越大；反之，权重越小；通过不断降低源域中与目标域数据相矛盾的样本数据权重，从而过滤源域中与目标域不匹配的数据，使得训练的模型更加适配于每位用户。

用户在首次使用该血压监测系统进行血压监测时，需提供其一天中8:00-10:00、14:00-16:00、18:00-20:00三个时段测量到的真实血压数据以及4分钟的生理信号数据进行个人模型校准，校准后的模型在进行血压预测时，只需获取手表同时采集的5秒生理信号，即可预测出用户当前的血压值。本发明提高了血压监测设备的便携性的同时也能确保血压监测的精度。

基于相同的发明构思，本实施例还提出一种基于ECG信号和PPG信号的血压监测方法，包括以下步骤：

步骤1、生理信号数据采集，同步采集用户的ECG信号和PPG信号；

步骤2、数据传输，将所述ECG信号、PPG信号打包传输到手机移动端，手机移动端将数据上传至服务器；

步骤3、数据预处理，对采集的生理信号数据进行数据切分、时间对齐、信号去噪、归一化、峰值检测、数据剔除操作等预处理；

步骤4、特征提取，对数据预处理后的生理信号数据提取ECG信号和PPG信号的特征；

步骤5、数据异常值检测，通过孤立森林算法根据提取的生理信号特征进行异常值检测，并剔除相关的异常样本，剔除不满足条件的数据片段；

步骤6、特征选择，选择出与血压相关的生理信号特征；

步骤7、血压回归预测，结合生理信号血压数据库，使用迁移学习Tradaboost算法针对每个用户进行血压预测模型训练，结合线性回归算法对选取的生理信号特征进行血压回归预测，输出收缩压和舒张压的预测值；调整血压预测模型参数，使得血压预测模型适配于每一位用户。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。