心音信号的EMD算法分析 摘要：一个新的无心电图参考的无监督、低复杂度的检测S1和S2心音成分的方法被提出。该心音分割算法将心音信号分为四个部分：第一心音，收缩期期，第二心音和舒张期。该算法使用经验模式分解（EMD）在时间域产生心音的强度包络。所提出的算法的性能已使用100个数字心音记录包括正常和异常心音的14000个心脏周期评估。该算法已显示出超过90%的正确率。 关键词：心音描记法，听诊，第一心音，第二心音，收缩期，舒张期杂音，经验模式

分解，分割。 一 引言 心血管系统造成的杂音和畸变的许多病理条件听起来多比之前他们反映的其他症状多，如心电图（ECG）信号的变化[ 1 ]。虽然心电图和超声检查广泛用于心血管病的防治诊断，首先由医生进行了听诊心脏声音分析的老年艺术来评价心脏的功能状态。然而，诊断很少单独基于听诊，由于该听诊是主观的而且高度依赖于解释器的技能的事实。然而普通医生还是要采用有效和客观的手段做一个简单的初步检查，在疑似病人可能被送到心脏科医生作进一步检查之前。 听诊观察者试图分开听和分析心音，然后合成心脏功能。一个心动周期中，应该确定的重要组成部分是：在第一心脏声音（S1），收缩期期间，第二心脏声音（S2），和舒张期在这个序列中的时间。重要特征应该被量化，包括节律、定时时刻和心脏声音分量的相对强度，S2的分裂，杂音或其他额外的声音的存在以及杂音或其他额外的声音的时间、强度和质量。虽然不同声音定性描述很多是可用的，但是仅通过听去量化其性能是很困难的。有人可能会得出这样的结论：客观有效的基于心音成分的定量特征可用于做出可靠的诊断。 心音分割算法在文献中大致可分为两种主要的方法：需要参照心电图同步分割和不需要这样的。后者可以进一步分为有监督和无监督的方法。有心电图参考的方法，首先分别检测QRS波和T波以定位SI段和S2段[2]。低质量的心电信号中，T波并不总是清晰可见。在这种情况下，S2的声音可以通过非监督分类[ 3分]。为了避免数据采集时额外的硬件要求和繁杂的接线安排，许多研究者试图通过几种信号处理和统计信息的方法找出S1 S2，而不用心电图作为参考。在这种情况下，已经提出了一些监管技术，如人工神经网络[4]和决策树[5]。另一类方法涉及如envelogram的无监督的方法[6]，频谱量化方法[7]，以及使用梅尔频率倒谱系数的自组织映射[8]。据观察，在实践中这些方法并不适用于所有类型的心音（如严重的杂音情况下）。众所周知，关于人工瓣膜产生的心音的分割，这些声音是取决于几个因素，如手术方法，植入的位置和人工瓣膜类型。 在实践中，上述的方法不能提供所需的不变性，为了能够适用于各种心脏瓣膜植入患者。 在当前国家的最先进的心脏声音的分割算法应用的最可靠的和不变的特点是含蓄的或明确的S1-S2间隔规律。 在本文中，针对S1和S2心音检测，提出了不依赖于心电图参考或病人的任何其他先验信息的新颖方法。 在该算法中，利用EMD算法将PCG信号分解为一些特定的功能，称为基本模式分量函数（IMF），提取一些PCG信号基本特征。 首先，运用EMD的IMF分量的包络将所有可用的心音分割成声音裂片，然后，划分S1 S2心音在每个心脏周期，通过检查两个音段之间的最近间隔为收缩间期（S1-S2期）。本文的结构如下：在第2节介绍了S1和S2了心音检测方法。第3部分和第4部分对结果进行讨论，最后在第5部分指出一些主要结论。 二、方法 A．数据收集 2010年至2011年，我们在德黑兰Modares医院的儿科门诊收集的心音。数据是从共100例，其中8例正常，年龄是1~26岁。报告出版和相应的图像的书面知情同意书从患者处获得。同意书的副本可供本期刊的首席主编审查。病人（除正常对照组）有心脏杂音历史，其中包括30例室间隔缺损（VSD），12例房间隔缺损（ASD），6例同时有ASD、VSD ，12例法洛四联症（TOF），6例肺动脉瓣狭窄（PS），34其他种类的疾病。在所有的情况下，经超声心动图诊断，一些病例的病史、诊断结果总结于表1。 数目 病人年龄 病人体重 心脏病 案例1 6 35 室间隔缺损 案例2 7 19 房间隔缺损 案例3 9 24 法洛四联症 案例4 12 26 肺动脉瓣狭窄 案例5 5 14 动脉导管未闭 案例6 13 35 主动脉瓣狭窄 案例7 24 85 健康 表1 一些研究患者的心音记录 B.预处理 电子听诊器记录的心音信号通常包含高频率噪声，因此，预处理是必不可少的。先滤波进行信号消噪，其次是归一化和分割。具体步骤如下： 滤波：由于心音信号主要是小于800 Hz，应用数字有限脉冲响应（FIR）设计巴特沃斯低通滤波器。 归一化：不同心音信号的振幅均归一化和限制在[-1 1]这个尺度。 归一化方程如下所示： [n][n]max|[n]|PCGnormPCGxxx

