麻烦大家也帮我看看这几个图

无水1324

那计算的时间与参数应该没有什么关系阿；

还有你画分岔图是用频闪法，还是？总感觉你那个分岔图有问题

jinqiankun

程序：还是从论坛上找的
global w;
wd=1.0;
range=[500:2:2100];
period=2*pi/wd;      
k=0;
step=2*pi/512;      
for w=range
    y0=[0.9;0.002;-0.4;0.002;0.9;0.002;-0.4;0.002;0.9;0.002;-0.4;0.002;-0.4;0.002];
    ww=w
    k=k+1;
    tspan=[0:step:50*period];
    [t,Y]=ode45('shongdongprocedureLastm3_100',tspan,y0);
    y0=(Y(end,:))';
    j=1;
    for i=50:100
        tspan=[i*period:step:(i+1)*period];
        [t,Y]=ode45('shongdongprocedureLastm3_100',tspan,y0);
        YY1(k,j)=Y(end,1);   % get the omega data from every period end
        YY2(k,j)=Y(end,13);
        j=j+1;               %取出每一个周期内的第一个解的最后一个值。
        y0=(Y(end,:))';
        ii=i
        
    end
   
end

[ 本帖最后由 jinqiankun 于 2007-12-14 10:38 编辑 ]

无水1324

周期应该与w有关系吧？

jinqiankun

无量纲后 周期都是2* pi
tao=w*t