首页 / 专利库 / 生物工程 / 脑机接口 / 便携式脑电采集设备及其在SSVEP和运动想象中的应用

便携式脑电采集设备及其在SSVEP和运动想象中的应用

阅读:597发布:2020-12-03

专利汇可以提供便携式脑电采集设备及其在SSVEP和运动想象中的应用专利检索,专利查询,专利分析的服务。并且一种便携式脑电采集设备及其在SSVEP和运动想象中的应用,设备包括:用于提供电源的系统供电 电路 ,及依次连接的脑 电极 帽转接线、PGA放大电路、AD转换器、STM32处理器和WIFI模 块 。将 压缩 感知 理论嵌入便携式脑电采集设备当中,在保证数据 精度 的前提下,很好地提升了数据的传输效率。在SSVEP和运动想象中的应用包括:用有限穿越可视图复杂网络对SSVEP脑电 信号 或运动想象EEG脑 电信号 进行数据分析,结合传统的 机器学习 以及 深度学习 理论,能够有效提高意念控制过程中的准确率,便携式脑电采集设备在更广泛领域的应用提供了可能。,下面是便携式脑电采集设备及其在SSVEP和运动想象中的应用专利的具体信息内容。

1.一种便携式脑电采集设备,其特征在于,包括,系统供电电路(1)、脑电极帽转接线(2)、PGA放大电路(3)、AD转换器(4)、STM32处理器(5)和WIFI模(6),其中,所述脑电极帽转接线(2)的输入端连接脑电极帽用于采集脑电信号,输出端依次连接PGA放大电路(3)、AD转换器(4)和STM32处理器(5),所述STM32处理器(5)分别连接PGA放大电路(3)和AD转换器(4)用于对PGA放大电路(3)和AD转换器(4)的工作状态进行控制,所述WIFI模块(6)连接STM32处理器(5)用于通过无线局域网使STM32处理器(5)与上位机进行通信,所述系统供电电路(1)分别连接PGA放大电路(3)、AD转换器(4)、STM32处理器(5)和WIFI模块(6)用于提供电源。
2.根据权利要求1所述的便携式脑电采集设备,其特征在于,其中,
所述的系统供电电路(1)提供3.3伏和5伏两种直流电源;
所述的电极帽转接线(2)采用国际标准1.5mm孔式插口;
所述的PGA放大电路(3)采用可编程增益放大器,能通过程序调节放大倍数,提高测量精度
所述的AD转换器(4)的转换精度为24位;
所述的STM32处理器(5)采用STM32F407ZET6芯片,STM32处理器(5)一方面接收AD转换器发送的数字信号,一方面对AD转换器的工作状态进行控制;
所述的WIFI模块(6)最大传输速率为40Mbps。
3.根据权利要求1所述的便携式脑电采集设备,其特征在于,所述的STM32处理器(5)基于数据压缩感知理论将接收到的脑电信号μ进行压缩,得到观测信号 将观测信号 通过WIFI模块(6)上传至上位机,然后上位机基于正交匹配追踪算法确定稀疏系数的估计,进而重构脑电信号估计 作为实际应用中的脑电信号。
4.根据权利要求3所述的便携式脑电采集设备,其特征在于,所述STM32处理器(5)基于数据压缩感知理论将接收到的脑电信号μ进行压缩,包括:
(1)选用一组稀疏基 B为选稀疏集的个数,对长度为J的脑电信号μ进
行稀疏表示: 系数向量η满足该系数向量η中非
零元素个数 远小于B,其中稀疏基选用正弦基、余弦基、小波基和curvelet基中的一种;
(2)取用大小为M×J的高斯观测矩阵 矩阵中任意元素均服从高斯分布N(0,1/M),观测信号表示为: 其中 是大小为M×B的传感矩阵,M是得到的观测
信号的维数,满足 观测矩阵 和稀疏矩阵 必须满足不相关性的要求。
5.根据权利要求3所述的便携式脑电采集设备,其特征在于,所述上位机基于正交匹配追踪算法确定稀疏系数的估计,进而重构脑电信号估计 包括:
(1)以观测信号 以及传感矩阵G为输入,并设定稀疏度 初始化残差 构建索引集 初始化迭代次数q=1;
(2)分别计算残差rq-1和传感矩阵Gq-1的所有列gj,j∈[1,B]的内积,G0=G,确定内积中最大值对应的标λq∈[1,B],即
(3)更新索引集Λq=Λq-1∪{λq},记录传感矩阵Gq-1中的重建原子集合
即每一次迭代将使得传感矩阵增多一列;
(4)确定重构脑电信号估计 并更新残差 q=q+1;
(5)判断迭代次数是否满足 若满足,则停止迭代;若不满足,返回第(2)步。
6.一种权利要求1所述的便携式脑电采集设备在SSVEP中的应用,其特征在于,包括如下步骤:
1)得到通过便携式脑电采集设备获取并重构的SSVEP脑电信号 L
为SSVEP脑电信号的长度,uh,l是第h个通道第l个数值;
2)通过典型相关分析方法将4个电极的脑电信号融合为一个脑电数据X;
3)通过快速傅里叶变换将融合后的脑电数据X转换为频域下的能量序列F:
F=ABS(FFT(X)),
其中ABS(·)为取绝对值的操作,FFT(·)为对括号中信号进行快速傅里叶变换;
4)对能量序列F建立有限穿越可视图复杂网络,提取所述有限穿越可视图复杂网络中对应视觉刺激界面中4个闪烁频率处的度值,即闪烁频率分别为6.67Hz、7.5Hz、8.57Hz和
10Hz处的节点度值;
5)将每个闪烁频率所对应的度值k1,k2,k3,k4放入到特征向量Kf中,得到Kf=[k1,k2,k3,k4],依据刺激图片的闪烁频率对特征向量Kf进行标记,得到4类特征向量,用特征向量作为输入训练Fisher分类器,得到优化后的Fisher分类器;
6)被试者注视任一个刺激图片,执行步骤1)~步骤5)得到新的特征向量,将新的特征向量输入优化后的Fisher分类器进行分类,分类结果为刺激图片的闪烁频率;
7)建立脑机接口,通过步骤6)的分类结果生成控制指令控制轮椅完成向前直行、向后直行、左转和右转四个动作,这四个动作分别对应实验刺激界面中的四个闪烁图片,实现运动想象轮椅控制,同时通过便携式脑电采集设备实时监测被试者的咬牙信号,若发现咬牙信号,则控制轮椅停止运动,否则轮椅一直进行运动。
7.根据权利要求6所述的便携式脑电采集设备在SSVEP中的应用,其特征在于,步骤1)包括:
(1)设定SSVEP脑电实验刺激界面,包含4个刺激图片,4个刺激图片的闪烁频率依次为
6.67Hz、7.5Hz、8.57Hz和10Hz;
(2)通过便携式脑电采集设备分别获取8个以上被试者中每一个被试者经4个刺激图片诱发的4种SSVEP脑电信号
8.根据权利要求6所述的便携式脑电采集设备在SSVEP中的应用,其特征在于,步骤2)包括:
(1)创建模板信号Y:
其中f1,f2,f3,f4分别表示4个刺激图片的闪烁频率,t为从1到L的序数;
(2)分别确定最优系数向量WU和WY,使得相关系数ρ取得最大值,即:
其中上角标T表示取向量转置的操作,E[·]表示取括号中元素平均值的操作;确定4个电极的SSVEP脑电信号 融合后的脑电数据X=UTWU。
9.根据权利要求6所述的便携式脑电采集设备在SSVEP中的应用,其特征在于,步骤4)中所述的建立有限穿越可视图复杂网络,包括:
(1)针对一元序列D,设置有限穿越视距为Z;
(2)找出一元序列D中所有不相邻的两个节点 和节点 之间的所有节点,对所有节点中的每一个节点 进行判断是否满足条件: 表示一元序列D的第a个节
点, 表示一元序列D的第b个节点, 表示一元序列D的第d个节点,当不满足所述条件的节点数小于等于Z个,则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点 和节点之间存在连边,否则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点 和节点 之间不存在连边;
(3)所有相邻的节点对之间都存在连边;
(4)得到有限穿越可视图复杂网络。
10.根据权利要求6所述的便携式脑电采集设备在SSVEP中的应用,其特征在于,步骤5)包括:
(1)每类特征向量选取90%作为训练样本,其余作为测试样本,设计 个分类器,K取4,得到6个分类器,利用Fisher判别函数在训练样本上训练6个分类器,得到最优投影方向向量,每一个分类器的训练过程包括:
计算各类样本的均值向量mα:
其中,ωα代表类别,Nα是ωα中样本数量,x是ωα中样本;
计算样本类内离散度矩阵Sα和总类内离散度矩阵Sw:
Sw=S1+S2,
计算该分类器的最优投影方向向量w*:
计算训练集内样本在最优投影方向向量w*上的投影点y,
y=(w*)Tx;
(2)每一个Fisher分类器的分类过程包括:
对于输入的样本x,计算在最优投影方向向量w*上的投影点y,
y=(w*)Tx
计算分割阈值y0:
其中,m′α为样本在投影后一维空间的均值,
确定样本类别:
y>y0,x∈ω1
y<y0,x∈ω2;
(3)把测试样本输入第(1)步训练后的6个分类器中进行分类,每个分类器输出一个频率值,并把该频率值下的标记数加1,统计得到标记数最多的频率值作为该样本最终的频率值。
11.一种权利要求1所述的便携式脑电采集设备在运动想象中的应用,其特征在于,包括如下步骤:
1)通过便携式脑电采集设备,开展运动想象脑电实验,获取重构后的被试者运动想象EEG脑电信号,并对运动想象EEG脑电信号进行预处理;
2)对预处理后的运动想象EEG脑电信号 其中,Q为每个信号的数据
长度,Xc,g表示第c个电极采集的信号中第g个数值,采用有限穿越可视图理论构建有限穿越可视图复杂网络Ac,得到大脑多层复杂网络;
3)对于每一个有限穿越可视图复杂网络Ac,分别提取网络指标数据,所述的网络指标数据包括:网络平均节点度、网络平均节点介数、网络平均聚集系数、网络全局聚集系数、网络聚集系数熵和网络平均最短路径;
4)对于任一被试者在每一个具体动作下的大脑多层复杂网络,获得所述网络指标数据,组成一维序列,通过具体动作编号为该一维序列设定标签,构建样本集,搭建初始深度卷积神经网络模型,使用样本集对初始深度卷积神经网络模型进行有监督的训练,得到能够用于运动想象EEG脑电信号有效分类、辨识的深度卷积神经网络模型及参数,其中,所述样本集中的每一个样本由各个网络指标数据组成的一维序列以及相应的标签构成;
5)被试者想象不同的具体的动作,通过便携式脑电采集设备获取运动想象EEG脑电信号并对运动想象EEG脑电信号进行预处理,通过执行步骤2)和步骤3)确定网络指标数据,组成一维序列,并输入训练得到的深度卷积神经网络模型,分类判断出运动想象EEG脑电信号对应的是何种运动想象,以分类结果为基础生成脑控指令,发送给机械外骨骼,控制机械外骨骼运动。
12.根据权利要求11所述的便携式脑电采集设备在运动想象中的应用,其特征在于,步骤1)中所述的运动想象EEG脑电信号是在被试者想象具体动作的同时,通过便携式脑电采集设备,分别获取被试者对应于P3,P4,O1,O2,C3,C4,Cz七个电极的运动想象EEG脑电信号,电极分布符合10~20国际标准导联;包括被试者分别想象四种具体动作:
动作1:通过右手写毛笔字;
动作2:通过左手端起杯喝水;
动作3:通过右脚完成踢足球
动作4:通过左脚完成连续单脚跳动。
13.根据权利要求11所述的便携式脑电采集设备在运动想象中的应用,其特征在于,步骤1)中所述的对运动想象EEG脑电信号进行预处理,是对运动想象EEG脑电信号进行带通滤波并去除伪迹,得到能够用于实现运动想象状态辨识的运动想象EEG脑电信号:
其中,Q为每个信号的数据长度,Xc,g表示第c个电极采集的信号中第g个数值。
14.根据权利要求11所述的便携式脑电采集设备在运动想象中的应用,其特征在于,步骤2)包括:
(1)对每一个电极的运动想象EEG脑电信号 采用有限穿越可视图理论得到有
限穿越可视图复杂网络 β,γ=1,2,...,Q;其中,
所述的建立有限穿越可视图复杂网络,包括:
(1.1)针对一元序列D,设置有限穿越视距为Z;
(1.2)找出一元序列D中所有不相邻的两个节点 和节点 之间的所有节点,对所有节点中的每一个节点 进行判断是否满足条件: 表示一元序列D的第a个
节点,表示一元序列D的第b个节点, 表示一元序列D的第d个节点,当不满足所述条件的节点数小于等于Z个,则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点 和节点 之间存在连边,否则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点 和节点 之间不存在连边;
(1.3)所有相邻的节点对之间都存在连边;
(1.4)得到有限穿越可视图复杂网络。
(2)对7个电极的运动想象EEG脑电信号 分别重复第(1)步,建立7个
有限穿越可视图复杂网络,得到大脑多层复杂网络{A1,A2,...,A7}。
15.根据权利要求11所述的便携式脑电采集设备在运动想象中的应用,其特征在于,步骤4)所述的初始深度卷积神经网络模型依次包括:
一个数据输入层:输入数据即样本为由大脑多层复杂网络的各个网络指标数据构成的一维序列及其对应标签,所述的输入数据是具有标签信息的数据;
三个一维卷积层:卷积核数量为8,卷积核大小为1×7,激活函数是一个线型修正单元,采用Relu激活函数
一个最大池化层:池化核大小为1×4,具体过程描述为 表示对当前池化
覆盖下的输入数据的元素提取最大值作为输出;
四个一维卷积层:卷积核数量为16,卷积核大小为1×7,采用Relu激活函数
一个最大池化层:池化核大小为1×4;
分类器:作为神经网络的输出层,输出分类结果,并用H个神经元表示,通过Soft max函数处理,Soft max函数本质上是归一化的指数函数,定义为 其中s=1…
H,e为自然对数值,zs为第s个神经元的输出,式子中的分母充当了正则项的作用,使得
16.根据权利要求11所述的便携式脑电采集设备在运动想象中的应用,其特征在于,步骤4)所述的训练包括:
(1)将样本集中80%的数据作为训练集,用于训练初始深度卷积神经网络模型,确定初始深度卷积神经网络模型参数,20%的数据作为测试集,用于测试已经训练好的深度卷积神经网络模型的泛化能
(2)将网络权值、阈值取接近于0的随机值,并初始化训练所需的超参数,所述的超参数包括学习率、权值衰减、学习率衰减;
(3)通过反向传播算法对初始深度卷积神经网络模型进行训练;
(4)在训练中通过梯度下降法对网络权值进行调整,公式为 是更新
前的网络权值, 是更新后的网络权值,ξ是学习率, 是误差对网络权值的偏导;
(5)当损失函数停止收敛或者网络训练迭代次数超过设定值M时,停止训练。

