一种EASI十二导联多功能Holter系统

　　心血管疾病，又称为循环系统疾病，是一系列涉及循环系统的疾病，循环系统指人体内运送血液的器官和组织，主要包括心脏、血管（动脉、静脉、微血管），可以细分为急性和慢性，一般都是与动脉硬化有关。这些疾病都有着相似的病因、病发过程及治疗方法。实际上，心血管疾病是由心脏学家、胸部外科医生、血管外科医生、神经学家及介入性放射学家来治疗。心电图指的是心脏在每个心动周期中，由起搏点、心房、心室相继兴奋，伴随着心电图生物电的变化，通过心电描记器从体表引出多种形式的电位变化的图形（简称ECG）。心电图是心脏兴奋的发生、传播及恢复过程的客观指标。

　　Holter即“动态心电图”，俗称“背盒子” 往往心电图难以扑捉到有效的诊断依据，但患者又有明显自觉症状，所以多数情况下医生建议她们来进行心电的动态监测。 由于Holter是个英文名字，不好理解，但是记录仪是背在身上的，很多人就把它形象的称其为“背包”。传统基于单片机的Holter能实现心电信号的显示、传输及简单的判病报警，但受限于单片机的运算速度，难以采用检测准确但耗费资源相对较多的算法实现对心电数据的在线分析。用TMS320VC5509A为处理系统，搭配具有高共模抑制比的EASI十二导联心电信号采集方案，实现一个具有实时检测、无线传输和24小时心电数据记录功能的多功能Holter系统。

　　1 EASI十二导联系统

　　EASI十二导联系统是临床医学诊断中常用的心电导联系统。EASI心电图是一种衍生的十二导联心电图，优点是可以通过5个心电电极产生三路双极导联数据，并推导出极近似标准十二导联的心电数据。电极贴放位置如图1所示。

　　EASI十二导联系统可以通过数学公式：Ld=a（A-I）+b（E-S）+c（A-S）转化后得出十二导联心电信息。其中A-I、E-S、A-S为三路双极导联数据，能分别反映心电向量在冠状面、矢状面和额状面的心电变化，a、b、c为经过大样本研究后得出的相关系数。

　　2 系统硬件设计

　　2.1 系统方案

　　系统框图设计如图2所示。

　　本Holter系统主要由EASI十二导联心电采集电路、TMS320VC5509A、外围大容量存储器、液晶显示模块、Zigbee模块和电源管理部分组成。电极从人体体表获取生物心电信号，低通滤波后接入带通放大电路进行两级放大及A/D转换，A/D转换后的数据除了能直接存储在MicroSD卡以便离线分析外，还可由TMS320VC5509A运行经二次微分小波变换算法对心电信号进行实时检测，或者通过Zigbee传输方式将心电信号传输到远程监护中心进行实时监护。

　　2.2 带通放大电路

　　心电信号幅度一般在0.5 mV～4 mV，频带通常为0.05 Hz～100 Hz，频带范围不宽，人体体表阻抗比较大，容易吸取较多共模噪声。为了使心电信号受到的干扰更小，本系统的心电放大电路采用两级放大的架构，体表两点的电位差为心电信号源，前级放大使用高输入阻抗的仪表放大器INA326，采用差分输入，把心电信号放大5倍；后级放大采用高运放OPA2335设计带通放大电路，把0.5 Hz～106 Hz的信号放大200倍。

　　2.3 右腿驱动电路

　　右腿驱动电路是专为克服50 Hz工频共模干扰、提高共模抑制比（CMRR）而设计的，原理是采用以人体为相加点的共模电压作并联负反馈，其方法是提取前级放大电路中的共模电压，经驱动电路倒相放大后再加回人体右腿上。

　　2.4 TMS320VC5509A系统

　　Holter是便携式设备，要求保证至少有24小时连续的医学应用，同时要运行心电检测等复杂算法。因而作为系统的DSP必须具备低功耗和高性能的特点。本文使用TMS320VC5509A作为系统运算和控制，系统配置DSP的McBSP接口与模数转换芯片ADS8325和CC2430 Zigbee模块进行通信，通过SD控制器与MicroSD卡连接，通过GPIO口与液晶模块相连接，设置3个由按键控制的外部中断，用于多功能切换。

　　3 软件设计

　　3.1 DSP程序设计

　　DSP系统初始化需要完成以下操作：初始化CPU、外设、引脚功能与中断，通过判断中断源确定是否存储在MicroSD卡上；然后开始A/D转换采集心电信号，根据此前所选择的方式通过DMA保存采集到的数据；同时数据累积到特征提取所需后对其进行处理，判断是否有异常。

　　3.2 算法设计

　　ECG信号中QRS波群的检测是对ECG信号自动分析诊断的前提和基础，只有在R波标定后才可做进一步分析。本系统采用墨西哥草帽小波变换进行R波检测。该小波母函数是高斯函数的二阶导数，心电信号的特征点与小波变换的模极大值点存在对应关系。式（1）为小波母函数函数表达式：

　　具体检测方法为：在小波变换尺度4上，分等长区间分别求模极大值，再对这组模极大值求均值，将该均值二分之一作为阈值，求出过阈值的连续区间中极大值为R波的相应位置，再修正时移。此时与尺度4相应时延为20点。

　　为了检查算法的有效性，本实验采用国际通用的MIT-BIH数据库进行测试，结果如表1所示。

　　经实验验证，由于噪声在小波变换的第3、第4尺度上已得到抑制，所以系统中所采用的方法可以有效地从噪声干扰中识别出心电信号中R波的位置，并且识别准确率达到了99.83%。

　　在完成R波识别流程后，分别以R波位置为起始点，向前在长度为0.04 s的区间中搜索模极小值点位置，以对Q波进行定位。对S波进行识别的基本流程与Q波相似，不同点是向R波后向检索，并且由于S波延续时间较Q波长，搜索区间长度为0.06 s。

　　3.3 上位机软件设计

　　本文的心电远程实时监护界面采用LabVIEW虚拟仪器编程语言设计，主要功能为实时从串口采集心电数据，切换显示十二导联数据，分析和存储等功能。

　　本系统中，Holter终端节点采集分析的数据应用ZigBee无线协议传输至网络协调器节点。网络协调器节点将接收到的用户信息数据进行融合处理，通过串口传送到PC机上。

　　4 系统测试

　　4.1 心电采集和预处理

　　按照图1的EASI导联系统电极位置所示，通过心电电极片连接人体和系统前端电路，采集心电信号。图3为实际采集的AI通道的心电信号并在CCSv3.3上显示的心电信号片段。

　　4.2 算法实现

　　在TMS320VC5509A采用小波变换的方法，在小波变换尺度4上对心电信号进行实时特征提取，并对QRS波群的各个特征点参数进行检测计算。小波变换尺度4如图4所示。

　　如图4所示，圆圈标记为R波在尺度4上对应的位置，三角形标记为Q波，矩形标记为S波。然后据此再对其他参数如P波、T波及其端点检测。

　　4.3 远程监护

　　在上位机运行心电远程实时监护界面，可同时实时显示用户的心电数据及根据EASI导联推导出的十二导联数据。

　　本文实现了一个基于DSP的多功能Holter系统。硬件系统、EASI十二导联心电采集电路、心电识别算法、Zigbee无线传输通信和上位机监护程序等都已调试成功。Holter系统能实现基于小波变换算法实时检测从心电模拟前端采集的心电信号，并通过Zigbee实时无线传输到PC机进行心电远程实时监护与十二导联心电数据的推导及实时显示，达到预期目标。