1、生产实习课程论文基于典型相关分析和小波变换的眼电伪迹去除研究学 院（系）： 电子信息与电气工程学部专 业： 生物医学工程 学 生 姓 名： 学 号： 指 导 教 师： 邱天爽 完 成 日 期： 2013.07.20 - II -目 录1 绪论.11.1 脑电信号处理的意义11.2 去除眼电伪迹方法的进展11.2.1 早期的人工处理.11.2.2 现代的相关去噪算法.11.3 wCCA 算法的提出 .22 wCCA 算法 22.1 基于典型相关分析的盲源分离方法22.2 小波阈值去噪32.3 基于 wCCA 的盲源分离方法去除眼电伪迹 33 程序说明.43.1 算法流程图43.2 相关 matl
2、ab 函数 .43.3 相关参数描述44 实验结果分析.64.1 仿真结果64.2 数据比较分析7结 论.8基 于 典 型 相 关 分 析 和 小 波 变 换 的 眼 电 伪 迹 去 除 研 究11 绪论1.1 脑电信号处理的意义脑电活动首次于 1924 年被德国精神教授测量并定名为 EEG。EEG 信号作为一种直接反应大脑内部状态的生物电信号，其中蕴含了大量的心理、生理和病理信息。目前被广泛运用于神经心理学、大脑意识及认知、脑部疾病的诊治、脑机接口等诸多研究领域中。与 EEG 信号研究紧密关系的另一种典型技术为脑机接口技术 BCI，BCI 可以完全不依赖于外围神经核肌肉的参与，直接实现人与计
3、算机之间或外部环境之间的通信。BCI 一般可以分为以自发脑电信号的 BCI 系统和使用命令的转换算法 。BCI 系统可以使有运动障碍的人通过 EEG 信号来与外界进行交流。提高 BCI 系统有效性的另一个重要手段就是提高采集到的微弱的 EEG 信号的信噪比。脑电信号一般通过放置于大脑头皮的电极进行采集，但是实际采集到的脑电信号非常微弱，只有微伏极。由于脑电信号是一种易变的非平稳信号，其在采集过程当中，会不可避免地混入非脑神经组织产生的各种伪迹信号，如眼电（眨眼或眼动），它的幅度比脑电信号大好几倍，所以如何对脑电信号进行预处理、去除各种伪迹成分，并从中提取出有效的脑电信号成分，是各国研究者关注的
4、重要问题，具有重大的理论和实践意义。1.2 去除眼电伪迹方法的进展1.2.1 早期的人工处理在早期，意思和研究者通过实验控制来处理无关的电生理伪迹成分，比如让患者和被试者避免或者减少眨眼、眼动、吞咽以及四肢运动等，这样会添加附加的实验任务，并且不易于控制，如当患者或被试者为儿童时，比较难以控制，因而会影响实验效果。一般情况下，EEG 信号伪迹去除的通用方法是去除含有伪迹成分的 EEG 信号【片断。比如，识别眼电伪迹（主要包括眨眼和眼动伪迹），通常通过检测眼电导联记录的电平超过一定的固定阈值，其他的伪迹成分或干扰的检测可以通过人工标记并去除，去除含有伪迹成分的 EEG 片断必然会引起有效的 EE
5、G 信号成分的大量损失。比如识别人物当中的眨眼可能就是识别任务的一种反应，若剔除就可能导致重要信息的丢失，另外，对于一些病人来说，剔除被污染的脑电数据就意味着病情的漏诊。1.2.2 现代的相关去噪算法自适应滤波法。自适应滤波器可以自动调节参数，在设计时无需任何关于信号和噪声的先验统计知识，在信号和噪声的先验知识未知的情况下，可以采用自适应滤波方法来去除噪声。基 于 典 型 相 关 分 析 和 小 波 变 换 的 眼 电 伪 迹 去 除 研 究2小波变换方法是 20 世纪 80 年代中期发展起来的一种时域分析方法。传统的傅里叶分析方法在处理平稳信号方面具有显著优势，经过其变换的信号具有最大的频率
6、分辨率，但是不具备时空定位信息，而小波变换由于其窗口可以根据频率分辨率的高低的而进行自适应调节，从而具有多分辨特性，小波变换在低频部分具有较高的频率分辨率和较低时间分辨率。而在高频部分具有较高的时间分辨率和较低的频率分辨率，这种多尺度特性适合于分析生物医学信号等非平稳信号。盲源分离是信号处理领域一个新的研究热点，它尝试在源信号和传输系统特性均未知的情况下对混合信号进行分离。盲源分离法将伪迹成分和 EEG 信号分解成不同的源信号成分，通过将与伪迹有关的源信号成分置零，可以得到去除伪迹后的信号。盲源分离问题可以采用许多不同的算法以及原则来估计源信号成分，通常采用基于二阶统计量和高阶统计量的方法来实
7、现盲源分离。1.3 wCCA 算法的提出针对脑电信号中眼电伪迹去除尚存在的问题，提出一种基于典型相关分析与小波变换的 (wavelet-enhanced canonical correlation analysis, wCCA)自动去除眼电伪迹的算法。首先，充分利用脑电信号和眼电伪迹的空间分布特征，将基于典型相关分析的盲源分离算法以一种全新的方式应用于混合信号中，从而保证典型相关分析分解得到的第一个典型相关变量（即左右脑区之间的最公共成分），就是与眼电伪迹相关的分量。其次为了恢复泄漏在该伪迹分量中的脑电成分，对伪迹分量进行小波阈值去噪，仅将小波系数高于某一阈值的分量置零。与其他三种基于盲源分离
8、去除眼电伪迹的方法相比较，该方法在有效地自动去除眼电伪迹的同时，很好地保留了潜在的脑电信号。 2 wCCA 算法2.1 基于典型相关分析的盲源分离方法讨论相关关系常用的一种方法是讨论第一组每个变量和第二组中每个变量的相关，得到 pq 个相关系数，用这些相关系数反映两组变量的关系。但这样做是不够的，既繁琐，又抓不住要领。另外一种方法类似于主分量分析，对每组变量做一个线性组合，称其为这组变量的综合变量，然后研究两组综合变量的相关，通过少数几个综合变量来反映两组变量的相关性质，这样可以抓住它们的主要关系，而且又简明。因此典型相关分析揭示了两组变量之间的内在关系，更深刻的反映了这两组随机变量之间的线性
9、相关情况。综合变量对间的相关强弱程度不同，就形成了不同的典型相关对.在实际中，往往只需重点研究相关关系较大的几对典型变量，因为它们反映了两组变量间相互关系的绝大部分信息。这就是典型相关分析的主要思想。基 于 典 型 相 关 分 析 和 小 波 变 换 的 眼 电 伪 迹 去 除 研 究3假设 X 和 Y 是两组观测信号。典型相关分析寻找 X 和 Y 的线性组合,即：使得产生的新变量 u 和 v 之间的相关程度最大。2.2 小波阈值去噪典型相关分析分解出的伪迹分量中还含有少量的高频脑电成分，若将该分量完全置零可能会造成感兴趣脑电成分的损失，故需要进一步恢复泄漏在伪迹分量中的脑电成分。小波阈值去噪
10、是同时基于时空域和频率的去噪方法。经过 CCA 分解得到的伪迹分量中，眼电成分的能量在小波域集中在一些大的小波系数中，而脑电成分的能量却分布于整个小波域内。因此经过小波分解后，伪迹成分的小波系数幅值要大于脑电成分的小波系数幅值。对伪迹分量进行小波阈值去噪，将小波系数高于某一阈值的分量置零，这样可以把大部分脑电成分对应的小波系数保留，而将伪迹成分系数置零。2.3 基于 wCCA 的盲源分离方法去除眼电伪迹由电极采集到的脑电信号是由不同来源的脑电和伪迹叠加而成的。不同的脑电信号在头皮上分布的区域是不同的，不同脑区采集到的脑电信号是不同的。另一方面，由于大脑是电的良导体，眼电信号从前额处后向传播，遍
11、历整个头部。充分利用脑电信号和眼电伪迹空间分布特点的不同，本文将 CCA 算法以一种全新的方式应用于混合信号中，令 X(t)为左侧脑区的脑电信号，Y(t)为右侧脑区的脑电信号。同时将垂直眼电信号分别加入到 X(t)和 Y(t)中，来提高得到的第一个典型相关变量中伪迹成分的百分比。CCA 分解得到的第一对典型相关变量之间的相关性最大，故该分量可以认为是 X(t)和 Y(t)之间最公共成分，即左右脑区之间最公共成分。而该成分是由与伪迹相关及少量高频脑电成分构成的。借助于这种方式的 CCA，巧妙的回避了基于盲源分离方法中面临的伪迹成1.kkTxxyywxywv基 于 典 型 相 关 分 析 和 小
12、波 变 换 的 眼 电 伪 迹 去 除 研 究4分人工识别的问题。当然若将该伪迹分量全部置零再重构的话，会造成其中脑电成分的损失。为此，可利用小波阈值去噪来去除该伪迹分量中眼电成分，保留脑电成分。3 程序说明3.1 算法流程图3.2 相关 matlab 函数本程序主要用的的函数有waverec 小波重构函数 、wavedec小波分解函数、appcoef 低频分量（尺度系数）提取函数、detcoef高频分量（小波系数）提取函数3.3 相关参数描述X1-将FP1 、F3、 C3 、O1 、垂直眼电信号组合作为一路51000Y1 -将FP2、F4、 C4 、O2 、垂直眼电信号组合作为一路51000
13、A、B -分别为X1 与Y1每一行的均值为 51的矩阵X、Y -经中心化处理的X1和Y1 信号WX、 Wy-由CAA算法处理得到的特征向量基 于 典 型 相 关 分 析 和 小 波 变 换 的 眼 电 伪 迹 去 除 研 究5u1、 v1-由典型相关向量得到的典型相关变量C1、 C2、L1、L2- 由小波分解得到的小波解向量和相应的记录长度cA5 c1A5-u v 的尺度系数（由函数appcoef得到的低频分量）cD1、cD2、cD3、cD4、cD5、c1D1、c1D2 、c1D3、c1D4、 c1D5-由detcoef得到的u、v各层小波系数k1k6、 p1p6-各个低频阈值和高频阈值U、V
14、 -经小波重构得到的信号x1、 y1-别对两路信号中的第一行处理后得到的信号基 于 典 型 相 关 分 析 和 小 波 变 换 的 眼 电 伪 迹 去 除 研 究64 实验结果分析4.1 仿真结果FP1 与 FP2 两路信号的处理前后结果对比图0 200 400 600 800 1000-5000500FP1信 信 信 信 信 信 信 信0 200 400 600 800 1000-5000500FP2信 信 信 信 信 信 信 信0 200 400 600 800 1000-5000500FP1信 信 信 信 信 信 信 信wCCA0 200 400 600 800 1000-5000500
15、FP2信 信 信 信 信 信 信 信wCCAF3 与 F4 两路信号的处理前后结果对比图0 500 1000-5000500F3信 信 信 信 信 信 信 信0 500 1000-5000500F4信 信 信 信 信 信 信 信0 500 1000-5000500F3信 信 信 信 信 信 信 信wCCA0 500 1000-5000500F4信 信 信 信 信 信 信 信wCCA基 于 典 型 相 关 分 析 和 小 波 变 换 的 眼 电 伪 迹 去 除 研 究7C3 与 C4 两路信号的处理前后结果对比图0 200 400 600 800 1000-5000500C3信信信信信信信信0
16、200 400 600 800 1000-5000500C4信信信信信信信信0 200 400 600 800 1000-5000500C3信信信信信信信信wCCA0 200 400 600 800 1000-5000500C4信信信信信信信信wCCAO1 与 O2 两路信号的处理前后结果对比图0 200 400 600 800 1000-5000500O1信信信信信信信信 0 200 400 600 800 1000-5000500O2信信信信信信信信0 200 400 600 800 1000-5000500O1信信信信信信信信wCCA0 200 400 600 800 1000-5000
17、500O2信信信信信信信信wCCA4.2 数据比较分析由实验结果截图可以看出，本次的实验程序基本上解决了脑电信号中眼电伪迹的去除工作，不同部位电极处的信噪比不同，可知不同部位的去除效果不一样，这可能和眼电传到不同部位的时间和距离有关，与眼电的空间分布特性和传输有关，FP1 与 FP2 处的处理效果不是很好，主要和前额两处电极受眼电干扰大的原因有关，O1 和 O2 处受眼电影响较小，这点由图可以看出来，可能是因为这两点和眼电产生部位的距离有关，从图中看不出来对这两个部位的处理效果，但是实验处理前后的数据是有变化的，所以说即使不能肉眼观察出来这两处的处理效果，但是实验数据能说明该算法还是对这两点起
18、到了去伪迹的效果。基 于 典 型 相 关 分 析 和 小 波 变 换 的 眼 电 伪 迹 去 除 研 究8反思与讨论在评价眼电伪迹去除效果的时候，需要综合考虑两个方面的目标：（1）考虑眼电伪迹去除的干净程度，即是否将眼电伪迹去除彻底；（2）考虑脑电信号的损失程度，即是否在去除眼电伪迹的同时也破坏了感兴趣的脑电信号。这样才能科学地评价一个伪迹去除算法的好坏。我们在这方面的研究与比较还有所欠缺参考文献典型相关分析算法理论及其在模式分类中的应用：西安电子科技大学 赵峰 2005.01基于典型相关分析和小波变换的眼电伪迹去除：大连理工大学 邱天爽实时脑电信号眼电伪差去除方法的研究：西安交通大学学报 2004.12基于脑电的脑_机接口_关键技术和应用前景：电子科技大学学报 2009.05MATLAB 小波分析-机械工程出版社脑电信号中眼电伪迹去除方法研究：电子科技大学 刘铁军 2008.12
