ruichard101 发表于 2008-4-8 19:04

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

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

octopussheng 发表于 2008-4-8 20:39

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

ruichard101 发表于 2008-4-8 21:03

回复 2楼 的帖子

谢谢oct，我试试！~~~

ruichard101 发表于 2008-4-9 16:49

回复 2楼 的帖子

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

无水1324 发表于 2008-4-9 18:13

回复 4楼 的帖子

对，应该是程序出了问题

xinwilliam 发表于 2008-4-9 18:45

回复 5楼 的帖子

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

无水1324 发表于 2008-4-10 14:04

回复 6楼 的帖子

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

octopussheng 发表于 2008-4-10 15:22

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

ruichard101 发表于 2008-4-13 18:48

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

ruichard101 发表于 2008-4-13 19:03

回复 5楼 的帖子

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=
for h=1:length(M)
    T=2*pi/(M(h)*eta*R*L*(R/C)^2*(L/(2*R))^2/(m*C));
    =ode45('youmofun',,x0,[],M(h));
    y(:,1)
    plot(M(h),y(100*180:100:end,1),'k.');hold on
end
初值取小了是空解，我从一试到了十，只有在4以后才出现解，而且还都是虚数:@L
到底是哪里出问题了呢？

ruichard101 发表于 2008-4-13 19:10

原帖由 octopussheng 于 2008-4-10 15:22 发表 http://www.chinavib.com/forum/images/common/back.gif
怀疑他的系统是不是不收敛了！
应该是收敛呀！分叉图都有了，只是我得不出那样的图:@(
哎！得到一个像样的分叉图咋这么难啊！:@o

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

superliu 发表于 2008-4-13 21:05

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

无水1324 发表于 2008-4-14 18:55

回复 11楼 的帖子

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

ruichard101 发表于 2008-4-14 18:57

回复 13楼 的帖子

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

无水1324 发表于 2008-4-14 19:06

回复 14楼 的帖子

神啊！

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

[ 本帖最后由 咕噜噜 于 2008-4-14 19:10 编辑 ]
页: [1] 2
查看完整版本: 我得到的分岔图怎么是这样？