[发明专利]基于信息增益的脑电波形时频特征提取方法

摘要

本发明提供了脑电波形识别技术领域的一种基于信息增益的脑电波形时频特征提取方法。脑电波形是脑电图中的基本元素，在进行医疗诊断时必须先识别出各类脑电波形，本发明是自动识别脑电波形技术中的特征提取方法。本发明的方法包括小波变换、离散化和信息增益选择三个主要步骤。小波变换把每一个时域样本变换为小波域样本，变换结果按照时间上从前到后频率上从高频到低频的次序排列。离散化把每个小波域样本从连续值向量转变为离散值向量。信息增益选择在离散化的小波域样本基础上依据信息增益从小波域特征中选取对分类作用较大的特征。本发明的优点在于：突出了各类脑电波形的相互差异，提升了脑电波形的识别效率。

主权项

1.基于信息增益的脑电波形时频特征提取方法，其特征在于，包括如下步骤：/n步骤S1：从脑电源信号中获取脑电波形，按设定长度截取各脑电波形得到脑电波形的样本向量，依据所述脑电波形的长度、类别信息以及样本向量构建脑电波形的样本集合；/n步骤S2：利用离散小波变换处理所述样本集合得到小波域样本集合；/n步骤S3：针对小波域样本向量的每个维度执行离散化过程，把每个小波域样本从连续值向量转变为离散值向量；/n步骤S4：对所述离散值向量的小波域样本集合执行信息增益选择过程，输出有利于分类的特征；/n步骤S5：依据所述特征构建特征向量，用于脑电波形的识别；/n所述步骤S3中的离散化过程具体包括：/n步骤S31：取出所有小波域样本向量在该维度上的值，按照脑电波形的类别分成不同的集合；/n步骤S32：计算各集合的均值和方差；/n步骤S33：依据所述均值和方差分别为各集合绘制一条正态分布曲线；/n步骤S34：求出相邻正态分布曲线的交点，利用所述交点把数轴分成多个取值区间，并给每个区间一个标记；/n步骤S35：把小波域样本向量在该维度上的连续值替换为对应的区间标记；/n所述步骤S4中的信息增益选择过程具体包括：/n步骤S41：设置信息增益阈值为G；设置check[i]＝false，i＝1,…,M；用node结构创建变量T，设置T.cur＝0，T.feature[i]＝0，i＝1,…,M，设置T.sample[i]＝true，i＝1,…,N；创建队列queue，将T插入queue；/n步骤S42：若queue为空，check[i]＝true所对应的特征为本过程所选择的特征，结束信息增益选择过程，返回主步骤；否则，进入步骤S43；/n步骤S43：从queue队头取出变量赋值给T；从F中删除T.feature[i]，i＝1,…,T.cur，赋值给subF；/n步骤S44：以T.sample[i]＝true为条件从小波域样本集合中筛选出一个子集D；依据D中的样本类别信息，针对subF中所有特征依据样本类别信息和区间的标记值计算信息增益，其中信息增益最大的为gMax，对应的特征标识为fMax；/n步骤S45：若gMax＞G，则进入步骤S46；否则，check[T.feature[i]]＝true,i＝1,…,T.cur,转入步骤S42；/n步骤S46：根据小波域样本第fMax个特征的不同值个数K，分别针对i＝1,…,K，用node结构创建newT[i],设置newT[i].feature＝T.feature，newT[i].cur＝T.cur+1,newT[i].feature[newT[i].cur]＝fMax，newT[i].sample＝T.sample，若样本的第fMax个特征不为第i个值则把newT[i].sample的相应项改为false，把newT[i]插入队列queue；转入步骤S42；/n其中，M为小波域样本的特征个数，N为小波域样本集合中样本的个数；用整数1,…,M标识特征，F为所有特征标识的集合；node结构包含的成员cur为整数类型，cur用于记录算法在当前位置已选出的特征数量；成员feature[M]为整型数组，它的前cur项用于记录算法在当前位置已选出的特征标识；sample[N]成员为布尔型数组，用于记录算法在当前位置已筛去了哪些样本，sample[i]等于false表示样本i已被筛去；check[M]为布尔型数组，用于记录算法所选择的特征，check[i]等于true表示第i个特征被算法选中。/n