基于运动想象的混合脑机接口特征提取方法与流程

1.本发明属于脑电肌电信号特征分析领域，本发明公开了一种基于运动想象的混合脑机接口特征提取方法。


背景技术：

2.由于脑血管疾病的高发病率和频繁发生的各种意外事故，中国的瘫痪患者数量与日俱增。针对瘫痪患者，在经过诸如外科手术等急性期的临床处理后，康复成为主要的一种治疗手段，它能够帮助患肢恢复运动功能，重新学习日常生活活动，其中手部动作康复训练过程尤为重要。本发明方法帮助病人康复训练，首次采用脑机协同的理论对动作意图进行识别，识别率相比单一脑电特征和肌电特征更高，泛化性较强。采取选取多特征融合互补的方式，可以提高动作意图识别的准确性、稳定性以及模型的鲁棒性。


技术实现要素：

3.本发明针对现有技术的不足，提出了一种基于运动想象的混合脑机接口特征提取方法。
4.本发明与单一脑电信号和肌电信号相比，具有如下优点和有益效果：本文利用肌肉协同之间协同元的相似性与不同性来进行特征提取；本文利用脑功能网络和肌肉协同这两种方法来进一步得出能综合反应受试者基于hbci的互补融合特征；
5.为实现以上目的，本发明方法主要包括以下步骤：
6.步骤(1)采集脑电数据。实验样本频率设置为5000hz。根据实验要求，受试者按照随机动作顺序进行曲腕、曲肘和握拳动作。
7.步骤(2)信号预处理。脑电信号质量的保证是准确有效研究的基础。在脑电信号采集过程中，不可避免地会受到各种信号和噪声的干扰。由于eeg信号非常微弱，因此很容易受到许多干扰源的干扰，包括眼球运动和眨眼伪影，心电图和工频干扰。因此，为了保证基于脑电信号的研究结果的可靠性，有必要对记录的信号进行相应的有效降噪步骤，以尽可能减少噪声对研究结果的干扰。
8.步骤(3)提取拓扑特征，脑电信号基于似然同步理论构建脑功能网络，选取合适阈值进行二值化处理，提取二值化脑功能网络的拓扑特征，肌电信号采用非负矩阵分解，协同矩阵选取合适的数据变异度寻找不同的协同元，提取协同元的特征。将脑电特征和肌电特征进行特征融合。
9.利用似然同步构建脑功能网络程序实现的理论依据如下：假设有m个同一时间段记录的信号序列x
k,i
，其中k表示信道数量(i＝1,2,
…
,m)，i表示信号记录的时间点(i＝1,2,
…
,n)，对每一道信号利用时间延迟重建内嵌向量为：
10.x
k,i
＝(x
k,i
,x
k,i l
,x
k,i 2l
,
…
,x
k,i (m-1)l
)
ꢀꢀ
(1)
11.式中，l和m分别表示时间延迟和嵌入的维数。定义任意两个内嵌向量x
k,i
和x
k,j
之间的欧氏距离小于ε的概率为：
[0012][0013]
式中，|
·
|操作符是求两个向量之间的欧氏距离；θ(x)是heaviside阶梯函数，即满足
[0014][0015]
式(2)中w1和w2是两个窗，w1是泰勒修正参数，用来消除由于自相关特性对向量间距离过近造成的影响；w2窗用来提高时间分辨率，避免信号向量之间的时间间隔太长。当前的研究中，对于窗口w1和w2的取值需要满足条件：w1《《w2《《n(《《表示远小于)，对于eeg数据的分析一般取w1＝100，w2取时间序列长度的1/10，对式(1)中的时间延迟参数l和嵌入的维数m尚没有唯一的取值标准，本文按照研究经验取l＝10，m＝10。综上可知，式(2)中对于每个k和每个时间点i对应的距离参数ε
k,i
[0016]
由概率所决定。假定m道信号中与内嵌向量x
k,i
之间的距离小于ε
k,i
,并满足窗口条件w1＜|i-j|＜w2的内嵌向量x
k,j
有h
i,j
个，即：
[0017][0018]
由实际意义可知h
i,j
的值介于0和m之间，表示i和j时刻非常接近或者同步的内嵌向量的数量。综上信息，对于k道信号和另外m-1道信号i和j时刻的似然同步值s
i,j,k
定义如下:
[0019][0020]
最后对所有时间点j求平均，得到k道信号和其他m-1道信号在i时刻的似然同步值为：
[0021][0022]
提取协同元的特征，具体如下：用非负矩阵分解的方法对肌电信号进行分解。如公式(1)所示，vm×
t
为待分解的表面肌电信号包络矩阵，m表示肌肉的数目，t表示数据的长度；wm×n＝[w
i1
,w
i2
,
…
,w
in
]为协同结构矩阵，表示肌肉间的线性组合模式i＝1，2，
…
，m；w
in
幅值大小表示肌肉在组合模式中的贡献程度；cn×
t
＝[c
1j
,c
2j
,
…
,c
nj
]为激活系数矩阵，代表协同结构如何按时间被调制的神经信号，同时表示每一个肌肉协同对一块肌肉整体激发的贡献，j＝1,2，
…
，t；em×
t
为噪声信号；n代表肌肉协同数目；
[0023]vm
×
t
＝wm×
ncn
×
t
em×
t
ꢀꢀ
(7)
[0024]
要使||e||最小，即求如式(8)的最大似然解。
[0025]
{w,c}＝argmax
w,c
p(v|w,c)＝argmin
w,c
[-logp(v|w,c)]
ꢀꢀ
(8)
[0026]
假设噪声服从高斯分布，那么其目标函数等于如下的欧式距离函数：
[0027][0028]
通过传统的梯度法，得到以下乘性迭代规则：
[0029][0030][0031]
为了确定最小的肌肉协同数目，通常计算每一个n值下由分解得到的w和c矩阵重构的肌电包络信号vrm×
t
，然后根据vm×
t
与vrm×
t
生成的数据变异度vaf，来确定肌肉协同数目n；确定n值方法为：
①
、vaf的大于设置阈值，
②
相邻vaf的差值小于设定阈值；
[0032]
vaf的具体计算公式如下：
[0033]
vaf＝1-(vm×
t-vrm×
t
)2/vm×
t2
ꢀꢀ
(12)
[0034]
步骤(4)使用极限学习机分类；极限学习机它包含着三大部分，分别是输入层、隐藏层和输出层。当输入为一组拥有m类总数为n的训练样本时，一个隐藏层极点数为l、激活函数为g(x)的单隐藏层反馈神经网络可以表示为：
[0035][0036]
其中xj＝[x
j1
,x
j2
,
…
,x
jn
]
t
为输入，tj＝[t
j1
,t
j2
,
…
,t
jn
]
t
为对应的期望输出，bi为隐藏极点i的激活阈值，wi＝[w
i1
,w
i2
,
…
,w
in
]
t
将输入层与隐藏层联系在一起，βi＝[β
i1
,β
i2
,
…
,β
im
]
t
将隐藏层与输出层联系在一起，oj是单隐藏层反馈神经网络的实际输出。单隐藏层反馈神经网络能够零误差拟合训练样本等价于因为elm可以随机生成隐藏层参数{wi,bi}。
[0037]
此时，将训练的标签值带入到神经网络的输出层中，得到公式
[0038]
t＝hβ (14)
[0039]
式中
[0040][0041]
通常情况下隐藏神经元的个数l是远远小于训练样本的个数n的，即l《《n，所以我们通过求最小值来确定输出权值，公式为：
[0042]
β＝arg min||t-hβ||
ꢀꢀ
(16)
[0043]
由极小范数解准则可知，式(16)可以由以下的公式求解得到：
[0044][0045]
其中h

