声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

楼主: chuandong418

[分形与混沌] LET工具箱如何使用?

[复制链接]
 楼主| 发表于 2007-8-24 18:51 | 显示全部楼层

回复 #30 hohoo 的帖子

出现了以下一系列错误,55555
??? Input argument "b" is undefined.

Error in ==> www at 23
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));

Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

Error in ==> findlyap at 164
      eval(IntegrationStr);

Error in ==> let at 242
  LineHandles=findlyap(MainHandles);

??? Error using ==> let('start')
Input argument "b" is undefined.

??? Error while evaluating uicontrol Callback
回复 支持 反对
分享到:

使用道具 举报

 楼主| 发表于 2007-8-24 18:52 | 显示全部楼层

回复 #30 hohoo 的帖子

第几个ff=zero(5,1)要去掉啊?
发表于 2007-8-24 18:53 | 显示全部楼层
你有qq吗 这样说太费劲了
 楼主| 发表于 2007-8-24 18:56 | 显示全部楼层

回复 #33 hohoo 的帖子

50665321
发表于 2007-8-24 18:56 | 显示全部楼层
function ff=www(t,x)
b=1;
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);

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(:)];
 楼主| 发表于 2007-8-25 17:13 | 显示全部楼层
这oct怎么还不回来哦,问题还解决不了呢!:@(
 楼主| 发表于 2007-8-28 20:42 | 显示全部楼层
难道没人用过?:@(
发表于 2008-4-8 17:36 | 显示全部楼层

可否传我个LET工具箱??

我的EMAIL: 150283138@qq.com
不甚感激!到处都没有下到:@(
发表于 2008-4-8 20:42 | 显示全部楼层
问题解决了?
发表于 2008-4-14 18:02 | 显示全部楼层
问题解决了么?我也遇到相同的问题了
发表于 2008-4-14 18:27 | 显示全部楼层

回复 39楼 的帖子

oct上面有很多问题等着你回答呢
发表于 2008-4-14 20:29 | 显示全部楼层
现在不知道回答哪一个了,呵呵!
发表于 2008-4-26 17:35 | 显示全部楼层

回复 42楼 的帖子

请otc回答let工具箱的问题。。。

我初学混沌,还在摸索,,

let工具箱也不太懂。。

谢谢!!!
发表于 2008-4-26 18:23 | 显示全部楼层
let工具箱其实用起来还是蛮简单的

1. 适用于连续系统,对时间序列的LE求解不适用
2. 在进行LE求解之前,需要注意应将非自治系统写成自治系统的形式,然后参考工具箱给出的Lorenz、Rossler系统的例子,将微分方程定义函数写成标准形式
3. 用let求解Lyapunov指数,在设置窗口中设置相关参数即可!
发表于 2008-4-27 17:17 | 显示全部楼层

let工具箱使用

你的程序是正确的。。
把b前边的%去掉就可以运行了。。。
我已经运行通过,,
结果如图所示:

[ 本帖最后由 wxl789 于 2008-4-27 17:22 编辑 ]
let 计算结果.bmp
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 04:10 , Processed in 0.081200 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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