一种基于相位同步的多类别运动想象脑电信号分类方法技术

本发明专利技术涉及一种基于相位同步的多类别运动想象脑电信号分类方法。本发明专利技术首先利用相位锁定值分别计算训练样本和测试样本的相位同步特征，然后计算两者的相关系数并对其去平均取绝对值后按从大到小的顺序排列，接着根据排序后的相关系数进行脑电信号的粗分类，再将粗分类后的脑电信号进行细分类，此过程涉及共空域模式特征提取方法和线性判别分析分类方法。本发明专利技术包括脑电信号采集、数据预处理、滤波、相位同步特征相关系数计算、特征提取和分类以及分类准确率计算。分类结果表明：采用基于相位同步的脑电信号分类方法，得到了较好的分类效果，其中基于相位同步的脑电信号的粗分类可以高效率、高准确率地确定测试样本大致所属的类别，减少计算量。

【技术实现步骤摘要】
一种基于相位同步的多类别运动想象脑电信号分类方法
本专利技术属于生物特征识别领域中的脑电信号分类领域，具体涉及一种基于脑电信号相位同步获取相位同步特征相关系数并运用于共空域模式算法进行多类别运动想象脑电信号分类的方法。
技术介绍
对于神经肌肉系统功能严重丧失的患者来说,一种新的与外界交流的手段是急需的。脑机接口(BCI)正是这样一种技术,它允许使用者直接通过意念而不是传统的神经肌肉通道来实现与外界环境的自主交流。近年来脑机接口(BCI)研究已经逐渐从两类别的模式识别发展为多类别的模式识别，从而实现多种意念与外界环境的自主交流，提升脑机接口技术的实用性。如何有效地对多类别脑电信号进行分类，提升分类的正确率和分类的速率是目前运动想象多类别脑电信号分类研究中必须面临并解决的一个关键问题。而传统的多类别脑电信号分类方法如基于“一对一”共空域模式方法、基于Jacobi算法的近似联合对角化方法等有计算量较大等缺点。
技术实现思路
本专利技术针对现有技术的不足和改善需要，提供了一种基于相位同步的多类别运动想象脑电信号分类方法，该方法是利用相位锁定值分别计算训练样本和测试样本的相位同步特征并根据两者的相关系数依次将多对相应两个类别的脑电信号运用共空域模式算法通过线性判别分析进行多次分类来实现的。本专利技术所采用的技术方案：1.脑电信号采集：设计实验，使用多通道脑电采集设备采集实验过程中被试者的脑电信号，完成被试者信息录入和脑电信号采集。2.数据预处理：对原始脑电信号进行数据预处理，目的是减少伪迹的干扰，提高信噪比，从而提高特征提取的准确性。预处理过程包括带通滤波。3.滤波：创建滤波器，将预处理后的脑电信号滤波到计算利用相位同步特征进行分类所需的频段。4.相位同步特征相关系数计算：本专利技术采用相位锁定值(PLV)来表示各次实验所得脑电信号每两个通道之间的相位关系，具体的PLV计算公式如下：PLV＝|<exp(j{Φx(t)-Φy(t)})>|；其中Φx(t)和Φy(t)分别为t时刻脑电信号x(t)和y(t)的瞬时相位。采用希尔伯特变换来计算信号的相位值，信号x(t)的希尔伯特变换x～(t)定义如下：公式中的P为柯西主值。由此可以定义x(t)的解析信号为：其中Ax(t)和Φx(t)分别为信号x(t)的瞬时振幅以及瞬时相位。同样地，可以定义信号y(t)的解析信号，并计算瞬时相位Φy(t)。本专利技术中，根据单次实验的时间流程和样本，采样频率等采用合适的时间窗口来计算单次实验特定频段的PLV值。对于训练数据集而言，分别计算各类运动想象脑电信号PLV均值：假设训练数据集中存在Ni段属于第i类运动想象脑电信号样本的不重叠时间片段，则这Ni段片段的PLV均值计算如下：其中ΔΦ为信号x(t)和y(t)的瞬时相位差。假设本专利技术选定的脑电通道数为M，利用两两通道构建不同的通道对，计算某一类运动想象单次实验样本所有通道对的PLV值，此时得到一个M×M的上三角矩阵K；该矩阵不仅包含了不同脑电通道两两之间的相位关系，还包含了脑电通道的空间信息。接下来，对于训练数据集而言，按上述方法计算属于这一类运动想象的其他单次实验样本所有通道对的PLV值，从而得到这一类别运动想象数据集的PLV均值，用同样的方法计算其他各类别运动想象训练数据集的PLV均值，而测试数据集只需计算单次实验样本所有通道对的PLV值。最后，分别计算以上测试样本所有通道对PLV值所组成的矩阵A与各类别运动想象数据集PLV均值所组成的矩阵B的相位同步特征相关系数，具体的相关系数r计算公式如下：其中m,n为各矩阵的行列数，从而形成相关系数列向量E；其中ri是相关系数，i是运动想象脑电信号类别数。5.特征提取和分类：本专利技术中，特征提取和分类的步骤大致如下：第一步，基于相位同步特征相关系数的测试样本粗分类：将上述过程计算所得列向量E中的相关系数去平均后取其绝对值按从大到小排序形成列向量Q，即其中ri'＝abs(ri-mean(E)),abs(ri-mean(E))表示相关系数去平均后的绝对值，mean(E)表示相关系数的平均值，i是运动想象脑电信号类别数。若原始脑电信号类别数为偶数，则把Q中相关系数按顺序依次进行无重叠的两两组合，即每一个相关系数只能被组合一次，并把各两两组合相关系数所对应的两个类别脑电数据集作为训练样本进行CSP特征提取，并对测试样本用LDA分类方法进行粗分类,将其归为两个类别中的某一类；若原始脑电信号类别数为奇数，则把Q中相关系数按顺序依次进行无重叠的两两组合，此时Q中最小的相关系数不参与组合，即最小的相关系数所对应的脑电信号类别不参与组合，并把各两两组合相关系数所对应的两个类别脑电数据集作为训练样本进行CSP特征提取，并对测试样本用LDA分类方法进行粗分类,将其归为两个类别中的某一类。第二步，利用上述相关系数两两组合方法所得到的测试样本可能会属于的新的脑电信号类别，对测试样本进行细分类，按不同的情况具体操作如下：①若原始脑电信号类别数为偶数，且新的脑电信号类别数(即为下一次分类的原始脑电信号类别数)也为偶数，则把新的脑电信号类别按顺序依次进行无重叠的两两组合，即每一个新的脑电信号类别只能被组合一次，并把各两两组合所对应的两个类别脑电数据集作为训练样本进行CSP特征提取并对测试样本用LDA分类方法进行分类,将其归为两个类别中的某一类，形成新的脑电信号类别；若新的脑电信号类别数(即为下一次分类的原始脑电信号类别数)为奇数，则把新的脑电信号类别按顺序依次进行无重叠的两两组合，此时最后的脑电信号类别不参与组合，并把各两两组合所对应的两个类别脑电数据集作为训练样本进行CSP特征提取并对测试样本用LDA分类方法进行分类,将其归为两个类别中的某一类，形成新的脑电信号类别。②若原始脑电信号类别数为奇数，则把新的脑电信号类别和上一步骤中不参与组合的脑电信号类别组合成又一新的脑电信号类别(即为下一次分类的原始脑电信号类别)，然后根据这一新的脑电信号类别数的奇偶性按本步骤①中对新的脑电信号类别的处理方法进行操作。第三步，根据第二步中形成的原始脑电信号类别数和新的脑电信号类别数的奇偶性，循环第二步的操作，直到测试样本被最终分类。对上述步骤中涉及的CSP特征提取方法和LDA分类方法做如下介绍：1)特征提取：本专利技术中，采用共空域模式CSP方法对两两组合的两个类别脑电数据集进行特征提取。假设X1和X2分别为两类想象运动任务下维数为N×T的多导联脑电信号,其中N为EEG信号的通道数,T为每一个通道采集点数。X1和X2归一化后的协方差矩阵R1和R2分别为:其中XT表示矩阵X的转置,trace(X)为矩阵X对角线上元素的总和。对多次实验得到的空间协方差R1和R2取平均值,得到平均协方差矩阵和构造混合协方差矩阵R′,并对R′做主成分分析可得:其中U是混合协方差矩阵R′的特征向量组成的矩阵,∑是R′的对角线元为相应特征值的对角矩阵。定义白化矩阵P为:P＝∑-1/2UT；将矩阵和分别通过白化矩阵P进行白化变换:然后再对S1和S2做主成分分析:可以证明矩阵S1的特征向量矩阵U1和矩阵S2的特征向量矩阵U2是相等的,即U1＝U2＝V；同时两个特征值的对角阵∑1和∑2之和为单位矩阵,即：∑1+∑2＝I；可知,当矩本文档来自技高网...

