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

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

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

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

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

y=;
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=;
=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

这么好的程序居然没人顶，表示很熵心{:3_47:}

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

{:{39}:}

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

dongxuxy 发表于 2012-6-10 20:00 static/image/common/back.gif


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

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

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

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

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

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

simber 发表于 2012-6-8 05:40 static/image/common/back.gif
这么好的程序居然没人顶，表示很熵心

我想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

{:{19}:}看不懂

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

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

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

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

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

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

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

为什么运行后感觉不太对啊？恒定的LE？
页: [1]
查看完整版本: 史上最短最快的李氏指数matlab程序