史上最短最快的李氏指数matlab程序

simber · 发表于 2012-6-7 06:33

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

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

x

本帖最后由牛小贱于 2015-3-15 15:15 编辑

本人相信简单就是美，今天花了一天把之前的程序简化了，其中用QR分解来代替smith正交归一化和向量长度的求解。实测速度比原来循环程序快5倍左右。

程序里用的是henon map，也可以求解微分方程的李氏指数，只需要把Q和y两行用求解微分方程和jacobi矩阵代码代替。

nit = 100000
n=2
y=[0;0];
a=1.4; b=0.3;
Q=eye(n);
cum=zeros(1,n);
Lexp=zeros(nit,n);
for i=1:nit
Q=[-2*a*y(1), 1;b, 0]*Q;
y=[1-a*y(1)^2+y(2); b*y(1)];
[Q R]=qr(Q);
cum=cum+log(diag(R)');
Lexp(i,:)=cum./(i);
end;
plot(Lexp);
Lys=sum(Lexp(end-1999:end,:))/2000

复制代码

simber · 发表于 2012-6-8 05:40

这么好的程序居然没人顶，表示很熵心

dongxuxy · 发表于 2012-6-10 20:00

simber · 发表于 2012-6-11 01:29

dongxuxy 发表于 2012-6-10 20:00

登录/注册后可看大图

发现个规律了，这个论坛越弱智越简单没意义的贴回复的人人越多，比如帖子里面要求回复邮箱发程序这种

chenshumeionly · 发表于 2012-6-19 18:05

看不太懂呢我最近正在找这类的程序我画的是硬弹簧Duffing系统的李亚普洛夫指数图你这个程序能用吗

傲视天下 · 发表于 2012-6-20 00:19

我来顶顶顶顶顶顶顶顶顶顶！

ChaChing · 发表于 2012-6-20 02:37

simber 发表于 2012-6-8 05:40

登录/注册后可看大图

这么好的程序居然没人顶，表示很熵心

我想LZ应该是分享学习成果, 希望有人欣赏并参与讨论精进才是!
并非一味没意义的顶吧:@)

Rainyboy · 发表于 2012-6-20 13:06

本帖最后由 Rainyboy 于 2012-6-20 13:08 编辑

术业有专攻，实际上科研很多事就是近乎孤芳自赏。
况且，仅凭一段程序就表示欣赏的一定是专业的同行了，而这样近的同行（还包括无法看到本贴的）又能有几个呢？我想这也就是科技论文的都要求有详尽的参考文献，摘要和综述（简介）的原因吧，不然能看懂的人真的太少。若有以此工作为基础发表的学术论文，可以发一个名称或者连接或者附件上来，也许大家更能明白你的工作的意义。

ejec · 发表于 2012-7-1 11:16

楼主能不能讲讲给咱扫扫盲我正好想用李氏指数不知道你的代码怎么用谢谢

我这可是注册了3天多才注册上的就想留这个言希望楼主能给回复一下万分感谢

ddlianyi · 发表于 2012-9-20 21:56

看不懂

165390292 · 发表于 2012-9-25 22:01

看不明白啊，楼主解释一下就更完美了。

redplum · 发表于 2012-10-5 21:01

有具体的数学模型和应用就可以

x589 · 发表于 2014-12-24 23:31

顶顶顶可是楼主能不能好事情做到底把程序完善下啊

寒域 · 发表于 2015-3-4 16:08

为什么运行后感觉不太对啊？恒定的LE？

账号		自动登录	找回密码
密码			我要加入

[计算数学] 史上最短最快的李氏指数matlab程序

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

点评

评分

点评

评分