CN113080993B - 一种基于改进多基分类器集成的自动心律失常分析方法 - Google Patents

Abstract

本发明公开了一种基于改进多基分类器集成的自动心律失常分析方法，包括以下步骤：预处理模块，包括心电信号降噪、心拍识别及截取；分类模块，包括用以对心电图数据进行心律失常自动分类的XGBoost分类模块、GBDT分类模块、BiGRU分类模块；融合模块，用以依据模型集成规则—Stacking，将分类模块产生的第一分类结果、第二分类结果及第三分类结果进行集成，得到集成数据；然后利用浅层神经网络对心律失常集成数据进行分类，得到最终的心率失常分类融合结果。如此，本发明将XGBoost、GBDT、BiGRU三种基分类器的分类结果利用Stacking集成规则进行融合，采用机器学习模型与深度学习模型相结合的方法，提升了心律失常整体分类性能及准确率。

Description

一种基于改进多基分类器集成的自动心律失常分析方法

技术领域

本发明涉及医学信号处理技术领域，具体为一种基于改进多基分类器集成的自动心律失常分析方法。

背景技术

心律指心跳的节奏。正常人的心脏跳动频率每分钟约为60～100次，并且健康的心律应该是十分均匀的，心脏病或心脏神经调节功能不正常时，可出现心律不齐或心律失常。心律失常是心血管疾病常见且最为严重的病症之一，可导致患者突然死亡，严重威胁人类健康。

如今，随着人工智能技术的发展，智能辅助诊断系统因其快速、可靠的特点已成为最具前景的临床诊断解决方案。由于心电远程监护技术的迅速发展和普及，以及当下可穿戴设备的流行我们可以很容易地得到病人的心电信号以及其他生理特征。传统上，医生可以通过观察心电信号中蕴含的信息来判断病人是否患心律失常，然而仅通过观察变化微弱的心电信号来推断病人状态的过程是费力的，并且容易因复杂的心电变化而产生错误的判断。因此，利用计算机来进行智能辅助的诊断系统可以克服对心电信号的评估限制，从而更好的辅助医生对病人心律失常的患病情况进行推断。

发明内容

(一)解决的技术问题

针对现有技术的不足，本发明提供了一种基于改进多基分类器集成的自动心律失常分析方法，提升心律失常的检测效果。

(二)技术方案

为实现以上目的，本发明通过以下技术方案予以实现：一种基于改进多基分类器集成的自动心律失常分析方法，包括以下步骤：

S1：数据预处理，读入12导联的心电信号的数据，采用小波变换滤除原始心电信号中的工频干扰，采用巴特沃斯带阻滤波器滤除肌电干扰，采用中值滤波滤除原始心电信号中的基线漂移；

S2：心拍数据提取，通过小波变换检测R波位置，进而对R波前后数据点进行固定长度截取，从而完成对心拍数据提取；

S3：基分类器模型构建，包括用以对心电图数据进行心律失常自动分类的XGBoost分类模块、GBDT分类模块、BiGRU分类模块；其中，XGBoost作为第一基分类模块、GBDT作为第二基分类模块、GRU作为第三基分类模块，通过对训练数据进行K折交叉的划分，使得XGBoost输出K个第一分类结果，GBDT输出K个第二分类结果，BiGRU输出K个第三分类结果；

S4：模型融合，依据模型集成规则—Stacking，将分类模块产生的第一分类结果、第二分类结果及第三分类结果进行集成，得到集成数据，然后搭建好浅层神经网络以对心律失常集成数据进行分类；

S5：训练模型参数，初始化所述集成模型的参数，将构建好的数据集拆分为训练集和测试集；利用训练集样本对集成模型进行训练，生成所述集成模型的参数并保存；

S6:模型预测，对测试集样本进行自动识别，获得心律失常的二分类结果。

优选的，步骤S2具体包括：

