大家看一下此Poincare图是否反应混沌运动？

octopussheng

做映射图时，选点的步长是保证图正确与否的标准啊！不能随便选择的！

quiet1016

quiet1016

quiet1016

function ff
options = odeset('RelTol',1e-5,'AbsTol',[1e-5 1e-5]);
sol= ode45(@ischao,[0 1000],[0 0],options);  %[T,Y] =
step=7.5;    %步长
x = 0:0.1:1000;
y1 = deval(sol,x,1);
y2 = deval(sol,x,2);
%plot(y1(100:step:end),y2(100:step:end),'*')     %y1 vs y2
subplot(1,2,1)
plot(y1(1:step:end-step),y1(1+step:step:end),'*');  %y1(t-1) vs y1(t)
subplot(1,2,2)
plot(y2(1:step:end-step),y2(1+step:step:end),'*');  %y2(t-1) vs y2(t)
function dy = ischao(t,y)
dy = zeros(2,1);    % a column vector
dy=[y(2);15.0*cos(1.0*t)-0.375*y(1)^3-y(1)+1.299*y(1)^2-0.0059*y(2)];

quiet1016

quiet1016

好像怎么都得不到楼主贴的那个POINCARE图

octopussheng

呵呵，确实啊！我也没算到！

16443

[quote]subplot(1,2,1)
plot(y1(1:step:end-step),y1(1+step:step:end),'*');  %y1(t-1) vs y1(t)
subplot(1,2,2)
plot(y2(1:step:end-step),y2(1+step:step:end),'*');  %y2(t-1) vs y2(t)
[quote]

这样做相图是什么意思啊 ？

无水1324

你这个步长选择的依据是什么？
Poincare图示根据什么原理做出来的？

16443

就是么，给介绍一下1

无水1324

根据现有的两种做Poincare方法，他这个都不能很好的解释。难道是一种新的理论算法，请quiet1016不吝赐教！

quiet1016

不是做Poincare，是y(t-x)与y(t)的关系，其中x＝7.5,我就是想看看楼主之前给的图究竟是怎么得到的

lzhyxcx

对不起，我的方程与第一个Poincare图不符，其中一个参数不一样。谢谢各位的讨论！

无水1324

不会吧。:@):@L

zhuoer

建议将系统的Lyapunov指数求出来，更能准确的判断是否进入了混沌