CN113180680A - 一种改进的基于奇异谱分析的心电信号降噪方法 - Google Patents
一种改进的基于奇异谱分析的心电信号降噪方法 Download PDFInfo
- Publication number
- CN113180680A CN113180680A CN202110531843.4A CN202110531843A CN113180680A CN 113180680 A CN113180680 A CN 113180680A CN 202110531843 A CN202110531843 A CN 202110531843A CN 113180680 A CN113180680 A CN 113180680A
- Authority
- CN
- China
- Prior art keywords
- component
- reconstruction
- signal
- electrocardiosignal
- components
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Power Engineering (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种改进的基于奇异谱分析的心电信号降噪方法。具体流程包括:对心电信号选择合适的窗口长度构造轨迹矩阵并进行奇异值分解;计算每一个奇异值所对应的信号重构分量并对第一个重构分量采用低阶多项式拟合的方法去除趋势项;采用陷波器去除各重构分量中的工频干扰;通过计算重构分量间的相关系数的绝对值并与阈值比较选择出有效的信号分量;最后将有效信号分量叠加便可得到降噪后的心电信号。本发明方法由数据驱动,无需参考输入,适用于高斯白噪声、基线漂移、肌电干扰、电力线等噪声的去除,可提高后续的信号分析效果。本方法不仅信号失真程度小,而且在较低信噪比情况下也能较好地降噪,可推广至单路或多路其他生理信号的降噪。
Description
技术领域
本发明涉及心电信号的降噪处理,具体涉及一种改进的基于奇异谱分析的心电信号降噪方法。
背景技术
体表心电图是通过体表电极和采集设备记录到的标准化导联上的心脏电活动的波形,包含着丰富的生理和病理信息,被认为是人体的重要生理信号。心电图能反映受试者的心脏各部位的生理状况,是评价心脏功能和预测心脏疾病的重要依据之一。
心电信号本身非常微弱,通常情况下最大幅值不超过5mV,在采集过程中极易受到各种干扰,如人体呼吸运动引起的基线漂移、肌肉收缩与颤动引起的肌电干扰、环境噪声、电力线噪声等等。这些干扰和噪声的存在使得信号的波形发生变化,影响研究者对于波形的形态、心跳间隔等特征的获取,导致误诊或者使得信号失去临床价值。因此对心电信号进行降噪的预处理极为重要,有利于为后续的进一步信号分析提供基础,从而降低仪器的误警率。
近年来,研究人员一直致力于开发鲁棒性高的降噪技术,常用的降噪方法包括数字滤波法、自适应滤波、基于小波变换的降噪方法等,都有一定的局限性。如,数字滤波方法无法直接去除频谱与心电信号有重叠的噪声,而且可能改变心电信号的局部特征波甚至引起信号失真;自适应滤波通常需要一个高质量的参考输入信号,这种参考信号通常难以获取;小波基函数的选择对于基于小波变换的降噪方法有很大的影响。
奇异谱分析是一种独立于信号模型由数据驱动的时间序列分解方法,它将时间序列分解为趋势、振荡和噪声分量,它们本质上独立且可加。基于奇异谱分析的方法已在气候学、生物医学图像处理等方面得到有效应用。但在心电信号降噪领域的应用还有待更深入的探索,如:对于基线漂移、工频干扰等噪声难以去除;对于如何选择用于叠加的重构分量的个数也是一个难题。
发明内容
本发明的目的在于提出了一种改进的基于奇异谱分析的心电信号降噪方法。本发明方法能够对与心电信号频谱有重叠的噪声进行去除,且能高效地去除多种噪声包括基线漂移、电力线干扰、肌电干扰等。本发明方法不依赖于函数或者模型,由数据驱动,无需参考输入。该方法不仅不会改变局部特征波,使得信号失真程度小,对有效信息损害小,而且在较低信噪比情况下也能较好地实现降噪。
本发明提出了一种改进的基于奇异谱分析的心电信号降噪方法,具体步骤如下:
(1)读取一段采样频率为fs、长度为N的单导联心电信号数据,得到一组一维心电信号片段;
(2)对于步骤(1)得到的一维心电信号片段X=[x1,x2,…,xN],设定窗口长度为L(L≤N/2),将数据长度为N的原始单导联心电信号映射成K个L滞后向量,并用这些滞后向量形成一个轨迹矩阵,该轨迹矩阵表示为XTraj:
其中,K=N–L+1,且K>>L;
对轨迹矩阵XTraj做奇异值分解,将得到的L个奇异值按大小降序排列,即:
λ1≥λ2≥…≥λL
轨迹矩阵XTraj可以写成上述L个奇异值λ1,λ2,…,λL对应的分量矩阵Zi相加的形式,i为1,2至L:
(3)对于步骤(2)得到的分量矩阵Zi,对它们所在逆对角线上的数值累加求平均,得到信号重构分量Ci(i为1,2至L),可以表示为:
(4)仅针对步骤(3)得到的信号重构分量Ci中的第一个重构分量C1,用低阶多项式拟合C1得到f_C1,再用C1减去f_C1,得到去除非线性趋势后的新的第一个重构分量C1’;
(5)对步骤(4)中得到的新的第一个重构分量C1’和信号重构分量Ci(i为2,3至L),采用截止频率为50Hz的陷波滤波器去除所有分量中的工频干扰,得到滤波后的重构分量FCi(i为1,2至L);
(6)对经过步骤(5)得到的滤波后的重构分量FCi,计算它们中任意两个滤波后的重构分量FCp和FCq(p,q均为1,2至L)之间的相关系数ρ,得到大小为L×L的相关系数矩阵;其中,相关系数ρ的计算公式为:
设循环次数为w(初始值为1,且w≤L),分别将第w个滤波后的重构分量FCw与第一个滤波后的重构分量FC1、第二个滤波后的重构分量FC2之间的相关系数的绝对值记为CORR1w和CORR2w;然后设定经验阈值为t(0<t<1),若相关系数CORR1w>t或CORR2w>t,则w=w+1,并继续该循环,否则循环结束;当循环结束时,设置用于叠加的重构分量个数r为r=w–1;
(7)将步骤(6)选择的前r个滤波后的重构分量(即前r个奇异值λ1,λ2,…,λr所对应的重构分量FC1,FC2,…,FCr)叠加,得到降噪后的心电信号片段Y=[y1,y2,…,yN],即:
本发明中,步骤(4)仅针对步骤(3)得到的信号重构分量Ci中的第一个重构分量C1,采用减去低阶多项式拟合f_C1的方法得到去除非线性趋势后的新的第一个重构分量C1’,该步骤可去除信号中的基线漂移噪声。
本发明中,步骤(5)中对步骤(4)中得到的新的第一个重构分量C1’和信号重构分量Ci(i为2,3至L),进行截止频率为50Hz的陷波滤波,该步骤可去除信号中的工频干扰。
本发明中,步骤(6)中使用相关系数阈值法来确定要选择的重构分量个数:计算任意两个滤波后的重构分量FCp和FCq(p,q均为1,2至L)之间的相关系数矩阵,并计算所有滤波后的重构分量FCw分别与第一个滤波后的重构分量FC1、第二个滤波后的重构分量FC2之间的相关系数的绝对值,将这两个值与设定的阈值进行循环比较(循环次数w的初始值为1,且w≤L);当且仅当两者都小于等于阈值时,循环结束;将循环结束时的w值减去1,作为后续用于叠加的重构分量的个数。
本发明具有以下有益效果:
1.本发明方法仅由数据驱动,无需依赖模型函数或者其他参考信号,适用于高斯白噪声、基线漂移、肌电干扰、电力线等噪声的去除,可提高后续的信号分析效果;
2.本发明方法通过多项式拟合的方法去除第一个重构分量C1中的非线性趋势,使得原始心电信号序列中的基线漂移得到有效去除,最后得到的降噪后信号序列中不含或只含非常少量的基线漂移噪声;
3.本发明方法通过让所有重构分量Ci都经过一个截止频率为50Hz的陷波滤波器,尽可能地去除原始信号中的工频干扰,将奇异谱分析方法与滤波方法有效结合,有利于噪声的进一步去除;
4.本发明方法通过滤波后的重构分量之间的相关系数阈值法来实现滤波后的重构分量个数的选择,在低信噪比情况下也能较好地实现分量个数的选择,对提高降噪算法的鲁棒性有一定意义;
5.本发明方法不仅信号失真程度小,而且在较低信噪比情况下也能较好地降噪,可推广至单路或多路其他生理信号的降噪。
附图说明
图1是对心电信号(ECG)应用本发明方法进行降噪的总体流程图。
图2是本实施例中原始干净信号、添加10dB高斯白噪声后的带噪信号、降噪后信号的对比图。其中(a)为10s原始干净信号片段的示意图,(b)为干净信号片段添加10dB高斯白噪声后得到的带噪信号的示意图,(c)为采用本发明方法对带噪信号进行降噪后的信号示意图。
图3是本实施例中带噪信号、降噪后信号与噪声的对比图。其中,(a)为带噪信号与降噪后信号的对比示意图,(b)为本实施例残差(去除掉的噪声)的示意图。
具体实施方式
下面结合附图和实施例对本发明方法及其应用做进一步详细说明。实施例仅用于对本发明进行说明,而非限定条件。在本发明技术方案的基础上,凡根据本发明原理对该实施例的各种变形或修正,均不应排除在本发明的保护范围之外。
实施例1:将本发明的一种改进的基于奇异谱分析的心电信号降噪方法应用于添加10dB高斯白噪声的心电信号(ECG),其流程如图1所示;本实施例的数据来自公开的数据库MIT-BIH Arrhythmia Database;工作流程如下:
(1)选择MIT-BIH Arrhythmia Database数据库中标号为100的记录中的某10s片段,其采样率为360Hz,使用滤波等方法对其进行预处理,使原始信号尽可能地干净;得到的干净ECG信号如图2(a)所示;对得到的原始干净ECG信号添加10dB的高斯白噪声,添加噪声后得到的带噪信号如图2(b)所示;
(2)对步骤(1)中得到的10s带噪信号片段,设定窗口长度L为90,将数据长度为3600的原始心电信号映射成3511个90-滞后向量,并用这些滞后向量形成一个轨迹矩阵,表示为XTraj:
对轨迹矩阵XTraj做奇异值分解,将得到的90个奇异值按大小降序排列,即:
λ1≥λ2≥…≥λ90
轨迹矩阵XTraj可以写成上述90个奇异值λ1,λ2,…,λ90对应的分量矩阵Zi相加的形式,i为1,2至90:
(3)对于步骤(2)得到的分量矩阵Zi,对它们所在逆对角线上的数值累加求平均,得到信号重构分量Ci(i为1,2至90),可以表示为:
(4)仅针对对步骤(3)得到的重构分量Ci中的第一个重构分量C1,用低阶多项式拟合C1得到f_C1,再用C1减去f_C1,得到去除非线性趋势后的新的重构分量C1’;
(5)对步骤(4)中得到的新的第一个重构分量C1’和信号重构分量Ci(i为2,3至90),采用截止频率为50Hz的陷波滤波器去除所有分量中的工频干扰,得到滤波后的重构分量FCi(i为1,2至90);
(6)经过步骤(5)得到的滤波后的重构分量FCi,计算它们中任意两个滤波后的重构分量FCp和FCq(p,q均为1,2至90)之间的相关系数,得到大小为90×90的相关系数矩阵;设循环次数为w(初始值为1,且w≤90),分别将第w个滤波后的重构分量FCw与第一个滤波后的重构分量FC1、第二个滤波后的重构分量FC2之间的相关系数的绝对值记为CORR1w和CORR2w;然后设定经验阈值为t=0.0015,若相关系数CORR1w>t或CORR2w>t,则w=w+1,并继续该循环,否则循环结束;当循环结束时,设置用于叠加的重构分量个数r为r=w–1;这里得到的r=14;
(7)将步骤(6)选择的前14个滤波后的重构分量叠加,得到降噪后的心电信号片段Y=[y1,y2,…,y3600],即:
降噪后心电信号如图2(c)所示;图2展示了原始干净信号、带噪信号、降噪后信号的对比图,可见噪声被明显地去除,且降噪后信号未有明显失真,有效信息保留完整;
表1展示了对原始干净心电信号分别添加10dB高斯白噪声、基线漂移、肌电干扰、电力线噪声后应用本发明方法进行降噪后得到的评价指标结果;
评价指标分别为:
(1)降噪后信噪比SNROUT:
(2)均方误差MSE:
(3)相关系数CORR:
表1
10dB噪声 | SNR<sub>OUT</sub> | MSE | CORR |
高斯白噪声 | 17.4148 | 2.3148e-04 | 0.9909 |
基线漂移 | 20.3121 | 1.1879e-04 | 0.9954 |
肌电干扰 | 13.8507 | 5.2591e-04 | 0.9804 |
电力线噪声 | 33.0667 | 6.2996e-06 | 0.9998 |
Claims (4)
1.一种改进的基于奇异谱分析的心电信号降噪方法,其特征在于具体步骤如下:
(1)读取一段采样频率为fs、长度为N的单导联心电信号数据,得到一组一维心电信号片段;
(2)对于步骤(1)得到的一维心电信号片段X=[x1,x2,…,xN],设定窗口长度为L(L≤N/2),将数据长度为N的原始单导联心电信号映射成K个L滞后向量,并用这些滞后向量形成一个轨迹矩阵,该轨迹矩阵表示为XTraj:
其中,K=N–L+1,且K>>L;
对轨迹矩阵XTraj做奇异值分解,将得到的L个奇异值按大小降序排列,即:
λ1≥λ2≥…≥λL
轨迹矩阵XTraj写成上述L个奇异值λ1,λ2,…,λL对应的分量矩阵Zi相加的形式,其中i为1,2至L:
(3)对于步骤(2)得到的分量矩阵Zi,对它们所在逆对角线上的数值累加求平均,得到信号重构分量Ci,其中i为1,2至L,表示为:
(4)仅针对步骤(3)得到的信号重构分量Ci中的第一个重构分量C1,用低阶多项式拟合C1得到f_C1,再用C1减去f_C1,得到去除非线性趋势后的新的第一个重构分量C1’;
(5)对步骤(4)中得到的新的第一个重构分量C1’和信号重构分量Ci(i为2,3至L),采用截止频率为50Hz的陷波滤波器去除所有分量中的工频干扰,得到滤波后的重构分量FCi(i为1,2至L);
(6)对经过步骤(5)得到的滤波后的重构分量FCi,计算它们中任意两个滤波后的重构分量FCp和FCq(p,q均为1,2至L)之间的相关系数ρ,得到大小为L×L的相关系数矩阵;其中,相关系数ρ的计算公式为:
设循环次数为w(初始值为1,且w≤L),分别将第w个滤波后的重构分量FCw与第一个滤波后的重构分量FC1、第二个滤波后的重构分量FC2之间的相关系数的绝对值记为CORR1w和CORR2w;然后设定经验阈值为t(0<t<1),若相关系数CORR1w>t或CORR2w>t,则w=w+1,并继续该循环,否则循环结束;当循环结束时,设置用于叠加的重构分量个数r为r=w–1;
(7)将步骤(6)选择的前r个滤波后的重构分量(即前r个奇异值λ1,λ2,…,λr所对应的重构分量FC1,FC2,…,FCr)叠加,得到降噪后的心电信号片段Y=[y1,y2,…,yN],即:
2.根据权利要求1所述的方法,其特征在于步骤(4)对步骤(3)得到的信号重构分量Ci中的第一个重构分量C1,采用减去低阶多项式拟合f_C1的方法得到去除非线性趋势后的新的重构分量C1’,该步骤可去除信号中的基线漂移噪声。
3.根据权利要求1所述的方法,其特征在于步骤(5)中对步骤(4)中得到的新的第一个重构分量C1’和信号重构分量Ci(i为2,3至L),进行截止频率为50Hz的陷波滤波,该步骤可去除信号中的工频干扰。
4.根据权利要求1所述的方法,其特征在于步骤(6)中使用相关系数阈值法来确定要选择的重构分量个数:计算任意两个滤波后的重构分量FCp和FCq(p,q均为1,2至L)之间的相关系数矩阵,并计算所有滤波后的重构分量FCw分别与第一个滤波后的重构分量FC1、第二个滤波后的重构分量FC2之间的相关系数的绝对值,将这两个值与设定的阈值进行循环比较,循环次数w的初始值为1,且w≤L;当且仅当两者都小于等于阈值时,循环结束;将循环结束时的w值减去1,作为后续用于叠加的重构分量的个数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110531843.4A CN113180680B (zh) | 2021-05-17 | 2021-05-17 | 一种改进的基于奇异谱分析的心电信号降噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110531843.4A CN113180680B (zh) | 2021-05-17 | 2021-05-17 | 一种改进的基于奇异谱分析的心电信号降噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113180680A true CN113180680A (zh) | 2021-07-30 |
CN113180680B CN113180680B (zh) | 2022-06-21 |
Family
ID=76982072
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110531843.4A Active CN113180680B (zh) | 2021-05-17 | 2021-05-17 | 一种改进的基于奇异谱分析的心电信号降噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113180680B (zh) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5175710A (en) * | 1990-12-14 | 1992-12-29 | Hutson William H | Multi-dimensional data processing and display |
US20070203418A1 (en) * | 2006-02-27 | 2007-08-30 | Cardiosoft L.L.P. | Multi-Channel System for Beat to Beat QT Interval Variability |
CN103190898A (zh) * | 2013-04-23 | 2013-07-10 | 何怡刚 | 心磁信号噪声自适应滤波消除设计方法 |
CN103431856A (zh) * | 2013-08-30 | 2013-12-11 | 深圳市理邦精密仪器股份有限公司 | 在多导联同步心电信号中选择心电导联的方法及装置 |
EP2686696A1 (en) * | 2011-03-17 | 2014-01-22 | Koninklijke Philips N.V. | Restriction of the imaging region for mri in an inhomogeneous magnetic field |
CN105266800A (zh) * | 2015-12-02 | 2016-01-27 | 广东工业大学 | 一种基于低信噪比条件下胎儿心电盲分离方法 |
US9289150B1 (en) * | 2012-08-17 | 2016-03-22 | Analytics For Life | Non-invasive method and system for characterizing cardiovascular systems |
CN109077715A (zh) * | 2018-09-03 | 2018-12-25 | 北京工业大学 | 一种基于单导联的心电信号自动分类方法 |
CN109165556A (zh) * | 2018-07-24 | 2019-01-08 | 吉林大学 | 一种基于grnn身份识别方法 |
CN109157220A (zh) * | 2018-09-11 | 2019-01-08 | 上海宏桐实业有限公司 | 一种基于多通道动态监测的呼吸指数提取系统及工作方法 |
-
2021
- 2021-05-17 CN CN202110531843.4A patent/CN113180680B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5175710A (en) * | 1990-12-14 | 1992-12-29 | Hutson William H | Multi-dimensional data processing and display |
US20070203418A1 (en) * | 2006-02-27 | 2007-08-30 | Cardiosoft L.L.P. | Multi-Channel System for Beat to Beat QT Interval Variability |
EP2686696A1 (en) * | 2011-03-17 | 2014-01-22 | Koninklijke Philips N.V. | Restriction of the imaging region for mri in an inhomogeneous magnetic field |
US9289150B1 (en) * | 2012-08-17 | 2016-03-22 | Analytics For Life | Non-invasive method and system for characterizing cardiovascular systems |
CN103190898A (zh) * | 2013-04-23 | 2013-07-10 | 何怡刚 | 心磁信号噪声自适应滤波消除设计方法 |
CN103431856A (zh) * | 2013-08-30 | 2013-12-11 | 深圳市理邦精密仪器股份有限公司 | 在多导联同步心电信号中选择心电导联的方法及装置 |
CN105266800A (zh) * | 2015-12-02 | 2016-01-27 | 广东工业大学 | 一种基于低信噪比条件下胎儿心电盲分离方法 |
CN109165556A (zh) * | 2018-07-24 | 2019-01-08 | 吉林大学 | 一种基于grnn身份识别方法 |
CN109077715A (zh) * | 2018-09-03 | 2018-12-25 | 北京工业大学 | 一种基于单导联的心电信号自动分类方法 |
CN109157220A (zh) * | 2018-09-11 | 2019-01-08 | 上海宏桐实业有限公司 | 一种基于多通道动态监测的呼吸指数提取系统及工作方法 |
Non-Patent Citations (3)
Title |
---|
JORGE BARRIOS-MURIEL: "A simple SSA-based de-noising technique to remove ECG interferencein EMG signals", 《BIOMEDICAL SIGNAL PROCESSING AND CONTROL》 * |
KAIYUE HE: "Recognition of Premature Ventricular Contraction Beat from 12-", 《IEEE》 * |
赵烨南: "心律失常的多重分形去趋势波动分析", 《北京生物医学工程》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113180680B (zh) | 2022-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jain et al. | Riemann Liouvelle fractional integral based empirical mode decomposition for ECG denoising | |
CN108158573B (zh) | 基于自适应阈值小波变换的心电信号降噪方法 | |
Sayadi et al. | Multiadaptive bionic wavelet transform: Application to ECG denoising and baseline wandering reduction | |
Smital et al. | Adaptive wavelet wiener filtering of ECG signals | |
Bahaz et al. | Efficient algorithm for baseline wander and powerline noise removal from ECG signals based on discrete Fourier series | |
Bhateja et al. | A novel approach for suppression of powerline interference and impulse noise in ECG signals | |
Rahman et al. | Noise cancellation in ECG signals using computationally simplified adaptive filtering techniques: Application to biotelemetry | |
Butt et al. | Denoising practices for electrocardiographic (ECG) signals: a survey | |
Malghan et al. | A review on ECG filtering techniques for rhythm analysis | |
Narwaria et al. | Removal of baseline wander and power line interference from ECG signal-a survey approach | |
Galiana-Merino et al. | Power line interference filtering on surface electromyography based on the stationary wavelet packet transform | |
Lu et al. | Model-based ECG denoising using empirical mode decomposition | |
Tang et al. | ECG de-noising based on empirical mode decomposition | |
Almalchy et al. | Noise removal from ECG signal based on filtering techniques | |
Rashmi et al. | ECG denoising using wavelet transform and filters | |
CN110292374B (zh) | 基于奇异谱分析和变分模态分解的心电信号去基线漂移方法 | |
Lin et al. | Discrete-wavelet-transform-based noise reduction and R wave detection for ECG signals | |
Vega-Martinez et al. | ECG baseline drift removal using discrete wavelet transform | |
CN113180680B (zh) | 一种改进的基于奇异谱分析的心电信号降噪方法 | |
Dubey et al. | Two-stage nonlocal means denoising of ECG signals | |
CN113598785B (zh) | 基于小波域稀疏特性的心电信号去噪方法 | |
Suchetha et al. | Empirical mode decomposition-based subtraction techniques for 50 Hz interference reduction from electrocardiogram | |
Mishra et al. | ECG signal denoising using time-frequency based filtering approach | |
Muppalla et al. | Performance evaluation of different denoising techniques for physiological signals | |
Kaur et al. | Adaptive wavelet thresholding for noise reduction in electrocardiogram (ECG) signals |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |