欢迎来到知嘟嘟! 联系电话:13336804447 卖家免费入驻,海量在线求购! 卖家免费入驻,海量在线求购!
知嘟嘟
我要发布
联系电话:13336804447
知嘟嘟经纪人
收藏
专利号: 2020101144570
申请人: 杭州电子科技大学
专利类型:发明专利
专利状态:已下证
专利领域: 医学或兽医学;卫生学
更新日期:2024-01-05
缴费截止日期: 暂无
价格&联系人
年费信息
委托购买

摘要:

权利要求书:

1.一种基于数据结构化分解的脑电信号分析方法,其特征在于:步骤1、采集受试者的脑电数据;对所得向量化一维脑电信号进行预处理,同时根据标准电极位置配置将一维脑电信号转成可以保留空间信息的二维脑电信号;

1‑1.对N个受试者均在M个不同的环境或时间下进行脑电数据采集,得到N·M组脑电数据;同一受试者在M次采集时大脑保持同样的被测任务类别;N个受试者在测试中对应的被测任务类别包含了所有的被测任务类别;

1‑2.对步骤1‑1所得的所有脑电数据进行预处理;

1‑3 .将各个一维的脑电数据 转化为二维脑电矩阵步骤2、对二维脑电数据进行分解;

2‑1.建立二维脑电数据集 二维脑电数据集 内的每个脑电数据Xi,j均为步骤1‑3中所得的其中一个二维脑电矩阵;脑电数据Xi,j由第i个受试者在第j种环境或时间下测得;

2‑2.将所有脑电信号Xi,j逐一分解为三个分量的组合,i=1,2,…,N,j=1,2,…,M,如式(1)所示;

式(1)中,Ci是表示脑电信号Xi,j的公共信息矩阵;Aj是脑电信号Xi,j的全局信息低秩矩阵;Ei,j是噪声矩阵;

2‑3.针对公共信息矩阵Ci、全局信息低秩矩阵Aj和噪声矩阵Ei,j建立目标函数 如式(2)所示;

式(2)中,||·||*表示核范式计算,||·||1表示l1范式计算;λi,j是控制稀疏度权重的标量;

为了将约束吸收到目标函数中,用增广的拉格朗日乘子将式(2)替换为式(3);

式(3)中,Yi,j为脑电数据Xi,j对应的拉格朗日乘子;μi,j是控制重构误差的标量;

步骤3、分别求解出全局信息低秩矩阵Aj、公共信息矩阵Ci、噪声矩阵Ei,j;

3‑1.使用奇异值阈值算法求解Aj;全局信息低秩矩阵Aj的求解问题表达为如式(4)所示;

式(4)中,||Aj||*表示核范式计算;||·||F表示F范式计算;<·,·>为两个向量的内积运算;

根据奇异值阈值算法,建立第一正交矩阵U、第二正交矩阵V和对角矩阵∑的表达式如式(5)所示;

根据式(4)和(5)重新建立全局信息低秩矩阵Aj的表达式如式(6)所示;

T

Aj=USτ(∑)V           (6)式(6)中,Sτ(∑)=sign(∑)·max(0,|∑|‑τ); sign(∑)表示矩阵∑的符号函数;max(·,·)表示两个数值中取最大值的运算;

3‑2.根据式(3)中Ei,j相关的部分,求解噪声矩阵Ei,j表达为如式(7)所示;

根据式(7)建立噪声矩阵Ei,j的表达式如式(8)所示式(8)中,

3‑3.求解公共信息矩阵Ci的问题表达为如式(9)所示;

式(9)中,

通过求解式(9)所述问题中关于Ci的偏导,并令其等于0,得到公共信息矩阵Ci的表达式如式(10)所示步骤4、构建样本子空间与识别方案;

利用各个受试者对应的公共信息矩阵Ci和全局信息低秩矩阵Aj分别构建各个受试者的子空间Si;i=1,2,…,N;子空间Si的集合表达式如式(11)所示;

式(11)中,ci、aj分别通过将公共信息矩阵Ci、全局信息低秩矩阵Aj向量化得到;wk表示第k个特征参数,k=1,2,…,M;

步骤5、使用步骤4得到的各个子空间Si,对被测人员的采集脑电数据时大脑执行的任务类别进行识别;

5‑1.对被测人员进行脑电数据的采集,并将被测人员的脑电数据根据步骤1中的方法预处理并转化为二维脑电矩阵Z;

5‑2.构建被测人员的子空间Sz;子空间Sz的集合表达式如式(12)所示;

式(12)中,z通过将二维脑电矩阵Z向量化得到;

5‑3.分别对比子空间Sz与步骤4得到的M个子空间Si的相似度;取与子空间Sz相似度最高的那个子空间Si作为目标子空间;取目标子空间对应的受试者对应的任务类别判断为被测人员的在脑电采集时大脑的任务类别。

2.根据权利要求1所述的一种基于数据结构化分解的脑电信号分析方法,其特征在于:一组脑电数据包括一次采集中获得的多个时刻的脑电数据。

3.根据权利要求1所述的一种基于数据结构化分解的脑电信号分析方法,其特征在于:步骤1‑1中所述的被测任务类别为使用者的情绪类别,包括积极、中性、消极三种情绪类别。

4.根据权利要求1所述的一种基于数据结构化分解的脑电信号分析方法,其特征在于:步骤1‑2中的预处理过程如下:首先使用带通滤波,从脑电数据中提取2‑40Hz频段的脑电成分;再计算滤波后的脑电数据各个通道不同时刻的脑电数据平均值,并将每个脑电通道的各个时刻脑电数据分别减去对应通道的平均值。