基于SCSP-LDA的脑电信号特征提取与分类方法技术

本发明专利技术提出一种基于SCSP‑LDA的脑电信号特征提取与分类方法，该方法首先，将脑电数据进行特征值分解，对其进行搜索并筛选得到新的特征空间，然后利用CSP提取其特征，最后，利用LDA对特征提取后的数据进行特征优化和数据分类，实现运动想象EEG信号解码。本发明专利技术将CSP算法转化为广义特征值的求解问题，并由广义特征值的求解问题结合稀疏搜索算法寻找最佳空间滤波器，引出一种稀疏共空间模式(Sparse Common Spatial Pattern,SCSP)算法，该SCSP是算法可以有效地提取特征最为明显的通道，从而实现通道稀疏，再结合特征分类算法LDA，可以完整地实现EEG的解码。

Feature extraction and classification of EEG based on scsp-lda

【技术实现步骤摘要】
基于SCSP-LDA的脑电信号特征提取与分类方法
本专利技术涉及生物信号处理和模式识别领域，尤其是涉及一种脑-机接口中基于稀疏共空间模式(SparseCommonSpatialPatterns,SCSP)-LDA的EEG特征提取与分类方法。
技术介绍
脑-机接口(brain-computerinterface,BCI)是使人在不需依赖外围神经系统和肌肉，只需采集大脑发出的信号便可实现与计算机或者其他设备通信的系统。其作为一种人脑直接与外围设备进行交流的途径，能够为那些因为脑部损伤而患有运动障碍的病人带来与外界进行交流的解决办法。脑电信号(electroencephalography，EEG)作为脑-机接口的主要研究方向，对EEG进行解码便是研究的重点。共空间模式(CommonSpatialPatterns，CSP)算法最早由Fukunage等人提出，之后Romeser和他的同事把它用在BCI中。CSP作为一种有效的空间滤波和特征提取算法，可以从多通道脑电数据中提取每个类的空间分布成分。但是，在一般的EEG采集中，为了获取更多的信息，往往会采集多通道的脑电数据，这样必然会存在多个通道脑电信号出现信号重叠的情况。
技术实现思路
本专利技术要解决的是脑电数据多通道采集时存在的多通道脑电信号重叠的问题。为解决上述技术问题，本专利技术提出一种基于SCSP-LDA的EEG特征提取与分类方法。首先利用稀疏共空间模式进行EEG的特征提取，并对EEG通道进行稀疏，然后通过LDA对提取到的特征进行降维并分类。该方法包括如下步骤：步骤1、将CSP算法转化为广义特征值的求解问题，具体步骤如下：步骤11、将CSP算法视为一种在广义瑞利商基础上的算法，其相对应的表达式可写为：式中，w为空间滤波矩阵，T为求转置，Xi(i＝1,2)表示将原始样本数据经过去均值预处理后得到的样本矩阵，Ci(i＝1,2)为样本数据的协方差矩阵；步骤12、在J(w)式中分子分母同乘以一非零常数k，其值不发生改变，即：J(kw)＝J(w)步骤13、利用所述步骤12得到的式子可以对所述步骤11中的式子进行简化，通过对常数k的调整，进而改变w的值，令wTC2w＝1，将对J(w)求极值的问题化简为对wTC1w求极值的问题；从而能够把以上问题转变为求拉格朗日算子的表达式，即：L(λ,w)＝wTC1w-λ(wTC2w-1)式中，λ为拉格朗日乘子；步骤14、对所述步骤13中的表达式中的w求偏导，同时令左边为0，故而有：通过以上步骤，CSP算法便转化为广义特征值的求解问题；其中，设则w为M关于特征值λ的广义特征向量；步骤2、由广义特征值的求解问题结合稀疏搜索算法寻找最佳空间滤波器，得到SCSP算法，用于所述SCSP算法对EEG的通道进行稀疏和特征提取；步骤3、提取新特征空间：步骤31、利用SCSP算法提取到的特征空间W，对原始数据进行滤波处理，得到滤波后的矩阵Z，即：ZN×T＝WN×N·EN×T式中，Z为滤波后的矩阵，W为利用SCSP算法求得的最优特征空间，E为原始数据中每一次测量得到的一组数据，其中N表示脑电数据的通道数，T表示每一次测量得到的数据点数；步骤32、从滤波后的矩阵ZN×T提取特征空间矩阵fp：式中，fp为特征空间矩阵，Zp由经过滤波的数组Z的前后各取k行构成；步骤4、利用LDA算法对所述特征空间矩阵fp进行特征降维与分类。优选的，所述步骤2包括如下具体步骤：步骤21、选择通道稀疏模式并设置稀疏度spl；步骤22、输入两类EEG协方差矩阵A、B；步骤23、运行选定的搜索模式；步骤24、定义通道存储数组ix；步骤25、以最大化两类EEG方差的差异，即以求解式A(ix,ix)X＝λB(ix,ix)X中特征值的大小为评价标准，迭代搜索的方式遍历所有通道后筛选出最佳通道，并且后一次循环是在前一次得到的最优解的基础上进行，以此保证每次得到结果都是最优解，直到通道数达到设定的稀疏度，即得到最佳通道组合，输出最终选择出的通道组合ix；步骤26、利用步骤25得到的最佳通道组合ix计算出最佳特征子空间。优选的，所述步骤21中所述选择通道稀疏模式包括：向前选择(forwardselection，FS)模式、向后消去(backwardelimination，BE)模式和贪婪搜索(greedysearch，GS)模式。优选的，所述步骤4包括如下具体步骤：步骤41、将两类EEG数据带入到LDA算法中，求得LDA的投影矩阵A*，有如下公式：式中，A*为LDA的投影矩阵，m1、m2分别为两类数据的均值向量，Sw两类数据的类内距离之和；步骤42、将由SCSP算法提取到的特征空间矩阵fp向法向量A*投影，实现特征降维，即：y＝(A*)T*fp式中，A*为LDA的投影矩阵，fp为特征空间矩阵；步骤43、根据得分y并设置合适的阈值实现对提取到的特征的分类。与现有技术相比，本专利技术具有如下有益效果：在本专利技术将CSP算法转化为广义特征值的求解问题，并由广义特征值的求解问题结合稀疏搜索算法寻找最佳空间滤波器，引出一种稀疏共空间模式(SparseCommonSpatialPattern,SCSP)算法，该SCSP是算法可以有效地提取特征最为明显的通道，从而实现通道稀疏，再结合特征分类算法LDA，可以完整地实现EEG的解码。附图说明图1为EEG特征提取与特征分类过程示意图；图2为EEG稀疏搜索算法程序流程图；以及图3为EEG分类结果比较柱状图。具体实施方式以下，参照附图对本专利技术的实施方式进行说明。本专利技术的一种基于SCSP-LDA的脑电信号特征提取与分类方法，首先，将脑电数据进行特征值分解，对其进行搜索并筛选得到新的特征空间，然后利用CSP提取其特征。最后，利用LDA对特征提取后的数据进行特征优化和数据分类，实现运动想象EEG信号解码。其总体流程图如图1所示，该方法的实现步骤如下：步骤1、将CSP算法转化为广义特征值的求解问题，具体步骤如下：步骤11、将CSP算法视为一种在广义瑞利商基础上的算法，其相对应的表达式可写为：式中，w为空间滤波矩阵，T为求转置，Xi(i＝1,2)表示将原始样本数据经过去均值预处理后得到的样本矩阵，Ci(i＝1,2)为样本数据的协方差矩阵；步骤12、在J(w)式中分子分母同乘以一非零常数k，其值不发生改变，即：J(kw)＝J(w)步骤13、利用步骤12得到的式子可以对步骤11中的式子进行简化，通过对常数k的调整，进而改变w的值，令wTC2w＝1，将对J(w)求极值的问题化简为对wTC1w求极值的问题；从而能够把以上问题转变为求拉格朗日算子的表达式，即：L(λ,w)＝wTC1本文档来自技高网...