说明书全文

便携式脑电采集设备及其在SSVEP和运动想象中的应用

技术领域

[0001] 本发明涉及一种便携式脑电采集设备。特别是涉及一种便携式脑电采集设备及其在 SSVEP和运动想象中的应用。

背景技术

[0002] 脑机接口技术是一种新型人机交互系统,是在人或动物的脑与外部设备间建立的直接连接通路,通过脑电采集设备读取脑中的信号,然后经过计算处理,把这种信号转化为相关的控制指令,并代替人的肢体或语言器官实现人与外界的交流以及对外部环境的控制。脑机接口技术涉及神经科学、信号检测、模式识别等多学科,在残疾人康复训练、生物医学、游戏娱乐等领域发挥着重要的作用。
[0003] 运动想象范式是在没有肢体运动的情况下,被试者自主进行思维活动,想象某一特定的运动场景,产生连续的EEG脑电信号,由脑电采集设备感知,分析系统对其特征进行分析的过程。相对于其他脑电实验范式而言,运动想象不需要附加刺激界面,在设备要求便携化、小型化时,仍可以保证较高的准确率,具有独特的优点。
[0004] 稳态视觉诱发电位(SSVEP),是指依托于一定频率的刺激在EEG脑电信号中产生响应的倍频响应,其中刺激主要是固定频率的闪烁,它可以可靠地应用于脑机接口系统。相对于其他信号而言,SSVEP具有更高的分类准确率、更快的信息传输率,系统和实验设计更加简便,而且需要的训练次数也比较少。
[0005] 复杂性科学作为当代科学发展的前沿领域之一,其发展广泛应用于社会学、物理统计学、经济学、控制学、工程学、生物医学等多个跨学科研究领域,该理论提出的非线性时间序列分析方法,把每个通道当作一个节点,把节点间的某种联系当作连边,这样该系统就被抽象成了复杂网络,把对的研究同一转化为对复杂网络的分析,利用网络指标可分析非线性时间序列的内在特性,准确把握其规律。