读取降噪后12导联的心电信号数据，对每个导联心电信号根据同一时刻R波峰值处的位置固定向前截取200个样本点，向后截取400个样本点，使得每个导联的每个心拍固定截取600个点的数据，对相同时刻R波顶点每个导联所截取的600个点的心电信号进行行级上的拼接，此时原始每个心拍的心电信号的数据形成12*600维的样本，作为基分类模型的输入数据X。

优选的，步骤S3具体包括：

XGBoost分类模块，其作为第一分类模快，通过对训练数据进行3折交叉的划分，使得XGBoost输出3个第一分类结果，其中，对每个心电样本的预测表示为：

其中x_i∈R^12*600是第i个心电训练样本，f_k表示第k棵树，

表示样本x_i的预测结果。在建模过程中不断优化目标函数，使其最小，表示为：

其中

是样本训练误差，Ω(f_i)是树的正则项，XGBoost通过将Obj^t使用二阶泰勒公式展开得：

其中，g_i、h_i表示在第i个样本下泰勒展开式的一阶导数、二阶导数，通过公式8可以得到Obj^t关于w_j的一元二次方程，由此找到使Obj^t取值最小的w_j，表示为：

其中，G_j＝∑g_i；H_j＝∑h_i，通过上述步骤完成第t棵树的建立，从而完成XGBoost模型的构建，并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^XGBoost、A₂ ^XGBoost、A₃ ^XGBoost和对应的针对测试集的预测结果B₁ ^XGBoost、B₂ ^XGBoost、B₃ ^XGBoost；

GBDT分类模块，其作为第二分类模块，通过对训练数据进行3折交叉的划分，使得GBDT输出3个第二分类结果，表示为：

使用负二项对数似然作为损失函数，表示为：

L(y,F)＝ylog(P(y＝1|x))+(1-y)log(1-P(y＝1|x)),y∈{0，1} (5)

计算第m-1棵树的预测函数的当前梯度值：

其中F_m-1(x)是第m-1棵树的预测函数，x_i∈R^12*600是第i个训练集心电样本，

为第i个样本的预测分类概率值，

将回归树作为基学习器，进行线性搜索以得到最优叶子节点值，使用Ne wton-Raphson近似：

其中R_m,j为第m棵树对应的叶子结点区域，第m棵树的预测概率函数如下：

F_m(x)为第m棵树的预测概率函数，预测样本的正分类概率为：

通过上述步骤完成第m棵树的建立，从而完成GBDT模型的构建，并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^GBDT、A₂ ^GBDT、A₃ ^GBDT和对应的针对测试集的预测结果B₁ ^GBDT、B₂ ^GBDT、B₃ ^GBDT；

BiGRU分类模块，其作为第三分类模块，通过对训练数据进行3折交叉的划分，使得BiGRU输出3个第三分类结果，具体的，所BiGRU深度神经网络采用的公式如下：

z_t＝σ(W_z·[h_t-1,x_t]) (10)

r_t＝σ(W_r·[h_t-1,x_t]) (11)

所述公式10-13中：

h_t-1表示上一个神经元的输出；x_t表示本次神经元的输入；W_z表示更新门的权重；σ表示sigmoid函数；z_t表示更新门；r_t表示重置门；w_r表示重置门的权重；

表示当前神经元的候选输出值；w表示输出状态的权重，tanh表示双曲正切函数；表示本次神经元的输出值，采用上述公式1-4分别计算t时刻BiGRU所对应的前向隐层状态h_t，和反向隐层状态h_t'，然后对h_t和h_t'加权求和得到t时刻的隐层状态h_t”，表示为：

h_t”＝w_th_t+v_th_t'+b_j (14)

其中，w_t、v_t分别表示t时刻h_t、h_t'对应的权重，b_j表示待训练的偏置，最后计算心律失常二分类结果，表示为：

其中，

表示为输出的心律失常预测概率，w_i表示为待训练的权重系数矩阵。过上述步骤完成BiGRU模型的构建，并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^BiGRU、A₂ ^BiGRU、A₃ ^BiGRU和对应的针对测试集的预测结果B₁ ^BiGRU、B₂ ^BiGRU、B₃ ^BiGRU。

