一种运动想象脑电信号分类方法及系统技术方案

本发明专利技术提供一种运动想象脑电信号分类方法及系统，涉及脑电信号分类识别研究技术领域。该运动想象脑电信号分类方法，包括以下步骤，S1：采集获取被试者不同运动想象状态下脑电波信号，S2：对采集的脑电波信号进行滤波、降噪预处理，S3：对预处理后的信号进行特征提取，构建联合特征矩阵，S4：对信号进行降维处理后，再分别进行训练和测试，获取最优分类器模型。通过对可编辑数字滤波器与期望信号产生的误差信号进行梯度计算和修正，改善了信号输出质量，构建时频域特征向量，避免了经验模态分解中的“模态混叠”问题，提高了脑电信号识别准确率，更好地满足处理需求。更好地满足处理需求。更好地满足处理需求。

【技术实现步骤摘要】
一种运动想象脑电信号分类方法及系统


[0001]本专利技术涉及脑电信号分类识别研究
，具体为一种运动想象脑电信号分类方法及系统。

技术介绍

[0002]脑机接口技术(BrainComputerInterface，BCI)是一种新的人机接口方式，通过分析提取人类的进行思维活动时脑电信号，实现脑意识与人的交流以及对外部设备的控制，在医疗康复、军事作战、教育娱乐和智能家居等领域有着广阔的发展前景，常见的BCI系统有基于稳态视觉诱发电位、基于P300成分和基于运动想象的BC系统等。
[0003]运动想象是一种脑机接口经典的范式，使用者可自主发出控制指令。其生理学基础是身体运动可在大脑运动感觉区产生mu(8
‑
12Hz)和beta(13
‑
30Hz)节律，出现事件相关同步化和事件相关去同步化现象。相较于其他范式，具有不需要外界刺激诱发的优点。基于运动想象的脑机接口系统可用于脑卒中康复、控制电动轮椅、文本拼写器和假肢等设备，除了能简化正常人与机器交互的流程之外，还能够帮助残障人士提高生活自理能力。运动想象脑电信号的分类准确率因受信号特性和用户认知程度等多种因素影响，难以达到较高值，而且不同用户的分类准确率差异较大。运动想象系统要实时地应用在实际问题中，就必须要求其拥有很高的分类精度以及很快的反应速度。
[0004]常见的分类方法有支持向量机、线性判别分析和K最近邻等，研究人员根据分类的样本特征和种类选择最合适的分类算法，通过寻优算法确定分类器超参数，训练出最优分类器，最终实现对运动想象脑电信号的识别分类，然而，大部分的研究都只是对脑电信号的特征进行单一的时域或频域的特征提取，会出现经验模态分解中的“模态混叠”问题，对脑电信号的识别准确率不高。

技术实现思路