发明内容

[0006] 本发明所要解决的技术问题是,提供一种能够实现多自由度意念控制的便携式脑电采集设备及其在SSVEP和运动想象中的应用。
[0007] 本发明所采用的技术方案是:一种便携式脑电采集设备,包括,系统供电电路、脑电极帽转接线、PGA放大电路、AD转换器、STM32处理器和WIFI模,其中,所述脑电极帽转接线的输入端连接脑电极帽用于采集脑电信号,输出端依次连接PGA放大电路、AD转换器和STM32处理器,所述STM32处理器分别连接PGA放大电路和AD转换器用于对PGA 放大电路和AD转换器的工作状态进行控制,所述WIFI模块连接STM32处理器用于通过无线局域网使STM32处理器与上位机进行通信,所述系统供电电路分别连接PGA放大电路、 AD转换器、STM32处理器和WIFI模块用于提供电源。其中,
[0008] 所述的系统供电电路提供3.3伏和5伏两种直流电源;
[0009] 所述的电极帽转接线采用国际标准1.5mm孔式插口;
[0010] 所述的PGA放大电路采用可编程增益放大器,能通过程序调节放大倍数,提高测量精度
[0011] 所述的AD转换器的转换精度为24位;
[0012] 所述的STM32处理器采用STM32F407ZET6芯片,STM32处理器一方面接收AD转换器发送的数字信号,一方面对AD转换器的工作状态进行控制;
[0013] 所述的WIFI模块最大传输速率为40Mbps。
[0014] 所述的STM32处理器(5)基于数据压缩感知理论将接收到的脑电信号μ进行压缩,得到观测信号 将观测信号 通过WIFI模块(6)上传至上位机,然后上位机基于正交匹配追踪算法确定稀疏系数的估计,进而重构脑电信号估计 作为实际应用中的脑电信号。
[0015] 所述STM32处理器(5)基于数据压缩感知理论将接收到的脑电信号μ进行压缩,包括:
[0016] (1)选用一组稀疏基 B为选稀疏集的个数,对长度为J的脑电信号μ进行稀疏表示: υ=1,2,…,B,系数向量η满足该系数向量η中
非零元素个数 远小于B,其中稀疏基选用正弦基、余弦基、小波基和curvelet基中的一种;
[0017] (2)取用大小为M×J的高斯观测矩阵 矩阵中任意元素均服从高斯分布N(0,1/M),观测信号表示为: 其中 是大小为M×B的传感矩阵,M是得到的观测信号的维数,满足 观测矩阵 和稀疏矩阵 必须满足不相关性的要求。
[0018] 所述上位机基于正交匹配追踪算法确定稀疏系数的估计,进而重构脑电信号估计包括:
[0019] (1)以观测信号 以及传感矩阵G为输入,并设定稀疏度 初始化残差 构建索引集 初始化迭代次数q=1;
[0020] (2)分别计算残差rq-1和传感矩阵Gq-1的所有列gj,j∈[1,B]的内积,G0=G,确定内积中最大值对应的标λq∈[1,B],即
[0021] (3)更新索引集Λq=Λq-1∪{λq},记录传感矩阵Gq-1中的重建原子集合即每一次迭代将使得传感矩阵增多一列;
[0022] (4)确定重构脑电信号估计 并更新残差 q=q+1;
[0023] (5)判断迭代次数是否满足 若满足,则停止迭代;若不满足,返回第(2)步。
[0024] 一种便携式脑电采集设备在SSVEP中的应用,包括如下步骤:
[0025] 1)得到通过便携式脑电采集设备获取并重构的SSVEP脑电信号 h=1,2,3,4, L为SSVEP脑电信号的长度,uh,l是第h个通道第l个数值;
[0026] 2)通过典型相关分析方法将4个电极的脑电信号融合为一个脑电数据X;
[0027] 3)通过快速傅里叶变换将融合后的脑电数据X转换为频域下的能量序列F:
[0028] F=ABS(FFT(X)),
[0029] 其中ABS(·)为取绝对值的操作,FFT(·)为对括号中信号进行快速傅里叶变换;
[0030] 4)对能量序列F建立有限穿越可视图复杂网络,提取所述有限穿越可视图复杂网络中对应视觉刺激界面中4个闪烁频率处的度值,即闪烁频率分别为6.67Hz、7.5Hz、8.57Hz和10 Hz处的节点度值;
[0031] 5)将每个闪烁频率所对应的度值k1,k2,k3,k4放入到特征向量Kf中,得到 Kf=[k1,k2,k3,k4],依据刺激图片的闪烁频率对特征向量Kf进行标记,得到4类特征向量,用特征向量作为输入训练Fisher分类器,得到优化后的Fisher分类器;
[0032] 6)被试者注视任一个刺激图片,执行步骤1)~步骤5)得到新的特征向量,将新的特征向量输入优化后的Fisher分类器进行分类,分类结果为刺激图片的闪烁频率;
[0033] 7)建立脑机接口,通过步骤6)的分类结果生成控制指令控制轮椅完成向前直行、向后直行、左转和右转四个动作,这四个动作分别对应实验刺激界面中的四个闪烁图片,实现运动想象轮椅控制,同时通过便携式脑电采集设备实时监测被试者的咬牙信号,若发现咬牙信号,则控制轮椅停止运动,否则轮椅一直进行运动。
[0034] 步骤1)包括:
[0035] (1)设定SSVEP脑电实验刺激界面,包含4个刺激图片,4个刺激图片的闪烁频率依次为6.67Hz、7.5Hz、8.57Hz和10Hz;
[0036] (2)通过便携式脑电采集设备分别获取8个以上被试者中每一个被试者经4个刺激图片诱发的4种SSVEP脑电信号 h=1,2,3,4。
[0037] 步骤2)包括:
[0038] (1)创建模板信号Y:
[0039]
[0040] 其中f1,f2,f3,f4分别表示4个刺激图片的闪烁频率,t为从1到L的序数;
[0041] (2)分别确定最优系数向量WU和WY,使得相关系数ρ取得最大值,即:
[0042]
[0043] 其中上角标T表示取向量转置的操作,E[·]表示取括号中元素平均值的操作;确定4个电极的SSVEP脑电信号 h=1,2,3,4融合后的脑电数据X=UTWU。
[0044] 步骤4)中所述的建立有限穿越可视图复杂网络,包括:
[0045] (1)针对一元序列D,设置有限穿越视距为Z;
[0046] (2)找出一元序列D中所有不相邻的两个节点 和节点 之间的所有节点,对所有节点中的每一个节点 进行判断是否满足条件: 表示一元序列D的第a 个节点,表示一元序列D的第b个节点, 表示一元序列D的第d个节点,当不满足所述条件的节点数小于等于Z个,则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点 和节点 之间存在连边,否则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点和节点 之间不存在连边;
[0047] (3)所有相邻的节点对之间都存在连边;
[0048] (4)得到有限穿越可视图复杂网络。
[0049] 步骤5)包括:
[0050] (1)每类特征向量选取90%作为训练样本,其余作为测试样本,设计 个分类器, K取4,得到6个分类器,利用Fisher判别函数在训练样本上训练6个分类器,得到最优投影方向向量,每一个分类器的训练过程包括:
[0051] 计算各类样本的均值向量mα:
[0052]
[0053] 其中,ωα代表类别,Nα是ωα中样本数量,x是ωα中样本;
[0054] 计算样本类内离散度矩阵Sα和总类内离散度矩阵Sw:
[0055]
[0056] Sw=S1+S2,
[0057] 计算该分类器的最优投影方向向量w*:
[0058]
[0059] 计算训练集内样本在最优投影方向向量w*上的投影点y,
[0060] y=(w*)Tx;
[0061] (2)每一个Fisher分类器的分类过程包括:
[0062] 对于输入的样本x,计算在最优投影方向向量w*上的投影点y,
[0063] y=(w*)Tx
[0064] 计算分割阈值y0:
[0065]
[0066] 其中,m′α为样本在投影后一维空间的均值, α=1,2;
[0067] 确定样本类别:
[0068] y>y0,x∈ω1
[0069] y<y0,x∈ω2;
[0070] (3)把测试样本输入第(1)步训练后的6个分类器中进行分类,每个分类器输出一个频率值,并把该频率值下的标记数加1,统计得到标记数最多的频率值作为该样本最终的频率值。
[0071] .一种便携式脑电采集设备在运动想象中的应用,包括如下步骤:
[0072] 1)通过便携式脑电采集设备,开展运动想象脑电实验,获取重构后的被试者运动想象 EEG脑电信号,并对运动想象EEG脑电信号进行预处理;
[0073] 2)对预处理后的运动想象EEG脑电信号 c=1,2,...,7,其中,Q为每个信号的数据长度,Xc,g表示第c个电极采集的信号中第g个数值,采用有限穿越可视图理论构建有限穿越可视图复杂网络Ac,得到大脑多层复杂网络;
[0074] 3)对于每一个有限穿越可视图复杂网络Ac,分别提取网络指标数据,所述的网络指标数据包括:网络平均节点度、网络平均节点介数、网络平均聚集系数、网络全局聚集系数、网络聚集系数熵和网络平均最短路径;
[0075] 4)对于任一被试者在每一个具体动作下的大脑多层复杂网络,获得所述网络指标数据,组成一维序列,通过具体动作编号为该一维序列设定标签,构建样本集,搭建初始深度卷积神经网络模型,使用样本集对初始深度卷积神经网络模型进行有监督的训练,得到能够用于运动想象EEG脑电信号有效分类、辨识的深度卷积神经网络模型及参数,其中,所述样本集中的每一个样本由各个网络指标数据组成的一维序列以及相应的标签构成;
[0076] 5)被试者想象不同的具体的动作,通过便携式脑电采集设备获取运动想象EEG脑电信号并对运动想象EEG脑电信号进行预处理,通过执行步骤2)和步骤3)确定网络指标数据,组成一维序列,并输入训练得到的深度卷积神经网络模型,分类判断出运动想象EEG脑电信号对应的是何种运动想象,以分类结果为基础生成脑控指令,发送给机械外骨骼,控制机械外骨骼运动;
[0077] 步骤1)中所述的运动想象EEG脑电信号是在被试者想象具体动作的同时,通过便携式脑电采集设备,分别获取被试者对应于P3,P4,O1,O2,C3,C4,Cz七个电极的运动想象 EEG脑电信号,电极分布符合10~20国际标准导联;包括被试者分别想象四种具体动作:
[0078] 动作1:通过右手写毛笔字;
[0079] 动作2:通过左手端起杯喝水;
[0080] 动作3:通过右脚完成踢足球
[0081] 动作4:通过左脚完成连续单脚跳动。
[0082] 步骤1)中所述的对运动想象EEG脑电信号进行预处理,是对运动想象EEG脑电信号进行带通滤波并去除伪迹,得到能够用于实现运动想象状态辨识的运动想象EEG脑电信号:
[0083]
[0084] 其中,Q为每个信号的数据长度,Xc,g表示第c个电极采集的信号中第g个数值。
[0085] 步骤2)包括:
[0086] (1)对每一个电极的运动想象EEG脑电信号 采用有限穿越可视图理论得到有限穿越可视图复杂网络 β,γ=1,2,...,Q;其中,
[0087] 所述的建立有限穿越可视图复杂网络,包括:
[0088] (1.1)针对一元序列D,设置有限穿越视距为Z;
[0089] (1.2)找出一元序列D中所有不相邻的两个节点 和节点 之间的所有节点,对所有节点中的每一个节点 进行判断是否满足条件: 表示一元序列D的第a 个节点,表示一元序列D的第b个节点, 表示一元序列D的第d个节点,当不满足所述条件的节点数小于等于Z个,则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点和节点 之间存在连边,否则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点和节点 之间不存在连边;
[0090] (1.3)所有相邻的节点对之间都存在连边;
[0091] (1.4)得到有限穿越可视图复杂网络。
[0092] (2)对7个电极的运动想象EEG脑电信号 c=1,2,...,7分别重复第(1)步,建立7个有限穿越可视图复杂网络,得到大脑多层复杂网络{A1,A2,...,A7}。
[0093] 步骤4)所述的初始深度卷积神经网络模型依次包括:
[0094] 一个数据输入层:输入数据即样本为由大脑多层复杂网络的各个网络指标数据构成的一维序列及其对应标签,所述的输入数据是具有标签信息的数据;
[0095] 三个一维卷积层:卷积核数量为8,卷积核大小为1×7,激活函数是一个线型修正单元,采用Relu激活函数
[0096] 一个最大池化层:池化核大小为1×4,具体过程描述为 表示对当前池化核覆盖下的输入数据的元素提取最大值作为输出;
[0097] 四个一维卷积层:卷积核数量为16,卷积核大小为1×7,采用Relu激活函数[0098] 一个最大池化层:池化核大小为1×4;
[0099] 分类器:作为神经网络的输出层,输出分类结果,并用H个神经元表示,通过Softmax函数处理,Softmax函数本质上是归一化的指数函数,定义为 其中s=1···H, e为自然对数值,zs为第s个神经元的输出,式子中的分母充当了正则项的作用,使得
[0100] 步骤4)所述的训练包括:
[0101] (1)将样本集中80%的数据作为训练集,用于训练初始深度卷积神经网络模型,确定初始深度卷积神经网络模型参数,20%的数据作为测试集,用于测试已经训练好的深度卷积神经网络模型的泛化能
[0102] (2)将网络权值、阈值取接近于0的随机值,并初始化训练所需的超参数,所述的超参数包括学习率、权值衰减、学习率衰减;
[0103] (3)通过反向传播算法对初始深度卷积神经网络模型进行训练;
[0104] (4)在训练中通过梯度下降法对网络权值进行调整,公式为 是更新前的网络权值, 是更新后的网络权值,ξ是学习率, 是误差对网络权值的偏导;
[0105] (5)当损失函数停止收敛或者网络训练迭代次数超过设定值M时,停止训练。
[0106] 本发明的便携式脑电采集设备及其在SSVEP和运动想象中的应用,具有如下有益效果:
[0107] (1)与被试者连接的便携式脑电采集设备的供电电路中的电压对人体均无伤害,保证了被试者的安全,与上位机采用WIFI通信,打破线缆的距离限制,增大被试者与计算机之间的空间,采用高精度AD转换器增加了转换精度,大大提高了数据的准确性,为后续分析提供基础数据支撑
[0108] (2)将压缩感知理论嵌入便携式脑电采集设备当中,在保证数据精度的前提下,很好地提升了数据的传输效率,为便携式脑电采集设备在更广泛领域的应用提供了可能;
[0109] (3)将有限穿越可视图复杂网络理论应用了SSVEP脑电信号以及运动想象EEG脑电信号分析当中,充分发挥有限穿越可视图复杂网络的抗噪能力,能够有效提升脑机接口中控制的准确率;
[0110] (4)结合有限穿越可视图复杂网络与深度卷积神经网络,有效挖掘运动想象EEG脑电信号中的潜在特征,提高分类准确率。附图说明
[0111] 图1是本发明便携式脑电采集设备的结构框图
[0112] 图2是本发明的便携式脑电采集设备在SSVEP中应用的流程图
[0113] 图3是本发明的便携式脑电采集设备在运动想象中应用的流程图;
[0114] 图4是有限穿越可视图复杂网络的原理图。

