请问这样的poincare映射图和lyapunov指数图能说明系统是混沌的吗？

ssaw0920

请问各位学长：
对于同一个动力系统得到这样的poincare映射图和lyapunov指数图能说明系统是混沌的吗？

siyanyu

从最大Lyapunov指数来看，可是说是混沌的

ssaw0920

但是第一、二个图近似为直线，不知道是不是准周期的？

liliangbiao

指数怎么这么大呢？？

wxl789

从lyapunov指数来看系统是混沌的

你的poincare截面范围好像太大了，把它缩小。然后再看看

ssaw0920

指数最小有500多
另外请问wxl789学长，怎样把截面范围变小呢？

liliangbiao

指数不可能这么大啊1！我不知道为什么 因为我始终觉得你的程序有问题？？
并且如果指数这么大，哪么这个系统的耗散度可想而之，并且轨线间的分离程度可想而之！

ssaw0920

请问如何判断我的系统的耗散程度呢？
我运行最长时间到20秒，计算20000个点，位移速度曲线近似为周期的，速度时间曲线接近于时间轴，这样能说明系统不耗散吗？或者是我算lyapunov指数时候运行的点数不够多呢？我计算了500个点。
谢谢！

ssaw0920

这是我编的画poincare截面的程序，请学长帮我指导一下看看哪里不太对，麻烦学长了！
%poincare截面绘制
clear;
options=odeset('RelTol',1e-10);
tall=1;
h=0.01;
tt=1;
k=tall/tt;
t=[];
x=[];
x1=[];
x2=[];
x3=[];
n=input('输入n值n=');%输入n的值,n要大于1，因为等于1时一个点都画不出来，因为T=3/20>t=1/10
i=1;
    [t1,x1]=ode15s(@dafen1,[(i-1)*tt:h:1/4*tt+(i-1)*tt],[0,0,0,0,0,0,0,0,0,0],[]);
    t=[t;t1];
    x=[x;x1];
    [t2,x2]=ode15s(@dafen2,[1/4*tt+(i-1)*tt+h:h:3/4*tt+(i-1)*tt],x1((1/4*tt)/h,:),[]);
    t=[t;t2];
    x=[x;x2];
    [t3,x3]=ode15s(@dafen1,[3/4*tt+(i-1)*tt+h:h:tt+(i-1)*tt],x2((1/2*tt)/h,:),[]);
    t=[t;t3]
    x=[x;x3];
    for i=2:n
        [t1,x1]=ode15s(@dafen1,[(i-1)*tt+h:h:1/4*tt+(i-1)*tt],x3((1/4*tt)/h,:),[]);
        t=[t;t1];
        x=[x;x1];
        [t2,x2]=ode15s(@dafen2,[1/4*tt+(i-1)*tt+h:h:3/4*tt+(i-1)*tt],x1((1/4*tt)/h,:),[]);
        t=[t;t2];
        x=[x;x2];
        [t3,x3]=ode15s(@dafen1,[3/4*tt+(i-1)*tt+h:h:tt+(i-1)*tt],x2((1/2*tt)/h,:),[]);
        t=[t;t3]
        x=[x;x3];
        x1;
        x2;
        x3;
    end
x;
T=3/20;
j=1;
ypoin=[];
for k=1:n*tt/T+1
    ypoin(k,:)=x(j,:);
    j=j+T/h;
end
ypoin;

subplot(3,3,1)
plot(ypoin(:,3),ypoin(:,5),'k.')
hold on
subplot(3,3,2)
plot(ypoin(:,4),ypoin(:,6),'k.')
hold on
subplot(3,3,3)
plot(ypoin(: ,7),ypoin(: ,9),'k.')
hold on
subplot(3,3,4)
plot(ypoin(:,8),ypoin(:,10),'k.')
hold on
subplot(3,3,5)
plot(x(:,3),x(:,5),'k')
hold on
subplot(3,3,6)
plot(x(:,4),x(:,6),'k')
hold on
subplot(3,3,7)
plot(x(: ,7),x(: ,9),'k')
hold on
subplot(3,3,8)
plot(x(:,8),x(:,10),'k')

ssaw0920

因为动力学方程比较长，所以只能用附件方式上传，不好意思麻烦学长了！！

ssaw0920

怎么没有反映，麻烦学长帮我看看，谢了先！！！