【技术保护点】
一种基于相位同步的多类别运动想象脑电信号分类方法，其特征在于该方法包括以下步骤：步骤1、脑电信号采集：使用多通道脑电采集设备采集实验过程中被试者的脑电信号，完成被试者信息录入和脑电信号采集；步骤2、数据预处理：对原始脑电信号进行数据预处理，具体为带通滤波；步骤3、滤波：创建滤波器，将预处理后的脑电信号滤波到计算利用相位同步特征进行分类所需的频段；步骤4、相位同步特征相关系数计算：采用相位锁定值PLV来计算各次实验所得脑电信号每两个通道之间的相位关系，具体的PLV计算公式如下：PLV＝|<exp(j{Φx(t)‑Φy(t)})>|其中Φx(t)和Φy(t)分别为脑电信号x(t)和y(t)的瞬时相位；采用希尔伯特变换来计算信号的相位值，信号x(t)的希尔伯特变换x定义如下：x~(t)=1πP∫-∞+∞x(τ)t-τdτ]]>公式中的P为柯西主值；由此可以定义x(t)的解析信号为：Zx(t)=x(t)+jx~(t)=Ax(t)ejΦx(t)]]>其中Ax(t)和Φx(t)分别为信号x(t)的瞬时振幅以及瞬时相位；同样地，可以定义信号y(t)的解析信号，并计算瞬时相位Φy(t)；根据单次实验的时间流程和样本采样频率采用时间窗口来计算单次实验特定频段的PLV值；对于训练数据集而言，分别计算各类运动想象脑电信号PLV均值：假设训练数据集中存在Ni段属于第i类运动想象脑电信号样本的不重叠时间片段，则这Ni段片段的PLV均值计算如下：PLVavgi=1Ni|Σn=1Ni<exp(jΔΦ)>|]]>其中ΔΦ为信号x(t)和y(t)的瞬时相位差；设选定的脑电通道数为M，利用两两通道构建不同的通道对，计算某一类运动想象单次实验样本所有通道对的PLV值，此时得到一个M×M的上三角矩阵ΚK=0k12k13···k1M00k23···k2M············000···k(M-1)M00000]]>该矩阵不仅包含了不同脑电通道两两之间的相位关系，还包含了脑电通道的空间信息；接下来，对于训练数据集而言，按上述方法计算属于这一类运动想象的其他单次实验样本所有通道对的PLV值，从而得到这一类别运动想象数据集的PLV均值，用同样的方法计算其他各类别运动想象训练数据集的PLV均值，而测试数据集只需计算单次实验样本所有通道对的PLV值；最后，分别计算以上测试样本所有通道对PLV值所组成的矩阵Α与各类别运动想象数据集PLV均值所组成的矩阵Β的相关系数，具体的相关系数(r)计算公式如下：r=ΣmΣn(Amn-A‾)(Bmn-B‾)(ΣmΣn(Amn-A‾)2)(ΣmΣn(Bmn-B‾)2)]]>其中m,n为各矩阵的行列数，从而形成相关系数列向量RR=r1r2···ri]]>其中ri是相关系数，i是运动想象脑电信号类别数；步骤5、特征提取和分类：特征提取和分类的步骤大致如下：第一步，基于相位同步特征相关系数的测试样本粗分类：将上述过程计算所得列向量R中的相关系数去平均后取其绝对值按从大到小排序形成列向量C，即C=r1′r2′···ri′,r1′>r2′>...>ri′]]>其中ri'＝abs(ri‑mean(R)),abs()表示相关系数去平均后的绝对值，mean()表示相关系数的平均值,i是运动想象脑电信号类别数；若原始脑电信号类别数为偶数，则把C中相关系数按顺序依次进行无重叠的两两组合，即每一个相关系数只能被组合一次，并把各两两组合相关系数所对应的两个类别脑电数据集作为训练样本进行CSP特征提取并对测试样本用LDA分类方法进行粗分类,将其归为两个类别中的某一类；若原始脑电信号类别数为奇数，则把C中相关系数按顺序依次进行无重叠的两两组合，此时C中最小的相关系数不参与组合，即最小的相关系数所对应的脑电信号...

