本发明专利技术公开了一种基于非线性动力学的运动想象脑电信号特征提取方法。该方法解决了运动想象脑电信号在实际应用环境中可能出现信噪比低从而导致分类正确率低和不稳定的问题。本发明专利技术具体步骤包括：(1)采集数据：(2)空间滤波；(3)基线校正；(4)带通滤波；(5)特征提取；(6)分类。本发明专利技术具有能够有效保证稳定的脑电信号分类正确率、有效的提供更好实时性的优点。
【技术实现步骤摘要】
本专利技术属于信息
，更一步涉及在生命科学领域中应用脑-机接口(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)特征提取 对两个不同频带的子信号在不同尺度下进行粗粒化处理，对本文档来自技高网...
【技术保护点】
基于非线性动力学的运动想象脑电信号特征提取方法，具体步骤如下：(1)采集数据：脑电信号采集系统通过受试者佩戴的电极帽，采集受试者想象单侧手指运动时的电极帽中多个电极的脑电信号；(2)空间滤波：采用共同平均参考的方法，对采集到的多个电极的脑电信号进行处理，得到共同平均参考后的脑电信号；(3)基线校正：将共同平均参考后的脑电信号减去基线，得到基线校正后的脑电信号；(4)带通滤波：利用有限脉冲响应滤波器，对基线校正后的脑电信号进行带通滤波，得到特定频带的子信号，重复执行带通滤波，得到两个不同频带的子信号；(5)特征提取：对两个不同频带的子信号在不同尺度下进行粗粒化处理，对粗粒化处理得到的粗粒化序列进行样本熵估计，将两个不同频带、不同尺度下得到的样本熵估计值组成特征向量；(6)分类：利用朴素贝叶斯分类器，对得到的特征向量进行分类，采取交叉验证方法对分类结果进行验证，得到分类正确率。
【技术特征摘要】
【专利技术属性】
技术研发人员：刘鹏，何嘉全，朱孟波，赵瑞霞，胡凯，朱振营，秦伟，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。