yiquncw 发表于 2009-12-12 13:16

帮忙分析一下分岔图结果

下图是用最大值法计算得到的分岔图，请问能得出什么结论？
图形看着周期与非周期运动区分不开，是程序的问题，还是其他的原因？
系统为在偏心质量、重力、陀螺力矩作用下的转子系统，无碰磨，横坐标为转角速度
请高手赐教！


[ 本帖最后由 yiquncw 于 2009-12-12 13:18 编辑 ]

无水1324 发表于 2009-12-12 15:47

应该是程序的问题了

yiquncw 发表于 2009-12-12 15:54

回复 沙发 无水1324 的帖子

function norub_bifur
q0=
options = odeset('RelTol',1e-6,'AbsTol',1e-6)
T=2*pi
ts=
for w=linspace(1,2501,250)
    =ode45('norub',ts,q0,options,w)
=getmax(z(:,1))
plot(w,Xmax,'b','markersize',1)
    hold on
    clear Xmax
end

这是我的程序，getmax子程序就是论坛里给出求最大值的方法，
是哪里的问题呢？

yiquncw 发表于 2009-12-12 16:06


另外补充一下，上图是用x=y截面图法得到的分岔图，也很不理想

yiquncw 发表于 2009-12-12 22:13

请斑竹们赐教啊，下面是算的弯摆耦合转子振动分岔图，在某转速下画的poincare图都是有奇异吸引子的，可分岔图没有对应的混沌信息啊


分岔图用的关键程序语句为：
for n=1:length(w)
    T=2*pi
    ts=
    =ode45('norub',ts,q0,[],w(n))
    q0=z(end,:)
    figure(1)
    subplot(2,1,1)
    plot(w(n),z(40000:50:end,1),'k.','markersize',1)
应该是频闪法吧。
请斑竹们找找原因，给点思路。

yiquncw 发表于 2009-12-13 09:45

等待指教...

yibanlixue 发表于 2009-12-13 11:35

plot(w(n),z(40000:50:end,1),'k.','markersize',1)

频闪法应该是plot(w(n),z(40000:500:end,1),'k.','markersize',1)

yiquncw 发表于 2009-12-13 13:36

我再算一下，有好多结果再贴出来
页: [1]
查看完整版本: 帮忙分析一下分岔图结果