是隐藏层矩阵h的moore-penrose增广逆，简称伪逆。对于分类问题，elm一般采用单对全标签编码将多分类问题转化为多输出函数回归。计算后查询样本的实际输
出，然后由实际输出的索引给出预测的类标签。
[0046]
本发明与现有技术相比，具有如下优点和有益效果：
[0047]
本发明利用肌肉协同之间协同元的相似性与不同性来进行特征提取；通过脑功能网络和肌肉协同这两种方法来进一步得出能综合反应受试者基于混合脑机接口的互补融合特征；相比单一脑电信号和肌电信号，动作识别率有显著提升。
附图说明
[0048]
图1为本发明的实施流程图；
[0049]
图2为本发明所进行实验中使用的电极分布图；
[0050]
图3为预处理脑电信号前后对比；
[0051]
图4为受试者1屈腕、屈肘、握拳生成的脑功能网络和对应二值化网络；
[0052]
图5为受试者1屈腕、屈肘对应的协同元分布比例。
具体实施方式
[0053]
下面结合附图对本发明的实施例作详细说明：本实施例在以本发明技术方案为前提下进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
[0054]
本发明的实施例采集了八位志愿者分别用左右手执行腕曲，肘曲和握拳共计6种动作的eeg信号。在数据预处理阶段使用了wt，脑电信号采用基于manhattan距离的sl构造加权fbn。二进制fbn由阈值和mst两种方法构建。肌电信号采用非负矩阵分解的方法构建肌肉协同模型。最后将两者特征融合，使用极限学习机对该动作进行动作识别。
[0055]
本文将介绍参与者、实验任务、数据集的信息与这项工作所涉及的理论，介绍方法的配置，并说明所提出方案的性能。最后与已有的相关研究进行比较，提出预期的工作。
[0056]
如图1所示，本产品包括如下步骤：
[0057]
步骤(1)采集脑电数据。实验样本频率设置为5000hz。根据实验要求，受试者按照随机动作顺序进行曲腕、曲肘和握拳动作。
[0058]
步骤(2)信号预处理。脑电肌电信号质量的保证是准确有效研究的基础。图2为本发明所进行实验中使用的电极分布图，在脑电肌电信号采集过程中，不可避免地会受到各种信号和噪声的干扰。由于eeg信号非常微弱，因此很容易受到许多干扰源的干扰，包括眼球运动和眨眼伪影，心电图和工频干扰。emg信号会受到肌肉疲劳干扰，因此，为了保证基于脑电信号的研究结果的可靠性，有必要对记录的信号进行相应的有效降噪步骤，以尽可能减少噪声对研究结果的干扰。对记录的信号进行wt去噪，以尽可能减少噪声对研究结果的干扰。如图3为志愿者1左手wf的去噪肌电信号。
[0059]
步骤(3)提取拓扑特征，脑电信号基于似然同步理论构建脑功能网络，选取合适阈值进行二值化处理，提取二值化脑功能网络的拓扑特征，图4为受试者1屈腕、屈肘、握拳生成的脑功能网络和对应二值化网络，图5为受试者1屈腕、屈肘对应的协同元分布比例，肌电信号采用非负矩阵分解，协同矩阵选取合适的数据变异度寻找不同的协同元，提取协同元的特征。将脑电特征和肌电特征进行特征融合。
[0060]
步骤(4)使用极限学习机对动作分类；
[0061]
最后8位志愿者的分类准确性通过表1单一脑电信号特征、表2单一肌电信号特征、表3融合特征的识别率所示；脑功能网络的拓扑特征比均方根和肌肉协同特征识别率高，而部分单一特征虽然相比其他单一特征识别率较低，但通过互补特征融合的方式弥补了自身缺陷。脑机协同理论主要是大脑皮层通过脊髓和周围神经控制肌肉组织的运动，使得肢体可以完成一定的运动功能，而肢体的运动又可以通过传入神经影响着大脑皮层的活动，两者相互协同作用。脑功能网络的拓扑特征与肌肉协同的特征融合有效的提高了动作意图识别的识别率和鲁棒性，且具有良好的泛化能力；
[0062]
表1
[0063][0064]
表2
[0065][0066]
表3
[0067]