基于dcg心电信号的r波检测新算法.docxVIP

基于dcg心电信号的r波检测新算法 r波的正确定位有助于预防和诊断心脏猝死的临床意义。动态心电图（cdg）作为一种诊断心脏疾病的工具，在临床上没有创伤性。心跳速率变化分析、心率减速法分析、宫理化率振动分析和t波电交替分析是基于dg的有效报警心律失常方法，这些分析方法取决于正确的r波检测。 非实时的R波检测方法包括小波分析法和神经网络等方法, 虽然检测效果较好, 但运算量大、原理复杂, 不易重现;而实时的R波检测方法, 其检测效率高、速度快、原理简单、易于重现, 已经成为近几年的研究热点. 文献提出了两种实时的R波检测算法——最大值双重搜索技术和差分运算方法, 两种算法的共同点是利用心电信号中R波斜率大的特点, 以相邻点之间的斜率为主要特征对R波进行检测.上述算法对噪声干扰比较敏感, 在检测R波之前需要对原始信号进行基线漂移校正, 对于采样频率较低的DCG心电信号的处理效果不够理想. 笔者提出了一种基于DCG心电信号的R波检测新算法, 利用了R波波峰斜率大、幅度高的特点, 提取了平均双向斜率和相对高度两项特征, 对MIT-BIH Long-Term ECG数据库中的全部7 个样本实现了R波波峰点坐标的定位.通过与另外两种实时的R波检测算法的对比, 验证了本文所提算法的可行性. 1 r波的准确定位 一个正常心电图的每个心动周期由P波、Q波、R波、S波、T波以及U波组成 (如图1 所示) .其中特征最为显著的是R波, 它充分反映了心室收缩时心脏的电行为, 因此解决心动周期提取问题的关键就是R波的准确定位. 观察发现, QRS波群中的QR波和RS波斜率高、幅度大.利用这两个特点, 对心电信号这两方面的特征进行提取可以突显R波波峰的位置, 从而达到准确定位R波的目的. 1.1 斜率步长的设定 在对心电信号的处理中, 定义斜率步长为用于求取斜率的两采样点间的采样点间隔. 设N (i) 为心电图中采样点i的纵坐标 ( i点对应的电压值) , f为心电信号的采样频率.采样点i和采样点j之间的斜率为 由于采样频率f是一个定值, 因此相邻采样点间的时间间隔相等, 于是式 (1) 可以简化为 在式 (2) 中, 采样点i和j的斜率步长就是i -j . 因为DCG心电信号的采样频率较小, 如果斜率步长设置得不够合理, 容易误将伪R波 (如图2 (a) 所示) 或畸形T波 (如图2 (b) 所示) 误识为R波.为了防止这两种情况的发生, 需要对斜率步长进行合理设定. 本文所采用的MIT-BIH Long-Term ECG数据库所有心电信号的采样频率均为128, Hz, 文献表明, 一般QRS波的时间是0.06, s左右, 也就是说QR斜率步长与RS斜率步长的和大约在7 个采样点左右 (0.06÷ (1÷128) =7.68) , 为突显QR斜率和RS斜率高的同时保持R波的对称性, QR坡的斜率步长和RS坡的斜率步长统一取3. 1.2 平均双向斜率的计算 为突显出R波两侧斜率高的特性, 本文定义了前向斜率、后向斜率以及平均双向斜率3 个概念. 设t为斜率步长.对采样点i而言, 点i与点i -t之间的斜率k1(i) 为点i的前向斜率, 即 而点i与点i +t之间的斜率的相反数k2(i) 为点i的后向斜率, 即 点i的前向斜率k1(i) 和后向斜率k2(i) 的平均值m (i) 为点i的平均双向斜率, 即 将式 (3) 和式 (4) 代入式 (5) , 得 由于对每个点i而言t都相等, 因此在计算平均双向斜率时, 式 (6) 可简化为 如图3 所示, 在平均双向斜率序列中, R波波峰的位置更加突出, 识别起来更加容易. 1.3 初始平均双向斜率阈值的确定 定位第1 个R波波峰是定位全部R波前需要解决的首要问题.只有确定了第1 个R波波峰的位置才能获得R波形态特征的阈值, 进而迭代地把R波检测进行下去. 第1 个R波波峰的定位依赖于初始平均双向斜率阈值的确定:首先在信号开头的至少含有2 个R波的信号段中求取最大平均双向斜率x ;然后选择一个合适的阈值r (如0.6) , 以xr为平均双向斜率的阈值重新搜索这一段信号, 可以定位第1 个最大平均双向斜率的位置loc.在该点附近 (loc-5 ~loc+5 ) 搜索到的最大峰值top (如式 (8) 所示) 就是第1 个R波波峰的位置. 1.4 相对高度序列 为突显R波幅度特性的同时又能克服某些具有基线漂移 (如图4 (a) 和4 (b) 所示) 的信号在R波检测中的干扰, 本文定义前向相对高度、后向相对高度以及相对高度的概念. 对采样点i而言, 点i与点i -t之间的纵坐标之差h1(i) 为点i的前向相对高度, 即 点i与点i +t之间的纵坐标之差h2(i) 为点i的后向相对高度, 即 点i的前向相对高度h1(