Ⅱ.临床癫痫患者脑电信号近似熵分析92

1、癫痫脑电信号分析研究指导教师: 夏阳 副教授 研 究 生: 刘衍素电子科技大学生命科学与技术学院1本文所要解决的问题将近似熵用于临床癫痫患者脑电信号分析，探讨理论分析方法在临床上的实用价值；用似然同步、近似熵分析痫样放电大鼠皮层和海马脑电信号，探讨癫痫大鼠在不同脑功能状态下、脑功能状态转换时的不同脑区及其之间的功能改变。 2内容提要.研究背景.临床癫痫患者脑电信号分析.癫痫大鼠脑电信号分析实验及数据预处理近似熵分析似然同步分析.工作总结和展望.致谢3.研究背景癫痫Epilepsy是由于脑部神经细胞群异常放电引起的一过性脑功能紊乱综合症，是一种严重且较顽固的慢性疑难性脑部疾病，其发作时表现为大脑

2、神经元群兴奋性增高以及过度同步化放电，并导致短暂性中枢神经系统功能失常 。 研究意义癫痫的发病率高约5 ，危害大 ；癫痫EEG 有其特征，但在临床上从EEG的中获得的信息有限，需要深入挖掘；动物模型可以做到有创分析，可以完成临床不能进行的一些研究 4 癫痫脑电信号中常见的特征波有5癫痫脑电信号的分析方法线性分析时域分析频域分析时频分析非线性分析关联维数Correlation dimension李亚普诺夫指数Lyapunov exponent复杂度Complexity算法复杂度近似熵 同步分析 似然同步Synchronization likelihood匹配追踪算法6近似熵Pincus，1991

3、)： 用一个非负数来表示一个时间序列的复杂性方法，越复杂的时间序列对应的近似熵越大 。算法优点计算所需数据短1005000点抗噪及抗野点能力强特别是对偶尔产生的瞬态强干扰适用于确定性信号及随机信号 本文采用洪波等1999提出的快速算法，其计算速度提高了5倍 （1） （3） （2） （4） （5）7似然同步Stam，2002) ：是一种一个时间序列和另一个或多个时间序列间动态相互作用的测量方法 算法优点既能测量线性又能测量非线性的相互作用 能够反映这种相互作用随时间的变化 ,适用于不稳定的时间序列 用似然同步分析各脑区、各节律间同步性差异及状态转化时的变化未见报道 （1） （2） （3） （4）

4、 （5）8.临床癫痫患者脑电信号近似熵分析数据来源：四川大学华西医学中心神经内科提供的临床癫痫患者的Video-EEG数据及相关临床信息 数据点及相似容限的选取 不同r时ApEn与N的对应关系 无痫样放电的EEG信号参数为： ( 表示标准差) ， 点的滑动窗，重叠240点 不同r时ApEn与N的对应关系 无痫样放电的EEG信号9分析所用数据的选取无痫样放电的EEG片断 相邻且有痫样放电的EEG片断 10分析结果采用上述参数对10例临床癫痫患者计算有无痫样放电脑电信号的相邻600个点3秒的近似熵平均值结果如图 。有无痫样放电EEG的ApEn比较 结果显示:无痫样放电时脑电信号的近似熵显著大于有痫

5、样放电时脑电信号的近似熵 p=0.0020.05 ;2不同脑功能状态间转换时没有痫样放电与连续性痫样放电 连续性痫样放电与周期性痫样放电 LC *ApEnnonApEncontinuous ApEncontinuous=ApEnperiod LH *ApEnnonApEncontinuous ApEncontinuous=ApEnperiod RC *ApEnnonApEncontinuous ApEncontinuous=ApEnperiod RH *ApEnnonApEncontinuous ApEncontinuous=ApEnperiod 注：*P0.01 18由表可见：不同脑功能状态

6、间转换时无痫样放电 连续性痫样放电连续性痫样放电 周期性痫样放电ApEn显著降低p0.05复杂的混沌状态 简单的同步活动1、该两种状态间神经元群的活动可能无差异；2、ApEn可能不能反映该状态间神经元群活动的差异。验证了临床癫痫患者脑电信号的ApEn分析结果193随时间演化的近似熵分析 三种脑功能状态下约45s脑电信号的近似熵随时间的演化20分频段分析不同脑功能状态间转换时无痫样放电 连续性痫样放电连续性痫样放电 周期性痫样放电左右ECoG的ApEn显著增加p0.01左右ECoG和EHG的高节律的ApEn显著降低p0.05左右ECoG和EHG的ApEn显著降低p0.05脑功能状态转换时的ApE