优选的，步骤S4，包括：

依据模型集成规则—Stacking，将各个分类模块产生的针对训练集上的第一分类结果、第二分类结果及第三分类结果进行行级别上的联结，得到集成数据，表示为：

XI＝[A₁ ^I,A₂ ^I，A₃ ^I]∈R¹，I＝{XGBoost，GBDT，BiLSTM} (15)

然后对三个分类模块预测结果的行级拼接数据进行列级别上的再次拼接，作为三个数据特征，表示为：

X′＝[X^XGBoost，X^GBDT，X^BiLSTM]∈R³ (16)

再把上述得到的三个数据特征，拼接真实的心律失常标签，作为之后的浅层神经网络的输入，表示为：

其中，

为拼接后的样本数据，负责训练之后的融合模型，Y是训练集样本对应的真实标签；

然后搭建好浅层神经网络以对心律失常集成数据进行分类，其中浅层神经网络模块为3层全连接神经网络，输入的数据为

输入3维X数据到3层全连接神经网络中；

第一个全连接层单元的神经元个数为5个，全连接层单元后的激励单元为Relu函数，经过第一层全连接层单元后的特征图维度为5；

第二个全连接层单元的神经元个数为5个，全连接层单元后的激励单元为Relu函数，经过第一层全连接层单元后的特征图维度为5；

第三个全连接层单元的神经元个数为1个，全连接层单元后的激励单元为Sigmoid函数，经过第一层全连接层单元后的特征图维度为1；

所述全连接神经网络采用的公式如下：

其中a^l是第l-1层的输出，a^l-1，b^l-1，W^l-1，f是第l-1层的激活值、偏置、权重和激活函数。

将12导联信号进行上述集成模型的操作后，输出心律失常的预测结果，输出结果的最终表示为：

其中，a³，b³，W³，f是第三层神经网络的输入特征、偏置、权重和激活函数。

优选的，步骤S5具体包括：

初始化所述宽深神经网络的参数，将构建好的数据集随机抽取70％数目的样本当作训练集，其他30％的样本视为测试集；利用三个基分类器，对训练集中的样本继续进行3折交叉验证，并将验证的结果作为融合模块的输入，训练融合模块。每迭代一次更新一次训练参数，直至最后三个基模型及最终的集成模型的损失函数稳定在某一数值附近，停止训练并保存当前网络的训练参数和模型结构信息。

优选的，步骤S6具体包括：将测试集样本输入到已训练好的集成模型中，进行自动识别，获得心律失常的二分类预测结果。

(三)有益效果

本发明提供了一种基于改进多基分类器集成的自动心律失常分析方法。具备以下有益效果：

本发明提出了基于改进多基分类器集成的自动心律失常分类方法，适用于临床上辅助医生进行心律失常的检测。传统的心律失常的检测方法仅仅利用了单一的检测模型，并没有利用Stacking策略集成其他模型一同进行检测，本发明提出的方法基于基于改进多基分类器集成的检测模型，通过对病人心电信号特征进行提取，缓解了心律失常检测难度，使得检测过程中可以利用诸多模型的优势，从而让心律失常检测变得更加准确。

附图说明

图1为基于改进多基分类器集成及诊断规则的自动心律失常分类系统及方法的总体示意图；

图2为Stacking模型集成策略原理示意图；

图3为浅层神经网络示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

请参阅图1-3，本发明提供一种技术方案：如附图1所示，一种基于改进多基分类器集成及诊断规则的自动心律失常分类系统及方法，包括以下步骤：

S1：数据预处理，读入12导联的心电信号的数据，采用小波变换滤除原始心电信号中的工频干扰，采用巴特沃斯带阻滤波器滤除肌电干扰，采用中值滤波滤除原始心电信号中的基线漂移；

