CN105748067B - 一种基于随机梯度自适应滤波的诱发脑电提取方法 - Google Patents

一种基于随机梯度自适应滤波的诱发脑电提取方法 Download PDF

Info

Publication number
CN105748067B
CN105748067B CN201610082340.2A CN201610082340A CN105748067B CN 105748067 B CN105748067 B CN 105748067B CN 201610082340 A CN201610082340 A CN 201610082340A CN 105748067 B CN105748067 B CN 105748067B
Authority
CN
China
Prior art keywords
signal
stimulation
time
moment
brain potential
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.)
Expired - Fee Related
Application number
CN201610082340.2A
Other languages
English (en)
Other versions
CN105748067A (zh
Inventor
李凌
李如峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Huinao Intelligent Technology Co ltd
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201610082340.2A priority Critical patent/CN105748067B/zh
Publication of CN105748067A publication Critical patent/CN105748067A/zh
Application granted granted Critical
Publication of CN105748067B publication Critical patent/CN105748067B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details 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)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Psychiatry (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Power Engineering (AREA)
  • Psychology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种基于随机梯度自适应滤波的诱发脑电提取方法,通过多道脑电测量系统采集多次重复刺激的脑电信号,其包含自发脑电、刺激施加后的自发脑电与诱发脑电的混合信号,并对各次混合信号叠加求平均,得到平均混合信号;以刺激施加时刻作为提取时间计数的起点,从平均混合信号中提取长度为N的一段信号作为第一主信号、对第一主信号延时一个时间间隔的信号作为第一参考信号,并输入LMS自适应滤波器,得到各通道的第一输出信号;再将第一输出信号作为第二参考信号、长度为N的一段自发脑电作为第二主信号,并输入LMS自适应滤波器,得到各通道的诱发脑电。本发明在没有足够先验知识的背景下,获得较高质量的诱发脑电,有效降低重复刺激次数。

Description

一种基于随机梯度自适应滤波的诱发脑电提取方法
技术领域
本发明属于神经信息技术领域,涉及一种诱发脑电提取方法,尤其涉及一种基于随机梯度(LMS,least mean square)自适应滤波的诱发脑电提取方法。
背景技术
自发脑电(自发脑电信号)是指利用精密电子设备在头皮未知无创伤地记录下来的大脑细胞群的自发性生物电活动;相对的,给予人们一定数量的重复刺激(视觉,听觉,体感等)则会产生诱发脑电(诱发脑电信号),诱发脑电被认为是研究人脑高级功能的新窗口。
然而,在实际应用中,由于诱发脑电总是淹没在较强的背景噪声(包括自发脑电、工频干扰、眼电、肌电、心电等)中,其幅值只有0.2~20uV,信噪比为0~10dB,而且脑电信号本身又具有随机性和非平稳性,诱发脑电和自发脑电在频谱上有相当大的重叠区,使得从观测的脑电信号中提取诱发脑电更加困难。传统的提取诱发脑电的方法是叠加平均法,该方法采用了如下假设:(1)诱发脑电和噪声为加法性的关系,且相互独立;(2)每次刺激后所得到的诱发脑电波形是一致的,即诱发脑电为确定性信号;(3)噪声与刺激无关,且是零均值的信号。在这样的假设前提下,通过对多次重复刺激测量结果求和平均,就可以得到减小噪声影响的诱发脑电信号。但这一方法要求同一采集通道重复刺激测量次数在100次以上,这需要耗费更多的时间来进行实验,而且并不是每次实验都会产生诱发脑电,因而实际处理结果要比理论值更小。另一方面,随着实验次数的增多,会使得受试者疲劳或者不适,影响实验结果。
为了解决传统提取方法中存在的问题,不同研究者提出了各种模型和方法,试图通过其他的方式来尽可能无损的提取诱发脑电信号。近年来广泛应用的有维纳滤波、小波、神经网络、时间序列分析法、独立分量分析法等。但由于不同方法适用的局限性,现在还没有一种公认的单次提取的方法得到广泛应用,同时各个类型的诱发脑电信号又具有不同的特征,在实际实验操作中,往往需要根据信号的特点来设计相应的使用方法,此时信号特点的获取就成为另一个问题。这也是这类方法发展所面临的困难之一。因此现在得到广泛认同的提取方法还是传统的叠加平均法。
发明内容
本发明的发明目的在于:为了克服现有技术中的缺陷,提高微弱信号提取的效率和精度,降低重复刺激的数量,提高应用度,提出了一种基于随机梯度(LMS,least meansquare)自适应滤波的诱发脑电提取方法。本发明能够在无法得到信号和噪声统计特性的先验知识的前提下,获得最佳的提取效果,能够在较少重复刺激次数的前提下,更准确地提取到诱发脑电信号。
本发明的一种基于随机梯度自适应滤波的诱发脑电提取方法,包括下列步骤:
脑电信号采集及预处理:
采用多通道脑电测量设备对K(K可以小于100)复刺激进行测量信号采集,所采集的测量信号包括刺激施加时刻、首次刺激施加时刻之前的脑电信号以及每次刺激施加时刻之后脑电信号,将首次刺激施加时刻之前的脑电信号定义为自发脑电信号,每次刺激施加时刻之后脑电信号定义为自发和诱发脑电的混合信号,其中每次刺激施加时刻之后的混合信号的采集长度相同;对每个通道的测量信号分别进行预处理:以每次刺激施加时刻作为起始点对齐K次混合信号,再进行叠加平均操作,得到平均混合信号,并记录平均混合信号的起始时刻;
分别获取每个通道的第一输出信号:
从平均混合信号中提取一段长度为N的信号作为第一主信号,其中第一主信号的起点对应平均混合信号的起始时刻;将第一主信号的起点延时一个时间间隔的时刻作为第一参考信号的起点,从平均混合信号中提取一段长度为N的信号作为第一参考信号;将第一主信号和第一参考信号输入LMS自适应滤波器进行滤波(自适应谱线增强处理),得到每个通道的第一输出信号;
分别获取每个通道的诱发脑电:
从自发脑电信号中提取一段长度为N的信号作为第二主信号,将第一输出信号作为第二参考信号;将第二主信号和第二参考信号输入随机梯度LMS自适应滤波器进行滤波(自适应噪声抵消处理),得到每个通道的诱发脑电。
综上所述,由于采用了上述技术方案,本发明的有益效果是:本发明可以在没有信号先验知识的前提下,有效地降低重复刺激次数,提取出高质量的诱发脑电。本发明通过对少数次重复刺激下的测量数据进行统计平均,然后利用LMS自适应滤波器对其依次进行自适应谱线增强、自适应噪声抵消处理,得到符合要求的诱发脑电,能有效的减少重复刺激次数,缩短时间,降低成本。同时,计算简单快捷,在样本较少的背景下具有更高的实用价值,同时自适应滤波方法的应用,进一步提高了诱发脑电信号提取的鲁棒性和稳定性。
附图说明
图1是实施例中,本发明对多通道脑电测量设备的其中一个通道的真实诱发脑电提取效果图;
图2是实施例中,本发明对64个通道的诱发脑电提取效果图;
图3是实施例中,基于AR参数模型的时间序列分析法对64个通路的诱发脑电提取效果图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合实施方式和附图,对本发明作进一步地详细描述。
实施例
本实施例中,以64通道的脑电测量设备为例对脑电信号进行采集,获取对应的测量信号。同时设置重复刺激次数为59次,采样率为1024Hz,时间长度为1.2秒,提取信号的长度N设置为200毫秒数据。通过执行下列步骤实现本发明的诱发脑电提取。
步骤A:脑电信号采集及预处理。
采用64通道脑电测量设备记录一个视觉重复刺激试验,得到64路测量信号。为了确保提取到长度为N的自发脑电,各路测量信号的数据记录初始时刻与首次刺激时刻的间隔应不低于200毫秒,本实施例中,记录首次刺激前的200毫秒数据,以及每次刺激施加时刻后的1000毫秒数据,因而获得一个三维信号:导联数(通道个数)*时间样本点*刺激重复次数。
以刺激施加时刻作为自发脑电和自发脑电+诱发脑电分割点,即将首次刺激施加时刻之前的脑电信号定义为自发脑电信号,每次刺激施加时刻之后脑电信号定义为自发和诱发脑电的混合信号。
对每个通道的测量信号分别进行预处理:以每次刺激施加时刻作为起始点对齐59次混合信号,再进行叠加平均操作,得到平均混合信号,并记录平均混合信号的起始时刻,如首次刺激施加时刻作为平均混合信号的起始时刻。
步骤B:分别获取各通道的第一输出信号。
以平均混合信号的起始时刻作为提取时间计数的零点,从平均混合信号提取一段长度为200毫秒的信号(200个数据点)作为LMS自适应滤波器的第一主信号。同时将第一主信号的起点延时一个时间间隔的时刻作为第一参考信号的提取时间计数的零点,从平均混合信号中提取一段长度为200毫秒的信号(200个数据点)作为LMS自适应滤波器的第一参考信号。
再将第一主信号、第一参考信号输入LMS自适应滤波器进行滤波,即进行自适应谱线增强处理,获得经过谱线增强处理后的200点长的平均样本,即各通道的第一输出信号。
步骤C:分别获取每个通道的诱发脑电。
将各通道的第一输出信号作为LMS自适应滤波器的第二参考信号,以各通道数据记录初始时刻开始到刺激施加时刻之间的数据点(200个)作为LMS自适应滤波器的第二主信号,再将第二主信号和第二参考信号输入LMS自适应滤波器进行滤波,即进行自适应噪声抵消处理,获得最后的输出信号,则得到每个通道要提取的诱发脑电。
对每个通道的测量信号分别执行步骤B、C后,得到每个通道的诱发脑电,则最后可获得一个诱发脑电矩阵数据(64*200)。
为了说明本发明诱发脑电提取的效果,将本发明提取的测量信号与传统的叠加平均的方法进行了比对,即59次重复刺激数据叠加后作为比对标准,用本发明的方法提取的诱发脑电与该比对标准求相关系数,获得全脑64个通道的提取效果。图1显示了本发明对其中一个通道的真实诱发脑电提取效果图,本发明与标准诱发脑电的相关系数达到了0.8099,而基于AR建模的时间序列分析法(利用AR模型构建白化滤波器对脑电信号中的自发脑电等非平稳信号进行白化滤除,然后进行逆变换得到诱发脑电信号。)与标准诱发脑电的相关系数仅为0.2582。说明了在没有足够先验知识的情况下,本发明能够有效的提取诱发脑电。图2和图3分别显示了本发明和基于AR参数模型的时间序列分析法对全脑64个通路59次重复刺激次数下的诱发脑电提取效果图。该图展现了在不同重复刺激次数的测量结果的输入下,不同方法与标准诱发脑电的相关系数结果的差异。相关系数的值越高,代表相关性越高,提取的效果越好。从图中可以清晰看出,本发明对所有通道都能取得很好的诱发脑电提取效果。在20次重复刺激次数下,本发明方法对全脑所有导联通道处理后的平均相关系数达到了0.8599,而在同样的重复次数刺激下,基于AR建模的时间序列分析法对全脑所有导联通道处理后的结果仅为0.5129。这进一步说明了本发明方法更具有普遍适用性,具有较好的实用价值。为了尽可能的用较少的重复次数获得更高质量的诱发脑电信号,通过记录5次、10次、15次下的重复刺激次数下的处理数据,利用本发明中所述方法进行处理,得到其平均相关系数分别为0.5893、0.7762、0.8255。随着测试次数的增加提取效果更加准确,20次的重复刺激次数的结果已经基本达到了应用要求。
以上所述,仅为本发明的具体实施方式,本说明书中所公开的任一特征,除非特别叙述,均可被其他等效或具有类似目的的替代特征加以替换;所公开的所有特征、或所有方法或过程中的步骤,除了互相排斥的特征和/或步骤以外,均可以任何方式组合。

Claims (4)

1.一种基于随机梯度自适应滤波的诱发脑电提取方法,其特征在于,包括下列步骤:
脑电信号采集及预处理:
采用多通道脑电测量设备对K次重复刺激进行测量信号采集,所采集的测量信号包括刺激施加时刻、首次刺激施加时刻之前的脑电信号以及每次刺激施加时刻之后脑电信号,将首次刺激施加时刻之前的脑电信号定义为自发脑电信号,每次刺激施加时刻之后脑电信号定义为自发和诱发脑电的混合信号,其中每次刺激施加时刻之后的混合信号的采集长度相同;
对每个通道的测量信号分别进行预处理:以每次刺激施加时刻作为起始点对齐K次混合信号,再进行叠加平均操作,得到平均混合信号,并记录平均混合信号的起始时刻;
分别获取每个通道的第一输出信号:
以平均混合信号的起始时刻作为提取时间计数的零点,从平均混合信号中提取一段长度为N的信号作为第一主信号,其中第一主信号的起点对应平均混合信号的起始时刻;
将第一主信号的起点延时一个时间间隔的时刻作为第一参考信号的提取时间计数的零点,从平均混合信号中提取一段长度为N的信号作为第一参考信号;
将第一主信号和第一参考信号输入随机梯度LMS自适应滤波器进行滤波,得到每个通道的第一输出信号;
分别获取每个通道的诱发脑电:
从自发脑电信号中提取一段长度为N的信号作为第二主信号,将第一输出信号作为第二参考信号;
将第二主信号和第二参考信号输入LMS自适应滤波器进行滤波,得到每个通道的诱发脑电。
2.如权利要求1所述的方法,其特征在于,K的取值小于100。
3.如权利要求2所述的方法,其特征在于,K的取值为20。
4.如权利要求1、2或3所述的方法,其特征在于,长度N的取值为200个数据点。
CN201610082340.2A 2016-02-05 2016-02-05 一种基于随机梯度自适应滤波的诱发脑电提取方法 Expired - Fee Related CN105748067B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610082340.2A CN105748067B (zh) 2016-02-05 2016-02-05 一种基于随机梯度自适应滤波的诱发脑电提取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610082340.2A CN105748067B (zh) 2016-02-05 2016-02-05 一种基于随机梯度自适应滤波的诱发脑电提取方法

Publications (2)

Publication Number Publication Date
CN105748067A CN105748067A (zh) 2016-07-13
CN105748067B true CN105748067B (zh) 2018-11-13

Family

ID=56330019

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610082340.2A Expired - Fee Related CN105748067B (zh) 2016-02-05 2016-02-05 一种基于随机梯度自适应滤波的诱发脑电提取方法

Country Status (1)

Country Link
CN (1) CN105748067B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106361328B (zh) * 2016-10-21 2019-07-19 电子科技大学 一种磁共振环境下的脑电信号提取方法
CN108836321A (zh) * 2018-05-03 2018-11-20 江苏师范大学 一种基于自适应噪声消除系统的脑电信号预处理方法
CN110327043B (zh) * 2019-08-05 2022-04-01 哈尔滨工业大学 一种基于稀疏建模的事件相关电位波形图谱求解方法
CN110833407B (zh) * 2019-09-17 2022-06-14 首都医科大学宣武医院 基于matlab的皮质-皮质间诱发电位数据处理方法
CN114041806A (zh) * 2021-11-09 2022-02-15 电子科技大学 一种自适应滤波的tms诱发电位提取装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5513649A (en) * 1994-03-22 1996-05-07 Sam Technology, Inc. Adaptive interference canceler for EEG movement and eye artifacts
CN1608581A (zh) * 2004-11-08 2005-04-27 中国科学院物理研究所 应用于心磁噪声抑制的自适应数字滤波方法
CN103405229A (zh) * 2013-07-12 2013-11-27 电子科技大学 一种基于靴带抽样的诱发脑电提取方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5513649A (en) * 1994-03-22 1996-05-07 Sam Technology, Inc. Adaptive interference canceler for EEG movement and eye artifacts
CN1608581A (zh) * 2004-11-08 2005-04-27 中国科学院物理研究所 应用于心磁噪声抑制的自适应数字滤波方法
CN103405229A (zh) * 2013-07-12 2013-11-27 电子科技大学 一种基于靴带抽样的诱发脑电提取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Application of adaptive noise cancellation with neural-network-based fuzzy inference system for visual evoked potentials estimation;Haie Yin et al;《Medical Engineering & Physics》;20040131;第26卷(第1期);第89-90页 *
基于LMS算法的自适应滤波器在水声信号处理中的应用;张艳喜 等;《科学技术与工程》;20070630;第7卷(第12期);第2830-2831页 *

Also Published As

Publication number Publication date
CN105748067A (zh) 2016-07-13

Similar Documents

Publication Publication Date Title
CN105748067B (zh) 一种基于随机梯度自适应滤波的诱发脑电提取方法
CN113288181B (zh) 基于稳态视觉诱发电位脑电信号识别的个体模板重构方法
CN110353702A (zh) 一种基于浅层卷积神经网络的情感识别方法及系统
CN101515200B (zh) 基于瞬态视觉诱发脑电的目标选择方法
Mc Laughlin et al. Cochlear implant artifact attenuation in late auditory evoked potentials: A single channel approach
US8484270B2 (en) System and method for adaptive stimulus-response signal filtering
CN106236080B (zh) 基于多通道的脑电信号中肌电噪声的消除方法
Yang et al. Noise characterization, modeling, and reduction for in vivo neural recording
CN104688220A (zh) 一种去除脑电信号中眼电伪迹的方法
CN109965869A (zh) 基于脑源域空间的mi-eeg识别方法
Luke et al. Kalman filter based estimation of auditory steady state response parameters
CN108403108A (zh) 基于波形优化的阵列式表面肌电信号分解方法
Hansen et al. Unmixing oscillatory brain activity by EEG source localization and empirical mode decomposition
Stevenson et al. A nonlinear model of newborn EEG with nonstationary inputs
CN113576498B (zh) 基于脑电信号的视听觉美学评价方法及系统
CN103989482A (zh) 用于幻听检测的声音刺激器及用于幻听检测的装置
David et al. Decoupling action potential bias from cortical local field potentials
CN106344011B (zh) 一种基于因子分析的诱发脑电信号提取方法
CN116269445A (zh) 一种用于ssvep短时窗信号的精准目标识别方法
CN109002798A (zh) 一种基于卷积神经网络的单导视觉诱发电位提取方法
CN103405229B (zh) 一种基于靴带抽样的诱发脑电提取方法
Nawrocka et al. 780. Biomedical signal identification and analysis
El Khadiri et al. A Comparison of the Denoising Performance Using Capon Time-Frequency and Empirical Wavelet Transform Applied on Biomedical Signal.
CN113545792A (zh) 一种基于tms-eeg的神经生物标记物的检测方法
CN107657278A (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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200807

Address after: 710000 25 / F, block D, Tsinghua Science Park, Keji 2nd Road, Zhangba Street office, hi tech Zone, Xi'an City, Shaanxi Province

Patentee after: Xi'an Huinao Intelligent Technology Co.,Ltd.

Address before: 611731 Chengdu province high tech Zone (West) West source Avenue, No. 2006

Patentee before: University of Electronic Science and Technology of China

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20181113