CN111789592B - 一种基于拓扑特征融合的脑电识别方法 - Google Patents
一种基于拓扑特征融合的脑电识别方法 Download PDFInfo
- Publication number
- CN111789592B CN111789592B CN202010620170.5A CN202010620170A CN111789592B CN 111789592 B CN111789592 B CN 111789592B CN 202010620170 A CN202010620170 A CN 202010620170A CN 111789592 B CN111789592 B CN 111789592B
- Authority
- CN
- China
- Prior art keywords
- fusion
- hidden layer
- electroencephalogram
- local
- topological
- 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
- 238000000034 method Methods 0.000 title claims abstract description 52
- 230000004927 fusion Effects 0.000 title claims abstract description 41
- 238000005457 optimization Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims abstract description 5
- 230000003925 brain function Effects 0.000 claims abstract 6
- 230000009471 action Effects 0.000 claims description 10
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000011160 research Methods 0.000 claims description 8
- 238000010606 normalization Methods 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 5
- 238000005452 bending Methods 0.000 claims description 4
- 210000000707 wrist Anatomy 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 3
- 238000004891 communication Methods 0.000 claims description 2
- 230000000295 complement effect Effects 0.000 claims description 2
- 230000007613 environmental effect Effects 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000004422 calculation algorithm Methods 0.000 abstract description 5
- 239000000284 extract Substances 0.000 abstract 1
- 210000004556 brain Anatomy 0.000 description 13
- 230000008569 process Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 230000001537 neural effect Effects 0.000 description 4
- 241000282414 Homo sapiens Species 0.000 description 3
- 230000004913 activation Effects 0.000 description 3
- 210000004247 hand Anatomy 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000012549 training Methods 0.000 description 3
- 230000007177 brain activity Effects 0.000 description 2
- 238000007635 classification algorithm Methods 0.000 description 2
- 238000002790 cross-validation Methods 0.000 description 2
- 238000000354 decomposition reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 210000003205 muscle Anatomy 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 206010033799 Paralysis Diseases 0.000 description 1
- 229910021607 Silver chloride Inorganic materials 0.000 description 1
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 210000005252 bulbus oculi Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000005465 channeling Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000013145 classification model Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 210000003414 extremity Anatomy 0.000 description 1
- 238000007499 fusion processing Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- HKZLPVFGJNLROG-UHFFFAOYSA-M silver monochloride Chemical compound [Cl-].[Ag+] HKZLPVFGJNLROG-UHFFFAOYSA-M 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
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/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/7253—Details of waveform analysis characterised by using transforms
- A61B5/726—Details of waveform analysis characterised by using transforms using Wavelet transforms
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Physiology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Psychiatry (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Evolutionary Computation (AREA)
- Fuzzy Systems (AREA)
- Mathematical Physics (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明公开了一种基于拓扑特征融合的脑电识别方法,本发明为了提高手部动作分类的性能,采用了一种LBELM和基于似然同步分析方法量化脑功能网络的新方法。与传统思想不同,本发明对二值FBN提取了两种最优的拓扑特征,并通过LBELM对两种最优的拓扑特征进行LBELM特征空间层的融合;同时本发明对LBELM做出了进一步的改进,增加了隐藏层参数优化以获得更高、更稳定的识别效果;使用了基于MSEPRESS的留一法优化算法,以获得最优的正则化系数和融合比例。
Description
技术领域
本发明属于脑电信号特征分析领域,公开了一种基于局部平衡极限学习机(locals-balanced Extreme Learning Machine,LBELM)进行拓扑特征融合的脑电识别方法。
背景技术
人类在运动想象(Motor Imagery,MI)任务中神经活动的表征和描述仍然是一个重要的挑战。目前,基于MI的神经康复技术和基于脑活动解释的脑-机接口(Brain-computer Interface,BCI)引起了人们的极大兴趣,试图为(部分或完全) 瘫痪患者提供一种潜在的新的脑肌沟通工具。最近的一些综述概述了BCI***的最新进展,并讨论了表征所记录的神经活动的常用方法。基于脑功能网络分析(Functional Brain Network,FBN)的功能连接性研究成为了研究热点,这些图论方法已被广泛应用并成功地描述人脑动力学。
为了利用脑电构造FBN,采用连通性测度来估计节点(电极)之间传递的信息值。文献中有许多线性和非线性连接性度量被用来构造FBN,如互信息、熵、相关性和Granger 因果关系。线性连接性测量通常无法识别大脑的非线性行为。因此,为了分析一个高度非线性的脑电信号,研究人员采用非线性方法来构建FBN。
目前,在如何最好地检测神经生理数据中的非线性相互依赖性方面还没有达成共识。事实上,不同的算法已经被证明可以检测大脑区域之间的非线性相互作用。如前所述,同步度量旨在量化多个耦合动力***之间的协调关系。文献中定义了不同类型的同步。三个主要类别包括完全同步、广义同步和相位同步。完全同步是当相同的***彼此有强耦合时实现的,而一般同步被定义为两个***的状态变量之间的函数关系,其中一个***作为驱动信号,另一个***充当响应信号。因此,一般同步放宽了完全同步所需的等状态准则,因为大多数实际***不满足强耦合这种假设。最后,相位同步被定义为当振幅保持混沌且弱相关时样本输出相位的全局夹带。
两个动力***相互作用的最普遍形式是广义同步。非线性非平稳***间相互关系的一种最有前途的量化方法是同步似然(Synchronization Likelihood,SL),它是基于嵌入向量(相似模式)在多个动态***中自回归的似然性。以往研究中还使用了其他脑连接性测量方法,如部分定向相关性和相位锁值来表征大脑活动。
虽然FBN是一种重要的工具,可以用来开发更有效的BCI,但它可能不足以将人类大脑动力学完全描述为一个复杂的解剖和功能网络中的信息流。
发明内容
本发明针对现有技术的不足,提出了一种基于拓扑特征融合的脑电识别方法。
本发明提出一种基于拓扑特征融合的脑电识别方法。与传统方法不同的是,本发明不是采用单一的拓扑特征,而是使用多种拓扑特征在LBELM的特征空间进行融合。同时在原LBELM的基础上,辅以优质隐藏层单元参数,以获得最佳的MI脑电识别率。
为实现以上目的,本发明方法主要包括以下步骤:
步骤(1)采集MI脑电。实验设备的采样频率设置为5000Hz。根据实验要求,受试者按照随机动作顺序进行左右手的曲腕、曲肘和握拳动作。
步骤(2)信号预处理。信号采集过程中,采集***的电极帽基于高质量Ag/AgCl 传感器的主动电极,结合新型集成噪声衰减电路,可以达到较低的噪声水平;同时做足完备的实验范式以防止声源、光源等外界干扰。信号采集结束后,由于眼电、心电和工频干扰,有必要对采集信号在进行滤波处理。
步骤(3)构建基于SL的FBN。脑电信号包含了数十亿个相互连接的神经元组成的复杂而密集的网络的信息。近年来,研究者们应用基于图论的方法成功地研究了这些复杂的脑网络。
步骤(4)构建二值FBN。使用度量量化大脑节点之间信息,量化出的连接往往会存在弱连接和伪连接的情况。为了防止这些连接对实验结果造成干扰,我需要选择合适的阈值来消除这些弱连接。
步骤(5)提取拓扑特征。实际网络都兼有确定和随机两大特征,确定新的法则或特征通常隐藏在统计性质中,因此,对复杂网络各种性质的统计描述十分重要。虽然目前提出的大多数网络统计特征仅描述网络拓扑特征,即只考虑节点间的位置关系而不考虑它们之间的形状的大小,但是由于网络节点的连边表示它们形形色色的相互作用,所以这种统计描述也包含动力学成分,具有非常重要的意义。
步骤(6)使用LBELM分类。ELM发展自单隐层神经网络(single hidden-layerfeedforward neural networks,SLFNs),它被认为是一个“广义”的SLFNs。对于一个 n维的,拥有N个样本,m个类别的训练样本L个隐藏单元,激活函数为gi(sj) 的SLFNs输出为
其中sj=[sj1,sj2,…,sjn]T为第j个输入,tj=[tj1,tj2,…,tjm]T为对应于第j个输入的期望输出;bi为第i个隐藏单元的激活阈值;ci=[ci1,ci2,…,cin]连接了输入层和第i个隐藏单元;βi=[βi1,βi2,…,βim]T连接了第i个隐藏单元和输出层;oj为SLFNs对应于第j个输入的实际输出。
SLFNs的所有参数需要使用梯度下降方法进行收敛,所以SLFNs的学***方误差意义下的逼近误差:
其中,
为了提高ELM的稳定性和泛化性能,本发明应用了正则化以防止过拟合,
这里的λ是正则化系数,正则化之后ELM的隐藏层输出矩阵为
其中I为单位矩阵。
特征在输入ELM之前需要归一化。由于篇幅的限制,本文只讨论了模值归一化和高斯分布生成ELM随机参数的情况。
假设样本s由多个互补特征构成,即s=[x1;x2;…;xq],其中k=1,2,…,q。他们的维数为ak、模值为||xk||2,k=1,2,…,q。同时对应的假设某个隐藏层单元的输入权重矩阵由多个局部隐藏层单元构成,即 c=[w1|w2|…|wq],其中k=1,2,…,q。wk的维数一一对应了xk的维数。ELM的隐藏层参数由高斯分布N(0,1)生成,所以 (事实上,模长会随机分布在上述值的附近)。
在单个隐藏层单元中,观察公式6中各个局部特征在ELM特征空间中的融合过程,鉴于的范围是[-1,1],可以发现局部特征的融合比例为这表明局部特征的融合比例只与输入局部特征的模和维数有关,各个局部特征在ELM特征空间中是无法自由融合的。
鉴于局部特征的融合在激活函数进行非线性化之前,所以不需要对ELM网络结构进行改动,只需要对输入数据进行局部的拆分和融合。这样就避免了对所有隐藏层单元“理解”进行局部的拆分和融合。
为了实现输入空间中局部的拆分与融合,相对于极限学***衡极限学习机LBELM在归一化时增加了一个变量B=[b1,b2,…,bq],样本归一化后变为其中归一化后s′仍然满足模值为一,与此同时可以发现矩阵H中的输入特征相关项
公式(7)中使用到了约等号,是因为局部隐藏单元wk是由N(0,1)随机产生,所以从单个隐藏层单元来分析,这个波动会使得该单元的融合比例在附近浮动,这并不会对该单元训练和测试过程中的融合产生影响;但是从多个隐藏层单元来分析,所有单元的融合比例都在附近浮动,这表示所有单元的融合比例都略微不一样。所以隐藏层的随机产生所造成的波动会对融合产生干扰。隐藏层局部输入权重间的模值比||w1||:||w2||:…:||wq||越接近则分类效果越好,即点E=[||w1||,||w2||,…,||wq||]和点的距离越近越好。。
为了克服随机参数造成的干扰,本文在生成隐藏层单元时设置阈值σ,只有当norm(E-F)≤σ时,将E对应的隐藏层单元加入到LBELM网络中。
步骤(7)使用留一法(leave-one-out cross validation,LOO)进行参数寻优。LOO方法是机器学***方和(predictedresidual sum of squares,PRESS)统计量,该统计量可有效计算线性模型(而不是多次再训练)的LOO方法的均方误差(mean square error,MSE)。
本发明与已有的诸多脑电信号的方法相比,具有如下特点:
对二值FBN提取了两种最优的拓扑特征,并通过LBELM进行特征空间层的融合;对LBELM做出了进一步的改进,增加了隐藏层参数优化以获得更高、更稳定的识别效果;使用了基于MSEPRESS的LOO优化算法,以获得最优的正则化系数和融合比例。
附图说明
图1为本发明的实施流程图;
图2为本发明所进行实验中使用的电极分布图;
图3为正则化ELM分类算法的伪代码;
图4为LBELM分类算法的伪代码;
图5为8位志愿者的动作识别率;
具体实施方式
下面结合附图对本发明的实施例作详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
本发明的实施例采集的受试者有8名志愿者组成,其中包括4名男性、4名女性,他们分别用左右手执行腕曲,肘曲和握拳共计6种动作。发明在数据预处理阶段使用了小波去噪,特征提取阶段采用基于Manhattan距离的SL构造加权FBN,二进制FBN由阈值法构建。模式识别阶段使用LBELM融合二值FBN的两种拓扑特征并进行分类任务。为了获得最佳参数,本发明还采用了留一法(leave-one-out cross validation,LOO) 来优化本发明中无法靠经验取值的参数。
下文将介绍参与者、实验任务、数据集的信息与这项工作所涉及的理论,介绍方法的配置,并说明所提出方案的性能。最后与已有的相关研究进行比较,提出预期的工作。
如图1所示,本实施例包括如下步骤:
步骤(1)采集脑电数据。实验室的环境设置必须能够有效控制住声音、光线以及温度等这些外界干扰。具体要求有:1)实验室的所有室外光源必须通过窗帘进行遮蔽,打开室内的稳定光源,以排除光线扰动;2)在数据采集过程中,利用空调将温度控制在25±2℃以内,以防止温度扰动;3)实验室采取隔音措施且采集过程中,采集人员不能随意走动和交流,以防止声音的干扰。实验采集使用的是actiCAP 128导的电极帽,其电极布局如图2 ( a )所示。实施列采集在保证目标研究区域覆盖的情况下,选择了其中如图2 ( b ) 所示的40个电极。实验前,受试者已被告知完整的实验过程和相关细节,以便熟悉整个实验过程。要求受试者避免身体和头部摇晃,以及所有超出实验计划的肢体运动,以获得最小程度的干扰。数据采集前,在受试者的电极帽内注入导电胶以降低阻抗。调整电极的接触条件,直到阻抗降到阈值以下,以确保采集信号的质量。
实验采样频率设置为5000Hz。根据实验要求,志愿者用左手或右手分别进行腕屈、肘屈和握拳的dMI动作。具体步骤如下:
1.采集开始前,受试者保持静止状态15秒。在此过程中,受试者必须坐在扶手椅上,放松手部肌肉,冷静下来,尽量减少干扰。
2.受试者根据要求同时想象和模仿相应的动作。动作执行时间要求为作5秒,做完一次动作休息15秒,然后再执行下一个动作。
3.一组动作包含上述6个动作的10次。将这60个动作的执行顺序随机打乱,以防止意识的惯性思维。
步骤(2)信号预处理。保证脑电信号质量是进行准确有效研究的基础。在脑电信号采集过程中,不可避免地会受到各种信号和噪声的干扰。由于脑电信号非常微弱,很容易受到多种干扰源的干扰,包括眼球运动和眨眼伪影、心电图和工频干扰。因此,为了保证基于EEG信号的研究成果的可靠性,有必要对信号做相应的去噪处理,以尽可能减少噪声对研究的干扰。
本研究选取小波阈值法进行去噪。具体步骤如下:
1.确定小波基和分解层,对提取的脑电信号进行分解。选择sym6作为小波基,将分解层数设为6,然后利用MATLAB中的wavedec函数实现第一步。
2.获取去噪的全局阈值。首先使用wnoisest函数来估计细节系数的标准差,其水平为1。在上述阶段之后,使用带有惩罚项的函数wbmpen来捕获全局阈值。
3.对处理后的信号进行去噪和重构。为了操作简单,本文选择了硬阈值法。利用保持逼近系数的函数wdencmp在MATLAB中实现去噪和重构。
步骤(3)通过SL构建加权FBN。C.J.Stam提出的SL分析方法是多元数据集广义同步的一种无偏测度。SL可以说是最流行的神经生理数据中估计GS的指标。该指数与广义互信息的概念密切相关,它依赖于同时发生的模式的检测,这些模式可以是复杂的,并且在信号之间有很大的不同。选取脑电采集的电极点为网络的节点,使用SL量化两两节点之间的SL值。
SL分析的参数设置为:嵌入维数d=10、嵌入滞后τ=10、泰勒窗[t1,t2]= [100,250]、参考概率Pref=0.05、运算速度speed=1500。
步骤(4)构建二值FBN。本发明采用阈值法作为权重FBN网络的二值化方法。选取合适的阈值Th,当节点间的SL值大于等于Th,则认为该节点间的网络连接处于联通状态,置为1;当节点间的SL值小于Th时,认为该节点间的网络连接为伪连接,置为0。
阈值法采用的阈值Th=0.17。
步骤(5)提取拓扑特征。本发明选取了度、聚类系数、介数中心度、核数作为候选拓扑特征。经过分析比较后采用了最富有信息的两个拓扑特征—度和核数。表1分别展示了阈值法二值FBN和通过构造MST法二值FBN四种候选拓扑特征的分类性能(统计量采用了ELM分类的MSEPRESS,数值越小表明该特征在分类中越优异);
表1
步骤(6)使用LBELM分类。LBELM对标准正则化ELM(Regularized ELM,RELM) 的局部特征融合缺陷进行了改进,使得LBELM拥有了按可控比例在ELM特征空间融合局部特征的能力。相对标准正则化ELM(伪代码如图3所示),LBELM(伪代码如图4所示) 除了需要提供正则化系数λ之外,还需要提供三个额外参数:
正则化系数λ。该参数可以控制隐藏层输出矩阵的稀疏程度,防止分类模型过拟合;参数λ设置为LOO算法输出的最优正则化系数;
输入特征的结构参数K=[a1,a2,…,aq]。该参数表示了输入特征的局部结构,从而把输入特征拆分成局部特征。对于本发明在该示例中的应用为两个40维的拓扑特征,所以该实例的参数K=[40,40];
局部特征融合比例B。该参数控制了局部特征之间的融合比例。参数B设置为LOO算法输出的最优融合比例;
隐藏层单元筛选阈值σ。该参数控制了生成的局部隐藏层单元间的模值比例。该实例的参数σ=0.01。
步骤(7)使用LOO进行参数寻优。基于MSEPRESS的参数寻优方法将会遍历所有的参数对,选择分类的MSEPRESS最小的参数对作为最优参数组合。融合比例B的候选区间为10[-1:0.1:1]/1;正则化系数候选区间为e[-4:0.1:6]。
该实例的动作分类正确率如图5所示,使用RELM进行融合无法确保融合的效果,而使用LBLEM就可以按照最优的融合比例进行融合,实现了识别率的较大提升。
Claims (6)
1.一种基于拓扑特征融合的脑电识别方法,其特征在于,包括如下步骤:
步骤(1)采集脑电数据;
步骤(2)信号预处理;
步骤(3)通过似然同步性构建加权脑功能网络;
步骤(4)构建二值脑功能网络;
步骤(5)提取拓扑特征;
提取步骤(4)中二值脑功能网络的两种拓扑特征;
步骤(6)使用局部平衡极限学习机实现特征融合并分类;
通过局部平衡极限学习机实现步骤(5)中两种拓扑特征在极限学习机特征空间的最优比例融合,以实现分类效果的提升,并通过最优的正则化系数获得最佳的分类泛化性能;
中局部“理解”之间的比值变为即通过改变融合比值B就能实现融合比例的调节;是隐藏层对第k个局部输入特征的“理解”;样本s由多个互补特征构成,即s=[x1;x2;…;xq],其中维数为ak、模值为||xk||2,其中表示第k个局部特征中第ak个元素;同时对应的假设某个隐藏层单元的输入权重矩阵由多个局部隐藏层单元构成,即c=[w1|w2|…|wq],其中wk的维数一一对应了xk的维数;ELM的隐藏层参数由高斯分布N(0,1)生成,所以
因为隐藏层的随机产生所造成的波动会对融合产生干扰,所以更加期待隐藏层局部输入权重间的模值比||w1||:||w2||:…:||wq||接近于即点E=[||w1||,||w2||,…,||wq||]和点的距离近;
为了克服随机参数造成的干扰,在生成隐藏层单元时设置阈值σ,只有当norm(E,F)≤σ时,将E对应的隐藏层单元加入到局部平衡极限学习机网络中;
步骤(7)使用留一法进行参数寻优;
设置融合比值B候选区间为10[-1:0.1:1];正则化系数λ候选区间为e[-4:0.1:6];通过留一法得到分类效果最佳的融合比值B和正则化系数λ。
2.根据权利要求1所述的一种基于拓扑特征融合的脑电识别方法,其特征在于:所述的采集脑电数据;具体为:
进行实验室的环境布置,必须能够有效控制住声音、光线以及温度干扰;在数据采集之前,将电极帽安装在受试者头上并完成电极阻抗测试,以确保收集信号的质量;实验样本频率设置为5000Hz,受试者按照实验范式的要求完成数据采集。
3.根据权利要求1所述的一种基于拓扑特征融合的脑电识别方法,其特征在于:采集脑电数据时,受试者按照随机动作顺序进行左右手的曲腕、曲肘和握拳动作。
4.根据权利要求1所述的一种基于拓扑特征融合的脑电识别方法,其特征在于:信号预处理;具体为:
对记录的信号进行小波阈值法去噪,以尽可能减少噪声对研究的干扰。
5.根据权利要求1所述的一种基于拓扑特征融合的脑电识别方法,其特征在于:通过似然同步性构建加权脑功能网络;具体为:
选取脑电采集的电极点为网络的节点,使用SL量化两两节点之间的SL值。
6.根据权利要求1所述的一种基于拓扑特征融合的脑电识别方法,其特征在于:构建二值脑功能网络;具体为:
选取阈值Th,当节点间的SL值大于等于Th,则认为该节点间的网络连接处于联通状态,置为1;当节点间的SL值小于Th时,则认为该节点间的网络连接为伪连接,置为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010620170.5A CN111789592B (zh) | 2020-06-30 | 2020-06-30 | 一种基于拓扑特征融合的脑电识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010620170.5A CN111789592B (zh) | 2020-06-30 | 2020-06-30 | 一种基于拓扑特征融合的脑电识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111789592A CN111789592A (zh) | 2020-10-20 |
CN111789592B true CN111789592B (zh) | 2022-04-05 |
Family
ID=72809803
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010620170.5A Active CN111789592B (zh) | 2020-06-30 | 2020-06-30 | 一种基于拓扑特征融合的脑电识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111789592B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113398422B (zh) * | 2021-07-19 | 2022-08-19 | 燕山大学 | 基于运动想象-脑机接口和虚拟现实的康复训练***及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109711278A (zh) * | 2018-12-07 | 2019-05-03 | 浙江大学 | 一种低复杂度的脑电波信号压缩和分类方法 |
CN109893126A (zh) * | 2019-03-21 | 2019-06-18 | 杭州电子科技大学 | 基于脑功能网络特征的癫痫发作预测方法 |
CN109948640A (zh) * | 2018-12-26 | 2019-06-28 | 杭州电子科技大学 | 基于双参数核优化型极限学习机的肌电信号分类方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9974463B2 (en) * | 2009-08-21 | 2018-05-22 | Beth Israel Deaconess Medical Center, Inc. | Hand-held device for electrical impedance myography |
-
2020
- 2020-06-30 CN CN202010620170.5A patent/CN111789592B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109711278A (zh) * | 2018-12-07 | 2019-05-03 | 浙江大学 | 一种低复杂度的脑电波信号压缩和分类方法 |
CN109948640A (zh) * | 2018-12-26 | 2019-06-28 | 杭州电子科技大学 | 基于双参数核优化型极限学习机的肌电信号分类方法 |
CN109893126A (zh) * | 2019-03-21 | 2019-06-18 | 杭州电子科技大学 | 基于脑功能网络特征的癫痫发作预测方法 |
Non-Patent Citations (4)
Title |
---|
An Improved Recognition Approach for Noisy Multispectral Palmprint by Robust L2 Sparse Representation with a Tensor-Based Extreme Learning Machine;Dongxu Cheng等;《sensors》;20190102;第1-25页 * |
Image Classification Using Low-Rank Regularized Extreme Learning Machine;Qin Li等;《IEEE Access》;20190104;第877-883页 * |
Using the Partial Directed Coherence to Assess Functional Connectivity in Electroencephalography Data for Brain–Computer Interfaces;Jorge Antonio Gaxiola-Tirado等;《 IEEE Transactions on Cognitive and Developmental Systems》;20171123;第776-783页 * |
基于极限学习机的大坝多测点模型分析;邵楠等;《城市勘测》;20200531(第2期);第169-172页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111789592A (zh) | 2020-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Seal et al. | DeprNet: A deep convolution neural network framework for detecting depression using EEG | |
US6546378B1 (en) | Signal interpretation engine | |
CN111209885A (zh) | 一种手势信息处理方法、装置、电子设备及存储介质 | |
CN113693613B (zh) | 脑电信号分类方法、装置、计算机设备及存储介质 | |
CN114533086B (zh) | 一种基于空域特征时频变换的运动想象脑电解码方法 | |
CN109497996B (zh) | 一种微状态eeg时域特征的复杂网络构建及分析方法 | |
CN110969108A (zh) | 一种基于自主运动想象脑电的肢体动作识别方法 | |
CN111657941B (zh) | 基于肌肉核心激活区域的电极校正及肌电模式识别方法 | |
CN110522412B (zh) | 基于多尺度脑功能网络分类脑电信号的方法 | |
Kumar et al. | OPTICAL+: a frequency-based deep learning scheme for recognizing brain wave signals | |
CN112861604A (zh) | 一种与用户无关的肌电动作识别与控制方法 | |
CN110013248A (zh) | 脑电张量模式识别技术与脑机交互康复*** | |
CN110333783A (zh) | 一种用于鲁棒肌电控制的无关手势处理方法及*** | |
CN112401905B (zh) | 一种基于源定位和脑网络的自然动作脑电识别方法 | |
Jinliang et al. | EEG emotion recognition based on granger causality and capsnet neural network | |
Liang et al. | Identity recognition using biological electroencephalogram sensors | |
CN113378650B (zh) | 一种基于脑电源成像和正则化共空间模式的情绪识别方法 | |
CN111789592B (zh) | 一种基于拓扑特征融合的脑电识别方法 | |
CN114027786B (zh) | 基于自监督式记忆网络的睡眠呼吸障碍检测方法及*** | |
Li et al. | Granger causal inference based on dual laplacian distribution and its application to MI-BCI classification | |
CN111931578B (zh) | 一种基于最小生成树和区域双层网络的脑电识别方法 | |
CN110738093B (zh) | 基于改进小世界回声状态网络肌电的分类方法 | |
CN113128384A (zh) | 一种基于深度学习的脑卒中康复***脑机接口软件关键技术方法 | |
Ren et al. | Extracting and supplementing method for EEG signal in manufacturing workshop based on deep learning of time–frequency correlation | |
Reddy et al. | Optimization Driven Spike Deep Belief Neural Network classifier: a deep-learning based Multichannel Spike Sorting Neural Signal Processor (NSP) module for high-channel-count Brain Machine Interfaces (BMIs) |
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 |