henon_lle_nicolis2.m

来自「实现 henon混沌仿真 实现 henon混沌仿真 实现 henon混沌仿真 实」· M 代码 · 共 27 行

% clear;
close all;

k=300;
for i=1:301
    a=(i-1)*5/1000;
    x(1)=0.6;
    y(1)=0.4;
    dx=[0.002;0.01];
    d0=norm(dx);
    e=1/d0*dx;
    le=0;
    for j=1:k
        J=[-2*a*x(j) 1;0.3 0];
        vec=J*e;
        d=norm(vec);
        e=1/d*vec;
        le=le+1/k*log(d);
        x(j+1)=1+y(j)-a.*x(j).*x(j);
        y(j+1)=0.3*x(j);
    end
    lle(i)=le;
end
a=0:0.005:1.5;
plot(a,zeros(1,301),'k');
hold on
plot(a,lle,'b');