7、n主要取决于较高节律ApEn的变化癫痫波发放终止可能与低节律脑电活动同步增加有关节律快波24.5Hz成分的相对含量显著增多p0.05能量分析脑功能状态转换时的ApEn主要取决于较高节律ApEn的变化节律的快波24.5Hz成分增多？21近似熵分析总结：0.530Hz的分析：证实了痫样放电时相关脑区神经元群活动从复杂的混沌状态趋向于简单的同步活动；随时间演化的ApEn分析进一步证明了近似熵可以鉴别背景脑电信号和痫样放电；分频段的分析不同脑功能状态间转换时 ：连续性痫样放电时ECoG中节律的快波成分显著增多；可以推测在脑功能状态发生转换时各脑区脑电信号ApEn的变化主要取决于高节律ApEn的变化；痫

8、样放电的终止，可能与低节律脑电活动的同步性增强有关；由、，我们推测节律在痫样放电过程中可能起着更重要的作用。 22似然同步分析参数取值： 窗口宽度数据长度根据数据处理的实际需要选取；并将计算结果在时间及空间上进行平均 ；0.530Hz脑电信号分析结果分频段脑电信号分析结果230.530Hz脑电信号分析结果无痫样放电与连续性痫样放电 连续性痫样放电与周期性痫样放电 LCLH *SynonSyperiod RCRH *SynonSycontinuous Sycontinuous=Syperiod LCRC *SynonSycontinuous Sycontinuous=Syperiod LHRH

9、*SynonSyperiod 注：*P0.01；*P0.052不同脑功能状态之间转换 ：脑区间的同步性有显著性差异 p0.05注：*P0.01；*P0.05 24由表可见：不同脑功能状态间转换时无痫样放电 连续性痫样放电连续性痫样放电 周期性痫样放电LCLH、LHRH p0.01 RCRH、LCRCp0.05的同步性均显著增强LCLH、LHRH的同步性显著降低p0.05反映了痫样放电的产生是由大量神经元群过度同步放电所致。 由匹罗卡品诱导的痫样放电活动中，左侧海马可能起着更重要的作用 左侧优势现象253 随时间演化的似然同步分析约33s三种脑功能状态的似然同步随时间的演化261相同脑功能状态

10、分频段脑电信号分析结果无痫样放电 连续性痫样放电 周期性痫样放电 与 SySy *SySy *SySy Sy=Sy Sy=Sy 与 *SySy Sy=Sy Sy=Sy 注：*P0.01；*P0.05 27由表可见：相同脑功能状态无痫样放电连续性痫样放电周期性痫样放电节律的同步性显著高于、节律的同步性节律的同步性显著低于、节律的同步性各节律的同步性间无显著性差异在痫样放电前节律可能已先于其它节律产生了类似于痫样放电的同步现象 痫样放电过程中较高节律的同步化可能占优势节律与其它节律的同步性差异可以区分不同的脑功能状态282不同脑功能状态间转换时无痫样放电与连续性痫样放电 连续性痫样放电与周期性痫样

11、放电 节律 *SynonSycontinuous Sycontinuous=Syperiod 节律 *SynonSyperiod 节律 *SynonSyperiod 节律 *SynonSyperiod 节律 *SynonSyperiod 注：*P0.01；*P0.05 29由表可见：不同脑功能状态间转换时无痫样放电 连续性痫样放电连续性痫样放电 周期性痫样放电左右皮层和海马脑电信号各节律的同步性均显著增强节律的同步性无显著性变化其余各节律的同步性均显著降低 低节律、段的同步性增强的程度低于高节律、段 在痫样放电过程中脑电信号高节律、同步性增强可能占主导地位 在痫样放电后期，高节律脑电活动的复杂

12、度增加早于痫样放电的终止 结合相同脑功能状态下的分析结果30似然同步分析总结：0.530Hz的分析痫样放电的产生是由大量神经元群过度同步放电所致；皮罗卡品诱发的痫样放电大鼠模型具有左侧优势现象；可以动态地刻画痫样放电过程中脑区间的同步性演化过程 ；分频段分析1相同脑功能状态下痫样放电前节律可能已先于其它节律产生了类似于痫样放电的同步现象；痫样放电过程中较高节律的同步化可能占优势；节律与其它节律的同步性差异可以区分不同的脑功能状态。 2不同脑功能状态间转换时痫样放电过程中脑电信号高节律、 同步性增强可能占主导地位；节律在痫样放电结束过程中可能起着重要作用；31 Paper1Yang Xia; H

13、ongmei Yan; Yansu Liu刘衍素; Dezhong Yao; Andrew C.N Chen.Hippocampus: the origin of pilocarpine-induced epilepsy. 12th ECCN congress.Stockholm. 2005.5. 录用2Xia Yang; Liu Yansu刘衍素; Yan Hongmei; Wu Hongyi; You Zili; Yao Dezhong. Synchronization analysis of EEG signal between cortex and hippocampus of pil

14、ocarpine-induced epileptic rat. 27th Annual Iternational Conference of the IEEE Engineering in Medicine and Biology Society (EMBS). ShangHai 2005.9. (Submit)3痫样放电大鼠颅内脑电信号的同步分析4癫痫大鼠颅内脑电信号不同节律的同步分析32本文的创新点不同频段脑电信号的近似熵分析；似然同步用于大鼠模型脑电信号分析；研究结果提示在该模型中，可能存在痫样放电的左侧优势现象，这一现象在癫痫动物模型研究中尚未见报道，值得进一步研究。33.工作总结和展

