首页 / 专利库 / 信号处理 / 小波变换 / 离散小波变换 / 一种基于离散小波变换和多尺度熵的脑电情感识别方法

一种基于离散小波变换和多尺度熵的脑电情感识别方法

阅读:787发布:2020-05-16

专利汇可以提供一种基于离散小波变换和多尺度熵的脑电情感识别方法专利检索,专利查询,专利分析的服务。并且本 发明 公开一种基于离散 小波变换 和多尺度熵的脑电 情感识别 方法,先用小波变换分解重构出四种节律波,并计算出四种节律波的小波 能量 以及小波熵值,然后又运用多尺度熵 算法 ,提取出十四种尺度下的多尺度熵值,最后用主成分分析法对这些特征量进行 降维 后放入经 遗传算法 参数寻优的SVM分类器中识别。本发明同时考虑脑电情感数据集中数据的非线性和非平稳性特性,保证了该方法的分类 精度 、准确率和执行速度。,下面是一种基于离散小波变换和多尺度熵的脑电情感识别方法专利的具体信息内容。

1.一种基于离散小波变换和多尺度熵的脑电情感识别方法,其特征在于按照如下步骤依次进行:
步骤一:对脑电信号进行小波变换,并计算脑电信号的小波能量和小波熵;
对脑电信号x(t)用db4小波做四层小波分解可得:
其中Ai为近似分量其对应的近似系数为CAi,Dj为不同尺度下的细节分量其对应的细节系数为CDj,小波能量是该尺度下小波系数的平方和,设i为分解层数,那么求得总的小波能量为: 根据香浓熵的定义需取相对小波能量 那么可求
得小波熵:
步骤二:计算脑电信号的多尺度熵;
计算基于样本熵的多尺度熵之前,首先对于长度为M的时间序列{Y}={y1,y2,L,yM}进行粗粒化处理,得到新的序列为: 其中序列长度为 t
为尺度,然后按照如下过程来计算新序列的样本熵值;
1)对粗粒化后的序列{X}重构出一组m维矢量,其中i=1,2,ΛN-m+1,Xi=xi,xi+1,L,xi+m-1;
2)定义Xi与Xj间的距离di,j为对应元素切比夫距离的最大值,其中i≠j,计算公式为di,j=max[|Xi+k-Xj+k|];
3)给定阈值r,对每个在取值范围内的i,计算Xi与其余矢量Xj的距离di,j,统计di,j小于r的数目Lm(i),以及Lm(i)与距离总数N-m的比值,记为
4)计算 的平均值,记为
m+1
5)增加维数为m+1,重复以上过程,得到C ,当N为有限值时,那么该序列长度为N的样本熵为:
6)重复以上过程,最终得到不同尺度下的样本熵值;
步骤三:对上述特征量进行主成分分析法降维
经过上述方法进行一系列的特征提取可以得到4个小波能量的特征值,1个小波熵的特征值,以及14个多尺度熵的特征值,共计19个特征参数,如果将这19个参数都放入分类器中进行分类将导致情感识别和程序运行速度上存在较大的难度,所以本发明采用主成分分析法对19个特征量进行降维处理,最终降维后只得到3个特征特征维数;
步骤四:将特征样本放入支持向量机进行训练分类;
将样本进行随机排序,抽出一半的样本输入到SVM分类器中进行训练,然后把剩下的样本输入到已经训练好的分类器中进行测试,为了得到SVM最佳的惩罚因子c和核函数参数g,本发明利用遗传算法进行参数寻优,设置g和c的区间均为[0 100],最大遗传代数为200,交叉概率为0.4,变异概率为0.01。

说明书全文

一种基于离散小波变换和多尺度熵的脑电情感识别方法

技术领域

[0001] 本发明涉及生物信息技术领域以及数据挖掘领域,具体涉及一种基于离散小波变换(Discrete Wavelet Transform,DWT)与多尺度熵(Multi-Scale Entropy,MSE)相结合的脑电情感识别方法。

背景技术

[0002] 近年来,随着计算机领域的飞速发展人工智能的研究已经达到了较高平,因此其研究内容也在进一步的延伸。机器对人的情感进行识别是人工智能领域重要的应用延伸之一,这方面的研究会大大促进拟人机器人、人性化的商品设计以及情感沟通机器人等领域的发展,具有十分重要的意义。目前,在情感识别研究领域中,利用脑电信号(Electroencephalogram,EEG)来进行情感识别是主要研究方向之一,其特点是不易伪装,并且相对于其他的生理信号有较高的识别准确度。
[0003] 在脑电情感识别领域中,特征量的提取是影响识别准确率的关键因素,本发明中选用小波能量、小波熵和多尺度熵作为特征参数,其中小波能量可以反映出四个节律波的能量分布状况,小波熵可以反映出信号谱在各个空间上的无序程度、多尺度熵可以反映脑电信号在各种尺度下的相关程度和复杂程度,这些特征可以从各个层面来反映脑电信号的特点,能够达到很好的情感识别效果。

