近似熵结合HHT在脑电信号处理中的应用.doc

Hilbert-Huang变换结合近似熵在脑电信号处理中的应用彭军强1* 吴平东1 殷罡1 （1北京理工大学机械与车辆工程学院，100081）摘 要：Hilbert-Huang Transform（HHT）方法和近似熵方法都是信号处理中的非线性处理方法，近年来在非线性的脑电信号处理中有了不少研究和应用，取得了一些研究成果。本文应用两种方法的结合提出了一种新的脑电信号处理方法，首先用HHT方法把脑电信号分解为多个内在的模式分量，然后求取其中能代表脑电非线性特性的2个分量的近似熵值，用其比值作为脑电的一个特征量。通过对正常静坐、正常驾驶、疲劳静坐、疲劳驾驶4种脑电的实际分析处理，证明本文提出的方法可以更好的表现出脑电的非线性特征。关键词：Hilbert Huang Transform（HHT）；近似熵；脑电中图分类号 R318 文献标识码 AApplication of Hilbert-Huang Transform combining Approximate Entropy in the EEG signal processingPENG Jun-qiang1, WU Ping-dong1, YIN Gang1（1School of Mechanical and Vehicular Engineering, Beijng Institute of Technology, Beijing 10081）Abstract：Hilbert-Huang Transform (HHT) and Approximate Entropy methods both are nonlinear methods in the signal processing. In recent years, the two methods have been applied to the EEG signal processing and some research work has given some useful resolutions. Combining the two methods, the paper proposes a new method in the EEG signal processing, first, decomposing the EEG signal into several IMFs(Intrinsic Mode Function), second, calculating the APEN(Approximate Entropy) values of the two weights which can represent the EEGs nonlinear characters, and the ratio of the two APEN values can be an eigenvector of the EEG. Applying this method to process four kinds of EEG signals, the results indicate that the method can abstract the nonlinear characters of the EEG signal more accurately.Key words：Hilbert-Huang Transform（HHT）；Approximate Entropy（APEN）；EEG0 引 言脑电信号作为人脑思维活动的一种外在表现，在对人脑的工作方式，思维和意识的产生等方面的研究中起着重要的作用。近年来针对脑电信号的应用研究也越来越多，例如脑机接口（BCI），疲劳驾驶检测等。在这些应用中，最关键的步骤就是对脑电信号的分析处理。通过信号处理方法，找出脑电信号中所包含的信息特征，然后进行应用。但是由于大脑的高度复杂性和非线性，脑电信号的微弱性，易受干扰性，从脑电信号中解读大脑的思维内容现在是不可能的，现在所进行的研究工作只是用各种信号处理方法找到脑电信号中能反映某种思维状态的信号特征，然后加以应用。传统的信号处理方法是各种时域方法和频域方法，还有小波分析等时频结合的方法。近年来随着非线性科学的发展，非线性方法开始在脑电信号中得到了应用，其中Hilbert Huang Transform（HHT）方法和近似熵方法是两种比较有潜力的非线性信号处理方法。HHT方法是一种非线性的信号变换方法，把非线性信号变换分解为多个单一模式信号。它同小波变换的最主要区别就是HHT变换是根据信号本身自适应的选取变换基底，而小波分析一旦固定小波基后就不能改变。因此HHT变换更能从本质上对脑电信号进行分解。而近似熵方法是用来评价非线性信号的复杂程度。本文结合这2种方法，对用HHT方法分解得到的特定模式的脑电信号进行近似熵计算来得到脑电信号的内在非线性特征。1 HHT和近似熵方法简介HHT(Hilbert-Huang Transform)技术1998年由NASA的Norden E Huang等提出的新的信号处理方法。该方法适用于非线性非平稳的信号分析, 被认为是近年来对以傅立叶变换为基础的线性和稳态谱分析的一个重大突破。目前HHT 技术已用于地球物理学和生物医学等领域的研究, 并取得了较好的结果。HHT 方法包含两个主要步骤:(1) 对原始数据进行分解。即先通过Experical Mode Decomposition (EMD) 方法, 把数据分解为满足Hilbert 变换要求的n 阶内在模式函数( IMFIntrinsic Mode Function) 。(2) 对分解出的每一阶IMF 进行Hilbert 变换, 得出各自的瞬时频率, 作出时频分布图。其中, EMD 处理后的数据, 即每一阶的IMF 应满足:(1) 数据的极值点和过零点交替出现, 且数目相等或最多相差一个。(2) 在任何点上, 由局部最大值和局部最小值定义的包络的均值必须是零。本文主要用HHT进行脑电信号的分解。近似熵( Approximate Entropy, APEN ) 是由Pincus于1991 年提出的一种度量序列复杂性和统计量化的非线性动力学参数。它用一个非负数来表示一个时间序列的复杂性, 反映时间序列中新信息发生率, 越复杂的时间序列对应的近似熵越大。研究表明, 近似熵能够表征人的生理状态的变化情况。通过实验验证了近似熵在应用中所表现的很好的性质: (1) 所需数据点数少, 可以达到实时观测的目的; (2) 抗噪声和抗干扰能力。特别是对产生的瞬态强干扰有较好的承受能力; (3) 适用于确定性和随机性的信号。最后这一特点对分析生物系统是很有利的, 因为生物信号往往既有确定性成分,也有随机成分。本研究的脑电信号就属于这样的信号, 因此选择近似熵方法作为一种分析工具。HHT和近似熵的具体算法参见相关资料。2 分析用脑电信号的获取实验用的脑电信号来源于正常驾驶和疲劳驾驶实验。正常驾驶实验：实验者在精力充沛时，在模拟驾驶仪上进行汽车驾驶的操作。疲劳驾驶实验：夜间，实验者在汽车模拟驾驶仪上连续开车，同时用日本光电的脑电仪记录实验者脑电。实验者为实验室研究生，身体健康。采用国际统一的10-20系统标准安放电极，参考电极为双耳电极A1和A2，记录位置为F3、F4、C3、C4、P3、P4、O1、O2、T3、T4。电极阻抗均保证达到R正常静坐(S1)疲劳驾驶(S4)疲劳静坐(S3)；而d4分量的近似熵值和d2分量正好相反，从大到小为：疲劳静坐(S3) 疲劳驾驶(S4) 正常静坐(S1) 正常驾驶(S2)；因此，取d2/d4比值作为一个脑电的特征量，可以更清晰的表现出四种信号的区别：正常驾驶(S2)正常静坐(S1)疲劳驾驶(S4)疲劳静坐(S3)。同时由于这种方法中只是使用到了HHT分解中的d2和d4分量，而一般情况下d1分量代表脑电的高频干扰，而d5分量及以后代表了脑电的低频干扰，如身体移动引起的干扰和眼电心电干扰，所以这种方法不必对原始脑电进行任何处理，这也是这种方法的一个优点。根据上面的分析，可以得到如下结果：疲劳时，脑电HHT分解的d2分量，也就是2540Hz左右的脑电的复杂度会减小，而d4分量，也就是515Hz左右的脑电的复杂度会增加；驾驶时，与其相反，脑电HHT分解的d2分量的复杂度会增加，而d4分量的复杂度会增加。5 结论和展望本文通过对正常静坐、正常驾驶、疲劳静坐、疲劳驾驶4种脑电信号用HHT变换结合近似熵计算的方法进行分析，得到了可以较好的分辨4种脑电的特征，说明这种方法可以较好的反映出脑电内在的非线性特征，但同时这种方法的计算量比较大，因此如何改进HHT变换和近似熵的算法减小计算量，是需要进一步研究的内容。参 考 文 献1 Norden E. Huang, Zheng Shen, Steven R. Long. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. The Royal Society, 1998, 454, 903-9952 赵治栋 唐向宏. 基于Hilbert-Huang Transform 的心音信号谱分析. 传感技术学报 2005.3 18(1)Spectrum Analysis of Heart Sound Signal Based on Hilbert-Huang Transform Chinese Journal of Sensors and Actuators3 江朝晖 高翠云 冯焕清. 脑电近似熵分析的思维分类识别. 应用科学学报2006.7 24(4)Mental Task Classif ication Based on Approximate Entropy of EEG. Journal of Applied Sciences4 徐进 郑崇勋 和卫星. 基于脑电近似熵分析的麻醉深度监测研究. 航天医学与医学工程. 2004.6 17(3)Monitoring of Anesthesia Depth Based on Approximate Ent ropy of EEG1XU J in , ZHENG Chong2xun , HE