15、望本文工作总结1. 将近似熵分析方法应用于临床癫痫患者的脑电信号分析，结果显示近似熵分析能够用于癫痫患者的背景脑电信号和痫样放电的鉴别检测；2. 成功建立了匹罗卡品诱发痫样放电大鼠模型，在本实验室首次探索了痫样放电大鼠模型的实验研究，为进一步深入开展癫痫电生理实验研究奠定了根底；3. 通过近似熵、似然同步方法分析痫样放电大鼠模型的脑电信号，提示这两种方法都可以刻画痫样放电过程中脑电信号的非线性动力学过程。 34对今后工作的展望1.节律脑电信号在痫样放电过程中的作用有待于进一步研究；2. 结合分子生物学研究手段，研究匹罗卡品诱发的痫样放电大鼠模型癫痫波发放的神经电生理机制与基因表达的关系；3.

16、基于动物模型，并与临床密切结合，进行癫痫预测与控制的新方法研究。 35与本研究有关的论文1WU Hongyi; XIA Yang; LAI Yongxiu; LIU Yansu刘衍素;Yao dezhong. Study of Epileptic EEG Using Bispectrum Analysis. Proceedings of International Conference on Neural Interface and Control. 26-28 may 2005. WuHan China. 119-122. 2Hongyi Wu; Yang Xia; Hongmei Yan;

17、 Yansu Liu刘衍素 and Dezhong Yao. Bispectrum analysis and its normalization of epileptic Rats EEG. 27th Annual Iternational Conference of the IEEE Engineering in Medicine and Biology Society (EMBS). ShangHai. 2005.9. (Submit) 36. 致谢首先衷心感谢我的导师夏阳副教授在学习、实验及论文写作中给予我的悉心指导和帮助。特别感谢尧德中教授给予我的关心、帮助和指导。感谢生命学院各位领导

18、和老师在我学习期间给予的无私帮助和指导。感谢本教研室的吴宏伊、王玲、徐鹏、赵丽娜等同学对我完本钱课题所提供的无私帮助。还要感谢我的家人对我的关心、理解、支持和鼓励。 37陈述完毕谢谢各位专家的评审！38RiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9Ld

19、OgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%

20、v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8J

21、bNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq

22、!t*w-z1D4G7JbMeQhThTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3

23、E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G8KbNeQiT

24、lWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSkVnYq$t*w-A

25、1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8K8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLd

26、OgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMePhSkWnZq$u*x+A2D5H8KbNfQiUlXo#s%

27、v(y0B3F6I9LdOgRjVmYq!t!t&w-z1C4G7JbMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTlWoZr%u(x+B2E6H9KcOfRiUmXp!s&v)z0C3F7IaMdPgSkVnYq$t*x-A1D5G8JbNeQiTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A2D5G

28、8KbNeQiTlXo#r%v(y+B3E6I9LcOgRjUmYp!t&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LdOgRjVmYp!t&w-z1C4G7JaMePhTkWnZr$u*x+A2E5H8KcNfQiUlXp#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H9KcNfRiUlXp#s&v)y0C3F6IaLdPgSjVnYq!t*w-A1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPdPgSk

29、VnYq$t*w-A1D5G8JbNeQhTlWo#r%u(y+B2E6H9LcOfRjUmXp!s&w)z0C4F7IaMdPhSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%v(y+B3E6H9LcOgRjUmYp!s&w)z1C4F7JaMdPhSkWnZq$u*x-A2D5H8KbNfQiTlXo#s%v(y0B3E6I9LcOgRjVmYp!t&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KcNfQiUlXo#s%v)y0B3F6I9LdOgSjVmYq!t&w-z1D4G7JbMePhTkWoZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0

30、C3F6IaLdOgSjVnYq!t*w-z1D4G8JbMeQhTkWoZr%u(x+B2E5H9KcOfRiUmXp#s&v)z0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTlWoZr%u(y+B2E6H9KcOfRjUmXp!s&v)z0C4F7IaMdPgSkVnZq$t*x-A1D5G8KbNeQiTlWo#r%u(y+B3E6H9LcOfRjUmYp!s&w)z0C4F7JaMdPhSkVnZq$u*x-A2D5G8KbNfQiTlXo#r%v(y0B3E6I9LcOgRjVmYp!t&w&w)z1C4G7JaMePhSkWnZr$u*x+A2D5H8KbNfQiUlXo#s%v(y0B3F6I9LdOgRjVmYq!t&w-z1C4G7JbMePhTkWnZr$u(x+A2E5H8KcNfRiUlXp#s%v)y0C3F6IaLdOgSjVmYq!t*w-z1D4G7JbMeQhTkWoZr$u(x+B2E5H9KcNfRiUmXp#s&v)y0C3F7IaLdPgSjVnYq$t*w-A1D4G8JbNeQhTl