LE问题求助！
小弟初次接触混沌，现在在算一个有干摩擦曲柄滑块机构的LE指数，用的是定义法求最大LE，但是结果不收敛，并且从仿真结果看运动是周期的，但是LE却是正的，不知各位高手有没有碰到这种情况，还请指点一二。求LE的代码如下：function LE=CalLE()
global K;
global LE;
d0=1e-9;
YY1=;
YY2=YY1+;
Endt=0.1;
n=1;
for k=120:120
K(n)=k;
temp=0;
for j=1:2500
=SysSimulation(YY1,50,50,k,Endt);
=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，横轴是迭代次数
页:
[1]