大家帮忙看一下我的图，这回传上来了

siyanyu

这是我做的相图和用LET做的最大Lyapunov图，可是却是矛盾的
相图说明是混沌的，可Lyapunov图却表明不是混沌的
为什么会出现这样的情况呢？






[ 本帖最后由 siyanyu 于 2008-7-28 16:33 编辑 ]

liliangbiao

怎么看不到呢？？你没有上传正确的格式吧？

siyanyu

这回能看到了吧:@)
用LET做完的图若是负值不就说明是非混沌吗
可是相图却是混沌的
这是怎么回事啊
帮下忙呵呵
多谢啦

zhailiangjun

你也是用的let工具箱啊，呵呵。咱们遇见的情况差不多。我得到的图是周期的解，用let算出来确是混沌的。期待有人来帮忙啊

liliangbiao

你的指数算的绝对不正确，请检查修改的let程序再算！

zhailiangjun

这个let的工具箱不知道liliangbiao大哥是不是经常使用呐，这个计算的可信度高吗。我计算的是哈密顿系统，所得的4个李雅普诺夫指数之和为0，对称的。但就是数值和别人算得不一样。不知道那地方出问题了。

liliangbiao

我以前用过 现在不用了，我怎么在图像中找不到你说的那个4个李雅普诺夫指数之和为0，对称的.至于你说的那个数值和别人算的不一样，我不知道，但是如果差别不大，近乎相等，即可！

siyanyu

不太明白你说的指数，怎么算？
我就是把方程带进去了，代入初值
就运行了
还望指教一下:handshake

liliangbiao

你可以发给我你的修改的let程序看看！如果方便的话，你可以贴上！lixf@mail.lzjtu.cn

siyanyu

好的，其实我就是把工具箱的方程改了一下，其他的都没有动
这是程序：

%Rearrange input data in desired format
%Note: the input data is a column vector
x=X(1); y=X(2);z=X(3);
Q=[X(4),X(7),X(10);
   X(5),X(8),X(11);
   X(6),X(9),X(12)];

%Duffing's equation
dx=y;
dy=-3.8*x-0.5*(y+10*x^2*y)-5*x^3-2*10.2*cos(2*z)*(-x+x^3/6);
dz=2;                %where z = t, this transformation is for changing
                %the non-autonomous system to a autonomous one

DX1=[dx;dy;dz];        %Output data

%Linearized system
J=[    0,    1,          0;
   -3.8-10*x*y-15*x^2-2*10.2*cos(2*z)*(-1+x^2/2), -0.5-5*x^2,  4*10.2*sin(2*z)*(-x+x^3/6);
        0,   0,          0];

%Variational equation
F=J*Q;

%Put output data in a column vector
OUT=[DX1;F(:)];

liliangbiao

你的Jacobi矩阵不正确啊！！

siyanyu

我同样还有一个也是一样的问题，相图混沌，可是LET的却不是

siyanyu

不正确吗？没发现哪里不对啊？

J=[    0,    1,          0;
-3.8-10*x*y-15*x^2-2*10.2*cos(2*z)*(-1+x^2/2),    -0.5-5*x^2,    4*10.2*sin(2*z)*(-x+x^3/6);
        0,   0,          0];
分别是对x,y,z求的导啊

liliangbiao

怎么对x求导还会有x^2？？？？？

liliangbiao

应该是=-3.8*x-0.5*(y+10*x^2*y)-5*x^3-2*10.2*cos(2*z)*(-x+x^3/6);对x求偏导 是

J=[-3.8-10*x*y-15*x^2-2*10.2*cos(2*z)*(-1+x^2/2]

[ 本帖最后由 liliangbiao 于 2008-7-29 12:01 编辑 ]