专利名称：基于非线性动力学的运动想象脑电信号特征提取方法
技术领域：
本发明属于信息技术领域，更一步涉及在生命科学领域中应用脑-机接口(Brain-Computer Interface, BCI)系统的基于非线性动力学的运动想象脑电信号特征提取方法。本发明用于对运动想象脑电信号特征的提取，通过分类器对特征进行分类，实现单侧手指运动想象的判别，最终可应用于轮椅、鼠标等运动想象BCI在线系统，从而提高有缺陷人士与外界的交流能力。
背景技术：
在准备和执行单侧手指运动想象时，人们大脑皮层的功能性连接发生改变，从而导致其对侧大脑运动感觉区mu和beta节律的脑电信号能量减弱，而其同侧大脑运动感觉区mu和beta节律的脑电信号能量增强。这种想象单侧手指运动时特定脑区特定频率脑电信号的能量变化，被称为事件相关去同步现象。该现象是判别左右手指运动想象脑电信号最根本的特征。因此，通过受试者运动想象时脑电信号的分析，判别受试者运动想象的方向，从而实现对外界装置的控制。北京工业大学拥有的专利技术“想象单侧肢体运动的脑电特征的提取方法”(申请号 200810056839. 1，申请目 2008. 01. 25，授权号 CN101219048B,授权日 2010. 06. 23)中提出一种想象单侧肢体运动的脑电特征的提取方法。该专利技术将共同空间模式(CommonSpatial Pattern, CSP)方法和线性判别分析(FDA)相结合，降低了输入矢量的维数，提高了分类器的推广性，在一定程度上提高了分类正确率。共同空间模式方法对于信噪比高的运动想象脑电信号的分类虽然能获得很好的正确率，但是该专利技术仍然存在的不足是，在实际应用中，通常不具备实验室较为理想的数据采集环境，所以实际应用采集的运动想象脑电信号往往与实验室采集的数据相比包含大量噪声，信噪比低。用共同空间模式方法分析处理实际应用采集的运动想象脑电信号的正确率不高，稳定性差，算法运算量大，不利于实际应用的推广。
发明内容
本发明的目的在于克服上述已有脑电信号特征提取方法的不足，提出一种基于非线性动力学的运动想象脑电信号特征的提取方法。该方法充分考虑运动想象时脑电信号在实际环境中信噪比低的特征，保证了一定分类正确率。实现本发明方法的主要思路是将采集的多通道脑电信号做空间滤波、基线校正、带通滤波处理，用非线性动力学的多尺度熵方法提取脑电信号特征向量；利用朴素贝叶斯分类器分类，采取交叉验证方法，得到分类正确率。根据上述主要思路，本发明方法的具体步骤如下(I)采集数据脑电信号采集系统通过受试者佩戴的电极帽，采集受试者想象单侧手指运动的脑电信号，受试者手指不能有实际上的动作。
⑵空间滤波采用共同平均参考的方法，将受试者电极帽上每个电极采集的脑电信号减去所有电极采集的脑电信号的均值，得到共同平均参考后的脑电信号。(3)基线校正将共同平均参考后的脑电信号减去基线，得到基线校正后的脑电信号。(4)带通滤波利用有限脉冲响应滤波器，对基线校正后的脑电信号进行带通滤波，得到不同频带的子信号。
(5)特征提取对所有不同频带的子信号在不同的尺度下进行粗粒化处理，对经过粗粒化之后的所有子信号进行样本熵估计，将所有得到的样本熵估计值组成特征向量。(6)分类利用朴素贝叶斯分类器，对得到的特征向量进行分类，采取交叉验证方法对分类结果进行验证，得到分类正确率。本发明与现有技术相比具有如下优点第一，本发明由于考虑了运动想象时脑电信号在实际环境中信噪比低的情况，利用多尺度熵作为特征向量，克服了现有技术中要求运动想象脑电信号信噪比高的、稳定性差的不足，使得本发明具有能够有效保证稳定的脑电信号分类正确率的优点。第二，本发明由于考虑了运动想象时脑电信号处理的实时性，利用多尺度熵作为特征向量，算法运算量小，克服了现有技术中算法运算量大的问题，使得本发明具有能够有效的提供更好实时性的优点。
图I为本发明的流程图；图2为本发明采集数据步骤中屏幕提示符示意图；图3为本发明采集数据步骤的示意图。
具体实施例方式下面结合附图I对本发明做进一步的描述。步骤I，采集数据脑电信号采集系统通过受试者佩戴的电极帽，采集想象单侧手指运动的脑电信号。脑电信号由佩戴在受试者头上的电极帽获取，并通过脑电放大器放大以及模/数转换器转换，输入计算机，以信号电压幅值的形式存储并显示。受试者佩戴电极帽，坐在椅子上平视距其Im左右的显示器。脑电信号采集系统的采样频率为250Hz，测试电极分别为C3和C4，脑电信号的波动值域为±100 μ V。在本发明采集数据步骤中需要给受试者运动想象提示，具体的提示种类参照图2。在图2中显示器的提示符有三种，显示器中标有十字叉的代表准备提示，显示器中标有向左的箭头代表想象左手运动提不，向右的箭头代表想象右手运动提不。在本发明采集数据实施例中图2所示的运动想象提示按设定的时间次序出现，参考图3。在图3中开始O秒至3秒时屏幕出现准备提示，并在第2秒出现短暂的提示音(lkHz,70ms)。随后，屏幕出现想象左手运动提示或想象右手提示，并持续I. 25秒。在第4秒，受试者开始想象相应的手指运动，并持续3秒。每次实施例后休息I. 5秒至2. 5秒。左右手指运动想象实验各120次，顺序上随机排列。步骤2，空间滤波采用共同平均参考的方法，对得到的脑电信号计算每个电极信号所有样本点的均值，将所有电极脑电信号的样本点的均值做平均得到总体样本点的均值，将所有脑电信号的每个样本点减去总体样本点的均值，从而得到共同平均参考之后的脑电信号。步骤3，基线校正以受试者想象单侧手指运动之前的200ms脑电信号为基线，将共同平均参考空间滤波后的脑电信号减去基线，得到基线校正后的脑电信号。步骤4，带通滤波利用有限脉冲响应滤波器，对基线校正后的脑电信号进行带通滤波，得到特定频带的子信号，重复执行带通滤波，得到频带为8-13HZ的子信号和18-26HZ子信号。步骤5，特征提取对采集的多个电极中的C3和C4电极，分别通过步骤4的带通滤波，得到两个电极、两个频带的4个子信号。将这4个子信号分别在尺度I、尺度2、尺度3下进行粗粒化处理，得到12个粗粒化序列。本发明实施例中，尺度I的数值取1，尺度2的数值取2，尺度3的数值取3。对得到的12个粗粒化序列分别计算在阈值r = O. 15XSD下的样本熵，其中SD为求样本熵的信号的标准差，在本发明的实施例中先通过matlab中的var函数取得信号的方差，将方差开平方得到该标准差。将得到的样本熵作为特征值，从而得到12个特征值，将这12个特征值组成一个特征向量。对两个不同电极、两个不同频带的子信号在不同尺度下进行粗粒化处理的具体步骤是第一步按照下式计算粗粒化样本点；
权利要求
1.基于非线性动力学的运动想象脑电信号特征提取方法，具体步骤如下 (1)采集数据 脑电信号采集系统通过受试者佩戴的电极帽，采集受试者想象单侧手指运动时的电极帽中多个电极的脑电信号； (2)空间滤波 采用共同平均参考的方法，对采集到的多个电极的脑电信号进行处理，得到共同平均参考后的脑电信号； (3)基线校正 将共同平均参考后的脑电信号减去基线，得到基线校正后的脑电信号； (4)带通滤波 利用有限脉冲响应滤波器，对基线校正后的脑电信号进行带通滤波，得到特定频带的子信号，重复执行带通滤波，得到两个不同频带的子信号； (5)特征提取 对两个不同频带的子信号在不同尺度下进行粗粒化处理，对粗粒化处理得到的粗粒化序列进行样本熵估计，将两个不同频带、不同尺度下得到的样本熵估计值组成特征向量； (6)分类 利用朴素贝叶斯分类器，对得到的特征向量进行分类，采取交叉验证方法对分类结果进行验证，得到分类正确率。
2.根据权利要求I所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于步骤(2)中所述的共同平均参考方法如下 第一步对步骤(I)得到的每个电极的脑电信号的所有样本点取均值，得到每个电极的均值； 第二步将第一步得到的每个电极的均值求和再除以电极数目，得到总体样本点的均值； 第三步将所有电极的脑电信号的每个样本点分别减去总体样本点的均值。
3.根据权利要求I所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于步骤(3)所述的基线为受试者想象单侧手指运动之前的200ms脑电信号。
4.根据权利要求I所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于步骤(4)所述的不同频带的子信号包括8-13HZ的子信号以及18-26HZ子信号。
5.根据权利要求I所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于步骤(5)中所述不同的尺度是指尺度的数量不少于3个，且相互不重复。
6.根据权利要求I所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于步骤(5)中所述的粗粒化的具体步骤是 第一步按照下式计算粗粒化样本点；
7.根据权利要求I所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于步骤(5)中信号的样本熵估计方法如下； 第一步按照下式将粗粒化序列的样本点组成一个矢量组 uDi(P) = iyP yP+i. · · · yP+m-J 其中um(p)为一个m维的矢量组;111为矢量组的维数；p为粗粒化序列的第P个样本点，其范围是I < P < N-m ；N为粗粒化序列的样本点总数;yp为粗粒化序列的样本点； 第二步按照下式计算两个m维矢量组之间的距离 d = max{Iy(k+a)-y (I+a) |} 其中d为两个m维矢量组之间的距离邱为矢量组的维数；max{}为求最大值的函数；y(k+a)为矢量组^仏)中第a个元素；y(l+a)为矢量组Um (I)中第a个元素；um(k)为第k个矢量组；um(l)为第I个矢量组，k和I互不相等； 第三步设定阈值r，r = O. 15XSD，其中SD为子信号的标准方差，统计距离d小于阈值r的矢量组的数目并记为模板匹配数，计算模板匹配数与距离总个数的比值； 第四步按照下式求模板匹配数与距离总个数的比值对矢量组总数的平均值
8.根据权利要求I所述的运动想象脑电信号特征的提取方法，其特征在于步骤(7)中所述的交叉验证方法是指，选择步骤(6)得到的任意一个脑电信号的特征向量作为测试数据，其余脑电信号的特征向量作为训练数据；分类器利用训练数据建立分类模型，将测试数据代入分类模型得到分类类别，比较分类类别与实际任务类别，得到正确分类或错误分类的结果；依次将每个脑电信号的特征向量作为一次测试数据，统计所有测试数据的分类结果，得到分类正确率。
全文摘要
本发明公开了一种基于非线性动力学的运动想象脑电信号特征提取方法。该方法解决了运动想象脑电信号在实际应用环境中可能出现信噪比低从而导致分类正确率低和不稳定的问题。本发明具体步骤包括(1)采集数据(2)空间滤波；(3)基线校正；(4)带通滤波；(5)特征提取；(6)分类。本发明具有能够有效保证稳定的脑电信号分类正确率、有效的提供更好实时性的优点。
文档编号G06K9/62GK102940490SQ20121043425
公开日2013年2月27日 申请日期2012年10月19日 优先权日2012年10月19日
发明者刘鹏, 何嘉全, 朱孟波, 赵瑞霞, 胡凯, 朱振营, 秦伟 申请人:西安电子科技大学