席朝岭;邱联奎
途语导航
【摘 要】为去除心电信号中的各种噪声,本文以小波变换的多分辨率分析为理论基础,利用自适应阈值调整小波变换系数,用调整后的系数进行心电信号重建.采用MIT-BIH数据库中的心电信号进行仿真、验证,有效地去除了噪声信号.与传统滤波器具去噪相比有明显的优越性.【期刊名称】王洪文传《中国医疗设备》
胶带机【年(卷),期】2010(025)011
PFOS
【总页数】4页(P19-21,24)
【关键词】ECG;小波变换;心电信号仿真;去噪
【作 者】席朝岭;邱联奎
【作者单位】河南科技大学,电信学院,河南,洛阳,471003;河南科技大学,电信学院,河南,洛阳,471003
【正文语种】中 文
【中图分类】R319
心电信号(ECG)记录了心脏的电生理现象，客观地反映了心脏各个部位的生理状况，为临床诊断提供了心脏功能及疾病的重要信息。由于心电信号比较微弱，幅值范围在50µV ～5mV之间，频率主要在0.05 ～100Hz范围之间，90%的心电信号频谱在0.25～35Hz之间[1]，在信号采集时，易受仪器、病人及环境等因素的影响，会产生基线漂移。受人体肌电、50Hz工频和其它噪声的干扰，心电信号往往会被淹没在较强的背景噪声之中[2]。传统的数字滤波方法大都针对某种噪声单独设计滤波器，对心电信号中的其它干扰不起作用。由于心电信号微弱、信噪比低、奇异点较多以及有用信号频带与噪声频带相互重叠等，使得传统的滤波去噪方法在处理心电这种短时低能量的瞬变信号时，信噪比的提高往往是以信号局部特征的损失为代价的。如何处理保留心电信号局部特征和有效抑制噪声之间存在的矛盾，是心电信号去噪的关键。近年来，随着小波理论研究和应用的不断深入，小波变换由于具有良好的时频局部化特性，它被逐渐运用到心电信号处理之中，成为缓解、处理上述问题的一种新手段。
小波的多分辨率分析特性能将信号在不同尺度下进行多分辨率分解，由于有用信号和噪声信号的频率存在差异，所以在不同分解尺度上的小波系数呈现不同特征。这样小波去噪就是利用多分辨率时频局部化特性，把信号中的不同频率成分分解到不同的尺度空间上，然后对各个尺度下的小波系数利用阈值进行处理，去掉由噪声引起的小波系数，保留由有用信号引起的小波系数，用处理过的小波系数进行反变换，得到去噪后的信号[3]。
1.1 小波基的特性
小波基的数学特性主要有：① 正交性。正交性是小波基的一个非常优良的性质，能够保证信号的精确重构。②正则性。正则性又称规则性,它用于描述信号在某点或某一区间内的光滑性和奇异性，主要影响着小波系数重构的稳定性。③ 消失矩。由数学分析的知识可知，一般光滑函数都能用多项式来描述(如Taylor展开)，因此小波的消失矩越高，光滑函数在小波展开式中的零元就越多(实际小波变换中，严格为零的小波系数也很少) ，消失距越高计算得到的零元越多，对数据压缩非常有用。④ 紧支性。紧支性是小波的重要性质，支集越短的小波，局部化能力越强，越有利于确定信号的突变点，而且支集越短的小波，产生的大幅值的小波系数就越少，这与实际中希望小波系数的幅值尽可能小相一致。此外，短支
集能提高计算速度，这对心电信号的实时处理也是十分重要的。⑤ 对称性。对称性是描述小波性能的一个基本特征。小波具有对称性意味着进行Mallat分解时的滤波器组具有线形相位,可以避免信号分解和重构时的相位失真，是刻画小波性能的一个基本特征[4]。
1.2 心电信号去噪中小波基的选择
小波变换的小波基不具备唯一性, 不同的小波基具有不同的时频特性，选取不同的小波基对同一心电信号进行处理，将得到不同的结果，所以，在小波变换去噪中，小波基的选取是影响小波去噪效果的一个主要因素。小波基的选择也是小波变换去噪中的一个难点，一般情况下，小波基的选择难以统一成一般原则。实际应用中，根据上文所述小波的数学特性，根据要处理信号的特点，针对具体问题进行具体分析，通过经验或不断的试验来选择小波基，用小波变换方法处理信号的结果与理论结果来判定小波基的好坏，由此确定最优的小波基。下面简单介绍几种常见的小波基：dx4
(1) Haar小波。它的优点是它是唯一一个具有对称性的紧支正交实数小波，支撑长度为1，做小波变换时，计算量很小。它的缺点就是光滑性太差，用它重构的信号，会出现“锯齿”现象。Haar小波的这种不连续性使其在实际的信号处理中失去了使用价值，因此不适于
处理心电信号。
(2) Marr小波和Morlet小波。它们的优点是具有清晰的函数表达式，且具有对称性，但是它们的尺度函数不存在，不具有正交性，这就不能够对分解后的信号进行重构，所以也不适于心电信号去噪。
(3) Meyer小波。它是紧支撑正交对称小波，无穷次连续可微，有无穷阶消失矩。但它没有快速算法，影响计算速度，对心电信号来说，实时处理、分析的数据量较大，所以从处理速度方面考虑，一般不采用Meyer小波。
(4) Biorthogonal小波。它是具有对称性的紧支双正交小波，可以用来实现心电信号去噪，但因为该小波系中的各小波基不具有正交性，只具有双正交性，相比具有同样消失矩阶数的正交小波来说，计算的简便性和计算时间可能会受到影响，所以应用时要合理选择滤波器的长度。
JY改造系统
(5) Daubechies小波。它是一类紧支正交小波，通常表示为dbN的形式，其中N对应了小波函数的消失矩的阶数，且支撑长度为2N-1，正则性随着N的增加而增加，这些对于心电信
号去噪来说都是很好的特性。但该类小波对称性很差，导致信号在分解与重构时相位失真严重，所以该类小波不太适合心电信号去噪。
(6) Smylet小波。它是近似对称的一类紧支正交小波函数，它具有Daubechies小波系的一切良好特性，在对称性方面进行的改进,使该小波系在处理心电信号时，可以很大程度避免不必要的失真，因此适合对心电信号去噪。
(7) coiflet小波。它也是一类具有近似对称性的紧支正交小波。消失距为N时，支撑长度为6N-1，理论上它可以用于心电信号去噪，而且coifliet小波symN小波的对称性要好一些，但这是以支撑长度的大幅度增加为代价的, coiflet具有与smy3N相同的支撑长度。对于心电信号来说，由于数据量较大，我们希望所选取小波的支撑范围尽量短，这样有利于节省处理时间，所以从计算时间上考虑，该类小波不如snytelt小波系。
综上所述，在这几类常用的小波中，只有Biohrtgoonal小波系和symlet小波系比较适合于心电信号去噪，且各具优缺点：Biohrtogonal小波系严格对称但不正交，symlet小波系正交但不严格对称。将信号作小波变换后，由于尺度系数反映的是信号的概貌信息，所以尺度函数的波形越接近于信号的形状，越有利于在不同的分解尺度下保留原始信号的形态。通过
比较，smyelt小波系中的各个小波基所对应的尺度函数更接近于心电信号中完整的P-QRS-T波，所以本文选择了symlet小波对心电信号进行去噪处理[5]。
1.3 小波分解层数的确定
小波分解层数(或称分解尺度)的确定非常重要，分解层数不同，去噪效果也不同，有时甚至差别很大。利用小波变换去噪时，通常利用信号和噪声的频率范围确定分解层数。根据心电信号及所含噪声的特点，本文使用symlet小波进行8尺度分解。
1.4 阈值函数的选取与阈值λ的确定
小波去噪的基本思想是：噪声与信号在各尺度(即各频带)上的小波系数具有不同的特点，这些特点是在小波域中，信号的能量相对集中在某几个位置上，而噪声的分布一般比较广。根据瞬时性的特点，信号表现为一些大的系数，而一些小的系数则更多的是由噪声和信号能量的突变所产生的。将各尺度上由噪声产生的小波分量去掉，这样保留下来的小波系数基本上就是原始信号的小波系数，再利用小波逆变换重构出原始信号。
由此可知：小波去噪的关键是如何选择阈值除去由噪声产生的小波分解系数分量。阈值的
选择同样非常关键，如果阈值过小，则方差偏大，数据欠平滑，去噪后仍留有噪声；但阈值过大，会使数据过于平滑，信号的奇异性可能丧失，重要的信号特征会被滤掉，引起偏差[6]。目前应用最广泛的是Donoho提出的硬阈值和软阈值降噪方法。根据心电信号中噪声的特点，阈值的选取不能单一，硬阈值由于阈值函数不连续性会引起较大的方差，且不稳定，对数据小的变化比较灵敏。软闭值法能提供较好的光滑性，具有较好的去噪效果，但有时会损失过多的有用信号信息，不易于保持信号的特征。所以本文采用软阈值σ方法处理，选用自适应阈值来克服这种缺点。阈值选取公式如下：λ=σ/ln(j+z)[7]，其中，N为心电信号采样点数，j为本级尺度，z为常数，本文实验中取z=1。
1.5 小波变换去噪的步骤
小波变换去噪法是一种原理简单的去噪方法，它可以简单的分为以下几个基本步骤：
(1)根据需要选择小波基，确定分解层数，然后将含噪声的信号在不同尺度下进行小波分解，分解到不同频带内，得到含噪信号的小波系数。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议