CN102940490A - 基于非线性动力学的运动想象脑电信号特征提取方法 - Google Patents

Abstract

本发明公开了一种基于非线性动力学的运动想象脑电信号特征提取方法。该方法解决了运动想象脑电信号在实际应用环境中可能出现信噪比低从而导致分类正确率低和不稳定的问题。本发明具体步骤包括：(1)采集数据：(2)空间滤波；(3)基线校正；(4)带通滤波；(5)特征提取；(6)分类。本发明具有能够有效保证稳定的脑电信号分类正确率、有效的提供更好实时性的优点。

Description

基于非线性动力学的运动想象脑电信号特征提取方法

技术领域

本发明属于信息技术领域，更一步涉及在生命科学领域中应用脑-机接口(Brain-Computer Interface，BCI)系统的基于非线性动力学的运动想象脑电信号特征提取方法。本发明用于对运动想象脑电信号特征的提取，通过分类器对特征进行分类，实现单侧手指运动想象的判别，最终可应用于轮椅、鼠标等运动想象BCI在线系统，从而提高有缺陷人士与外界的交流能力。

背景技术

在准备和执行单侧手指运动想象时，人们大脑皮层的功能性连接发生改变，从而导致其对侧大脑运动感觉区mu和beta节律的脑电信号能量减弱，而其同侧大脑运动感觉区mu和beta节律的脑电信号能量增强。这种想象单侧手指运动时特定脑区特定频率脑电信号的能量变化，被称为事件相关去同步现象。该现象是判别左右手指运动想象脑电信号最根本的特征。因此，通过受试者运动想象时脑电信号的分析，判别受试者运动想象的方向，从而实现对外界装置的控制。

北京工业大学拥有的专利技术“想象单侧肢体运动的脑电特征的提取方法”(申请号200810056839.1，申请目2008.01.25，授权号CN101219048B，授权日2010.06.23)中提出一种想象单侧肢体运动的脑电特征的提取方法。该专利技术将共同空间模式(Common Spatial Pattern，CSP)方法和线性判别分析(FDA)相结合，降低了输入矢量的维数，提高了分类器的推广性，在一定程度上提高了分类正确率。共同空间模式方法对于信噪比高的运动想象脑电信号的分类虽然能获得很好的正确率，但是该专利技术仍然存在的不足是，在实际应用中，通常不具备实验室较为理想的数据采集环境，所以实际应用采集的运动想象脑电信号往往与实验室采集的数据相比包含大量噪声，信噪比低。用共同空间模式方法分析处理实际应用采集的运动想象脑电信号的正确率不高，稳定性差，算法运算量大，不利于实际应用的推广。

发明内容

本发明的目的在于克服上述已有脑电信号特征提取方法的不足，提出一种基于非线性动力学的运动想象脑电信号特征的提取方法。该方法充分考虑运动想象时脑电信号在实际环境中信噪比低的特征，保证了一定分类正确率。

实现本发明方法的主要思路是：将采集的多通道脑电信号做空间滤波、基线校正、带通滤波处理，用非线性动力学的多尺度熵方法提取脑电信号特征向量；利用朴素贝叶斯分类器分类，采取交叉验证方法，得到分类正确率。

根据上述主要思路，本发明方法的具体步骤如下：

(1)采集数据：

脑电信号采集系统通过受试者佩戴的电极帽，采集受试者想象单侧手指运动的脑电信号，受试者手指不能有实际上的动作。

(2)空间滤波：

采用共同平均参考的方法，将受试者电极帽上每个电极采集的脑电信号减去所有电极采集的脑电信号的均值，得到共同平均参考后的脑电信号。

(3)基线校正：

将共同平均参考后的脑电信号减去基线，得到基线校正后的脑电信号。

(4)带通滤波：

利用有限脉冲响应滤波器，对基线校正后的脑电信号进行带通滤波，得到不同频带的子信号。

(5)特征提取：

对所有不同频带的子信号在不同的尺度下进行粗粒化处理，对经过粗粒化之后的所有子信号进行样本熵估计，将所有得到的样本熵估计值组成特征向量。

(6)分类：

利用朴素贝叶斯分类器，对得到的特征向量进行分类，采取交叉验证方法对分类结果进行验证，得到分类正确率。

本发明与现有技术相比具有如下优点：

第一，本发明由于考虑了运动想象时脑电信号在实际环境中信噪比低的情况，利用多尺度熵作为特征向量，克服了现有技术中要求运动想象脑电信号信噪比高的、稳定性差的不足，使得本发明具有能够有效保证稳定的脑电信号分类正确率的优点。

