技术领域
本发明涉及脑机接口(brain-computer interface,BCI)系统想象动作脑电 信号特征的提取方法,特别涉及利用离散小波变换和BP神经网络提取想象动 作脑电特征的方法。
背景技术
脑-机接口(brain-computer interface,BCI)是指一种不依赖于脑的正常输出 通路(即外周神经和肌肉)的脑-机(计算机或其它装置)通讯系统,它是一种全新 的通讯和控制方法。BCI系统通常由四个部分组成:即信号采集系统、信号 处理系统、
模式识别系统和控制装置系统,如图1所示。信号分析和处理环 节是每个BCI系统的核心部分,其功能是将输入的脑
电信号转换成控制外界 装置的
输出信号。
想象运动是在动作未发生的情况下,对将要发生的动作的一种预测。当 人们想象单侧肢体运动时,大脑对侧的运动感觉区的mu节律和beta节律
能量 减小,而同侧的运动感觉区mu节律和beta节律能量增大,这是大脑神经元突 触后电位相互削弱和增强的结果。这种现象也被称为事件相关去同步和事件 相关同步。也就是说,这两种现象反映了特定
频率段上的脑电信号,在能量 上的减小或增加变化。它成为判断左右手想象运动最根本的依据。但是利用
现有技术,如
叠加平均法、快速
傅立叶变换法,
自回归模型谱估计,独立分 量分析等方法,从含有大量噪声的脑电信号中提取出特征信号,存在着低信 噪比、识别准确率不高,识别速度慢等缺点。
发明内容
为了改善现有技术的不足,本发明的目的是以想象动作思维引起的能量 变化作为区分左右手想象运动的特征,提供一种基于
离散小波变换和BP神经 网络相结合的方法,对含有大量噪声的原始脑电信号进行特征提取及脑电分 类,以提高
信噪比,分类准确率和通讯速度。
为了实现上述目的,本发明采用了如下技术方案:
通过
电极帽采集脑电信号,受试者进行N次想象动作试验,分别包含想 象左手动作和想象右手动作试验各次。数据由差分电极从国际标准的10~ 20
导联系统的C3,C4通道获得(如图1所示)。采集到的脑电信号经过放大, 模数(A/D)转换等前置处理后传送到计算机中,并以信号
电压幅值形式存储于
存储器中。利用计算机实现小波变换和BP神经网络相结合组成的小波神经网 络对脑电信号进行特征提取和分类。N的取值范围是140~300;该方法依次 包括下列步骤:
1.脑电信号预处理
1)计算想象左手动作的
平均功率按照平均功率公式
计算次想象左手动作从C3、C4通道 获取的脑电信号(以下简称左手C3和左手C4)在0~9s内所有
采样点对应的平 均功率PLC3、PLC4,式中为第j个脑电数据所有次试验的平均功率,N 为实验次数,xf(i,j) 2为第i次第j个脑电数据。对求出的平均功率设置3.5~8s 的
时间窗;处理过程如图3所示。
2)计算想象右手动作的平均功率
同理,按照平均功率公式
,计算次想象右手动作从C3、 C4通道获取的脑电信号(以下简称右手C3和右手C4)在0~9s内所有采样点 对应的平均功率PRC3、PRC4。对求出的平均功率设置3.5~8s的时间窗。 处理过程如图4所示。
2.离散小波变换
1)左手C3、C4的离散小波变换
对左手C3、C4在
加窗段内的平均功率进行离散二进小波变换。选用 Daubechies类db5小波,采用6层分解。分解结果如图5所示,可以看出左手C3 在尺度6上的逼近信号L3a6高于左手C4在尺度6上的逼近信号L4a6。
2)右手C3、C4的离散小波变换
同理,对右手C3、C4在加窗段内的平均功率进行离散二进小波变换。同 样选用Daubechies类db5小波,采用6层分解。分解结果如图6所示,可以看出 右手C3在尺度6上的逼近信号R3a6低于右手C4在尺度6上的逼近信号R4a6。
经上述分析可知,对于尺度6上的逼近信号a6,左右手存在明显差异。该 差异表现为左手C3在尺度6上的逼近信号L3a6高于左手C4在尺度6上的逼近信 号L4a6,而右手C3在尺度6上的逼近信号R3a6则低于右手C4在尺度6上的逼近 信号R4a6。因此需要按照公式La6=L3a6-L4a6,Ra6=R3a6-R4a6计算出左右 手C3、C4在尺度6上的逼近信号a6的差值La6和Ra6。选取La6和Ra6作为区分 左右手的特征量。
其中,所述的离散二进小波变换,一维快速分解与重建
算法为:
上式中h(n)和g(n)称为共轭镜像
滤波器组,分别对应于低通和
高通滤波器 的单位脉冲响应。式(1)为小波分解公式,信号f(t)可从j尺度到j+1尺度进行 逐步分解。原始信号
,其小波分解为{AJ df,(Djf)1≤j≤J},J为某一整数, Aj df是在尺度2j下分解得到的逼近信号,Djf是尺度2j下分解得到的细节信 号。式(2)为小波重建公式。
3.利用BP神经网络对特征量进行分类
所述的
基于小波变换和BP神经网络的脑电特征提取方法,其特征在于将 上述离散小波变换与BP神经网络相结合组成小波BP神经网络。
本发明具有如下优点:
本文采用离散小波变换和BP神经网络提取想象运动脑电特征的方法。由 于小波变换具有多
分辨率的特点,在低频时,
时间分辨率较低但频率分辨率 较高;在高频时,时间分辨率较高但频率分辨率较低,它的这种变焦距特性, 容易将类别间差距最大的部分突出表示,从而将不同类之间的差异“放大”, 有助于提高想象动作电位的信噪比和识别正确率。另外,小波变换是一种线 性变换,计算速度快,适合于在线分析。
附图说明
图1为BCI系统基本结构示意图
图2为想象左右手动作实验基本构成示意图
图3为左手C3、C4在0~9s内所有采样点对应的平均功率
图4为右手C3、C4在0~9s内所有采样点对应的平均功率
图5为左手C3、C4经离散小波分解后各尺度的细节信号和逼近信号
其中:d1-d6为尺度1-6上的细节信号,a6为尺度6上的逼近信号;其中 L3a6为左手C3在尺度6上的逼近信号,L4a6为左手C4在尺度6上的逼近 信号;
图6为想象右手C3、C4经离散小波分解后各尺度的细节信号和逼近信 号,其中:d1-d6为尺度1-6上的细节信号,a6为尺度6上的逼近信号;其中 R3a6为右手C3在尺度6上的逼近信号,R4a6为右手C4在尺度6上的逼近 信号;
图7为想象左右手C3、C4在尺度6上的逼近信号a6的差值
其中:La6为左手C3与C4在尺度6上逼近信号a6的差值,Ra6为右手 C3与C4在尺度6上逼近信号a6的差值。
具体实施方式
下面对本发明作进一步详细说明:
想象运动是在动作未发生的情况下,对将要发生的动作的一种预测。人 体大脑的C3,C4
位置包含了想象对侧手运动时最为丰富的信息,即手的运动感 觉区。本发明研究的是想象运动的脑电分类,并针对运动感觉敏感的C3,C4 通道的脑电信号,采用离散小波变换和BP神经网络方法对原始信号进行特征 提取、分类,以达到区分左右手想象运动的目的。
1.实验设计方案
实验是由一个含有反馈信号控制的在线
脑机接口系统完成,所执行的意 识任务是通过想象左右手运动来控制反馈的
光标。因此,在该实验中信号处 理的任务就是在一定的时间内识别出该意识任务是想象左手运动还是右手运 动,实现对测试数据的在线分类和分析。试验次数为140次,分别包含想象左 手动作和想象右手动作试验各70次。实验的电极位置见图1,实验
进程示意见 图2。
所有的N次实验都是在同一天执行的,每两次实验中间均有几分钟的间 隔,实验时间为9s/次。在起初的2s,受试者保持休息状态,t=2s时,显示器上 出现1个持续1s的十字光标,同时会伴随1个
声音信号提示实验开始(受试者准 备开始想象任务)。在t=3~9s时,有1个指示左右方向的箭头代替十字光标作 为指令,同时,要求受试者按照指令所提示的方向通过想象左右手运动控制
进度条向箭头指示的方向移动。通过位于C3,C4前后(前“+”后“-”)各2.5cm 位置的2
对电极记录脑电信号。脑电信号的
采样频率是128Hz。
2.脑电信号的预处理
1)计算想象左手动作的平均功率
按照平均功率公式
,计算70次左手C3、C4在0~9s内所 有采样点对应的平均功率PLC3、PLC4,式中为第j个脑电数据所有次试 验的平均功率,N为实验次数,xf(i,j) 2为第i次第j个脑电数据。对求出的平均 功率设置3.5~8s的时间窗。处理过程如图3所示。加窗后左手C3、C4对应 的平均功率分别为PLC3′、PLC4′。
2)计算想象右手动作的平均功率
同理,按照平均功率公式
,计算70次右手C3、C4在0~9s内 所有采样点对应的平均功率PRC3、PRC4。对求出的平均功率设置3.5~8s的 时间窗。处理过程如图4所示。加窗后右手C3、C4对应的平均功率分别为PRC3 ′、PRC4′。
3.离散小波变换
离散二进小波变换的一维快速分解与重建算法为:
上式中h(n)和g(n)称为共轭镜像
滤波器组,分别对应于低通和高通滤波器 的单位脉冲响应。式(1)为小波分解公式,信号f(t)可从j尺度到j+1尺度进行逐步 分解。原始信号
,其小波分解为{AJ df,(Djf)1≤j≤J},J为某一整数,aj df是 在尺度2j下分解得到的逼近信号,Djf是尺度2j下分解得到的细节信号。式(2) 为小波重建公式。
1)左手C3、C4的离散小波变换
对左手C3、C4在加窗段内的平均功率PLC3’、PLC4’分别进行离散二进小 波变换。选用Daubechies类db5小波,采用6层分解。分解结果如图5所示,可 以看出左手C3在尺度6上的逼近信号L3a6高于左手C4在尺度6上的逼近信号 L4a6。
2)右手C3、C4的离散小波变换
同理,对右手C3、C4在加窗段内的平均功率PRC3’、PRC4’分别进行离散 二进小波变换。同样选用Daubechies类db5小波,采用6层分解。分解结果如图 6所示,可以看出右手C3在尺度6上的逼近信号R3a6低于右手C4在尺度6上的 逼近信号R4a6。
对于小波类型和分解层数的选取,曾先后尝试过Daubechies类db4、db5、 db6三种小波分解类型以及4、5、6层分解。考虑小波基函数与待分析信号波 形的相似性,最终选用Daubechies类db5小波,采用6层分解。并发现了左右手 C3、C4在尺度6上的逼近信号a6存在明显差异。该差异表现为左手C3在尺度6 上的逼近信号L3a6高于左手C4在尺度6上的逼近信号L4a6,而右手C3在尺度6 上的逼近信号R3a6则低于右手C4在尺度6上的逼近信号R4a6。因此需要按照公 式La6=L3a6-L4a6,Ra6=R3a6-R4a6计算出左右手C3、C4在尺度6上的逼近 信号a6的差值La6和Ra6。选取La6和Ra6作为区分左右手的特征量。
4.利用BP神经网络对特征量进行分类
选用BP神经网络作为分类器。BP模型是一个多层
感知机构,是由
输入层、
中间层(隐层)和
输出层构成的前馈网络。相对于其它神经网络模型来说, BP神经网络具有自适应功能、泛化功能以及很强的容错能
力。BP网络的过程 由前向过程和误差反向传播过程组成,其中
输入信号经输入层和隐层神经元 逐层处理,前向传输到输出层输出结果;若输出层的输出值与样本值有误差, 则该误差沿原来的连接通道反向传播,经
修改各层神经元连接
阈值与权值, 缩小误差并反复
迭代,当误差小于允许值时,网络训练过程结束。
其前向计算过程如下:
(1)输入层结点i的输出Oi等于其输入Xi;
(2)隐层
节点j的输入
,输出
,式 中ωji为隐层节点j与输入层结点i之间的连接权;θj为隐层节点j的
阀值,f为非 线形Sigmoid传递函数。
(3)输出层结点l的输入
,输出
式中ωlj为输出层节点l与隐含层结点j之间的连接权;θl为隐层节点1的阀值。
对给定的训练样本集(xp1,xp2,...,xpn)→(lp1,lp2,..lpn),p=1,2,3,...,P为样本号, 网络运算结果与训练样本目标之间的均方误差和表示为:
对于输出层与隐层之间的权值ωli有:
ωlj(k+1)=ωlj(k)+η·δj·oi
δl=f(netl)·(tl-ol)
对于输入层与隐层之间的权值ωji有:
ωji(k+1)=ωji(k)+η·δj·oi
式中k为迭代次数,η为学习率;0<η<1。
网络训练学习的过程就是通过调节网络内部连接权使网络误差最小。BP 网络内部连接权的调整过程也就是误差的反向传播过程。
将上述离散小波变换与BP神经网络相结合组成小波BP神经网络提取想象 运动脑电特征的方法,实际实施效果体现在可以提高信噪比,处理后
诱发电位 波形明显,容易识别。同时缩短了提取视觉诱发电位所需的时间,有助于提高 脑机接口的通讯速度。与传统方法相比,该项发明具有准确率高,识别速度 快的优点。