我计算Lyapunov指数遇到的问题

11m

请高手帮忙，下面是我计算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;

11m

正确的LE2=-1.621，我的计算结果是0.218是个正值，具体问题在哪里呢？

11m

为什么我得不到负值李指数？