S2：心拍数据提取，读取降噪后12导联的心电信号数据，对每个导联心电信号根据同一时刻R波峰值处的位置固定向前截取200个样本点，向后截取400个样本点，使得每个导联的每个心拍固定截取600个点的数据，对相同时刻R波顶点每个导联所截取的600个点的心电信号进行行级上的拼接，此时原始每个心拍的心电信号的数据形成12*600维的样本，作为宽深神经网络模型的输入数据X；

S3：基分类器模型构建，XGBoost分类模块，其作为第一分类模快，通过对训练数据进行3折交叉的划分，使得XGBoost输出3个第一分类结果。其中，对每个心电样本的预测表示为：

其中x_i∈R^12*600是第i个心电训练样本，f_k表示第k棵树，

表示样本x_i的预测结果。在建模过程中不断优化目标函数，使其最小，表示为：

其中

是样本训练误差，Ω(f_i)是树的正则项。XGBoost通过将Obj^t使用二阶泰勒公式展开得：

其中，g_i、h_i表示在第i个样本下泰勒展开式的一阶导数、二阶导数，通过公式8可以得到Obj^t关于w_j的一元二次方程，由此找到使Obj^t取值最小的w_j，表示为：

其中，G_j＝∑g_i；H_j＝∑h_i，通过上述步骤完成第t棵树的建立，从而完成XGBoost模型的构建。并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^XGBoost、A₂ ^XGBoost、A₃ ^XGBoost和对应的针对测试集的预测结果B₁ ^XGBoost、B₂ ^XGBoost、B₃ ^XGBoost；

GBDT分类模块，其作为第二分类模块，通过对训练数据进行3折交叉的划分，使得GBDT输出3个第二分类结果，表示为：

使用负二项对数似然作为损失函数，表示为：

L(y,F)＝ylog(P(y＝1|x))+(1-y)log(1-P(y＝1|x)),y∈{0，1} (5)

计算第m-1棵树的预测函数的当前梯度值：

其中F_m-1(x)是第m-1棵树的预测函数，x_i∈R^12*600是第i个训练集心电样本，

为第i个样本的预测分类概率值，

将回归树作为基学习器，进行线性搜索以得到最优叶子节点值，使用Newton-Raphson近似：

其中R_m,j为第m棵树对应的叶子结点区域，第m棵树的预测概率函数如下：

F_m(x)为第m棵树的预测概率函数，预测样本的正分类概率为：

过上述步骤完成第m棵树的建立，从而完成GBDT模型的构建。并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^GBDT、A₂ ^GBDT、A₃ ^GBDT和对应的针对测试集的预测结果B₁ ^GBDT、B₂ ^GBDT、B₃ ^GBDT；

BiGRU分类模块，其作为第三分类模块，通过对训练数据进行3折交叉的划分，使得BiGRU输出3个第三分类结果，具体的，所BiGRU深度神经网络采用的公式如下：

z_t＝σ(W_z·[h_t-1,x_t]) (10)

r_t＝σ(W_r·[h_t-1,x_t]) (11)

所述公式10-13中：

h_t-1表示上一个神经元的输出；x_t表示本次神经元的输入；W_z表示更新门的权重；σ表示sigmoid函数；z_t表示更新门；r_t表示重置门；w_r表示重置门的权重；

表示当前神经元的候选输出值；w表示输出状态的权重，tanh表示双曲正切函数；表示本次神经元的输出值，采用上述公式1-4分别计算t时刻BiGRU所对应的前向隐层状态h_t，和反向隐层状态h_t'，然后对h_t和h_t'加权求和得到t时刻的隐层状态h_t”，表示为：

h_t”＝w_th_t+v_th_t'+b_j (14)

其中，w_t、v_t分别表示t时刻h_t、h_t'对应的权重，b_j表示待训练的偏置，最后计算心律失常二分类结果，表示为：

其中，

表示为输出的心律失常预测概率，w_i表示为待训练的权重系数矩阵。过上述步骤完成BiGRU模型的构建。并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^BiGRU、A₂ ^BiGRU、A₃ ^BiGRU和对应的针对测试集的预测结果B₁ ^BiGRU、B₂ ^BiGRU、B₃ ^BiGRU；

S4：模型融合，依据模型集成规则—Stacking，如附图2，将各个分类模块产生的针对训练集上的第一分类结果、第二分类结果及第三分类结果进行行级别上的联结，得到集成数据，表示为：

XI＝[A₁ ^I,A₂ ^I，A₃ ^I]∈R¹，I＝{XGBoost，GBDT，BiLSTM} (15)

然后对三个分类模块预测结果的行级拼接数据进行列级别上的再次拼接，作为三个数据特征，表示为：

X′＝[X^XGBoost，X^GBDT，X^BiLSTM]∈R³ (16)

再把上述得到的三个数据特征，拼接真实的心律失常标签，作为之后的浅层神经网络的输入，表示为：

其中，

为拼接后的样本数据，负责训练之后的融合模型，Y是训练集样本对应的真实标签。

然后搭建好浅层神经网络以对心律失常集成数据进行分类，如附图3，其中浅层神经网络模块为3层全连接神经网络，输入的数据为

输入3维X数据到3层全连接神经网络中；

第一个全连接层单元的神经元个数为5个，全连接层单元后的激励单元为Relu函数，经过第一层全连接层单元后的特征图维度为5；

第二个全连接层单元的神经元个数为5个，全连接层单元后的激励单元为Relu函数，经过第一层全连接层单元后的特征图维度为5；

第三个全连接层单元的神经元个数为1个，全连接层单元后的激励单元为Sigmoid函数，经过第一层全连接层单元后的特征图维度为1。

所述全连接神经网络采用的公式如下：

其中a^l是第l-1层的输出，a^l-1，b^l-1，W^l-1，f是第l-1层的激活值、偏置、权重和激活函数。

将12导联信号进行上述集成模型的操作后，输出心律失常的预测结果，输出结果的最终表示为：

其中，a³，b³，W³，f是第三层神经网络的输入特征、偏置、权重和激活函数；

S5：训练模型参数，初始化所述宽深神经网络的参数，将构建好的数据集随机抽取70％数目的样本当作训练集，其他30％的样本视为测试集；利用三个基分类器，对训练集中的样本继续进行3折交叉验证，并将验证的结果作为融合模块的输入，训练融合模块。每迭代一次更新一次训练参数，直至最后三个基模型及最终的集成模型的损失函数稳定在某一数值附近，停止训练并保存当前网络的训练参数和模型结构信息；

S6：模型预测，将测试集样本输入到已训练好的集成模型中，进行自动识别，获得心律失常的二分类预测结果。

通过上述方法步骤，本发明实施例具有以下有益效果：

1、本发明实施例改进了现有的单一模型在心律失常检测中的应用方法，适用于在临床上辅助心律失常的检测。传统的心律失常检测方法仅仅利用了单一的检测模型，并没有融合其他模型的有益信息，本发明提出的方法基于Stacking策略，通过对病人心电信号特征进行提取，充分利用三个基分类器模型的优点，缓解了心律失常检测难度，使得算法更贴近临床检测实际，从而让心律失常检测变得更加准确。

2、本发明实施例使用Stacking策略作为模型的融合策略，相对于传统的心律失常检测算法中仅单一检测，可以多角度的得挖掘心电信号特征，将深度学习模型与机器学习模型相结合，对于因心电信号变化复杂而导致准确度达不到临床要求的检测可以提供更为可靠的实施效果。

需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个......”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

尽管已经示出和描述了本发明的实施例，对于本领域的普通技术人员而言，可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由所附权利要求及其等同物限定。

Claims (3)

1.一种基于改进多基分类器集成的自动心律失常分析方法，其特征在于，包括以下步骤：

S1：数据预处理，读入12导联的心电信号的数据,采用小波变换滤除原始心电信号中的工频干扰，采用巴特沃斯带阻滤波器滤除肌电干扰，采用中值滤波滤除原始心电信号中的基线漂移；

S2：心拍数据提取，通过小波变换检测R波位置，进而对R波前后数据点进行固定长度截取，从而完成对心拍数据提取；

S3：基分类器模型构建，包括用以对心电图数据进行心律失常自动分类的XGBoost分类模块、GBDT分类模块、BiGRU分类模块；其中，XGBoost作为第一基分类模块、GBDT作为第二基分类模块、GRU作为第三基分类模块，通过对训练数据进行K折交叉的划分，使得XGBoost输出K个第一分类结果，GBDT输出K个第二分类结果，BiGRU输出K个第三分类结果；

S4：模型融合，依据模型集成规则—Stacking，将分类模块产生的第一分类结果、第二分类结果及第三分类结果进行集成，得到集成数据，然后搭建好浅层神经网络以对心律失常集成数据进行分类；

S5：训练模型参数，初始化所述集成模型的参数，将构建好的数据集拆分为训练集和测试集；利用训练集样本对集成模型进行训练，生成所述集成模型的参数并保存；

S6:模型预测，对测试集样本进行自动识别，获得心律失常的二分类结果；

XGBoost分类模块，其作为第一分类模快，通过对训练数据进行3折交叉的划分，使得XGBoost输出3个第一分类结果，其中，对每个心电样本的预测表示为：

其中x_i∈R^12*600是第i个心电训练样本，f_k表示第k棵树，

表示样本x_i的预测结果；在建模过程中不断优化目标函数，使其最小，表示为：

其中

是样本训练误差，Ω(f_i)是树的正则项，XGBoost通过将Obj^t使用二阶泰勒公式展开得：

其中，g_i、h_i表示在第i个样本下泰勒展开式的一阶导数、二阶导数，通过公式可以得到Obj^t关于w_j的一元二次方程，由此找到使Obj^t取值最小的w_j，表示为：

其中，G_j＝∑g_i；H_j＝∑h_i，通过上述步骤完成第t棵树的建立，从而完成XGBoost模型的构建，并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^XGBoost、A₂ ^XGBoost、A₃ ^XGBoost和对应的针对测试集的预测结果B₁ ^XGBoost、B₂ ^XGBoost、B₃ ^XGBoost；

GBDT分类模块，其作为第二分类模块，通过对训练数据进行3折交叉的划分，使得GBDT输出3个第二分类结果，表示为：

使用负二项对数似然作为损失函数，表示为：

L(y,F)＝ylog(P(y＝1|x))+(1-y)log(1-P(y＝1|x)),y∈{0，1} (5)

计算第m-1棵树的预测函数的当前梯度值：

其中F_m-1(x)是第m-1棵树的预测函数，x_i∈R^12*600是第i个训练集心电样本，

为第i个样本的预测分类概率值，

将回归树作为基学习器，进行线性搜索以得到最优叶子节点值，使用Newton-Raphson近似：

其中R_m,j为第m棵树对应的叶子结点区域，第m棵树的预测概率函数如下：

F_m(x)为第m棵树的预测概率函数，预测样本的正分类概率为：

通过上述步骤完成第m棵树的建立，从而完成GBDT模型的构建，并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^GBDT、A₂ ^GBDT、A₃ ^GBDT和对应的针对测试集的预测结果B₁ ^GBDT、B₂ ^GBDT、B₃ ^GBDT；

BiGRU分类模块，其作为第三分类模块，通过对训练数据进行3折交叉的划分，使得BiGRU输出3个第三分类结果，具体的，所BiGRU深度神经网络采用的公式如下：

z_t＝σ(W_z·[h_t-1,x_t]) (10)

r_t＝σ(W_r·[h_t-1,x_t]) (11)

所述公式10-13中：

h_t-1表示上一个神经元的输出；x_t表示本次神经元的输入；W_z表示更新门的权重；σ表示sigmoid函数；z_t表示更新门；r_t表示重置门；w_r表示重置门的权重；

