CN115062711A - 基于多源域适应字典学习和稀疏表示的脑电情感识别方法 - Google Patents
基于多源域适应字典学习和稀疏表示的脑电情感识别方法 Download PDFInfo
- Publication number
- CN115062711A CN115062711A CN202210716694.3A CN202210716694A CN115062711A CN 115062711 A CN115062711 A CN 115062711A CN 202210716694 A CN202210716694 A CN 202210716694A CN 115062711 A CN115062711 A CN 115062711A
- Authority
- CN
- China
- Prior art keywords
- class
- sparse
- source domain
- dictionary
- sparse representation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 69
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 20
- 230000008909 emotion recognition Effects 0.000 title claims abstract description 19
- 230000006870 function Effects 0.000 claims abstract description 36
- 238000005457 optimization Methods 0.000 claims abstract description 17
- 238000012360 testing method Methods 0.000 claims abstract description 14
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- 239000011159 matrix material Substances 0.000 claims description 29
- 238000012549 training Methods 0.000 claims description 16
- 239000000126 substance Substances 0.000 claims description 13
- 239000013598 vector Substances 0.000 claims description 11
- 238000013507 mapping Methods 0.000 claims description 7
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 238000009795 derivation Methods 0.000 claims description 2
- 238000010626 work up procedure Methods 0.000 claims description 2
- 239000007788 liquid Substances 0.000 claims 1
- 238000013508 migration Methods 0.000 abstract description 5
- 230000005012 migration Effects 0.000 abstract description 5
- 238000001182 laser chemical vapour deposition Methods 0.000 abstract 1
- 230000037007 arousal Effects 0.000 description 17
- PIZHFBODNLEQBL-UHFFFAOYSA-N 2,2-diethoxy-1-phenylethanone Chemical compound CCOC(OCC)C(=O)C1=CC=CC=C1 PIZHFBODNLEQBL-UHFFFAOYSA-N 0.000 description 15
- 238000002474 experimental method Methods 0.000 description 11
- 238000012706 support-vector machine Methods 0.000 description 5
- 238000009826 distribution Methods 0.000 description 4
- 230000006978 adaptation Effects 0.000 description 3
- 238000010276 construction Methods 0.000 description 2
- 230000002996 emotional effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/16—Devices for psychotechnics; Testing reaction times ; Devices for evaluating the psychological state
- A61B5/165—Evaluating the state of mind, e.g. depression, anxiety
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
- A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Physics & Mathematics (AREA)
- Psychiatry (AREA)
- Public Health (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Fuzzy Systems (AREA)
- Signal Processing (AREA)
- Mathematical Physics (AREA)
- Evolutionary Computation (AREA)
- Physiology (AREA)
- Child & Adolescent Psychology (AREA)
- Developmental Disabilities (AREA)
- Educational Technology (AREA)
- Hospice & Palliative Care (AREA)
- Psychology (AREA)
- Social Psychology (AREA)
- Image Analysis (AREA)
Abstract
本发明涉及EEG脑电技术领域，尤其涉及基于多源域适应字典学习和稀疏表示的脑电情感识别方法，包括：使用LCSVD算法得到字典和稀疏的初始值；计算源域和目标域上类内、类间稀疏重建误差；根据最小化、最大化类内稀疏重建误差构建判别项，计算字典学习和稀疏表示的目标函数；采用参数迭代优化方法对目标函数求解；通过未标注的测试样本，计算测试样本类别。本发明共有字典的学习准则是最小化类内稀疏重建误差和最大化类间稀疏重建误差，保证学习到的稀疏表示具有识别能力；为避免负迁移的发生，每个源域对应一个领域适应权重，且可以在自适应学习中得到权重最佳值；模型参数的求解通过参数交替优化的方法，所有的参数可同时达到最优解。
Description
技术领域
本发明涉及EEG脑电技术领域，尤其涉及基于多源域适应字典学习和稀疏表示的脑电情感识别方法。
背景技术
现有的基于EEG的情感识别需要在训练阶段花费大量的带标签数据，一种比较简单而直接的方法是重复利用之前采集的EEG数据训练一个通用分类器，而不考虑个体间的差异。但是，传统的机器学习算法建立在训练数据和测试数据独立同分布的假设之上；而这个假设对情感BCI往往难以成立，这是因为EEG信号具有非线性和非平稳性的固有特性，针对相同情绪状态从不同人群中提取的特征集，如脑电通道的谱带功率，并没有表现出很强的相关性；因此，脑电通道谱带功率及其导数的训练和测试数据往往具有不同的分布，也就是说，当同一分类器应用于其他受试者的脑电数据或从各种其他数据集中提取的数据时，性能会显著降低。
另外，现有方法中只考虑单个源领域，将所有相关的辅助脑电图数据视为一个源领域，没有涉及多源域学习问题。
发明内容
针对现有算法的不足，本发明将多个源域和目标域的EEG样本投影到一个共享的投影子空间中，在共享子空间中学习一个域不变字典，域不变字典的学习准则是最小化类内稀疏重建误差和最大化类间稀疏重建误差，学习到的稀疏特征具有很强的识别能力，并满足类内分散度小和类间分离大的特性；本发明为每个源域分配一个领域适应权重，且在模型学习中自适应学习到每个源领域权重最佳值，能有效避免负迁移的发生；多源域领域适应字典学习和稀疏表示(MDA-DLSR)目标函数的求解采取参数交替优化的方法，保证所有的参数同时达到最优解。
本发明所采用的技术方案是：一种基于多源域适应字典学习和稀疏表示的脑电情感识别方法包括以下步骤：
首先，构建字典学习和稀疏表示的训练样本，使用LC-SVD算法得到字典和稀疏表示矩阵的初始值；
DLSR是为无监督学习提出，其优化问题是在词典空间中最小化原始信号与重构信号之间的重构误差，设训练集X＝[x1,x2,...,xn]∈Rm×n，其中，m和n分别表示数据的维数和样本数，通过优化经验函数
DLSR使用字典D∈Rm×K的有限原子对数据进行分解，A＝[a1,...,an]∈RK×n是稀疏系数，l()表示损失函数，K是字典原子数；在DLSR中，原始信号和重构信号之间的重构误差均方是最常见的损失函数，通常由l1范数将稀疏性引入目标函数，因此，这一思想表示成：
为了避免D的值任意大，从而避免ai的值任意小，需要对字典原子或稀疏系数进行额外的约束，以限制其l2范数的值。式(1)中添加该约束后表示为：
其中，T0表示字典原子中非零个数阈值，显然，原始DLSR是无监督学习，其在优化问题中没有考虑类别的标签信息；因此，在监督学习中，需要将标签类别信息加入到目标式中，由此，产生了一系列有监督DLSR方法，标签一致性LC-SVD算法使用线性分类器来构建模型：
其中，W∈RC×K是分类参数,T0表示字典原子的阈值，α是正则化参数，D表示字典，A表示稀疏表示矩阵。
然后，使用投影矩阵将源域投影到低维子空间，计算源域上的类内稀疏重建误差和类间稀疏重建误差；
然后，使用投影矩阵将目标域投影到低维子空间，计算目标域上的类内稀疏重建误差和类间稀疏重建误差；
然后，根据最小化类内稀疏重建误差和最大化类间稀疏重建误差的构建判别项，并计算字典学习和稀疏表示的目标函数，采用参数迭代优化方法对目标函数求解；
最后，通过未标注的测试样本，根据得到的最优投影矩阵和字典，计算测试样本类别。
进一步的，具体步骤包括：
步骤一、构建字典学习和稀疏表示的训练样本，从DEAP数据集采集EEG信号，构成源域样本集，设sC个源域样本集
si＝{1,2,...,sC}和一个目标域
且每个领域的数据都不重叠，均带有类别标签，sC个源域和目标域的特征空间一致，但它们的边缘概率分布和条件概率分布可以不同，第si个源域Xsi和目标域Xt中的样本分别表示为
和
样本的个数分别是nsi和nt，全部源域样本的个数为ns，
步骤二、使用投影矩阵Psi∈Rm×d，将第si个源域Xsi投影到低维子空间，基于字典学习和稀疏表示理论计算类内稀疏重建误差：
其中，
表示Xsi中第c类的第j个样本，
是
的稀疏表示，D表示字典，
是源域Xsi上的类内散度，
函数返回
的稀疏表示系数的K维映射向量；
通过下式计算得到
其中，
函数返回
的稀疏表示系数的K维映射向量，
的稀疏表示
中第c类系数保留原值的，其它值均为零，tr(·)表示矩阵的秩操作。
投影空间中第si个源域Xsi上，基于字典学习和稀疏表示理论计算类间稀疏重建误差：
其中，
是源域Xsi上的类间散度，
函数返回
的稀疏表示系数的K维映射向量；
通过下式计算得到
其中，
函数返回
的稀疏表示系数的K维映射向量，与
不同的是，其对应的
的稀疏表示
中第c类系数为零的，其它值保留。
步骤三、使用投影矩阵Pt∈Rm×d，将目标域Xt投影到低维子空间，目标域Xt上的类内稀疏重建误差
和类间稀疏重建误差
计算式如下：
其中，
是目标域Xt上的类内散度，通过下式计算得到
其中，
是目标域Xt上的类间散度，
在共享子空间中通过公共字典来学习多源域和目标域EEG数据的共有判别知识，一方面，多源域EEG信号的稀疏表示和分类器在子空间中是独立的；另一方面，域不变字典用于在各个域之间建立潜在的联系，将判别信息从多源域传递到目标域；
步骤四、根据最小化类内稀疏重建误差和最大化类间稀疏重建误差构建判别项保证了模型的判别能力，为了有效利用多个源领数据的判别知识，同时避免负迁移，MDA-DLSR赋予每个源域自适应权重，基于以上思想，MDA-DLSR的目标函数表示为：
其中，向量λ＝[λ1,λ2,…,λsC]是权重向量，元素λsi表示第si个源领域在目标函数中的重要程度，r(r>1)是平衡指数；在传统的权重求解问题中，λ的值需人工设置，常通过网络搜索法获得最优，但MDA-DLSR中λ的值作为一个模型参数自适应得到最优值；
为简化目标函数的计算，将式(10)中两项的分子式和分母式进行合并，同时，定义矩阵
式(10)可以重写为：
进一步的，公式(11)是一个非凸问题，涉及三个参数{P,D,λ}，采用参数迭代优化方法进行求解，
进一步的，参数迭代优化方法进行求解，具体包括：
1、固定参数{D,λ}，求解P，为了避免过度拟合，将一个正则化项加入到式(11)的分母中，根据跟踪比优化(trace ratio optimization,TRO)策略，对于投影矩阵P，必有一个最大值q*，使得：
其中，μ为一个正数；
因此得到：
tr(PTAP)-q*tr(PT(B+μI)PT)≤0 (13)
经过整理，式(13)进一步写成以下形式：
tr(PT(A-q*(B+μI)P)≤0 (14)
为求解式(14)，定义以下关于q的函数f(q)：
f(q)具有两个性质：(1)f(q)是一个递减函数；(2)f(q)＝0当且仅当q＝q*，也就是说，q的最优值总是存在的，通过迭代交替更新P和q能得到相应的最优解；因此，对(15)进行q的一阶求导，可得f'(q)＝-tr(P(q)T(B+μI)P(q))，令f'(q)＝0，得到q的极值，考虑目标函数关于矩阵P的约束PTP＝I，P的最优值可通过下式求解：
显然，式(16)是一个特征值分解问题：
(A-qB-qμI)P＝γP (17)其最优解是求解式(11)得到的关于P的前d个最大特征值。
2、固定参数{P,λ}，求解D，目标函数关于每个类别的子字典表示为：
根据公式(5)和(7)，以及矩阵A和B的定义，矩阵Ac和Bc分别对应矩阵A和B在第c类样本上的运算结果；
子字典Dc采用梯度上升法进行求解，Dc＝Dc+αJ'(Dc)，其中，α是学习率，J'(Dc)则可由下式计算得到：
3、固定参数{P,D}，求解λ，在式(10)上引入拉格朗日乘子β并忽略常数项，得到下式：
λ的极值的必要条件为
同时考虑
得到：
由式(21)得到：
将上式化简，并代入式(10)，得λsi的解析解：
步骤五、对于一个未标注的测试样本z，根据得到的最优投影矩阵Pt和字典D，使用下式计算其在第j类上的类别：
其中，
是第j类子字典Dj的伪逆；
最后，使用投票法来得到样本z的类标签，即
其中，Δj表示第j类的票数。
本发明的有益效果：
1、通过多个源领域和目标领域寻找一个投影子空间，在子空间中通过共有字典建立多个源域和目标领域之间的桥梁，为充分利用源领域数据的识别能力，共有字典的学习准则是最小化类内稀疏重建误差和最大化类间稀疏重建误差；同时，通过源领域自适应地权重学习可以有效避免负迁移的发生；
2、在真实脑电情感识别数据集DEAP上，对比本发明方法与现有技术的四种方法，结果表明本发明方法在DEAP数据集上的Arousal和Valence识别精度最高。
附图说明
图1是本发明的基于多源域适应字典学习和稀疏表示的脑电情感识别方法流程图；
图2是本发明的方法与现有方法在32名个体上的Arousal识别精度比较；
图3是本发明方法与现有方法在32名个体上的Valence识别精度比较；
图4是本发明的参数m和d在DEAP数据集上的Arousal识别精度；
图5是本发明的参数m和d在DEAP数据集上的Valence识别精度。
具体实施方式
下面结合附图和实施例对本发明作进一步说明，此图为简化的示意图，仅以示意方式说明本发明的基本结构，因此其仅显示与本发明有关的构成。
使用DEAP数据集验证本发明方法有效性，DEAP数据集来自英国伦敦玛丽皇后大学等4所大学，记录了32名志愿者在观看40段音乐视频时的EEG信号和视频信号，对每位志愿者进行了40次数据采集，每次采集过程中播放一段特定的63秒视频，这些视频由参与者根据唤醒(Arousal)和效价(Valence)等指标进行评分，评分范围为1-9；实验中采用的类别标签设定为：将评分划分成2个带有阈值5的二元分类问题：高/低Arousal和高/低Valence(低：≤5，高：>5)；EEG信号按照国际10/20系统放置32个电极采集得到；数据预处理时将EEG数据下采样至128Hz，去除信号的伪影并删除3秒基线数据，并应用4.0-45.0Hz的频率实施带通滤波；在EEG信号的众多特征提取方法中，微分熵(differential entropy,DE)被认为能反映EEG特征的连续随机变量的复杂性；对于固定长度的EEG片段，DE特征可以等价于特定频带内的对数能量谱；实验使用256采样点的短时傅里叶变换和1s的非重叠窗口来提取EEG信号的五个频带(δ:1-3Hz,θ:4-7Hz,α:8-13Hz,β:14-30Hz,γ:31-50Hz)；然后计算每个频带的DE特征，由于每个频带信号有62个通道，一个样本可提取310维的DE特征。
如图1所示，一种基于多源域适应字典学习和稀疏表示的脑电情感识别方法包括以下步骤：
首先，构建字典学习和稀疏表示的训练样本，使用LC-SVD算法得到字典和稀疏表示矩阵的初始值；
然后，使用投影矩阵将源域投影到低维子空间，计算源域上的类内稀疏重建误差和类间稀疏重建误差；
然后，使用投影矩阵将目标域投影到低维子空间，计算目标域上的类内稀疏重建误差和类间稀疏重建误差；
然后，根据最小化类内稀疏重建误差和最大化类间稀疏重建误差的构建判别项，并计算字典学习和稀疏表示的目标函数；采用参数迭代优化方法对目标函数求解；
最后，通过一个未标注的测试样本，根据得到的最优投影矩阵和字典，计算测试样本类别。
进一步的，具体步骤包括：
步骤一、构建字典学习和稀疏表示的训练样本：
步骤二、使用投影矩阵Psi∈Rm×d，将第si个源域Xsi投影到低维子空间，基于字典学习和稀疏表示理论计算类内稀疏重建误差：
步骤三、使用投影矩阵Pt∈Rm×d，将目标域Xt投影到低维子空间，目标域Xt上的类内稀疏重建误差
和类间稀疏重建误差
步骤四、根据最小化类内稀疏重建误差和最大化类间稀疏重建误差的构建判别项保证了模型的判别能力，MDA-DLSR赋予每个源域自适应权重，构建MDA-DLSR的目标函数；
步骤五、对于一个未标注的测试样本z，根据得到的最优投影矩阵Pt和字典D，使用下式计算其在第j类上的类别。
MDA-DLSR方法的训练步骤表示如下：
为验证所提方法的性能，实验与两类方法进行了比较：一类是基线方法：高斯核支持向量机(support vector machine,SVM)和标签一致LC-KSVD方法；另一类是领域适应方法：自适应子空间特征匹配(adaptive subspace feature matching,ASFM)方法,最大独立域自适应(maximum independence domain adaptation,MIDA)方法，稳健主成分分析(robust principal component analysis,RPCA)；具体的参数设置如下：
高斯核的核参在网格{10-3,10-2,...,103}中搜索，正则化参数在网格{2-6,2-5,...,26}中搜索；LC-KSVD的字典原子数设置为训练样本个数的1/3；MIDA的投影空间维度在网格{20,30,...,100}中搜索；ASFM的阈值参数设置为0.45；RPCA的平衡参数设置为样本数和特征数最大值的平方根；本发明方法的子空间维度是通过搜索网格{20,30,...,100}来确定的；每个类的原子数在{10,15,20,25,30,35}中选择，实验中源领域和目标领域的设置策略采取留一法，即数据集中的每个个体均可作为目标领域，其它的个体均作为源域。因此，DEAP数据集上进行了32轮实验，其中31名个体作为多个源域，因为训练集样本较多，实验中在每名个体样本中随机选择1/2数量的EEG数据进行训练，目标域中随机选择20个样本作为训练集，剩余的目标域样本作为测试集；这一实验过程进行了10次，记录了每个方法的分类精度。所有算法均在MATLAB 2019a中实现。
对比实验：
实验比较了MDA-DLSR方法在DEAP数据集上的Arousal和Valence的识别精度；各方法在32名个体上的平均识别精度如表1所示。为了充分显示各方法在每名个体上的具体结果，图2-3显示了各方法在32名个体上的Arousal和Valence的识别精度；从实验结果可以看到：
1)首先，基线分类方法SVM和LC-KSVD在DEAP数据集上的跨个体EEG情感识别任务中无法获得令人满意的Arousal和Valence识别精度；因为本质上它们并不是为了解决领域适应问题，SVM和LC-KSVD在模型将所有的源域数据和目标域训练数据混合在一起，此时，由于不同领域数据分布上的差别使得源领域对于目标域上数据识别的辅助效果是有限的。
2)其次，在对比的领域适应方法中，本发明方法比使用其它对比的方法表现更好，主要原因是从多个源领域习得的共享字典可以使用更多的判别知识来辅佐目标域上分类器的建立，而且每个源领域可以自适应地学习到对应的权重，可以有效防止负迁移的发生，保证迁移学习的效果，此外，MDA-DLSR方法的参数迭代学习策略保证了所有参数同时达到最优解。
3)从表1结果可以看出，MDA-DLSR方法不仅Arousal和Valence的识别精度是最优的，其对应的标准差也是较小的，说明方法具有较好的稳定性。而所有方法中，非领域适应方法SVM和LC-KSVD的标准差较大，说明这2个方法在处理跨个体EEG情感识别任务易发生过拟合的现象。
4)从图2-3结果可以看出，32名个体的Arousal和Valence识别精度有较大的区别，有些个体的识别精度达到80％以上，有些个体的识别精度略高于50％，这也许是因为在EEG信号采样中受到外界和自身干扰，获得的EEG数据质量较低。
表1 DEAP数据集的Arousal和Valence的识别精度(标准差)比较
参数分析
本发明方法中需要寻优的参数有平衡指数γ、子空间的维数m和字典原子个数d；其中，平衡指数γ的寻优范围是{1,1.5,...,3}，子空间的维数m的寻优范围是{10,20,...,100}和字典原子个数d的寻优范围是{30,40,...,120}；表2列出了参数γ在DEAP数据集上的Arousal和Valence识别精度；图4-5分别显示了参数m和d在DEAP数据集上的Arousal和Valence识别精度，从实验结果可以看出：
(1)平衡指数γ的作用是调节各源域在目标函数的权重，因为各源域的权重是介于[0,1]之间的实数，γ的值越小，各源域的权重比就越大；相反，γ的值越小，各源域的权重比就越接近；根据表2中的结果，平衡指数γ对应的Arousal和Valence识别精度的变化是温和的，实验中可设置γ＝2。
表2参数γ在DEAP数据集上的Arousal和Valence识别精度
(2)投影空间的维数m和字典原子个数d对MDA-DLSR方法的性能起到决定性作用，不同m和d值对应不同的Arousal和Valence识别精度；从图4-5可以看出，MDA-DLSR方法在m和d较小时就可以得到较高的Arousal和Valence，当m大于50且d大于80时，MDA-DLSR方法的性能可以达到稳定；根据图4-5的结果，实验可以设置m和d值分别设置为60和80。
本发明方法为多个源领域和目标领域寻找一个投影子空间，在子空间中通过共有字典建立多个源域和目标领域之间的桥梁，为充分利用源领域数据的识别能力，共有字典的学习准则是最小化类内稀疏重建误差和最大化类间稀疏重建误差。同时，通过源领域自适应地权重学习可以有效避免负迁移的发生，在真实脑电情感识别数据集DEAP上的实验证明了本发明方法的有效性。
以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。
Claims (7)
1.基于多源域适应字典学习和稀疏表示的脑电情感识别方法，其特征在于，包括以下步骤：
步骤一、构建字典学习和稀疏表示的训练样本，使用LC-SVD算法得到字典和稀疏表示矩阵的初始值；
步骤二、使用投影矩阵将源域投影到低维子空间，计算源域上的类内稀疏重建误差和类间稀疏重建误差；
步骤三、使用投影矩阵将目标域投影到低维子空间，计算目标域上的类内稀疏重建误差和类间稀疏重建误差；
步骤四、根据最小化类内稀疏重建误差和最大化类间稀疏重建误差构建判别项，并计算字典学习和稀疏表示的目标函数，采用参数迭代优化方法对目标函数求解；
步骤五、通过未标注的测试样本，根据得到的最优投影矩阵和字典，计算测试样本类别。
2.根据权利要求1所述的基于多源域适应字典学习和稀疏表示的脑电情感识别方法，其特征在于，步骤一的具体步骤如下：
构建字典学习和稀疏表示的训练样本，设sC个源域样本集
si＝{1,2,...,sC}和一个目标域
第si个源域Xsi和目标域Xt中的样本分别表示为
和
样本的个数分别是nsi和nt，全部源域样本的个数为ns，
3.根据权利要求1所述的基于多源域适应字典学习和稀疏表示的脑电情感识别方法，其特征在于，步骤二的具体步骤如下：
使用投影矩阵Psi∈Rm×d，将第si个源域Xsi投影到低维子空间，计算类内稀疏重建误差：
其中，
表示Xsi中第c类的第j个样本，
是
的稀疏表示，D表示字典，
是源域Xsi上的类内散度，
函数返回
的稀疏表示系数的K维映射向量；
通过下式计算得到
其中，
函数返回
的稀疏表示系数的K维映射向量，
的稀疏表示
中第c类系数保留原值，其它值均为零，tr(·)表示矩阵的秩操作；
投影空间中第si个源域Xsi上，计算类间稀疏重建误差：
其中，
是源域Xsi上的类间散度，
函数返回
的稀疏表示系数的K维映射向量；
通过下式计算得到
其中，
函数返回
的稀疏表示系数的K维映射向量。
4.根据权利要求3所述的基于多源域适应字典学习和稀疏表示的脑电情感识别方法，其特征在于，步骤三的具体步骤如下：
使用投影矩阵Pt∈Rm×d，将目标域Xt投影到低维子空间，目标域Xt上的类内稀疏重建误差
和类间稀疏重建误差
计算式如下：
其中，
是目标域Xt上的类内散度；
通过下式计算得到
其中，
是目标域Xt上的类间散度，
5.根据权利要求4所述的基于多源域适应字典学习和稀疏表示的脑电情感识别方法，其特征在于，所述步骤四的具体步骤如下：
根据最小化类内稀疏重建误差和最大化类间稀疏重建误差构建判别项，并赋予每个源域自适应权重，计算字典学习和稀疏表示的目标函数，公式为：
其中，向量λ＝[λ1,λ2,…,λsC]是权重向量，元素λsi表示第si个源领域在目标函数中的重要程度，r(r＞1)是平衡指数；
为简化目标函数的计算，将式(10)中两项的分子式和分母式进行合并，同时，定义矩阵
式(10)可以重写为：
公式(11)涉及三个参数{P,D,λ}，采用参数迭代优化方法进行求解。
6.根据权利要求5所述的基于多源域适应字典学习和稀疏表示的脑电情感识别方法，其特征在于，参数迭代优化方法进行求解，具体包括：
1、固定参数{D,λ}，求解P，将一个正则化项加入到式(11)的分母中，根据跟踪比优化策略，对于投影矩阵P，必有一个最大值q*，使得：
其中，μ为一个正数；
得到：
tr(PTAP)-q*tr(PT(B+μI)PT)≤0 (13)
经过整理，式(13)进一步写成以下形式：
tr(PT(A-q*(B+μI)P)≤0 (14)
为求解式(14)，定义以下关于q的函数f(q)：
通过迭代交替更新P和q能得到相应的最优解；因此，对(15)进行q的一阶求导，可得f'(q)＝-tr(P(q)T(B+μI)P(q))，令f'(q)＝0，得到q的极值，考虑目标函数关于矩阵P的约束PTP＝I，P的最优值可通过下式求解：
显然，式(16)是一个特征值分解问题：
(A-qB-qμI)P＝γP, (17)其最优解是求解式(11)得到的关于P的前d个最大特征值；
2、固定参数{P,λ}，求解D，目标函数关于每个类别的子字典表示为：
根据公式(5)和(7)，以及矩阵A和B的定义，矩阵Ac和Bc分别对应矩阵A和B在第c类样本上的运算结果；
子字典Dc采用梯度上升法进行求解，Dc＝Dc+αJ'(Dc)，其中α是学习率，J'(Dc)则可由下式计算得到：
3、固定参数{P,D}，求解λ，在式(10)上引入拉格朗日乘子β并忽略常数项，得到下式：
λ的极值的必要条件为
同时考虑
得到：
由式(21)得到：
将上式化简，并代入式(10)，得λsi的解析解：
7.根据权利要求6所述的基于多源域适应字典学习和稀疏表示的脑电情感识别方法，其特征在于，计算测试样本类别的公式为：
其中，
是第j类子字典Dj的伪逆。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210716694.3A CN115062711A (zh) | 2022-06-23 | 2022-06-23 | 基于多源域适应字典学习和稀疏表示的脑电情感识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210716694.3A CN115062711A (zh) | 2022-06-23 | 2022-06-23 | 基于多源域适应字典学习和稀疏表示的脑电情感识别方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115062711A true CN115062711A (zh) | 2022-09-16 |
Family
ID=83202988
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210716694.3A Withdrawn CN115062711A (zh) | 2022-06-23 | 2022-06-23 | 基于多源域适应字典学习和稀疏表示的脑电情感识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115062711A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116702018A (zh) * | 2023-04-19 | 2023-09-05 | 东南大学 | 基于ga-pdpl算法的跨被试脑电情感识别方法及装置 |
-
2022
- 2022-06-23 CN CN202210716694.3A patent/CN115062711A/zh not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116702018A (zh) * | 2023-04-19 | 2023-09-05 | 东南大学 | 基于ga-pdpl算法的跨被试脑电情感识别方法及装置 |
| CN116702018B (zh) * | 2023-04-19 | 2024-03-01 | 东南大学 | 基于ga-pdpl算法的跨被试脑电情感识别方法及装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Que et al. | Back to the future: Radial basis function networks revisited | |
| Chen et al. | Model Metric Co-Learning for Time Series Classification. | |
| CN110472649B (zh) | 基于多尺度分析和集成树模型的脑电情感分类方法及系统 | |
| CN112101152A (zh) | 一种脑电情感识别方法、系统、计算机设备、可穿戴设备 | |
| CN113191225B (zh) | 一种基于图注意力网络的情绪脑电识别方法及系统 | |
| CN112465069B (zh) | 一种基于多尺度卷积核cnn的脑电情绪分类方法 | |
| Köster et al. | A two-layer model of natural stimuli estimated with score matching | |
| CN111310656A (zh) | 基于多线性主成分分析的单次运动想象脑电信号识别方法 | |
| CN115062711A (zh) | 基于多源域适应字典学习和稀疏表示的脑电情感识别方法 | |
| CN115105076A (zh) | 基于动态卷积残差多源迁移的脑电情绪识别方法及系统 | |
| CN115221969A (zh) | 基于emd数据增强和并行scn的运动想象脑电信号识别方法 | |
| CN113191206B (zh) | 基于黎曼特征迁移的脑磁图信号分类方法、装置及介质 | |
| CN113592030A (zh) | 基于复值奇异谱分析的图像检索方法和系统 | |
| Kaleem et al. | Empirical mode decomposition based sparse dictionary learning with application to signal classification | |
| Jiang et al. | Picture-induced EEG signal classification based on CVC emotion recognition system | |
| Puri et al. | Wavelet packet sub-band based classification of alcoholic and controlled state EEG signals | |
| CN116919422A (zh) | 基于图卷积的多特征情感脑电识别模型建立方法及装置 | |
| CN110991364A (zh) | 脑电信号分类方法及系统 | |
| CN115758118A (zh) | 一种基于脑电互信息的多源流形嵌入特征选择方法 | |
| CN111008590B (zh) | 基于多维特征提取无监督eeg信号分类方法 | |
| Nuo et al. | A Fusion Transfer Learning Method of Motor Imagery EEG Signals Based on Riemannian Space | |
| Nguyen et al. | Machine Learning in Classification of Parkinson’s Disease Using Electroencephalogram with Simon’s Conflict | |
| Pham et al. | Feature-reduction fuzzy co-clustering algorithm for hyperspectral image segmentation | |
| Han et al. | EEG signals classification based on wavelet packet and ensemble Extreme Learning Machine | |
| Zhao et al. | Multilinear generalization of common spatial pattern |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WW01 | Invention patent application withdrawn after publication | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20220916 |