基于变分模态分解的心电信号QRS波检测方法与流程

本发明涉及生物医学信号提取技术领域，具体的说，是一种基于变分模态分解的心电信号qrs波检测方法。

背景技术：

心电信号是心脏电生理活动在人体体表的表现，提供心脏功能等生理状况的有重要价值的临床医学信息，是临床心脏病智能诊断的基础。

qrs波群是心电信号中最重要的部分，在临床医学上qrs波是心脏信息的重要载体，能够反映出很多心脏问题。

qrs波的检测算法包括差分阈值法、双阈值检测算、经验模态分解和小波变换法。差分阈值法在理论上是最简单且速度最快的方法，但由于心电信号存在高频噪声和基线漂移等噪声干扰，使得差分阈值法在实际应用中效果不理想。双阈值检测算法检测效率较高，但是由于该算法没有确定的数学模型，当心电信号出现异常时该算法会出现较多的误检和漏检。小波变换在心电信号检测中具有良好的时频特性，能够对心电信号的特征进行有效的检测。但是在应对实际检测到的心电信号时，由于心电信号存在一些干扰，虽然小波变换相较于其它算法有一定的优势但是也会出现一些误检和漏检得情况。

技术实现要素：

针对现有技术的不足，本发明提供了一种心电信号qrs波的检测方法，能够有效地检测出心电信号的qrs波。

为达到上述目的，本发明提供了一种基于变分模态分解心电信号qrs波的检测方法，具体步骤如下。

步骤（1）：对心电信号进行变分模态分解，选取最优分解层数c。对一个心电信号添加三个参考信号，得到仿真的心电信号。对仿真心电信号依次进行k=5、6、7、8的变分模态分解，对每一层模态分解进行频谱分析，将能够完全分解出三个参考信号的对应k值，认为是变分模态分解的最有分解层。x为心电信号，经过k层变分模态分解得到k个模态分量；每个模态分量的频率中心表示为：，在迭代求解过程中不断更新模态分量及其频率中心，更新步骤如下：

（a）初始化和n为0；

（b）根据下面两式更新和，为惩罚因子（取值2000），λ为拉格朗日乘法算子；

（c）更新λ，其中τ为噪声因子，τ=0；

（d）如果或n>k，则停止迭代，否则换回执行步骤（b），ε表示判断约束条件，ε=0.01。

步骤（2）：将心电信号进行c层变分模态分解，得到c个模态分量。

步骤（3）：确定包含心电信号qrs波的模态分量。对每一层模态分量进行频谱分析，根据qrs波的频率范围选择对应的模态分量进行重构。

步骤（4）：对重构的心电信号求取模极值对。通过小波变换的奇异值检测原理确定心电信号的奇异值进而求取心电信号的模极值对。

步骤（5）：根据信号的模极值对定位r波峰值位置，进而确定qrs波的范围。

正负模极值对之间的零点，就是r波的位置。r波为中心在100ms范围内定位qrs波的起始点。在距离r波的一定范围内定位q波峰，与pr段电位比较，如果大于或等于pr段电位就把q点认为是qrs波的起始点。同样以r波作为中心在100ms范围内定位qrs波末位置。由ecg可知在s波后会出现一个斜率转折点，能找到这个转折点，则这一点认为是心电信号qrs波的终点位置。

附图说明

图1为基于变分模态分解的心电信号qrs波的检测方法的流程框图。

图2为正常心电信号。

图3为仿真心电信号。

图4为心电信号变分模态分解中参考信号分布情况图。

图5为仿真心电信号7层模态分解图。

图6为待分解心电信号。

图7为心电信号7层模态分解。

图8为心电信号小波分解4尺度分解图。

图9为心电信号qrs波检测结果。

具体实施方法

下面结合附图和实例对本发明做进一步说明,用于实施的硬件环境是普通计算机，软件环境是：matlabr12和windows10。我们用matlab软件实现了本发明提出的方法。被分析信号来自美国麻省理工学院与bethisrael医院联合建立的mit-bih心电数据库中的心电信号。

参照图1，图1是本发明方法，基于变分模态分解的心电信号qrs波的检测方法的流程框图，具体实施还需要包括以下步骤。

步骤（1）：对心电信号进行变分模态分解，选取最优分解层数c。对如附图2所示的心电信号添加高频、中频和低频三个参考信号，根据心电信号的频率特点添加的三个参考信号频率分别为10hz、30hz、50hz；参考信号幅值分别为0.2mv、0.05mv、0.01mv；得到仿真的心电信号，如附图3所示。对仿真心电信号依次进行k=5、6、7、8的变分模态分解，对每一层模态进行频谱分析，将能够完全分解出三个参考信号的对应k值，认为是变分模态分解的最有分解层。经频谱分析后k=5、6、7、8的变分模态分解模态中三个参考信号情况如附图4所示当心电信号分解层数k为7时，分解得到的模态分量包含所添加的低、中、高三个参考信号如附图5所示，因此心电信号的最佳分解层数为c=7。

步骤（2）：将如附图6所示的心电信号进行c层变分模态分解，得到c个模态分量。变分模态分解的参数设置如下：c=7，α=2000，τ=0，tol=1e-6。逐步更新μk，ωk和λ，直到满足停止迭代的精度条件，得到7个模态分量。

步骤（3）：确定包含心电信号qrs波的模态分量。对心电信号经过变分模态分解后得到的模态分量做频谱分析，由如附图7所示的模态分量的频谱范围可知心电信号的qrs波特性主要集中在模态3、4、和5上。因此将模态3、模态4和模态5重构为新的待检测心电信号。

步骤（4）：对重构的心电信号求取模极值对。通过小波变换的奇异值检测原理确定心电信号的奇异值进而求取心电信号的模极值对如附图8所示。

步骤（5）：根据信号的模极值对定位r波峰值位置，进而确定qrs波的范围。模极值点连线的过零点位置是心电信号的峰值点，以此确定心电信号r波峰值位置。r波为中心在100ms范围内定位qrs波的起始点。在距离r波的一定范围内定位q波峰，与pr段电位比较，如果大于或等于pr段电位就把q点认为是qrs波的起始点。同样以r波作为中心在100ms范围内定位qrs波末位置。由ecg可知在s波后会出现一个斜率转折点，能找到这个转折点，则这一点认为是心电信号qrs波的终点位置。检测结果如附图9所示。