CN111680785B - 基于稀疏特性与对抗神经网络相结合的ecg信号处理方法 - Google Patents
基于稀疏特性与对抗神经网络相结合的ecg信号处理方法 Download PDFInfo
- Publication number
- CN111680785B CN111680785B CN202010482112.0A CN202010482112A CN111680785B CN 111680785 B CN111680785 B CN 111680785B CN 202010482112 A CN202010482112 A CN 202010482112A CN 111680785 B CN111680785 B CN 111680785B
- Authority
- CN
- China
- Prior art keywords
- signal
- noise
- neural network
- ecg signal
- network
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- 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/7235—Details of waveform analysis
-
- 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/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- 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/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Biophysics (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Signal Processing (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Surgery (AREA)
- Mathematical Physics (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Data Mining & Analysis (AREA)
- Physiology (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Fuzzy Systems (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Complex Calculations (AREA)
Abstract
一种基于稀疏特性与对抗神经网络相结合的ECG信号处理方法,利用稀疏提取出ECG信号中深层次特征,对传统的对抗神经网络进行改进,利用对抗网络的特性不断学习优化,高精度实现了去除ECG信号中的噪声干扰的同时减少了运算时间。考虑到不同人体之间存在的个体差异特征,针对远程医院背景下的信号噪声多,干扰大等因素,利用大数据特征优势,引入深度学习,针对ECG信号的时域特性通过发生器和鉴别器之间的博弈,利用对抗性思维不断积累ECG信号噪声分布的知识,利用支持向量机(SVM)算法对去噪后的信号进行定性评估。
Description
技术领域
本发明涉及心电信号处理技术领域,具体涉及一种基于稀疏特性与对抗神经网络相结合的ECG信号处理方法。
背景技术
ECG信号包含六种不同类型的波形,分别是P、Q、R、S、T、U。在心脏状况良好的情况下,会出现有规律性的ECG信号曲线,但ECG信号采集过程中往往会受噪声的干扰,使得心电图质量不佳,因此能够快速准确的滤出噪声还能尽可能的保留信号信息的完整性对ECG信号的检测尤为重要。
目前ECG信号降噪方法大多是单靠某一种数字滤波器去除噪声,由于噪声种类繁多,其实用效果并不理想,有的降噪后波形严重失真,丢失了信号的医学特性。由于信号中部分噪声频域分布特征复杂,传统的方法很难滤出干净。在许多实际的去噪问题中,噪声信号包含大量的稀疏信息,这有助于去噪。然而,常用的去噪方法如低通滤波法和小波去噪法有许多局限性,如投影空间刚性或高频成分损失。随着深度学习的发展,基于深度学习的降噪方法虽然取得了一些成就,但也存在着一些问题。例如模糊神经网络和小波神经网络只能去除ECG信号中的一种噪声,不能反映噪声的复杂性和多样性。改进的小波神经网络和小波神经网络能够去除ECG信号中的三种经典噪声,但其结果会失去原始重要信息。
发明内容
本发明为了克服以上技术的不足,提供了一种利用稀疏提取出ECG信号中深层次特征,利用对抗网络的特性不断学习优化,高精度去除ECG信号中的噪声的ECG信号处理方法。
本发明克服其技术问题所采用的技术方案是:
一种基于稀疏特性与对抗神经网络相结合的ECG信号处理方法,包括:
a)从MIT-BIH噪声压力测试数据库中选择EM、BM和MA噪声记录作为噪声数据v;
b)深度学习的对抗神经网络中生成网络输入的信号为含有噪声数据v的信号y,将含有噪声数据v的信号y重新构建干净的原始信号y*实现信号y的降噪,得到降噪后的信号y′,将原始信号y*与降噪后的信号y′作为输入信号输入进对抗神经网络中的判别器,使用支持向量机来评价降噪后的ECG信号质量;
c)通过对抗神经网络的学习法将步骤b)中的生成网络模型进行训练,训练后的生成网络模型将心电信号进行高保真降噪;
d)通过公式lD=log(1-D(y*))+log(D(G(y)))计算对抗神经网络中的判别器的损失函数lD,式中G(y)为以含有噪声数据v的信号y作为输入信号并通过生成网络去噪后的信号,D(G(y)为判别网络判别信号是来自原始信号而不是生成网络去噪后的信号的概率;
e)当判别器无法判别输入的信号是原始信号还是生成网络生成的信号时,保存模型。
进一步的,步骤a)中从MIT-BIH噪声压力测试数据库中分别选择编号为103、105、111、116、122、205、213、219、223、230的噪声记录。
进一步的,步骤a)中通过将信噪比分别设置为15dB、20dB和25dB在计ECG信号中加入噪声数据v,每种噪声训练样本数量为54000个,测试样本数量为5940个。
进一步的,步骤b)包括如下步骤:
b-1)通过公式Y={y=Ax+v}计算不确定系数线性模型中采样的一组数据Y,其中A∈Rn×m,A为宽为n高为m的线性变换矩阵,R为矩阵,x为稀疏编码,x∈Rm,将x的位置是非零的概率定义为γ,x非零值服从均匀分布U(0,1),噪声数据v从正态分布N(0,SNR)中采样,
b-2)设置m=100,n=50,将γ设为0.1;
b-3)选择递归神经网络作为编码部分,编码结构方程为h=0,ht+1=f(Uht+Wy+b),式中ht为稀疏编码x的估计值,h为初始隐藏状态,ht+1是编码部分的第t+1个状态,Wy与b为学习参数,设置4层递归网络作为编码部分,隐藏层单元的数量与x的维数相同,隐藏层每一层都将前一层的输出作为输入;
b-4)通过公式x=hT,r(x)=f(Wrx+br)+x计算残差块r(x),式中hT是编码部分的第T个状态,Wr与br为学习参数;
b-5)通过公式y′=W0r(x)+b0计算降噪后的信号y′,式中W0与b0为学习参数。
进一步的,步骤c)包括如下步骤:
c-3)通过公式lG=lg(y,y*)+αldist+βldist-max计算生成器的损失函数lG,α=0.7,β=0.2,ldist-max为ldist的最大值。
本发明的有益效果是:本发明克服了传统去噪方法对投影空间的限制及高频分量损失的局限性,提出了一种基于稀疏特性与对抗神经网络相结合的ECG信号处理方法。利用稀疏提取出ECG信号中深层次特征,对传统的对抗神经网络进行改进,将生成网络设计为去噪模型,包括编码部分、去噪部分和线性恢复部分。利用对抗网络的特性不断学习优化,高精度实现了去除ECG信号中的噪声干扰的同时减少了运算时间。考虑到不同人体之间存在的个体差异特征,针对远程医院背景下的信号噪声多,干扰大等因素,对生成网络的损失函数进行创新。同时,为了保持数据的原始特征,利用鉴别器网络来帮助生成网络中的去噪模型学习,利用对抗性思维不断积累ECG信号噪声分布的知识。最后利用支持向量机(SVM)算法对去噪后的信号进行定性评估。
具体实施方式
下面对本发明做进一步说明。
一种基于稀疏特性与对抗神经网络相结合的ECG信号处理方法,包括:
a)从MIT-BIH噪声压力测试数据库中选择EM、BM和MA噪声记录作为噪声数据v;
b)深度学习的对抗神经网络中生成网络输入的信号为含有噪声数据v的信号y,将含有噪声数据v的信号y重新构建干净的原始信号y*实现信号y的降噪,得到降噪后的信号y′,将原始信号y*与降噪后的信号y′作为输入信号输入进对抗神经网络中的判别器,使用支持向量机来评价降噪后的ECG信号质量;
c)通过对抗神经网络的学习法将步骤b)中的生成网络模型进行训练,训练后的生成网络模型将心电信号进行高保真降噪;
d)通过公式lD=log(1-D(y*))+log(D(G(y)))计算对抗神经网络中的判别器的损失函数lD,式中G(y)为以含有噪声数据v的信号y作为输入信号并通过生成网络去噪后的信号,D(G(y)为判别网络判别信号是来自原始信号而不是生成网络去噪后的信号的概率;
e)当判别器无法判别输入的信号是原始信号还是生成网络生成的信号时,保存模型。
利用稀疏提取出ECG信号中深层次特征,对传统的对抗神经网络进行改进,利用对抗网络的特性不断学习优化,高精度实现了去除ECG信号中的噪声干扰的同时减少了运算时间。考虑到不同人体之间存在的个体差异特征,针对远程医院背景下的信号噪声多,干扰大等因素,利用大数据特征优势,引入深度学习,针对ECG信号的时域特性通过发生器和鉴别器之间的博弈,利用对抗性思维不断积累ECG信号噪声分布的知识,利用支持向量机(SVM)算法对去噪后的信号进行定性评估。
优选的,步骤a)中从MIT-BIH噪声压力测试数据库中分别选择编号为103、105、111、116、122、205、213、219、223、230的噪声记录。步骤a)中通过将信噪比分别设置为15dB、20dB和25dB在计ECG信号中加入噪声数据v,每种噪声训练样本数量为54000个,测试样本数量为5940个。
进一步的,步骤b)包括如下步骤:
b-1)通过公式Y={y=Ax+v}计算不确定系数线性模型中采样的一组数据Y,其中A∈Rn×m,A为宽为n高为m的线性变换矩阵,R为矩阵,x为稀疏编码,x∈Rm,将x的位置是非零的概率定义为γ,x非零值服从均匀分布U(0,1),噪声数据v从正态分布N(0,SNR)中采样,
b-2)设置m=100,n=50,将γ设为0.1;
b-3)选择递归神经网络作为编码部分,编码结构方程为h=0,ht+1=f(Uht+Wy+b),式中ht为稀疏编码x的估计值,h为初始隐藏状态,ht+1是编码部分的第t+1个状态,Wy与b为学习参数,设置4层递归网络作为编码部分,隐藏层单元的数量与x的维数相同,隐藏层每一层都将前一层的输出作为输入;
b-4)通过公式x=hT,r(x)=f(Wrx+br)+x计算残差块r(x),式中hT是编码部分的第T个状态,Wr与br为学习参数;
b-5)通过公式y′=W0r(x)+b0计算降噪后的信号y′,式中W0与b0为学习参数。
步骤c)包括如下步骤:
c-3)通过公式lG=lg(y,y*)+αldist+βldist-max计算生成器的损失函数lG,α=0.7,β=0.2,ldist-max为ldist的最大值。
Claims (4)
1.一种基于稀疏特性与对抗神经网络相结合的ECG信号处理方法,其特征在于,包括:
a)从MIT-BIH噪声压力测试数据库中选择EM、BM和MA噪声记录作为噪声数据v;
b)深度学习的对抗神经网络中生成网络输入的信号为含有噪声数据v的信号y,将含有噪声数据v的信号y重新构建干净的原始信号y*实现信号y的降噪,得到降噪后的信号y′,将原始信号y*与降噪后的信号y′作为输入信号输入进对抗神经网络中的判别器,使用支持向量机来评价降噪后的ECG信号质量;
c)通过对抗神经网络的学习法将步骤b)中的生成网络模型进行训练,训练后的生成网络模型将心电信号进行高保真降噪;
d)通过公式lD=log(1-D(y*))+log(D(G(y)))计算对抗神经网络中的判别器的损失函数lD,式中G(y)为以含有噪声数据v的信号y作为输入信号并通过生成网络去噪后的信号,D(G(y)为判别网络判别信号是来自原始信号而不是生成网络去噪后的信号的概率;
e)当判别器无法判别输入的信号是原始信号还是生成网络生成的信号时,保存模型;
步骤b)包括如下步骤:
b-1)通过公式Y={y=Ax+v}计算不确定系数线性模型中采样的一组数据Y,其中A∈Rn ×m,A为宽为n高为m的线性变换矩阵,R为矩阵,x为稀疏编码,x∈Rm,将x的位置是非零的概率定义为γ,x非零值服从均匀分布U(0,1),噪声数据v从正态分布N(0,SNR)中采样,
b-2)设置m=100,n=50,将γ设为0.1;
b-3)选择递归神经网络作为编码部分,编码结构方程为h=0,ht+1=f(Uht+Wy+b),式中ht为稀疏编码x的估计值,h为初始隐藏状态,ht+1是编码部分的第t+1个状态,Wy与b为学习参数,设置4层递归网络作为编码部分,隐藏层单元的数量与x的维数相同,隐藏层每一层都将前一层的输出作为输入;
b-4)通过公式x=hT,r(x)=f(Wrx+br)+x计算残差块r(x),式中hT是编码部分的第T个状态,Wr与br为学习参数;
b-5)通过公式y′=W0r(x)+b0计算降噪后的信号y′,式中W0与b0为学习参数。
2.根据权利要求1所述的基于稀疏特性与对抗神经网络相结合的ECG信号处理方法,其特征在于:步骤a)中从MIT-BIH噪声压力测试数据库中分别选择编号为103、105、111、116、122、205、213、219、223、230的噪声记录。
3.根据权利要求1所述的基于稀疏特性与对抗神经网络相结合的ECG信号处理方法,其特征在于:步骤a)中通过将信噪比分别设置为15dB、20dB和25dB在计ECG信号中加入噪声数据v,每种噪声训练样本数量为54000个,测试样本数量为5940个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010482112.0A CN111680785B (zh) | 2020-05-29 | 2020-05-29 | 基于稀疏特性与对抗神经网络相结合的ecg信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010482112.0A CN111680785B (zh) | 2020-05-29 | 2020-05-29 | 基于稀疏特性与对抗神经网络相结合的ecg信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111680785A CN111680785A (zh) | 2020-09-18 |
CN111680785B true CN111680785B (zh) | 2021-09-24 |
Family
ID=72453765
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010482112.0A Active CN111680785B (zh) | 2020-05-29 | 2020-05-29 | 基于稀疏特性与对抗神经网络相结合的ecg信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111680785B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112487914B (zh) * | 2020-11-25 | 2021-08-31 | 山东省人工智能研究院 | 一种基于深度卷积生成对抗网络的ecg降噪方法 |
CN112807000B (zh) * | 2021-02-04 | 2023-02-28 | 首都师范大学 | 鲁棒性脑电信号的生成方法及装置 |
CN113081002B (zh) * | 2021-03-31 | 2023-03-31 | 灵犀医学科技(北京)有限公司 | 脑电信号的伪差去除方法、装置及电子设备 |
CN113935378A (zh) * | 2021-10-14 | 2022-01-14 | 山东省人工智能研究院 | 一种基于对抗性深度全卷积网络的ecg降噪方法 |
CN114781445B (zh) * | 2022-04-11 | 2022-11-18 | 山东省人工智能研究院 | 一种基于可解释性的深度神经网络的心电信号降噪方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108113647A (zh) * | 2016-11-28 | 2018-06-05 | 深圳先进技术研究院 | 一种心电信号分类装置及方法 |
CN110223254A (zh) * | 2019-06-10 | 2019-09-10 | 大连民族大学 | 一种基于对抗生成网络的图像去噪方法 |
CN110840445A (zh) * | 2019-12-12 | 2020-02-28 | 山东省计算中心(国家超级计算济南中心) | 一种动态心电信号的自动降噪方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107832737B (zh) * | 2017-11-27 | 2021-02-05 | 上海优加利健康管理有限公司 | 基于人工智能的心电图干扰识别方法 |
CN110974217B (zh) * | 2020-01-03 | 2022-08-09 | 苏州大学 | 基于卷积自编码器的双阶段心电信号降噪方法 |
-
2020
- 2020-05-29 CN CN202010482112.0A patent/CN111680785B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108113647A (zh) * | 2016-11-28 | 2018-06-05 | 深圳先进技术研究院 | 一种心电信号分类装置及方法 |
CN110223254A (zh) * | 2019-06-10 | 2019-09-10 | 大连民族大学 | 一种基于对抗生成网络的图像去噪方法 |
CN110840445A (zh) * | 2019-12-12 | 2020-02-28 | 山东省计算中心(国家超级计算济南中心) | 一种动态心电信号的自动降噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111680785A (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111680785B (zh) | 基于稀疏特性与对抗神经网络相结合的ecg信号处理方法 | |
CN111814656B (zh) | 一种基于对抗生成网络的心电信号降噪方法 | |
CN112946749B (zh) | 基于数据增广训练深度神经网络压制地震多次波的方法 | |
CN108416740A (zh) | 一种消除椒盐噪声的迭代自适应中值滤波算法 | |
CN113723171B (zh) | 基于残差生成对抗网络的脑电信号去噪方法 | |
CN112884149B (zh) | 一种基于随机敏感度st-sm的深度神经网络剪枝方法及系统 | |
CN111650654B (zh) | 联合emd与wt算法的地面磁共振信号尖峰噪声剔除方法 | |
CN113642484B (zh) | 一种基于bp神经网络的大地电磁信号噪声压制方法及系统 | |
He et al. | Noise reduction of ECG signals through genetic optimized wavelet threshold filtering | |
CN115877461A (zh) | 一种基于多尺度注意力交互网络的沙漠地震噪声抑制方法 | |
CN113486750A (zh) | 一种基于改进vmd算法和小波包去噪算法的油浸式变压器局部放电信号去噪方法 | |
CN110706186B (zh) | 基于压缩感知及Renyi熵的水下信息抗干扰方法 | |
CN114648048A (zh) | 基于变分自编码和PixelCNN模型的心电信号降噪方法 | |
CN114169368A (zh) | 一种基于信号降噪自编码器sde的信号降噪方法 | |
CN116720059A (zh) | CEEMDAN-Grey-SSA弱信号提取方法 | |
CN111047537A (zh) | 一种图像去噪中恢复细节的系统 | |
CN115561817A (zh) | 一种基于多重注意力机制的沙漠地震去噪方法 | |
CN116184502A (zh) | 一种基于子空间投影注意力网络的井中das噪声抑制方法 | |
CN113205491A (zh) | 一种基于深度学习的织物起毛起球评级方法 | |
CN117743768B (zh) | 基于去噪生成对抗网络和扩散模型的信号去噪方法及系统 | |
CN115659150B (zh) | 一种信号处理方法、装置及设备 | |
Han et al. | Genetic algorithm assisted wavelet noise reduction scheme for chaotic signals | |
CN115017964B (zh) | 基于注意力机制稀疏表示的大地电磁信号去噪方法及系统 | |
CN111414808B (zh) | 基于平移不变分数阶小波稀疏表示的机械故障诊断方法 | |
CN114754857B (zh) | 一种两段式光纤传感水声信号补偿方法和装置 |
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 |