【技术保护点】
1.一种基于SCSP-LDA的脑电信号特征提取与分类方法，其特征在于，该方法包括如下步骤：/n步骤1、将CSP算法转化为广义特征值的求解问题，具体步骤如下：/n步骤11、将CSP算法视为一种在广义瑞利商基础上的算法，其相对应的表达式J可写为：/n

【技术特征摘要】
1.一种基于SCSP-LDA的脑电信号特征提取与分类方法，其特征在于，该方法包括如下步骤：
步骤1、将CSP算法转化为广义特征值的求解问题，具体步骤如下：
步骤11、将CSP算法视为一种在广义瑞利商基础上的算法，其相对应的表达式J可写为：



式中，w为空间滤波矩阵，T为求转置，Xi(i＝1,2)表示将原始样本数据经过去均值预处理后得到的样本矩阵，Ci(i＝1,2)为样本数据的协方差矩阵；
步骤12、所述步骤11中所述表达式的分子分母同乘以一非零常数k，其值不发生改变，即：
J(kw)＝J(w)；
步骤13、利用所述步骤12得到的式子对所述步骤11中的式子进行简化，通过对常数k的调整，进而改变w的值，令wTC2w＝1，将对所述表达式J求极值的问题化简为对wTC1w求极值的问题；从而能够把以上问题转变为求拉格朗日算子的表达式，即：
L(λ,w)＝wTC1w-λ(wTC2w-1)
式中，λ为拉格朗日乘子；
步骤14、对所述步骤13中的表达式中的w求偏导，同时令左边为0，故而有：



通过以上步骤，CSP算法便转化为广义特征值的求解问题；其中，设则w为M关于特征值λ的广义特征向量；
步骤2、由广义特征值的求解问题结合稀疏搜索算法寻找最佳空间滤波器，得到SCSP算法，用于所述SCSP算法对EEG的通道进行稀疏和特征提取；
步骤3、提取新特征空间：
步骤31、利用SCSP算法提取到的特征空间W，对原始数据进行滤波处理，得到滤波后的矩阵Z，即：
ZN×T＝WN×N·EN×T
式中，Z为滤波后的矩阵，W为利用SCSP算法求得的最优特征空间，E为原始数据中每一次测量得到的一组数据，其中N表示脑电数据的通道数，T表示每一次测量得到的数据点数；
步骤32、从滤波后的矩阵ZN×T提取特征空间矩阵fp：

...

【专利技术属性】
技术研发人员：付荣荣，田永胜，王世伟，
申请(专利权)人：燕山大学，
类型：发明
国别省市：河北;13