请教关于分岔程序

gghhjj

原帖由 无水1324 于 2007-7-12 10:06 发表

                               
登录/注册后可看大图

这个大概的意思是：周期运动的话，会出现相同的最大值，其实道理上与 闪频差不多，这是这种方法在系统的解出现跳跃的时候，是不准确的。个人的想法

个人认为这种方法是不可取的
1. 正如无水所说，这个方法只能用于个别系统，对于未知系统根本无法判断结果是否正确
2. 这种方法效率十分低

octopussheng

确实是有一定的问题，今天我也用一下无水上面的方法试试，呵呵

octopussheng

:@L :@L :@L :@L :@L

无水，我用你前面写的程序画出来的分岔图怎么是这样的？？？

为避免你翻页看麻烦，把程序贴出来了
function xdot=dddd(t,x,flag,omega)
eta=0.1 ; a= -48.704; b= 24.35e4;gamma= 1;
xdot=[x(2);a*x(1)-b*x(1)^3-eta*x(2)+gamma*cos(omega*t)];


clear;clc;close all;
omega=1:0.01:2;
for h=1:length(omega)
    T=2*pi/omega(h);
    [t,x]=ode45('dddd',[0:T/100:220*T],[0 0],[],omega(h));
    plot(omega(h),x(10000:100:end,2),'k.');hold on
end

无水1324

是这样的也很正常，你再去看一下这类系统出现混沌是在那个区域。将omega调整到那个区域就可以了

octopussheng

呵呵，方法是可以算的，和分岔参数的设置有关系，我再算算吧！

octopussheng

楼主看看这个分岔图吧，算了一个小时左右，看看效果如何～

mechanic05

前辈你好！你能关心这个问题我很高兴！！谢谢！
就你画出的分岔图来说代表什摸？其实我的问题是a的前面有负号．如果去掉的话系统应该是周期运动，也就是说你画出的分岔图应该代表周期运动，不会发生分岔．你的图怎摸解释？麻烦再把它换成－a试试．谢谢！:@)

无水1324

a= -48.704,他的这个方程很奇怪的

mechanic05

前辈好！
我的问题是一个梁的大挠度问题．其中a的值在大于零时是稳定的，不会屈曲．a在小于零时出现屈曲．
前辈，我画出的分岔图瞬态响应太大．就原来的程序我把 plot(omega(h),x(10000:100:end,2),'k.');hold on　修改成 plot(omega(h),x(100000:100:end,2),'k.');hold on为什摸老出现矩阵长度不符的问题？:@L

咕噜噜

哈，这个图好像人脸啊
不过这个方程很奇怪

咕噜噜

omega与x矩阵大小看一看是不是一致，主要是后面画图的时候

mechanic05

前辈好！
我把１００００改成１０００００，应该可以减小瞬态响应吧？:@)

无水1324

是的
但是你计算的总点数都没有这么多哈

mechanic05

:@L :@L :@L
我刚装上７．０，怎么同样的程序６．５好象更快一些？
我用７．０已经算了一个小时了，还在算．:lol

[ 本帖最后由 无水1324 于 2007-7-16 11:52 编辑 ]

无水1324

原帖由 mechanic05 于 2007-7-16 11:34 发表

                               
登录/注册后可看大图

:@L :@L :@L
我刚装上７．０，怎么同样的程序６．５好象更快一些？
我用７．０已经算了一个小时了，还在算．:lol

其实6.5的计算功能还要好一些，只是一些功能、工具箱有限