第二，本发明由于考虑了运动想象时脑电信号处理的实时性，利用多尺度熵作为特征向量，算法运算量小，克服了现有技术中算法运算量大的问题，使得本发明具有能够有效的提供更好实时性的优点。

附图说明

图1为本发明的流程图；

图2为本发明采集数据步骤中屏幕提示符示意图；

图3为本发明采集数据步骤的示意图。

具体实施方式

下面结合附图1对本发明做进一步的描述。

步骤1，采集数据：

脑电信号采集系统通过受试者佩戴的电极帽，采集想象单侧手指运动的脑电信号。脑电信号由佩戴在受试者头上的电极帽获取，并通过脑电放大器放大以及模/数转换器转换，输入计算机，以信号电压幅值的形式存储并显示。

受试者佩戴电极帽，坐在椅子上平视距其1m左右的显示器。脑电信号采集系统的采样频率为250Hz，测试电极分别为C3和C4，脑电信号的波动值域为±100μV。

在本发明采集数据步骤中需要给受试者运动想象提示，具体的提示种类参照图2。在图2中显示器的提示符有三种，显示器中标有十字叉的代表准备提示，显示器中标有向左的箭头代表想象左手运动提示，向右的箭头代表想象右手运动提示。

在本发明采集数据实施例中图2所示的运动想象提示按设定的时间次序出现，参考图3。在图3中开始0秒至3秒时屏幕出现准备提示，并在第2秒出现短暂的提示音(1kHz，70ms)。随后，屏幕出现想象左手运动提示或想象右手提示，并持续1.25秒。在第4秒，受试者开始想象相应的手指运动，并持续3秒。每次实施例后休息1.5秒至2.5秒。左右手指运动想象实验各120次，顺序上随机排列。步骤2，空间滤波：

采用共同平均参考的方法，对得到的脑电信号计算每个电极信号所有样本点的均值，将所有电极脑电信号的样本点的均值做平均得到总体样本点的均值，将所有脑电信号的每个样本点减去总体样本点的均值，从而得到共同平均参考之后的脑电信号。

步骤3，基线校正：

以受试者想象单侧手指运动之前的200ms脑电信号为基线，将共同平均参考空间滤波后的脑电信号减去基线，得到基线校正后的脑电信号。

步骤4，带通滤波：

利用有限脉冲响应滤波器，对基线校正后的脑电信号进行带通滤波，得到特定频带的子信号，重复执行带通滤波，得到频带为8-13Hz的子信号和18-26Hz子信号。

步骤5，特征提取：

对采集的多个电极中的C3和C4电极，分别通过步骤4的带通滤波，得到两个电极、两个频带的4个子信号。将这4个子信号分别在尺度1、尺度2、尺度3下进行粗粒化处理，得到12个粗粒化序列。本发明实施例中，尺度1的数值取1，尺度2的数值取2，尺度3的数值取3。对得到的12个粗粒化序列分别计算在阈值r＝0.15×SD下的样本熵，其中SD为求样本熵的信号的标准差，在本发明的实施例中先通过matlab中的var函数取得信号的方差，将方差开平方得到该标准差。将得到的样本熵作为特征值，从而得到12个特征值，将这12个特征值组成一个特征向量。

对两个不同电极、两个不同频带的子信号在不同尺度下进行粗粒化处理的具体步骤是：

第一步：按照下式计算粗粒化样本点；

$y_j^{\left(\mathrm{\τ}\right)}=\frac{1}{\mathrm{\τ}}\underset{i}{\overset{\mathrm{j\τ}}{\mathrm{\Σ}}}{x}_i$

其中，

为粗粒化样本点；τ为尺度的值；j为粗粒化样本点下标，代表第j个粗粒化样本点，其范围是1≤j≤N\τ；N为子信号的样本点总数；i为子信号的第i个样本点，其范围是1≤i≤N，其初值为(j-1)τ+1；x_i为子信号的样本点；

第二步：按照粗粒化样本点下标从小到大的顺序，将得到的粗粒化样本点构建成粗粒化序列。

对粗粒化处理得到的粗粒化序列进行样本熵估计方法如下：

第一步：按照下式将粗粒化序列的样本点组成一个m维矢量组：

u_m(p)＝{y_p,y_p+1,...y_p+m-1}，1≤p≤N-m

