CN103156599A - 一种心电信号r特征波检测方法 - Google Patents
一种心电信号r特征波检测方法 Download PDFInfo
- Publication number
- CN103156599A CN103156599A CN201310116216XA CN201310116216A CN103156599A CN 103156599 A CN103156599 A CN 103156599A CN 201310116216X A CN201310116216X A CN 201310116216XA CN 201310116216 A CN201310116216 A CN 201310116216A CN 103156599 A CN103156599 A CN 103156599A
- Authority
- CN
- China
- Prior art keywords
- electrocardiosignal
- wavelet
- carried out
- point
- maximum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种心电信号R特征波检测方法,通过对采集的心电信号进行改进阈值的小波分解方法进行滤波,然后进行突出QRS波群的小波重构,作由时间域分析到能量域分析的能量窗变换,然后选取极大值点并进行优化,最后根据R波选取逻辑进行R波选取。本发明利用小波分析的多分辨率特性,将信号按频率分解到不同的尺度上,从而有针对性的对不同尺度上的信号进行处理,增加了算法的灵活性;采用能量窗变换的方法,将信号由时间域变换到能量域,有效抑制了高频噪声的干扰,提高了算法的稳固性;通过进行小波重构,将QRS波群提取出来,而将P波、T波当作噪声剔除,有效避免了高大的P波、T波在检测中造成的误检,提高了检测的精度。
Description
技术领域
本发明涉及心电信号自动检测与分析技术领域,特别涉及一种心电信号R特征波检测方法。
背景技术
心脏病具有隐蔽性和潜伏性,不发病的时候很难在心电图上表现出来,发病时又是短暂的,来不及观察心电图。为此需要给病人携带24小时Holter,进行24小时心电信号采集。但是这样就会造成数据量巨大,医生需要大量时间来一一查看心电图,寻找异常点,极大的增加了医生的负担。同时,诊断过程中会受到个人认知和情感的影响,使心脏病诊断具有主观性。应用电信号自动分析技术可以纠正这一偏差。现有的分析算法对于R波的分析算法相对完善,但对高频噪声和QRS波群形态改变上的处理效果不尽人意。
发明内容
本发明的目的是提供一种心电信号R特征波检测方法,以解决现有分析算法在应对高频噪声和QRS波群形态改变上的处理较差的问题。
本发明的目的是这样实现的:本发明所提供的心电信号R特征波检测方法,包括以下步骤:
a)信号采集:以250Hz的采集频率采集人体心电信号,并存储为TXT文档的数据形式,然后用Matlab软件将所述TXT文档存储的心电原始信号数据读取到电脑中。
b)对所述心电原始信号数据进行滤波处理:
b-1)对所述心电原始信号进行小波分解:选择DB6小波,对信号进行8层分解,得到各个尺度上的小波系数di;
b-2)采用改进的计算阈值方法,求取各尺度的阈值,对小波系数进行阈值化处理:
b-3)采用软阈值方法对信号进行阈值化处理:在不同尺度选取不同的阈值进行阈值化处理,得到滤波后的心电信号。
c)根据QRS波群的频率分布范围特征,选择第3、4尺度进行小波重构,得到重构后的心电信号S';
选择进行小波重构的尺度时,其主要考虑QRS波群与P波、T波频率分布范围的不同,综合QRS波群频率分布所在尺度,且P波、T波频率分布较少的尺度进行小波重构。
d)对经过小波重构的心电信号进行能量窗变换,并选取极大值点:
d-1)能量窗变换:按下式,将经过小波重构的心电信号S'由时间域分析变换到能量域分析,得到心电信号能量曲线:
其中,En表示第n个采样点的能量值;N为所选的窗口长度;M为总的采样点数;S'n表示所述小波重构后的心电信号S'的第n个数据;
所述窗口长度,是按如下方法确定的:根据采样频率的大小,选择最少的能将QRS波群覆盖的偶数个采样点的个数作为窗长;
d-2)选取极大值点:将所得到的心电信号能量曲线进行硬阈值化处理,即:
其中,Th为所选取的阈值,取Th=0.3*median(En),
然后选择经过硬阈值化处理后的心电信号能量曲线的波峰位置作为极大值点。
e)优化极大值点:设定2个时间阈值t1和t2,且t1<t2当任意两个极大值点的时间间隔小于t1时,就去掉这两个极大值点之间幅值较小的那个;当任意两个极大值点的时间间隔大于t2时,就在这两个极大值点之间寻找另一未被识别的极值点;当任意两个极大值点的时间间隔既大于等于t1,又小于等于t2,则该两个极大值点均保留,最终得到的经优化的极大值点,并且每一个所述经优化的极大值点对应一个QRS波群。
f)根据步骤e)中每个极大值点所在的时间点,在步骤b)中所述滤波后的心电信号上相应的时间点左右各7个采样点的范围内搜寻信号幅值最大的点,做为检测到的R波。
步骤b-3)中软阈值化处理的表达式为:
步骤c)中小波重构的表达式为:
步骤d-1)中所述窗口长度为26。
步骤e)中t1和t2是按如下方法确定的:
以Et作为步骤d-2)所得到的所有极大值点的时间间隔的平均值,规定t1=0.5×Et,t2=1.5×Et。
本发明对信号进行小波分解过程中选取特定的小波基函数及小波分解层数,同时在对信号进行阈值化处理过程中采用改进的阈值方法,使得滤波后的信号在滤除心电信号中夹杂的肌电干扰,基线漂移和工频干扰的同时,尽可能多的保留了有用的信息,改善了通用阈值过度平滑的现象。本发明通过进行小波重构,将QRS波群提取出来,而将P波、T波当作噪声剔除,有效避免了高大的P波、T波在检测中造成的误检,提高了检测的精度。另外,本发明采用了能量窗变换方法,将信号变换到能量域去分析从而解决了时域分析中,信号易受高频噪声的影响,且在滤波中不能被全部滤除的问题。在能量窗变换中,本发明充分考虑了心电信号的时域特征及QRS波群的跨越时间,进行窗长的选择。实验结果证明只有窗长为26时,噪声的伪波峰产生的多检及低幅值的QRS波造成的漏检现象才能最有效的避免。
本发明的方法解决了现有技术中在应对高频噪声和QRS波群形态改变方面所存在问题,可实现对心电信号R特征波快速、精确检测。
附图说明
图1为本发明的流程图。
图2为心电原始信号。
图3为进行滤波后的心电信号。
图4为进行小波重构后的心电信号。
图5为进行硬阈值化处理后的选取极大值点的心电信号。
图6为对极大值点进行优化的流程图。
图7为所检测到的R波。
具体实施方式
实施例1:
本实施例在Intel Pentium Dual E22002.20GHz,内存为3.00GB,操作***为Window XP的计算机中实现,整个R波检测算法采用Matlab语言编写。
本发明的实施过程如图1所示:
a)心电信号采集并输出:利用心电信号采集设备(北京迈创通元电子仪器有限公司的七导同步心电图模块(ECG CB),其采样频率为250HZ,心电信号的分析时间间隔长度设定为10s)采集人体的心电信号并以TXT的形式存储。在Matlab软件中,将TXT文档中的心电数据读取到环境中,如图2,得到心电原始信号。
b)利用改进阈值的小波分解方法对心电原始信号进行滤波处理:
b-1)选择Daubechies小波系列的DB6小波,进行8层分解,如表1所示:
表1:在250Hz的采样频率下对DB6小波进行8层分解
Scale | 250HZ |
1 | 62.5-125 |
2 | 31.25-62.5 |
3 | 15.625-31.25 |
4 | 7.8125-15.625 |
5 | 3.90625-7.8125 |
6 | 1.953125-3.90625 |
7 | 0.9765625-1.953125 |
8 | 0.48828125-0.9765625 |
然后提取各尺度上的小波系数di,进行阈值的计算,以得到经改进的阈值,即:
其中i表示小波分解层数,Ti为改进的阈值,e是自然常数,n表示采样点数,σi为小波系数绝对值的均值,其表达式为
相比现有的固定阈值方法和极大极小阈值法,按式(1)对阈值的算法进行改进之后,改进后的阈值具有边带自适应性,且保持了良好的去噪重构特性。
b-2)采用软阈值方法,在各尺度上选取相应经改进的阈值,按式(2)对心电信号进行阈值化处理,即:
其中j=i;
从而得到滤波后的心电信号,如图3所示。
通过采用改进阈值的小波分解方法进行滤波后,使得滤波后的心电信号,尽可能多的保留了有用的信息,改善了通用阈值过度平滑的现象,使滤波效果更加稳定。
c)进行突出QRS波群的小波重构:正常的心电信号QRS波群的频率分布范围是5-45Hz,从表1中可以看出,其主要集中于3、4尺度上,而P波和T波的频率分布范围为0.05到10Hz,在3、4尺度上没有或仅有少量分布,因此,根据QRS波群和P波、T波的频率分布范围的差别,选择QRS波群与P波、T波分布频率重叠最少的3、4尺度对经滤波后的心电信号进行小波重构,即:
经过小波重构以后,所得到的心电信号主要为QRS波群的信息,起到了凸显QRS波群的作用,如图4所示。
d)对小波重构后的心电信号进行由时间域分析到能量域分析的能量窗变换,并选取极大值点:
d-1)能量窗变换:按如下式(4),将经过小波重构的心电信号由时间域分析变换到能量域分析,得到心电信号能量曲线:
其中,N为所选的窗口长度(N=26),M为总的采样点数,S'n表示步骤c)小波重构后的心电信号S'的第n个数据。
能量窗变换中,窗长度的选取是一个关键,其直接决定R波检测算法是否有效。选择窗长时,根据采样频率的大小,选择最少的正好能将QRS波群覆盖的偶数个采样点的个数作为窗长。本实施例中,充分考虑了心电信号的时域特征及QRS波群的跨越时间,其N值的选择是按如下方法确定的:心电信号的采样频率为250Hz,而正常的QRS波群一般不超过0.1s,为25个采样点,我们选取窗长为偶数,为26。实验结果证明只有窗长为26时,噪声的伪波峰产生的多检及低幅值的QRS波造成的漏检现象才能最有效的避免。
时域分析中,信号易受高频噪声的影响,且在滤波中不能被全部滤除,针对此问题,采用能量床变换的方法将时间域分析变换到能量域分析。能量域相比于时间域分析,对噪声具有更好的鲁棒性。如图5,经能量窗变换以后,信号的QRS波群位置点变得更加突出,心拍与心拍之间的间隔更加明显,高频噪声的影响也相应的变弱。
d-2)通过设定阈值以选取极大值点:将得到的信号能量曲线进行硬阈值化处理:
其中,Th为所选取的阈值,取Th=0.3*median(En);
然后选取经过硬阈值化处理后的心电信号能量曲线的波峰位置作为极大值点,如图5所示。
e)优化极大值点:如图6所给出的流程图,设定2个时间阈值t1和t2,且t1<t2当任意两个极大值点的时间间隔小于t1时,就去掉这两个极大值点之间幅值较小的那个;当任意两个极大值点的时间间隔大于t2时,就在这两个极大值点之间寻找另一未被识别的极值点;如两个极大值点的时间间隔既大于等于t1,又小于等于t2,则该两个极大值点均保留,如此最终得到的经优化的每个极大值点都对应一个QRS波群。
图6中,Et表示步骤d-2)所得到的所有极大值点的时间间隔的平均值,t1=0.5×Et,t2=1.5×Et。
f)R波选取:根据步骤e)中所确定的每个极大值点所在的时间点,在步骤b)中滤波后心电信号上相应的时间点左右各7个采样点的范围内搜寻信号幅值最大的点,即为检测到的R波(图7)。
Claims (5)
1.一种心电信号R特征波检测方法,其特征是,包括以下步骤:
a)信号采集:以250Hz的采集频率采集人体心电信号,并存储为TXT文档的数据形式,然后用Matlab软件将所述TXT文档存储的心电原始信号数据读取到电脑中;
b)对所述心电原始信号数据进行滤波处理:
b-1)对所述心电原始信号进行小波分解:选择DB6小波,对信号进行8层分解,得到各个尺度上的小波系数di;
b-2)采用改进的计算阈值方法,求取各尺度的阈值,对小波系数进行阈值化处理:
b-3)采用软阈值方法对信号进行阈值化处理:在不同尺度选取不同的阈值进行阈值化处理,得到滤波后的心电信号;
c)根据QRS波群的频率分布范围特征,选择第3、4尺度进行小波重构,得到重构后的心电信号S';
d)对经过小波重构的心电信号进行能量窗变换,并选取极大值点:
d-1)能量窗变换:按下式,将经过小波重构的心电信号S'由时间域分析变换到能量域分析,得到心电信号能量曲线:
其中,En表示第n个采样点的能量值;N为所选的窗口长度;M为总的采样点数;S'n表示所述小波重构后的心电信号S'的第n个数据;
所述窗口长度,是按如下方法确定的:根据采样频率的大小,选择最少的能将QRS波群覆盖的偶数个采样点的个数作为窗长;
d-2)选取极大值点:将所得到的心电信号能量曲线进行硬阈值化处理,即:
其中,Th为所选取的阈值,取Th=0.3*median(En),
然后选择经过硬阈值化处理后的心电信号能量曲线的波峰位置作为极大值点;
e)优化极大值点:设定2个时间阈值t1和t2,且t1<t2,当任意两个极大值点的时间间隔小于t1时,就去掉这两个极大值点之间幅值较小的那个;当任意两个极大值点的时间间隔大于t2时,就在这两个极大值点之间寻找另一未被识别的极值点;当任意两个极大值点的时间间隔既大于等于t1,又小于等于t2,则该两个极大值点均保留,最终得到的经优化的极大值点,并且每一个所述经优化的极大值点对应一个QRS波群;
f)根据步骤e)中每个极大值点所在的时间点,在步骤b)中所述滤波后的心电信号上相应的时间点左右各7个采样点的范围内搜寻信号幅值最大的点,做为检测到的R波。
2.根据权利要求1所述的一种心电信号R特征波检测方法,其特征是,
步骤b-3)中软阈值化处理的表达式为:
4.根据权利要求1所述的一种心电信号R特征波检测方法,其特征是,步骤d-1)中所述窗口长度为26。
5.根据权利要求1所述的一种心电信号R特征波检测方法,其特征是,步骤e)中t1和t2是按如下方法确定的:
以Et作为步骤d-2)所得到的所有极大值点的时间间隔的平均值,规定t1=0.5×Et,t2=1.5×Et。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310116216.XA CN103156599B (zh) | 2013-04-03 | 2013-04-03 | 一种心电信号r特征波检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310116216.XA CN103156599B (zh) | 2013-04-03 | 2013-04-03 | 一种心电信号r特征波检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103156599A true CN103156599A (zh) | 2013-06-19 |
CN103156599B CN103156599B (zh) | 2014-10-15 |
Family
ID=48580499
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310116216.XA Active CN103156599B (zh) | 2013-04-03 | 2013-04-03 | 一种心电信号r特征波检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103156599B (zh) |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103565427A (zh) * | 2013-11-19 | 2014-02-12 | 深圳邦健生物医疗设备股份有限公司 | 准周期生理信号特征点的检测 |
CN103705234A (zh) * | 2013-12-05 | 2014-04-09 | 深圳先进技术研究院 | 动态心电信号数据中的检波方法和装置 |
CN104161510A (zh) * | 2014-08-27 | 2014-11-26 | 电子科技大学 | 一种多级导联心电信号qrs波形识别方法 |
CN104523266A (zh) * | 2015-01-07 | 2015-04-22 | 河北大学 | 一种心电信号自动分类方法 |
CN105748066A (zh) * | 2016-03-03 | 2016-07-13 | 深圳竹信科技有限公司 | 一种心电信号波形特征点的提取方法及装置 |
CN105899268A (zh) * | 2015-06-23 | 2016-08-24 | 中国科学院深圳先进技术研究院 | 基于gpu的并行心电信号分析方法 |
CN106562782A (zh) * | 2016-05-20 | 2017-04-19 | 彭慧敏 | 一种儿科护理专用心电监护仪 |
JP2017127398A (ja) * | 2016-01-18 | 2017-07-27 | 学校法人慶應義塾 | 情報処理装置、情報処理システム、情報処理方法及びプログラム |
CN107260163A (zh) * | 2017-07-31 | 2017-10-20 | 广东南方电信规划咨询设计院有限公司 | 无线远程心电监测***及方法 |
CN108158573A (zh) * | 2017-12-26 | 2018-06-15 | 智慧康源(厦门)科技有限公司 | 基于自适应阈值小波变换的心电信号降噪方法 |
CN108272451A (zh) * | 2018-02-11 | 2018-07-13 | 上海交通大学 | 一种基于改进小波变换的qrs波识别方法 |
CN108537100A (zh) * | 2017-11-17 | 2018-09-14 | 吉林大学 | 一种基于pca和lda分析的心电信号身份识别方法及*** |
CN108601543A (zh) * | 2016-11-03 | 2018-09-28 | 华为技术有限公司 | 一种ecg信号处理方法及装置 |
CN108697334A (zh) * | 2016-03-01 | 2018-10-23 | 圣犹达医疗用品心脏病学部门有限公司 | 用于标测心脏活动的方法和*** |
CN109117729A (zh) * | 2018-07-11 | 2019-01-01 | 上海夏先机电科技发展有限公司 | 心电图室性逸搏实时判断方法、装置、***及存储介质 |
CN110169767A (zh) * | 2019-07-08 | 2019-08-27 | 河北大学 | 一种心电信号的检索方法 |
CN110420022A (zh) * | 2019-07-29 | 2019-11-08 | 浙江大学 | 一种基于双密度小波变换的p波检测方法 |
CN110547786A (zh) * | 2018-06-01 | 2019-12-10 | 山东省计算中心(国家超级计算济南中心) | 基于多尺度小波变换的心电图波形的形态识别方法 |
CN111345815A (zh) * | 2020-02-11 | 2020-06-30 | 广州视源电子科技股份有限公司 | 一种心电信号中qrs波的检测方法、装置、设备和存储介质 |
CN111543981A (zh) * | 2020-03-16 | 2020-08-18 | 浙江好络维医疗技术有限公司 | 一种基于分段modwt及自适应阈值的动态心电图实时滤波方法 |
CN111920407A (zh) * | 2020-07-27 | 2020-11-13 | 广东省医疗器械研究所 | 基于小波变换的心电特征提取方法、***、装置及介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1951320A (zh) * | 2006-11-16 | 2007-04-25 | 上海交通大学 | 基于超完备特征的异常心电识别方法 |
US20080167567A1 (en) * | 2006-09-19 | 2008-07-10 | The Cleveland Clinic Foundation | Prediction and prevention of postoperative atrial fibrillation in cardiac surgery patients |
CN101828916A (zh) * | 2010-05-07 | 2010-09-15 | 深圳大学 | 心电信号处理*** |
US20110270111A1 (en) * | 2010-04-28 | 2011-11-03 | Maxime Cannesson | Method and apparatus for assessment of fluid responsiveness |
CN102626310A (zh) * | 2012-04-23 | 2012-08-08 | 天津工业大学 | 基于提升小波变换和改进近似包络的心电特征检测算法 |
-
2013
- 2013-04-03 CN CN201310116216.XA patent/CN103156599B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080167567A1 (en) * | 2006-09-19 | 2008-07-10 | The Cleveland Clinic Foundation | Prediction and prevention of postoperative atrial fibrillation in cardiac surgery patients |
CN1951320A (zh) * | 2006-11-16 | 2007-04-25 | 上海交通大学 | 基于超完备特征的异常心电识别方法 |
US20110270111A1 (en) * | 2010-04-28 | 2011-11-03 | Maxime Cannesson | Method and apparatus for assessment of fluid responsiveness |
CN101828916A (zh) * | 2010-05-07 | 2010-09-15 | 深圳大学 | 心电信号处理*** |
CN102626310A (zh) * | 2012-04-23 | 2012-08-08 | 天津工业大学 | 基于提升小波变换和改进近似包络的心电特征检测算法 |
Non-Patent Citations (3)
Title |
---|
ROMERO LEGARRETA,ETC: "R-wave Detection Using Continuous Wavelet Modulus Maxima", 《IEEE:BROWSE CONFERENCE PUBLICATIONS > COMPUTERS IN CARDIOLOGY》, 31 December 2003 (2003-12-31), pages 565 - 568 * |
Y.FERDI,ETC: "R wave detection using fractional digital differentiation", 《ITBM-RBM》, vol. 24, 31 December 2003 (2003-12-31), pages 273 - 280 * |
汪振兴: "心电信号特征提取和ST段识别算法研究", 《万方数据知识服务平台》, 30 November 2012 (2012-11-30) * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103565427A (zh) * | 2013-11-19 | 2014-02-12 | 深圳邦健生物医疗设备股份有限公司 | 准周期生理信号特征点的检测 |
CN103705234A (zh) * | 2013-12-05 | 2014-04-09 | 深圳先进技术研究院 | 动态心电信号数据中的检波方法和装置 |
CN103705234B (zh) * | 2013-12-05 | 2015-08-26 | 深圳先进技术研究院 | 动态心电信号数据中的检波方法和装置 |
CN104161510A (zh) * | 2014-08-27 | 2014-11-26 | 电子科技大学 | 一种多级导联心电信号qrs波形识别方法 |
CN104523266A (zh) * | 2015-01-07 | 2015-04-22 | 河北大学 | 一种心电信号自动分类方法 |
CN104523266B (zh) * | 2015-01-07 | 2017-04-05 | 河北大学 | 一种心电信号自动分类方法 |
CN105899268A (zh) * | 2015-06-23 | 2016-08-24 | 中国科学院深圳先进技术研究院 | 基于gpu的并行心电信号分析方法 |
US10258250B2 (en) | 2015-06-23 | 2019-04-16 | Shenzhen Institutes Of Advanced Technology Chinese Academy Of Sciences | GPU-based parallel electrocardiogram signal analysis method, computer readable storage medium and device |
CN105899268B (zh) * | 2015-06-23 | 2019-02-15 | 中国科学院深圳先进技术研究院 | 基于gpu的并行心电信号分析方法 |
JP2017127398A (ja) * | 2016-01-18 | 2017-07-27 | 学校法人慶應義塾 | 情報処理装置、情報処理システム、情報処理方法及びプログラム |
CN108697334A (zh) * | 2016-03-01 | 2018-10-23 | 圣犹达医疗用品心脏病学部门有限公司 | 用于标测心脏活动的方法和*** |
CN105748066A (zh) * | 2016-03-03 | 2016-07-13 | 深圳竹信科技有限公司 | 一种心电信号波形特征点的提取方法及装置 |
WO2017148452A1 (zh) * | 2016-03-03 | 2017-09-08 | 深圳竹信科技有限公司 | 一种心电信号波形特征点的提取方法及装置 |
CN105748066B (zh) * | 2016-03-03 | 2018-03-16 | 深圳竹信科技有限公司 | 一种心电信号波形特征点的提取方法及装置 |
CN106562782A (zh) * | 2016-05-20 | 2017-04-19 | 彭慧敏 | 一种儿科护理专用心电监护仪 |
CN108601543A (zh) * | 2016-11-03 | 2018-09-28 | 华为技术有限公司 | 一种ecg信号处理方法及装置 |
US11337638B2 (en) | 2016-11-03 | 2022-05-24 | Huawei Technologies Co., Ltd. | ECG signal processing method and apparatus |
CN107260163A (zh) * | 2017-07-31 | 2017-10-20 | 广东南方电信规划咨询设计院有限公司 | 无线远程心电监测***及方法 |
CN108537100A (zh) * | 2017-11-17 | 2018-09-14 | 吉林大学 | 一种基于pca和lda分析的心电信号身份识别方法及*** |
CN108158573A (zh) * | 2017-12-26 | 2018-06-15 | 智慧康源(厦门)科技有限公司 | 基于自适应阈值小波变换的心电信号降噪方法 |
CN108272451B (zh) * | 2018-02-11 | 2021-01-22 | 上海交通大学 | 一种基于改进小波变换的qrs波识别方法 |
CN108272451A (zh) * | 2018-02-11 | 2018-07-13 | 上海交通大学 | 一种基于改进小波变换的qrs波识别方法 |
CN110547786A (zh) * | 2018-06-01 | 2019-12-10 | 山东省计算中心(国家超级计算济南中心) | 基于多尺度小波变换的心电图波形的形态识别方法 |
CN110547786B (zh) * | 2018-06-01 | 2020-05-05 | 山东省计算中心(国家超级计算济南中心) | 基于多尺度小波变换的心电图波形的形态识别方法 |
CN109117729A (zh) * | 2018-07-11 | 2019-01-01 | 上海夏先机电科技发展有限公司 | 心电图室性逸搏实时判断方法、装置、***及存储介质 |
CN110169767A (zh) * | 2019-07-08 | 2019-08-27 | 河北大学 | 一种心电信号的检索方法 |
CN110169767B (zh) * | 2019-07-08 | 2021-09-21 | 河北大学 | 一种心电信号的检索方法 |
CN110420022A (zh) * | 2019-07-29 | 2019-11-08 | 浙江大学 | 一种基于双密度小波变换的p波检测方法 |
CN111345815A (zh) * | 2020-02-11 | 2020-06-30 | 广州视源电子科技股份有限公司 | 一种心电信号中qrs波的检测方法、装置、设备和存储介质 |
CN111543981A (zh) * | 2020-03-16 | 2020-08-18 | 浙江好络维医疗技术有限公司 | 一种基于分段modwt及自适应阈值的动态心电图实时滤波方法 |
CN111920407A (zh) * | 2020-07-27 | 2020-11-13 | 广东省医疗器械研究所 | 基于小波变换的心电特征提取方法、***、装置及介质 |
CN111920407B (zh) * | 2020-07-27 | 2023-09-22 | 广东省医疗器械研究所 | 基于小波变换的心电特征提取方法、***、装置及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN103156599B (zh) | 2014-10-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103156599A (zh) | 一种心电信号r特征波检测方法 | |
CN102697493B (zh) | 一种快速的脑电信号中眼电伪迹自动识别和去除的方法 | |
CN104523266B (zh) | 一种心电信号自动分类方法 | |
Seena et al. | A review on feature extraction and denoising of ECG signal using wavelet transform | |
CN103405227B (zh) | 基于双层形态学滤波的心电信号预处理方法 | |
CN102626310A (zh) | 基于提升小波变换和改进近似包络的心电特征检测算法 | |
CN105147248A (zh) | 基于生理信息的抑郁症评估***及其评估方法 | |
CN104473631B (zh) | 一种基于非负盲分离胎儿心电瞬时心率识别方法及*** | |
CN204931634U (zh) | 基于生理信息的抑郁症评估*** | |
CN102258368B (zh) | 胎心电检测的时域稀疏性线性混叠盲分离模型的判别方法 | |
CN105997043B (zh) | 一种基于腕式可穿戴设备的脉率提取方法 | |
CN103610461A (zh) | 基于双密度小波邻域相关阈值处理的脑电信号消噪方法 | |
CN105054926A (zh) | 心电信号特征信息的提取方法及装置 | |
CN106333658A (zh) | 一种光电容积脉搏波检测仪及检测方法 | |
Stevenson et al. | A nonparametric feature for neonatal EEG seizure detection based on a representation of pseudo-periodicity | |
CN103750835A (zh) | 一种心电信号特征检测算法 | |
CN107361753A (zh) | 基于脉搏波形特征点的人体健康状态监护方法 | |
CN105266800A (zh) | 一种基于低信噪比条件下胎儿心电盲分离方法 | |
CN106137185A (zh) | 一种基于模板构造小波的癫痫特征波检测方法 | |
CN101919704A (zh) | 一种心音信号定位、分段方法 | |
CN103190901A (zh) | 基于极值域均值模式分解和改进Hilbert包络的R波检测算法 | |
CN102240208A (zh) | 一种可在单片机中实现的心电去噪小波算法 | |
Elbuni et al. | ECG parameter extraction algorithm using (DWTAE) algorithm | |
Wu et al. | An improved method for ECG signal feature point detection based on wavelet transform | |
CN116584960A (zh) | 一种膈肌肌电信号降噪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |