摘要本文采用从单导到多导的检测方法，首先利用小波变换实现单导QRS波的检测，在此基础上，利用位置相关法进行多导QRS波的检测，并利用心电信号的2¹尺度小波变换的平方值来确定QRS波的起止点。经过大量数据的检测证明取得了很好的效果。

Identification of QRS Complexes and Their Boundaries in Multi-Lead Synchronized ECG

Yu Dekuang, Yin Bingsheng, Yang Yi, Li Benfu and Nong Debin

(The Fisrt Military Medical University,Guangzhou 510515)

Abstract

In this paper we present an algorithm for detection of QRS complexes from singlelead to multi-lead in ECG.We first detect QRS complexes in single-lead ECG using wavelet-transformation of signals on the bano of ahich we then detect QRS complexes in multi-lead ECG according to their position correlation; and we determine onsets and ends of QRS complexes with the square value of 2¹ scale wavelet-transformation of the ECG signals.Tests of large quantities of data have proven the algorithm satisfactory.

Key words:Multi-lead synchronized ECG, Identification of QRS complexes, Wavelet transformation, Onsets and ends of QRS complexes

0前言

随着计算机在医学中的广泛应用，心电图的计算机分析已成为心电图发展的一个重要方向。心电图的计算机分析系统主要由波形参数测量和诊断分类两个部分组成。其中的关键部分是波形的定位及其参数的测量，尤其是QRS波的定位和起止点的确定，其检测的准确率直接影响计算机分析的结果。QRS波是心电图中最具特征的波形，波形的检测都由QRS波开始，心率的计算和分析以及T波和P波的检测都基于QRS波的精确定位和范围的确定。人们对单导QRS波检测的研究较多，取得了比较大的进展。目前同步多导心电图在临床上的应用越来越广泛，多导同步心电图的波形检测也日益受到人们的重视。

1多导同步心电图的QRS波检测原理和方法

心电图的QRS波代表心室的去极化皮，由动作电位的峰电位所产生，时间宽度不大，是心电信号的高频部分^［1］。在同步十二导心电图中，各导的心电波形不相同，差异可能很大，有的比较明显，幅度较大，变化大；有的QRS波比较低平，和T波与P波的频率相近，幅度也接近，这样就给同步多导心电波形的识别带来一定难度。这里我们采用从单导检测到多导检测的思想^［2］，利用同步多导心电波形的位置相关性，在单导检测基础上进行多导的进一步判断。过程如下：

如图1所示，多导QRS波检测过程包括预处理和小波变换、单导QRS波判断、多导QRS波判断和多导QRS波位置的修正。

图1多导QRS波检测框图

(1)预处理：利用简单整系数滤波^［3］进行低通和带阻滤波。心电信号易受肌电、呼吸和工频干扰的影响，去除这些干扰有利于更好地进行波形检测。(信号采样率为250点/秒)

①带阻滤波(或称陷波)的z变换为：

其信号计算公式为：

y(n)=2y(n-5)-y(n-10)+［1600x(n-195)-3198x(n-200)+

1600x(n-205)-x(n-400)-x(n)］/1600

式中x(n)为数字心电信号，y(n)为滤波后信号。用此式不但能够去除50Hz干扰，而且对基线飘移也有抑制作用。

②低通滤波的z变换为：

其信号计算公式为：

y(n)=(3x(n)-3x(n-3)+3x(n-6)-x(n-9))/27+3y(n-1)-3y(n-2)+y(n-3)

式中x(n)为数字心电信号，y(n)为滤波后信号。

(2)信号变换与单导QRS波判断

此处我们对滤波后的数据进行二进小波变换，用Mallat算法实现，小波函数选用的是三次B样条小波^［6］。其实现公式如下：

其h_k和g_k值如下：

h₀=0.3750, h₁=02500, h₂=0.0625, h₃=h₄=0,

g₀=0, g₁=0.59261, g₂=0.10872, g₃=0.01643, g₄=0.00008,

h_n=g_n=0(n＞4) h_n=h_-n g_n=-g_-n

S₂₀f(n)=dn为数字心电信号，代入系数h和g的值，就可以求出信号的多尺度的小波变换。

利用上面公式求出2¹至2³的小波变换，再对2²和2³尺度的小波变换进行幅度阈值判断，判定QRS波^［4，5］。

(3)多导QRS波判断

由于同步十二导心电图QRS波的位置有一个对应关系，所有导联的某一确定QRS波的位置有所偏差，但在某一窗口内是比较集中的，针对这个特点我们把十二导信号的QRS波判断结果转换到一个序列上(0表示此处未被判为QRS波，〉=1表示此处被判为QRS波)，选择固定的窗口在此序列内进行搜索，如果在此窗口内有大于一定导联数的信号被判定为QRS波，则确定该窗口内有一个QRS波，并把该窗内能检测出QRS波的导联所检测出的QRS波的位置的平均值确定为QRS波的位置。具体步骤如下：

①设各导联的QRS波检测结果为p_i(k)(k=1…n_i)，p_i(k)为检测出的QRS波的位置，n_i为一导检测出的QRS波的个数，i为某一导联(i=1…12)。

②把各个导联的结果转换到位置序列上：

0表示此处未被检测为QRS波，1表示此处被检测为QRS波，i表示第i导联。

③把所有十二导的QRS波检测结果转换到同一位置序列上：

经过上面换算，所有的十二导检测的结果在一条位置序列上得到反映，R(n)为0表示该处未被检测为QRS波，为1表示有该处有一导检测为QRS波，为m表示有m个导联检测该处为QRS波(m为整数)。

④用固定的窗口Ts(取120mS)去搜索该位置序列：

若G(n)值大于某一阈值(我们取为8)，则定该窗口内有一QRS波，并把在该窗口内能检测出的QRS波的位置的平均值作为十二导的QRS波位置。

⑤对各导联的QRS波检测结果进行修正：根据④所确定的十二导的QRS波的位置pos(n)，分别在各个单导的窗口(pos(n)-Ts/2,pos(n)+Ts/2)中搜寻，若窗口中存在有单导的检测结果，则此单导检测结果就是该导联的某一QRS波的检测结果；若窗口中没有单导的检测结果，则在该窗口内找出小波变换值最大的点，找出和这点组成极大－极小值对的过零点，把此过零点定为该导联上的该QRS波的检测结果。

2QRS波的起止点的确定的方法和步骤

三次B样条小波具有很好的时频局部化特性，我们利用2¹尺度的小波变换的平方值进行QRS波的起止点的确定，由于小尺度的小波变换所突出的是信号的相对高频的信息，即QRS波的部分，而相对低频的部分，如T波和P波，则被压低，对于P－R段和ST段，其小尺度的小波变换值很小，针对小波变换值有正有负，并且利用此小波变换值确定QRS波的范围往往偏大的不足，我们选用其2¹尺度的小波变换的平方值，一方面实现对信号的放大，增大了变换结果的对比，另一方面，把信号由正负皆有变成全部的正值，便于判断。

设尺度为1的小波变换结果为W₂₁f(n):

<