其中：u_m(p)为一个m维的矢量组；m为矢量组的维数；p为粗粒化序列的第p个样本点，其范围是1≤p≤N-m；N为粗粒化序列的样本点总数；y_p为粗粒化序列的样本点；

第二步：按照下式计算两个m维矢量组之间的距离：

d＝max{|y(k+a)-y(1+a)|：0≤a≤m-1，k≠1}

其中：d为两个m维矢量组之间的距离；m为矢量组的维数；max{}为求最大

值的函数；y(k+a)为矢量组u_m(k)中第a个元素；y(1+a)为矢量组u_m(l)中第a个

元素；u_m(k)为第k个矢量组；u_m(l)为第1个矢量组；

第三步：根据阈值r，对每个p≤N-m的值，统计距离d小于阈值r的数目(称为

模板匹配数)及此数目与距离总数N-m-1的比值，该比值记作

$B_r^m\left(c\right)=N^m\left(p\right)/(N-m-1)$

其中：

为模板匹配数与其距离总个数的比值；c为第c个模板匹配数与距离总个数的比值，其范围是1≤c≤N-m；N为粗粒化序列的样本点总数；m为矢量组的维数；N^m(p)为模板匹配数；p为第p个矢量组；

第四步：按照下式求模板匹配数与距离总个数的比值对矢量组总数的平均值

$B^m\left(r\right)={(N-m)}^{-1}\underset{c=1}{\overset{N-m}{\mathrm{\Σ}}}{B}_r^m\left(c\right)$

其中：B^m(r)为模板匹配总数与距离总个数的比值对矢量组总数的平均值；r为阈值；m为矢量组的维数；N为粗粒化序列的样本点总数；

为模板匹配数与其距离总个数的比值；c为第c个模板匹配数与距离总个数的比值，其范围是1≤c≤N-m；

第五步：增大维数m重复第一步、第二步、第三步、第四步，得到m+1维的模板匹配数与距离总个数的比值对矢量组总数的平均值B^m+1(r)；

$B^{m+1}\left(r\right)={(N-m)}^{-1}\underset{i=1}{\overset{N-m}{\mathrm{\Σ}}}{B}_r^{m+1}\left(c\right)$

其中：B^m+1(r)为模板匹配总数与距离总个数的比值对矢量组总数的平均值；r为阈值；m为矢量组的维数；N为粗粒化序列的样本点总数；

为模板匹配数与其距离总个数的比值；c为第c个模板匹配数与距离总个数的比值，其范围是1≤c≤N-m；

第六步：熵估计值为

SampEn(m，r，N)＝-In[B^m+1(r)/B^m(r)]

其中：SampEn(m，r，N)为熵估计值的表达式；m为矢量组的维数；r为阈值；N为粗粒化序列的样本点总数；In()为对数函数；B^m+1(r)为m+1维的模板匹配数与距离总数的比值；Bm(r)为m维的模板匹配数与距离总数的比值。

步骤6，特征分类：

利用朴素贝叶斯分类器，对特征向量分类，采取交叉验证方法，对分类结果进行交叉验证，从而得到分类正确率。

选择步骤5得到的任意一个脑电信号的特征向量作为测试数据，其余脑电信号的特征向量作为训练数据；朴素贝叶斯分类器利用训练数据建立分类模型，将测试数据代入分类模型得到分类类别，比较分类类别与实际任务类别，得到正确分类或错误分类的结果；依次将每个脑电信号的特征向量作为一次测试数据，统计所有测试数据的分类结果，得到分类正确率。

Claims (8)

1.基于非线性动力学的运动想象脑电信号特征提取方法，具体步骤如下：

(1)采集数据：

脑电信号采集系统通过受试者佩戴的电极帽，采集受试者想象单侧手指运动时的电极帽中多个电极的脑电信号；

(2)空间滤波：

采用共同平均参考的方法，对采集到的多个电极的脑电信号进行处理，得到共同平均参考后的脑电信号；

(3)基线校正：

将共同平均参考后的脑电信号减去基线，得到基线校正后的脑电信号；

(4)带通滤波：

利用有限脉冲响应滤波器，对基线校正后的脑电信号进行带通滤波，得到特定频带的子信号，重复执行带通滤波，得到两个不同频带的子信号；

(5)特征提取：

对两个不同频带的子信号在不同尺度下进行粗粒化处理，对粗粒化处理得到的粗粒化序列进行样本熵估计，将两个不同频带、不同尺度下得到的样本熵估计值组成特征向量；

(6)分类：

利用朴素贝叶斯分类器，对得到的特征向量进行分类，采取交叉验证方法对分类结果进行验证，得到分类正确率。

