octopussheng 发表于 2007-8-19 11:58

回复 #75 无水1324 的帖子

从理论上看确实是这样,不过如果chuandong418 的系统做非自治的分岔图不好做的话,变换成自治的也是可行的啊,毕竟也是一条路~

无水1324 发表于 2007-8-19 12:04

回复 #76 octopussheng 的帖子

不可能非自治不好做(对于他的系统)

chuandong418 发表于 2007-8-19 13:03

回复 #69 yzsldj 的帖子

下面这个图是按照yzldj所说的方法做出的以b为分岔参数的分岔图,但是感觉很不对劲,因为偶做的poincare显示在b在3到20之间出现了一个周期运动,一个概周期运动, 还有混沌运动,感觉这个分岔图不正确,请大家帮忙分析!而且从它上面看不出任何多周期运动!

咕噜噜 发表于 2007-8-19 14:22

回复 #78 chuandong418 的帖子

是有点奇怪哈
其实我觉得非自治系统的分叉应该更好做才是

yzsldj 发表于 2007-8-19 14:32

回复 #78 chuandong418 的帖子

在b=6~9之间作图,相当于将上图局部放大,应该能看到多倍周期运动。

无水1324 发表于 2007-8-19 14:53

回复 #79 咕噜噜 的帖子

赞成!

chuandong418 发表于 2007-8-19 15:10

回复 #79 咕噜噜 的帖子

呵呵,但是如果按非自治系统去计算,matlab根本没反应!不知道为什么,只有将其按yzsldj的说法,将其转化为自治系统后才能运行!很奇怪,也希望yzsldj解释下这是为啥?

yzsldj 发表于 2007-8-19 15:21

回复 #82 chuandong418 的帖子

说实话,你放在#1楼的那个程序我没看懂,我以前也没碰到过非自治系统,我只是觉得转成自治系统来处理(孰优孰劣暂且不论),至少方法应该是可行

咕噜噜 发表于 2007-8-19 15:27

回复 #82 chuandong418 的帖子

你66楼的程序有错误

chuandong418 发表于 2007-8-19 17:38

回复 #84 咕噜噜 的帖子

修改正确后才去做的,呵呵,谢谢啊咕:loveliness:

咕噜噜 发表于 2007-8-19 18:32

哦?修改后的程序贴上来看看
我帮你运行一下

chuandong418 发表于 2007-8-20 17:14

回复 #86 咕噜噜 的帖子

b=2.0:0.1:20;
options = odeset('RelTol',1e-6,'AbsTol',);
forn=1:length(b);
w=2.0;
T=2*pi/w;
ts=;
x0=;
=ode45('www',ts,x0,options,b(n));
figure(1)
plot(b(n),X(5000:100:end,1),'k.');
xlabel('\fontsize{18}\omega');
ylabel('\fontsize{18}x');
hold on
end
grid



function ff=www(t,x,flag,b)
ff=zeros(5,1);
g=9800;
E=0.12;
f=0.12;
a=0.5;
%b=15.0;
q=0.16;
f0=25;
u=0.136/q;
w=2.0;
e=sqrt(x(1)^2+x(2)^2);
G=g/((2*pi*f0)^2*q);
ff=zeros(5,1);
dx(1)=x(3);
dx(2)=x(4);
dx(3)=-2*E*x(3)-x(1)-a*(x(1)^2+x(2)^2)*x(1)-b*(1-1/e)*(x(1)-f*x(2))+u*w^2*cos(w*x(5));
dx(4)=-2*E*x(4)-x(2)-a*(x(1)^2+x(2)^2)*x(2)-b*(1-1/e)*(f*x(1)+x(2))+u*w^2*sin(w*x(5))-G;
dx(5)=1;
ff=;

octopussheng 发表于 2007-8-20 20:09

帮你算了一下,怎么和上面的图一样啊,有没有修改??

chuandong418 发表于 2007-8-20 20:42

回复 #88 octopussheng 的帖子

就是修改后做出来的图形啊!
对了,你用过let来算李亚扑懦夫指数把?
我用得时候设置initial conditions好象提示输入格式不对?
这个位置的格式是什么样的?能否告知@谢谢哦@:loveliness:

chuandong418 发表于 2007-8-20 20:46

回复 #88 octopussheng 的帖子

还有let工具箱中为什么没有分岔参数的指定啊?
还有下面的设置intial conditions 的说明,看不懂,帮指点下!谢谢了
The number of initial conditions must be equal to the total number
%       of state variables of the system.Users need not to specify the
%       initial conditions for the variational equation since the program
%       can generate the initial conditions for the variational equation as
%       long as the users provide the correct number of linearized ODEs.
页: 1 2 3 4 5 [6] 7
查看完整版本: 下面的分岔图是否正确