基于matlab的滤波器设计第三章.doc

第三章 ECG检测方法 本章将主要介绍本文所采用的ECG信号的检测方法的基本原理。首先要对ECG信号进行预处理：采用滤波技术消除各种噪声，采用Tompkins算法检测出QRS波的峰值，并对周期性的ECG信号进行分割。然后介绍采用AR建模技术提取ECG信号的特征的基本原理，以及如何采用相关系数及信噪比两个指标来确定建模阶次P，最后介绍本文采用的BP神经网络分类法和基于二次判别函数分类方法的基本原理及相关的计算公式。 3．1信号的预处理 本文中所使用的ECG信号取自MIT-BIH数据库：NSR信号取自“MIT-BIH arrhythmia database”，其采样频率是360Hz；VT信号和VF信号取自“MIT-BIH Ventricular Arrhythmia database”，其采样频率是250Hz。信号的频率都统一转化为360Hz。在提取信号特征之前，需要对ECG信号进行预处理：采用滤波技术消除实时测量时可能存在的噪音，采用Tompkins算法检测出QRS波的峰值，并对周期性的ECG信号进行分割，选取有效的窗口数据。 1.信号的滤波 在实时测量得到的ECG信号通常带有噪音，噪声的来源是多种多样的，病人的呼吸 ，电极的移动，电源的工频干扰，肌肉收缩引起的高频噪音等，这些噪音和干扰会对ECG信号检测准确性带来很大的影响。所以首先要消除这些噪音，才能进行进一步的分析。据有关资料显示，电源引起的噪音频率约(50-60Hz)，呼吸引起的噪音频率约0.2Hz左右，电极移动引起的低频噪音约0.3HZ肌肉收缩引起的高频噪音。本文采用带通滤波器（BPF）（由一个低通滤波器和一个高通滤波器）组成 ，对ECG信号进行滤波。其上下边带截止频率是1Hz和50Hz，能有效地消除各种噪音[24]。 低通滤波器（LFP）的系统函数如下： ( 3-1) 2. QRS波峰值检测 本文中的ECG信号QRS波的检测算法最初由Pan Tompkins提出，Hamilton 和Tompkins对该算法进行了进一步的研究，QRS波检测包括QRS波位置、宽度、面积的检测ECG 图3-1 QRS波检测系统框图 ECG信号通过一个带通滤波器（由一个低通滤波器和一个高通滤波器串联组成）。再进行微（差）分、平方、滑动窗口平均、最后通过自适应阈值的检波方法，检测到QRS波的峰值。采样的心电信号经过数字带通滤波器；微分后再平方，平方过程可突出微分后的频响曲线的斜率，并有助于减少由于T波具有更多高频分量而引起的；再做移动窗口(150ms)积分，积分所得的结果既包含QRS波群的斜率，又包含QRS波群的宽度的信息，然后应用两种阈值进行检测，一是滤波后信号的阈值，另一是由窗口积分后的阈值。阈值的大小随信号的波动不断调整，这样可提高检测的可靠性，特别可以防止噪声引起的。 Hz，所以我们用一带通滤波器（用一个上限频率为15 Hz的低通滤波器与一个下限频率为5 Hz的高通滤波器串联构成）将QRS波与T波、U波等分离开来。 低通滤波器（LFP）的系统函数如下： （3-2） 低通滤波器的输入、输出关系还可以用差分方程表示为： （3-3） 式中 x(nT) :是低通滤波器的输入信号， y(nT) 是低通滤波器的输出信号 一阶低通滤波器的传递函数由下式给出： （3-4） LPF的差分方程表示为： （3-5） 高通滤波器（HFP）的传递函数如下： （3-6） 高通滤波器（HFP）的常系数线性差分方程表示为： （3-7） 式中x(nT)是HPF的输入信号，p(nT)是HPF的输出信号，y(nT)是LPF的输出。 2)求导数 经过滤波后，我们要开始检测QRS波的峰值，首先对信号求导数，求导数的作用是求出波形各部分的斜率。 求导过程的传递函数如下： （3-7） 求导运算也可由下面的公式得到： （3-8） 式中x(nT)是求导函数的输入， y(nT)是求导函数的输出。 3)平方运算 求导后是进行平方运算，将求导的结果进行平方，使斜率的差值进一步扩大，便于找出峰值点。 （3-8） 式中x(nT)是平方运算的输入，y(nT)是平方运算的输出。 4)峰值的检测 峰值的检测通过一个自适应的阈值，当信号在某一个固定时间间期内变化时，峰值是一个确定的值，通过下述算法，阈值不断调整以检测峰值，峰值检测到后，还要鉴定一下是信号峰值还是噪音的峰值。 移动窗口积分可用下式表示： （3-9） 式中：N是窗口中的信号样本点数，x(nT) 是积分运算的输入，y(nT)是积分运算的输出。 窗口的大小以及N的选择很关键，