yuxiaobai 发表于 2008-1-15 13:31

用过LET工具箱的帮忙看一下吧,不胜感激!!

LET工具箱不是可以使用自己编写的混沌系统的程序吗?我自己改写了原系统中的lorenz混沌系统,把它变成了六个常微分方程的,也就是两个振子的同步问题。然后另存为lorenzeqsys.m然后使用LET工具箱时进行了如下设置:
ODE function:lorenzeqsys
initial Conditions:1 1 1 2 2 2
No. of linearized ODEs:36
下面是改完之后的程序:

SIGMA = 10;
RHO = 28;
BETA = 8/3;

%Rearrange input data in desired format
%Note: the input data is a column vector
x=X(1);y=X(2);z=X(3);x1=X(4);y1=X(5);z1=X(6);
Q=[ X(7), X(13), X(19), X(25), X(31), X(37);
    X(8), X(14), X(20), X(26), X(32), X(38);
    X(9), X(15), X(21), X(27), X(33), X(39);
    X(10), X(16), X(22), X(28), X(34), X(40);
    X(11), X(17), X(23), X(29), X(35), X(41);
    X(12), X(18), X(24), X(30), X(36), X(42)];
%Lorenz equation
dx=SIGMA*(y-x);
dy=-x*z+RHO*x-y;
dz=x*y-BETA*z;

dx1=SIGMA*(y1-x1);
dy1=-x*z1+RHO*x-y1;
dz1=x*y1-BETA*z1;


DX1=;        %Output data

%Linearized system
J=[-SIGMA, SIGMA,   0,       0,      0,0;
   RHO-z,    -1,    -x;       0,      0,0;
         y,   x, -BETA      0,      0,0;
         0,   0,   0,-SIGMA,SIGMA,0;
    RHO-z1,   0,   0,       0,   -1, -x;
      y1,   0,   0,       0,      x, -BETA];

%Variational equation   
F=J*Q;

%Output data must be a column vector
OUT=;

可是总是出错,说我设置的参数不对之类的,那位遇到过,帮忙看看那里错了吧。

yuxiaobai 发表于 2008-1-15 13:33

:'( 还有我已经搜过原来有关LET的帖子,可是没有特别具体的,所以我还是不太清楚,应该怎么办,麻烦各位指点一下吧。。。。。。。

octopussheng 发表于 2008-1-15 14:04

把弹出的错误提示窗口贴一下吧,看看先
页: [1]
查看完整版本: 用过LET工具箱的帮忙看一下吧,不胜感激!!