声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1386|回复: 2

[分形与混沌] 用过LET工具箱的帮忙看一下吧,不胜感激!!

[复制链接]
发表于 2008-1-15 13:31 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
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=[dx;dy;dz,dx1,dy1,dz1];        %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=[DX1; F(:)];

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

使用道具 举报

 楼主| 发表于 2008-1-15 13:33 | 显示全部楼层
:'( 还有我已经搜过原来有关LET的帖子,可是没有特别具体的,所以我还是不太清楚,应该怎么办,麻烦各位指点一下吧。。。。。。。
发表于 2008-1-15 14:04 | 显示全部楼层
把弹出的错误提示窗口贴一下吧,看看先
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-9-21 10:51 , Processed in 0.050491 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表