具体实施方式

[0115] 下面结合实施例和附图对本发明的便携式脑电采集设备及其在SSVEP和运动想象中的应用做出详细说明。
[0116] 如图1所示,本发明的一种便携式脑电采集设备,包括,系统供电电路1、脑电极帽转接线2、PGA放大电路3、AD转换器4、STM32处理器5和WIFI模块6,其中,所述脑电极帽转接线2的输入端连接脑电极帽用于采集脑电信号,输出端依次连接PGA放大电路3、 AD转换器4和STM32处理器5,所述STM32处理器5分别连接PGA放大电路3和AD转换器4用于对PGA放大电路3和AD转换器4的工作状态进行控制,所述WIFI模块6连接STM32处理器5用于通过无线局域网使STM32处理器5与上位机进行通信,所述系统供电电路1分别连接PGA放大电路3、AD转换器4、STM32处理器5和WIFI模块6用于提供电源。其中,
[0117] 所述的系统供电电路1提供3.3伏和5伏两种直流电源;
[0118] 所述的电极帽转接线2采用国际标准1.5mm孔式插口,能与各种进口和国产脑电极帽配套使用;
[0119] 所述的PGA放大电路3采用可编程增益放大器,是一种通用性很强的放大器,能通过程序调节放大倍数,提高测量精度;
[0120] 所述的AD转换器4的转换精度为24位,将模拟信号转换为数字信号;
[0121] 所述的STM32处理器5采用STM32F407ZET6芯片,具有高性能、低成本、低功耗的特点;STM32处理器5一方面接收AD转换器发送的数字信号,一方面对AD转换器的工作状态进行控制;
[0122] 所述的WIFI模块6最大传输速率为40Mbps,该WIFI模块与STM32处理器相连接,接收STM32处理器传输的所有数据,通过无线局域网将所有数据上传至上位机。
[0123] 所述的STM32处理器5基于数据压缩感知理论将接收到的脑电信号μ进行压缩,得到观测信号 将观测信号 通过WIFI模块6上传至上位机,然后上位机基于正交匹配追踪算法确定稀疏系数的估计,进而重构脑电信号估计 作为实际应用中的脑电信号。
[0124] 所述STM32处理器5基于数据压缩感知理论将接收到的脑电信号μ进行压缩,包括:
[0125] (1)选用一组稀疏基 B为选稀疏集的个数,对长度为J的脑电信号μ进行稀疏表示: υ=1,2,…,B,系数向量η满足该系数向量η中
非零元素个数 远小于B,其中稀疏基选用正弦基、余弦基、小波基和curvelet基中的一种;
[0126] (2)取用大小为M×J的高斯观测矩阵 矩阵中任意元素均服从高斯分布N(0,1/M),观测信号表示为: 其中 是大小为M×B的传感矩阵,M是得到的观测信号的维数,满足 观测矩阵 和稀疏矩阵 必须满足不相关性的要求。
[0127] 所述上位机基于正交匹配追踪算法确定稀疏系数的估计,进而重构脑电信号估计包括:
[0128] (1)以观测信号 以及传感矩阵G为输入,并设定稀疏度 初始化残差 构建索引集 初始化迭代次数q=1;
[0129] (2)分别计算残差rq-1和传感矩阵Gq-1的所有列gj,j∈[1,B]的内积,G0=G,确定内积中最大值对应的角标λq∈[1,B],即
[0130] (3)更新索引集Λq=Λq-1∪{λq},记录传感矩阵Gq-1中的重建原子集合即每一次迭代将使得传感矩阵增多一列;
[0131] (4)确定重构脑电信号估计 并更新残差 q=q+1;
[0132] (5)判断迭代次数是否满足 若满足,则停止迭代;若不满足,返回第(2)步。
[0133] 本发明的便携式脑电采集设备在SSVEP中的应用,如图2所示,包括如下步骤:
[0134] 1)得到通过便携式脑电采集设备获取并重构的SSVEP脑电信号 h=1,2,3,4, L为SSVEP脑电信号的长度,uh,l是第h个通道第l个数值;包括:
[0135] (1)设定SSVEP脑电实验刺激界面,包含4个刺激图片,4个刺激图片的闪烁频率依次为6.67Hz、7.5Hz、8.57Hz和10Hz;
[0136] (2)通过便携式脑电采集设备分别获取8个以上被试者中每一个被试者经4个刺激图片诱发的4种SSVEP脑电信号 h=1,2,3,4。
[0137] 2)通过典型相关分析方法将4个电极的脑电信号融合为一个脑电数据X;包括:
[0138] (1)创建模板信号Y:
[0139]
[0140] 其中f1,f2,f3,f4分别表示4个刺激图片的闪烁频率,t为从1到L的序数;
[0141] (2)分别确定最优系数向量WU和WY,使得相关系数ρ取得最大值,即:
[0142]
[0143] 其中上角标T表示取向量转置的操作,E[·]表示取括号中元素平均值的操作;确定4个电极的SSVEP脑电信号 h=1,2,3,4融合后的脑电数据X=UTWU。
[0144] 3)通过快速傅里叶变换将融合后的脑电数据X转换为频域下的能量序列F:
[0145] F=ABS(FFT(X)),
[0146] 其中ABS(·)为取绝对值的操作,FFT(·)为对括号中信号进行快速傅里叶变换;
[0147] 4)对能量序列F建立有限穿越可视图复杂网络,提取所述有限穿越可视图复杂网络中对应视觉刺激界面中4个闪烁频率处的度值,即闪烁频率分别为6.67Hz、7.5Hz、8.57Hz和10 Hz处的节点度值;所述的建立有限穿越可视图复杂网络,如图4所示,包括:
[0148] (1)针对一元序列D,设置有限穿越视距为Z;
[0149] (2)找出一元序列D中所有不相邻的两个节点 和节点 之间的所有节点,对所有节点中的每一个节点 进行判断是否满足条件: 表示一元序列D的第a 个节点,表示一元序列D的第b个节点, 表示一元序列D的第d个节点,当不满足所述条件的节点数小于等于Z个,则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点 和节点 之间存在连边,否则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点和节点 之间不存在连边;
[0150] (3)所有相邻的节点对之间都存在连边;
[0151] (4)得到有限穿越可视图复杂网络。包括:
[0152] (1)每类特征向量选取90%作为训练样本,其余作为测试样本,设计 个分类器, K取4,得到6个分类器,利用Fisher判别函数在训练样本上训练6个分类器,得到最优投影方向向量,每一个分类器的训练过程包括:
[0153] 计算各类样本的均值向量mα:
[0154]
[0155] 其中,ωα代表类别,Nα是ωα中样本数量,x是ωα中样本;
[0156] 计算样本类内离散度矩阵Sα和总类内离散度矩阵Sw:
[0157]
[0158] Sw=S1+S2,
[0159] 计算该分类器的最优投影方向向量w*:
[0160]
[0161] 计算训练集内样本在最优投影方向向量w*上的投影点y,
[0162] y=(w*)Tx;
[0163] (2)每一个Fisher分类器的分类过程包括:
[0164] 对于输入的样本x,计算在最优投影方向向量w*上的投影点y,
[0165] y=(w*)Tx
[0166] 计算分割阈值y0:
[0167]
[0168] 其中,m′α为样本在投影后一维空间的均值, α=1,2;
[0169] 确定样本类别:
[0170] y>y0,x∈ω1
[0171] y<y0,x∈ω2;
[0172] (3)把测试样本输入第(1)步训练后的6个分类器中进行分类,每个分类器输出一个频率值,并把该频率值下的标记数加1,统计得到标记数最多的频率值作为该样本最终的频率值。
[0173] 5)将每个闪烁频率所对应的度值k1,k2,k3,k4放入到特征向量Kf中,得到 Kf=[k1,k2,k3,k4],依据刺激图片的闪烁频率对特征向量Kf进行标记,得到4类特征向量,用特征向量作为输入训练Fisher分类器,得到优化后的Fisher分类器;
[0174] 6)被试者注视任一个刺激图片,执行步骤1)~步骤5)得到新的特征向量,将新的特征向量输入优化后的Fisher分类器进行分类,分类结果为刺激图片的闪烁频率;
[0175] 7)建立脑机接口,通过步骤6)的分类结果生成控制指令控制轮椅完成向前直行、向后直行、左转和右转四个动作,这四个动作分别对应实验刺激界面中的四个闪烁图片,实现运动想象轮椅控制,同时通过便携式脑电采集设备实时监测被试者的咬牙信号,若发现咬牙信号,则控制轮椅停止运动,否则轮椅一直进行运动。
[0176] 本发明的便携式脑电采集设备在运动想象中的应用,如图3所示,包括如下步骤:
[0177] 1)通过便携式脑电采集设备,开展运动想象脑电实验,获取重构后的被试者运动想象 EEG脑电信号,并对运动想象EEG脑电信号进行预处理;其中,
[0178] 所述的运动想象EEG脑电信号是在被试者想象具体动作的同时,通过便携式脑电采集设备,分别获取被试者对应于P3,P4,O1,O2,C3,C4,Cz七个电极的运动想象EEG脑电信号,电极分布符合10~20国际标准导联;包括被试者分别想象四种具体动作:
[0179] 动作1:通过右手写毛笔字;
[0180] 动作2:通过左手端起水杯喝水;
[0181] 动作3:通过右脚完成踢足球射门;
[0182] 动作4:通过左脚完成连续单脚跳动。
[0183] 所述的对运动想象EEG脑电信号进行预处理,是对运动想象EEG脑电信号进行带通滤波并去除伪迹,得到能够用于实现运动想象状态辨识的运动想象EEG脑电信号:
[0184]
[0185] 其中,Q为每个信号的数据长度,Xc,g表示第c个电极采集的信号中第g个数值。
[0186] 2)对预处理后的运动想象EEG脑电信号 c=1,2,...,7,其中,Q为每个信号的数据长度,Xc,g表示第c个电极采集的信号中第g个数值,采用有限穿越可视图理论构建有限穿越可视图复杂网络Ac,得到大脑多层复杂网络;包括:
[0187] (1)对每一个电极的运动想象EEG脑电信号 采用有限穿越可视图理论得到有限穿越可视图复杂网络 β,γ=1,2,...,Q;其中,
[0188] 所述的建立有限穿越可视图复杂网络,如图4所示,包括:
[0189] (1.1)针对一元序列D,设置有限穿越视距为Z;
[0190] (1.2)找出一元序列D中所有不相邻的两个节点 和节点 之间的所有节点,对所有节点中的每一个节点 进行判断是否满足条件: 表示一元序列D的第a 个节点,表示一元序列D的第b个节点, 表示一元序列D的第d个节点,当不满足所述条件的节点数小于等于Z个,则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点和节点 之间存在连边,否则认定在有限穿越视距为Z的有限穿越可视图复杂网络中节点和节点 之间不存在连边;
[0191] (1.3)所有相邻的节点对之间都存在连边;
[0192] (1.4)得到有限穿越可视图复杂网络。
[0193] (2)对7个电极的运动想象EEG脑电信号 c=1,2,...,7分别重复第(1)步,建立7个有限穿越可视图复杂网络,得到大脑多层复杂网络{A1,A2,...,A7}。
[0194] 3)对于每一个有限穿越可视图复杂网络Ac,分别提取网络指标数据,所述的网络指标数据包括:网络平均节点度、网络平均节点介数、网络平均聚集系数、网络全局聚集系数、网络聚集系数熵和网络平均最短路径;其中,
[0195] (3.1)网络平均节点度O:任意节点θ的度oθ表示与该节点相连的连边个数,则其中<·>表示取平均值的操作,l′为网络中节点个数;
[0196] (3.2)网络平均节点介数 表示节点介数,其中σλκ为连接节点λ和节点κ的最短路径数目,σλκ(θ)为经过节点θ的连接节点λ和节点κ的最短路径数目;
[0197] (3.3)网络平均聚集系数 表示节点聚集系数,其中τθ,Δ表示复杂网络中包含节点θ的封闭三角形的数目,τθ表示复杂网络中从节点θ出发至少有两个边的三角形的数目;
[0198] (3.4)网络全局聚集系数
[0199] (3.5)网络聚集系数熵
[0200] (3.6)网络平均最短路径 其中节点λ和节点κ互不相同,且Uλκ表示节点λ和节点κ之间的最短路径长度。
[0201] 4)对于任一被试者在每一个具体动作下的大脑多层复杂网络,获得所述网络指标数据,组成一维序列,通过具体动作编号为该一维序列设定标签,构建样本集,搭建初始深度卷积神经网络模型,使用样本集对初始深度卷积神经网络模型进行有监督的训练,得到能够用于运动想象EEG脑电信号有效分类、辨识的深度卷积神经网络模型及参数,其中,所述样本集中的每一个样本由各个网络指标数据组成的一维序列以及相应的标签构成;
[0202] 所述的初始深度卷积神经网络模型依次包括:
[0203] 一个数据输入层:输入数据即样本为由大脑多层复杂网络的各个网络指标数据构成的一维序列及其对应标签,所述的输入数据是具有标签信息的数据;
[0204] 三个一维卷积层:卷积核数量为8,卷积核大小为1×7,激活函数是一个线型修正单元,采用Relu激活函数
[0205] 一个最大池化层:池化核大小为1×4,具体过程描述为 表示对当前池化核覆盖下的输入数据的元素提取最大值作为输出;
[0206] 四个一维卷积层:卷积核数量为16,卷积核大小为1×7,采用Relu激活函数[0207] 一个最大池化层:池化核大小为1×4;
[0208] 分类器:作为神经网络的输出层,输出分类结果,并用H个神经元表示,通过Softmax函数处理,Soft max函数本质上是归一化的指数函数,定义为 其中s=1···H, e为自然对数值,zs为第s个神经元的输出,式子中的分母充当了正则项的作用,使得
[0209] 所述的训练包括:
[0210] (1)将样本集中80%的数据作为训练集,用于训练初始深度卷积神经网络模型,确定初始深度卷积神经网络模型参数,20%的数据作为测试集,用于测试已经训练好的深度卷积神经网络模型的泛化能力;
[0211] (2)将网络权值、阈值取接近于0的随机值,并初始化训练所需的超参数,所述的超参数包括学习率、权值衰减、学习率衰减;
[0212] (3)通过反向传播算法对初始深度卷积神经网络模型进行训练;
[0213] (4)在训练中通过梯度下降法对网络权值进行调整,公式为 是更新前的网络权值, 是更新后的网络权值,ξ是学习率, 是误差对网络权值的偏导;
[0214] (5)当损失函数停止收敛或者网络训练迭代次数超过设定值M时,停止训练。
[0215] 5)被试者想象不同的具体的动作,通过便携式脑电采集设备获取运动想象EEG脑电信号并对运动想象EEG脑电信号进行预处理,通过执行步骤2)和步骤3)确定网络指标数据,组成一维序列,并输入训练得到的深度卷积神经网络模型,分类判断出运动想象EEG脑电信号对应的是何种运动想象,以分类结果为基础生成脑控指令,发送给机械外骨骼,控制机械外骨骼运动。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