[0005]针对现有技术的不足，本专利技术提供了一种运动想象脑电信号分类方法及系统，通过矩阵乘以一个随机矩阵来实现联合特征矩阵的降维处理，避免经验模态分解中的“模态混叠”问题，对脑电信号的识别准确率更高。
[0006]为实现以上目的，本专利技术通过以下技术方案予以实现：一种运动想象脑电信号分类方法，包括以下步骤：
[0007]S1：采集获取被试者不同运动想象状态下脑电波信号；
[0008]S2：对采集的脑电波信号进行滤波、降噪预处理；
[0009]S3：对预处理后的信号进行特征提取，构建联合特征矩阵；
[0010]S4：对信号进行降维处理后，再分别进行训练和测试，获取最优分类器模型。
[0011]优选的，采用对可编辑数字滤波器与期望信号产生的误差信号进行梯度计算和修正，输出最终降噪后的信号，包括以下步骤：
[0012]S21：脑电波信号X(n)输送至可编辑数字滤波器中，经过算法处理和补偿后得到滤
波器输出信号y(n)，同时期望输出信号Y(n)也输入至可编辑数字滤波器中；
[0013]S22：计算Y(n)和y(n)的差值得到误差信号e(n)，随后将每次的输入信号和误差信号同时输入到梯度计算与系数修正模块，计算出下一轮滤波所使用的系数；
[0014]S23：通过可编程数字滤波器不断更新滤波器的系数，用滤波器模拟得到的反馈量和实际反馈量之间的误差，并最终输出得到最终输出信号L(n)。
[0015]优选的，对每导联信号进行互补集合经验模态分解，并利用皮尔逊相关系数筛选每个导联脑电信号的本征模态函数分量，构建时频域特征向量，包括以下步骤：
[0016]S31：信号经预处理后，对每个导联信号进行互补集合经验模态分解，并利用皮尔逊相关系数筛选每个导联脑电信号的本征模态函数分量；
[0017]S32：提取筛选后每个导联信号的时域和频域特征向量；
[0018]S33：通过共同空间模式算法获得空间滤波器，提取每个导联脑电信号的空域特征向量，构建联合特征矩阵。
[0019]优选的，通过矩阵乘以一个随机矩阵来实现联合特征矩阵的降维处理，包括以下步骤：
[0020]记描述该运动想象特征向量m的矩阵为D∈IRd
×
p，其中，d为关键点模块的特征向量长度，p为关键点模块的数量；
[0021]将维数为d的向量投影到一个维数为n的子空间中，其中n＜＜d；
[0022]通过矩阵D乘以一个随机矩阵R来实现特征向量矩阵的降维处理，将特征向量矩阵维数降低为n
×
p：
[0023]Y＝RD
[0024]其中，Y∈IRn
×
p为降维后的特征向量矩阵，R为随机投影矩阵，R∈IRn
×
d，均值为0，方差为1。
[0025]优选的，所述步骤S4中，通过训练和测试获取最优分类器模型的方法为将处理后的信号分为训练组和测试组，训练组将脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维和类别直方图计算后，得到训练数据的分类器模型最优参数，测试组将脑电信号作为训练数据进行格式转换、特征点提取、特征向量计算、降维和类别直方图计算后，利用训练组中的分类器模型最优参数进行稀疏表达，得到最优分类器模型，随后将测试组数据输入训练好的分类器模型中进行分类，完成对多类运动想象的识别分类。
[0026]具体的，对训练数据进行格式转换的方法为将采集的运动想象脑电信号由二维信号转换为三维信号：将第一脑电信号存储模块中存储的脑机接口系统采集的运动想象脑电信号X∈IRn
×
m
×
T构建二维空间矩阵。
[0027]对特征点提取和特征向量计算的方法为：(1)关键点检测，第一脑电信号存储模块中存储的脑机接口系统采集的运动想象脑电信号X∈IRn
×
m
×
T表示为X＝[X1,X2,...,XT]，对于X中的每个值利用Harris关键点检测器检测关键点，针对检测到的所有关键点提取大小为(η
×
η
×
T)的模块，(2)特征值提取，针对(1)中的模块所有值去均值，针对模块里面的每一个像素沿着时间方向计算其方差M2，偏态M3和峰度M4，组成与模块对应的时间矩阵Mr,r＝{2,3,4}，其中Mr＝[mij],i,j＝1,2,...,η，(3)特征向量计算，将(2)中的每一个时间矩阵转化为一个长度为η2的向量，时间矩阵Mr,r＝{2,3,4}中r的三个值组成特征向量m∈IRd，其中d＝3η2为关键点模块的特征向量长度。
[0028]一种运动想象脑电信号分类系统，包括：
[0029]脑电波信号采集模块，用于采集获取被试者不同运动想象状态下脑电波信号；
[0030]预处理模块，用于对采集的脑电波信号进行滤波、降噪预处理；
[0031]特征提取模块，对预处理后的信号进行特征提取，构建联合特征矩阵；
[0032]特征处理模块，对信号进行降维处理，所述特征处理模块包括训练模块和测试模块，用于对脑电信号数据进行训练和测试，获取结果数据，所述特征处理模块还本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种运动想象脑电信号分类方法，其特征在于，包括以下步骤：S1：采集获取被试者不同运动想象状态下脑电波信号；S2：对采集的脑电波信号进行滤波、降噪预处理；S3：对预处理后的信号进行特征提取，构建联合特征矩阵；S4：对信号进行降维处理后，再分别进行训练和测试，获取最优分类器模型。2.根据权利要求1所述的一种运动想象脑电信号分类方法，其特征在于，采用对可编辑数字滤波器与期望信号产生的误差信号进行梯度计算和修正，输出最终降噪后的信号，包括以下步骤：S21：脑电波信号X(n)输送至可编辑数字滤波器中，经过算法处理和补偿后得到滤波器输出信号y(n)，同时期望输出信号Y(n)也输入至可编辑数字滤波器中；S22：计算Y(n)和y(n)的差值得到误差信号e(n)，随后将每次的输入信号和误差信号同时输入到梯度计算与系数修正模块，计算出下一轮滤波所使用的系数；S23：通过可编程数字滤波器不断更新滤波器的系数，用滤波器模拟得到的反馈量和实际反馈量之间的误差，并最终输出得到最终输出信号L(n)。3.根据权利要求1所述的一种运动想象脑电信号分类方法，其特征在于，对每导联信号进行互补集合经验模态分解，并利用皮尔逊相关系数筛选每个导联脑电信号的本征模态函数分量，构建时频域特征向量，包括以下步骤：S31：信号经预处理后，对每个导联信号进行互补集合经验模态分解，并利用皮尔逊相关系数筛选每个导联脑电信号的本征模态函数分量；S32：提取筛选后每个导联信号的时域和频域特征向量；S33：通过共同空间模式算法获得空间滤波器，提取每个导联脑电信号的空域特征向量，构建联合特征矩阵。4.根据权利要求1所述的一种运动想象脑电信号分类方法，其特征在于，通过矩阵乘以一个随机矩阵来实现联合特征矩阵的降维处理，包括以下步骤：(1)记描述该运动想象特征向量m的矩阵为D∈IRd...

【专利技术属性】
技术研发人员：裴昕，
申请(专利权)人：东北电力大学，
类型：发明
国别省市：