声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1683|回复: 2

[分形与混沌] 我计算Lyapunov指数遇到的问题

[复制链接]
发表于 2008-11-16 23:32 | 显示全部楼层 |阅读模式

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

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

x
请高手帮忙,下面是我计算Henon的李指数的MATLAB程序,用的是离散方法。但是计算出来的结果是LE1=log(max(abs(eig(F))))=0.4163   LE2=log(min(abs(eig(F))))=0.2183    实际应该是 LE1 = 0.418, LE2 = -1.621  其中我的LE1计算是正确的,但是LE2的计算结果不正确。请高手帮忙看问题出在哪里。谢谢大家。
%             x(n+1) = 1 - a*x(n)^2 + y(n)
%             y(n+1) = b*x(n)
a=1.4;
b=0.3;
x=zeros(401,1);y=zeros(401,1);
x(1)=0;y(1)=0;
for i=1:100
    x(i+1) = 1 - a*x(i)^2 + y(i);
    y(i+1 )= b*x(i);
end
F=[1 0;0 1];
for i=101:300
    x(i+1) = 1 - a*x(i)^2 + y(i);
    y(i+1 )= b*x(i);
    J=[-2*a*x(i) 1;b 0];
    F=F*J;
end
z=log(max(abs(eig(F))));
LE=z/200;

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

 楼主| 发表于 2008-11-17 10:45 | 显示全部楼层
正确的LE2=-1.621,我的计算结果是0.218是个正值,具体问题在哪里呢?
 楼主| 发表于 2008-11-18 15:39 | 显示全部楼层
为什么我得不到负值李指数?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-26 00:04 , Processed in 0.059803 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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