专利名称：基于小波变换和bp神经网络的脑电特征提取方法
技术领域：
本发明涉及脑机接口 (brain-computer interface, BCI)系统想象动作脑电 信号特征的提取方法，特别涉及利用离散小波变换和BP神经网络提取想象动 作脑电特征的方法。
背景技术：
脑-机接口(brain-computer interface, BCI)是指一种不依赖于脑的正常输出 通路(即外周神经和肌肉)的脑-机(计算机或其它装置)通讯系统，它是一种全新 的通讯和控制方法。BCI系统通常由四个部分组成即信号采集系统、信号 处理系统、模式识别系统和控制装置系统，如图1所示。信号分析和处理环 节是每个BCI系统的核心部分，其功能是将输入的脑电信号转换成控制外界 装置的输出信号。
想象运动是在动作未发生的情况下，对将要发生的动作的一种预测。当 人们想象单侧肢体运动时，大脑对侧的运动感觉区的mu节律和beta节律能量 减小，而同侧的运动感觉区mu节律和beta节律能量增大，这是大脑神经元突 触后电位相互削弱和增强的结果。这种现象也被称为事件相关去同步和事件 相关同步。也就是说，这两种现象反映了特定频率段上的脑电信号，在能量 上的减小或增加变化。它成为判断左右手想象运动最根本的依据。但是利用 现有技术，如叠加平均法、快速傅立叶变换法，自回归模型谱估计，独立分 量分析等方法，从含有大量噪声的脑电信号中提取出特征信号，存在着低信噪比、识别准确率不高，识别速度慢等缺点。
发明内容
为了改善现有技术的不足，本发明的目的是以想象动作思维引起的能量
变化作为区分左右手想象运动的特征，提供一种基于离散小波变换和BP神经 网络相结合的方法，对含有大量噪声的原始脑电信号进行特征提取及脑电分 类，以提高信噪比，分类准确率和通讯速度。
为了实现上述目的，本发明采用了如下技术方案-
通过电极帽釆集脑电信号，受试者进行N次想象动作试验，分别包含想
象左手动作和想象右手动作试验各^iV次。数据由差分电极从国际标准的10
2
20导联系统的C3,C4通道获得(如图1所示)。采集到的脑电信号经过放大， 模数(A/D)转换等前置处理后传送到计算机中，并以信号电压幅值形式存储于 存储器中。利用计算机实现小波变换和BP神经网络相结合组成的小波神经网 络对脑电信号进行特征提取和分类。N的取值范围是140 300;该方法依次 包括下列步骤
1.脑电信号预处理
l)计算想象左手动作的平均功率
按照平均功率公式i^ 计算^iV次想象左手动作从C3、 C4通道
获取的脑电信号(以下简称左手C3和左手C4)在0 9s内所有采样点对应的平
均功率PLC3、 PLC4，式中/[力为第j个脑电数据所有次试验的平均功率，N 为实验次数，《,.,，为第i次第j个脑电数据。对求出的平均功率设置3. 5 8s的时间窗；处理过程如图3所示。
2)计算想象右手动作的平均功率
同理，按照平均功率公式^)=去|>/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作为区分左右手的特征量。 '其中:，所述的离散二进小波变换， 一维快速分解与重建算法为<formula>formula see original document page 7</formula> (2)上式中h(n)和g(n)称为共轭镜像滤波器组，分别对应于低通和高通滤波器 的单位脉冲响应。式(l)为小波分解公式，信号f(t:)可从j尺度到j+l尺度进行 逐步分解。原始信号/(/)=《/,其小波分解为"7，(iVU,J为某一整数，是在尺度下分解得到的逼近信号，是尺度2卩下分解得到的细节信 号。式(2)为小波重建公式。3.利用BP神经网络对特征量进行分类所述的基于小波变换和BP神经网络的脑电特征提取方法，其特征在于将 上述离散小波变换与BP神经网络相结合组成小波BP神经网络。 本发明具有如下优点本文采用离散小波变换和BP神经网络提取想象运动脑电特征的方法。由 于小波变换具有多分辨率的特点，在低频时，时间分辨率较低但频率分辨率 较高；在高频时，时间分辨率较高但频率分辨率较低，它的这种变焦距特性， 容易将类别间差距最大的部分突出表示，从而将不同类之间的差异"放大"， 有助于提高想象动作电位的信噪比和识别正确率。另外，小波变换是一种线 性变换，计算速度快，适合于在线分析。
图1为BCI系统基本结构示意2为想象左右手动作实验基本构成示意3为左手C3、 C4在0 9s内所有釆样点对应的平均功率图4为右手C3、 C4在0 9s内所有采样点对应的平均功率图5为左手C3、 C4经离散小波分解后各尺度的细节信号和逼近信号其中dl-d6为尺度l-6上的细节信号，a6为尺度6上的逼近信号；其中L3a6为左手C3在尺度6上的逼近信号，L4a6为左手C4在尺度6上的逼近信号；图6为想象右手C3、 C4经离散小波分解后各尺度的细节信号和逼近信 号，其中dl-d6为尺度l-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次。实验的电极位置见图l，实验进程示意见 图2。所有的N次实验都是在同一天执行的，每两次实验中间均有几分钟的间 隔，实验时间为9s/次。在起初的2s，受试者保持休息状态，t-2s时，显示器上 出现l个持续ls的十字光标，同时会伴随l个声音信号提示实验开始(受试者准 备开始想象任务)。在t-3 9s时，有l个指示左右方向的箭头代替十字光标作 为指令，同时，要求受试者按照指令所提示的方向通过想象左右手运动控制 进度条向箭头指示的方向移动。通过位于C3， C4前后(前"+"后"-")各2. 5cm 位置的2对电极记录脑电信号。脑电信号的采样频率是128Hz。2. 脑电信号的预处理 l)计算想象左手动作的平均功率按照平均功率公式<formula>formula see original document page 9</formula>
计算70次左手C3、 C4在0 9s内所有采样点对应的平均功率PLC3、 PLC4，式中&)为第j个脑电数据所有次试 验的平均功率，N为实验次数，x;力为第i次第j个脑电数据。对求出的平均功率设置3. 5 8s的时间窗。处理过程如图3所示。加窗后左手C3、 C4对应的平均功率分别为PLC3' 、 PLC4'。2)计算想象右手动作的平均功率同理，按照平均功率公式4)=丄|>〗(,、,.)，计算70次右手C3、 C4在0 9s内所有采样点对应的平均功率PRC3、 PRC4。对求出的平均功率设置3.5 8s的 时间窗。处理过程如图4所示。加窗后右手C3、 C4对应的平均功率分别为PRC3 '、PRC4'。 3.离散小波变换离散二进小波变换的一维快速分解与重建算法为- * (1)</ = Z - 2《/ + !>("- (2) 上式中h(n)和g(n)称为共轭镜像滤波器组，分别对应于低通和高通滤波器 的单位脉冲响应。式(l)为小波分解公式,信号柳可从j尺度到j+l尺度进行逐步分解。原始信号/(/)=《/,其小波分解为(^/,(",/U， J为某一整数，^/是在尺度2'下分解得到的逼近信号，",/是尺度2'下分解得到的细节信号。式(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的输出Q等于其输入Z,;(2) 隐层节点j的输入"",Zc^+A，输出0 _~^——，式J 7 ■' '' l + exp(—""》中^,为隐层节点j与输入层结点i之间的连接权；《.为隐层节点j的阀值，伪非 线形Sigmoid传递函数。(3) 输出层结点1的输入一 =^x。,+《，输出0/ = ——^——， l + exp(-"W,)式中^.为输出层节点l与隐含层结点j之间的连接权；《为隐层节点l的阀值。 对给定的训练样本集(、1 ,...^ )^(//^2，...，^)， p-l,2,3,…，P为样本号，网络运算结果与训练样本目标之间的均方误差和表示为-J/2对于输出层与隐层之间的权值％有 fij". (A +1) =(A) + /7.《.A《=X(n)对于输入层与隐层之间的权值％有 . (A:'+1) =(A:) + 77 《 o,.式中k为迭代次数，ll为学习率；0<T1<1。网络训练学习的过程就是通过调节网络内部连接权使网络误差最小。BP网络内部连接权的调整过程也就是误差的反向传播过程。将上述离散小波变换与BP神经网络相结合组成小波BP神经网络提取想象 运动脑电特征的方法，实际实施效果体现在可以提高信噪比，处理后诱发电位 波形明显，容易识别。同时縮短了提取视觉诱发电位所需的时间,有助于提高 脑机接口的通讯速度。与传统方法相比，该项发明具有准确率高，识别速度 快的优点。
权利要求
1、基于小波变换和BP神经网络的脑电特征提取方法，通过电极帽采集脑电信号，受试者进行N次想象动作试验，分别包含想象左手动作和想象右手动作试验各次；采集到的脑电信号经过放大，模数转换前置处理后传送到计算机中，并以信号电压幅值形式存储于存储器中；其特征在于，所述的N的取值范围是140～300；依次包括下列步骤1)脑电信号预处理i)计算想象左手动作的平均功率按照平均功率公式<math-cwu><![CDATA[<math> <mrow><mover> <msub><mi>P</mi><mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo></mrow> </msub> <mo>‾</mo></mover><mo>=</mo><mfrac> <mn>1</mn> <mi>N</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>N</mi></munderover><msubsup> <mi>x</mi> <mrow><mi>f</mi><mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo></mrow> </mrow> <mn>2</mn></msubsup> </mrow></math>]]></math-cwu><!--img id="icf0002" file="S2008100568387C00012.gif" wi="120" he="43" img-content="drawing" img-format="tif"/-->计算次想象左手动作从C3、C4通道获取的脑电信号在0～9s内所有采样点对应的平均功率PLC3、PLC4，式中为第j个脑电数据所有次试验的平均功率，N为实验次数，xf(i，j)2为第i次第j个脑电数据；对PLC3、PLC4设置3.5～8s的时间窗；记加窗后左手C3、C4通道对应的平均功率分别为PLC3′、PLC4′；ii)计算想象右手动作的平均功率按照平均功率公式<math-cwu><![CDATA[<math> <mrow><mover> <msub><mi>P</mi><mrow> <mo>(</mo> <mi>j</mi> <mo>)</mo></mrow> </msub> <mo>‾</mo></mover><mo>=</mo><mfrac> <mn>1</mn> <mi>N</mi></mfrac><munderover> <mi>Σ</mi> <mrow><mi>i</mi><mo>=</mo><mn>1</mn> </mrow> <mi>N</mi></munderover><msubsup> <mi>x</mi> <mrow><mi>f</mi><mrow> <mo>(</mo> <mi>i</mi> <mo>,</mo> <mi>j</mi> <mo>)</mo></mrow> </mrow> <mn>2</mn></msubsup> </mrow></math>]]></math-cwu><!--img id="icf0005" file="S2008100568387C00015.gif" wi="120" he="43" img-content="drawing" img-format="tif"/-->，计算次想象右手动作从C3、C4通道获取的脑电信号在0～9s内所有采样点对应的平均功率PRC3、PRC4；对PRC3、PRC4设置3.5～8s的时间窗；加窗后右手C3、C4对应的平均功率分别为PRC3′、PRC4′；2)离散小波变换i)左手C3、C4的离散小波变换对左手C3、C4在加窗段内的平均功率PLC3′、PLC4′进行离散二进小波变换；选用Daubechies类db5小波，采用六层分解；ii)右手C3、C4的离散小波变换对右手C3、C4在加窗段内的平均功率PRC3′、PRC4′分别进行离散二进小波变换；同样选用Daubechies类db5小波，采用六层分解；按照公式La6＝L3a6-L4a6，Ra6＝R3a6-R4a6计算出左右手C3、C4在尺度6上的逼近信号a6的差值La6和Ra6；其中L3a6、L4a6分别为左手C3、C4在尺度6上的逼近信号，R3a6、R4a6分别为右手C3、C4在尺度6上的逼近信号；选取La6和Ra6作为区分左右手的特征量；3)利用BP神经网络对特征量进行分类按照建立BP神经网络、对特征量进行训练、仿真的步骤来完成这一过程。
全文摘要
本发明是脑机接口系统想象动作脑电信号特征的提取方法，特别是基于小波变换和BP神经网络的脑电特征提取方法。本发明是以想象动作思维引起的能量变化作为区分左右手想象运动的特征，按照平均功率公式分别计算出想象左右手动作从C3、C4通道获取的脑电信号(以下简称左右手C3、C4)在0～9s内所有采样点对应的平均功率。设置时间窗，对加窗段的数据进行离散二进小波变换，选取第六尺度上的逼近信号a6，作为信号特征；以BP神经网络作为分类器进行分类。本文采用离散小波变换和BP神经网络提取想象运动电位的方法有助于提高想象动作电位的信噪比和识别正确率；另外，小波变换是一种线性变换，计算速度快，适合于在线分析。
文档编号G06F17/00GK101221554SQ200810056838
公开日2008年7月16日 申请日期2008年1月25日 优先权日2008年1月25日
发明者刘净瑜, 亮 孙, 左国玉, 李明爱, 蕊 王, 郝冬梅, 阮晓钢 申请人:北京工业大学