CN111563581A - 一种基于小波相干的脑肌功能网络构建方法 - Google Patents
一种基于小波相干的脑肌功能网络构建方法 Download PDFInfo
- Publication number
- CN111563581A CN111563581A CN202010460960.1A CN202010460960A CN111563581A CN 111563581 A CN111563581 A CN 111563581A CN 202010460960 A CN202010460960 A CN 202010460960A CN 111563581 A CN111563581 A CN 111563581A
- Authority
- CN
- China
- Prior art keywords
- wavelet
- signals
- network
- muscle function
- brain
- 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
-
- 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/06—Physical realisation, i.e. hardware implementation of neural networks, neurons or parts of neurons
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/369—Electroencephalography [EEG]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/389—Electromyography [EMG]
-
- 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
- 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)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Psychiatry (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Evolutionary Computation (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Computational Linguistics (AREA)
- Neurology (AREA)
- Fuzzy Systems (AREA)
- Psychology (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及一种基于小波相干的脑肌功能网络构建方法。在人体运动过程中,大脑皮层通过脊髓和周围神经控制肌肉组织的运动,使得肢体可以完成一定的运动功能,而肢体的运动又可以通过传入神经反馈给大脑皮层。脑肌电功能耦合可以从生理上反映皮层与肌肉之间的关系。为了更好的探索人体行为感知,应用脑电与肌电信号研究运动过程中不同通道信号之间的脑肌功能耦合强度,并以此作为复杂网络加权边的值,构建脑肌功能网络。将两两信号之间的小波相干值作为脑肌功能耦合强度,同时为了增大不同动作网络模型的差异,应用阈值法和固定加权边法简化网络模型。实验结果证明,构建的脑肌功能网络具有小世界特性,应用网络特征,可对人体行为动作进行识别。
Description
技术领域
本发明属于信号处理领域,涉及一种基于小波相干的脑肌功能网络的构建方法。
背景技术
脑电信号(EEG)与表面肌电信号(sEMG)分别是大脑活动与肌肉活动产生的一种非线性、非平稳信号。EEG广泛应用于神经系统疾病、脑机接口、运动想象等领域。而sEMG广泛应用于康复训练、医学传感器、机械控制等领域。
在人类自主运动过程中,除了大脑不同区域间的协同作用外,运动皮层、大脑体感皮层和运动肌肉组织的神经细胞之间也会自动发生同步。大脑皮层、运动神经和肌肉组织之间的相互作用构成了脑肌功能耦合(CMC)。近年来,越来越多的研究者对脑电、肌电的耦合分析进行了研究。研究脑肌电耦合的方法主要有相干性分析、相位同步指标、广义同步指标、格兰杰因果关系以及信息理论等。在这些方法中,相干分析是应用最广泛的。
复杂网络理论在脑电信号的研究中得到了广泛的应用,尤其是脑功能网络的研究已经有二十多年的历史。网络的拓扑性质指的是在复杂网络中不体现具体的节点位置与连接情况就能呈现的特征。近些年来,伴随着脑影像功能网络的研究深入,许多测度指标被设定来衡量网络整体或者局部特征。这些测度中最主要的有:度、特征路径长度、聚类系数等。而常见的网络模型有随机网络、规则网络、小世界网络和无标度网络。在大部分真实网络具有小世界特性、大型网络具有无标度特性这两个发现得到验证后,脑功能网络进入了快速发展时期。脑功能网络被广泛应用于不同脑功能或功能障碍的脑科学研究。
发明内容
本发明涉及一种基于小波相干的脑肌功能网络的构建方法。首先,选取合适的脑电信号与肌电信号,然后计算脑电信号与肌电信号两两通道之间的小波相干值,以此作为复杂网络的加权边值,得到相应的邻接矩阵。接着将得到的邻接矩阵分为三部分,EEG-EEG、EMG-EMG、EEG-EMG,将三个部分的小波相干值分别归一化,得到归一化邻接矩阵。最后运用阈值法去除归一化邻接矩阵中较弱的加权边,即将邻接矩阵中较弱的边置为零。由此得到的邻接矩阵即为简化后的脑肌功能网络。实验过程中,我们将两两信号之间的小波相干值作为脑肌功能耦合强度,同时为了增大不同动作网络模型的差异,我们分别应用阈值法简化网络模型。最终的实验结果证明,我们构建的脑肌功能网络具有小世界特性,应用网络特征,可对人体行为动作进行识别。
为了实现以上目的,本发明方法主要包括以下步骤:
步骤(1).首先选取与屈腕、伸腕、握拳三个行为动作紧密相关的8通道脑电信号与4通道肌电信号并进行预处理,8通道脑电信号为C1、C2、C3、C4、C5、FC1、FC3、CP1,4通道肌电信号采集自尺侧伸腕肌、桡侧屈腕肌、屈指浅肌、伸指肌四块肌肉;
步骤(2).计算步骤(1)中脑电信号与肌电信号两两通道之间的小波相干值,以此作为复杂网络的加权边值,得到相应的邻接矩阵;
步骤(3).将步骤(2)得到的邻接矩阵分为三部分,EEG-EEG、EMG-EMG、EEG-EMG,将三个部分的小波相干值分别归一化,得到归一化邻接矩阵;
步骤(4).运用阈值法去除归一化邻接矩阵中较弱的加权边,即将邻接矩阵中较弱的边置为零。由此得到的邻接矩阵即为简化后的脑肌功能网络;
所述的小波相干值计算方法具体如下:
1)定义信号x的CWT的范数平方为小波能量,即
WPx(t,f)=||CWTx(t,f)||2
WP是时间t与小波中心频率f的函数,CWT为连续小波变换;
2)信号x与y的交叉小波变换定义如下:
3)与传统的相干性估计类似,我们将小波相干值进行平滑。基于时间轴的平滑函数定义如下:
λ=t/a,c1为归一化常数,∧为卷积运算。基于尺度轴的平滑函数定义如下:
Sscale(CWTx(t,f))=CWTx(t,f)∧c2Π(0.6a)
c2为归一化常数,Π为矩形函数。矩形函数的长度是由尺度去相关决定的。
实验时我们设定为0.6。
4)小波相干的定义如下:
a-1为尺度的逆,用于归一化XWT。Schwartz不等式确保WCxy∈[0,1]。
平滑过程S定义如下:
S(w)=Sscale[Stime(W)]
本发明与已有的诸多肌电信号的消噪算法相比,具有如下特点:
本发明将肌电信号与脑电信号同时加入复杂网络,构建皮层肌肉功能网络,应用复杂网络理论的数学方法对构建网络进行研究,通过与规则网络、随机网络的对比验证了所构建网络的小世界特性,并提取网络特征成功对不同行为动作进行分类,这说明大脑皮层与肌肉组织之间确实存在信息交互,并且以所构建的网络探索更加高效的基于脑肌电信号的人体行为感知方法是可行的。
附图说明
图1为本发明的实施流程图;
图2为本发明实施例中不同行为动作不同频段构建的皮层肌肉功能网络邻接矩阵图;
图3为本发明实施例中加权聚类系数与加权特征路径长度随阈值的变化曲线;
图4为本发明实施例中网络模型图;
图5为本发明实施例中规则网络、应用阈值法简化得到的脑肌功能网络、随机网络的对比;
图6为本发明实施例中特征散点图(左为阈值法所得,右为固定加权边法所得)。
具体实施方式
下面结合附图对本发明的实施例作详细说明:本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,本实施例包括如下步骤:
步骤一,首先选取与屈腕、伸腕、握拳三个行为动作紧密相关的8通道脑电信号与4通道肌电信号并进行预处理,8通道脑电信号为C1、C2、C3、C4、C5、FC1、FC3、CP1,4通道肌电信号采集自尺侧伸腕肌、桡侧屈腕肌、屈指浅肌、伸指肌四块肌肉。
计算脑电信号与肌电信号两两通道之间的小波相干值,以此作为复杂网络的加权边值,得到相应的邻接矩阵。分别归一化EEG-EEG、EMG-EMG、EEG-EMG三个部分的小波相干值,得到归一化邻接矩阵。图2为不同行为动作不同频段构建的皮层肌肉功能网络邻接矩阵图。
步骤二,运用阈值法去除归一化邻接矩阵中较弱的加权边,即将邻接矩阵中较弱的边置为零。由此得到的邻接矩阵即为简化后的脑肌功能网络。如图3为加权聚类系数与加权特征路径长度随阈值的变化曲线。图4为构建脑肌功能网络后的网络模型示意图。
步骤三,在得到了简化后的脑肌功能网络之后,验证其小世界特性。将脑肌功能网络、规则网络、随机网络的平均加权聚类系数、平均加权特征路径进行比较。小世界网络具有较高的聚类系数和较小的特征路径长度。图5为规则网络、应用阈值法简化得到的脑肌功能网络、随机网络的对比。
步骤四,将脑肌功能网络的加权聚类系数和加权特征路径长度作为特征来识别屈腕、伸腕、握拳这三种行为动作。表1给出了应用Fisher线性判别作为分类器时不同简化网络模型方法下每两个动作的识别率。可以发现,应用经阈值法简化的皮层肌肉功能网络进行动作分类具有较高的精度。图6为特征散点图。
表1应用Fisher线性判别作为分类器时不同简化模型方法下的识别率
Claims (2)
1.一种基于小波相干的脑肌功能网络构建方法,其特征在于,该方法包括如下步骤:
步骤(1).首先选取与屈腕、伸腕、握拳三个行为动作紧密相关的8通道脑电信号与4通道肌电信号并进行预处理,8通道脑电信号为C1、C2、C3、C4、C5、FC1、FC3、CP1,4通道肌电信号分别采集自尺侧伸腕肌、桡侧屈腕肌、屈指浅肌、伸指肌四块肌肉;
步骤(2).计算步骤(1)中脑电信号与肌电信号两两通道之间的小波相干值,以此作为复杂网络的加权边值,得到相应的邻接矩阵;
步骤(3).将步骤(2)得到的邻接矩阵分为三部分,脑电与脑电信号之间、肌电与肌电信号之间、脑电与肌电信号之间,将三个部分的小波相干值分别归一化,得到归一化邻接矩阵;
步骤(4).运用阈值法去除归一化邻接矩阵中较弱的加权边,即将邻接矩阵中较弱的边置为零;由此得到的邻接矩阵即为简化后的脑肌功能网络。
2.根据权利要求1所述的一种基于小波相干的脑肌功能网络构建方法,其特征在于,所述的小波相干值计算方法具体如下:
1)定义信号x的CWT的范数平方为小波能量,即
WPx(t,f)=||CWTx(t,f)||2
WP是时间t与小波中心频率f的函数,CWT为连续小波变换;
2)信号x与y的交叉小波变换定义如下:
3)与传统的相干性估计类似,将小波相干值进行平滑;基于时间轴的平滑函数定义如下:
λ=t/a,c1为归一化常数,∧为卷积运算;基于尺度轴的平滑函数定义如下:
Sscale(CWTx(t,f))=CWTx(t,f)∧c2Π(0.6a)
c2为归一化常数,Π为矩形函数;矩形函数的长度是由尺度去相关决定的;
4)小波相干的定义如下:
a-1为尺度的逆,用于归一化XWT;Schwartz不等式确保WCxy∈[0,1];平滑过程S定义如下:
S(w)=Sscale[Stime(W)]。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010460960.1A CN111563581B (zh) | 2020-05-27 | 2020-05-27 | 一种基于小波相干的脑肌功能网络构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010460960.1A CN111563581B (zh) | 2020-05-27 | 2020-05-27 | 一种基于小波相干的脑肌功能网络构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111563581A true CN111563581A (zh) | 2020-08-21 |
CN111563581B CN111563581B (zh) | 2023-08-18 |
Family
ID=72069755
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010460960.1A Active CN111563581B (zh) | 2020-05-27 | 2020-05-27 | 一种基于小波相干的脑肌功能网络构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111563581B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112971786A (zh) * | 2021-02-05 | 2021-06-18 | 郑州大学 | 一种基于脑肌电信号小波相干系数的卒中康复评估方法 |
CN113558639A (zh) * | 2021-07-05 | 2021-10-29 | 杭州电子科技大学 | 基于格兰杰因果关系和图论的运动意图脑肌网络分析方法 |
CN113558638A (zh) * | 2021-07-05 | 2021-10-29 | 杭州电子科技大学 | 基于Vine Copula的脑肌耦合模型构建方法 |
CN113974652A (zh) * | 2021-10-26 | 2022-01-28 | 杭州电子科技大学 | 基于皮层肌肉功能网络模型的肌肉控制精准度确定方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722727A (zh) * | 2012-06-11 | 2012-10-10 | 杭州电子科技大学 | 基于脑功能网络邻接矩阵分解的脑电特征提取方法 |
CN107126193A (zh) * | 2017-04-20 | 2017-09-05 | 杭州电子科技大学 | 基于滞后阶数自适应选择的多变量因果关系分析方法 |
CN109497999A (zh) * | 2018-12-20 | 2019-03-22 | 杭州电子科技大学 | 基于Copula-GC的脑肌电信号时频耦合分析方法 |
CN109567799A (zh) * | 2018-12-26 | 2019-04-05 | 杭州电子科技大学 | 基于平滑小波相干性的肌电信号特征提取方法 |
CN109770924A (zh) * | 2019-01-24 | 2019-05-21 | 五邑大学 | 一种基于广义一致性构建脑功能网络与相关向量机的疲劳分类方法 |
CN110464343A (zh) * | 2019-08-16 | 2019-11-19 | 杭州电子科技大学 | 一种基于自主手部动作的增强型脑肌相干方法 |
CN110495880A (zh) * | 2019-08-16 | 2019-11-26 | 杭州电子科技大学 | 基于经颅电刺激脑肌耦合的运动障碍皮层可塑性管理方法 |
-
2020
- 2020-05-27 CN CN202010460960.1A patent/CN111563581B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722727A (zh) * | 2012-06-11 | 2012-10-10 | 杭州电子科技大学 | 基于脑功能网络邻接矩阵分解的脑电特征提取方法 |
CN107126193A (zh) * | 2017-04-20 | 2017-09-05 | 杭州电子科技大学 | 基于滞后阶数自适应选择的多变量因果关系分析方法 |
CN109497999A (zh) * | 2018-12-20 | 2019-03-22 | 杭州电子科技大学 | 基于Copula-GC的脑肌电信号时频耦合分析方法 |
CN109567799A (zh) * | 2018-12-26 | 2019-04-05 | 杭州电子科技大学 | 基于平滑小波相干性的肌电信号特征提取方法 |
CN109770924A (zh) * | 2019-01-24 | 2019-05-21 | 五邑大学 | 一种基于广义一致性构建脑功能网络与相关向量机的疲劳分类方法 |
CN110464343A (zh) * | 2019-08-16 | 2019-11-19 | 杭州电子科技大学 | 一种基于自主手部动作的增强型脑肌相干方法 |
CN110495880A (zh) * | 2019-08-16 | 2019-11-26 | 杭州电子科技大学 | 基于经颅电刺激脑肌耦合的运动障碍皮层可塑性管理方法 |
Non-Patent Citations (4)
Title |
---|
CAO CHANG-QING 等: "Application of telemetry EEG and EMG in conscious free moving dogs and monkeys for safety assessment and translational science" * |
F LOTTE 等: "A review of classification algorithms for EEG-based brain–computer interfaces: a 10 year update" * |
XIAOBI CHEN 等: "The cortical-muscular functional coupling performance evaluation based on harmonic wavelet and symbolic phase transfer entropy" * |
高云园 等: "基于相干性的多频段脑肌电信号双向耦合分析" * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112971786A (zh) * | 2021-02-05 | 2021-06-18 | 郑州大学 | 一种基于脑肌电信号小波相干系数的卒中康复评估方法 |
CN113558639A (zh) * | 2021-07-05 | 2021-10-29 | 杭州电子科技大学 | 基于格兰杰因果关系和图论的运动意图脑肌网络分析方法 |
CN113558638A (zh) * | 2021-07-05 | 2021-10-29 | 杭州电子科技大学 | 基于Vine Copula的脑肌耦合模型构建方法 |
CN113974652A (zh) * | 2021-10-26 | 2022-01-28 | 杭州电子科技大学 | 基于皮层肌肉功能网络模型的肌肉控制精准度确定方法 |
CN113974652B (zh) * | 2021-10-26 | 2024-04-26 | 杭州电子科技大学 | 基于皮层肌肉功能网络模型的肌肉控制精准度确定方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111563581B (zh) | 2023-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111563581B (zh) | 一种基于小波相干的脑肌功能网络构建方法 | |
CN110969108B (zh) | 一种基于自主运动想象脑电的肢体动作识别方法 | |
Wang et al. | Classification of surface EMG signals using optimal wavelet packet method based on Davies-Bouldin criterion | |
Xie et al. | Fuzzy approximate entropy analysis of chaotic and natural complex systems: detecting muscle fatigue using electromyography signals | |
Naik et al. | Applications of ICA and fractal dimension in sEMG signal processing for subtle movement analysis: a review | |
Dai et al. | Prediction of individual finger forces based on decoded motoneuron activities | |
CN111860410A (zh) | 基于多特征融合cnn的肌电手势识别方法 | |
Wang et al. | The Analysis of Surface EMG Signals with the Wavelet‐Based Correlation Dimension Method | |
CN112488002B (zh) | 基于n170的情绪识别方法及识别系统 | |
Khezri et al. | Neuro-fuzzy surface EMG pattern recognition for multifunctional hand prosthesis control | |
Powar et al. | A novel pre-processing procedure for enhanced feature extraction and characterization of electromyogram signals | |
CN113558639A (zh) | 基于格兰杰因果关系和图论的运动意图脑肌网络分析方法 | |
CN109674445B (zh) | 一种结合非负矩阵分解和复杂网络的肌间耦合分析方法 | |
Sun et al. | One-channel surface electromyography decomposition for muscle force estimation | |
Mahdavi et al. | Surface electromyography feature extraction based on wavelet transform | |
CN115640827B (zh) | 对电刺激数据处理的智能闭环反馈网络方法及系统 | |
Zhang et al. | Muscle force estimation based on neural drive information from individual motor units | |
CN116340824A (zh) | 基于卷积神经网络的肌电信号动作识别方法 | |
Yang et al. | Fast removal of ocular artifacts from electroencephalogram signals using spatial constraint independent component analysis based recursive least squares in brain-computer interface | |
Yu et al. | Wrist torque estimation via electromyographic motor unit decomposition and image reconstruction | |
Khodadadi et al. | Nonlinear analysis of biceps surface EMG signals for chaotic approaches | |
Wang et al. | Deep convolutional neural network for decoding EMG for human computer interaction | |
CN113558638A (zh) | 基于Vine Copula的脑肌耦合模型构建方法 | |
Arjunan et al. | Fractal properties of surface electromyogram for classification of low-level hand movements from single-channel forearm muscle activity | |
CN113128384A (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 |