基于Pan-Tompkins算法ECG波形检测方法研究.pdfVIP

与波形的检测241

可以使用一种简单的峰值搜索算法来检测每一个ECG波形。为了获得相应的

QRS波位置，滤波器引入的净延迟应从检测到的峰值位置中减去。

请注意，不能直接在ECG信号上进行峰值搜索：QRS波并不总是一个心

动周期中的最高波形，而且干扰信号很容易影响搜索过程。同时请注意，本

图示中的ECG信号在计算导数之前已被滤波至有限带宽90Hz，并且没有

基线漂移。

解决方案2：Murthy和Rangaraj[201]提出了一种基于加权平方一阶导数算子和移动

平均（MA）滤波器的QRS波检法。在该方法中，定义了一种滤波导算子，其形式为

N

X

2

g(n)|x(n−i+1)−x(n−i)|(N−i+1),(4.4)

1

i=1

其中x(n)为ECG信号，N为计算一阶差分的窗口宽度，差分值被平方并由因

子(N−i+1)加权。该加权因子从当前差分到时间上早N个样本的差分之间

线性递减，起到平滑作用。为了进一步平滑结果，还使用了对M个点进行处

理的MA滤波器，以获得

M−1

1X

g(n)

g(n−j).(4.5)

M1

j=0

在采样率为100Hz的情况下，滤波器窗口宽度被设置为M=N=8个样本。

该算法为每个QRS波群一个峰值，并抑制P波和T波。

在诸如g(n)这样的处理信号中寻找峰值，可以通过一个简单的峰值搜索算法

来实现，步骤如下：

1.扫描信号g(n)中可能包含峰值的一部分，并确定最大值gmax；或者让

gmax表示g(n)在整个可用时间段内的最大值。

2.将阈值定义为最大值的一个分数，例如Th=0.5gmax。

3.对于所有g(n)Th，选择那些对应的g(n)值大于某个预定义数量M的g

(n)前后样本的样本，即：

{p}{n|[g(n)Th]AND

[g(n)g(n−i),i1,2,...,M]AND

[g(n)g(n+i),i1,2,...,M]}.(4.6)

集合{p}按上述方式定义，包含g(n)中峰值的索引。

DETECTIONOFEVENTSANDWAVES241

Asimplepeak-searchingalgorithmmaybeusedtodetecteachECGbeat.Thenet

delayintroducedbythefiltersshouldbesubtractedfromthedetected