求高手帮忙画一个分岔图
最近需要画一个系统的分岔图, 现学了一点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 的变化关系
回复 1 # dorje 的帖子
时间不长太小,有可能是方程或积分本身的问题,比如方程刚性。最好你看看matcont里面的help 我之前也看过一点matcont,尝试画了下分岔图,但是看不懂怎么回事,不过到没遇到你说的这种错误,你是怎么画的呢,我是按着里面画的,不知对你是否有用
Matcont : A Matlab package for dynamical systems with applications to neural activity
页:
[1]