2.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于：步骤(2)中所述的共同平均参考方法如下：

第一步：对步骤(1)得到的每个电极的脑电信号的所有样本点取均值，得到每个电极的均值；

第二步：将第一步得到的每个电极的均值求和再除以电极数目，得到总体样本点的均值；

第三步：将所有电极的脑电信号的每个样本点分别减去总体样本点的均值。

3.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于：步骤(3)所述的基线为受试者想象单侧手指运动之前的200ms脑电信号。

4.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于：步骤(4)所述的不同频带的子信号包括8-13Hz的子信号以及18-26Hz子信号。

5.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于：步骤(5)中所述不同的尺度是指尺度的数量不少于3个，且相互不重复。

6.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于：步骤(5)中所述的粗粒化的具体步骤是：

第一步：按照下式计算粗粒化样本点；

$y_j^{\left(\mathrm{\τ}\right)}=\frac{1}{\mathrm{\τ}}\underset{i}{\overset{\mathrm{j\τ}}{\mathrm{\Σ}}}{x}_i$

其中，

为粗粒化样本点；τ为尺度的值；j为粗粒化样本点的序号，代表第j个粗粒化样本点，其范围是1≤j≤N\τ；N为子信号的样本点总数；i为子信号的第i个样本点，其范围是1≤i≤N，其初值为(j-1)τ+1；x_i为子信号的样本点；

第二步：按照粗粒化样本点的序号从小到大的顺序，将得到的粗粒化样本点构建成粗粒化序列。

7.根据权利要求1所述的基于非线性动力学的运动想象脑电信号特征提取方法，其特征在于：步骤(5)中信号的样本熵估计方法如下；

第一步：按照下式将粗粒化序列的样本点组成一个矢量组：

u_m(p)＝{y_p,y_p+1,...y_p+m-1}

其中：u_m(p)为一个m维的矢量组；m为矢量组的维数；p为粗粒化序列的第p个样本点，其范围是1≤p≤N-m；N为粗粒化序列的样本点总数；y_p为粗粒化序列的样本点；

第二步：按照下式计算两个m维矢量组之间的距离：

d＝max{|y(k+a)-y(1+a)|}

其中：d为两个m维矢量组之间的距离；m为矢量组的维数；max{}为求最大值的函数；y(k+a)为矢量组u_m(k)中第a个元素；y(1+a)为矢量组u_m(l)中第a个元素；u_m(k)为第k个矢量组；u_m(l)为第1个矢量组，k和1互不相等；

第三步：设定阈值r，r＝0.15×SD，其中SD为子信号的标准方差，统计距离d小于阈值r的矢量组的数目并记为模板匹配数，计算模板匹配数与距离总个数的比值；

第四步：按照下式求模板匹配数与距离总个数的比值对矢量组总数的平均值

$B^m\left(r\right)={(N-m)}^{-1}\underset{c=1}{\overset{N-m}{\mathrm{\Σ}}}{B}_r^m\left(c\right)$

其中：B^m(r)为模板匹配总数与距离总个数的比值对矢量组总数的平均值；r为阈值；m为矢量组的维数；N为粗粒化序列的样本点总数；

为模板匹配数与其距离总个数的比值；c为第c个模板匹配数与距离总个数的比值，其范围是1≤c≤N-m；

第五步：增大维数m重复第一步、第二步、第三步、第四步，得到m+1维的模板匹配数与距离总个数的比值对矢量组总数的平均值B^m+1(r)；

第六步：根据下式计算熵估计值：

SampEn(m，r，N)＝-In[B^m+1(r)/B^m(r)]

其中：SampEn(m，r，N)为熵估计值的表达式；m为矢量组的维数；r为阈值；N为粗粒化序列的样本点总数；In()为对数函数；B^m+1(r)为m+1维的模板匹配数与距离总数的比值；B^m(r)为m维的模板匹配数与距离总数的比值。

8.根据权利要求1所述的运动想象脑电信号特征的提取方法，其特征在于：步骤(7)中所述的交叉验证方法是指，选择步骤(6)得到的任意一个脑电信号的特征向量作为测试数据，其余脑电信号的特征向量作为训练数据；分类器利用训练数据建立分类模型，将测试数据代入分类模型得到分类类别，比较分类类别与实际任务类别，得到正确分类或错误分类的结果；依次将每个脑电信号的特征向量作为一次测试数据，统计所有测试数据的分类结果，得到分类正确率。

Images