表示当前神经元的候选输出值；w表示输出状态的权重，tanh表示双曲正切函数；表示本次神经元的输出值，采用上述公式1-4分别计算t时刻BiGRU所对应的前向隐层状态h_t，和反向隐层状态h_t'，然后对h_t和h_t'加权求和得到t时刻的隐层状态h_t”，表示为：

h_t”＝w_th_t+v_th_t'+b_j (14)

其中，w_t、v_t分别表示t时刻h_t、h_t'对应的权重，b_j表示待训练的偏置，最后计算心律失常二分类结果，表示为：

其中，

表示为输出的心律失常预测概率,w_i表示为待训练的权重系数矩阵；过上述步骤完成BiGRU模型的构建，并通过3折交叉验证的方式，分别产生针对训练集的预测结果A₁ ^BiGRU、A₂ ^BiGRU、A₃ ^BiGRU和对应的针对测试集的预测结果B₁ ^BiGRU、B₂ ^BiGRU、B₃ ^BiGRU；

步骤S4中，依据模型集成规则—Stacking，将各个分类模块产生的针对训练集上的第一分类结果、第二分类结果及第三分类结果进行行级别上的联结，得到集成数据，表示为：

X_I＝[A₁ ^I,A₂ ^I，A₃ ^I]∈R¹，I＝{XGBoost，GBDT，BiLSTM} (15)

然后对三个分类模块预测结果的行级拼接数据进行列级别上的再次拼接，作为三个数据特征，表示为：

X′＝[X^XGBoost，X^GBDT，X^BiLSTM]∈R³ (16)

再把上述得到的三个数据特征，拼接真实的心律失常标签，作为之后的浅层神经网络的输入，表示为：

其中，

为拼接后的样本数据，负责训练之后的融合模型，Y是训练集样本对应的真实标签；

然后搭建好浅层神经网络以对心律失常集成数据进行分类，其中浅层神经网络模块为3层全连接神经网络，输入的数据为

输入3维X数据到3层全连接神经网络中；

第一个全连接层单元的神经元个数为5个，全连接层单元后的激励单元为Relu函数，经过第一层全连接层单元后的特征图维度为5；

第二个全连接层单元的神经元个数为5个，全连接层单元后的激励单元为Relu函数，经过第一层全连接层单元后的特征图维度为5；

第三个全连接层单元的神经元个数为1个，全连接层单元后的激励单元为Sigmoid函数，经过第一层全连接层单元后的特征图维度为1；

所述全连接神经网络采用的公式如下：

其中a^l是第l-1层的输出，a^l-1，b^l-1，W^l-1，f是第l-1层的激活值、偏置、权重和激活函数；

将12导联信号进行上述集成模型的操作后，输出心律失常的预测结果，输出结果的最终表示为：

其中，a³，b³，W³，f是第三层神经网络的输入特征、偏置、权重和激活函数；

步骤S5中，初始化宽深神经网络的参数，将构建好的数据集随机抽取70％数目的样本当作训练集，其他30％的样本视为测试集；利用三个基分类器，对训练集中的样本继续进行3折交叉验证，并将验证的结果作为融合模块的输入，训练融合模块；每迭代一次更新一次训练参数，直至最后三个基模型及最终的集成模型的损失函数稳定在某一数值附近，停止训练并保存当前网络的训练参数和模型结构信息。

2.如权利要求1所述的一种基于改进多基分类器集成的自动心律失常分析方法，其特征在于，步骤S2，具体包括：

读取降噪后12导联的心电信号数据，对每个导联心电信号根据同一时刻R波峰值处的位置固定向前截取200个样本点，向后截取400个样本点，使得每个导联的每个心拍固定截取600个点的数据，对相同时刻R波顶点每个导联所截取的600个点的心电信号进行行级上的拼接，此时原始每个心拍的心电信号的数据形成12*600维的样本，作为基分类模型的输入数据X。

3.如权利要求1所述的一种基于改进多基分类器集成的自动心律失常分析方法，其特征在于，步骤S6，具体包括：将测试集样本输入到已训练好的集成模型中，进行自动识别，获得心律失常的二分类预测结果。

Images