使用LabVIEW进行心电信号处理.docVIP

使用LabVIEW进行心电信号处理 心电图是一种记录心脏产生的生物电流的技术。临床医生可以利用心电图对患者的心脏状况进行评估，并做出进一步诊断。ECG记录是通过对若干电极(导联)感知到的生物电流进行采样获得的。图1中显示了典型的单周期心电图波形。 图1 典型的单周期心电图波形 通常说来，记录的心电信号会被噪声和人为引入的伪影所污染，这些噪声和伪影在我们感兴趣的频段内，并且与心电信号本身有着相似的特性。为了从带有噪声的心电信号中提取出有用的信息，我们需要对原始的心电信号进行处理。 从功能上来说，心电信号的处理可以大致分为两个阶段：预处理和特征提取(如图2所示)。预处理阶段消除和减少原始心电信号中的噪声，而特征提取阶段则从心电信号中提取诊断信息。 图2 典型的心电信号处理流程图 使用LabVIEW和相关工具箱，如高级信号处理工具箱(ASPT)和数字滤波器设计工具箱(DFDT)等，用户可以方便地创建针对两个阶段的信号处理应用，包括消除基线漂移、清除噪声、QRS综合波检测、胎儿心率检测等。本文着重讨论使用LabVIEW进行典型的心电信号处理的方法。 1. 心电信号预处理 心电信号预处理可以帮助用户去除心电信号中的污染。广义上讲，心电信号污染可以分为如下几类： ? 电源线干扰 ? 电极分离或接触噪声 ? 病人电极移动过程中人为引入的伪影 ? 肌电（EMG）噪声 ? 基准漂移 在这些噪声中，电源线干扰和基准漂移是最为重要的，可以强烈地影响心电信号分析。除了这两种噪声，其它噪声由于可能是宽频带的且复杂的随机过程，也会使心电信号失真。电源线干扰是以60 Hz (或 50 Hz)为中心的窄带噪声，带宽小于1Hz。通常，心电信号的采集硬件可以消除电源线干扰。但是，基准漂移和其它宽带噪声通过硬件设备很难抑制。而软件设计则成为更为强大而可行的离线式心电信号处理方法。用户可以使用以下方法来消除基准漂移和其它宽带噪声。 消除基准漂移 基准漂移的产生通常源于呼吸，频率在0.15 到 0.3 Hz之间，可以通过使用高通数字滤波器进行抑制。用户还可以使用小波变换通过消除心电信号的趋势来消除基准漂移。 1. 数字滤波器方法 LabVIEW 的DFDT工具箱提供了图形化和交互式的方法，用于快速而高效地设计和实现有限冲击响应(FIR)或无限冲击响应(IIR)滤波器。例如，用户可以使用Classical Filter Design Express VI设计Kaiser窗FIR高通滤波器消除基准漂移。图3显示了一个高通滤波器的实例，用户可以使用这个实例VI的程序框图来消除基准漂移。 图3 设计并使用高通滤波器消除基准漂移 2．小波变换方法 除了数字滤波器，小波变换也是一种消除指定频带内信号的有效方法。LabVIEW ASPT工具箱提供了WA Detrend VI，它可以消除信号的低频趋势。图4显示了使用WA Detrend VI消除基准漂移的实例。 图4?使用WA Detrend VI消除基准漂移 这个实例使用了Daubechies6 (db06)小波，因为这种小波与实际的心电信号相似。在这个实例中，心电信号的采样周期为60秒，共1200个采样点。这样根据如下公式，趋势级别为0.5： ???????????????????????????????????????????????????????????????? 其中，t是采样周期而N为采样点数。 图5显示了原始的心电信号，以及由基于数字滤波器方法和基于小波变换方法处理后得到的心电信号。你可以看到处理后的心电信号在保持原有心电信号主要特性的同时，几乎不含有基准漂移信息。你还可以注意到基于小波变换的方法更具优势，因为这种方法不引入延时，而且比基于数字滤波器方法的失真更小。 消除宽带噪声 在消除了基准漂移后，得到的心电信号比原来的信号更加清晰和稳定。但是，其它类型的噪声仍然会影响心电信号的特征提取。这些噪声是宽频带的复杂随机过程，所以不能使用传统的数字滤波器。为了消除宽带噪声，需要使用Wavelet Denoise Express VI。 这种Express VI首先应用小波变换将心电信号分解到各个子带，然后利用阈值或收缩功能调整小波系数，最后重建出消除噪声后的信号。下图显示了对心电信号使用非抽样小波变换(UWT)的实例。 图6 利用UWT消除心电信号中的宽带噪声 UWT相比于离散小波变换(DWT)，在平滑度和精度之间有着更好的折中。通过将去噪后的心电信号和未经去噪的心电信号相比较，如图7所示，可以发现宽带噪声被极大地抑制而心电信号的细节则保持不变。 图