本发明专利技术公开了一种基于提升小波变换和改进近似包络的心电信号特征检测算法，属于微弱生物电信号处理技术领域。目前临床应用的心电信号检测技术多数不能兼顾检测精度与实时性要求，本发明专利技术提出了一种基于提升小波改进半软阈值去噪和改进近似包络的心电信号进行预处理算法以及基于斜率阈值的心电特征检测算法，根据心电信号的波形特点及时域分布特性设定检测准则，对心电信号分别进行R波、QRS波群起止点、P波及T波位置的检测。本发明专利技术检测算法简单，速度快，适合并行处理，占用内存空间少，便于DSP芯片实现。该算法即使在强噪声和大P/T波干扰的情况下也能准确检测出R点位置，与MIT-BIH数据库标注比较，对包含严重噪声干扰的105号数据的R波误检率仅为0.27％。
【技术实现步骤摘要】
本专利技术属于微弱生物电信号处理
，特别是受严重噪声干扰的微弱心电信号的去噪及其特征点的实时检测 算法。
技术介绍
由于心电图(electrocardiogram, ECG)诊断心脏病和心血管疾病的方法具有无创性，方便易行使心电图诊断在临床上的到了广泛的应用。无创的体表电极采集到的心电信号比较微弱，仅为毫伏级，所以ECG信号极易受到外界环境干扰。这些干扰包括工频干扰、基线漂移、肌电干扰、电极接触噪声、电极极化噪声、运动干扰和放大电路内部噪声。为了抑制各种干扰，增强心电信号中的有效成分，提高心电特征点检测的准确率，除了对心电记录仪的硬件抗干扰能力有较高的要求外，对采集到的ECG信号经模数转换后的数字滤波处理也是非常重要的。目前对ECG信号的检测方法有很多，大致可以分为三大类，即经典的数字滤波器技术、自适应滤波器技术和以小波变换、数学形态学、神经网络及基于经验模式分解(Empirical Mode Decomposition, EMD)的方法等为代表的现代信号处理技术。这些算法在有其处理优越性的同时也存在着不足，如滤波法，设计思想简单，处理速度快，大多可实现实时检测，但检测精度不够高。神经网络法虽然检测速度比较快，精度也比较高，但需要依靠人工网络的学习训练，不仅对训练样本的广泛性和代表性有较高的要求，而且训练学习要花费较多的时间，实际应用困难。数学形态学算法对信号前期的预处理要求高，并且计算量较大。EMD算法与小波变换法类似，虽然能够达到较好的检测精度，但算法的计算量非常大，且目前还没有相应的快速算法。小波变换法以其优越的时频局域特性在心电信号处理方面得到广泛的研究与应用，基于小波变换的检测算法虽然能够达到很高的检测率，但实时性差，即使使用Mallet快速分解算法，计算量仍比较大，实现实时检测比较困难。1996年Sweldens成功建立了一种新的实现小波变换的结构-提升框架，这种提升格式构造出来的小波称之为第二代小波。它既保留了第一代小波的时频局域特性，同时也克服了其平移和伸缩的不变性。此外，这种基于提升格式的小波变换的所有运算均在时域进行，无需Fourier分析，运算速度相比第一代小波变换大大提高，适合并行处理，占用内存空间少，便于DSP芯片实现。本专利技术采用提升小波改进半软阈值法对心电信号进行去噪处理，为了提高R波检测准确率，本专利技术通过一种改进的近似包络法在增强R波能量的同时能够很好的抑制噪声和大P/T波的干扰，该算法在确保准确率的同时运算速度快，占用内存少，可用于临床上心电图仪采集的心电信号处理与检测。
技术实现思路
本专利技术的目的在于解决目前存在的心电特征检测算法不能同时兼顾检测精度与实时性的问题，提供一种检测精度高，实时性强的基于提升小波改进半软阈值与改进近似包络的心电信号预处理算法和基于斜率阈值的检测算法。本专利技术通过以下技术方案实现基于提升小波改进半软阈值和改进包络预处理算法以及基于斜率阈值心电信号特征检测算法的实现步骤如下(I)从MIT-BIH数据库中提取心电信号x(t)，确定小波基函数以及提升小波分解层数N，对信号X (t)进行N层分解；(2)设定提升小波高频系数的阈值Thl，对提取的高频提升小波系数进行改进半软阈值去噪处理；(3)确定重构心电信号的尺度信号，对心电信号进行重构得重构信号y(t)；(4)对重构心电信号y(t)进行Hilbert变换并取改进的近似包络，完成心电信号的预处理得信号z (t)；(5)设定R波检测的斜率阈值Th2，依照斜率阈值检测策略，对信号z⑴进行R波位置的检测；(6)采用误检策略以及漏检回溯策略进行误检与漏检的查询；(7)根据临床心电信号波形特征点的时域分布特性，以检测到的R波位置为基准，检测QRS波群起止点、P波及T波的位置。在上述步骤(I)中，根据小波基函数特性以及心电图信号的时频分布确定小波基函数和分解层数N。在上述步骤⑵中，设定提升小波系数系数阈值m = Ioge(IengtKdj)),对提取的提升小波高频系数进行改进半软阈值去噪处理。在上述步骤(3)中，重构心电信号时，舍掉高尺度与低尺度信号，选择中间尺度对心电信号进行重构。在上述步骤(4)中，本专利技术提出一种改进的近似包络，能够在增强R波能量的同时很好的抑制各种噪声以及大P/T的干扰，比传统包络更容易实现R波的检测。在上述步骤(5)中，由于预处理后心电信号中R波斜率明显增大，所以本专利技术提出基于斜率阈值的特征检测算法。在上述步骤(7)中，在对其它特征点进行检测时，由于取改进近似包络后其它特征波形受到抑制，所以在检测其它特征波时，以检测到的R波位置为基准，在重构的信号进行其它特征点的检测。本专利技术的心电信号特征点检测算法的效果和益处是本专利技术基于提升小波变换和改进近似包络，提高了心电信号的去噪处理的速度。经过改进近似包络处理，提高了 R波的检测精度，也使得其它特征点的检测精度有所提高。该处理算法适合并行处理，速度快，精度高，占用内存空间少，便于DSP芯片实现。附图说明图I是本专利技术算法总体流程图；图2是心电信号预处理及R波检测流程图3是向心电信号中加入噪声后，用本专利技术的预处理算法去噪的效果图；图4是特征点检测结果图；图5是本专利技术算法成功判定室性期前收缩位置图；图6是本专利技术算法成功去除伪迹干扰图；图7是本专利技术算法成功去除强噪声干扰图；图8是本专利技术算法在存在基线漂移的情况下成功检测为弱R波图；图9是本专利技术算法在大T波干扰下成功检测R波位置图。 具体实施例方式以下结合技术方案和附图详细叙述本专利技术的具体实施例。本专利技术基于提升小波变换、改进近似包络和斜率阈值的心电图信号中R波、QRS波群起止点、P波以及T波位置的检测算法，图2是心电信号预处理及R波检测的具体流程图，具体实施步骤为I.提升小波改进半软阈值去噪根据心电信号特征波形以及噪声的频域分布特性，利用提升小波小波变换的时频局域特性，选用sym8小波作为小波基函数对心电信号进行三层提升分解，提取高频小波系数并进行改进半软阈值处理，选择三尺度信号对心电信号进行重构。2.改进近似包络处理，首先对重构的心电信号y(t)依据下式进行Hilbert变换，H\y{t)\ = yh ( ) = y(t)*h(t) = —(I) Tt由传统包络定义可得包络信号E1 (t)为El(t) = ^y{tf+yh(t)2⑵为了增强QRS能量并减弱噪声和不需要的P/T波。本专利技术提出一种改进的近似包络，定义如下E2 (t) = y(t) |3+|yh(t) I4⑶图3是向提取的心电信号中加入噪声后，用以上步骤预处理算法处理的效果图。3. R波波峰位置检测策略由于心电信号随机性强，同一被检者的R波的幅值和间隔可能在短时间内发生剧烈的变化，所以本专利技术提出了依据斜率为阈值的特征检测方法，并设定平均R-R间隔值根据当前检测R波位置自适应变化。此外本专利技术还采用自适应阈值、回溯检漏和跳检的检测策略，更好的防止了漏检、误检的发生，R波的具体检测步骤如下(I)设定预处理后信号的起始点i，求该点斜率，并取绝对值dif (i)；(2)判断dif (i)与设定阈值(本文在检测R波时取阈值Th2 = 0. I)的大小，若dif(i) > Thr且其后连续四点的斜率值递增并满足阈值条件，则在i点前后各本文档来自技高网...
【技术保护点】
【技术特征摘要】
【专利技术属性】
技术研发人员：李鸿强，王小飞，
申请(专利权)人：天津工业大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。