求高手帮忙画一个分岔图

dorje

最近需要画一个系统的分岔图， 现学了一点matcont, 但是画到分岔图时老提示错误，请高手帮忙给看一下：

系统为：
x'=(a*(1-x)/(e+z))-(b*z*x/(d+x))
y'=x*x-y                        
z'=y-(c*(1-x)*z/(e+z))

其中参数取值： a=0.007; b=954.5;c=0.035;d=0.029;e=0.00002

我算出两个霍普夫分岔点为：
first point found
tangent vector to first point found
label = H , x = ( 0.073637 0.005422 0.003067 0.005892 )
First Lyapunov coefficient = -7.371038e+001
label = H , x = ( 0.332767 0.110734 0.002296 0.167406 )
First Lyapunov coefficient = -3.494008e+001

elapsed time  = 15.9 secs
npoints curve = 300

可是画分岔图的时候总是报错： time step too small。 烦请高手指点一下， 最好能给画一下，我需要的是x,z 以c为控制参数的分岔图， 并且需要系统周期随c 的变化关系

hsfy919

回复 1 # dorje 的帖子

时间不长太小，有可能是方程或积分本身的问题，比如方程刚性。最好你看看matcont里面的help

kangarooli

我之前也看过一点matcont，尝试画了下分岔图，但是看不懂怎么回事，不过到没遇到你说的这种错误，你是怎么画的呢，我是按着里面画的，不知对你是否有用
Matcont : A Matlab package for dynamical systems with applications to neural activity