下面的分岔图是否正确

chuandong418

恩，偶页想不明白，是不是运行程序有问题，还望各位高手指点，谢谢了

无水1324

原帖由 无水1324 于 2007-8-15 12:54 发表

                               
登录/注册后可看大图

应该是书写错误，必须要ff=[dx(1);dx(2);dx(3);dx(4)]的

我都想不通你们为什么说上面这个ff可以省略？

yzsldj

以w为分岔参数的分岔图

chuandong418

谢谢yzslj这位大侠！

[ 本帖最后由 咕噜噜 于 2007-8-17 18:51 编辑 ]

咕噜噜

开始的那些应该是瞬态吧，你过滤掉瞬态没有？

yzsldj

已经过滤掉瞬态了，共迭代30000次，只取了最后的1500次的结果，
关键是把原来的非自治系统先改为自治系统。

octopussheng

原帖由 yzsldj 于 2007-8-18 08:30 发表

                               
登录/注册后可看大图

已经过滤掉瞬态了，共迭代30000次，只取了最后的1500次的结果，
关键是把原来的非自治系统先改为自治系统。

关键是把原来的非自治系统先改为自治系统？

这个是什么意思，能不能解释一下啊？

咕噜噜

:lol 我也想知道为什么，yzsldj 快出来
为什么要这么做

yzsldj

将原4阶的非自治系统改为如下5阶的自治系统
F =[x(3);
       x(4);
        -2*E*x(3)-x(1)-a*(x(1)^2+x(2)^2)*x(1)-b*(1-1/sqrt(x(1)^2+x(2)^2))*(x(1)-f*x(2))+u*w^2*cos(w*x(5));
        -2*E*x(4)-x(2)-a*(x(1)^2+x(2)^2)*x(2)-b*(1-1/sqrt(x(1)^2+x(2)^2))*(f*x(1)+x(2))+u*w^2*sin(w*x(5))-G;
        1];  
再考虑分岔图，就简单了

咕噜噜

能不能说一说为什么这么做

yzsldj

这样做的目的是为了便于运用局部最大值方法作分岔图，
因为我不知道对于非自治系统有什么好的方法来作分岔图，
我一直是先把非自治系统改为自治系统，然后再运用局部最大值方法作分岔图。

无水1324

有点意思，第一次看到这种说法。

不过，可能对于多频激励的情况可能用起来比较好

yzsldj

原帖由 无水1324 于 2007-8-18 18:07 发表

                               
登录/注册后可看大图

有点意思，第一次看到这种说法。

不过，可能对于多频激励的情况可能用起来比较好

其实，在很多教材上都有类似的介绍。

无水1324

是吗？
能否推荐几本看看

octopussheng

这个想法确实不错，在有些书上确实有介绍，不过以前还真不想到这方面去

在张琪昌的那本书上将中心流形定理的时候，推导中心流形的时候就有说过，将时间t看成一个变量，将它增加到系统方程里面去，相当于给系统多加了一个自由度，无水你看看！