脑电信号中去除眼电成分..docVIP

脑电信号中眨眼眼电成分的提取 摘要：眨眼伪迹是脑电信号采集过程中的常见噪声,严重影响其有用信息的提取。该文尝试采用独立成分分析中的快速算法分离脑电信号中的各个独立分量,并通过相关性分析自动识别独立分量中的眨眼伪迹干扰并去除。研究结果表明该方法能有效识别和去除眨眼伪迹,在脑电信号的处理中有重要应用价值,值得深入研究和推广。 目的利用独立分量分析方法(ICA) 将混合在观测信号中相互独立的源信号分出来。方法记录3个正常人自然眨眼和水平扫视条件下7道脑电信号和2道眼电信号，选取7道脑电信号进行处理，2道眼电信号用来指示干扰源的情况。使用扩展相似对角化算法( JADE) 将脑电信号分解成多个独立分量，同时利用伪迹脑地形图特征，判断出与眼电伪迹相关分量并将其去除。结果存在于前额电极的眼电干扰被消除，同时其他电极上的信号细节成分较好地保留下来。独立分量分析方法成功去除了脑电信号中的眼电伪迹。 本文针对脑电信号的眼电伪迹去除的问题，运用ICA（独立分量分析）和小波去噪两种方法实现了眼电伪迹去除，并比较分析了两种方法各自的优点和缺点。 关键词：脑电信号 眨眼眼电 ICA 小波去噪 1 引言 脑电( electroencephalogram ，EEG) 信号是一种微弱( μV 级) 的电生理信号，同时具有很强的随机性，极易受其他电生理信号干扰。其中，眼电伪迹是一种最主要的干扰成分。它产生于人体自身，当眨眼( blink) 或是眼球运动( eye movement)时，会在测量电极处引起较大的电位变化形成眼电( electro-oculogram，EOG) 。在采集EEG 时，EOG 从其源发出，弥散到整个头皮，导致采集到的EEG 信号产生明显畸变，形成伪迹，其幅度可达到100 mV。为减少EOG 伪迹影响，要求受试者长时间控制自己的眼部运动。但这通常会引起眼部不适，尤其是部分特定人群( 如患有多动症的儿童、精神分裂症患者等) 的无意运动难以控制。采集到的EEG 信号中会包含EOG 伪迹。 本文采用了ICA算法和小波去噪算法这两种方法来进行脑电信号中眨眼眼电成分的提取，下文将分别介绍这两种算法的原理，实现思想、流程图，以及实验结果。并通过实验结果的分析来比较讨论两种算法的优缺点。 2 ICA 2.1 ICA基础理论 独立分量分析(Independence Component Analyse，ICA)是由盲源信号分离技术发展来的多道信号处理方法【5】。其基本含义是将多道观察信号按照统计独立的原则通过优化算法分解为若干独立成分。从而实现信号的增强和分析。因为采集到的EEG信号是自发脑电信号与各种伪迹的线形混合，满足信号源独立的条件。从而适用ICA方法来除去伪迹。ICA方法的判别依据根据度量各分量之间独立程度的判据不同，有多种形式．如：互信息极小、负熵最大、极大似然估计等。 ICA可描述如下：是M维观测信号，是产生观测信号的N个相互统计独立的源信号，且观测信号X(t)是源信号S(t)经过未知矩阵A线性混合而产生的，即x(t)=As(t)。ICA的目的就是：在混合矩阵A和源信号s(t)未知的情况下，仅利用观测信号x(t)和源是统计独立这一假设，尽可能真实的分离出源信号，即s(t)=Wx(t)，其中W 是分离矩阵，W =A-1。 2.2 FastICA 算法【6】 ICA算法的研究可分为基于信息论准则的迭代估计方法和基于统计学的代数方法两大类，从原理上来讲，它们都是利用了源信号的独立性和非高斯性。基于信息论的方法研究中，主要有基于最大熵、最小互信息、最大似然估计和负熵最大化等角度的一系列估计算法。如FastICA算法、Infomax算法和最大似然估计算法。本文采用极大化非高斯性的ICA估计方法中的FastICA算法。它对任何类型的数据都适用，以其收敛速度快、分离效果好被广泛应用于信号处理领域。FastlCA算法可分为两部分：信号预处理和独立分量提取。 (1)信号预处理包括去均值和白化。从观测信号中减去信号的均值向量，使得观测信号成为零均值变量。该预处理只是为了简化ICA算法。对观测信号白化处理可去除各观测信号之间的相关性，从而简化后续独立分量的提取过程。 (2)独立分量的提取即是分离矩阵不断优化的过程。该过程是一个迭代逼近的过程，通过重复的迭代来寻找W；X(t)的非高斯性最大值。衡量非高斯性的公式为： i=1,2,...,n 式中，K是正常数；G()是非二次函数，E()表示数学期望；为具有和s相同方差的高斯变量。寻找Wi ，使得N(si)达到极大，就可以获得一个独立分量si(t)。分离过程中，对Wi进行调整，当相邻两次的Wi无变化或变化很小时，对一个独立分量的