用过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=;
可是总是出错,说我设置的参数不对之类的,那位遇到过,帮忙看看那里错了吧。 :'( 还有我已经搜过原来有关LET的帖子,可是没有特别具体的,所以我还是不太清楚,应该怎么办,麻烦各位指点一下吧。。。。。。。 把弹出的错误提示窗口贴一下吧,看看先
页:
[1]