发明内容

[0004] 发明是为了解决现有技术存在的上述技术问题,提供一种可提高分类精度的基于小波变换和多尺度熵的脑电情感识别方法。
[0005] 为了达到上述的目的,本发明的技术方案为:
[0006] 一种基于离散小波变换和多尺度熵的脑电情感识别方法,其特征在于按照如下步骤依次进行:
[0007] 步骤一:对脑电信号进行小波变换,并计算脑电信号的小波能量和小波熵。
[0008] 对脑电信号x(t)用db4小波做四层小波分解可得:其中Ai为近似分量其对应的近似系数为CAi,Dj为不同
尺度下的细节分量其对应的细节系数为CDj。小波能量是该尺度下小波系数的平方和,设i为分解层数,那么求得总的小波能量为: 根据香浓熵的定义需
取相对小波能量 那么可求得小波熵:
[0009] 步骤二:计算脑电信号的多尺度熵。
[0010] 计算基于样本熵的多尺度熵之前,首先对于长度为M的时间序列{Y}={y1,y2,L,yM}进行粗粒化处理,得到新的序列为: 其中序列长度为 1≤k≤N,t为尺度。然后按照如下过程来计算新序列的样本熵值。
[0011] 1)对粗粒化后的序列{X}重构出一组m维矢量,其中i=1,2,ΛN-m+1,Xi=xi,xi+1,L,xi+m-1。
[0012] 2)定义Xi与Xj间的距离di,j为对应元素切比夫距离的最大值,其中i≠j,计算公式为:di,j=max[|Xi+k-Xj+k|]。
[0013] 3)给定阈值r,对每个在取值范围内的i,计算Xi与其余矢量Xj的距离di,j,统计di,j小于r的数目Lm(i),以及Lm(i)与距离总数N-m的比值,记为
[0014] 4)计算 的平均值,记为
[0015] 5)增加维数为m+1,重复以上过程,得到Cm+1。当N为有限值时,那么该序列长度为N的样本熵为:
[0016] 6)重复以上过程,最终得到不同尺度下的样本熵值。
[0017] 步骤三:对上述特征量进行主成分分析法降维
[0018] 经过上述方法进行一系列的特征提取可以得到4个小波能量的特征值,1个小波熵的特征值,以及14个多尺度熵的特征值,共计19个特征参数。如果将这19个参数都放入分类器中进行分类将导致情感识别和程序运行速度上存在较大的难度,所以本发明采用主成分分析法对19个特征量进行降维处理,最终降维后只得到3个特征特征维数。
[0019] 步骤四:将特征样本放入支持向量机进行训练分类。
[0020] 将样本进行随机排序,抽出一半的样本输入到SVM分类器中进行训练,然后把剩下的样本输入到已经训练好的分类器中进行测试。为了得到SVM最佳的惩罚因子c和核函数参数g,本发明利用遗传算法进行参数寻优,设置g和c的区间均为[0 100],最大遗传代数为200,交叉概率为0.4,变异概率为0.01。
[0021] 本发明公开一种基于离散小波变换和多尺度熵的脑电情感识别方法,先用小波变换分解重构出四种节律波,并计算出四种节律波的小波能量以及小波熵值,然后又运用多尺度熵算法,提取出十四种尺度下的多尺度熵值,最后用主成分分析法对这些特征量进行降维后放入经遗传算法参数寻优的SVM分类器中识别。本发明同时考虑脑电情感数据集中数据的非线性和非平稳性特性,保证了该方法的分类精度、准确率和执行速度。附图说明
[0022] 图1是本发明实施例流程图
[0023] 图2是不同尺度下四种情感的多尺度熵值变化曲线。
[0024] 图3是主成分累计贡献率的变化曲线。
[0025] 图4是遗传算法对分类器进行参数优化的适应度曲线。
[0026] 具体实时方式
[0027] 本发明的基于一种基于离散小波变换和多尺度熵的脑电情感识别方法,如图1所示,按照如下步骤依次进行:
[0028] 步骤一:对脑电信号进行小波变换,并计算脑电信号的小波能量和小波熵。
[0029] 对脑电信号x(t)用db4小波做四层小波分解可得:其中Ai为近似分量其对应的近似系数为CAi,Dj为不同
尺度下的细节分量其对应的细节系数为CDj。小波能量是该尺度下小波系数的平方和,设i为分解层数,那么求得总的小波能量为: 根据香浓熵的定义需
取相对小波能量 那么可求得小波熵:
[0030] 步骤二:计算脑电信号的多尺度熵。
[0031] 计算基于样本熵的多尺度熵之前,首先对于长度为M的时间序列{Y}={y1,y2,L,yM}进行粗粒化处理,得到新的序列为: 其中序列长度为 1≤k≤N,t为尺度。然后按照如下过程来计算新序列的样本熵值。
[0032] 1)对粗粒化后的序列{X}重构出一组m维矢量,其中i=1,2,ΛN-m+1,Xi=xi,xi+1,L,xi+m-1。
[0033] 2)定义Xi与Xj间的距离di,j为对应元素切比雪夫距离的最大值,其中i≠j,计算公式为:di,j=max[|Xi+k-Xj+k|]。
[0034] 3)给定阈值r,对每个在取值范围内的i,计算Xi与其余矢量Xj的距离di,j,统计di,j小于r的数目Lm(i),以及Lm(i)与距离总数N-m的比值,记为
[0035] 4)计算 的平均值,记为
[0036] 5)增加维数为m+1,重复以上过程,得到Cm+1。当N为有限值时,那么该序列长度为N的样本熵为:
[0037] 6)重复以上过程,最终得到不同尺度下的样本熵值。
[0038] 步骤三:对上述特征量进行主成分分析法降维。
[0039] 经过上述方法进行一系列的特征提取可以得到4个小波能量的特征值,1个小波熵的特征值,以及14个多尺度熵的特征值,共计19个特征参数。如果将这19个参数都放入分类器中进行分类将导致情感识别和程序运行速度上存在较大的难度,所以本发明采用主成分分析法对19个特征量进行降维处理,最终降维后只得到3个特征特征维数。
[0040] 步骤四:将特征样本放入支持向量机进行训练分类。
[0041] 将样本进行随机排序,抽出一半的样本输入到SVM分类器中进行训练,然后把剩下的样本输入到已经训练好的分类器中进行测试。为了得到SVM最佳的惩罚因子c和核函数参数g,本发明利用遗传算法进行参数寻优,设置g和c的区间均为[0 100],最大遗传代数为200,交叉概率为0.4,变异概率为0.01。
[0042] 本发明实施例进行情感识别的数据集是DEAP数据集。DEAP数据集是由Koelstra等人构建的一个多模态情感数据库,它包含了32个受试者观看40个时长为1分钟的音乐视频,记录了32个通道的脑电信号,以及8个通道的外围生理信号,采样频率为128Hz。每个受试者观看完一个视频后,都会对视频从唤醒度、效价、喜爱度、熟悉度等四个维度进行评分。
[0043] 本发明实施例仅对DEAP数据集中的脑电信号进行处理,数据集中包含32×40=1280个数据样本。为提高识别效率,每个样本仅选择了每个视频中间的9秒数据(第34~42秒的数据),因此每个样本共128×9×3=3456个数据点。
[0044] 目前最常见的情感分类方法是基于Lange的二维度情感分类模型,该模型由两个基本情感维度构成,分别是愉悦度(Valence)和唤醒度(Arouse)。Valence的值越大代表情绪越愉悦,变化范围由消极状态逐渐过渡到积极状态。Arouse的值越大代表越兴奋,变化范围由平静状态逐渐过渡到兴奋状态。本文利用该二维模型取愉快(HVHA)、放松(HVLA)、悲伤(LVLA)、愤怒(LVHA)这四种情绪来进行分类。
[0045] 图2为不同尺度下四种情感的多尺度熵值变化曲线,可以看出在整体趋势上其结果区别相对明显。为了不丢失信号中的有效特征,本发明选用1到14尺度下的多尺度熵值作为特征参数。图3为主成分累计贡献率的变化曲线,可以看出前3个成分的累计功效率达到了98.67%,他们已经包含了原特征参数的主要信息,所以选择用PCA处理后的前3个特征量作为分类器的输入量。图4为遗传算法对分类器进行参数优化的适应度曲线,可以看出前15周期收敛迅速,后面区域平缓进而实现了参数的寻优。最后将所有的分类结果汇总到表1中。
[0046] 表1本发明方法额区电极情感识别结果(%)
[0047]
[0048] 由表1可以看出该方法下额区的脑电信号均具有较高的识别率,其中最高为91.64%,最低为78.73%。从所有额区通道的识别率均值来看,悲伤(LVLA)的识别率是最高的均值为90.51%而愉快(HVHA)的识别率最低,为82.72%。从左右脑电极的平均准确率来看,左脑为83.85%右脑为86.64%,右脑的平均准确率较高,这也与神经生理学中右脑控制情绪左脑控制逻辑思维相吻合。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