（1） 一例正常心音的滤波和归一化记录如图1所示

图1 滤波和归一化后的一个正常心音 C.经验模式分解（EMD） 心音信号首先被分解成的IMF分量。为了获得IMF分量，先找到信号的局部最大和最小值。局部最小值和最大值的包络线分别由三次样条拟合。让1(t)m的表示为这两个包络线的平均值，由原始信号减去1(t)m作为如下： (t)x-1(t)m=1(t)h （2） 把1(t)h作为一个新的信号并且重复上述，直到1(t)h满足IMF的标准[9]。所得到的信号是第一IMF分量，定义为1(t)c，剩余信号是1(t)r，其中 1(t)r=(t)x-1(t)c （3）

把1(t)r视为一个新的原始信号，并且在迭代过程被再次执行以提取IMF分量，直到第k个残差信号(t)kr成为一个单一的功能，这意味着没有更多的IMF可以被提取。因此(t)x的表达式为：

(t)x=1(t)(t)kkiicr （4）

其中1(t)c到(t)kc是信号的k个IMF分量。 EMD方法会使得信号更加对称，通过自适应消除骑行波和分解信号。 在本实验中，IMF5被提取。该IMF显著减少了杂音和噪音的影响。其余的IMF分量对检测主要声音并不是很有效的方法。所提取的IMF5和它的绝对值示于图2。

图2. IMF5 （a）正常PCG （b）EMD提取的IMF5 （c）IMF5的绝对值 D.分割 在这部分中，心音信号应划分为周期。四个条款在此步骤中需确认：第一心音（S1），第二心音（S2），收缩期和舒张期。为了检测这些条款与噪声的干扰抑制，IMF5的香农能量（SE）的使用和计算方法如下[10-14]： 2255[n]=-[]log[]imfimfSEnnxx [5]

该方程表明，香农能量的特征是抑制信号的低幅度分量[11]。因此主要成分S1 S2将被分出。与其他算法相比，如绝对值或香农熵，香农能量对分割心音杂音更好[11]。 IMF5的香农能量及其包络计算和绘制在图3。

图3 （a）IMF5的香农能量 （b）包络 心音杂音很容易被提取，影响S1和S2的识别结果。然而，更高的阈值对于处理杂音会更容易重新检测到丢失的S2。S1和S2由以下三个步骤确认[13]： 1.如果在50毫秒内检测到两个高于阈值的峰值，则低能量的一个会被淘汰。 2.对于峰之间的每个间隔，长度较短的间隔不超过300毫秒被表示一个收缩期，而另一个超过300毫秒的间隔是舒张期。不确定的间隔进行了注明。 3.对于那些不确定的时间间隔，第二阈值设置为了发现可能没被识别的S1和S2。 这些步骤用于确保S1和S2可以被正确识别，然后确认收缩期和舒张期。 然后运用迭代识别方法检测可能被丢失的组件。迭代识别方法检查检测点是否周期性发生。对于S1S2应该被检测的位置，该方法可以确定是否有丢失的点，检测到的S1S2如图4所示。 图4 （a）心电图作为参考评价该方法的准确性 （b）提出的方法检测到S1S2 三．结果 从德黑兰的Modares医院收集的患者心音记录用于验证该方法的识别准确率。该数据库包含了100条记录，信息表提供在表1。 利用迭代识别患有各种心脏病的7个案例如下。 临床病例1 图5显示一个健康的24岁男子无心脏杂音的心音样本。为了评估该方法的准确性，提供了心音及相应的心电图。注意检测到的S1和S2。 图5 健康、无杂音的心电图和心音图 （a）心电图 （b）检测到的S1 和S2心音 临床病例2 图6给出一个患有室间隔缺损的6岁女孩的心音样本。为了评估该方法的准确性，提供了心音及相应的心电图。值得注意的是，尽管有很严重的收缩期杂音，S1S2也被检测到了。

图6 VSD病例的心音样本 （a）心电图 （b）检测到的S1 和S2心音 临床病例3 图7给出了一患有房间隔缺损和轻度室间隔缺损的7岁女孩的心音样本。为了评估该方法的准确性，提供了心音及相应的心电图。请注意，在收缩期存在房间隔缺损杂音，然而，还是检测到了S1和S2。 图7 ASD案例心音图的心音样本 （a）心电图 （b）检测到的S1 和S2心音 临床病例4 图8提出的心音样本是来自一个患有主动脉瓣狭窄的13岁女孩。为了评估该方法的准确性，提供了心音及相应的心电图。请注意，尽管收缩期有杂音，S1和S2也被检测到了。

图8 AS案例心音图的心音样本 （a）心电图 （b）检测到的S1 和S2心音 临床案例5