|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
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(:)];
可是总是出错,说我设置的参数不对之类的,那位遇到过,帮忙看看那里错了吧。 |
|