基于多参数公共空间频率模式算法的脑电特征提取

摘  要： 关键词： 脑机接口；脑电图；CSSP算法；特征提取；支持向量机
       多年来，人们一直希望直接利用大脑中的电生理信号建立一种向外界传递信息和发送命令的通道，即所谓的脑机接口系统[1]。脑机接口BCI(Brain Computer Interface)是一种不依赖于大脑的外周神经与肌肉正常输出通路的通讯和控制系统[1]。由于非植入式BCI所获取的脑电信号EEG(Electro Encephalo Gram)较微弱，同时受到心电、肌电等信号的干扰，使得提取脑电信号的特征存在较大的困难，所以寻求一种有效且适用性强的信号特征提取方法决定着BCI系统能否迅速发展及广泛应用。
     目前，脑电信号的特征提取通常有时域、频域和空域方式，提取方法主要有以下几种:(1)AR(Autoregressive)模型谱估计。该方法获取了频域上的谱信息但损失了时域上的信息，在对时间敏感的信号应用上，得不到好的效果[2]。(2)时域分析法。时域分析法的一个优点是它能够获取时域和频域上的特征，但是算法比较复杂，计算量比较大，不符合BCI系统的实时性要求。(3)公共空间模式CSP(Common Spatial Pattern)。CSP算法是基于两个协方差矩阵的同时对角化来获取空间滤波器。该空间滤波器对两种类型的信号进行空间滤波时起到相反的效果，当其中一类经过滤波后的方差越大时，另一类反而越小，这使得两类的特征存在比较大的差异，对它们进行分类就较容易[3]。但是CSP算法的不足之处就是无法操作频域上的信息。(4)公共空间频率模式CSSP(Common Spatio-Spectral Pattern)。CSSP算法原理上与CSP算法类似，但是CSSP算法优于CSP算法之处是它不仅可以获取空域上的信息，同时还可以对频域上的信息进行操作，使得提取的特征更加明显。

1.2 多参数的CSSP算法
     本文研究的数据集的每一个序列是从3个电极(C3、Cz、C4)中采集得到的，所以输入信号是一个三维的矩阵，而不是一维的向量。在使用CSSP算法对脑电信号进行特征提取时，所引入的延时因子τ，只是一个单一的变量值[5-6]，对所有从不同电极所采集的信号进行相同时间的延时。考虑到从不同电极所采集到的信号反映不同类别的特征所在的频段会存在一定的偏差，所以用单一的延时因子τ对所有序列进行延时，并不能将两类别差别最大的特征提取出来。因此，本文提出将单一的延时因子τ替换为多维变量[τ]，即对不同电极采集的信号采用不同的延时因子。按此方法，多参数的CSSP算法的公式演变为：

2 实验数据特征提取
     本文研究的数据集来自于2003年BCI竞赛中运动想象的标准数据集，由奥地利工业大学生物工程学的医学信息部门提供。
     该数据集采自于想象左右手运动的方法来控制一个运动杆的任务。该数据集包含280次实验，从中随机抽取140次实验作为训练集，另外的140次作为测试集。每次实验所经历的时间为9 s。前2 s为准备阶段，受试者精神处于放松状态，无任何动作。在第2 s出现一个声音的刺激信号，提示受试者做好准备，从第2 s到第3 s，屏幕上出现一个十字形的持续提示信号，从第3 s开始，屏幕上随机出现向左或向右的箭头，提示实验者进入想象左右手运动的实验阶段。
2.1 数据预处理
     单边的肢体运动或想象运动时，大脑同侧产生事件相关同步电位ERS(Event-Related Synchronization)，大脑对侧产生事件相关去同步电位ERD(Event-Ralated Desynchronization)[7]，这两种电位主要位于运动感觉区并且主要反应在mu节律和beta节律这两个波段。由于这两个波的频率带主要集中在8 Hz～30 Hz，所以需对运动想象脑电信号进行8 Hz～30 Hz带通滤波的预处理。本文通过加窗滤波的方式对数据集的每一通道进行带通滤波。采用的窗函数为矩形窗，阶数为60。
2.2 特征提取
     本文所采用的数据集的每一次实验所用的时间为9 s，而真正的实验阶段为3 s～9 s，所以只将3 s～9 s的数据用于特征提取。在参考文献[2]中，实验者在整个实验过程中，在第4 s～5 s对左右手想象运动表现得最为明显，可以认为这一时间段是实验者脑活动最为活跃的阶段，因此，用该时间段的数据进行特征提取。
     设经过预处理得到的第4 s～5 s的两类样本数据为Tr_13×128×140和 Tr_23×128×140(三个下标分别代表电极数、每秒的采样数、试验次数)，按照三维的延时参数[τ]对样本数据进行延时操作，得到δ[τ]Tr_13×128×140和δ[τ]Tr_23×128×140，将这两类数据各自按行连接起来，构成X_13×17920和 X_23×17920，利用式(4)构造出两类的输入信号，利用多参数的CSSP算法求出空间滤波器W。本文只采用能量最大w1与最小w2的两个方向组成空间频率滤波器，再运用式(5)求出输出信号Z。将经两个方向投影后的输出信号Z分别做方差运算作为脑电分类的特征，即特征f为：
     f=(var(z1)  var(z2))                              (9)
3 实验步骤及数据结果分析
     本文采用支持向量机的分类方法。支持向量机是通过某个内核函数将输入信号映射到一个高维特征空间，进而在高维特征空间寻找一个最优的分类面。本文采用高斯核作为核函数。
3.1 实验步骤
     (1)给定延时参数[τ]3×1以及初始值、支持向量机参数c和g的范围并选择c和g的步长。本文给定的[τ]3×1范围为[1 1 1]′～[6 6 6]′(采样频率为128 Hz时，每一个单位代表1/128 s)，初始值为[1 1 1]′(参数为[0 0 0]′时即不采用延时参数，本质上即CSP算法)。c和g参数的范围为20～212，步长为2。
     (2)利用多参数的CSSP算法以及给定的[τ]3×1值提取出140个训练样本特征。利用网格搜索法与五折交叉验证法，在给定的支持向量机参数范围内对训练集进行训练，求出使得分类正确率最高的支持向量机参数，得到最佳的分类模型。
     (3)在所给的延时参数[τ]3×1范围内，计算出每一个[τ]3×1所对应的测试集的特征，利用已获得的最佳分类模型对这些测试集特征进行分类。
     (4)在给定的范围内，赋予[τ]3×1新的值，返回第二步骤，直到[τ]3×1参数范围内所有值都用来对训练集进行训练为止。
3.2 数据结果分析
     利用上面所阐述的方法，通过寻找训练集与测试集各自所对应的最佳延时参数，可以提取出训练集与测试集不同类别间差别较大的特征，对于不同的测试者产生的脑电信号可以训练出一个最佳的分类模型。令C3、Cz、C4这三个电极采集的数据序列的延时参数[τ]3×1为τ1、τ2、τ3三个变量。在训练集采用最佳的延时参数[3 4 5]的情况下，将参数τ1固定，改变参数τ2、τ3得到测试集的正确率，如图1所示。