基于小波变换的一种QRS波检测方法研究.pdfVIP

基于小波变换的一种QRS波检测方法宰 李芳胡志忠 南京航空航天大学信息科学与技术学院，江苏南京210016 摘要：本文将小波变换应用于心电信号(EcG)QRS波检测，选取双正交样条小波函数对心电信号 按Mallat算法进行快速变换；从等效滤波器的角度分析了信号奇异点(R波峰值点)与其小波变换模极 大值对的零交叉点的关系；在检测中还运用了一系列策略以增强算法的抗干扰能力。经MIT／BIH标准 心律失常数据库验证，这种方法具有很好的效果。 关键词：双正交样条小波心电信号QRs波Mallat算法 l引 言 心电图(ECG)在临床诊断方面具有蘑要的意义，尤其是对心血管疾病的诊断。QRS波是心电信号中最 具有特征的波形，波形的榆测都是由QRS波检测开始的，心率的计算和分析以及P波、T波的检测，都是 测算法已发展多年，但以往人们丰要采用滤波器法和模板I』‘配法研究心电信号ll之J。但是，’前者抗干扰能力 差而后者对抗高频和基线漂移不明显，且计算时间较长。因此，它们的局限性很大。 小波变换是一种窗口大小(即窗口面积)固定但其形状町改变的，时间窗和频率窗都可改变的时频局 部化分析方法，正是这种特性，使得小波变换被广泛应用于工程技术领域中。其中一个重要应用就是检测 信号的特征点。近年来，随着小波分析技术的成熟，国内外不少学者把小波变换用于心电信号的波形检测 和特征提取的研究中。西安交大李翠微、郑崇勋等13】和印度Sahambi等人以二次样条小波和高斯函数的一 次导数作为小波，提高了QRs波的检测率。为了进一步优化和完善现有的各种QRS波榆测算法，本文采 用双正交样条小波变换对QRs波进行检测，实验结果表明，该方法在朱经任何预处理的条件F，能较为准 确的检测QRS波群特征。 2小波变换 2．1小波变换和Mauat算法 由小波分析理论知道，非平稳信号厂(x)的连续小波变换定义为： 彬厂(x)=厂(x)·虮(z)=三e／(，渺(等矽 其中s称为尺度；虬(x)=去沙(兰)是母小波y(工)在尺度s上的伸缩。 对于数字信号／(栉)进行二进小波变换，应用Mallat算法求得，其数字滤波器的表示形式为 s：’厂(刀)=∑吃是川／∽一2川后) 士，主要研究方向电路与系统，信号检测及处理． W2’厂(刀)=∑g。蹦。1／o一2川七) (3) ^Eo 其中，是o／(刀)是要处理的EcG信号，嘎’厂(甩)是信号／(胛)的二进制小波变换，是。称为光滑算子，魄 和g。分别是正交的数字低通滤波器H(国)和高通滤波器G(国)的系数。 2．2小波基的选取 ‘由式(2)和(3)可见，小波基函数的选取不是惟一的，选用不同的小波基对信号进行小波变换，可以突出 不同特点的信号特征的具体形式，由参考文献【41知道样条小波同时具有线性相位、短支集、对称和反对称 性、及半正交性等优良性质。它的滤波器系数个数有限，且系数都是有理数。本文选用二次B样条小波， 它是一平滑函数的一阶导数，使用文献‘51中的滤波器得到小波分解滤波器系数为：ho=l／4，hl=3／4， 92=3／4， h2=3／4，h3=l／4；go=一1／4，gI=一3／4， 93=l／4。 3QRS波检测的原理 心电信号是一种微弱的牛物电信号，是有一系列波组成的，包括R波、QRS复合波、P波、T波和u 波组合而成的。QRS波代表左、右两心室除极过程的电位变化：QRS波的第一个向下的波为Q波；继Q 波之后的向。上的波为R波；R波后的向下波为S波。 小波变换w|厂(x)在S充分接近于零时，其模的最大值点是信号变化的最激烈点。对于二进制的小波变换， 信号的奇异点在较小的尺度上对应小波变换的模极大值点。由于S=2，(j∈Z，且．．0)不能真正趋于零，对于 数字信号的二进小波变换，信号的奇异点对应小波变换在s=21尺度上正极大值、负极大值对的零交叉点。 而在小波变换的