CN114159081A - 基于时变混合Copula互信息的肌电耦合方法 - Google Patents
基于时变混合Copula互信息的肌电耦合方法 Download PDFInfo
- Publication number
- CN114159081A CN114159081A CN202111522190.XA CN202111522190A CN114159081A CN 114159081 A CN114159081 A CN 114159081A CN 202111522190 A CN202111522190 A CN 202111522190A CN 114159081 A CN114159081 A CN 114159081A
- Authority
- CN
- China
- Prior art keywords
- time
- varying
- copula
- mutual information
- function
- 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.)
- Pending
Links
- 241000039077 Copula Species 0.000 title claims abstract description 108
- 238000010168 coupling process Methods 0.000 title claims abstract description 37
- 230000008878 coupling Effects 0.000 claims abstract description 21
- 238000005859 coupling reaction Methods 0.000 claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 14
- 238000005315 distribution function Methods 0.000 claims abstract description 13
- 230000006870 function Effects 0.000 claims description 53
- 238000004364 calculation method Methods 0.000 claims description 13
- 230000003183 myoelectrical effect Effects 0.000 claims description 11
- 239000000126 substance Substances 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 241000932075 Priacanthus hamrur Species 0.000 claims description 2
- 210000003205 muscle Anatomy 0.000 abstract description 18
- 230000007246 mechanism Effects 0.000 abstract description 6
- 238000007781 pre-processing Methods 0.000 abstract 1
- 230000008569 process Effects 0.000 description 6
- 210000003169 central nervous system Anatomy 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 210000002310 elbow joint Anatomy 0.000 description 3
- 210000002161 motor neuron Anatomy 0.000 description 3
- 210000000653 nervous system Anatomy 0.000 description 3
- 210000001364 upper extremity Anatomy 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000008602 contraction Effects 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- 230000001054 cortical effect Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007659 motor function Effects 0.000 description 2
- 210000001087 myotubule Anatomy 0.000 description 2
- 210000005036 nerve Anatomy 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 210000000225 synapse Anatomy 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 208000001840 Dandruff Diseases 0.000 description 1
- 208000012661 Dyskinesia Diseases 0.000 description 1
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 1
- 206010060820 Joint injury Diseases 0.000 description 1
- 206010061296 Motor dysfunction Diseases 0.000 description 1
- 206010049565 Muscle fatigue Diseases 0.000 description 1
- 206010056720 Muscle mass Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 210000000852 deltoid muscle Anatomy 0.000 description 1
- 230000004064 dysfunction Effects 0.000 description 1
- 239000012636 effector Substances 0.000 description 1
- 238000002567 electromyography Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 239000004519 grease Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 210000001611 motor endplate Anatomy 0.000 description 1
- 230000004118 muscle contraction Effects 0.000 description 1
- 230000004220 muscle function Effects 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 230000002232 neuromuscular Effects 0.000 description 1
- 210000000715 neuromuscular junction Anatomy 0.000 description 1
- 210000002569 neuron Anatomy 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 230000009466 transformation Effects 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/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/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]
- A61B5/397—Analysis of electromyograms
-
- 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/725—Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Signal Processing (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Psychiatry (AREA)
- Physiology (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
技术领域
本发明属于神经系统运动控制机制研究领域,具体涉及基于时变混合Copula互信息的肌电耦合方法。
背景技术
肌电(electromyography,EMG)信号是神经肌肉活动产生的信号,肌肉在收缩过程中所产生的生理电信号,是人体运动控制系统中运动执行末端的信息表征。起主导作用的神经系统是通过运动神经元(motor neuron,MN)连接,每个运动神经元细胞的突触分出数十根分支直达运动终板效应器与肌肉纤维耦合,形成一对一的神经肌肉接头支配多条肌纤维收缩,它们整合在一起形成了运动单元(motor units,MU)。人体肌肉中包含了难以计数的运动单元,当肌肉运动时,大脑运动皮层区就会产生刺激信号或运动指令向下传导,中枢神经系统中的运动神经元发出控制肌肉收缩的信号到相对应的肌肉块。通过表面电极测量所得的肌电信号称为表面肌电信号(surface EMG,sEMG),它作为无创采集的人体生理电信号,是浅层肌肉以及神经干上电活动在皮肤表面的综合表现,其中蕴含着丰富的信息。在运动过程中,肌肉之间是彼此联系和相互作用的,并且受到运动神经系统的调控,肌电信号间的耦合可以反映肌肉的运动功能状态和不同的运动范式下中枢神经系统的不同控制策略,探索内在的运动功能控制机制,为理解运动障碍的病理机制提供了理论基础。
Conway于1995年发现了在运动过程中脑电信号与相关的肌电之间的皮层肌肉相干性,随后,Brown等发现在运动过程中多通道的肌电信号之间也具有一定的关系,研究者开始采用不同的相关性分析方法获取肌电信号之间的功能联系,进一步探索运动神经控制机理。Nguyen等发现不同运动范式下的肌电信号之间的耦合关系有所差异。
近年来,相干性、互信息等方法被广泛应用到肌电耦合分析中。杜义浩等利用相干性分析,发现运动致疲劳过程中协同肌受中枢神经系统控制以更加同步的方式活动;陈玲玲等利用互信息,分析上肢各肌电通道间的关联特性,构建肌肉功能网络。相干性耦合方法只能描述线性的耦合关系,互信息在研究中应用于非线性耦合分析。然而,互信息的估计严重依赖联合概率密度的精确表示。而通过Copula函数估计互信息,不仅能有效避免对联合密度函数的估计,而且计算复杂度低。Copula函数种类繁多,不同的Copula函数所得到的Copula互信息估计值不相同,吴亚婷等利用Copula互信息进行肌间耦合网络分析,发现肱三头肌与三角肌中束、三角肌后束的耦合强度相对较高,肌间功能联系紧密;而肱二头肌在该运动下独立于其他肌肉。Copula互信息应用的是参数静态的常相关Copula函数对互信息进行估计,这种静态Copula函数不能实时反映相关结构的动态变化,具有一定的局限性,少数应用参数是动态的Copula函数,但多数采用单个Copula函数,考虑到实际情况中的复杂结构,仅用一种Copula函数无法与数据分布类型很好地匹配,可能会出现失真的情况,不能很好地描述各信号间的功能耦合关系。
发明内容
针对现有技术的不足,本发明提出了基于时变混合Copula互信息的肌电耦合方法,结合多种时变Copula函数构建时变混合Copula模型,计算肌电信号的耦合特性。
基于时变混合Copula互信息的肌电耦合方法,具体包括以下步骤:
步骤一:同步采集多通道sEMG信号,经过去噪处理后作为计算样本。
作为优选,所述去噪处理为:对sEMG信号进行去均值、去基线漂移,再利用无限冲激响应陷波滤波器抑制50Hz工频干扰,然后用带通滤波器对其进行0.5~200HZ的带通滤波。
步骤二:选择步骤一去除噪声后的两通道计算样本X={X1,X2,...Xn}和Y={Y1,Y2,...Yn},X、Y分别为来自连续分布函数u和v的同分布样式,n为计算样本的长度,使用计算样本的经验分布函数估计,得到对应的估计值序列和
步骤三:选择多个时变二元Copula函数,构建得到如下时变混合Copula模型:
其中,s为时变混合Copula模型中时变二元Copula函数的数量,Cj(u,v;θjt)为第j个时变二元Copula函数的分布函数,λjt、θjt分别为第j个时变二元Copula函数的权重和时变参数,且∑λjt=1,λjt≥0。θcopula为时变混合Copula模型的时变参数。
对时变混合Copula模型求偏导,得到时变Copula密度函数为:
其中,cj(u,v;θjt)为第j个时变二元Copula密度函数。
步骤四:将步骤二得到的计算样本的估计值序列和代入对数似然函数中进行求解,然后利用期望最大化算法和拟牛顿算法估计得到时变混合Copula模型的估计时变参数定义时变混合Copula互信息MI(X,Y)为:
步骤五:重复步骤二~四,计算每两个通道计算样本之间的时变混合Copula互信息,得到肌电信号间的耦合强度。
本发明具有以下有益效果:
采用多个不同的时变Copula函数构建时变混合Copula模型,定量分析了运动过程中肌电信号间的耦合特性,可以很好地描述特定运动下的肌电信号之间的相关性,反映特定运动下的相关肌肉块之间的相互作用以及中枢神经系统的控制策略,弥补了传统的肌电耦合方法并不能有效地描述sEMG信号的耦合强度大小的缺点,为进一步探究运动功能障碍的病理机制提供一定的理论依据。
附图说明
图1为多通道肌电耦合方法流程图;
图2为多通道sEMG信号同步采集示意图;
图3为经过去噪处理
具体实施方式
以下结合附图对本发明作进一步的解释说明;
如图1所示,基于时变混合Copula互信息的肌电耦合方法,具体包括以下步骤:
步骤一:本实施例招募了5名健康受试者,均为无上肢运动功能障碍或关节损伤历史的健康人,年龄为21~25岁,体重为50~76kg,身高为165~175cm,所有受试者均为右利手。如图2所示,采用美国Delsys公司生产的TrignoTM Wireless sEMG信号采集装置,在2000NHz的采样频率下同步采集受试者与肘关节屈曲动作密切相关的肱桡肌(Brachioradialis,B)、肱二头肌(Biceps Brachii,BB)、肱三头肌(TricepsBrachii,TB)的3通道sEMG信号。首先用医用酒精擦拭采集部位,去除皮肤表面的油脂和皮屑。然后由受试者手持哑铃根据节拍器匀速地进行肘关节运动,运动10s后,为防止肌肉疲劳,休息5s,每位受试者进行三次实验。对采集得到的sEMG信号进行去均值、去基线漂移,再利用无限冲激响应陷波滤波器抑制50Hz工频干扰,然后用带通滤波器对其进行0.5~200HZ的带通滤波。经过去噪处理后的计算样本频域图如图3所示。
步骤二:选择步骤一去除噪声后的两通道计算样本X={X1,X2,...Xn}和Y={Y1,Y2,...Yn},X、Y分别为来自连续分布函数u和v的同分布样式,n为计算样本的长度,使用计算样本的经验分布函数估计,得到对应的估计值序列和
其中,I为示性函数。
步骤三:本实施例中选择了4个时变二元Copula函数,分别为时变二元高斯Copula函数、时变二元Gumbel Copula函数、时变二元Clayton Copula函数和时变SJC Copula函数。
(1)时变二元高斯Copula函数的分布函数为:
其中,Φ-1(·)为标准正态分布的逆函数,ρt为高斯时变相关系数,演化方程为:
其中,Λ(·)为Logistic变换函数,Λ(x)=(1-e-x)(1+e-x)-1,保证ρt∈(-1,1),ωn、βn和αn为时变方程参数,时变高斯Copula函数无法捕捉尾部相关变化。
(2)时变二元Gumbel Copula函数的分布函数为:
其中,θ1t为Gumbel时变相关系数,演化方程为:
其中,Λ(x)=1+e-x,保证θ1t∈(0,1),时变Gumbel Copula函数是非对称的,可以较好地捕捉到上尾部相关变化。
(3)时变二元Clayton Copula函数的分布函数为:
其中,θ2t为Clayton时变相关系数,演化方程为
其中,Λ(x)=e-x,保证θ2t∈(0,1),时变Clayton Copula函数是非对称的,可以较好地捕捉到下尾部相关变化。
(4)时变Symmetrized Joe-Clayton(SJC)Copula函数的分布函数为:
其中:
将上述4个时变二元Copula函数混合,构建时变混合Copula模型:
其中,λ1t,λ2t,λ3t,λ4t分别表示4个时变二元Copula函数的权重,λ1t+λ2t+λ3t+λ4t=1,θcopula为时变混合Copula模型的时变参数。
对时变混合Copula模型求偏导,得到时变Copula密度函数为:
通过添加一个惩罚项进行求解:
其中,γn为光滑参数,控制模型复杂度,为SCAD惩罚函数,再利用期望最大化算法和拟牛顿算法估计时变混合Copula模型的参数。然后利用期望最大化算法和拟牛顿算法估计得到时变混合Copula模型的估计时变参数
采用Akaike信息准则(AIC)衡量构建的时变混合Copula模型对肌间耦合关系的拟合优度:
AIC=-2ln(L)+2k (15)
其中,k为模型中参数的个数,L为似然函数值,AIC的值越小,模型拟合程度越好。由于运动过程中肌电信号间的耦合突出表现在β(15-30Hz)与γ(30-60Hz)频段,采用切比雪夫Ⅱ型带通滤波器提取这两个频段的肌电信号。表1给出了第1位受试者第1次肘屈时的时变Copula、混合Copula、时变混合Copula在β、γ频段的AIC指标:
表1
由表1可知,时变混合Copula的AIC值比单一时变Copula与混合Copula的AIC值低,这说明在第1位受试者的第1次上肢肘关节运动中,时变混合Copula对肌间结构的拟合要比单一时变Copula、混合Copula更好。
为比较时变Copula、混合Copula、时变混合Copula的性能,统计所有受试者根据AIC最小原则在每次肘屈动作、每个特征频段上所选择的Copula函数类型,如表2所示:
表2
由表2可知,大多数情况下都选择了时变混合Copula,少部分选择了混合Copula。这说明时变混合Copula函数更适合描述肌电之间的耦合关系。
步骤五:根据互信息与Copula熵的等价关系,定义时变混合Copula互信息MI(X,Y)为:
其中,为时变混合Copula熵,E[·]表示求取期望,MI(X,Y)为计算样本X和Y的时变混合Copula互信息。根据公式计算每两个通道计算样本之间的时变混合Copula互信息,得到肌电信号间的耦合强度。
以上所述的实施例仅是对本发明的优选实例方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。
Claims (6)
1.基于时变混合Copula互信息的肌电耦合方法,其特征在于:具体包括以下步骤:
步骤一:同步采集多通道sEMG信号,经过去噪处理后作为计算样本;
步骤二:选择去噪处理后的两通道计算样本X={X1,X2,...Xn}和Y={Y1,Y2,...Yn},X、Y分别为来自连续分布函数u和v的同分布样式,n为计算样本的长度,使用计算样本的经验分布函数估计,得到对应的估计值序列和
步骤三:选择多个时变二元Copula函数,构建得到如下时变混合Copula模型:
其中,s为时变混合Copula模型中时变二元Copula函数的数量,Cj(u,v;θjt)、λjt、θjt分别为第j个时变二元Copula函数的分布函数、权重、时变参数,∑λjt=1,λjt≥0;θcopula为时变混合Copula模型的时变参数;
对时变混合Copula模型求偏导,得到时变Copula密度函数为:
其中,cj(u,v;θjt)为第j个时变二元Copula密度函数;
步骤四:将步骤二得到的计算样本的估计值序列和代入对数似然函数中进行求解,然后利用期望最大化算法和拟牛顿算法估计得到时变混合Copula模型的估计时变参数定义时变混合Copula互信息MI(X,Y)为:
步骤五:重复步骤二~四,计算每两个通道计算样本之间的时变混合Copula互信息,得到肌电信号间的耦合强度。
2.如权利要求1所述基于时变混合Copula互信息的肌电耦合方法,其特征在于:所述去噪处理为:对sEMG信号进行去均值、去基线漂移,再利用无限冲激响应陷波滤波器抑制50Hz工频干扰,然后用带通滤波器对其进行0.5~200HZ的带通滤波。
4.如权利要求1所述基于时变混合Copula互信息的肌电耦合方法,其特征在于:步骤三中,选择时变二元高斯Copula函数、时变二元Gumbel Copula函数、时变二元ClaytonCopula函数和时变SJC Copula函数构建时变混合Copula模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111522190.XA CN114159081A (zh) | 2021-12-13 | 2021-12-13 | 基于时变混合Copula互信息的肌电耦合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111522190.XA CN114159081A (zh) | 2021-12-13 | 2021-12-13 | 基于时变混合Copula互信息的肌电耦合方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114159081A true CN114159081A (zh) | 2022-03-11 |
Family
ID=80486140
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111522190.XA Pending CN114159081A (zh) | 2021-12-13 | 2021-12-13 | 基于时变混合Copula互信息的肌电耦合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114159081A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060224054A1 (en) * | 2005-03-30 | 2006-10-05 | Kabushiki Kaisha Toshiba | Pulse wave detecting device and method therefor |
CN109497999A (zh) * | 2018-12-20 | 2019-03-22 | 杭州电子科技大学 | 基于Copula-GC的脑肌电信号时频耦合分析方法 |
CN111931606A (zh) * | 2020-07-23 | 2020-11-13 | 杭州电子科技大学 | 基于混合Copula互信息的肌间耦合分析方法 |
CN111931129A (zh) * | 2020-07-23 | 2020-11-13 | 杭州电子科技大学 | 基于高斯Copula传递熵的肌间耦合网络分析方法 |
CN112130668A (zh) * | 2020-09-27 | 2020-12-25 | 杭州电子科技大学 | 一种R藤Copula互信息的肌间耦合分析方法 |
CN112232301A (zh) * | 2020-11-16 | 2021-01-15 | 杭州电子科技大学 | 基于多尺度Copula互信息的肌间耦合网络分析方法 |
CN112509689A (zh) * | 2021-02-08 | 2021-03-16 | 杭州电子科技大学 | 基于时变Copula互信息的肌间耦合分析方法 |
-
2021
- 2021-12-13 CN CN202111522190.XA patent/CN114159081A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060224054A1 (en) * | 2005-03-30 | 2006-10-05 | Kabushiki Kaisha Toshiba | Pulse wave detecting device and method therefor |
CN109497999A (zh) * | 2018-12-20 | 2019-03-22 | 杭州电子科技大学 | 基于Copula-GC的脑肌电信号时频耦合分析方法 |
CN111931606A (zh) * | 2020-07-23 | 2020-11-13 | 杭州电子科技大学 | 基于混合Copula互信息的肌间耦合分析方法 |
CN111931129A (zh) * | 2020-07-23 | 2020-11-13 | 杭州电子科技大学 | 基于高斯Copula传递熵的肌间耦合网络分析方法 |
CN112130668A (zh) * | 2020-09-27 | 2020-12-25 | 杭州电子科技大学 | 一种R藤Copula互信息的肌间耦合分析方法 |
CN112232301A (zh) * | 2020-11-16 | 2021-01-15 | 杭州电子科技大学 | 基于多尺度Copula互信息的肌间耦合网络分析方法 |
CN112509689A (zh) * | 2021-02-08 | 2021-03-16 | 杭州电子科技大学 | 基于时变Copula互信息的肌间耦合分析方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102835955B (zh) | 一种无需设定阈值的脑电信号中眼电伪迹自动去除的方法 | |
CN106709469B (zh) | 基于脑电和肌电多特征的自动睡眠分期方法 | |
CN109674445B (zh) | 一种结合非负矩阵分解和复杂网络的肌间耦合分析方法 | |
CN111931606B (zh) | 基于混合Copula互信息的肌间耦合分析方法 | |
CN102697493A (zh) | 一种快速的脑电信号中眼电伪迹自动识别和去除的方法 | |
CN113229831B (zh) | 基于肌电和肌氧信号的运动功能监测管理方法 | |
CN106236080B (zh) | 基于多通道的脑电信号中肌电噪声的消除方法 | |
CN112130668B (zh) | 一种R藤Copula互信息的肌间耦合分析方法 | |
CN103258120A (zh) | 一种基于脑电信号的中风康复程度指标计算方法 | |
CN112232301A (zh) | 基于多尺度Copula互信息的肌间耦合网络分析方法 | |
CN111708978B (zh) | 多尺度时频肌间耦合分析方法 | |
Li et al. | ECG denoising method based on an improved VMD algorithm | |
CN110638444A (zh) | 基于MEMD-rTVgPDC的皮层肌肉耦合分析方法 | |
CN112541415B (zh) | 基于符号传递熵和图论的脑肌功能网络运动疲劳检测方法 | |
Wu et al. | EMGdi signal enhancement based on ICA decomposition and wavelet transform | |
CN112509689A (zh) | 基于时变Copula互信息的肌间耦合分析方法 | |
CN113723557A (zh) | 一种基于多频带时空卷积网络的抑郁脑电分类系统 | |
Islam et al. | Probability mapping based artifact detection and wavelet denoising based artifact removal from scalp EEG for BCI applications | |
CN117281479A (zh) | 一种基于表面肌电信号多维度特征融合的人体下肢慢性疼痛判别方法、存储介质及装置 | |
Guo et al. | Analysis and feature extraction of EEG signals induced by anesthesia monitoring based on wavelet transform | |
CN114159081A (zh) | 基于时变混合Copula互信息的肌电耦合方法 | |
CN116584960A (zh) | 一种膈肌肌电信号降噪方法 | |
CN116172576A (zh) | 一种基于多模块神经网络的脑电信号伪迹去除方法 | |
Altamirano | EMG pattern prediction for upper limb movements based on wavelet and Hilbert-Huang transform | |
El Khadiri et al. | A Comparison of the Denoising Performance Using Capon Time-Frequency and Empirical Wavelet Transform Applied on Biomedical Signal. |
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 |