小波阈值去噪技术在ECG信号处理中的应用.PDF

小波阈值去噪技术在ECG信号处理中的应用 1、 引言 心脏细胞除极和复极的电生理现象，是心脏运行的基础。心电信号记录了心脏细 胞的除极和复极过程，在一定程度上客观反映了心脏个部位的生理状况，因而在临床医学中 有重要意义。[1] 由于人体心电信号随着检测状态及时间的变化，一方面具有较明显的非平 稳性特点，另一方面还包含了许多干扰，如工频干扰、肌电干扰、呼吸干扰等，以加性或乘 性方式与心电信号混杂，引起心电信号的畸变，从而掩盖了原始心电波形中的特征信息，使 整个心电波形模糊不清，难以进行识别诊断。传统去除干扰的方法有 RC 滤波、数字滤波补 偿基线漂移以及基线拟和等方法，但都存在一些不足之处。本文借助近几年来一种新的信号 时频分析理论——小波变换理论（WT：Wavelet Transforms）运用于心电信号的测量，利用 小波变换多尺度多分辨的特点，将心电信号进行分解，不同频带的信号便显现在小波分解的 不同尺度上，进行信号重构时，去除高频干扰及基线漂移所在尺度的信息，使重构后信号不 再含有干扰成分，以便正确估计心电信号的各特征参数并检出所期望的心电波形，进而提取 有诊断价值的信息。 2、 小波阈值去噪原理 2.1 小波变换 傅立叶分析是将信号分解成一系列不同频率的正弦波的叠加，同样小波分析是将 信号分解成一系列小波函数的叠加，而这些小波函数都是一个母小波函数经过平移与尺度伸 缩得来的。小波变换的定义是把某一被称为基本小波（也叫母小波）的函数做位移 b 后， 再在不同尺度 a 下与待分析的信号 x(t)做内积：其中设 x(t)是平方可积函数（记作 ），ϕ （t）是基本小波或母小波（MW）函数，且满足容许条件 则 称为 x(t)的小波变换。式中： Wx(a, b) 是 x(t)的小波变换式， a 0 是尺度 因子；b 反映位移，其值可正可负，上标*代表共轭 ， 是基 本小波的位移与尺度伸缩。式（1）不但是连续变量，而且 a 和 b 也是连续变换，因此称为 连续小波变换（CWT）。式（1）的等效频域表示为： 式中 X (ω ), Ψ* (aω ) 分别是 x (t )和Ψ (t ) 的傅立叶变换。 从式（1）（2）可知，如果 x (t )为信号函数，则小波变换是信号与小波函数的 内积，是对信号满足一定附加条件的滤波，这种附加条件反映在小波函数及小波因子选择上。 小波变换提供了很好的局部化特性，它既可在时域，也可在频域局部化定位观测。高频时使 用小尺度 a值，低频时使用大尺度 a 值，分析频率有高有低，但各分析频段内分析的品质因 数却保持一致[3] ，如果希望在时域上观测的愈细致，就愈要压缩观察范围，并提高分析频 率。利用小波变换所具有的这种数学显微镜特点和频域带通特性，可以把所需的信号分离出 来，进行分析研究。 2.2 小波阈值消噪算法分析 所谓阈值降噪，就是按照一定的预设阈值压缩信号的小波变化系数，然后用被压 缩后的系数重构以达到降噪的目的。目前应用最广泛的是 Donoho 提出的硬阈值和软阈值降 h 噪方法。因为在小波域中，信号的能量相对集中在某几个位置上，而噪声的分布一般比较广， t 根据瞬时性的特点，信号表现为一些大的系数，而一些小的系数则更多的是由噪声和信号能 t 量的突变所产生的，所以小波阈值去噪主要是利用了有效信号和噪声信号在小波变换下奇异 p 性截然不同的表现特征来去除噪声，保留有效信号。心电信号的主要频率成分在 100Hz 以 : 下，而肌电干扰噪声在 5~2000Hz，所以相对于 ECG 信号来说，肌电信号是一种高频干扰。 / 所以先通过小波分析多分辨率分析方法