llp.m

Lorenz系统的李雅谱诺夫指数

function llp
global a;
global b;
global r;
a=10;b=8/3;r=30;
y=[1;1;1;1;0;0;0;1;0;0;0;1];
lp=0;
for k=1:100 
    [T,Y] = ode45('Lorenz2',1, y); 
    y = Y(size(Y,1),:);
    y0 = [y(4) y(7) y(10);y(5) y(8) y(11);y(6) y(9) y(12)];
    y0=GS(y0);
    mod(1)=norm(y0(:,1));
    mod(2)=norm(y0(:,2));
    mod(3)=norm(y0(:,3));
    lp = lp+log(abs(mod));
    y0(:,1)=y0(:,1)/mod(1);
    y0(:,2)=y0(:,2)/mod(2);
   y0(:,3)=y0(:,3)/mod(3);
   y(4:12) = y0';
end
lp=lp/100