发明名称 基于脑电特征的情绪状态识别方法
摘要 一种基于脑电特征的情绪状态识别方法,有:数据采集阶段,是在国际情绪图片诱发条件下,提取被试在不同愉悦度等级图片诱发下的64导脑电数据;数据预处理阶段,对采集到的64导脑电数据进行改变参考电位、降采样、带通滤波、去除眼电四个步骤的处理;特征提取阶段,是将预处理后的信号用共空间模式算法进行滤波后提取时域特征;特征识别,使用支持向量机分类器对特征进行识别,将不同情绪状态区分开。本发明尝试的OVR共空间模式算法,可去除背景信号的干扰,用于多类情绪诱发脑电的信号增强。去除背景信号干扰后,不同类别的情绪脑电间差异增大,用时域方差特征进行识别,被试者的识别正确率较为理想,可以准确地区分开不同愉悦度的情绪。
申请公布号 CN102715911B 申请公布日期 2014.05.28
申请号 CN201210199052.7 申请日期 2012.06.15
申请人 天津大学 发明人 张迪;明东;陈龙;李南南;柯余峰;许敏鹏;綦宏志;万柏坤
分类号 A61B5/16(2006.01)I;A61B5/0476(2006.01)I 主分类号 A61B5/16(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 杜文茹
主权项 1.一种基于脑电特征的情绪状态识别方法,其特征在于,包括如下阶段:(1)数据采集阶段,所述的数据采集是在国际情绪图片诱发条件下,提取被试在不同愉悦度等级图片诱发下的64导脑电数据,包括进行如下步骤:1)实验设计:从国际情绪图片库中选择了45张图片,按愉悦度范围分为3个等级,每个等级15张图片;等级1、2、3分别对应代表消极、中性、积极情绪图片;包含45次任务,每次刺激任务时间为14秒,有三个时段,即图片呈现时段A持续6s、休息时段B持续6s和提醒时段C持续2s;2)脑电采集:采集EEG的实验程序是在stim2平台下设计的,脑电采集仪器使用的是scan4.3,使用的脑电电极帽是10-20标准64导联;(2)数据预处理阶段,对采集到的64导脑电数据进行改变参考电位、降采样、带通滤波、去除眼电四个步骤的处理;(3)特征提取阶段,所述的特征提取阶段是将预处理后的信号用共空间模式算法进行滤波后提取时域特征,包括如下步骤:1)共空间模式滤波:三类任务原始的脑电信号为64×384维的矩阵:设为<img file="2012101990527100001DEST_PATH_IMAGE001.GIF" wi="20" he="22" />,<img file="283648DEST_PATH_IMAGE002.GIF" wi="21" he="22" />和<img file="2012101990527100001DEST_PATH_IMAGE003.GIF" wi="21" he="22" />,其中64为导联通道数,384为每个通道在刺激任务下的前3秒的数据点数,则每一类信号的归一化协方差矩阵分别为:<img file="532095DEST_PATH_IMAGE004.GIF" wi="266" he="69" />这里的<img file="DEST_PATH_IMAGE005.GIF" wi="30" he="25" />表示对<img file="816446DEST_PATH_IMAGE006.GIF" wi="18" he="17" />的转置,<i>trace</i>表示矩阵的迹;构造合成的空间协方差矩阵为:<img file="DEST_PATH_IMAGE007.GIF" wi="126" he="29" /><img file="166656DEST_PATH_IMAGE008.GIF" wi="21" he="22" />可分解为:<img file="DEST_PATH_IMAGE009.GIF" wi="94" he="31" /><img file="381605DEST_PATH_IMAGE010.GIF" wi="25" he="28" />和<img file="DEST_PATH_IMAGE011.GIF" wi="23" he="24" />分别为其特征矢量矩阵和特征值对角阵,白化变换可使方差均匀化,白化矩阵为:<img file="687822DEST_PATH_IMAGE012.GIF" wi="107" he="34" />为了更好地阐述这个算法,先考虑如何得到愉悦度等级为1的诱发脑电的空间滤波器,设:<img file="DEST_PATH_IMAGE013.GIF" wi="94" he="31" />令:<img file="877495DEST_PATH_IMAGE014.GIF" wi="102" he="32" /><img file="DEST_PATH_IMAGE015.GIF" wi="102" he="31" />这里<img file="904881DEST_PATH_IMAGE016.GIF" wi="24" he="29" />是对应<img file="DEST_PATH_IMAGE017.GIF" wi="22" he="29" />的白化矩阵;如果<img file="408675DEST_PATH_IMAGE018.GIF" wi="20" he="25" />可以分解成:<img file="DEST_PATH_IMAGE019.GIF" wi="94" he="31" />则<img file="897294DEST_PATH_IMAGE020.GIF" wi="19" he="42" />可被分解为:<img file="DEST_PATH_IMAGE021.GIF" wi="94" he="31" />且有:<img file="195551DEST_PATH_IMAGE022.GIF" wi="89" he="33" />将式<img file="769621DEST_PATH_IMAGE014.GIF" wi="102" he="32" />,<img file="139422DEST_PATH_IMAGE015.GIF" wi="102" he="31" />,<img file="967701DEST_PATH_IMAGE019.GIF" wi="94" he="31" />,<img file="764756DEST_PATH_IMAGE021.GIF" wi="94" he="31" />,<img file="560542DEST_PATH_IMAGE022.GIF" wi="89" he="33" />综合,即得到:<img file="DEST_PATH_IMAGE023.GIF" wi="373" he="35" />这里的<img file="671718DEST_PATH_IMAGE024.GIF" wi="23" he="28" />的列向量是矩阵<img file="DEST_PATH_IMAGE025.GIF" wi="22" he="28" />的特征向量,可以看出,经变换后信号协方差矩阵的特征矢量对应的特征值之和为1,所以在<img file="872279DEST_PATH_IMAGE001.GIF" wi="20" he="22" />方差值最大的方向上其余的信号<img file="840235DEST_PATH_IMAGE002.GIF" wi="21" he="22" />和<img file="874050DEST_PATH_IMAGE003.GIF" wi="21" he="22" />的方差值就很小,所以选取<img file="585654DEST_PATH_IMAGE026.GIF" wi="22" he="27" />中与最大的特征值对应的特征向量定义为<img file="DEST_PATH_IMAGE027.GIF" wi="27" he="29" />,则对应于<img file="637793DEST_PATH_IMAGE001.GIF" wi="20" he="22" />的空间滤波器,即信号的投影方向为:<img file="2012101990527100001DEST_PATH_IMAGE002.GIF" wi="156" he="39" />这里<img file="297761DEST_PATH_IMAGE016.GIF" wi="24" he="29" />是对应<img file="DEST_PATH_IMAGE029.GIF" wi="26" he="29" />求得的<img file="672111DEST_PATH_IMAGE017.GIF" wi="22" he="29" />的白化矩阵,<img file="657385DEST_PATH_IMAGE024.GIF" wi="23" he="28" />中与最大的特征值对应的特征向量定义为<img file="154094DEST_PATH_IMAGE030.GIF" wi="31" he="33" />;在相应的模式下的投影即为:<img file="DEST_PATH_IMAGE031.GIF" wi="90" he="30" />同理可获得<img file="162501DEST_PATH_IMAGE002.GIF" wi="21" he="22" />、<img file="215908DEST_PATH_IMAGE003.GIF" wi="21" he="22" />的空间滤波器分别为:<img file="233849DEST_PATH_IMAGE032.GIF" wi="161" he="34" />这里<img file="DEST_PATH_IMAGE033.GIF" wi="26" he="29" />是对应<img file="652192DEST_PATH_IMAGE034.GIF" wi="28" he="29" />求得的<img file="DEST_PATH_IMAGE035.GIF" wi="24" he="29" />的白化矩阵,<img file="131584DEST_PATH_IMAGE036.GIF" wi="25" he="28" />中与最大的特征值对应的特征向量定义为<img file="DEST_PATH_IMAGE037.GIF" wi="33" he="33" />;<img file="926365DEST_PATH_IMAGE038.GIF" wi="155" he="34" />这里<img file="DEST_PATH_IMAGE039.GIF" wi="25" he="28" />是对应<img file="807602DEST_PATH_IMAGE040.GIF" wi="28" he="29" />求得的<img file="DEST_PATH_IMAGE041.GIF" wi="25" he="31" />的白化矩阵,<img file="396846DEST_PATH_IMAGE042.GIF" wi="25" he="28" />中与最大的特征值对应的特征向量定义为<img file="DEST_PATH_IMAGE043.GIF" wi="31" he="33" />;以上得到的空间滤波器是64×64维的,每一行向量被称为空间滤波子,对应的空间滤波器与对应的信号相乘即得到滤波后的信号,经过共空间模式滤波后的情绪脑电信号,去除了背景信号的干扰;2)为了更好地进行后续的分类识别,根据滤波后的信号<img file="629113DEST_PATH_IMAGE044.GIF" wi="31" he="32" />做如下变化后作为脑电特征值,<img file="DEST_PATH_IMAGE045.GIF" wi="301" he="72" />(4)特征识别特征提取后,使用支持向量机分类器对特征进行识别,将不同情绪状态区分开。
地址 300072 天津市南开区卫津路92号