心音信号自动分析与识别的研究现状吴玉春，郭兴明，肖守中重庆大学生物工程学院，重庆（400030）E-mail：wychun690@摘要：分析了当前心音信号自动分析与识别的研究现状，通过分类比较了谱分析、小波分析、神经网络方法以及数学形态学方法等几种用于心音的研究方法各自的特点，讨论了对心音信号进行自动分析、识别的研究状况、研究重点和发展方向。
关键词：心音,谱分析；小波变换；神经网络中图分类号：R318．041. 引言心音信号是包含有重要生理、病理信息的医学信号之一，是临床评估心脏功能状态的一种基本方法，某些心血管疾病发展到出现明显的病理特征之前，会在心音相应组成成分上有所反映，如冠状动脉狭窄会在心音的舒张期产生异常杂音。
但是，由于目前心音的研究中存在以下几个问题，心音信号的自动分析与识别还不能很好的应用于临床，第一，心音的产生机制尚无定论，目前的研究表明，心音是由心脏的瓣膜和大血管在血流冲击下形成的振动，以及心脏内血流的加速与减速形成的湍流与涡流以及对心脏瓣膜、心房、室壁的作用所产生的振动，再加上心肌在周期性的心电活动作用下起刚性的迅速增加和减小形成的振动，经过心胸传导系统到达体表形成的体表心音[1]；第二，目前对于心音自动分析的研究大都处在理论研究阶段，临床上对于心音的诊断大多还是靠医生的经验和主观判断，用听诊器采集心音音频信号，这种诊断主要依赖于医生的经验水平，往往只能是有经验的心脏病专家，才能做出正确的诊断结论；第三，目前心音信号的自动分析技术还不成熟，不全面，一般只能对正常的心音信号做出比较准确的识别，而对于异常心音信号，还不能用先进的信号处理手段进行自动分析及识别。
2.常用方法下面就目前用于心音信号自动识分析及识别的一些常用方法进行分类介绍：2.1 谱分析上世纪90年代以前，心音分析工作者主要采用FFT、AR（自回归模型）和ARMA（自回归滑动平均模型）等谱分析，而且其研究对象也局限在动物心音。
有研究[1]表明：心胸声传导系统是一个时变模型。
因此，采用FFT或参数模型（LPC、AR和ARMA）都是基于平稳性好的分析方法，难以反映出心音信号的动态变化过程，另外，功率谱只能反映幅频特性，而往往忽视了蕴含着大量信息的相频特性。
后来Sava HP等[2]采用改进的前后预算的Prony 方法对用机械膜修复后的心脏心音的谱分布进行了研究，Prony方法可以直接求得信号的幅值、相位、阻尼因子和频率4个要素，能够处理衰减信号、估计阻尼系数。
何莹等[3]将高阶统计量方法应用于心音信号分析，建立了心音信号的AR双谱模型，获得了心音的双谱幅度图，并采用模型参数作为特征参量对心音信号进行了二类模式识别，其实验结果表明，心音信号本质上为非高斯信号，高阶统计量在心音信号分析和处理中有一定的应用价值。
李彬彬等[4]采用改进的HHT算法对心音信号进行分析，李等在EMD筛选过程中引入基于两个参数θ1和θ2的判断标准，来缓解判断条件设置不当引起的误差。
模态幅值：（1）评估函数：（2）筛选停止条件：1）评估函数不大于θ1，其目的是为了保证整个数据段中波动的平均值足够小；2）评估函数不大于θ2，其目的是为了限制EMD分解时局部数据大的漂移。
这两个参数在保证数据整体性的同时，也考虑到局部数据的影响。
利用这种改进的HHT对正常心音信号进行分析，从时频特征、时域特征、能量分布等多角度提取了心音的特征值，为进一步的心音分析提供了有用的方法和数据。
朱冰莲等[5]采用短时傅立叶变化（STFT，又称窗口傅立叶变换）对心音信号进行分析。
通过适当地调整STFT 的参数(窗口的类型和宽度) ,可以比较准确地得到心音信号的时频表示。
当选用窄窗进行分析时,时间分辨率高,可得到心音信号的时域能量包络图,由此可以准确地量度心音成分的时域特征参数；而选用宽窗进行分析时,便得到每个心音成分( S1 和S2) 的频谱,由此可以以高的频率分辨率识别出它们的组成。
朱等首先用窄时间窗进行分析，频谱图具有高时间分辨率，得到了心动周期等时域特征参数。
进而逐渐加宽时间窗，最后得到高频率分辨率的频谱图，可以看出各心音成分的频谱特征。
实验结果表明，不同时频尺度的STFT分析可以很好地描述正常心音信号的时域和频域特征。
2.2 小波分析因为小波变换中的时-频窗具有自适应性，所以可以得到更能反映心音病例特征的信息。
心音信号的小波变换谱在不同尺度上的分布相应于心脏不同部分的振动信息。
尺度较小的部分相应于信号的高频成分，可以表示心脏瓣膜的振动信息、血流本身及其对其它部分相互作用的信息。
而尺度较大的部分相应于信号的低频部分，这是心肌与大血管等心脏的质量较大的部分在瓣膜血流和肌肉的作用下产生的振动。
近年来，国内生物医学工作者应用连续小波变换对心音的频率范围进行估计和对心音做三维分析。
M.Nadir Kurnaz,Tamer Olmez[6]应用小波变换对一个窗口中两个时间段的心音信号进行自动分段并提取特征参数，在这个窗口中实现了两个分析目的，一个是S1和S2的分段，另一个是特征提取。
分段以后，用小波变换分解心音信号，并由第六层分解的细节（高频）系数得到特征向量，从而使用动态设计分析得到最优的特征参数。
D.Kumar等[7]采用小波变换的高频特征量来检测第一心音和第二心音，由于当前的心音自动分析的算法存在受噪音影响大和不适于在无节律状况下进行心音分段的问题，他们提出了一种基于高频标示器的方法来解决这些问题，用快速小波分解来提取高频特征参数来估计瞬时心率，这个标示器是由横穿心脏瓣膜的主要压力差来推动工作的。
通过对人工心脏瓣膜和正常健康瓣膜两种情况下心音的分析，证实了算法的有效性。
李勇[8]等将连续小波变换应用于心音信号的分析中，在基于Mellin变化的连续小波快速算法的基础上提出了同步展宽法，解决了在分析尺度范围大，小波函数和信号中心频率较高条件下引起的计算量大的问题。
但是由于连续小波变换存在“频带”的重叠现象，胡晓、全海燕、王威廉等[9][10]利用离散二进制小波变换及其相应重构原始信号的公式对心音信号进行处理，他们采用具有双正交性和比Daubechies系更好对称性的coiflets小波系，对心音信号进行多分辨分析，将其分解成不同的频带，然后按时间顺序分段计算频段的规一化香农能量，能量大于某一合适阈值的时刻作为S1和S2的候选点，最后结合心音的特点确定S1和S2。
进行信号分析时，在有些情况下，选用不同的小波基函数，会得出截然不同的结果。
选取合适的小波基函数，会使计算出的小波系数能够最好地反映信号的特征。
2.3 人工神经网络人工神经网络（ANN）是对人脑的抽象、简化和模拟，反映人脑的基本特性。
它是由简单的处理单元所组成的大量并行分布的处理机，这种处理机具有存储和应用经验知识的自然特性，其自学习自适应特点可用于心音信号的自动识别。
网络输入特征矢量的确定，在神经网络模式识别方法中具有重要的作用。
输入特征的选取直接关系着网络识别方法及结果的准确性，对于生物医学信号而言，常用的网络输入特征有时域波形、线性预测系数、预测误差、带通滤波器组、短时富里叶变换谱、倒谱以及小波变换谱。
Zamri Mohd Zin [11]等采用积分小波变换和离散小波变换对心音信号进行处理和讨论，将通过积分小波变换得到的心音信号特征参数作为人工神经网络的输入向量，对多层感知神经网络进行训练和识别，进而得到比较好的识别效果。
姚晓帅等[12][13]提出了一种基于概率神经网络 PNN 的心音信号识别算法，首先采用小波变换对信号进行预处理，然后采用平均Shannon 能量变换提取心音信号包络波形，进而对包络进行第一心音与第二心音波峰的提取与定位，最后定义心音信号特征矢量，并将其作为概率神经网络的输入矢量，进行基于概率神经网络的S1、S2识别。
叶学松[14]采用小波变换的方法对心音在时域和频域上展开了局部化分析,把不同频段上整周期心音信号的平均能量与舒张期的平均能量的能量比值以及人体的个体特征参数(如年龄、抽烟与否、血压水平等) 作为输入参量，输入到径向基函数网络(RBF 网络) 进行训练和识别.2.4 数学形态学数学形态学可以用来简化图像数据，保持图像的基本形态特征，同时去掉图像终于研究无关的部分，其基本运算有四种：膨胀、腐蚀，开运算和闭运算，通常是用在二维图形的处理，但近年来，国内外许多学者将其应用于一维信号的处理中。
作为一维信号处理的工具，数学形态学方法已经用于心电信号处理。
HenryChu 将数学形态学用于心电信号的噪声滤除[15]。
Trahanias 将其用于 QRS 波群的检测[16]。
郭兴明，姚晓帅，陈剑等[17][18]将数学形态学用于心音包络的提取与识别研究，他们首先利用数学形态学滤波和全波整流对原始心音进行预处理；然后利用形态学闭运算提取心音包络；进而应用形态学开运算来消除噪声包络；最后依据生理学知识对心音信号进行自动识别。
这种方法提取的信号包络毛刺较多且边缘不够光滑，但可以通过优化结构元素来改善；另外，包络信号与原始信号的外形差别比较大，但是心音的时域特征没有大的改变，因此可以用来对心音信号进行识别。
2.5 其它方法ADRicke等[19]采用隐马尔可夫模型来对心音信号进行自动分段，他们首先对原始心音信号进行带通滤波，然后进行规一化并在连续的0.04秒以内计算平均香农能量，其中每个时间段有0.02秒的重叠，接下来提取心音信号的光谱特性，进而通过特征提取计算一系列回归系数，最后用机器状态不可见的概率性状态机—隐马尔可夫模型来拟合心音信号，隐马尔可夫模型可以拟合心音图的四个心音状态：第一心音S1，收缩期，第二心音S2，舒张期。
D Gill[20]等则首先采用同态滤波器来获得平滑的心音图包络，这利于突出自动监测心音信号中人们感兴趣的部分，而从这些部分中提取的特征序列被用于隐马尔可夫模型的观察值。
研究表明，主要心音成分检测和识别的任务，可以在不使用任何同步采集信号的情况下，自动的完成。
王文辉等[21]提出阶梯值比较法并结合经验数据实现心音自定位，但这种方法受经验参数的制约，定位的准确性不高。
陈萌辉等[22]通过使用双门限、迭代等方法,改进了基于信号能量的分段算法,并首次引入短时过零率以更准确地定位分段边界，对心音信号进行自动识别。
王衍文等[23]提出一种基于Choi-Williams分布的心音信号检测方法，利用Choi—Williams分布分析心音信号的时频特性，这种方法不仅准确地判定了Sl和S2的起始位置，而且对心音信号的时域、频域和能量分布作了全面的描述。
周静等[24]，杨群清等[25]，童基均等[26]都将非线性时间序列分析的方法-----相关维应用于心音时间序列的分析，通过计算心音信号的关联维数，获取冠状动脉的特征信息，来区别正常心音，得到一种对冠状动脉疾病进行无损检测的方法。