无水1324 发表于 2007-12-13 22:59

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

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

jinqiankun 发表于 2007-12-14 10:37

回复 #16 无水1324 的帖子

程序：还是从论坛上找的
global w;
wd=1.0;
range=;
period=2*pi/wd;      
k=0;
step=2*pi/512;      
for w=range
    y0=;
    ww=w
    k=k+1;
    tspan=;
    =ode45('shongdongprocedureLastm3_100',tspan,y0);
    y0=(Y(end,:))';
    j=1;
    for i=50:100
      tspan=;
      =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 发表于 2007-12-14 10:41

周期应该与w有关系吧？

jinqiankun 发表于 2007-12-14 10:45

回复 #18 无水1324 的帖子

无量纲后 周期都是2* pi
tao=w*t
页: 1 [2]
查看完整版本: 麻烦大家也帮我看看这几个图