LE问题求助！

buaalqq

小弟初次接触混沌，现在在算一个有干摩擦曲柄滑块机构的LE指数，用的是定义法求最大LE，但是结果不收敛，并且从仿真结果看运动是周期的，但是LE却是正的，不知各位高手有没有碰到这种情况，还请指点一二。求LE的代码如下：
function LE=CalLE()
global K;
global LE;
d0=1e-9;
YY1=[0;pi/2];
YY2=YY1+[0;d0];
Endt=0.1;
n=1;
for k=120:120
   K(n)=k;
   temp=0;   
    for j=1:2500
       [Time1,Y1]=SysSimulation(YY1,50,50,k,Endt);
       [Time2,Y2]=SysSimulation(YY2,50,50,k,Endt);
        YY1=Y1(end,:);
        YY2=Y2(end,:);
       d=sqrt((YY1-YY2)*(YY1-YY2)');
       YY2=YY1+(YY2-YY1)*d0/d;
       if j>600
           temp=temp+log(d/d0);
           LE(j-600)=temp/(Endt*(j-600));
       end
       disp(['k=',num2str(k),',j=',num2str(j)]);
    end
   
%  n=n+1;
end
为了验证对错，只算了k=120的情况，另外家里电脑性能不好，迭代次数不多，不知是不是这个原因。
结果如下，纵轴是最大LE，横轴是迭代次数

纵轴是最大LE，横轴是迭代次数