我得到的分岔图怎么是这样？

ruichard101

我分析的是一个简单的卧式jeffcott转子模型，同时考虑了油膜压力和剪力，得到的分叉图怎么是这样呢？我改了很多次初值，基本上都是这个形状。请问一下我问题出在哪了？

octopussheng

这个图建议你在5~50这个区域好好的算一下，可能会找到一些现象的！

ruichard101

谢谢oct，我试试！~~~

ruichard101

我算过了，和以前的基本没什么区别:@L
无论我怎么取初值，得到的解都是虚解。可虚部到后来显示的是0.0000i:@( ，我总觉得我的主程序有问题:@o

无水1324

对，应该是程序出了问题

xinwilliam

因为了转子模型，有可能是磨擦性或碰撞性分岔!

无水1324

有点奇怪的是从坐标的幅值很大10^12

octopussheng

怀疑他的系统是不是不收敛了！

ruichard101

算w的时候公式出了错误:loveliness:
最后出来的图如下：

ruichard101

clc;
clear;
global M w
m=2;
C=0.0001;
eta=0.018;
R=0.025;
L=0.01;
w=M*eta*R*L*(R/C)^2*(L/(2*R))^2/(m*C);
M=10:300;
x0=[4;4;4;4]
for h=1:length(M)
    T=2*pi/(M(h)*eta*R*L*(R/C)^2*(L/(2*R))^2/(m*C));
    [t y]=ode45('youmofun',[0:T/100:200*T],x0,[],M(h));
    y(:,1)
    plot(M(h),y(100*180:100:end,1),'k.');hold on
end
初值取小了是空解，我从一试到了十，只有在4以后才出现解，而且还都是虚数:@L
到底是哪里出问题了呢？

ruichard101

原帖由 octopussheng 于 2008-4-10 15:22 发表

                               
登录/注册后可看大图

怀疑他的系统是不是不收敛了！

应该是收敛呀！分叉图都有了，只是我得不出那样的图:@(
哎！得到一个像样的分叉图咋这么难啊！:@o

[ 本帖最后由 ruichard101 于 2008-4-13 19:12 编辑 ]

superliu

这可以肯定不是分岔图，一个分岔都没有

无水1324

弱问一下，你的这个M是频率吗？

ruichard101

不是频率,是无量纲质量,它是关于频率的函数.

无水1324

神啊！

你不是频率就不能用这种方法做分岔图，你要理解频闪法作分岔图的物理意义哈，是分频采样。

[ 本帖最后由 咕噜噜 于 2008-4-14 19:10 编辑 ]