[发明专利]一种心电信号R特征波检测方法

说明书

技术领域

本发明涉及心电信号自动检测与分析技术领域，特别涉及一种心电信号R特征波检测方法。

背景技术

心脏病具有隐蔽性和潜伏性，不发病的时候很难在心电图上表现出来，发病时又是短暂的，来不及观察心电图。为此需要给病人携带24小时Holter，进行24小时心电信号采集。但是这样就会造成数据量巨大，医生需要大量时间来一一查看心电图，寻找异常点，极大的增加了医生的负担。同时，诊断过程中会受到个人认知和情感的影响，使心脏病诊断具有主观性。应用电信号自动分析技术可以纠正这一偏差。现有的分析算法对于R波的分析算法相对完善，但对高频噪声和QRS波群形态改变上的处理效果不尽人意。

发明内容

本发明的目的是提供一种心电信号R特征波检测方法，以解决现有分析算法在应对高频噪声和QRS波群形态改变上的处理较差的问题。

本发明的目的是这样实现的：本发明所提供的心电信号R特征波检测方法，包括以下步骤：

a）信号采集：以250Hz的采集频率采集人体心电信号，并存储为TXT文档的数据形式，然后用Matlab软件将所述TXT文档存储的心电原始信号数据读取到电脑中。

b）对所述心电原始信号数据进行滤波处理：

b-1）对所述心电原始信号进行小波分解：选择DB6小波，对信号进行8层分解，得到各个尺度上的小波系数d_i；

b-2）采用改进的计算阈值方法，求取各尺度的阈值，对小波系数进行阈值化处理：

Ti=σi2log(n)ei-1,i=1,...,8]]>

其中，T_i为改进的阈值，i表示小波分解层数，e是自然常数，n表示采样点数，σ_i为小波系数绝对值的均值:

b-3）采用软阈值方法对信号进行阈值化处理：在不同尺度选取不同的阈值进行阈值化处理，得到滤波后的心电信号。

c）根据QRS波群的频率分布范围特征，选择第3、4尺度进行小波重构，得到重构后的心电信号S'；

选择进行小波重构的尺度时，其主要考虑QRS波群与P波、T波频率分布范围的不同，综合QRS波群频率分布所在尺度，且P波、T波频率分布较少的尺度进行小波重构。

d）对经过小波重构的心电信号进行能量窗变换，并选取极大值点：

d-1）能量窗变换：按下式，将经过小波重构的心电信号S'由时间域分析变换到能量域分析，得到心电信号能量曲线：