dorje 发表于 2011-1-22 03:32

求高手帮忙画一个分岔图

最近需要画一个系统的分岔图, 现学了一点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 发表于 2011-1-22 11:44

回复 1 # dorje 的帖子

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

kangarooli 发表于 2011-2-23 10:02

我之前也看过一点matcont,尝试画了下分岔图,但是看不懂怎么回事,不过到没遇到你说的这种错误,你是怎么画的呢,我是按着里面画的,不知对你是否有用
Matcont : A Matlab package for dynamical systems with applications to neural activity
页: [1]
查看完整版本: 求高手帮忙画一个分岔图