[发明专利]一种ECG信号联合基线校正及降噪的方法

说明书

一种ECG信号联合基线校正及降噪的方法，使得经过基线校正和降噪后恢复的信号，能够保持原始信号的平滑，并且恢复的信号基本保留了原始信号的细节信息特征。有效地进行了基线校正和噪声抑制，明显改善了传统滤波算法中存在的ECG峰值欠估计的问题，保证了恢复ECG信号的真实性。

技术领域

本发明涉及心电图信号处理技术领域，具体涉及一种ECG信号联合基线校正及降噪的方法。

背景技术

心电图(electrocardiogram,ECG)是有关心脏活动相关的电位变化图，是可以检测健康或病变心脏的电信号。然而，在ECG信号采集过程中，不可避免地会受到人为干扰和噪声的影响，其中包括随机噪声，基线漂移，肌肉伪影及周围电场的干扰。基线漂移主要由记录者呼吸及身体移动所产生的人为干扰，其频率较低，属于一种低频干扰。由于基线漂移的频率与ECG信号ST段会有重叠部分，会抬高ST段，进而影响判断。与此同时，由于心电信号本身比较微弱，易受到周围环境噪声的影响，且噪声的频谱与ECG信号重叠，使得ECG信号难与噪声分离，无法采用传统滤波方法有效消除。因此，ECG信号的基线校正和降噪具有重要意义。

目前传统的小波算法虽然也可以同时达到基线校正和降噪的目的，但是由于传统的滤波会对ECG原始信号进行压缩，因此恢复后的ECG信号会存在极大地峰值欠估计问题，造成ECG原始信号细节的丢失。

发明内容

本发明为了克服以上技术的不足，提供了一种在充分抑制基线漂移和噪声的同时，保留更多ECG原始细节信息，改善传统滤波中存在的峰值欠估问题的ECG 信号联合基线校正及降噪的方法。

本发明克服其技术问题所采用的技术方案是：

一种ECG信号联合基线校正及降噪的方法，包括：

a)计算机加载ECG信号数据y，提取ECG信号数据长度为N；

b)设置LTI滤波器，输入LTI滤波器的截止频率为fc、滤波阶数d、拆分矩阵的阶数K；

c)通过公式P＝A^T*A，Q＝A*A^T+B*B^T计算LTI滤波器所需的稀疏带状矩阵P、Q，其中A为d阶差分矩阵，B为d-1阶差分矩阵，A^T为A的转置矩阵，B^T为B的转置矩阵；

d)进行LTI滤波与稀疏结合优化，引入非对称惩罚函数θ(x；r)，其中式中r为非对称惩罚因子，x为原始ECG信号；

e)引入优化参数b、d，其中d＝Q^TQP^-1y-λP^Tb，式中λ为差分参数， P^-1为矩阵P的逆矩阵，P^T为矩阵P的转置矩阵，Q^T为矩阵Q的转置矩阵；

f)设置差分矩阵D，当K＝0时，D为N阶单位矩阵，当K≥1 时，N阶单位矩阵中每一行位置对应的数值减去下一行位置相对应的数值得到 (N-K)*N阶矩阵，当K＝1时，其为(N-1)*N阶差分矩阵；

g)设置生成的N阶矩阵为主对角线全为1，其余均为0 的矩阵，引入第三方矩阵M，其中M＝2*Gamma+D^-1*Gamma*D；

h)设定算法迭代最大次数Nit，对最终优化x进行迭代循环计算，其中 x＝P(Q^T*Q+P^TMP)d；

i)当x前后两次的差值小于0.005或达到最大迭代循环次数Nit时停止循环，输出经过基线校正和降噪后的恢复信号x。