应用QR法求Lyapunov指数

laozhou

下面是我用QR法画Duffing方程的Lyapunov指数，但结果分析其中一指数一直大于零，与实际不符，望各位高手给与指点！
--------------------------------------------------------------------------------------------------------
clear;clc;
tic;
w0=1;a=1;b=-1;c=1;
T=2*pi/w0;
F=0:0.01:1;
for i=1:length(F)
    [t,x]=ode45('my_duffing',[0:T/600:200*T],[0,0.0001],[],a,b,c,F(i),w0);
   N=length(x(:,1));
   L1=0;
   L2=0;
   Q=eye(2,2);
   for j=1:N
       A=[0,1;-b-3*c*x(j,1).^2,-a];
       [Q,R]=qr(A*Q);
       R=diag(R);
       L1=L1+log(abs(R(1)));
       L2=L2+log(abs(R(2)));
   end
   L11(i)=L1/N;
   L22(i)=L2/N;
   end
plot(F,L),grid on;
------------------------------------------------------------
M函数

function dx=my_duffing(t,x,flag,a,b,c,F,w0)
dx=[x(2);-a*x(2)-b*x(1)-c*x(1).^3+F*cos(w0*t)];
---------------------------------------------------------------

lushuang8

实际结果是什么样的？