计算相位响应曲线的方波扰动直接算法

通过相位响应曲线可对具有极限环周期运动的动力系统的性质有更为深入的理解.神经元是一个典型的动力系统,因此相位响应曲线提供了一种研究神经元重复周期放电行为的新思路.本文提出一种求解相位响应曲线的方法,即方波扰动的直接算法,通过Hodgkin-Huxley,Fitz Hugh-Nagumo,Morris-Lecar和Hindmarsh-Rose神经元模型验证该算法可计算周期峰放电、周期簇放电的相位响应曲线.该算法克服了其他算法在运用过程中的局限性.利用该算法计算结果表明:周期峰放电的相位响应曲线类型是由其分岔类型所决定;在Morris-Lecar模型中发现一种开始于Hopf分岔终止于鞍点同宿轨道分岔的阈上周期振荡,其相位响应曲线属于第二类型.通过大量的相位响应曲线的计算发现相位响应的相对大小及正负性仅取决于扰动所施加的时间,而且周期簇放电的相位响应曲线比周期峰放电的相位响应曲线更为复杂.