[发明专利]基于Copula-GC的脑肌电信号时频耦合分析方法

权利要求书

1.基于Copula-GC的脑肌电信号时频耦合分析方法，其特征在于：该方法包括如下步骤：

步骤一，获取人体上肢动作时的头皮脑电信号和表面肌电信号样本数据；

步骤二，脑肌电信号预处理；

步骤三，基于Copula-GC的时域分析；

具体为：EEG与EMG信号之间的Copula-GC值大小反映大脑皮层和运动神经元之间信息交换的数量；对于预处理后的EEG信号X＝{x_t}与EMG信号Y＝{y_t}之间的cGC通过下式计算得到：

其中，E表示样本空间的期望；F和G分别表示Y和X的条件边缘分布；c表示条件Copula的密度函数；和分别为X和Y的历史信息，m和n为变量X和Y的滞后阶次；利用贝叶斯信息准则进行回归分析确定m和n；EEG→EMG的cGC值代表大脑皮层传递到相应肌肉的信息量，EMG→EEG的cGC值代表肌肉细胞反馈回大脑皮层的信息量；

步骤四，基于子带分解和Copula-GC的频域耦合强度统计；

具体为：利用FIR滤波器将1-50Hz的EEG和EMG信号分隔为49个间隔为1Hz的子带信号，然后利用式(1)计算相同子带的EEG和EMG信号之间的cGC值，记为cGC_X→Y(f)和cGC_Y→X(f)，其中f代表子带频率；cGC值越大，说明在该频段内信息的传递量越大；

为进一步定量刻画EEG-EMG信号在不同频段的双向耦合特征差异，基于EEG-EMG同步研究中一致性分析的相干面积指标，定义EEG-EMG的耦合强度指标CS_X→Y和CS_Y→X分别如式(2)和(3)所示：

在式(2)和(3)中，Δf表示子带分辨率，cGC_X→Y(f)和cGC_Y→X(f)分别表示频率f处的Copula-GC值；CS的数值越大，说明该方向上的功能耦合越强；

步骤五，采用步骤三和步骤四的计算结果从时频域上分析脑肌电间的双向功能耦合影响。

2.根据权利要求1所述的基于Copula-GC的脑肌电信号时频耦合分析方法，其特征在于：步骤一中，脑肌电信号采集时，使用采集系统来获取位于初级运动区域的C3/C4通道的EEG信号，以及指浅屈肌和桡腕屈肌的EMG信号；在放置电极之前，洗净头皮，确保阻抗处于低水平，并用酒精清洁皮肤表面；EEG和EMG信号的采样频率设置为1000Hz。

3.根据权利要求1所述的基于Copula-GC的脑肌电信号时频耦合分析方法，其特征在于：步骤二的脑肌电信号预处理；

具体为：脑电信号与肌电信号采集过程中利用BrainAmp DC设备进行了工频干扰、基线漂移、溢出和眼动伪迹的去除，然后利用独立成分分析法对EEG中的眼电伪迹进行了去除，再分别对EEG和EMG信号进行小波软阈值消噪处理。