【技术特征摘要】
1.一种基于相位同步的多类别运动想象脑电信号分类方法，其特征在于该方法包括以下步骤：步骤1、脑电信号采集：使用多通道脑电采集设备采集实验过程中被试者的脑电信号，完成被试者信息录入和脑电信号采集；步骤2、数据预处理：对原始脑电信号进行数据预处理，具体为带通滤波；步骤3、滤波：创建滤波器，将预处理后的脑电信号滤波到计算利用相位同步特征进行分类所需的频段；步骤4、相位同步特征相关系数计算：采用相位锁定值PLV来计算各次实验所得脑电信号每两个通道之间的相位关系，具体的PLV计算公式如下：PLV＝|<exp(j{Φx(t)-Φy(t)})>|；其中Φx(t)和Φy(t)分别为t时刻脑电信号x(t)和y(t)的瞬时相位；采用希尔伯特变换来计算信号的相位值，信号x(t)的希尔伯特变换定义如下：公式中的P为柯西主值；由此可以定义x(t)的解析信号为：其中Ax(t)和Φx(t)分别为信号x(t)的瞬时振幅以及瞬时相位；同样地，可以定义信号y(t)的解析信号，并计算瞬时相位Φy(t)；根据单次实验的时间流程和样本采样频率，采用时间窗口来计算单次实验特定频段的PLV值；对于训练数据集而言，分别计算各类运动想象脑电信号PLV均值：假设训练数据集中存在Ni段属于第i类运动想象脑电信号样本的不重叠时间片段，则这Ni段片段的PLV均值计算如下：其中ΔΦ为信号x(t)和y(t)的瞬时相位差；设选定的脑电通道数为M，利用两两通道构建不同的通道对，计算某一类运动想象单次实验样本所有通道对的PLV值，此时得到一个M×M的上三角矩阵K；该矩阵不仅包含了不同脑电通道两两之间的相位关系，还包含了脑电通道的空间信息；接下来，对于训练数据集而言，按上述方法计算属于这一类运动想象的其他单次实验样本所有通道对的PLV值，从而得到这一类别运动想象数据集的PLV均值，用同样的方法计算其他各类别运动想象训练数据集的PLV均值，而测试数据集只需计算单次实验样本所有通道对的PLV值；最后，分别计算以上测试样本所有通道对PLV值所组成的矩阵A与各类别运动想象数据集PLV均值所组成的矩阵B的相关系数，具体的相关系数r计算公式如下：其中m,n为各矩阵的行列数，从而形成相关系数列向量E；其中ri是相关系数，i是运动想象脑电信号类别数；步骤5、特征提取和分类：特征提取和分类的步骤大致如下：第一步，基于相位同步特征相关系数的测试样本粗分类：将上述过程计算所得列向量E中的相关系数去平均后取其绝对值按从大到小排序形成列向量Q，即：其中ri'＝abs(ri-mean(E)),abs(ri-mean(E))表示相关系数去平均后的绝对值，mean(E)表示相关系数的平均值,i是运动想象脑电信号类别数；若原始脑电信号类别数为偶数，则把Q中相关系数按顺序依次进行无重叠的两两组合，即每一个相关系数只能被组合一次，并把各两两组合相关系数所对应的两个类别脑电数据集作为训练样本进行CSP特征提取，并对测试样本用LDA分类方法进行粗分类,将其归为两个类别中的某一类；若原始脑电信号类别数为奇数，则把Q中相关系数按顺序依次进行无重叠的两两组合，此时Q中最小的相关系数不参与组合，即最小的相关系数所对应的脑电信号类别不参与组合，并把各两两组合相关系数所对应的两个类别脑电数据集作为训练样本进行CSP特征提取，并对测试样本用LDA分类方法进行粗分类,将其归为两个类别中的某一类；第二步，利用上述相关系数两两组合方法所得到的测试样本可能会属于的新的脑电信号类别，对测试样本进行细分类，按不同的情况具体操作如下：①若原始脑电信号类别...

【专利技术属性】
技术研发人员：孔万增，徐飞鹏，周凌霄，徐思佳，任银芝，戴国骏，
申请(专利权)人：杭州电子科技大学，
类型：发明
国别省市：浙江;33