LET工具箱如何使用?

hohoo

你跟踪调试一下

chuandong418

octopussheng

明天回家了，今天和老板参加了一个研讨会，所以没有上线，先道个歉！

其实在LET工具箱里面有个readme文档，还有help，我建议你好好看看先，其实也不难的！

就对你图里面的那几个问题做一下回答吧！
（1）和（2），这两个其实应该是参考文献里面所用的初始条件和得到的结果，你可以从下面显示的参考文献里面找到，这个就不多说了！

（3）这个就是你求解时需要输入的微分方程的初始条件，呵呵，这个应该没有问题的吧

（4）这个其实是你的线性化方程的数目，从微分方程的维数上来解释，有如下的对应关系
2－－4
3－－9
4－－16
5－－25
6－－36
7－－49
这个就是微分方程子程序中的F＝J×Q的结果
你可以参考一下LET里面的Lorenzq程序，应该也不难理解的！

（5）这个是在当前路径下指定一个输出文件，文件名可以自己设定的！

chuandong418

呵呵，我都看了好几遍了，有几个细节问题不能解决,谢谢OCT.:loveliness:

chuandong418

基本上设置是没问题的，但是为什么会出现#17的错误提示?>

hohoo

界面上 没有问题 看看是不是你的函数有问题

octopussheng

这个主要是我上面回答的第4个问题，你检查一下Jacobian矩阵对不对！

chuandong418

我把函数贴上来，麻烦帮看下，那里出错了,谢谢
function ff=www(t,x,flag,b)
ff=zeros(5,1);
g=9800;
E=0.12;
f=0.12;
a=0.5;
%b=16.0;
q=0.16;
f0=25;
u=0.136/q;
w=2;
e=sqrt(x(1)^2+x(2)^2);
G=g/((2*pi*f0)^2*q);
ff=zeros(5,1);
X=x(1);
Y=x(2);
Z=x(3);
P=x(4);
O=x(5);
dX=x(3);
dY=x(4);
dZ=-2*E*x(3)-x(1)-a*(x(1)^2+x(2)^2)*x(1)-b*(1-1/e)*(x(1)-f*x(2))+u*w^2*cos(w*x(5));
dP=-2*E*x(4)-x(2)-a*(x(1)^2+x(2)^2)*x(2)-b*(1-1/e)*(f*x(1)+x(2))+u*w^2*sin(w*x(5))-G;
dO=1;
J=[0,0,1,0,0;
   0,0,0,1,0;
   -1-(2*a*x(1)^2+a*e^2)-b*(1-1/e),-2*a*x(1)*x(2)+b*f*(1-1/e),-2*E,0,-u*w^3*sin(w*x(5));
   -2*a*x(1)*x(2)-b*f*(1-1/e),-1-3*a*x(2)^2-b*(1-1/e),0,-2*E,u*w^3*cos(w*x(5));
   0,0,0,0,0];
Q=[x(6),x(11),x(16),x(21),x(26);
   x(7),x(12),x(17),x(22),x(27);
   x(8),x(13),x(18),x(23),x(28);
   x(9),x(14),x(19),x(24),x(29);
   x(10),x(15),x(20),x(25),x(30)];
F=J*Q;
ff=[dx(1);dx(2);dx(3);dx(4);dx(5);F(:)];

liliangbiao

原帖由 shenyongjun 于 2007-8-22 17:44 发表

                               
登录/注册后可看大图

这个软件可以计算自己编写的程序啊，点击Run Let Main program，然后选择setting，输入自己编辑的函数文件（按照软件要求的格式），同时进行各种参数设置即可进行计算

能的，确实能自己编写自己的程序，我把这个工具箱和另外一个工具箱 demo chaos 闹混了，谢谢 Shen!

octopussheng

function ff=www(t,x,flag,b)
g=9800;
E=0.12;
f=0.12;
a=0.5;
q=0.16;
f0=25;
u=0.136/q;
w=2;
e=sqrt(x(1)^2+x(2)^2);
G=g/((2*pi*f0)^2*q);
ff=zeros(5,1);
dX=x(3);
dY=x(4);
dZ=-2*E*x(3)-x(1)-a*(x(1)^2+x(2)^2)*x(1)-b*(1-1/e)*(x(1)-f*x(2))+u*w^2*cos(w*x(5));
dP=-2*E*x(4)-x(2)-a*(x(1)^2+x(2)^2)*x(2)-b*(1-1/e)*(f*x(1)+x(2))+u*w^2*sin(w*x(5))-G;
dO=1;
J=[0,0,1,0,0;
   0,0,0,1,0;
   -1-(2*a*x(1)^2+a*e^2)-b*(1-1/e),-2*a*x(1)*x(2)+b*f*(1-1/e),-2*E,0,-u*w^3*sin(w*x(5));
   -2*a*x(1)*x(2)-b*f*(1-1/e),-1-3*a*x(2)^2-b*(1-1/e),0,-2*E,u*w^3*cos(w*x(5));
   0,0,0,0,0];
Q=[x(6),x(11),x(16),x(21),x(26);
   x(7),x(12),x(17),x(22),x(27);
   x(8),x(13),x(18),x(23),x(28);
   x(9),x(14),x(19),x(24),x(29);
   x(10),x(15),x(20),x(25),x(30)];
F=J*Q;
ff=[dX;dY;dZ;dP;dO;F(:)];

你算算试试！

chuandong418

还是提示同样的错误,错误和#17一样.555

chuandong418

难道就没人用这个吗?有用过的指点下,谢谢了!:loveliness:

hohoo

你程序www(t,x,flag,b)中的b没有定义，这个值应该在你程序中定义，let界面上没有b参数的输入，你可以把b赋个值看看

chuandong418

还是回出现#17中出现的错误！:'(

hohoo

不是吧 我运行过了 哦 你这一句ff=zeros(5,1);没去掉吧