论文部分内容阅读
从20世纪80年代,随着非线性动力学理论的发展,非线性动力学被广泛地用于分析癫痫脑电信号。然而,由于基于混沌理论的非线性动力学分析方法需要假设脑电是低维混沌信号,对脑电信号的数据长度及平稳度要求较高,且对脑电信号中的噪声很敏感。为了克服传统非线性动力学方法的缺点,本文立足于癫痫脑电信号分析的要求,提出了新的方法用于刻画癫痫发作各阶段脑电信号的动力学特征变化。首先,结合传统递归图和排序递归图提出了混杂递归图方法分析失神发作脑电信号。该方法的新颖性在于定义时间序列的递归状态时,不仅考虑局部相空间距离而且考虑局部排序模式分布结构。仿真分析表明,基于混杂递归图的对角线结构分布的确定性测度DET,能够刻画模型参数的变化,并且抗噪声干扰能力更强。该方法用于分析大鼠失神发作间隙期、发作前期和发作期脑电信号的确定性特征,发现发作前期脑电信号的DET均值显著地大于发作间隙期,但显著地小于发作期。结果表明大鼠失神性癫痫发作前期脑电的确定性特性程度比发作间隙期脑电更高,但是比发作期脑电更低。其次,为了进一步调查失神性癫痫发作各阶段脑电信号隐藏的非线性动力学特征,基于脑电信号的排序模式分布特性,提出了一个新颖的相异性指数方法。该方法通过计算排序模式分布的距离来分析两段脑电信号的相异性。由神经元群模型模拟生成脑电信号,仿真分析相异性指数的性能。结合移动窗口技术,该方法用于分析110段大鼠失神发作脑电信号,其中58段失神发作脑电信号在其发作前相异性指数有显著地增加,能够成功地检测到发作前期状态,表明该方法可以用于检测失神发作脑电信号的动力学特征变化。再次,针对脑电同步分析在癫痫发作研究中的重要意义,基于排序模式分析提出了一个新颖的互信息估计方法。该方法通过对时间序列排序模式进行分类,来实现复杂的概率分布估计,从而直接估计出时间序列的信息量。通过耦合Henon映射模型和耦合神经元群模型模拟生成耦合时间序列,仿真分析排序互信息方法的性能。相比于传统的直方图方法,基于排序模式的互信息估计方法能够更好地刻画模型耦合系数的变化,而且抗噪声干扰能力更强。该方法应用于分析癫痫发作脑电信号,发现脑电信号间的互信息随着癫痫发作的开始而逐渐增加,在完全发作时达到最大。最后,提出了排序自互信息法用于刻画癫痫发作间隙期、发作前期和发作期脑电的动力学特征。仿真分析发现噪声和混沌序列的自互信息AMI(δ)值随着时滞偏移δ增大而衰减的速率不同,表明自互信息方法能够区别噪声和混沌序列。该方法用于分析癫痫发作各阶段脑电信号,发现发作间隙期、发作前期和发作期脑电的自互信息值分别在时滞偏移δ大约为5-6、7-8和9-10时达到稳定值。进一步结合线性判别分析,表明排序自互信息能够识别不同发作阶段的脑电信号。