摘要本文研究了一种准确实时提取接受反搏治疗病人心电信号特征值方法，分析了其干扰噪声的形式，设计了数字滤波器对干扰进行抑制，研究了心电QRS波检测匹配模板，运用模板匹配技术完成QRS波检测，结合子波多尺度分析对心电信号进行特征提取，根据提取的特征完成匹配模板修正，完善了反搏治疗病人心电信号特征提取，将有利于保证体外反搏治疗效果。

A Real-Time Approach to Extract Parameters Form a

Patient's ECG on EECP

Li Haiyun, Duan Yuanmin

(Dept. BME, Capital University of Medical Science, Beijing 100054)

Zheng Zhensheng

(Dept. BME, Sun Yat-sen University of Medical Science, Guangzhou 510089)

Abstract

A precise and real-time approach to extract parameters from a patient's ECG under the condition of the Enhanced External Counterpulsation(EECP) was developed. Different types of noises were analyzed. A digital filter was designed to suppress noises. A match themplate was designed to detect the ECG QRS. A new method has proposed to extract the feafures of ECG by wavelet multiresolution analysis. According to extracted features, the template was modified. The approach can improve extraction of parameters from a patient's real-time ECG, which is helpful to improve the results of treatment of EECP.

Key words:EECP;ECG;Wavelet

0引言

体外反搏系统是一种应用计算机控制的无创伤的机、电、气辅助循环装置。通过包裹在患者的下肢和臀部的密封气囊以特定模式对肢体加压，改变正常血液流场分布，保证在肢体、臀部受压引起血液返流回主动脉瓣处时正好是主动脉瓣关闭的瞬间，并充分利用心脏舒张期的整个时间，使施加于血管的压力保持足够长，最大限度提高舒张压，保证反搏血流以最充分的时间向缺血器官进行灌注，而又不加重心脏的负荷。另一方面，反搏排气时间的确定保证在下一个心脏收缩期之前解除对血管的压迫，使收缩压降至最低，最大限度地减轻心脏射血期阻力，从而达到最佳的反搏效应^［1］。实施体外反搏必须使充排气时序跟踪反搏期间心动周期的变化。以保证心动周期变化时，反搏的时序控制始终符合生理要求。由图1可知：心电信号T波结尾一般可视作心脏收缩期即将结束，与主动脉瓣关闭时刻非常接近^［2］。T波末端与主动脉瓣关闭时刻之间最大偏差＜40ms^［3］。

图1心脏的力学过程与电学

过程的时序关联

(图1中1=右心房收缩开始；2=左心房收缩开始；3=左心室收缩开始；4=右心室收缩开始；5=右心室射血开始；6=左心室射血开始；7=左心室射血终止；8=右心室射血终止；交叉线区域表示等容收缩，画点区域表示心室射血。)

故可以QRS波为基准，计算QT时间(为心电信号中从QRS波的R波到T波终点的一段时间)，反搏控制系统根据不断检测得到的QT间期自动跟踪调节反搏参数。从而保证充排气时序可以在反搏期间不断跟踪心率的变化来进行调整，以保证心率变化时，反搏的时序控制始终以QT作为反搏控制参数，可获得反搏的最佳效果而又不危害病人的健康。因此准确检出每一个QRS波是实施体外反搏的关键。比一般情况下的心电检测分析增加了不少难度 。我们综合了数字滤波、模板匹配、子波分析技术完善了反搏治疗病人心电信号实时处理。

1整系数心电信号数字滤波

采集到的反搏病人的心电信号常常伴有各种干扰，主要来自工频、呼吸运动、反搏引起的电极与皮肤的相对运动产生的噪声，必须进行滤波处理，在体外反搏治疗实施过程中，要求病人的心电信号的检测处理及反搏控制均在一个心动周期内完成，因此对数字滤波器的要求：精度适中，占时尽可能短。以适应体外反搏治疗装置的在线工作。选用整系数线性相位递归滤波器能满足要求。权系数为1的动平均滤波器具有低通滤波器特性，其系统函数为：

(1)

但其旁瓣较大，为减小其旁瓣的幅度，我们采用高阶的归一化形式，设定其系统函数为：

(2)

其差分方程为：

(3)

滤波结果见图2。

图2心电数字滤波

2心电信号的QRS波检测及子波分析

2.1QRSB波的检测

设计一个三参数匹配模板TLV=［A_T,D_T，Q_T］，其中

A_T：QRS波的幅度阈值D_T:QRS波上升沿梯度阈值Q_T：QRS波的QT时间阈值。

运用TLV，通过不断比较检出QRS波，但有可能将干扰和各种早搏信号以及其他一些心律失常信号误检为可以实施反搏的正常心动周期的QRS，这在体外反搏治疗过程中应严格避免的，会误触发反搏装置工作，给病人的健康造成损害。为克服这一不足，我们采用子波多尺度分析离线处理心电信号，剔除干扰和和各种早搏信号以及心律失常信号的影响，准确提取正常心动周期QT时间(由于体外反搏在每一心动周期内是一严格的实时控制系统，应用子波分析在线提取心电信号特征实施反搏控制没有可行性)作为一个动态跟踪标准，跟踪更新TLV中的Q_T参数，去识别检测到的各种早搏信号以及心律失常信号，一经识别后，中止反搏器的治疗，待信号回复正常再实施反搏。以保证病人治疗过程中的安全。

2.2心电信号的子波分析

子波变换的极值点和零交叉点可用来检测信号的奇异点位置。而且有两个突出特点：受高频噪声影响小，选择适当尺度下的子波分析，可以充分突出信号成分最大衰减高频噪声，即使在强噪声背景下也能准确定位信号的奇异点位置；其二，由于子波变换本质上的带通特性，以低频成分为主的基漂，对信号奇异点定位的准确性影响不大^［4］。

二进离散子波变换为：

(3)

构造ψ(t)使ψ_m,n(t)成为正交子波。可以得到一个函数的子波级数展开式：

(4)

令

(5)

基函数φ_m,n(t)为尺度函数，其定义为：

(6)

所以

(7)

上式以尺度的观点分析心电信号时，以尺度m₀为界限将心电信号分成两部分，m₀以下各尺度作为细化特征的近似；m₀以上各尺度用于基本特征的提取。选用不同的子波，就产生不同的效果。选取平滑函数的一阶或二阶导数可以作为子波，在这样两种情况下，可分别以子波变换的极值点或零交叉点来表征信号，而且用子波变换的极值点表征信号比用零交叉点表征信号有更多优点。我们采用信号子波变换的局部极值来抽取心电信号的特征点，也就是说，所选用的子波ψ(x)应等于平滑函数θ(x)的一阶导数，即要求：