我计算Lyapunov指数遇到的问题
请高手帮忙，下面是我计算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=;
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; 正确的LE2=-1.621，我的计算结果是0.218是个正值，具体问题在哪里呢？ 为什么我得不到负值李指数？
页:
[1]