本发明属于脑电信号处理的技术领域，尤其涉及一种基于层次分析的脑电信号特征分类的二分类方法。
背景技术：
脑机接口(braincomputerinterface,bci)是指能够使人不依赖周围神经系统和肌肉，而与外界进行通信或控制的设备。脑机接口技术作为一种全新的通信和控制技术，可以为思维正常但有严重运动障碍的患者提供语言交流和环境控制手段。此外，脑机接口技术不仅应用于未患者提供语言交流和环境控制，在自动控制、军事认知等科学领域也有潜在的应用价值。鉴于其巨大的应用前景，脑机接口已经引起国际科学界的高度重视，称为脑科学、康复工程、生物医学工程及人机控制领域的一个研究热点。
在所有能够被监测到的反映大脑活动的信号中，由于脑电信号(electroencephalogram，eeg)具有较好的时间分辨率，监测仪器较简单，非入侵等优点，被大部分脑机接口系统所采纳。
目前，对脑电信号特征分类的算法有神经网络、支持向量机、线性判别式分析和朴素贝叶斯等等分类器。这些分类器可以较为准确的根据脑电信号分类出结果。但也存在一些问题，例如，每次只会得到简单的分类结果而没有分类结果的可信度。
不同时间段或者不同频率段提取的特征对于分类准确度的贡献不同。而目前分类结果一般仅仅参考单一时间段或者频率，或者以不同时间段或者频率段的分类结果做专家投票机制，得出最后的分类结果。
综上所述，目前脑电信号在特征分类方面存在如下问题：分类依据特征单一或无法较好地综合多特征，分类精度不高；分类结果可信度未知，即只能简单地得出分类结果而不能得出对应结果的可信性。针对脑电信号特征分类的问题，本发明提出了一种基于层次分析的脑电信号特征分类的二分类方法，该方法可以提高分类精度和得到分类结果的可信度。
技术实现要素：
本发明为了解决上述问题，克服现有技术中脑电信号在特征分类方面存在的分类依据特征单一或无法较好地综合多特征、分类精度不高、分类结果可信度未知的问题，提供一种基于层次分析的脑电信号特征分类的二分类方法，该方法可以提高分类精度和得到分类结果的可信度。
为了实现上述目的，本发明采用如下技术方案：
一种基于层次分析的脑电信号特征分类的二分类方法，该方法的具体步骤包括：
(1)脑机接口系统采集运动想象脑电信号；
(2)将脑机接口系统采集的运动想象脑电信号进行格式转换，并进行数据预处理；
(3)采用特征提取算法分别提取不同时间段和频率段的脑电信号特征；
(4)将特征提取后的脑电信号征和其对应的标记按照合适比例分为训练集和测试集；
(5)根据训练集的脑电信号数据采用二分类算法得到预测分类模型；
(6)采用预测分类模型计算测试集的测试集分类结果，并且分别得出不同时间段和频率段结果的正确率；
(7)对不同时间段或频率段分类结果的正确率与实际结果进行两两比较，构造判断矩阵a；
(8)对判断矩阵a进行一致性检验，并判断是否通过一致性检验，若通过，进入步骤(9)，否则，返回步骤(7)；
(9)确定不同时间段或者不同频率段判断矩阵a的权重w；
(10)用权重w对不同时间段或者频率段训练结果加权得到最后结果及其对应的可信度。
进一步的，所述步骤(1)中脑机接口系统采集运动想象脑电信号的具体步骤为：
脑机接口系统对用户进行不同心理作业时的运动想象脑电信号进行采集，分别包括想象左手运动时的运动想象脑电信号和想象右手运动时的运动想象脑电信号，并将采集的运动想象脑电信号进行存储。
进一步的，所述步骤(2)中将脑机接口系统采集的运动想象脑电信号进行格式转换并进行数据预处理的具体步骤为：
(2-1)将采集的运动想象脑电信号由二维信号转换为三维信号：将步骤(1)中存储的脑机接口系统采集的运动想象脑电信号构建二维空间矩阵；
(2-2)对步骤(2-1)格式转换后的运动想象脑电信号进行空间滤波、频带的带通滤波、去除眼电和基线矫正预处理。
进一步的，所述步骤(3)中采用特征提取算法分别提取不同时间段和频率段的脑电信号特征的具体步骤为：
采用快速傅里叶变换、共空域模式或小波包变换方法提取不同时间段和频率段脑电信号特征。
进一步的，所述步骤(6)中分类结果分别为：
想象左手运动时的运动想象脑电信号对应的结果为1，想象右手运动时的运动想象脑电信号对应的结果为-1。
进一步的，所述步骤(7)中构造判断矩阵a的具体步骤为：
根据对不同时间段或频率段分类结果的正确率与实际结果进行两两比较得到的正确率大小的不同，引用数字1-9及其倒数作为标度构造出判断矩阵a。
进一步的，所述步骤(7)中，
数字1作为标度代表的含义为两个因素相比，具有相同的重要性，
数字3作为标度代表的含义为两个因素相比，前者比后者稍重要，
数字5作为标度代表的含义为两个因素相比，前者比后者明显重要，
数字7作为标度代表的含义为两个因素相比，前者比后者强烈重要，
数字9作为标度代表的含义为两个因素相比，前者比后者极端重要，
数字2、数字4、数字6和数字8分别表示上述相邻数字标度判断的中间值。
进一步的，所述步骤(8)中，对判断矩阵a进行一致性检验的具体步骤为：
(8-1)计算判断矩阵a的一致性指标ci：
其中，λmax为判断矩阵最大特征值，n为阶数；
(8-2)按照平均随机一致性指标表查找一致性指标ri；
(8-3)根据一致性指标ci与一致性指标ri的商计算一致性比例cr。
进一步的，所述步骤(9)中，采用几何平均法、算数平均法、特征向量法、或最小二乘法确定不同时间段或者不同频率段判断矩阵a的权重w。
进一步的，所述步骤(10)中，可信度分为非常好、较好、较差和非常差四个等级。
本发明的有益效果：
1.提高脑电信号在特征分类精度：本发明的一种基于层次分析的脑电信号特征分类的二分类方法，采用特征提取算法分别提取不同时间段和频率段的脑电信号特征，并采用二分类算法根据训练集的脑电信号数据得到预测分类模型，大大提高了分类精度；
2.得到分类结果性对应的可信度：本发明的一种基于层次分析的脑电信号特征分类的二分类方法，对不同时间段或频率段分类结果的正确率与实际结果进行两两比较，构造判断矩阵a；对判断矩阵a进行一致性检验，确定不同时间段或者不同频率段判断矩阵a的权重w；用权重w对不同时间段或者频率段训练结果加权得到最后结果及其对应的可信度。有效的得到了分类结果的可信度。
附图说明
图1是本发明的方法流程图。
具体实施方式：
应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。下面结合附图与实施例对本发明作进一步说明。
实施例1：
正如背景技术所介绍的，现有技术中脑电信号在特征分类方面存在的分类依据特征单一或无法较好地综合多特征、分类精度不高、分类结果可信度未知的问题，提供一种基于层次分析的脑电信号特征分类的二分类方法，该方法可以提高分类精度和得到分类结果的可信度。
为了实现上述目的，本发明采用如下技术方案：
如图1所示，
一种基于层次分析的脑电信号特征分类的二分类方法，该方法的具体步骤包括：
(1)脑机接口系统采集运动想象脑电信号；
(2)将脑机接口系统采集的运动想象脑电信号进行格式转换，并进行数据预处理；
(3)采用特征提取算法分别提取不同时间段和频率段的脑电信号特征；
(4)将特征提取后的脑电信号征和其对应的标记按照合适比例分为训练集和测试集；
(5)根据训练集的脑电信号数据采用二分类算法得到预测分类模型；
(6)采用预测分类模型计算测试集的测试集分类结果，并且分别得出不同时间段和频率段结果的正确率；
(7)对不同时间段或频率段分类结果的正确率与实际结果进行两两比较，构造判断矩阵a；
(8)对判断矩阵a进行一致性检验，并判断是否通过一致性检验，若通过，进入步骤(9)，否则，返回步骤(7)；
(9)确定不同时间段或者不同频率段判断矩阵a的权重w；
(10)用权重w对不同时间段或者频率段训练结果加权得到最后结果及其对应的可信度。
所述步骤(1)中脑机接口系统采集运动想象脑电信号的具体步骤为：
脑机接口系统对用户进行不同心理作业时的运动想象脑电信号进行采集，分别包括想象左手运动时的运动想象脑电信号和想象右手运动时的运动想象脑电信号，并将采集的运动想象脑电信号进行存储。
在本实施例中，用户佩戴eeg电极帽进行不同心理作业的运动想象，想象左手运动和想象右手运动，采集相应的脑电信号为x＝[xl:xr]。想象左手运动时的运动想象脑电信号对应的结果为1，想象右手运动时的运动想象脑电信号对应的结果为-1。
所述步骤(2)中将脑机接口系统采集的运动想象脑电信号进行格式转换并进行数据预处理的具体步骤为：
(2-1)将采集的运动想象脑电信号由二维信号转换为三维信号：将步骤(1)中存储的脑机接口系统采集的运动想象脑电信号构建二维空间矩阵；
为了体现eeg信号的空间信息，采集到的eeg信号存储为x∈irn×m×t,其中m＝n×m是为了体现eeg空间信息而构建的二维空间矩阵，其依据是国际10－20电极位置图。若是标准的10－20电极位置图，则构建的二维空间矩阵m为一个5*5的二维矩阵，相应位置上的值为该电极上在不同时间点测量到的脑电信号电压值，为0的位置则表示该点没有测量信号，该点的值始终为0。t是采集的时间点，若采集时间为1秒，采集频率为128,则t＝128。
经过此种数据转换，测量到的脑电信号则为的三维信号，也可以看作一个的图像随着时间的变化而变化的录像。
在对运动想象脑电信号的分析训练中，建立能够描述脑电信号空间信息的二维矩阵，将采集的运动想象脑电信号由二维信号转换为三维信号，使得脑电信号描述过程中，不仅包含时间信息，也能够包含空间信息。
(2-2)对步骤(2-1)格式转换后的运动想象脑电信号进行空间滤波、频带的带通滤波、去除眼电和基线矫正预处理。
当大脑进行不同的运动想象时，不同电极位置上的脑电信号会随着时间的推移进行变化。同一种运动想象任务导致的这种变化是相似的，特征提取的任务就是针对某一种运动想象任务找到能够描述这种相似性的东西。
在本实施例中，所述步骤(3)中采用特征提取算法分别提取不同时间段和频率段的脑电信号特征的具体步骤为：
采用快速傅里叶变换(fft)、共空域模式(csp)或小波包变换(wpd)方法提取不同时间段和频率段脑电信号特征。
在本实施例中，所述步骤(4)中用适合的比率将数据随机分为训练集和测试集，70％的训练集和30％的测试集。
在本实施例中，所述步骤(5)中用支持向量机(svm)二分类算法根据70％的训练集得出预测分类模型。
在本实施例中，所述步骤(6)中分类结果分别为：
想象左手运动时的运动想象脑电信号对应的结果为1，想象右手运动时的运动想象脑电信号对应的结果为-1。
在本实施例中，所述步骤(7)中构造判断矩阵a的具体步骤为：
在决策者心目中，不同时间段和频率段由于正确率不同，他们对最后分类贡献度不同。
根据对不同时间段或频率段分类结果的正确率与实际结果进行两两比较得到的正确率大小的不同，引用数字1-9及其倒数作为标度构造出判断矩阵
a＝(aij)n×m。
在本实施例中，所述步骤(7)中，
如表1所示，
表1判断矩阵标度定义
数字1作为标度代表的含义为两个因素相比，具有相同的重要性，
数字3作为标度代表的含义为两个因素相比，前者比后者稍重要，
数字5作为标度代表的含义为两个因素相比，前者比后者明显重要，
数字7作为标度代表的含义为两个因素相比，前者比后者强烈重要，
数字9作为标度代表的含义为两个因素相比，前者比后者极端重要，
数字2、数字4、数字6和数字8分别表示上述相邻数字标度判断的中间值。
在本实施例中，所述步骤(8)中，对判断矩阵a进行一致性检验的具体步骤为：
(8-1)计算判断矩阵a的一致性指标ci：
其中，λmax为判断矩阵最大特征值，n为阶数；
(8-2)按照平均随机一致性指标表查找一致性指标ri；
表2平均随机一致性指标
(8-3)根据一致性指标ci与一致性指标ri的商计算一致性比例cr
当cr<0.10时，认为判断矩阵的一致性可以接受，否则判断矩阵需要修正，则回到步骤6，重新构造判断矩阵。
在本实施例中，所述步骤(9)中，采用几何平均法、算数平均法、特征向量法、或最小二乘法确定不同时间段或者不同频率段判断矩阵a的权重w。
1)几何平均法(方根法)
计算步骤：①a的元素按行相乘得一新的向量：
②将新的向量的每个分量开n次方；
③将所得向量归一化即为权重向量。
2)算术平均法(求和法)
由于判断矩阵a中的每一列都近似的反映了权值的分配情景，故可用全部向量的算术平局值来估计权向量，即：
计算步骤：①a的元素按列归一化，即求aij/
②将归一化后的各列相加；
③将相加后的向量除以n即得权重向量。
3)特征向量法
将权重向量w右乘权重比矩阵a，有：
aw＝λmaxw
同上，λmax为判断矩阵的最大特征值，存在且唯一，w的分量均为正分量。最后，将求得的权重向量作归一化处理即为所求。
4)最小二乘法
用拟合的方法确定权重向量，使残差和平方最小。即求解以下模型：
wi＞0,i＝1,2,…n.
在本实施例中，所述步骤(10)中，可信度分为非常好、较好、较差和非常差四个等级。用权重w对不同时间段或者频率段训练结果加权得到最后结果及可信度，加权的的结果对应可信度如表3：
表3
本发明的有益效果：
1.提高脑电信号在特征分类精度：本发明的一种基于层次分析的脑电信号特征分类的二分类方法，采用特征提取算法分别提取不同时间段和频率段的脑电信号特征，并采用二分类算法根据训练集的脑电信号数据得到预测分类模型，大大提高了分类精度；
2.得到分类结果性对应的可信度：本发明的一种基于层次分析的脑电信号特征分类的二分类方法，对不同时间段或频率段分类结果的正确率与实际结果进行两两比较，构造判断矩阵a；对判断矩阵a进行一致性检验，确定不同时间段或者不同频率段判断矩阵a的权重w；用权重w对不同时间段或者频率段训练结果加权得到最后结果及其对应的可信度。有效的得到了分类结果的可信度。
以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。