llp.m

基于定义法计算Lyapunov指数。只要在matlab命令窗口中直接输入文件中函数的名称即可。

function llp
global a;
global b;
global r;
a=16;b=4;r=45.92;
y=[1;1;1;1;0;0;0;1;0;0;0;1];
lp=0;
for k=1:300
    [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/300