兰些量型塾墨竺型！！-计算物理.PDF

篓丝 竺竺!!兰些量型塾墨竺型!! 鳖』鳖 [文章编号】 100l一246x(2004)04—034l-04 相空间中脑电近似熵和信息熵的计算 游荣义1’2， 陈 忠2 (1．集美大学计算科学与应用物理系，福建厦门361021；2．厦门大学物理系，福建厦门361005) [摘 要] 提出一种基于相空间重构脑电信号来计算脑电近似熵和信息熵的新方法．实验计算结果表明，癫痫 患者脑电和正常人脑电的近似熵和信息熵随相空间嵌入维数的变化有明显的不同． [关键词] 脑电；相空间；近似熵；信息熵 [中图分类号] 04l；R318．04 [文献标识码】 A O 引言 序列，它是大量脑神经元电活动的总体表现．脑电通常用头皮电极或头盖内电极采集得到，临床上作为脑病 诊断的依据．由于脑电形态的复杂性和多样性，加上脑电信号常常混有各种伪迹或噪声，因此准确地分析有 一定的困难．长期以来，人们在探索脑电分析方法，特别是自动诊断方面做了大量的工作，相继也出现了大量 的分析方法心’31，但可用于临床上的成果依然有限．脑电的非线性被普遍承认以后，非线性动力学分析方法近 年已成为脑电分析的热门方法．其中，计算脑电的关联维数、Lyapunov指数以及近似熵和信息熵等是常见的 一些方法．本文提出一种基于相空间重构脑电信号来计算脑电近似熵和信息熵的新方法，实验计算结果显 示，癫痫患者脑电和健康者脑电的近似熵和信息熵随相空间嵌入维数的变化表现出完全不同的变化规律，为 脑病的I临床诊断分析提供一个新的参考依据． 1计算脑电的近似熵 近似熵(Approximate 示时间序列的复杂性的大小．我们曾用文[4]介绍的方法计算了100例癫痫患者脑电和50例健康者脑电的 近似熵．数据长度从500点到1000点不等．计算结果显示，近似熵的大小和样本数据长度关系不大．癫痫患 间变化，平均值为0．90，二者差别不是很大．可见，从近似熵的大小上看，癫痫患者的脑电和健康者的脑电具 有相当的复杂程度，从而难以区分．在上述方法的基础上，本文提出一种基于相空间重构脑电信号来计算脑 电近似熵的新方法． 1．1相空间重构脑电信号 延时扩展后获得一组由n个状态向量构成的相空间，代表相空间中n个坐标点，形成一条动力系统的运动 轨迹．第忌个向量为 x^=[z(后)，z(后+r)，并(矗+2r)，…，石(尼+(m一1)r)]7，七=l，2，…，n，n：￡一(m一1)r． (1) 1．2近似熵计算方法 近似熵的计算步骤如下： [收稿日期]2003—06—02；[修回日期]2003一11一18 [基金项目]福建省自然科学基金计划(批准号：c0310028)及福建省教育厅重点项目(闽教2001高188)资助项目 [作者简介]游荣义(1957一)，男，福建，副教授，博士，从事生物信号处理、神经网络应用等方面的研究． 342 计 算 物 理 第21卷 ①给定m和r(m取2，3，4，…，M=15分别计算，r取1或2)，按式(1)进行重构． ②计算任意向量x。与其余向量之间的相对欧氏距离如=DH／m(_『=1，2，…，_『≠I|})． 目并求出该数目与向量总数n的比，记为c7． ④将c?取对数，再对所有的七求平均值，记为P“，即 R Pm：土ylnc7． (2) n茸 ⑤m增加l，重复①，④的过程，求得P“，Pm“，P～2，…，∥． 维到高维的熵变化． 本计算方法之所以采用步骤⑤和⑥的过程，其目的在于观察脑电信号的近似熵是如何随相空间嵌入维 数的变化而变化的． 2计算脑电的信息熵 为了进一步了解脑电信源的信息量特征，我们同样采用基于相空间重构脑电信号来计算脑电信息熵 (Info瑚ation Entropy，记为InEn)的方法．其目的是考察脑电信息熵随相空间嵌入维数的变化规律