ssaw0920 发表于 2008-7-28 11:06

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

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

siyanyu 发表于 2008-7-28 11:19

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

ssaw0920 发表于 2008-7-28 15:47

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

liliangbiao 发表于 2008-7-28 16:17

指数怎么这么大呢？？

wxl789 发表于 2008-7-29 10:25

回复 4楼 的帖子

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

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

ssaw0920 发表于 2008-7-29 16:13

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

liliangbiao 发表于 2008-7-30 09:46

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

ssaw0920 发表于 2008-7-30 11:04

回复 7楼 的帖子

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

ssaw0920 发表于 2008-7-30 11:18

这是我编的画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;
    =ode15s(@dafen1,[(i-1)*tt:h:1/4*tt+(i-1)*tt],,[]);
    t=;
    x=;
    =ode15s(@dafen2,,x1((1/4*tt)/h,:),[]);
    t=;
    x=;
    =ode15s(@dafen1,,x2((1/2*tt)/h,:),[]);
    t=
    x=;
    for i=2:n
      =ode15s(@dafen1,[(i-1)*tt+h:h:1/4*tt+(i-1)*tt],x3((1/4*tt)/h,:),[]);
      t=;
      x=;
      =ode15s(@dafen2,,x1((1/4*tt)/h,:),[]);
      t=;
      x=;
      =ode15s(@dafen1,,x2((1/2*tt)/h,:),[]);
      t=
      x=;
      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 发表于 2008-7-30 11:26

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

ssaw0920 发表于 2008-8-1 17:57

怎么没有反映，麻烦学长帮我看看，谢了先！！！
页: [1]
查看完整版本: 请问这样的poincare映射图和lyapunov指数图能说明系统是混沌的吗？