基于小波变换心电信号特征检测算法.docVIP

基于小波变换心电信号特征检测算法 　　摘 要：本文提出了基于连续小波变换的心电信号QRS波识别算法。采用高斯函数的一阶导数作为小波基函数，其中QRS波中R波顶点范围的确定利用小波变换相应层次中的模极大值对的位置； R波顶点位置确定后，利用已有的平均心电周期，实现了P、T波提取的一种新方法―搜索窗口宽度自适应法，在这个新的窗口中利用原信号的差分值，找到P、T波的一些关键点，最终检测提取出完整的心电波形。 　　关键词：小波变换；QRS波识别；自适应窗口；心电波形 　　DOI：10.16640/j.cnki.37-1222/t.2016.17.104 　　1 引言 　　对心电信号进行检测分析时，首先要对提取到的心电信号进行去噪处理，包括去除基线漂移、去除50Hz工频干扰等。然后提取识别心电信号的各种波形，其主要任务是完成波形特征点参数值的确定，将模拟的心电信号进行数字化描述[1]。近些年来，对心电信号特征参数的提取涌现出了各种各样的算法，本文选择基于小波分析的方法对心电信号特征检测进行研究，该研究包括对QRS波群、P、T波等心电信号的典型波形特征和形态进行分析识别。 　　2 QRS波群检测研究 　　本文利用小波变换的方法检测信号的奇异点，其实质就是用小波变换对采集到的心电信号进行分解处理，得到不同频率成分的信号。而这个过程中首先要做的也是最重要的任务就是选择合适的小波函数。 　　（1）选择合适的小波函数 　　小波函数的选取方法多种多样，目前用于心电信号特征提取的小波函数有db5、haar、mexican-hat和二次样条[2]。本文在参考这些方法的基础上选用了具有一阶消失矩的小波函数，经过该小波函数的变换后，QRS波群正好对应相应的模极大值对，以此检测出QRS波。为了让模极值对的幅值更明显，小波函数应尽量选取平滑函数的一阶导数（即具有一阶消失矩），为满足上述要求本文选取高斯函数的一阶导数作为小波函数。 　　（2）R波峰值点的检测研究 　　利用具有一阶消失矩的小波函数分解变换心电信号后，得到一系列小波系数。在得到的小波系数中找到模极大值对的位置，同时找到原心电信号中峰值最高的那个点但必须保证这个点在模极大值对的位置区间内，将这个点作为R波顶点。若极值对中先出现极小值再出现极大值，则说明R波是正向的，若极值对先出现极大值再出现极小值，则说明R波是倒置的。确定好R波后，跳过200毫秒（心电信号不应期时间），继续分析检查后面的小波系数以此得到下一个R波。 　　（3）QRS波群起点检测研究 　　确定QRS波的起点和终点的常用方法有差分法，小波分析结合差分法[3]等，为了解决这些算法中存在的对干扰比较敏感，关键点判断时间过长等缺点，本文利用小波变换剔除干扰信号，除去大部分肌电干扰和基线漂移。通过采用小波分析的模极大值对判断方法获得Q波的顶点位置，然后根据心电信号的时序特征，在R波顶点位置的附近，找到各波形的区域范围，最后利用差分分析法准确得到QRS波波形的关键点。 　　（4）P波检测研究 　　P波在时间轴上出现在QRS波群中Q波的左侧，一般利用P波的宽度在心电周期（RR间期）的宽度中所占的比列，得到P波的搜索范围[4]。通常将R波顶点坐标值RP减去心电波形中PR段的最小宽度（即RP-0.06RT）作为搜索P波的起点，将R波顶点坐标值RP减去心电波形中PR段的长度（RP-0.205RT）作为搜索P波的终点，其中，RT为心电周期的平均值。确定完P波顶点后开始寻找P波起始点，具体思路是，从确定的P波顶点位置开始，向时间轴的左侧搜索大概一个P波宽度的区域，搜索过程中出现的第一个拐点就是P波起点，P波终点也用同样的方法得到，只是搜索方向为从P波顶点开始向时间轴右侧搜索。 　　（5）T波检测研究 　　首先要找到T波顶点，然后在T波顶点的基础上搜索T波的起点和终点。 　　T波搜索方法几乎与P波搜索方法相同。首先确定T波顶点位置的区间，主要依据T波信号的时宽性质和心电周期，在T波可能出现的大概范围内找到小波系数中的模极值对。然后在模极值对的相应范围内确定T波顶点位置，但要区分T波是倒置的还是正向的两种情况。最后根据心电的QT间期和QRS间期的时间长度，得到T波顶点范围大概为RP+0.1 RT ～RP+0.5 RT。查找到T波顶点后，在T波终点开始向S波方向搜索，找到的第一个拐点为T波起点，然后假设T波对称，找到T波终点。 　　3 小结 　　本文利用小波变换及小波对奇异点的检测特性，选择的小波函数为具有一阶消失矩的高斯函数的一阶导数。利用该函数对采集到的心电数据进行连续的小波分析处理，再根据原信号的差分值和小波变换得到的一系列小波函数，完成了基于连续小波变换分析的方法对QRS波群、P波和T波关键点的检测。 　　参考文献